Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cost function”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Employing a Hardware-in-the-Loop Approach to Realize a Fully Homomorphic Controller for a Small Modular Advanced High Temperature Reactor

This paper addresses the cybersecurity challenges of advanced nuclear reactors by integrating fully homomorphic encryption (FHE) into their control systems, enabling encrypted processing of control signals without compromising functionality. Advanced nuclear reactors, including Small Modular Reactors (SMRs) and microreactors, aim to achieve autonomous and remote operations, reducing costs and enhancing competitiveness. However, these advancements expand the attack surface for cyberattacks, particularly in autonomous and remote operation scenarios. Cyberattacks can exploit vulnerabilities to manipulate physical processes, causing shutdowns, asset damage, or public harm. Such attacks begin with passive reconnaissance, where adversaries intercept communications or observe behaviors to gather information, which is then leveraged to execute cyber-physical attacks by injecting malicious commands. Nuclear power must adopt cybersecurity protection measures to secure the integrity and availability of their digital control systems. This paper demonstrates the application of FHE to secure operations by enabling encrypted processing of sensitive signals and parameters -- ensuring privacy without exposing data. FHE supports secure mathematical operations on encrypted data without requiring decryption. Using a hardware-in-the-loop (HIL) approach, this paper implements an FHE-integrated controller on a BeagleBone Black (BBB) controlling a simulation of the Small Modular Advanced High Temperature Reactor (SmAHTR). By doing so, the encrypted controller protects the integrity of critical set points and control signals during transmission and processing. Thus, FHE-integrated controllers enhance secure operations of advanced nuclear reactors while maintaining functionality.

control systems↗
DOE OSTI2024

High-throughput genetics enables identification of nutrient utilization and accessory energy metabolism genes in a model methanogen

Archaea are widespread in the environment and play fundamental roles in diverse ecosystems; however, characterization of their unique biology requires advanced tools. This is particularly challenging when characterizing gene function. Here, we generate randomly barcoded transposon libraries in the model methanogenic archaeon Methanococcus maripaludis and use high-throughput growth methods to conduct fitness assays (RB-TnSeq) across over 100 unique growth conditions. Using our approach, we identified new genes involved in nutrient utilization and response to oxidative stress. We identified novel genes for the usage of diverse nitrogen sources in M. maripaludis including a putative regulator of alanine deamination and molybdate transporters important for nitrogen fixation. Furthermore, leveraging the fitness data, we inferred that M. maripaludis can utilize additional nitrogen sources including $\tiny{L}$-glutamine, $\tiny{D}$-glucuronamide, and adenosine. Under autotrophic growth conditions, we identified a gene encoding a domain of unknown function (DUF166) that is important for fitness and hypothesize that it has an accessory role in carbon dioxide assimilation. Finally, comparing fitness costs of oxygen versus sulfite stress, we identified a previously uncharacterized class of dissimilatory sulfite reductase-like proteins (Dsr-LP; group IIId) that is important during growth in the presence of sulfite. When overexpressed, Dsr-LP conferred sulfite resistance and enabled use of sulfite as the sole sulfur source. The high-throughput approach employed here allowed for generation of a large-scale data set that can be used as a resource to further understand gene function and metabolism in the archaeal domain.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Solid Oxide Cell Manufacturing Cost Tool

NETL has developed a detailed solid oxide cell (SOC) cost modeling tool. The tool provides a fully customizable spreadsheet that aids in cost estimation of large volume manufacturing of SOC (both fuel cell and electrolysis cell). The tool can be used to assist commercial developers and other end users to evaluate the costs of producing cells and stack components as SOC systems increase in scale. The manufacturing materials, labor, component costs, and specific equipment can be tailored to accurately reflect the desired manufacturing process. The tool is paired with a detailed user manual to assist in the understanding of the full functionality. This document is the spreadsheet tool. To access the associated user guide, please visit: <a href="https://netl.doe.gov/energy-analysis/details?id=f5823130-0d20-4a04-9c53-341682e58ac1" rel="noopener noreferrer">Energy Analysis | netl.doe.gov</a>

cost modeling↗
DOE OSTI2020

Cooperative Research and Development Agreement between National Energy Technology Laboratory and Birla Carbon U.S.A., Inc. [Abstract]

NETL researchers will be collaborating with Birla Carbon, U.S.A., Inc. (Birla Carbon) to leverage NETL's patented catalyst work on methane decomposition to commercialize an efficient and a cost-effective process for production of hydrogen and valuable carbon from natural gas. This catalytic methane decomposition is a promising technology, however, finding an appropriate catalyst to perform all the functions has been a challenge. NETL has developed a suitable catalyst for the process which has shown promising results. The heat required for this endothermic catalytic decomposition process is provided by combustion of a portion of the hydrogen. This process is cost effective, energy efficient, and environmentally friendly with minimal CO 2 mitigation costs. The objectives of the project include evaluation of the effect of feed gas composition on the performance, reduction of the catalyst production cost, and conduction of an in-depth large-scale cost analysis for a process suitable for applications at Birla Carbon.

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

A light weight regularization for wave function parameter gradients in quantum Monte Carlo

The parameter derivative of the expectation value of the energy, ∂E/∂p, is a key ingredient in variational Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve estimate of this derivative suffers from an infinite variance, which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator, which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias, which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of ∂E/∂p for a multi-Slater–Jastrow trial wave function on the LiH molecule and in the optimization of a multi-Slater–Jastrow trial wave function on the CuO molecule. This regularized estimator is a simple and efficient estimator of ∂E/∂p for VMC optimization techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Abstract We present a way to dramatically accelerate Gaussian process models for interatomic force fields based on many-body kernels by mapping both forces and uncertainties onto functions of low-dimensional features. This allows for automated active learning of models combining near-quantum accuracy, built-in uncertainty, and constant cost of evaluation that is comparable to classical analytical models, capable of simulating millions of atoms. Using this approach, we perform large-scale molecular dynamics simulations of the stability of the stanene monolayer. We discover an unusual phase transformation mechanism of 2D stanene, where ripples lead to nucleation of bilayer defects, densification into a disordered multilayer structure, followed by formation of bulk liquid at high temperature or nucleation and growth of the 3D bcc crystal at low temperature. The presented method opens possibilities for rapid development of fast accurate uncertainty-aware models for simulating long-time large-scale dynamics of complex materials.

Chemistry↗
DOE OSTI2025

Defect formation and microstructure tuning via proton irradiation to control electrochemical and phase reversibility in layered battery materials

The reversibility of phase transformation influences the functionality of electrode materials in batteries. In many battery materials, nanosized grains favor phase reversibility but at the cost of cyclability due to aggravated side reactions with the electrolyte. Here, in this study, we present a novel approach to enhance the phase transformation reversibility of layered oxide cathodes, exemplified by Na 2/3 Fe 1/2 Mn 1/2 O 2 through proton irradiation. In addition to forming defects, proton irradiation at sufficiently high doses can subdivide single grains into multiple nanodomains without physically rupturing them. Hence, the single grains of the material assume a pseudo-secondary particle nature without reducing the overall grain size. Preserving the grain size is advantageous, as it reduces side reactions, which is not possible with conventional grain size reduction methods. While chemical transformations and defect formation induced through proton irradiation can influence the stability of battery materials, it is expected that structural reorganization due to cycling-induced phase transformation will be contained within these nanodomains. Such confinement of phase transformation is potentially responsible for enhancing the reversibility of layered oxide materials in our study. Thus, our study suggests that grain subdivision could become an effective microstructure tuning strategy for managing electrochemical cycling-induced phase changes in battery electrodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Invariant Molecular Representations for Heterogeneous Catalysis

Catalyst screening is a critical step in the discovery and development of heterogeneous catalysts, which are vital for a wide range of chemical processes. In recent years, computational catalyst screening, primarily through density functional theory (DFT), has gained significant attention as a method for identifying promising catalysts. However, the computation of adsorption energies for all likely chemical intermediates present in complex surface chemistries is computationally intensive and costly due to the expensive nature of these calculations and the intrinsic idiosyncrasies of the methods or data sets used. This study introduces a novel machine learning (ML) method to learn adsorption energies from multiple DFT functionals by using invariant molecular representations (IMRs). To do this, we first extract molecular fingerprints for the reaction intermediates and later use a Siamese-neural-network-based training strategy to learn invariant molecular representations or the IMR across all available functionals. Our Siamese network-based representations demonstrate superior performance in predicting adsorption energies compared with other molecular representations. Notably, when considering mean absolute values of adsorption energies as 0.43 eV (PBE-D3), 0.46 eV (BEEF-vdW), 0.81 eV (RPBE), and 0.37 eV (scan+rVV10), our IMR method has achieved the lowest mean absolute errors (MAEs) of 0.18 0.10, 0.16, and 0.18 eV, respectively. These results emphasize the superior predictive capacity of our Siamese network-based representations. The empirical findings in this study illuminate the efficacy, robustness, and dependability of our proposed ML paradigm in predicting adsorption energies, specifically for propane dehydrogenation on a platinum catalyst surface.

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

Deep Learning Potential For Investigation of Chemical Conversion of CO2 at Air-Reactive Deep Eutectic Solvent Interface

Capturing CO2 from air is urgent to battle against the climate change. Traditional CO2 capture by alkaline solution or aqueous amine has high regeneration cost, equipment corrosion, and amine leakage. Deep eutectic solvents (DESs) are promising alternatives with low corrosion, non-toxicity, and biodegradable nature. Functionalized DESs can capture and convert CO2 at low partial pressures. Developing reactive DESs for efficient CO2 capture needs a foundational understanding of the reaction-transport coupling of CO2 at the air-DES interface. In this study, we develop a deep learning potential (DLP) to investigate the absorption and reaction of CO2 at the air-DES interface (DES: 1-ethyl-3-methylimidazolium 2-cyanopyrrolide ([Emim][2-CNpyr]) and ethylene glycol (EG)). The DLP model allows us to simulate chemical reactions at a lower computational cost with ab initio accuracy. By analyzing reaction free energy surfaces, molecular interactions among CO2, DES components and reaction products, and transport of CO2 through the interface, we identify three main reaction pathways of CO2 at the interface: forming carboxylate with [Emim], carbamate with [2-CNpyr], and carbonate with EG. The mechanistic understanding of CO2 chemisorption at the interface will facilitate the development of novel direct air capture technology using DESs.

Nguyen, Manh Tien↗
DOE OSTI2022

Modified Calcium Manganites for Thermochemical Energy Storage Applications

CaAl 0.2 Mn 0.8 O 3-δ (CAM28) and CaTi 0.2 Mn 0.8 O 3-δ (CTM28) are perovskite metal oxides developed for high-temperature thermochemical energy storage (TCES) applications, e.g., in support of air Brayton power generation. Previous reports for these compounds focus on the equilibrium non-stoichiometry (δ) as a function of temperature and oxygen partial pressure (pO 2 ) and the endotherm (or exotherm) accompanying changes in δ resulting from thermal reduction (or re-oxidation). Herein, we report results for elemental substitution and doping (Al, Co, Fe, La, Sr, Ti, Y, Zn, and Zr) of calcium manganites (CM) that establish the preference for CAM28 and CTM28. Techniques employed include conventional (screening and equilibrium) and ballistically heated multi-cycle thermogravimetric analysis (TGA), conventional and high temperature ( in-situ ) X-ray diffraction (XRD), and differential scanning calorimetry (DSC). Forward-looking results for A-site Y-doped materials, e.g., Ca 0.9 Y 0.1 MnO 3-δ (CYM910), establish a route to increasing the reduction enthalpy relative to CAM28 and CTM28, albeit at the expense of increased reduction temperatures and raw materials costs. A thermodynamic model presented for CAM28, but extendable to related materials, provides values for the reaction enthalpy and extent of reduction as a function of temperature and oxygen partial pressure for use in design efforts. Taken as a whole, the results support the choice of Al-doped CaMnO 3-δ as a low-cost material for TCES in a high temperature air Brayton application, but point the way to achieving higher stored energy densities that could lead to overall cost savings.

Miller, James E.↗
DOE OSTI2025

Calibration of RAFM Micromechanical Model for Creep Using Bayesian Optimization for Functional Output

A Bayesian optimization procedure is presented for calibrating a multimechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on Fe(8–9)%Cr are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600°C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework (Huang et al., 2021, “Bayesian optimization of functional output in inverse problems,” Optim. Eng., 22, pp. 2553–2574). Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. In conclusion, validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel.

42 ENGINEERING↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Commercial building HVAC demand flexibility with model predictive control: Field demonstration and literature insights

Model Predictive Control (MPC) for building Heating Ventilation and Air Conditioning (HVAC) systems is beginning to gain traction in the market, with a few controls companies incorporating it into their product offerings. However, it remains difficult to assess whether the energy cost savings are enough to justify the cost of MPC implementation for a particular building, given the limited number of reported demonstrations. For small commercial and residential buildings with relatively uniform systems, standardized approaches can help lower implementation costs. In contrast, for large buildings or district systems, the potential magnitude of cost savings could justify more customized solutions. Estimating the cost-effectiveness of MPC becomes more challenging for medium and large commercial buildings, where a one-size-fits-all solution may not be suitable, and the potential energy cost savings may be insufficient to justify a customized solution. To make MPC technology more appealing, incorporating additional value streams beyond energy efficiency alone can significantly increase its attractiveness. One such revenue stream is demand flexibility, in response to dynamic electricity prices, where MPC can leverage the thermal mass of the building to shift the load and support the grid. Building on an extensive literature review of MPC field studies focused on cost savings and demand flexibility, this paper presents the results of implementing MPC control in a large office building HVAC system in Berkeley, CA. Four different dynamic electricity price profiles were integrated into the MPC objective function to shift building demand while maintaining comfort, and field testing was performed with each price profile across four seasons. The results show potential for 40–65 % demand decrease percentage and up to 61 % annual cost savings compared to the existing rule-based control strategy, under the tested dynamic price scenarios. This paper also presents a sensitivity analysis on the cost savings with respect to the price profile variability, discusses the implementation effort for the price-responsive MPC, and compares the cost savings found in this study to those found in literature on the basis of dynamic price variability, or so-called Electricity Price Relative Standard Deviation.

Zanetti, Ettore↗
DOE OSTI2022

Weight function procedure for reduced order fracture analysis of arbitrary flaws in cylindrical pressure vessels

Fracture mechanics calculations using the finite element method can be computationally expensive, which makes them challenging to use in engineering evaluations of flaws in pressure vessels. The weight function (WF) technique is a reduced order fracture modeling method that is widely used to greatly reduce these computational costs. Its computational efficiency is essential for use in probabilistic fracture mechanics evaluations of embrittled nuclear reactor pressure vessels (RPVs), due to the large number of sampled flaws that must be evaluated. Although the WF technique is general, it is typically used only for axis-aligned flaws. Recent discoveries of off-axis flaws in operating nuclear reactors have necessitated detailed simulations of such flaws and the interactions between neighboring flaws. This paper presents a generalized WF approach applicable for analyzing arbitrary flaw geometries in thick-walled cylindrical vessels, including surface-breaking and subsurface flaws, which can either axis-aligned or off-axis, and can account for interactions with other flaws. Herein, this approach is demonstrated on representative simulations of multiple flaw geometries in a RPV subjected to transient loading conditions. In all cases, the WF technique gives good comparison with benchmark results from direct simulation with greatly reduced computational effort.

42 ENGINEERING↗
DOE OSTI2022

Multi-objective optimization of root phenotypes for nutrient capture using evolutionary algorithms

Root phenotypes are avenues to the development of crop cultivars with improved nutrient capture, which is an important goal for global agriculture. The fitness landscape of root phenotypes is highly complex and multidimensional. It is difficult to predict which combinations of traits (phene states) will create the best performing integrated phenotypes in various environments. Brute force methods to map the fitness landscape by simulating millions of phenotypes in multiple environments are computationally challenging. Evolutionary optimization algorithms may provide more efficient avenues to explore high dimensional domains such as the root phenotypic space. We coupled the three-dimensional functional–structural plant model, SimRoot, to the Borg Multi-Objective Evolutionary Algorithm (MOEA) and the evolutionary search over several generations facilitated the identification of optimal root phenotypes balancing trade-offs across nutrient uptake, biomass accumulation, and root carbon costs in environments varying in nutrient availability. Our results show that several combinations of root phenes generate optimal integrated phenotypes where performance in one objective comes at the cost of reduced performance in one or more of the remaining objectives, and such combinations differed for mobile and non-mobile nutrients and for maize (a monocot) and bean (a dicot). Functional–structural plant models can be used with multi-objective optimization to identify optimal root phenotypes under various environments, including future climate scenarios, which will be useful in developing the more resilient, efficient crops urgently needed in global agriculture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

A Multiple Filter Based Neural Network Approach to the Extrapolation of Adsorption Energies on Metal Surfaces for Catalysis Applications

Computational catalyst discovery involves the development of microkinetic reactor models based on estimated parameters determined from density functional theory (DFT). For complex surface chemistries, the number of reaction intermediates can be very large, and the cost of calculating the adsorption energies by DFT for all surface intermediates even for one active site model can become prohibitive. In this paper, we have identified appropriate descriptors and machine learning models that can be used to predict a significant part of these adsorption energies given data on the rest of them. Furthermore, our investigations also included the case when the species data used to train the predictive model are of different size relative to the species the model tries to predict—this is an extrapolation in the data space which is typically difficult with regular machine learning models. Due to the relative size of the available data sets, we have attempted to extrapolate from the larger species to the smaller ones in the current work. Here, we have developed a neural network based predictive model that combines an established additive atomic contribution based model with the concepts of a convolutional neural network that, when extrapolating, achieves a statistically significant improvement over the previous models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗