Search NASA⌕ Search

SEARCH · Search NASA

Results for “CCSD Computations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

GPU acceleration of rank-reduced coupled-cluster singles and doubles

Here, we have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSD introduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSD tensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSD equations for up to 400 electrons and 1550 basis functions—roughly 50% larger than the largest canonical CCSD computations that have been performed on any hardware. In addition to increased scalability, the RR-CCSD computations are faster than the corresponding CCSD computations for all but the smallest molecules. We test the accuracy of RR-CCSD for a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1 kcal mol -1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSD doubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSD computation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.

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

On the functional dependence of transition-potential coupled cluster

Orbital relaxation of the core region is a primary source of error in the computation of core ionization and core excitation energies. Recently, Transition-Potential Coupled Cluster (TP-CC) methods have been used to explicitly treat orbital relaxation using non-variational molecular orbitals determined by reoccupation of orbitals optimized for a fractional core occupation. The amount of fractional occupation is governed by parameter λ, and recommended values for accurate TP-CCSD and XTP-CCSD computations of carbon, nitrogen, oxygen, and fluorine K edges were previously determined. Herein, we explore the performance of several density functionals for generating the fractionally occupied orbitals used in TP-CCSD. These functionals include HF, BP86, BH&HLYP, B3LYP, M06-2X, and ωB97m-V. The fractionally occupied orbitals computed across the various functionals were subsequently employed as the initial orbitals for our TP-CCSD calculations of organic K-edge x-ray absorption and photoelectron spectra. Regardless of the functional used to generate the fractionally occupied orbitals, the TP-CCSD calculations yield accurate and comparable core ionization energies, core excitation energies, and oscillator strengths.

Coupled-cluster methods↗
DOE OSTI2014

Communication: Resolving the three-body contribution to the lattice energy of crystalline benzene: Benchmark results from coupled-cluster theory

Coupled-cluster theory including single, double, and perturbative triple excitations [CCSD(T)] has been applied to trimers that appear in crystalline benzene in order to resolve discrepancies in the literature about the magnitude of non-additive three-body contributions to the lattice energy. In this paper, the present results indicate a non-additive three-body contribution of 0.89 kcal mol -1 , or 7.2% of the revised lattice energy of -12.3 kcal mol -1 . For the trimers for which we were able to compute CCSD(T) energies, we obtain a sizeable difference of 0.63 kcal mol -1 between the CCSD(T) and MP2 three-body contributions to the lattice energy, confirming that three-body dispersion dominates over three-body induction. Taking this difference as an estimate of three-body dispersion for the closer trimers, and adding an Axilrod-Teller-Muto estimate of 0.13 kcal mol -1 for long-range contributions yields an overall value of 0.76 kcal mol -1 for three-body dispersion, a significantly smaller value than in several recent studies.

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

tert -Butyl peroxy radical: ground and first excited state energetics and fundamental frequencies

Alkylperoxy radicals (RO 2 ˙) are key intermediates in combustion and atmospheric oxidation processes. As such, reliable detection and monitoring of these radicals can provide a wealth of information about the underlying chemistry. The tert-butyl peroxy radical is the archetypal tertiary peroxy radical, yet its vibrational spectroscopy is largely unexplored. To aid in future experimental investigations, we have performed high-level theoretical studies of the fundamental vibrational frequencies of the ground- and first excited states. A conformer search on both electronic surfaces reveals single minimum-energy structures. We predict an à 2 A' ← X 2 A" adiabatic excitation energy of 7738 cm -1 via focal point analysis, approximating the CCSDT(Q)/CBS level of theory. This excitation energy agrees to within 17 cm -1 of the most accurate experimental measurement. We compute CCSD(T) fundamental vibrational frequencies via second-order vibrational perturbation theory (VPT2), using a hybrid force field in which the quadratic (cubic/quartic) force constants are evaluated with the ANO1 (ANO0) basis set. Anharmonic resonance polyads are treated with the VPT2 + K effective Hamiltonian approach. Among the predicted fundamental frequencies, the ground state O–O stretch, excited state O–O stretch, and excited state C–O–O bend fundamentals are predicted at 1138, 959, and 490 cm -1 , respectively. Basis set sensitivity is found to be particularly great for the O–O stretches, similar to what has already been noted in smaller, unbranched peroxy radicals. Exempting these O–O stretches, agreement with the available experimental fundamentals is generally good (±10 cm -1 ).

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

Next Generation, Explicitly Correlated Electronic Structure Methods

All common wave function methods of electronic structure theory have a fundamental and troublesome flaw: the inability to correctly describe the exact mathematical cusp behavior of many-electron wave functions in the vicinity of coalescence points (Coulomb singularities), and hence to fully account for instantaneous, short-range correlation among electrons. To achieve high thermochemical accuracy, next-generation methodologies were needed that solve the electron cusp problem by explicitly incorporating interelectronic variables (r 12 ) into molecular wave functions. We pursued a long-term research program to develop practical explicitly-correlated methods, to create attendant, freely available, state-of-the-art computer codes, and to undertake chemical applications of unprecedented size. Our work allows conventional CCSD(T) computations to be conjoined with explicitly-correlated second-order perturbation methods within a focal point extrapolation scheme for pinpointing the ab initio limit of electronic structure theory.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗
DOE OSTI2015

The benzene⋯naphthalene complex: A more challenging system than the benzene dimer for newly developed computational methods

High-level coupled cluster singles, doubles, and perturbative triples [CCSD(T)] computations with up to the aug-cc-pVQZ basis set (1924 basis functions) and various extrapolations toward the complete basis set (CBS) limit are presented for the sandwich, T-shaped, and parallel-displaced benzene⋯naphthalene complex. Using the CCSD(T)/CBS interaction energies as a benchmark, the performance of some newly developed wave function and density functional theory methods has been evaluated. The best performing methods were found to be the dispersion-corrected PBE0 functional (PBE0-D3) and spin-component scaled zeroth-order symmetry-adapted perturbation theory (SCS-SAPT0). The success of SCS-SAPT0 is very encouraging because it provides one method for energy component analysis of π-stacked complexes with 200 atoms or more. Most newly developed methods do, however, overestimate the interaction energies. The results of energy component analysis show that interaction energies are overestimated mainly due to the overestimation of dispersion energy.

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

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

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

Scattering study of the Ne + NeH{sup +}(v{sub 0} = 0, j{sub 0} = 0) → NeH{sup +} + Ne reaction on an ab initio based analytical potential energy surface

Initial state selected dynamics of the Ne + NeH{sup +}(v{sub 0} = 0, j{sub 0} = 0) → NeH{sup +} + Ne reaction is investigated by quantum and statistical quantum mechanical (SQM) methods on the ground electronic state. The three-body ab initio energies on a set of suitably chosen grid points have been computed at CCSD(T)/aug-cc-PVQZ level and analytically fitted. The fitting of the diatomic potentials, computed at the same level of theory, is performed by spline interpolation. A collinear [NeHNe]{sup +} structure lying 0.72 eV below the Ne + NeH{sup +} asymptote is found to be the most stable geometry for this system. Energies of low lying vibrational states have been computed for this stable complex. Reaction probabilities obtained from quantum calculations exhibit dense oscillatory structures, particularly in the low energy region and these get partially washed out in the integral cross section results. SQM predictions are devoid of oscillatory structures and remain close to 0.5 after the rise at the threshold thus giving a crude average description of the quantum probabilities. Statistical cross sections and rate constants are nevertheless in sufficiently good agreement with the quantum results to suggest an important role of a complex-forming dynamics for the title reaction.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Optimization of the Coupled Cluster Implementation in NWChem on Petascale Parallel Architectures

Coupled cluster singles and doubles (CCSD) algorithm has been optimized in NWChem software package. This modification alleviated the communication bottleneck and provided from 2- to 5-fold speedup in the CCSD iteration time depending on the problem size and available memory. Sustained 0.60 petaflop/sec performance on CCSD(T) calculation has been obtained on NCSA Blue Waters. This number included all stages of the calculation from initialization till termination, iterative computation of single and double excitations, and perturbative accounting for triple excitations. In the section of perturbative triples alone, the computation maintained 1.18 petaflop/sec performance level. CCSD computations have been performed on Guanine-Cytosine deoxydinucleotide monophosphate (GC-dDMP) to probe the conformational energy difference in DNA single strand in A- and B-conformations. The computation revealed significant discrepancy between CCSD and classical force fields in prediction of relative energy of A- and B-conformations of GC-dDMP.

calculations↗
DOE OSTI2018

Gas Phase Hydrolysis and Oxo-Exchange of Actinide Dioxide Cations: Elucidating Intrinsic Chemistry from Protactinium to Einsteinium

Gas-phase bimolecular reactions of metal cations with water provide insights into intrinsic characteristics of hydrolysis. For the actinide dioxide cations, actinyl(V) AnO 2 + , melding of experiment and computation provides insights into trends for hydrolysis, as well as for oxo-exchange between actinyls and water that proceeds by a hydrolysis pathway. Here this line of inquiry is further extended into the actinide series with CCSD(T) computations of potential energy surfaces, for the reaction pathway for oxo-exchange through hydrolysis of nine actinyl(V) ions, from PaO 2 + to EsO 2 + . The computed surfaces are in accord with previous experimental results for oxo-exchange, and furthermore predict spontaneous exchange for CmO 2 + , BkO 2 + , CfO 2 + and EsO 2 + , but not for AmO 2 + . Natural Bond Order analysis of the species involved in both hydrolysis and oxo-exchange reveals an inverse correlation between the barrier to hydrolysis and the charge on the actinide centre, q(An). Based on this correlation, it can be concluded that hydrolysis, and related phenomena such as oxo-exchange, become less favourable as the charge on the metal centre decreases. Finally, the new results provide a straightforward rationalization of trends across a wide swathe of the actinide series.

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

Implementation of analytic gradients for CCSD and EOM-CCSD using Cholesky decomposition of the electron-repulsion integrals and their derivatives: Theory and benchmarks

We present a general formulation of analytic nuclear gradients for the coupled-cluster with single and double substitution (CCSD) and equation-of-motion (EOM) CCSD energies computed using Cholesky decomposition (CD) representations of the electron repulsion integrals. By rewriting the correlated energy and response equations such that the storage of the largest four-index intermediates is eliminated, CD leads to a significant reduction in disk storage requirements, reduced I/O penalties, and an improved parallel performance. CD thus extends the scope of the systems that can be treated by (EOM-)CCSD methods, although analytic gradients in the framework of CD are needed to extend the applicability of (EOM-)CCSD methods in the context of geometry optimizations. This paper presents a formulation of analytic (EOM-)CCSD gradient within the CD framework and reports on the salient details of the corresponding implementation. The accuracy and the capabilities of analytic CD-based (EOM-)CCSD gradients are illustrated by benchmark calculations and several illustrative examples.

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

Convergent Concordant Mode Approach for Molecular Vibrations: CMA-2

The concordant mode approach (CMA) is a promising new scheme for dramatically increasing the system size and level of theory achievable in quantum chemical computations of molecular vibrational frequencies. Here, we achieve advances in the CMA hierarchy by computations targeting CCSD(T)/cc-pVTZ (coupled cluster singles and doubles with perturbative triples using a correlation-consistent polarized-valence triple-ζ basis set) benchmarks within the G2 molecular test set, executing a statistical analysis for 1501 frequencies from 111 compounds and then separately solving the refractory case of pyridine. First, MP2/cc-pVTZ (second-order Møller–Plesset perturbation theory with the same basis set) proves to be an excellent and preferred choice for generating the underlying (Level B) normal modes of the CMA scheme. Utilizing this Level B within the CMA-0A method reproduces the 1501 benchmark frequencies with a mean absolute error (MAE) of only 0.11 cm –1 and an attendant standard deviation of 0.49 cm –1 . Second, a convergent CMA-2 method is constituted that allows efficient computation of higher level (Level A) frequencies to any reasonable accuracy threshold by using only Hartree–Fock (HF) and MP2 or density functional theory (DFT) data to generate ξ parameters, which select the sparse off-diagonal force field elements for explicit evaluation at Level A. When Level B = MP2/cc-pVTZ, a cutoff of ξ = 0.02 provides an average maximum absolute error per molecule of only 0.17 cm –1 by incurring merely a 33% increase in average cost over CMA-0A. This CMA-2 method also eradicates the 4 problematic CMA-0A outliers of pyridine with even less effort (ξ = 0.04, 22% increase). Finally, the newly developed CMA procedures are shown to be highly successful when applied to 1-(1H-pyrrol-3-yl)ethanol, a new test molecule with diverse types of vibration.

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

Remarkably High Stability of Late Actinide Dioxide Cations: Extending Chemistry to Pentavalent Berkelium and Californium

Actinyl chemistry is extended beyond Cm to BkO 2 + and CfO 2 + through transfer of an O atom from NO 2 to BkO + or CfO + , establishing a surprisingly high lower limit of 73 kcal mol -1 for the dissociation energies, D[O-(BkO + )] and D[O-(CfO + )]. CCSD(T) computations are in accord with the observed reactions, and characterize the newly observed dioxide ions as linear pentavalent actinyls; these being the first Bk and Cf species with oxidation states above IV. Computations of actinide dioxide cations AnO 2 + for An=Pa to Lr reveal an unexpected minimum for D[O-(CmO + )]. For CmO 2 + , and AnO 2 + beyond EsO 2 + , the most stable structure has side-on bonded η 2 -(O 2 ), as An III peroxides for An=Cm and Lr, and as An II superoxides for An=Fm, Md, and No. It is predicted that the most stable structure of EsO 2 + is linear [O=Es V =O] + , einsteinyl, and that FmO 2 + and MdO 2 + , like CmO 2 + , also have actinyl(V) structures as local energy minima. Finally, the results expand actinide oxidation state chemistry, the realm of the distinctive actinyl moiety, and the non-periodic character towards the end of the periodic table.

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

Improved Accuracy in Semi-Experimental Structure Determination by Resolving Problems Associated with Rotation of Principal Inertial Axes of Isotopologues: Structures of 1,3-Oxazole ( c -C 3 H 3 NO)

The rotational spectrum of the normal isotopologue of 1,3-oxazole (c-C 3 H 3 NO) was observed from 43 to 750 GHz. Over 3900 transitions for the ground vibrational state are measured, assigned, and least-squares fit to sextic centrifugally distorted-rotor Hamiltonians. The measured frequencies and resulting spectroscopic constants from this extended spectral range, combined with previous measurements of the nuclear quadrupole coupling constants, will facilitate astronomical searches for oxazole across the majority of the range of modern radiotelescopes. Spectra for a set of 30 oxazole isotopologues, which include multiple isotopic substitutions of each atom, are used to determine the first semi-experimental equilibrium ($r$$^{SE}_{e}$) structure and semi-experimental substitution structure ($r$$^{SE}_{e}$), each using CCSD(T) computed values for the vibration–rotation interaction and electron-mass corrections. The large number of isotopologues, including 21 isotopologues observed for the first time, and the redundant substitutions of each atom provide sufficient spectroscopic information to determine the $r$$^{SE}_{e}$ structure with the expected high level of accuracy and precision (0.0001 or 0.0002 Å in bond distances and 0.013 to 0.025° in bond angles). In the course of this study, we analyzed a known issue for some $r$$^{SE}_{e}$ structure determinations of near-oblate asymmetric tops in which inclusion of individual isotopologues degrades the structure determination. We demonstrate that this problem primarily arises from the difference in the values of the computed vibration–rotation interaction corrections as evaluated at the computed re geometry vs the $r$$^{SE}_{e}$ geometry of the “real” molecule. Our solution to this problem substantially improves the $r$$^{SE}_{e}$ structure of oxazole and likely can be generalized to many other molecules.

Chemical structure↗
DOE OSTI2014

Anion A– • HX Clusters with Reduced Electron Binding Energies: Proton vs Hydrogen Atom Relocation Upon Electron Detachment

Clustering an anion with one or more neutral molecules is a stabilizing process that enhances the oxidation potential of the complex relative to the free ion. Several hydrogen bond clusters (i.e., A— • HX, where A— = H2PO4— and CF3CO2— and HX = MeOH, PhOH, and Me2NOH or Et2NOH) are examined by photoelectron spectroscopy and M06-2X and CCSD(T) computations. Remarkably, these species are experimentally found to have adiabatic detachment energies that are smaller than those for the free ion and reductions of 0.47 to 1.87 eV are predicted computationally. Hydrogen atom and proton transfers upon vertical photodetachment are two limiting extremes on the neutral surface in a continuum of mechanistic pathways that account for these results, and the whole gamut of possibilities are predicted to occur.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2024

High-resolution infrared spectroscopy of jet cooled cyclobutyl in the α-CH stretch region: large-amplitude puckering dynamics in a 4-membered ring radical

We report gas-phase cyclobutyl radical (c-C 4 H 7 ) is generated at a rotational temperature of T rot = 26(1) K in a slit-jet discharge mixture of 70% Ne/30% He and 0.5–0.6% cyclobromobutane (c-C 4 H 7 Br). A fully rovibrationally resolved absorption spectrum of the α-CH stretch fundamental band between 3062.9 cm −1 to 3075.7 cm −1 is obtained and analyzed, yielding first precision structural and dynamical information for this novel radical species. The α-CH stretch band origin is determined to be 3068.7887(4) cm −1 , which implies only a modest (≈0.8 cm −1 ) blue shift from rotationally unresolved infrared spectroscopic studies of cyclobutyl radicals in liquid He droplets [ A. R. Brown, P. R. Franke and G. E. Douberly, J. Phys. Chem. A , 2017, 121 , 7576–7587]. Of particular dynamical interest, a one-dimensional potential energy surface with respect to the ring puckering coordinate is computed at CCSD(T)/ANO2 level of theory and reveals a double minimum C s puckered geometry, separated by an exceedingly shallow planar C 2v transition state barrier (E barr ≈ 1 cm −1 ). Numerical solutions on this double minimum potential yield a zero-point energy for the ground state (E zero-point ≈ 27 cm −1 ) greatly in excess of the interconversion barrier. This is indicative of highly delocalized, large amplitude motion of the four-membered ring structure, for which proper vibrationally averaging of the moment of inertia tensor reproduces the experimentally determined inertial defect remarkably well. Finally, intensity alternation in the experimental spectrum due to nuclear spin statistics upon exchange of three indistinguishable H atom pairs (IH = ½) matches K a + K c = even : odd = 36 : 28 predictions, implying that the unpaired electron in the radical center lies in an out-of-plane p π orbital. Thus, high-resolution infrared spectroscopy provides first experimental confirmation of a shallow double minimum ring puckering potential with a highly delocalized ground state wave function peaked at a planar C 2v transition state geometry consistent with a cyclobutyl π radical.

Chemistry↗
DOE OSTI2014

An exploration of the ozone dimer potential energy surface

The (O{sub 3}){sub 2} dimer potential energy surface is thoroughly explored at the ab initio CCSD(T) computational level. Five minima are characterized with binding energies between 0.35 and 2.24 kcal/mol. The most stable may be characterized as slipped parallel, with the two O{sub 3} monomers situated in parallel planes. Partitioning of the interaction energy points to dispersion and exchange as the prime contributors to the stability, with varying contributions from electrostatic energy, which is repulsive in one case. Atoms in Molecules analysis of the wavefunction presents specific O⋯O bonding interactions, whose number is related to the overall stability of each dimer. All internal vibrational frequencies are shifted to the red by dimerization, particularly the antisymmetric stretching mode whose shift is as high as 111 cm{sup −1}. In addition to the five minima, 11 higher-order stationary points are identified.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗