Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONTRACTS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Guided Ion Beam and Quantum Chemical Investigation of the Thermochemistry of Thorium Dioxide Cations: Thermodynamic Evidence for Participation of f Orbitals in Bonding

Kinetic energy dependent reactions of ThO + with O 2 are studied using a guided ion beam tandem mass spectrometer. The formation of ThO 2 + in the reaction of ThO + with O 2 is observed to be slightly endothermic and also exhibits two obvious features in the cross section. These kinetic energy dependent cross sections were modeled to determine a 0 K bond dissociation energy of D 0 (OTh + - O) = 4.94 ± 0.06 eV. This value is slightly larger but within experimental uncertainty of less precise previously reported experimental values. The higher energy feature in the ThO 2 + cross section was also analyzed and suggests formation of an excited state of the product ion lying 3.1 ± 0.2 eV above the ground state. Additionally, the thermochemistry of ThO 2 + was explored by quantum chemical calculations, including a full Feller-Peterson-Dixon (FPD) composite approach with correlation contributions up to CCSDT(Q) and four-component spin-orbit corrections, as well as more approximate CCSD(T) calculations including semiempirical estimates of spin-orbit energy contributions. The FPD approach predicts D 0 (OTh + -O) = 4.87 ± 0.04 eV, in good agreement with the experimental value. Analogous FPD results for ThO+, ThO, and ThO 2 are also presented, including ionization energies for both ThO and ThO 2 . The ThO 2 + bond energy is larger than those of its transition metal congeners, TiO 2 + and ZrO 2 + , which can be attributed partially to an actinide contraction, but also to contributions from the participation of f orbitals on thorium that are unavailable to the transition metal systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Quantum-Based Molecular Dynamics Simulations Using Tensor Cores

Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics

Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. Furthermore, the situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version

An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin–orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S o → S 1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. Here, a UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗
DOE OSTI2020

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗
DOE OSTI2020

f-Orbital Mixing in the Octahedral f 2 Compounds UX 6 2– [X = F, Br, Cl, I] and PrCl 6 3–

Understanding how interactions between the f orbitals and ligand orbitals in lanthanide and actinide systems affect their physical properties is the central issue in f-element chemistry. A wide variety of approaches including both theoretical and experimental tools have been used to study these relationships. Among the most widely used tools has been crystal field theory (CFT), which bridges theory and experiment in that it is a model based largely on atomic theory that is parametrized using experimental data. Crystal field theory is quite accurate for the lanthanides, due in part to the highly contracted nature of the 4f orbitals. For actinides, crystal field theory is less accurate, potentially due to the treatment of orbital mixing. In CFT, orbital mixing is handled implicitly by allowing the electron repulsion parameters (Slater F k parameters) and the spin–orbit coupling constant to vary. As a result, orbital mixing in CFT is isotropic in that the F k parameters and the spin–orbit coupling constant affect all f orbitals equally. This approximation works well for the lanthanides due to the limited degree of orbital mixing in these complexes. In actinide complexes, the 5f orbitals have greater overlap with the ligand orbitals, and this approximation is less accurate than in the lanthanides. In this work, we report a modification of CFT that includes the effect of orbital mixing on electron repulsion and spin–orbit coupling for each f orbital. The model is applied to the tetravalent uranium hexahalide dianions and PrCl 6 3– for which the energies of many low-lying excited states are known. The new model generally fits the data as well the traditional CFT although with fewer parameters. However, the new model does not fit the data better than the more complex CFT models of Faucher and co-workers. The results of the model show in detail how changes in overlap and orbital energies influence the energies of the bonding and antibonding orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Cryogenic Vibrationally Resolved Photoelectron Spectroscopy of OH–(H2O): Confirmation of Multidimensional Franck-Condon Simulation Results for the Transition State of the OH + H2O Reaction

We present a transition state spectroscopic study of the OH + H2O reaction using the experimental technique of cryogenic Negative Ion Photoelectron Spectroscopy (NIPES). The recorded NIPE spectrum at 193 nm exhibits multiple vibrational progressions that include excitations to the shared H atom antisymmetric stretching mode with an interval of 0.32 eV as well as other progressions, mainly involving the H bending and O…O symmetric stretching modes. The Vertical Detachment Energy (VDE) was measured at 3.53 eV, whereas an upper limit for the Adiabatic Detachment Energy (ADE) was estimated at 2.90 eV. These values are in excellent agreement with the theoretically computed values of 3.51 eV and 2.83 eV respectively, obtained at the CCSD(T)/aug-cc-pV5Z level of theory. The recorded NIPE spectrum is in very good agreement when compared to the one recently reported from four-dimensional Franck-Condon simulations, in which a similar spectral profile was predicted. Besides observing the ground state, a charge-transfer excited state in the form of [OH–(H2O)+] is identified with a relative energy of 1.39 eV, well matching the previous prediction of 1.36 eV. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer. This research also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Cao, Wenjin↗
DOE OSTI2022

Experimental X-ray Charge-Density Studies–A Suitable Probe for Superconductivity? A Case Study on MgB 2

Case studies of 1T-TiSe 2 and YBa 2 Cu 3 O 7-δ have demonstrated that X-ray diffraction (XRD) studies can be used to trace even subtle structural phase transitions which are inherently connected with the onset of superconductivity in these benchmark systems. However, the utility of XRD in the investigation of superconductors like MgB 2 lacking an additional symmetry-breaking structural phase transition is not immediately evident. Nevertheless, high-resolution powder XRD experiments on MgB 2 in combination with maximum entropy method analyses hinted at differences between the electron density distributions at room temperature and 15 K, that is, below the T c of approx. 39 K. The high-resolution single-crystal XRD experiments in combination with multipolar refinements presented here can reproduce these results but show that the observed temperature-dependent density changes are almost entirely due to a decrease of atomic displacement parameters as a natural consequence of a reduced thermal vibration amplitude with decreasing temperature. Our investigations also shed new light on the presence or absence of magnesium vacancies in MgB 2 samples–a defect type claimed to control the superconducting properties of the compound. Here, we propose that previous reports on the tendency of MgB 2 to form non-stoichiometric Mg 1–x B 2 phases (1 – x ~ 0.95) during high-temperature (HT) synthesis might result from the interpretation of XRD data of insufficient resolution and/or usage of inflexible refinement models. Indeed, advanced refinements based on an Extended Hansen–Coppens multipolar model and high-resolution X-ray data, which consider explicitly the contraction of core and valence shells of the magnesium cations, do not provide any significant evidence for the formation of non-stoichiometric Mg 1–x B 2 phases during HT synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Introducing GPU Acceleration into the Python-Based Simulations of Chemistry Framework

We introduce the first version of GPU4P Y SCF, a module that provides GPU acceleration of methods in P Y SCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using the Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree–Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of 2 orders of magnitude with respect to the multithreaded CPU Hartree–Fock code of P Y SCF and the performance comparable to other open-source GPU-accelerated quantum chemical packages, including GAMESS and QUICK, on a single NVIDIA A100 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Implementation of McMurchie–Davidson Algorithm for Gaussian AO Integrals Suited for SIMD Processors

We report an implementation of the McMurchie− Davidson evaluation scheme for 1- and 2-particle Gaussian AO integrals designed for processors with Single Instruction Multiple Data (SIMD) instruction sets. Like in our recent MD implementation for graphical processing units (GPUs) [Asadchev, A.; Valeev, E. F.. J. Chem. Phys. 2024, 160, 244109.], variable-sized batches of shellsets of integrals are evaluated at a time. By optimizing for the floating point instruction throughput rather than minimizing the number of operations, this approach achieves up to 50% of the theoretical hardware peak FP64 performance for many common SIMD-equipped platforms (AVX2, AVX512, NEON), which translates to speedups of up to 30 over the state-of-the-art one-shellset-at-a-time implementation of Obara−Saika-type schemes in Libint for a variety of primitive and contracted integrals. As with our previous work, we rely on the standard C++ programming language such as the std::simd standard library feature to be included in the 2026 ISO C++ standard without any explicit code generation to keep the code base small and portable. The implementation is part of the open source LibintX library freely available at https://github.com/ValeevGroup/libintx.

Basis sets↗
DOE OSTI2026

Unimolecular Dynamics of Partially Deuterated Hydroperoxyalkyl Intermediates (•QOOD) in Cyclohexane and Cyclopentane Oxidation

The hydroperoxyalkyl radicals (•QOOH) formed in cyclohexane and cyclopentane oxidation with carbon radical center at the β site relative to the hydroperoxy group have recently been observed through their infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and cyclic ether products. Partial deuteration shifts the IR transitions associated with OH/OD motion to lower frequencies, providing a new window on the IR spectroscopy and unimolecular dynamics of these transient intermediates. Specifically, the overtone OD stretch (2ν OD ) has been identified at 5215.0 cm –1 and 5216.0 cm –1 for the β-QOOD intermediates in cyclohexane and cyclopentane oxidation, respectively. Unimolecular decay rates of the •QOOD intermediates at these energies have been obtained from the time-resolved appearance of OD products. The experimental rates are compared with statistical microcanonical rates evaluated using RRKM theory, including heavy-atom tunneling associated with simultaneous O–O bond elongation and C–O–O angle contraction along the reaction pathway. The experimental rates provide a further test of the previously determined transition state (TS) barriers, which were computed utilizing a benchmarking approach that builds on higher-level reference calculations for the smaller ethane oxidation system. Here, the experimental rate measurements for both •QOOD intermediates agree well with the computed rates after accounting for small changes in zero-point energy and a minor empirical adjustment of the TS barriers in both systems.

Ethers↗
DOE OSTI2024

Ultrafast X-ray Absorption Spectroscopy Reveals Excited-State Dynamics of B 12 Coenzymes Controlled by the Axial Base

Polarized time-resolved X-ray absorption spectroscopy at the Co K-edge is used to probe the excited-state dynamics and photolysis of base-off methylcobalamin and the excited-state structure of base-off adenosylcobalamin. For both molecules, the final excited-state minimum shows evidence for an expansion of the cavity around the Co ion by ca. 0.04 to 0.05 Å. The 5-coordinate base-off cob(II)alamin that is formed following photodissociation has a structure similar to that of the 5-coordinate base-on cob(II)alamin, with a ring expansion of 0.03 to 0.04 Å and a contraction of the lower axial bond length relative to that in the 6-coordinate ground state. Here these data provide insights into the role of the lower axial ligand in modulating the reactivity of B 12 coenzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Structural, Spectroscopic, and Kinetic Insight into the Heating Rate Dependence of Studtite and Metastudtite Dehydration

Uranyl peroxide minerals studtite (UO 2 O 2 ·4H 2 O) and metastudtite (UO 2 O 2 ·2H 2 O) are important materials in the nuclear fuel cycle. When heated, they dehydrate and transform to amorphous uranium oxides (UO x ), yet phase stability and heating rate dependence of these transformations are poorly understood. This information is critical to proper management of fuel cycle materials. In this work, we use in situ powder X-ray diffraction (PXRD), Raman spectroscopy, and thermogravimetric analysis (TGA) to monitor the dehydration of studtite and metastudtite. Strong linear correlation between the heating rate and phase transition temperature is observed. Geometric contraction and diffusion-related kinetic models describe studtite dehydration at slow heating rates, whereas Avrami–Erofeev or reaction order models become more accurate for faster thermal treatments. A second order model describes the transition from metastudtite to UO x regardless of the heating rate. Water retention during studtite dehydration is indicated by PXRD, Raman spectroscopy, and TGA. Here, we observe mixed-phase UO x dehydration products of metastudtite with a likely formation mechanism involving conversion of some uranyl centers from hexagonal to pentagonal bipyramidal coordination units via peroxide liberation. Our observations clarify over 100 years of measurements on these materials and represent an advancement in understanding the chemical behavior of nuclear fuel cycle materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Formation of Vertically Oriented Channels during Calcination of Surfactant-Templated Titania-Doped Mesoporous Silica Thin Films

Here, the mesoporous structure development during calcination of titania-doped mesoporous silica thin films is characterized by in situ grazing-incidence small angle x-ray scattering (GISAXS). Varying amounts of titania (≤ 6 mol%) are incorporated in cetyltrimethylammonium bromide (CTAB) templated film by two methods in this study: dispersion of titania in the silica matrix, and fixing titania at the pore surface by complexation of the titania precursor with a sugar-based co-surfactant, dodecyl maltoside (C 12 G 2 ). Before calcination, films templated with only CTAB display a mesoporous structure with a mixture of in-plane 2D hexagonal and randomly oriented porous structure, while C 12 G 2 -complexed film possessed randomly oriented mesostructure only. Both methods yield a final structure of hexagonal close-packed pores orthogonally oriented to the substrate after calcination at 500 °C. This is attributed to weak Si-O-Ti bonds which allow sintering to occur at 500 °C, combined with unidirectional thermal contraction of the film in the vertical direction. Anisotropic stress and annealing of the films allows the randomly oriented pores to merge vertically and form channels with pore diameter of 2.3 ± 0.3 nm. Titania doping greater than 1% is required for this transformation in films with no C 12 G 2 surfactant used while the transformation was observed for films with C 12 G 2 complexation even at low (0.005%) titania doping, perhaps accentuated by orientational disorder introduced by the sugar surfactant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗