Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent density functional theory”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Cavity Ring-Down Spectroscopy of Anthracene, 9-Methylanthracene, and 2 Methylanthracene in Supersonic Expansion

The S0–S1 absorption spectra of anthracene (C14H10), 9-methylanthracene (C15H12), and 2-methylanthracene (C15H12) are measured in the ultraviolet region between 330 and 375 nm (26,666 to 30,303 cm–1) with cavity ring-down spectroscopy in supersonic free-jet expansions of argon. The associated vibronic band systems and their spectroscopic assignments are discussed and compared to previous studies performed using fluorescence excitation and dispersed fluorescence techniques. Density functional theory (DFT) calculations were carried out to study the structures and evaluate the vibrational transitions of the ground and excited states. Time-dependent DFT calculations of the first electronic excited states and Franck–Condon factor calculations were carried out to assist in the assignment of the experimentally measured vibronic bands. The vibronic spectra obtained in absorption agree well with fluorescence excitation spectra in terms of peak positions but exhibit different relative band intensities. We find a very good match between experimentally obtained vibronic line positions and the peak positions of the quantum chemically calculated Franck–Condon excitation lines.

Salma Bejaoui↗
DOE OSTI2025

Direct observation of key aluminum hydroxide prenucleation oligomers for gibbsite nucleation and crystallization in sodium aluminate solution by liquid ToF-SIMS

The mechanism of gibbsite (aluminum hydroxide) crystallization from highly alkaline solutions such as Bayer liquors remains poorly understood, where aluminum (Al) transforms from largely tetrahedrally coordinated aluminate monomers in sodium aluminate solutions into a network of octahedra in gibbsite crystals. A variety of traditional analytical approaches applied to this system do not readily reveal the presence of higher-order oligomeric intermediates. To overcome this limitation, we employed in-situ liquid Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) to examine the Al species present in concentrated sodium aluminate solutions favoring crystallization of gibbsite or sodium aluminate. A complex mixture of Al oligomers was found. By comparing the change in the relative concentration of Al oligomers with +1 and -1 charge, we were able to identify three major Al oligomer candidates, iso-tetramers, iso-pentamers, and cyclic-hexamers, for the nucleation and crystallization of gibbsite. The concentrations of iso-tetramers and iso-pentamers significantly surpass those of cyclic-hexamers. Time-dependent in-situ Raman spectroscopy analysis indicated that the appearance of gibbsite coincided with the peak concentration of these oligomers. The Density-functional theory (DFT) calculation suggests that the formation of iso-oligomers is more favorable than that of cyclic-hexamers. The combined results suggest that iso-tetramers and iso-pentamers play the most substantial role in the nucleation and growth of gibbsite in the sodium aluminate solutions. Our findings also suggest that the oligomers that promote gibbsite crystallization are more stable in dilute sodium aluminate solutions, making these solutions particularly suitable for efficient gibbsite crystallization. In conclusion, our study fills a major knowledge gap in understanding Al speciation that leads to the nucleation and crystallization of gibbsite in concentrated sodium aluminate solutions.

Bayer liquor↗
NASA NTRS2013

Association Mechanisms of Unsaturated C2 Hydrocarbons with Their Cations: Acetylene and Ethylene

The ion-molecule association mechanism of acetylene and ethylene with their cations is investigated by ab initio quantum chemical methods to understand the structures, association energies, and the vibrational and electronic spectra of the products. Stable puckered cyclic isomers are found as the result of first forming less stable linear and bridge isomers. The puckered cyclic complexes are calculated to be strongly bound, by 87, 35 and 56 kcal/mol for acetylene-acetylene cation, ethylene-ethylene cation and acetylene-ethylene cation, respectively. These stable complexes may be intermediates that participate in further association reactions. There are no association barriers, and no significant inter-conversion barriers, so the initial linear and bridge encounter complexes are unlikely to be observable. However, the energy gap between the bridged and cyclic puckered isomers greatly differs from complex to complex: it is 44 kcal/mol in C4H4 +, but only 6 kcal/mol in C4H8 +. The accurate CCSD(T) calculations summarized above are also compared against less computationally expensive MP2 and density functional theory (DFT) calculations for structures, relative energies, and vibrational spectra. Calculated vibrational spectra are compared against available experiments for cyclobutadiene cation. Electronic spectra are also calculated using time-dependent DFT.

ethylene↗
DOE OSTI2026

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗
DOE OSTI2024

A Cellular Automaton Simulation for Predicting Phase Evolution in Solid-State Reactions

New computational tools for solid-state synthesis recipe design are needed in order to accelerate the experimental realization of novel functional materials proposed by high-throughput materials discovery workflows. This work contributes a cellular automaton simulation framework for predicting the time-dependent evolution of intermediate and product phases during solid-state reactions as a function of precursor choice and amount, reaction atmosphere, and heating profile. The simulation captures the effects of reactant particle spatial distribution, particle melting, and reaction atmosphere. Reaction rates based on rudimentary kinetics are estimated using density functional theory data from the Materials Project and machine learning estimators for the melting point and the vibrational entropy component of the Gibbs free energy. The resulting simulation framework allows for the prediction of the likely outcome of a reaction recipe before any experiments are performed. We analyze five experimental solid-state recipes for BaTiO 3 , CaZrN 2 , and YMnO 3 found in the literature to illustrate the performance of the model in capturing reaction selectivity and reaction pathways as a function of temperature and precursor choice. This simulation framework offers an easier way to optimize existing recipes, aid in the identification of intermediates, and design effective recipes for yet unrealized inorganic solids in silico .

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

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Predicting Bone Mechanical State During Recovery After Long-Duration Skeletal Unloading Using QCT and Finite Element Modeling

During long-duration missions at the International Space Station, astronauts experience weightlessness leading to skeletal unloading. Unloading causes a lack of a mechanical stimulus that triggers bone cellular units to remove mass from the skeleton. A mathematical system of the cellular dynamics predicts theoretical changes to volume fractions and ash fraction in response to temporal variations in skeletal loading. No current model uses image technology to gather information about a skeletal site s initial properties to calculate bone remodeling changes and then to compare predicted bone strengths with the initial strength. The goal of this study is to use quantitative computed tomography (QCT) in conjunction with a computational model of the bone remodeling process to establish initial bone properties to predict changes in bone mechanics during bone loss and recovery with finite element (FE) modeling. Input parameters for the remodeling model include bone volume fraction and ash fraction, which are both computed from the QCT images. A non-destructive approach to measure ash fraction is also derived. Voxel-based finite element models (FEM) created from QCTs provide initial evaluation of bone strength. Bone volume fraction and ash fraction outputs from the computational model predict changes to the elastic modulus of bone via a two-parameter equation. The modulus captures the effect of bone remodeling and functions as the key to evaluate of changes in strength. Application of this time-dependent modulus to FEMs and composite beam theory enables an assessment of bone mechanics during recovery. Prediction of bone strength is not only important for astronauts, but is also pertinent to millions of patients with osteoporosis and low bone density.

Chang, Katarina L.↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗
NASA NTRS2004

Photoionization of Benzene and Small Polycyclic Aromatic Hydrocarbons in Ultraviolet-Processed Astrophysical Ices: A Computational Study

We employed density functional theory (DFT) calculations to model the photoionization behavior of benzene and small polycyclic aromatic hydrocarbons when they are embedded in a matrix of water ice in order to investigate issues raised by recent experimental work by Gudipati and Allamandola. The ionization energies of benzene, naphthalene, anthracene, and pyrene were found to be lowered by 1.5-2.1 eV in water ice. Low-lying vertical electronic excitation energies were computed with time-dependent DFT for both neutral and ionized species and are found in both cases to be remarkably unaffected by the ice matrix. Chemical behavior in ultraviolet-photoprocessed ices is also discussed, with a focus on electron recombination and pathways leading to phenol and analogous products.

Woon, D. E.↗
DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Eigenfunctions of the magnetospheric radial-diffusion operator

The applicability of eigenfunction methods to the theory of magnetospheric radial diffusion is investigated analytically, with a focus on the operator (Lambda) for radiation-belt and ring-current particles. It is shown that the eigenfunctions of Lambda can be expressed in terms of a linear combination of ordinary Bessel functions of the first and second kinds, and that the solution for the drift-averaged phase-space density can be expressed as the sum of a quasi-static solution and a time-dependent superposition of the eigenfunctions of Lambda. Results for illustrative and practical sample problems are presented in extensive tables and graphs and discussed in detail.

Schulz, Michael↗
DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗
DOE OSTI2024

Electronic structure and optical properties of halide double perovskites from a Wannier-localized optimally-tuned screened range-separated hybrid functional

Halide double perovskites are a chemically diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and ab initio many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here, in this study, we apply the WOT-SRSH functional to five halide double perovskites and compare the results with those obtained from other known functionals and previous GW calculations. We also use the approach as a starting point for GW calculations and we compute the band structures and optical absorption spectrum for Cs 2 AgBiBr 6 , using both time-dependent DFT and the GW-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Shock Waves in a Bose-Einstein Condensate

A paper presents a theoretical study of shock waves in a trapped Bose-Einstein condensate (BEC). The mathematical model of the BEC in this study is a nonlinear Schroedinger equation (NLSE) in which (1) the role of the wave function of a single particle in the traditional Schroedinger equation is played by a space- and time-dependent complex order parameter (x,t) proportional to the square root of the density of atoms and (2) the atoms engage in a repulsive interaction characterized by a potential proportional to | (x,t)|2. Equations that describe macroscopic perturbations of the BEC at zero temperature are derived from the NLSE and simplifying assumptions are made, leading to equations for the propagation of sound waves and the transformation of sound waves into shock waves. Equations for the speeds of shock waves and the relationships between jumps of velocity and density across shock fronts are derived. Similarities and differences between this theory and the classical theory of sound waves and shocks in ordinary gases are noted. The present theory is illustrated by solving the equations for the example of a shock wave propagating in a cigar-shaped BEC.

Kulikov, Igor↗
DOE OSTI2026

Mixed stochastic-deterministic density functional theoretic decomposition of Kubo–Greenwood conductivities in the projector augmented wave formalism

Pairing the accuracy of the Kohn–Sham density-functional framework with the efficiency of a stochastic algorithmic approach, mixed stochastic-deterministic density functional theory (mDFT) achieves a favorable computational scaling with system sizes and electronic temperatures. We employ the recently developed mDFT formalism to investigate the dynamic charge-transport properties of systems in the warm dense matter regime. The optical conductivity spectra are computed for single- and multi- component mixtures of carbon, hydrogen, and beryllium using two complementary approaches: Kubo–Greenwood in the mDFT picture and real-time time-dependent mDFT. We further devise a decomposition of the Onsager coefficients leading up to the Kubo–Greenwood spectra to exhibit contributions from the deterministic, stochastic, and mixed electronic state transitions at different incident photon energies.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Physics and Chemistry of Te and HgTe-based Ternary Melts

Historically, the theoretical treatment of the liquid phases has always been more difficult and complicated than that for solid and gas phases. A liquid has no lattice structure as crystalline solids and the atoms/molecules in the liquid can migrate through it relatively rapidly. On the other hand, it is also interacting with many other atoms/molecules so that the simplifications of the kinetic theory of gases cannot be employed. For more complicated liquids, such as the liquids of high ionicity and those containing hydrogen bonds and electric dipoles, the understanding is far from complete. At the same time, accurate information on the physics and chemistry of semiconductor melts is needed for the quantitative descriptions of the process of crystal growth from melt. The pre-crystallization phenomena in the liquid phase are critical because the properties of the grown crystals depend on the state and structure of the melt as well as the thermal history of the melt during solidification process. However, the data on the liquid phase, such as thermophysical properties of semiconductor melts are scarce, especially for the HgTe-based II-VI ternary compound semiconductors because of their high vapor pressure and extreme toxicity. Analysis of the thermophysical properties of the melt can provide information about structural transitions of the melt during the solidification process. From a broader point of view, the structure of liquids is much more complicated than the crystalline solids, especially the relaxation behavior through different thermal histories. The theory of hetero-phase fluctuations of liquids is applicable to any many-body systems including condensed-matter physics, field theory, physics of nuclear-matter, cosmology, biology and even sociology. This book summarizes the physics and chemistry from the experimental measurements and theoretical analyses of phase diagram, thermodynamic properties, density, thermal conductivity, viscosity, and electrical conductivity on the binary, pseudo-binary and ternary melts of the most advanced IR-detector material systems of HgCdTe and HgZnTe as well as the analyses of these results. The main objectives of this study are: (1) to provide the phase diagrams and thermodynamic properties of Hg-Cd-Te and Hg-Zn-Te systems through quantitatively fitting the experimental data by assuming an associated solution model for the liquid phase, (2) to experimentally measure the thermophysical properties of the Hg-Cd-Te and Hg-Zn-Te melts, including density, viscosity, electrical conductivity and thermal conductivity as functions of temperature and composition and (3) to enhance the fundamental knowledge of hetero-phase fluctuations and relaxation phenomena in the melts and extend our understanding of the solidification process in order to interpret the experimental results of crystal growth so as to improve the melt growth processes of the compound semiconductor. The physics and chemistry of Te and HgTe-based ternary melts were explored through the studies of the structural transformation during melting, the supercooling during solidification, the relaxation phenomena after rapid cooling of the melts and the metal-semiconductor transition in the melts through the analyses of electrical conductivity and Lorenz number. An in-depth study on the thermophysical properties and their time-dependent structural dynamic processes taking place in the vicinity of the solid-liquid phase transition of the narrow homogeneity range HgTe-based ternary semiconductors as well as the analysis of the homogenization process in the melt will also be presented.

Thermophysical properties↗
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗