Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid structure simulation”

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

Control of Permanent Porosity in Type 3 Porous Liquids via Solvent Clustering

Porous liquids (PLs) are an exciting new class of materials for carbon capture due to their high gas adsorption capacity and ease of industrial implementation. They are composed of sorbent particles suspended in a nonadsorbed solvent, forming a liquid with permanent porosity. While PLs have a vast number of potential compositions based on the number of solvents and sorbent materials available, most of the research has been focused on the selection of the sorbent rather than the solvent. Therefore, PL design criteria on the supramolecular structures of the solvent are explored to create a fundamental understanding of how the solvent enables PL formation for rapid discovery of new PL compositions. Atomistic molecular dynamics simulation of eight solvents with a range of molecular sizes, shapes, and intramolecular bonding was performed, identifying that the shape and size of molecular clusters formed in the solvent are the driving predictor of PL formation rather than the size of the individual solvent molecule. The results demonstrate a significant departure from common approaches to PL formation based on the steric exclusion of solvent molecules from the sorbent via the size of the pore aperture. A modeling and experimental validation study further supports these findings. In conclusion, through this computational material design study, a previously unexplored mechanism in PL formation, solvent–solvent clustering, is identified as a critical factor for the accelerated discovery of liquid phase carbon capture materials.

Carbon capture↗

Fluid management technology: Liquid slosh dynamics and control

Flight experiments were defined for the Cryogenic On-Orbit Liquid Depot Storage, Acquisition and Transfer Satellite (COLD-SAT) test bed satellite and the Shuttle middeck to help establish the influence of the gravitational environment on liquid slosh dynamics and control. Several analytical and experimental studies were also conducted to support the experiments and to help understand the anticipated results. Both FLOW-3D and NASA-VOF3D computer codes were utilized to simulate low Bond number, small amplitude sloshing, for which the motions are dominated by surface forces; it was found that neither code provided a satisfactory simulation. Thus, a new analysis of low Bond number sloshing was formulated, using an integral minimization technique that will allow the assumptions made about surface physics phenomena to be modified easily when better knowledge becomes available from flight experiments. Several examples were computed by the innovative use of a finite-element structural code. An existing spherical-pendulum analogy of nonlinear, rotary sloshing was also modified for easier use and extended to low-gravity conditions. Laboratory experiments were conducted to determine the requirements for liquid-vapor interface sensors as a method of resolving liquid surface motions in flight experiments. The feasibility of measuring the small slosh forces anticipated in flight experiments was also investigated.

Dodge, Franklin T.↗

Configurational-Bias Monte Carlo Back-Mapping Algorithm for Efficient and Rapid Conversion of Coarse-Grained Water Structures into Atomistic Models

Coarse-grained molecular dynamics (MD) simulations represent a powerful approach to simulate longer time scale and larger length scale phenomena than those accessible to all-atom models. The gain in efficiency, however, comes at the cost of atomistic details. The reverse transformation, also known as back-mapping, of coarse grained beads into their atomistic constituents represents a major challenge. Most existing approaches are limited to specific molecules or specific force-fields and often rely on running a long time atomistic MD of the back-mapped configuration to arrive at an optimal solution. Such approaches are problematic when dealing with systems with high diffusion barriers. Here, we introduce a new extension of the configurational-bias-Monte-Carlo (CBMC) algorithm, which we term the crystalline-configurational-bias-Monte-Carlo (C-CBMC) algortihm, that allows rapid and efficient conversion of a coarse-grained model back into its atomistic representation. Although the method is generic, we use a coarse-grained water model as a representative example and demonstrate the back-mapping or reverse transformation for model systems ranging from the ice-liquid water interface to amorphous and crystalline ice configurations. A series of simulations using the TIP4P/Ice model are performed to compare the new CBMC method to several other standard Monte Carlo and Molecular Dynamics based back-mapping techniques. In all the cases, the C-CBMC algorithm is able to find optimal hydrogen bonded configuration many thousand evaluations/steps sooner than the other methods compared within this paper. For crystalline ice structures such as a hexagonal, cubic, and cubic-hexagonal stacking disorder structures, the C-CBMC was able to find structures that were between 0.05 and 0.1 eV/water molecule lower in energy than the ground state energies predicted by the other methods. Detailed analysis of the atomistic structures show a significantly better global hydrogen positioning when contrasted with the existing simpler back-mapping methods. The error in the RDFs of backmapped relative to reference configuration for the CCBMC, MD, and MC were found to be 6.9, 8.7, and 12.9 respectively for the hexagonal system. For the cubic system the relative error of the RDFs for the CCBMC, MD, and MC were found to be 18.2, 34.6, and 39.0 respectively. Finally, our results demonstrate the efficiency and efficacy of our new back-mapping approach, especially for crystalline systems where simple force-field based relaxations have a tendency to get trapped in local minima.

74 ATOMIC AND MOLECULAR PHYSICS↗

Proton rings from late-forming ballistic sheath fields

Many laser-driven ion experiments have seen ring-like patterns in the proton angular distribution across a wide array of laser and target parameters. These rings can impede measurement due to the small acceptance angle of detectors and often inhibit potential applications. A myriad of explanations for their formation have been proposed, yet most studies attribute them to some aspect of the laser–plasma interaction. Using 3D particle-in-cell simulations, we show that late-forming strong radial electric fields can arise due to charge separation while the beam is in flight, long after the laser–plasma interaction. These fields can accelerate ions to significant divergences (≈ 10°) as they propagate away from the target. We compare our results to a recent experiment where a high intensity, short pulse laser (I 0 ≈ 10 21 W/cm 2 , τ ≈ 30 fs) was incident upon thin (≈ 1 μm) liquid crystal targets. Our simulations capture all the main features of the experimental results—namely, robust ring formation and larger rings for higher energy protons. In addition, we show that rings do not form for sufficiently short preplasma scale lengths. Finally, we develop a phenomenological model to describe the spatiotemporal structure of the radial electric field and use this to explain the proton rings' energy and preplasma dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Liquid Structure of Shock-Compressed Hydrocarbons at Megabar Pressures

We present here results for the ionic structure in hydrocarbons (polystyrene, polyethylene) that were shock compressed to pressures of up to 190 GPa, inducing rapid melting of the samples. The structure of the resulting liquid is then probed using in situ diffraction by an x-ray free electron laser beam, demonstrating the capability to obtain reliable diffraction data in a single shot, even for low-Z samples without long range order. The data agree well with ab initio simulations, validating the ability of such approaches to model mixed samples in states where complex interparticle bonds remain, and showing that simpler models are not necessarily valid. While the results clearly exclude the possibility of complete carbon-hydrogen demixing at the conditions probed, they also, in contrast to previous predictions, indicate that diffraction is not always a sufficient diagnostic for this phenomenon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nanoscale view of assisted ion transport across the liquid–liquid interface

During solvent extraction, amphiphilic extractants assist the transport of metal ions across the liquid–liquid interface between an aqueous ionic solution and an organic solvent. Investigations of the role of the interface in ion transport challenge our ability to probe fast molecular processes at liquid–liquid interfaces on nanometer-length scales. Recent development of a thermal switch for solvent extraction has addressed this challenge, which has led to the characterization by X-ray surface scattering of interfacial intermediate states in the extraction process. Here, we review and extend these earlier results. We find that trivalent rare earth ions, Y(III) and Er(III), combine with bis(hexadecyl) phosphoric acid (DHDP) extractants to form inverted bilayer structures at the interface; these appear to be condensed phases of small ion–extractant complexes. The stability of this unconventional interfacial structure is verified by molecular dynamics simulations. The ion–extractant complexes at the interface are an intermediate state in the extraction process, characterizing the moment at which ions have been transported across the aqueous–organic interface, but have not yet been dispersed in the organic phase. In contrast, divalent Sr(II) forms an ion–extractant complex with DHDP that leaves it exposed to the water phase; this result implies that a second process that transports Sr(II) across the interface has yet to be observed. Calculations demonstrate that the budding of reverse micelles formed from interfacial Sr(II) ion–extractant complexes could transport Sr(II) across the interface. Our results suggest a connection between the observed interfacial structures and the extraction mechanism, which ultimately affects the extraction selectivity and kinetics.

Liang, Zhu↗

CO 2 in Ionene–Ionic Liquid Composite Membranes

Abstract Ionene – ionic liquid (IL) composites are promising materials for CO 2 separation, yet a molecular‐level understanding of their structure and its impact on CO 2 speciation, solubility, rotation, and diffusivity remains unclear. Herein, using multimodal nuclear magnetic resonance (NMR), time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS), atomic force microscopy (AFM), and molecular dynamics (MD) simulations, we reveal that the composites contain IL‐rich domains extending across hundreds of nanometres within the ionene matrix, and these bicontinuous domains span the entire membrane depth. CO 2 also absorbs into the ionene matrix, with the distribution between two CO 2 species varying with temperature and time. The rotational correlation times of these two species are on the timescale of 0.1 and 1 ns, respectively. As IL content increases, the ionic domains expand, resulting in higher CO 2 solubility due to enhanced molecular dynamics and increased free volume in both ionene backbones and IL‐rich regions. Although CO 2 diffusion in the membranes is an order of magnitude slower than in bulk IL, the activation energy for CO 2 diffusion remains comparable. Ionene‐IL composites represent a promising platform for designing CO 2 separation membranes, offering enhanced CO 2 diffusion and selectivity through IL‐rich domains, and increased CO 2 solubility and mechanical integrity from the ionene matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗

Characterization of the molten salt FMgNaK through ab initio molecular dynamics and experimental density measurements

The selection of a salt species for use in a molten salt reactor MSR is a key part of any MSR design. However, many salts have sparse or no thermophysical property data sets, especially those with higher melting points. One such salt is the eutectic mixture of (NaF)0.345(KF)0.59(MgF 2 )0.065, or FMgNaK, and is explored here through ab initio molecular dynamics (AIMD) simulation (1023-1273K) as well as experimental measurement of the liquid density using the Archimedean method (973-1223K). Predicted densities from AIMD are similar to experimentally-measured densities, though with more noise, suggesting that the small scale necessitated by AIMD simulations may be problematic for simulating FMgNaK in particular. Here the coefficient of thermal expansion is predicted from simulation, and salt structure is characterized. Mg-F-Mg chaining is observed in the salt network, though the low concentration of Mg inhibits chaining on the scale that is observed in (LiF) 0.67 (BeF 2 ) 0.33 , or FLiBe.

36 MATERIALS SCIENCE↗

Computational Study of Short-Pulse Laser-Induced Generation of Crystal Defects in Ni-Based Single-Phase Binary Solid–Solution Alloys

Single-phase concentrated solid–solution alloys exhibit enhanced mechanical characteristics and radiation damage resistance, making them promising candidate materials for applications involving an exposure to rapid localized energy deposition. We use large-scale atomistic modeling to investigate the mechanisms of the generation of vacancies, dislocations, stacking faults, and twin boundaries in Ni, Ni 50 Fe 50 , Ni 80 Fe 20 , and Ni 80 Cr 20 targets irradiated by short laser pulses in the regime of melting and resolidification. The decrease in the thermal conductivity and strengthening of the electron–phonon coupling due to the intrinsic chemical disorder in the solid-solution alloys are found to have important implications on localization of the energy deposition and generation of thermoelastic stresses. The interaction of the laser-induced stress waves with the melting front is found to play a key role in roughening of the crystal–liquid interface and generation of dislocations upon the solidification. A common feature revealed in the structural analysis of all irradiated targets is the presence of high vacancy concentrations exceeding the equilibrium values at the melting temperature by about an order of magnitude. On the basis of the results of molecular dynamics simulations of solidification occurring at fixed levels of undercooling, the generation of vacancies is correlated with the velocity of the solidification front, and the processes responsible for creating the strong vacancy supersaturation are revealed. The suppression of the vacancy generation in the solid–solution alloys is also revealed and related to combined effect of enhanced vacancy mobility and higher energy of the vacancy formation in the alloy systems. The analysis of the first atomic shells surrounding the vacancy sites in Ni–Fe alloys uncovers the preference for the vacancy sites to be surrounded by Fe atoms and suggests that atomic-scale chemical heterogeneities may play an important role in defining the behavior and properties of the single-phase concentrated solid-solution alloys.

36 MATERIALS SCIENCE↗

Latent Heating Retrievals Using the TRMM Precipitation Radar: A Multi-Seasonal Study

Rainfall is a key link in the hydrologic cycle and is a primary heat source for the atmosphere. The vertical distribution of latent heat release, which is accompanied by rainfall, modulates the large-scale circulations of the tropics and in turn can impact midlatitude weather. This latent heat release is a consequence of phase changes between vapor, liquid, and solid water. Present largescale weather and climate models can simulate latent heat release only crudely, thus reducing their confidence in predictions on both global and regional scales. This paper represents the first attempt to use NASA Tropical Rainfall Measuring Mission (TRMM) rainfall information to estimate the four-dimensional structure of global monthly latent heating profiles over the global tropics from December 1997 to October 2000. The Goddard Convective-Stratiform. Heating (CSH) algorithm and TRMM precipitation radar data are used for this study. We will examine and compare the latent heating structures between 1997-1998 (winter) ENSO and 1998-2000 (non-ENSO). We will also examine over the tropics. The seasonal variation of heating over various geographic locations (i.e., oceanic vs continental; Indian oceans vs west Pacific; Africa vs S. America) will be also examined and compared. In addition, we will examine the relationship between latent heating (max heating level) and SST. The period of interest also coincides with several TRMM field campaigns that recently occurred over the South China Sea in 1998 (SCSMEX), Brazil in 1999 (TRMM-LBA), and in the central Pacific in 1999 (KWAJEX). Sounding diagnosed Q1 budgets from these experiments could provide a means of validating the retrieved profiles of latent heating from the CSH algorithm.

Tao, Wei-Kuo↗

Transient cavities in liquids and the nature of the hydrophobic effect

The size distributions of transient cavities in water and organic liquids, obtained from computer simulations, have provided a new means to analyze the nature of the hydrophobic effect and to evaluate the adequacy of different analytical models of this effect. The poor solubility of non-polar solutes in water is attributed to a low probability of finding in water cavities of atomic and molecular size. It has been shown that water applies more force per unit area of cavity surface than do hydrocarbon liquids. Models that successfully capture the main characteristics of the hydrophobic effect must at least include information about the density and the radial distribution of oxygen atoms in liquid water. One such model, quantitatively accurate for molecular solutes of arbitrary shape, is presented.

NASA Center ARC↗

Atomic-level mechanisms of short-circuit diffusion in materials

This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.

36 MATERIALS SCIENCE↗

Molecular Investigation of Oxidized Graphene: Anatomy of the Double-Layer Structure and Ion Dynamics

We investigated the influence of surface oxidization of planar graphene electrodes on charge storage and ion dynamics of supercapacitors. Our approach compared two distinct ionic liquid (IL) electrolytes: tetraethylammonium tetrafluoroborate solvated in acetonitrile (TEA-BF4/ACN) and 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide solvated in acetonitrile (EMIm-TFSI/ACN). Both experimental electrochemical tests and molecular dynamics (MD) simulations showed positive, electrolyte-specific influences of hydroxyl-free electrode interfaces on capacitance. In the EMIm-TFSI/ACN system, the hydroxylated surface, because of its strong interaction with anions, impeded surface charge storage. On the other hand, in the case of TEA-BF4/ACN, the distribution and orientation of ACN across the system exerted vital influence on the capacitance, especially on the positive hydroxyl-free electrode. Furthermore, MD simulations of ion mobility with respect to the electrode surface in the lateral and perpendicular directions revealed significantly slower diffusion performance on the oxidized surface. Our efforts enhanced the level of fundamental understanding of the effects of hydroxyl groups on electrode–electrolyte interfaces and resulting supercapacitor performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reconfigurable Multicompartment Emulsion Drops Formed by Nematic Liquid Crystals and Immiscible Perfluorocarbon Oils

Liquid crystalline (LC) oils provide the basis of stimuli-responsive LC-in-water emulsions. Although past studies have explored the properties of single-phase LC emulsions, few studies have focused on complex multicompartment emulsions containing co-existing isotropic and LC domains. In this paper, we report a study of multiphase emulsions using LCs and immiscible perfluoroalkanes dispersed in water or glycerol (the latter continuous phase is used to enable characterization). We found that the nematogen 4'-pentyl-4-biphenylcarbonitrile (5CB) anchors homeotropically (perpendicularly) and weakly at liquid perfluorononane (F9) interfaces, consistent with the smectic layering of 5CB molecules. The proposed role of smectic layering is supported by experiments performed with 4-(trans-4-pentylcyclohexyl)benzonitrile, a nematogen that possesses a cyclohexyl group that frustrates the smectic packing and leads to tilted orientations at the F9 interface. By employing perfluorocarbon and hydrocarbon surfactants in combination with multiphase 5CB and F9 emulsion droplets dispersed in a continuous water or glycerol phase, we observe a range of emulsion droplet morphologies to form, including core–shell and Janus structures, with internal organizations that reflect an interplay of interfacial (anchoring energies; F9 and glycerol) and elastic energies within the confines of the geometry of the emulsion droplet. By comparing experimental observations to simulations of the LC–perfluorocarbon droplets based on a Landau-de Gennes model of the free energy, we place bounds on the orientation-dependent interfacial energies that underlie the internal ordering of these complex emulsions. Additionally, by forming core–shells emulsion droplets from 5CB (shell) and perfluoroheptane (cores), we demonstrate how a liquid-to-vapor phase transition in the perfluorocarbon core can be used to actuate the droplet and rapidly thin the nematic shell. Overall, the findingss reported in this paper demonstrate that multiphase LC emulsions formed from mixtures of perfluoroalkanes and LCs provide new opportunities to engineer hierarchical and stimuli-responsive emulsion systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

In this study, a method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob–Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Infrared Spectroscopy of Li + Solvation in EmimBF 4 and in Propylene Carbonate: Ab Initio Molecular Dynamics and Experiment

Infrared (IR) spectra of solutions of the lithium salt LiBF 4 in the ionic liquid 1-ethyl-3-methylimidazolium tetrafluoroborate (EmimBF 4 ) and in the organic solvent propylene carbonate (PC) are studied via infrared spectroscopy and ab initio molecular dynamics (AIMD) simulations. The measurements show that the major effects of LiBF 4 in both solutions, compared to their neat counterparts, are the appearance of a new broad band in the 300-450 cm -1 frequency region and a broadening of the IR structure in the 900-1200 cm -1 region with the development of a new peak at 980 cm -1 . Computational analysis indicates that hindered translational motions of Li + in its solvation cage are mainly responsible for the former, while the latter is due to Li + -induced structural changes and accompanying vibrational frequency shifts of constituent ions and molecules of the solutions. In addition, molecular motions in these and lower frequency regions are generally correlated, disclosing the collective nature of the vibrational dynamics, which involve multiple ions/molecules. Herein, detailed analysis of these features via AIMD simulations of the spectrum and its components arising from auto- and cross-correlations of motions of constituent molecular species, combined with generalized normal modes of the solutions and normal modes of small Li + -containing clusters, is presented. Furthermore, other minor spectral changes caused by the lithium salt as well as the interaction-induced effect on IR spectra are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗