Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive thermodynamics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗
DOE OSTI2022

What Is the Smallest Zeolite That Could Be Synthesized?**

Zeolites with a few unit cells are promising as catalyst and adsorbents. The quest to synthesize the smallest zeolites has recently resulted in 4 to 8 nm nanozeolites, about 2 to 4 unit cells. These findings pose the question of what is the smallest zeolite that could be obtained by hydrothermal synthesis. We address this question using molecular simulations and thermodynamic analysis. The simulations predict that amorphous precursors as small as 4 nm can crystallize zeolites, in agreement with the experiments. We find that interfacial forces dominate the structure of smaller particles, resulting in size-dependent compact isomers that have ring and pore distributions different from open framework zeolites. The instability of zeolites smaller than 3±0.5 nm precludes a classical mechanism of nucleation from solution or through assembly of small nanoslabs.

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

Leveraging slow $\mathrm{DOTA}$ f-element complexation kinetics to enable separations by kinetic design

The design of metal-concerned solvent extraction systems frequently leverages thermodynamically derived differences in selectivity. An alternative approach, leveraging kinetic control, has been considered much less seriously. Our recent manuscript describing DOTA (1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) complexation kinetics across the lanthanide series shows the observed rate constant steadily increases across the series, with some non-monotonic behavior observed at terbium and thulium. This contrasts the thermodynamic stability constant trend, where lanthanide-DOTA stability constants initially increase and then plateau as a function of ionic radii after samarium. To leverage the kinetic differences of DOTA with the lanthanides across the series, kinetically based separations must be utilized. Since DOTA has very slow complexation kinetics, a separations system must expedite DOTA-metal complexation to allow a separation approximating practical application. Here in this report, a DOTA-based solvent extraction system, where DOTA is the aqueous holdback reagent and bis-2,4,4-trimethylpentylphosphinic acid (Cyanex 272) is the organic phase extractant, is demonstrated and compares the separations chemistry of Nd, Eu and Am. The slowness of DOTA complexation was addressed by heating the system. Results showed, in general, separation between metals is better during early phase contact, and diminishes under longer contact times. Under all conditions, separations are better than would be predicted based on a thermodynamic basis. This report suggests that while slowly complexing ligands classically used for biological applications may not be appropriate for thermodynamically designed metal separations, their use for kinetically based systems may be appropriate and enable a new design basis for f-element separations.

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

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Subtle metastability of the layered magnetic topological insulator MnBi 2 Te 4 from weak interactions

Abstract Layered quantum materials can host interesting properties, including magnetic and topological, for which enormous computational predictions have been done. Their thermodynamic stability is much less visited computationally, which however determines the existence of materials and can be used to guide experimental synthesis. MnBi 2 Te 4 is one of such layered quantum materials that was predicted to be an intrinsic antiferromagnetic topological insulator, and later experimentally realized but in a thermodynamically metastable state. Here, using a combined first-principles-based approach that considers lattice, charge, and spin degrees of freedom, we investigate the metastability of MnBi 2 Te 4 by calculating the Helmholtz free energy for the reaction Bi 2 Te 3 + MnTe → MnBi 2 Te 4 . We identify a temperature range (~500–873 K) in which the compound is stable with respect to the competing binary phases, consistent with experimental observation. We validate the predictions by comparing the calculated specific heats contributed from different degrees of freedom with experimental results. Our findings indicate that the degrees of freedom responsible for the van der Waals interaction, lattice vibration, magnetic coupling, and nontrivial band topology in MnBi 2 Te 4 not only enable emergent phenomena but also play a crucial role in determining its thermodynamic stability. This conclusion lays the foundation for the future computational material synthesis of novel layered systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

What Is the Smallest Zeolite That Could Be Synthesized?**

Abstract Zeolites with a few unit cells are promising as catalyst and adsorbents. The quest to synthesize the smallest zeolites has recently resulted in 4 to 8 nm nanozeolites, about 2 to 4 unit cells. These findings pose the question of what is the smallest zeolite that could be obtained by hydrothermal synthesis. Here we address this question using molecular simulations and thermodynamic analysis. The simulations predict that amorphous precursors as small as 4 nm can crystallize zeolites, in agreement with the experiments. We find that interfacial forces dominate the structure of smaller particles, resulting in size‐dependent compact isomers that have ring and pore distributions different from open framework zeolites. The instability of zeolites smaller than 3±0.5 nm precludes a classical mechanism of nucleation from solution or through assembly of small nanoslabs.

Dhabal, Debdas↗
DOE OSTI2020

Polymorphic Metastability in Colloidal Semiconductor Nanocrystals

Metastable polymorphs of inorganic solids often possess material properties not present in the corresponding thermodynamic polymorphs, making them targets for the development of new functional materials. In contrast with isolating metastable bulk materials, syntheses of metastable polymorphs on the nanoscale are aided by fast non-equilibrium reaction kinetics and the favorable thermodynamic influence of surface energies, giving rise to greater ease of access to metastable high-temperature polymorphs and, in some cases, new polymorphs that do not exist in the bulk. The syntheses of metastable semiconductor nanocrystals are of interest for their potentially unique optoelectronic and physicochemical properties. However, in many material systems, synthesizing nanocrystalline products away from thermodynamic equilibrium in a predictable manner remains an outstanding challenge. This review outlines direct synthetic methodologies that have been developed to enable control over the nucleation and growth of metastable polymorphs of semiconductor nanocrystals by tailoring reaction conditions, precursor kinetics, ligand and surface effects, and other synthetic levers. The case studies reviewed herein expound on the direct syntheses of metastable ZnSe, Cu 2 SnSe 3 , CuInSe 2 , Ag 2 Se, and AgInSe 2 nanocrystals, and although there remain numerous examples of metastable nanocrystal syntheses outside of these metal chalcogenide systems, the concepts discussed are of general utility to the field of metastable nanocrystal syntheses as a whole. Explicit examples in which new functional properties are afforded by metastable polymorphs of the aforementioned material systems are presented within the context of applications for solar cells, photonics, and optical sensing. Finally, the factors that affect the kinetic persistence of metastable nanocrystalline polymorphs are discussed at length for these material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗
DOE OSTI2021

Investigation of $\gamma$-irradiated polyvinylidene fluoride and its acute toxicity

Polyvinylidene fluoride (PVDF) is an industrial polymer with many applications. The EPR spectrum of polyvinylidene fluoride (PVDF) irradiated by γ-rays at low temperature (-196°C) provides evidence for the presence of radicals formed by loss of a fluorine atom from PVDF, and polyenyl radicals. Polyenyl radicals dominate in the polymer irradiated at room temperature. Peroxide radicals formed as a result of oxidation of primary free radicals by atmospheric oxygen are observed in the EPR spectrum of a polymer irradiated in air, or when air is introduced to a polymer irradiated in a vacuum. Carbonization and oxidation of γ-irradiated PVDF were observed in the XPS spectrum. A model for polymer degradation based on correlated G3(MP2) molecular orbital theory electronic structure calculations is described in this report. The presence of HF in the products and a low yield of products with the composition CxFyHz indicate the absence of noticeable thermal destruction of the main polymer chain upon heating to 220°C of both the initial and γ-irradiated PVDF. Although any HF released during the PVDF dehydrofluorination would be toxic, a study of the acute toxicity of γ-irradiated PVDF showed that the polymer does not exhibit toxicity after a single intragastric route of administration to mice.

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

Equilibrium spin polarization arising from chirality

Chirality-induced spin selectivity (CISS) describes how chiral molecules and materials generate spin polarization even at thermal equilibrium. This observation has challenged established principles of microscopic reversibility and Onsager reciprocity. We resolve this paradox by formulating a pseudo-Hermitian quantum framework that separates thermodynamic equilibrium from time-reversal symmetry. Within this approach, structural chirality and electron correlations, irrespective of their microscopic origin, are sufficient to produce CISS observables. Chirality enters through a non-local metric η that couples spin and spatial motion, leading to real spectra, unitary evolution, and thermodynamic consistency. The framework predicts a chirality-induced spin magnetic ordering characterized by a spin-displacement order, which reconciles equilibrium spin polarization with detailed balance and explains the persistence of CISS in materials composed of light elements. We derive generalized Onsager-Casimir relations that respect the observed CISS symmetry, i.e., parity $(\mathscr{P})$-odd and time-reversal $(\mathscr{T})$-odd, but exhibiting $\mathscr{PT}$-even symmetry. This approach establishes a coherent foundation for equilibrium CISS and provides a route to link chemical chirality with measurable spin-to-charge conversion effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effect of initial γ -irradiation on infrared laser ablation of poly(vinyl alcohol) studied by infrared spectroscopy

We report the effect of initial γ-irradiation of poly(vinyl alcohol) (PVA) on the kinetics of polymer ablation under CO 2 laser irradiation was studied using infra-red spectroscopy. The rate of laser ablation rapidly decreases with increasing time of laser ablation up to 10 s. Initial γ-irradiation of a dose up to 100 kGy reduces the rate of the laser ablation, just like the effect of laser exposure. The number of OH and CH groups decrease in PVA both with laser irradiation and with γ-irradiation. The number of OH groups decreases more under laser irradiation. Unsaturated bonds are formed in γ–irradiated PVA on the surface of the crater formed by laser ablation of the polymer that has undergone radiolysis. Increasing the dose of initial γ-irradiation from 100 to 2300 kGy leads to an increase in the intensity of the IR bands of the unsaturated bonds. The number of IR bands in the spectra of the laser ablation crater of PVA does not change with the dose of initial γ-irradiation. A dose of radiolysis of 3500 kGy significantly changes the IR spectra of both the ablation crater and the coating formed from the ablation products. Many bands characteristic for PVA are no longer present and the intensities of bands corresponding to unsaturated bonds are increased. A computational model for loss of water and ketone formation was developed for PVA and PVA including a diol using G3(MP2). The highest activation energy barrier step in the decomposition reaction was formation of the enol species and water was shown to catalyze all decomposition steps.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Joint influence of γ–irradiation and high temperature shear grinding on the IR spectra and surface–energy properties of polyethylene

The optical and surface-energy properties of polyethylene (PE) subjected to γ-irradiation and post-radiation, high temperature, shear grinding are presented. Comparison of the IR spectra of the initial and γ-irradiated PE shows that oxygen-containing groups and unsaturated bonds in the polymer macromolecules are formed and accumulate during their radiolysis in air, increasing with γ-irradiation dose. A redistribution of the products of radiation-chemical conversion on the surface of the polymer granules throughout the entire volume of the powder occurs after high-temperature shear milling to a powder as observed in the IR spectrum. Polar functional groups due to oxidation reactions are observed on the surface layer of the sample after post-radiation grinding. Furthermore, the initial polymer surface is weakly basic with an acidity parameter of -0.34. The value of the acidity parameter of the sample made of ground irradiated PE has an acidic surface with an acidity parameter of 0.65. γ irradiation of the secondary processing of PE waste with high-temperature shear grinding is a promising approach to obtain a secondary PE powder with specified parameters and different degrees of functionalization.

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

Predicting short-range order evolution in WTaCrVHf refractory high-entropy alloys

Short-range order (SRO) in multicomponent concentrated alloys affects their mechanical response. Hence, is paramount to understand how composition modifies the chemical ordering in the system to design materials with optimal properties. Here, in this work, we present a methodology to predict the SRO and thermodynamic properties in chemically complex systems and apply it to the WTaCrVHf quinary alloy. We observe that the addition of Hf significantly modifies the SRO, mainly at intermediate to low temperatures, matching experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prediction of Solution Behavior via Calorimetric Measurements Allows for Detailed Elucidation of Polyoxometalate Transformation

The solution behavior of a polyoxometalate cluster, LiNa-U 24 Pp 12 (Li 24 Na 24 [(UO 2 O 2 ) 24 (P 2 O 7 ) 12 ]) that consists of 24 uranyl ions, peroxide groups, and 12 pyrophosphate linkers, was successfully predicted based on new thermodynamic results using a calorimetric method recently described for uranyl peroxide nanoclusters (UPCs), molybdenum blues, and molybdenum browns. The breakdown of LiNa-U 24 Pp 12 and formation of U 24 (Li 24 [UO 2 O 2 OH] 24 ) was monitored in situ via Raman spectroscopy using a custom heating apparatus. A combination of analytical techniques confirmed the simultaneous existence of U 24 Pp 12 and U 24 midway through the conversion process and U 24 as the single end product. The application of a molecular weight filter resulted in a complete and successful separation of UPCs from solution and, in conjunction with DOSY results, confirmed the presence of large intermediate cluster building blocks.

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

CHARMM-GUI Nanomaterial Modeler for Modeling and Simulation of Nanomaterial Systems

Molecular modeling and simulation are invaluable tools for nanoscience that predict mechanical, physicochemical, and thermodynamic properties of nanomaterials and provide molecular-level insight into underlying mechanisms. However, building nanomaterial-containing systems remains challenging due to the lack of reliable and integrated cyberinfrastructures. Here we present Nanomaterial Modeler in CHARMM-GUI, a web-based cyberinfrastructure that provides an automated process to generate various nanomaterial models, associated topologies, and configuration files to perform state-of- the-art molecular dynamics simulations using most simulation packages. The nanomaterial models are based on the interface force field, one of the most reliable force fields (FFs). The transferability of nanomaterial models among the simulation programs was assessed by single-point energy calculations, which yielded 0.01% relative absolute energy differences for various surface models and equilibrium nanoparticle shapes. Three widely used Lennard-Jones (LJ) cutoff methods are employed to evaluate the compatibility of nanomaterial models with respect to conventional biomolecular FFs: simple truncation at r = 12 Å (12 cutoff), force-based switching over 10 to 12 Å (10-12 fsw), and LJ particle mesh Ewald with no cutoff (LJPME). The FF parameters with these LJ cutoff methods are extensively validated by reproducing structural, interfacial, and mechanical properties. As a result, we find that the computed density and surface energies are in good agreement with reported experimental results, although the simulation results increase in the following order: 10-12 fsw <12 cutoff < LJPME. Nanomaterials in which LJ interactions are a major component show relatively higher deviations (up to 4% in density and 8% in surface energy differences) compared with the experiment. Nanomaterial Modeler's capability is also demonstrated by generating complex systems of nanomaterial-biomolecule and nanomaterial-polymer interfaces with a combination of existing CHARMM-GUI modules. We hope that Nanomaterial Modeler can be used to carry out innovative nanomaterial modeling and simulations to acquire insight into the structure, dynamics, and underlying mechanisms of complex nanomaterial-containing systems.

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

Pathways to High-Performance Salt Hydrate Thermochemical Energy Storage Materials and Systems

Thermochemical materials (TCMs) based on salt hydrates are promising for thermal energy storage as they combine high energy densities with low reaction temperatures. However, their adoption is hindered by poor structural integrity and degradation under hygrothermal cycling. Storage performance is governed not only by the chemical reaction, but also by the coupled thermo-chemo-mechanical behavior that evolves with cycling. Understanding and controlling this coupling across length scales (material-to-reactor) is necessary to improve TCM stability and lifetime. In this perspective, we discuss the shortcomings of current characterization approaches and emphasize the need for measuring transport properties and structural transformations using in situ techniques that capture the dynamic evolution of these materials. We also outline opportunities for multiscale modeling frameworks that link thermodynamics and mechanics, enabling predictive evaluation of composite architectures designed for cycling stability. We conclude by identifying research questions that must be addressed to transform TCMs into viable energy storage technologies.

Composites↗
DOE OSTI2025

Multimodal In Situ Characterization Uncovers Unexpected Stability of a Cobalt Electrocatalyst for Acidic Sustainable Energy Technologies

An accelerated development of durable and affordable sustainable energy technologies is often hindered by a limited understanding of how non-precious materials within these systems degrade. In acidic proton exchange membrane fuel cells and water electrolyzers, metallic cobalt (Co) is considered an unstable component that is often combined with precious metals or other stabilizers. To understand the mechanisms behind Co instability, we employ an experimental platform that quantifies dissolution with on-line inductively coupled plasma mass spectrometry and product formation with electrochemical mass spectrometry during electrochemical testing, along with ex- situ characterization. Under varied conditions (electrocatalysis, time, gas-type saturation, and ion concentration), windows of Co stability are observed that are different than predicted with classical chemical thermodynamics, suggesting new stabilization and degradation mechanisms than previously understood. Notably, Co is active for the hydrogen evolution reaction (HER), with prolonged stability that is ~300 mV different than thermodynamically projected. Additionally, in an oxygenated environment, Co concurrently performs the HER and oxygen reduction reaction (ORR) yet undergoes different morphology changes and dissolution mechanisms. Interestingly, in the absence of electrocatalysis, there is a 22x decrease in dissolution in an oxygen-free environment, proposing a route to decrease Co losses during device shutdown protocols. Lastly, under more extreme operating conditions, Co becomes stable after a substantial amount of dissolution, suggesting that high concentrations of Co 2+ ions in the microenvironment induce the formation of a stable CoHO 2 surface. Altogether, these results can be leveraged to improve the design and development of more robust and cost-effective sustainable energy technologies, as well as promote strategic strategies for prolonged material utilization.

Dissolution↗
DOE OSTI2022

Diamond and methane formation from the chemical decomposition of polyethylene at high pressures and temperatures

Abstract Polyethylene (C 2 H 4 ) n was compressed to pressures between 10 and 30 GPa in a diamond anvil cell (DAC) and laser heated above 2500 K for approximately one second. This resulted in the chemical decomposition of the polymer into carbon and hydrocarbon reaction products. After quenching to ambient temperature, the decomposition products were measured in the DAC at pressures ranging from ambient to 29 GPa using a combination of x-ray diffraction (XRD) and small angle x-ray scattering (SAXS). XRD identified cubic diamond and methane as the predominant product species with their pressure–volume relationships exhibiting strong correlations to the diamond and methane equations of state. Length scales associated with the diamond products, obtained from SAXS measurements, indicate the formation of nanodiamonds with a radius of gyration between 12 and 35 nm consistent with 32–90 nm diameter spherical particles. These results are in good agreement with the predicted product composition under thermodynamic and chemical equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗