Ethylene Dehydroaromatization over Ga-ZSM-5 Catalysts: Nature and Role of Gallium Speciation
Explore the source record for details and available documents.
99 GENERAL AND MISCELLANEOUS↗
SEARCH · Search NASA
Results for “optimized energy density functionals”
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.
Triplet Tuning: A Novel Family of Non-Empirical Exchange–Correlation Functionals
In the framework of density functional theory (DFT), the lowest triplet excited state (T1) can be evaluated using multiple formulations, the most straightforward of which are unrestricted density functional theory (UDFT) and time-dependent density functional theory (TDDFT). Assuming the exact exchange-correlation (XC) functional is applied, UDFT and TDDFT provide identical energies for T1 ( ET), which is also a constraint that we require our XC functionals to obey. However, this condition is not satisfied by most of the popular XC functionals, leading to inaccurate predictions of low-lying, spectroscopically and photochemically important excited states, such as T1 and the lowest singlet excited state (S1). Inspired by the optimal tuning strategy for frontier orbital energies [T. Stein, L. Kronik, and R. Baer, J. Am. Chem. Soc. 2009, 131, 2818], we proposed a novel and nonempirical prescription of constructing an XC functional in which the agreement between UDFT and TDDFT in ET is strictly enforced. Referred to as "triplet tuning", our procedure allows us to formulate the XC functional on a case-by-case basis, using the molecular structure as the exclusive input, without fitting to any experimental data. The first triplet tuned XC functional, TT-ωPBEh, is formulated as a long-range-corrected (LRC) hybrid of Perdew-Burke-Ernzerhof (PBE) and Hartree-Fock (HF) functionals [M. A. Rohrdanz, K. M. Martins, and J. M. Herbert, J. Chem. Phys. 2009, 130, 054112] and tested on four sets of large organic molecules. Compared to existing functionals, TT-ωPBEh manages to provide more accurate predictions for key spectroscopic and photochemical observables, including but not limited to ET, the optical band gap ( ES), the singlet-triplet gap (Δ EST), and the vertical ionization potential ( I⊥), as it adjusts the effective electron-hole interactions to arrive at the correct excitation energies. This promising triplet tuning scheme can be applied to a broad range of systems that were notorious in DFT for being extremely challenging.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials
Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride
Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.
Sassi, Michel J. [Pacific Northwest National Labor↗
Transitioning NWChem to the Next Generation of Manycore Machines
The NorthWest chemistry (NWChem) modeling software is a popular molecular chemistry simulation software that was designed from the start to work on massively parallel processing supercomputers [1-3]. It contains an umbrella of modules that today includes self-consistent eld (SCF), second order Møller-Plesset perturbation theory (MP2), coupled cluster (CC), multiconguration self-consistent eld (MCSCF), selected conguration interaction (CI), tensor contraction engine (TCE) many body methods, density functional theory (DFT), time-dependent density functional theory (TDDFT), real-time time-dependent density functional theory, pseudopotential plane-wave density functional theory (PSPW), band structure (BAND), ab initio molecular dynamics (AIMD), Car-Parrinello molecular dynamics (MD), classical MD, hybrid quantum mechanics molecular mechanics (QM/MM), hybrid ab initio molecular dynamics molecular mechanics (AIMD/MM), gauge independent atomic orbital nuclear magnetic resonance (GIAO NMR), conductor like screening solvation model (COSMO), conductor-like screening solvation model based on density (COSMO-SMD), and reference interaction site model (RISM) solvation models, free energy simulations, reaction path optimization, parallel in time, among other capabilities [4]. Moreover, new capabilities continue to be added with each new release.
Bylaska, Eric J.↗
Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen
Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.
first-principles calculations↗
Active Site Engineering via Optimizing the Heterogeneous Support Structure for Single-Atom Catalysis
Supported single-atom catalysts show a large range of activities and selectivities that depend on the local environment of the catalytic sites. Here a theory-based optimization strategy is presented that is based on a density functional theory determination of the transition states and intermediates for a low-dimensional coordinate representation of the heterogeneity of the active sites. The approach is applied to a vanadium catalyst on an amorphous SiO 2 support that involves a large kinetic network described using a full chemistry model. Without assuming a priori scaling relations or mechanism reduction, the optimal state of heterogeneity is found to lie at atomic configurations where the activation energies for two distinct key chemical processes are equal. It is found a posteriori that the behavior of the system is consistent with linear free energy scaling relations in the randomness parameters. The energetic span theory proves quite useful in reducing the full chemistry model to a small number of key reactions. The use of a nonlinear optimization algorithm in combination with energetic span theory provides significant simplification in treating disordered systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Possibility of controlling the chemical pattern of He and Ar “guided streamers” by means of N{sub 2} or O{sub 2} additives
Plasma chemistry at atmospheric pressure and low temperature has met an increased interest due to the effective production of reactive species (mainly reactive nitrogen species and reactive oxygen species) which are demanded in many applications like in biomedicine. Τhe generation of such emissive reactive species is herein studied. A single parameter, i.e., the gas composition, is tuned and abundant chemical species are obtained in a quite controllable manner. The system refers to a small dielectric barrier discharge based reactor fed with He–N{sub 2}/O{sub 2} or Ar–N{sub 2}/O{sub 2} gases, which provides plasma in the form of “guided streamers.” The plasma is sustained by positive high voltage pulses, and the emissive transitions versus the gas composition is determined by UV-VIS optical emission spectroscopy. The relative intensities of emissive species are recorded in the ambient air where the streamers propagate. The evolution of dominant species as a function of the gas composition is mapped, and optimal conditions in terms of species production are clearly revealed. The results show that additives do not necessarily enhance the density of all species, implying the need for plasma chemistry optimization in respect to every application. The study is extended inside the reactor, supporting the above statements. Finally, the rotational and vibrational distributions of critical probe molecules are recorded, testing the dependence of the gas temperature and energy transfer, respectively, on the gas composition.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Self-assembling of formic acid on the partially oxidized p (2×1) Cu(110) surface reconstruction at low coverages
Carbon dioxide (CO 2 ) reduction for synthetic fuel generation could be an integral part of a sustainable energy future. Copper (Cu) is the leading electrocatalyst for CO 2 reduction to produce multiple C-containing products such as C1 and C2 hydrocarbons and oxygenates. Understanding the mechanisms leading to their production could help optimize these pathways further. Adsorption studies of the many possible intermediates on well-characterized surfaces are crucial to elucidating these mechanisms. In this work, we explore the adsorption configurations of formic acid (HCOOH) on the surface of the partially oxidized p (2 × 1) reconstruction of the Cu(110) surface, using low-temperature scanning tunneling and atomic force microscopy, in conjunction with density functional theory modeling. We find that HCOOH adsorbs favorably on the CuO chain comprising the reconstruction. The adsorption interactions involve dative bonding of the carbonyl O to the oxidized Cu and hydrogen bonding of the OH group to the surface O or to an adjacently adsorbed HCOOH molecule. Furthermore, cooperative adsorption of the molecules occurs, forming two- to three-molecule-long oligomer chains, facilitated by intermolecular hydrogen bonding and mutual polarization of the CuO acid-base adsorption sites.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Electron Injection from Copper Diimine Sensitizers into TiO 2 : Structural Effects and Their Implications for Solar Energy Conversion Devices
Copper(I) diimine complexes have emerged as low cost replacements for ruthenium complexes as light sensitizers and electron donors, but their shorter metal-to-ligand-charge-transfer (MLCT) states lifetimes and lability of transient Cu(II) species impede their intended functions. Two carboxylated Cu(I) bis-2,9-diphenylphenanthroline (dpp) complexes [Cu(I)(dpp-O(CH2CH2O)(5))(dpp-(COOH)(2))](+) and [Cu(I)(dpp-O(CH2CH2O)(5))(dpp-(F-COOH)(2))](+) (F = tolyl) with different linker lengths were synthesized in which the MLCT-state solvent quenching pathways are effectively blocked, the lifetime of the singlet MLCT state is prolonged, and the transient Cu(II) ligands are stabilized. Aiming at understanding the mechanisms of structural influence to the interfacial charge transfer in the dye-sensitized solar cell mimics, electronic and geometric structures as well as dynamics for the MLCT state of these complexes and their hybrid with TiO2 nanoparticles were investigated using optical transient spectroscopy, X-ray transient absorption spectroscopy, time-dependent density functional theory, and quantum dynamics simulations. The combined results show that these complexes exhibit strong absorption throughout the visible spectrum due to the severely flattened ground state, and a long-lived charge-separated Cu(II) has been achieved via ultrafast electron injection (<300 fs) from the 1MLCT state into TiO2 nanoparticles. The results also indicate that the TiO2-phen distance in these systems does not have significant effect on the efficiency of the interfacial electron-transfer process. The mechanisms for electron transfer in these systems are discussed and used to develop new strategies in optimizing copper(I) diimine complexes in solar energy conversion devices.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structural anisotropy in Sb thin films
Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.
36 MATERIALS SCIENCE↗
Shape Matters: Understanding the Effect of Electrode Geometry on Cell Resistance and Chemo-Mechanical Stress
Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities when compared to a conventional geometry, i.e. a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher active material loading, while keeping ion transport distances small. However, the relationship between electrical and mechanical performance of shaped electrodes remains poorly understood. Many electrode designs have been explored, where the electrodes are individually shaped or intertwined, and advances in manufacturing and shape/topology optimization have made such designs a reality. Here, we explore sinusoidal half cells and interdigitated full cells. First, we use a simple electrostatics model to understand the cell resistance as a function of shape. Here we focus on low-temperature conditions, where the electrolyte conductivity decreases relative to that of the electrode; here, LiPF 6 EC:DMC electrolyte and MnO 2 electrode are considered. Next, we use a chemo-mechanics model to examine the stress that arises due to intercalation-driven volume expansion. We show that shaped electrodes provide a significant reduction in resistance in low-temperature conditions, however, they exhibit unfavorable stress concentrations. Overall, we find that the fully interdigitated electrodes may provide the best balance with respect to this resistance-stress trade-off.
25 ENERGY STORAGE↗
Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics
Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty
The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.
58 GEOSCIENCES↗
Investigating the solid-liquid phase transition of water nanofilms using the generalized replica exchange method
The generalized Replica Exchange Method (gREM) was applied to study a solid-liquid phase transition in a nanoconfined bilayer water system using the monatomic water (mW) model. Exploiting optimally designed non-Boltzmann sampling weights with replica exchanges, gREM enables an effective sampling of configurations that are metastable or unstable in the canonical ensemble via successive unimodal energy distributions across phase transition regions, often characterized by S-loop or backbending in the statistical temperature. Extensive gREM simulations combined with Statistical Temperature Weighted Histogram Analysis Method (ST-WHAM) for nanoconfined mW water at various densities provide a comprehensive characterization of diverse thermodynamic and structural properties intrinsic to phase transitions. Graph representation of minimized structures of bilayer water systems determined by the basin-hopping global optimization revealed heterogeneous ice structures composed of pentagons, hexagons, and heptagons, consistent with an increasingly ordered solid phase with decreasing density. Apparent crossover from a first-order solid-liquid transition to a continuous one in nanoconfined mW water with increasing density of the system was observed in terms of a diminishing S-loop in the statistical temperature, smooth variation of internal energies and heat capacities, and a characteristic variation of lateral radial distribution functions, and transverse density profiles across transition regions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Non-reactive facet specific adsorption as a route to remediation of chlorinated organic contaminants
The present work quantifies metal-contaminant interactions between palladium substrates and three salient chlorinated organic contaminants, namely trichloroethylene 1,3,5-trichlorobenzene (TCB), and 3,3′,4,4′-tetrachlorobiphenyl (PCB77). Given that Pd is one of the conventional catalytically active materials known for contaminant removal, maximizing catalytic efficiency through optimal adsorption dynamics reduces the cost of remediation of contaminants that are persistent water pollutants chronically affecting public health. Adsorption efficiency analyses from all-atom molecular dynamics (MD) simulations advance the understanding of reaction mechanisms available from density functional theory (DFT) calculations to an extractable feature scale that can fit the parametric design of supported metal catalytic systems and feed into high throughput catalyst selection. Data on residence time, site-specific adsorption, binding energies, packing geometries, orientation profiles, and the effect of adsorbate size show the anomalous behaviour of organic contaminant adsorption on the undercoordinated {110} surface as compared to the {111} and {100} surfaces. The intermolecular interaction within contaminants from molecular dynamics simulation exhibits refreshing results than ordinary single molecule density functional theory calculation. Since complete adsorption and dechlorination is an essential step for chlorinated organic contaminant remediation pathways, the presented profiles provide essential information for designing efficient remediation systems through facet-controlled palladium nanoparticles.
Guo, Hao↗
A charge-density-based general cation insertion algorithm for generating new Li-ion cathode materials
Future lithium (Li) energy storage technologies, in particular solid-state configurations with a Li metal anode, opens up the possibility of using cathode materials that do not necessarily contain Li in its as-made state. To accelerate the discovery and design of such materials, we develop a general, chemically, and structurally agnostic methodology for identifying the optimal Li sites in any crystalline material. For a given crystal structure, we attempt multiple Li insertions at symmetrically in-equivalent positions by analyzing the electronic charge density obtained from first-principles density functional theory. In this report, we demonstrate the effectiveness of this procedure in successfully identifying the positions of the Li ion in well-known cathode materials using only the empty host (charged) material as guidance. Furthermore, applying the algorithm to over 2000 candidate cathode empty host materials we obtain statistics of Li site preferences to guide future developments of novel Li-ion cathode materials, particularly for solid-state applications.
25 ENERGY STORAGE↗
Efficient Oxygen Electrocatalysis by Nanostructured Mixed-Metal Oxides
Oxygen electrocatalysis plays a critical role in the efficiency of important energy conversion and storage systems. While many efforts have focused on designing efficient electrocatalysts for these processes, optimal catalysts that are inexpensive, active, selective, and stable are still being searched. Nonstoichiometric, mixed-metal oxides present a promising group of electrocatalysts for these processes due to the versatility of the surface composition and fast oxygen conducting properties. Herein, we demonstrate, using a combination of theoretical and experimental studies, the ability to develop design principles that can be used to engineer oxygen electrocatalysis activity of layered, mixed ionic-electronic conducting Ruddlesden– Popper (R–P) oxides. We show that a density function theory (DFT) derived descriptor, O 2 binding energy on a surface oxygen vacancy, can be effective in identifying efficient R–P oxide structures for oxygen reduction reaction (ORR). Using a controlled synthesis method, well-defined nanostructures of R–P oxides are obtained, which along with thermochemical and electrochemical activity studies are used to validate the design principles. This has led to the identification of a highly active ORR electrocatalyst, nanostructured Co-doped lanthanum nickelate oxide, which when incorporated in solid oxide fuel cell cathodes significantly enhances the performance at intermediate temperatures (~550 °C), while maintaining long-term stability. Here, the reported findings demonstrate the effectiveness of the developed design principles to engineer mixed ionic-electronic conducting oxides for efficient oxygen electrocatalysis, and the potential of nanostructured Co-doped lanthanum nickelate oxides as promising catalysts for oxygen electrocatalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗