Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster method”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Converging high-level coupled-cluster energetics via adaptive selection of excitation manifolds driven by moment expansions

A novel approach to rapidly converging high-level coupled-cluster (CC) energetics in an automated fashion is proposed. The key idea is an adaptive selection of excitation manifolds defining higher--than--two-body components of the cluster operator inspired by CC(P;Q) moment expansions. Further, the usefulness of the resulting methodology is illustrated by molecular examples where the goal is to recover the electronic energies obtained using the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT) when the noniterative triples corrections to CCSD fail.

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

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Benchmarking quantum trial wavefunctions for phaseless auxiliary-field quantum Monte Carlo

The phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method is a stochastic imaginary-time projection technique for computing ground-state properties of strongly correlated quantum systems, with accuracy that depends critically on the choice of trial wavefunction. Here, we investigate ph-AFQMC with trial states prepared using parameterized quantum circuits. In this work, we present a comprehensive benchmarking study of quantum trial wavefunctions spanning unitary coupled-cluster, Hamiltonian-informed, Jastrow-inspired, and adaptively constructed ansatze. The benchmarking evaluates accuracy, expressibility, and scalability of these ansatze within the QC-AFQMC framework. We test these ansatze on linear hydrogen chains under bond stretching and find that several ansatz families produce chemically accurate ph-AFQMC energies across the dissociation curve. We have performed simulations using the CUDA-Q quantum development platform on the GPU partition of the Perlmutter supercomputer. When comparing ansatze at similar numbers of variational parameters, we find that different ansatz families yield comparable ph-AFQMC results despite exhibiting substantially different variational energies, optimization costs, and circuit depths. Our results indicate that the variational energy of an ansatz is not always a reliable indicator of its quality for ph-AFQMC and reveal instances of over-parameterization. In the strongly correlated regime, trial wavefunctions obtained from adaptive ansatze, exemplified here by ADAPT-VQE with the UCCSD operator pool, can outperform their fixed-ansatz counterparts (UCCSD) in terms of projected energies while using substantially more compact circuits, providing a flexible route to optimize quantum resources within the ph-AFQMC framework.

Rofougaran, Rod [LBNL, Berkeley; Columbia U.; PNL,↗
DOE OSTI2021

Probing Clustering Dynamics between Silicon and PAA or LiPAA Slurries under Processing Conditions

This work explores the complex interplay between slurry aggregation, agglomeration, and conformation (i.e., shape) of poly(acrylic acid) (PAA)- and lithiated PAA-based silicon slurries as a function of the shear rate and the resulting slurry homogeneity. These values were measured by small-angle neutron scattering (SANS) and rheology-coupled ultra-SANS at conditions relevant to battery electrode casting. Different binder solution preparation methods, either a ball milling (BM) process or a planetary centrifugal mixing (PCM) process, dramatically modify the resulting polymer dynamics and organization around a silicon material. This is due to the different energy profiles of mixing where the more violent and higher energy PCM causes extensive breakdown and reformation of the binder, which is now likely in a branched conformation, while the lower energy BM results in simply lower-molecular weight linear polymers. In conclusion, the breakdown and reorganization of the polymer structure affect silicon slurry homogeneity, which affects subsequent electrode architecture.

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

High accuracy barrier heights, enthalpies, and rate coefficients for chemical reactions

Abstract Quantitative chemical reaction data, including activation energies and reaction rates, are crucial for developing detailed kinetic mechanisms and accurately predicting reaction outcomes. However, such data are often difficult to find, and high-quality datasets are especially rare. Here, we use CCSD(T)-F12a/cc-pVDZ-F12// ω B97X-D3/def2-TZVP to obtain high-quality single point calculations for nearly 22,000 unique stable species and transition states. We report the results from these quantum chemistry calculations and extract the barrier heights and reaction enthalpies to create a kinetics dataset of nearly 12,000 gas-phase reactions. These reactions involve H, C, N, and O, contain up to seven heavy atoms, and have cleaned atom-mapped SMILES. Our higher-accuracy coupled-cluster barrier heights differ significantly (RMSE of ∼5 kcal mol −1 ) relative to those calculated at ω B97X-D3/def2-TZVP. We also report accurate transition state theory rate coefficients $${k}_{\infty }(T)$$ k ∞ ( T ) between 300 K and 2000 K and the corresponding Arrhenius parameters for a subset of rigid reactions. We believe this data will accelerate development of automated and reliable methods for quantitative reaction prediction.

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

Is the Trotterized UCCSD Ansatz Chemically Well-Defined?

The variational quantum eigensolver (VQE) has emerged as one of the most enticing near-term quantum algorithms that can be used to simulate many-body systems such as molecular electronic structures. Serving as an attractive ansatz in the VQE algorithm, unitary coupled cluster (UCC) theory has seen a renewed interest in recent literature. Yet, unlike the original classical UCC theory, implementation on a quantum computer requires a finite-order Suzuki-Trotter decomposition to separate the exponentials of the large sum of Pauli operators. While previous literature has recognized the non-uniqueness of different orderings of the operators in the Trotterized form of UCC methods, the question of whether or not different orderings matter at the chemical scale has not been addressed. In this letter, we explore the effect of operator ordering on the Trotterized UCCSD ansatz, as well as the much more compact k-UpCCGSD ansatz recently proposed by Lee et al. We observe a significant, system-dependent variation in the energies of Trotterizations with different operator orderings. The energy variations occur on a chemical scale, sometimes on the order of hundreds of kcal/mol. This letter establishes the need to define not only the operators present in the ansatz, but also the order in which they appear. This is necessary for adhering to the quantum chemical notion of a ``model chemistry'', in addition to the general importance of scientific reproducibility. As a final note, we suggest a useful strategy to select out of the combinatorial number of possibilities, a single well-defined and effective ordering of the operators.

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

Identifying Sample Provenance From SEM/EDS Automated Particle Analysis via Few-Shot Learning Coupled With Similarity Graph Clustering

Automated particle analysis (APA) provides a vast amount of compositional data via energy-dispersive X-ray spectroscopy along with size and shape data via scanning electron microscopy for individual particles in a sample. In many instances, APA data are leveraged to support identification of the source of a sample based on the detection of particles of a specific composition. Often, the particles that provide context make up a minuscule portion of the sample. Additionally, the interpretation of complex samples can be difficult due to the diversity of compositions both in the mixture and within a particle. In this work, we demonstrate a method to compute and cluster similarity graphs that describe inter-particle relationships within a sample using a multi-modal few-shot learning neural network. Here, as a proof-of-concept, we show that samples known to have been exposed to gunshot residue can be distinguished from samples occasionally mistaken for gunshot residue. Our workflow builds upon standard APA techniques and data processing methods to unveil additional information in a readily interpretable and quantitatively comparable format.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Scalable Heterogeneous Execution of a Coupled-Cluster Model with Perturbative Triples

The CCSD(T) coupled-cluster model with perturbative triples is considered a gold standard for computational modeling of the correlated behavior of electrons in molecular systems. A fundamental constraint is the relatively small global-memory capacity in GPUs compared to the main-memory capacity on host nodes, necessitating relatively smaller tile sizes for high-dimensional tensor contractions in NWChem's GPU-accelerated implementation of the CCSD(T) method. A coordinated redesign is described to address this limitation and associated data movement overheads, including a novel fused GPU kernel for a set of tensor contractions, along with inter-node communication optimization and data caching. The new implementation of GPU-accelerated CCSD(T) improves overall performance by 3.4x. Finally, we discuss the trade-offs in using this fused algorithm on current and future supercomputing platforms.

Kim, Jinsung↗
DOE OSTI2020

Distributed non-negative matrix factorization with determination of the number of latent features

The holistic analysis and understanding of the latent (that is, not directly observable) variables and patterns buried in large datasets is crucial for data-driven science, decision making and emergency response. Such exploratory analyses require devising unsupervised learning methods for data mining and extraction of the latent features, and non-negative matrix factorization (NMF) is one of the prominent such methods. NMF is based on compute-intense non-convex constrained minimization, which, for large datasets requires fast and distributed algorithms. However, current parallel implementations of NMF fail to estimate the number of latent features. In practice, identifying these features is both difficult and significant for pattern recognition and latent feature analysis, especially for large dense matrices. Here, we introduce a distributed NMF algorithm coupled with distributed custom clustering followed by a stability analysis on dense data, which we call DnMFk, to determine the number of latent variables. The results on synthetic data and the classical Swimmer data set demonstrate the accuracy of model determination while scaling nearly linearly across multiple processors for large data. Further, we employ DnMFk to determine the number of hidden features from a terabyte matrix.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Gaunt and Breit two-electron contributions to mean-field transformations and fine structure splitting

Materials utilized by novel energy systems are often studied using weakly correlated mean-field theories. However, if these systems incorporate heavy elements, relativistic effects must be included. Therefore, a Kramers unrestricted coupled cluster with singles and doubles excitation formalism within a molecular mean-field exact two-component framework (X2C mmf ) using a four-component Dirac–Hartree–Fock (DHF) reference state is presented. The exact X2C mmf transformed normal-order Hamiltonian incorporates all one-electron and two-electron (2e) contributions from the Coulomb, Gaunt, and Breit operators and is used with the equation of motion method to calculate the excitation energies of the alkali group of elements. Using this framework, the effects of 2e Gaunt and Breit integrals are studied. Results demonstrate growing contributions from these integrals to the generated X2C mmf mean-fields and electronic fine structure calculations with increasing atomic number. Overall, this paper outlines the method, its effect within the X2C mmf approach, and lays the foundation for future theoretical development of relativistic calculations within this framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

27 Al NMR diffusometry of Al 13 Keggin nanoclusters

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al 13 , [Al 13 O 4 (OH) 24 (H 2 O) 12 ] 7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, in this work, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO 4 core of the ϵ-Al 13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO 4 center of the ϵ-Al 13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T 1 and T 2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes–Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al 13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered.

27Al↗
DOE OSTI2022

Sequential Bond Dissociation Energies of Th + (CO) x , x = 3–6: Guided Ion Beam Collision-Induced Dissociation and Quantum Computational Studies

Collision-induced dissociation (CID) of [Th,xC,xO] + , x = 3–6, with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). Products are formed exclusively by the loss of CO ligands. Analyses of the kinetic energy-dependent CID product cross sections yield bond dissociation energies (BDEs) of (CO) x–1 Th + –CO at 0 K as 1.09 ± 0.05, 0.82 ± 0.07, 0.63 ± 0.05, and 0.70 ± 0.05 eV, respectively. Different structures of [Th,xC,xO] + were explored using various electronic structure methods, and BDEs for CO ligand loss from precursor [Th,xC,xO] + complexes were computed. Both experimental and theoretical results corroborate that the structures of [Th,xC,xO] + , x = 3–6, formed experimentally are homoleptic thorium cation carbonyl complexes, Th + (CO) x . The nonmonotonic trend in experimental BDEs is reproduced theoretically, although ambiguities in the spin states of the x = 4–6 complexes (doublet or quartet) remain. BDEs calculated at the coupled cluster with single, double, and perturbative triple excitations (CCSD(T))/cc-pVXZ//B3LYP/cc-PVXZ (X = T and Q) level and a complete basis set (CBS) extrapolation agree reasonably well with the experimental values for all complexes. Thorium oxide ketenylidene carbonyl cations, OTh + CCO(CO) y , y = 1–4, were calculated to be the most stable structures of [Th,xC,xO] + , x = 3–6, respectively; however, these are not observed in our experiment. Potential energy profiles (PEPs) having either quartet or doublet spin calculated at the B3LYP/cc-pVQZ level suggest that the failure to observe OTh + CCO(CO) y , y = 1–4, is the result of a barrier corresponding to the C–C bond formation, making the formation of OTh + CCO(CO) y inaccessible kinetically under the present experimental conditions.

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

Main and Satellite Features in the Ni 2p XPS of NiO

The origin and assignment of the complex main and satellite XPS features of the cations in ionic compounds has been the subject of extensive theoretical studies using different methods. There is agreement that within a molecular orbital model one needs to take into account different types of configurations. Specifically, those where a core electron is removed but no other configuration changes are made and those where in addition to ionization there are also shake or charge transfer changes to the ionic configuration. However, there are strong disagreements about the assignment of XPS features to these configurations. The present work is directed toward resolving the origin of main and satellite features for the Ni 2p XPS of NiO based on ab initio molecular orbital wavefunctions for a cluster model of NiO. A major problem in earlier ab initio XPS studies of ionic compounds has been the use of a common set of orbitals that was not able to properly describe all the ionic configurations that contribute to the full XPS spectra. This is resolved in the present work by using orbitals that are optimized for averages of the occupations of the different configurations that contribute to the XPS. The approach of using State Averaged orbitals is validated through comparisons between different averages and through use of higher order excitations in the wavefunctions for the ionic states. It represents a major extension of our earlier work on the main and satellite features of the Fe 2p XPS of Fe 2 O 3 and proves the reliability and the generality of the assignments of the character and origin of the different features of the XPS obtained with orbitals optimized for State Averages. These molecular orbital methods permit the characterization of the ionic states in terms of the importance of shake excitations and of the coupling of ionization of 2p 1/2 and 2p 3/2 spin-orbit split sub shells. Further, the work lays the foundation for definitive assignments of the character of main and satellite XPS features and points to their origin in the electronic structure of the material.

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

Bacterial genome editing by coupling Cre-lox and CRISPR-Cas9 systems

The past decade has been a golden age for microbiology, marked by the discovery of an unprecedented increase in the number of novel bacterial species. Yet gaining biological knowledge of those organisms has not kept pace with sequencing efforts. To unlock this genetic potential there is an urgent need for generic (i.e. non-species specific) genetic toolboxes. Recently, we developed a method, termed chassis-independent recombinase-assisted genome engineering (CRAGE), enabling the integration and expression of large complex gene clusters directly into the chromosomes of diverse bacteria. Here we expand upon this technology by incorporating CRISPR-Cas9 allowing precise genome editing across multiple bacterial species. To do that we have developed a landing pad that carries one wild-type and two mutant lox sites to allow integration of foreign DNA at two locations through Cre-lox recombinase-mediated cassette exchange (RMCE). The first RMCE event is to integrate the Cas9 and the DNA repair protein genes RecET, and the second RMCE event enables the integration of customized sgRNA and a repair template. Following this workflow, we achieved precise genome editing in four different gammaproteobacterial species. We also show that the inserted landing pad and the entire editing machinery can be removed scarlessly after editing. We report here the construction of a single landing pad transposon and demonstrate its functionality across multiple species. The modular design of the landing pad and accessory vectors allows design and assembly of genome editing platforms for other organisms in a similar way. We believe this approach will greatly expand the list of bacteria amenable to genetic manipulation and provides the means to advance our understanding of the microbial world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Are Marine Low Cloud Droplet Concentrations Buffered by Entrained Aitken‐Mode Aerosol (Final Technical Report)

During the summertime, the high-latitude oceans come to life with green phytoplankton, which gain their energy from sunlight and are food for sea creatures small and large. When the phytoplankton are eaten or die, sulfur-rich gases are released under the ocean surface and mix into the air. Observations suggest that, over the Southern Ocean, frequent storms lift this air high into the atmosphere while raining out particulates like salt. As a result, the sulfur-rich air then spawns high concentrations of small ‘Aitken-mode’ aerosol particles. We hypothesize that these particles work their way down into the marine boundary layer, where they can replenish the supply of cloud-condensation nuclei scavenged by frequent precipitation. Further, this process maintains high concentrations of liquid cloud droplets in austral summer, promoting more sunlight to be reflected to space. We call this ‘Aitken buffering’. The primary objective of this project has been to document and test what role Aitken-mode aerosols play in clouds over the Southern Ocean and elsewhere. This effort has included three main components: 1) developing a computer model that realistically simulates the aerosol processes and the small-scale turbulent air motions that move aerosols around and create the clouds, 2) using that model to interpret and extend these observations for process understanding, by allowing different factors that contribute to the aerosol budget, such as surface wind speed, precipitation, surface gas exchange, etc. to be separated, and 3) studying Aitken-mode aerosol and its variability with a focus over the Southern Ocean and Antarctica, using data from Atmospheric Radiation Measurement (ARM) sites and other available observations. Initial computer studies in more idealized conditions found that elevated concentrations of Aitken-mode aerosols above the clouds could help prevent the breakup of those clouds by acting as cloud-condensation nuclei after they were entrained into the cloudy boundary layer. Simulations of a day during the ACE-ENA field campaign showed that Aitken-mode aerosols could also prevent cloud breakup under more realistic conditions. A new method that extracts information about Aitken-mode aerosols from measurements of aerosols onto which cloud droplets can form finds that Aitken-mode aerosols do vary seasonally over the Southern Ocean, with a peak in summertime, as described above. Other work during this project has focused on understanding how patterns of water vapor, clouds and precipitation are coupled within low-lying clouds over the oceans, and also on how cloud droplets cluster within clouds and how the distribution of cloud droplet sizes change as dry air is mixed into clouds, with the latter studies also using observations from ACE-ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation of Relativistic Coupled Cluster Theory for Massively Parallel GPU-Accelerated Computing Architectures

In this paper, we report reimplementation of the core algorithms of relativistic coupled cluster theory aimed at modern heterogeneous high-performance computational infrastructures. The code is designed for parallel execution on many compute nodes with optional GPU coprocessing, accomplished via the new ExaTENSOR back end. The resulting ExaCorr module is primarily intended for calculations of molecules with one or more heavy elements, as relativistic effects on the electronic structure are included from the outset. In the current work, we thereby focus on exact two-component methods and demonstrate the accuracy and performance of the software. The module can be used as a stand-alone program requiring a set of molecular orbital coefficients as the starting point, but it is also interfaced to the DIRAC program that can be used to generate these. We therefore also briefly discuss an improvement of the parallel computing aspects of the relativistic self-consistent field algorithm of the DIRAC program.

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

Enabling Efficient and Accurate Computational Studies of MOF Reactivity via QM/MM and QM/QM Methods

Electronic structure calculations can provide unique insight into metal–organic framework (MOF) reactivity and defect formation. Such calculations can be broadly categorized as utilizing either periodic or cluster models, each with their respective advantages and disadvantages. In the present study, we demonstrate how multiscale methods can leverage the advantages of both approaches to enable high levels of accuracy and computational efficiency in studies of MOF reactivity. Using defect formation in a zeolitic imidazolate framework (ZIF) as a prototypical example, we benchmark a quantum mechanics/molecular mechanics (QM/MM) scheme that enables highly efficient cluster-based calculations on MOFs. We demonstrate the importance of correctly accounting for the influence of both dative bond cleavage in the QM cluster and long-range mechanical coupling to the bulk to achieving accurate QM/MM studies of MOFs. We subsequently leverage these cluster models in a QM/QM scheme that goes beyond standard DFT to yield “gold-standard” correlated wave function results on MOFs at modest computational cost. Crucially, we find several cases in which the incorporation of these correlated corrections yields qualitatively important corrections over conventional DFT values.

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

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗