Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled cluster 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

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New reactions of diazene and related species for modelling combustion of amine fuels

Potential energy surfaces for reactions involving N 2 H 2 isomers of diazene (diimide) have been explored using density functional theory, with energies based on coupled-cluster theory. A focus is on processes that create or consume these species, and isomerisation between the E (trans) and Z (cis) forms of HNNH. Furthermore, these include isomerisation and dissociation pathways for HNNH, addition of H atoms to form N 2 H 3 , abstraction by H atoms yielding short-lived NNH, and abstraction reactions of H with N 2 H 3 . Transition state and capture theories are applied for high-pressure-limiting behaviour, while low-pressure and falloff regions are characterised via the methods of Troe and coworkers. Rate constants and thermochemistry are provided to improve models of diamine chemistry, relevant to the combustion of NH 3 especially at high concentrations, high pressures or under reducing conditions. Results indicate that amine radical recombination mainly yields the E HNNH isomer, while H-abstraction from N 2 H 3 results in E HNNH and H 2 NN. However, at elevated temperature E → Z isomerisation becomes competitive, and Z HNNH, being more reactive, acts to enhance the diazene consumption rate.

10 SYNTHETIC FUELS↗

Charge Radius of the Short-Lived Ni 68 and Correlation with the Dipole Polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus 68 Ni at the N = 40 subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α$ D has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus 48 Ca. Three-particle–three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of 68 Ni.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^(8) determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

Quantum Computing↗

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^(8) determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

Quantum Computing↗

Coupled-cluster techniques for computational chemistry: The CFOUR program package

An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Lastly, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined in the appendix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal trend of charge radii of even-even Ca–Zn nuclei

Radii of nuclear charge distributions carry information about the strong and electromagnetic forces acting inside the atomic nucleus. Whereas the global behavior of nuclear charge radii is governed by the bulk properties of nuclear matter, their local trends are affected by quantum motion of proton and neutron nuclear constituents. The measured differential charge radii δ$\langle$$r^{2}_{c}$$\rangle$ between neutron numbers N = 28 and N = 40 exhibit a universal pattern as a function of n = N – 28 that is independent of the atomic number. In this work, we analyze this remarkable behavior in even-even nuclei from calcium to zinc using two state-of-the-art theories based on quantified nuclear interactions: the ab initio coupled cluster theory and nuclear density functional theory. Both theories reproduce the smooth rise of differential charge radii and their weak dependence on the atomic number. By considering a large set of isotopic chains, we show that this trend can be captured by just two parameters: the slope and curvature of δ$\langle$$r^{2}_{c}$$\rangle$(n). We demonstrate that these parameters show appreciable model dependence, and the statistical analysis indicates that they are not correlated with any single model property, i.e., they are impacted by both bulk nuclear properties as well as shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

From closed shells to open shells: Coupled-cluster calculations of atomic nuclei

Coupled-cluster theory is a powerful tool for first-principles calculations of atomic nuclei, enabling accurate predictions of nuclear observables across the Segrè chart. While coupled-cluster computations are especially efficient at shell closures, extensions have been developed to tackle open-shell nuclei, by exploiting the equation-of-motion method or by expanding the coupled-cluster wave function on top of a symmetry-breaking (either deformed or superfluid) reference state. In this study, we provide a comprehensive comparison of these different formulations applied to the calcium and nickel isotopes using nuclear two-and three-body interactions from chiral effective field theory. Here, based on ground-state energies, two-neutron separation energies, and two-neutron shell gaps, different coupled-cluster computations—based on symmetry-broken reference states and equationof-motion techniques— offer consistent descriptions of bulk properties across medium-mass isotopic chains.

Marino, Francesco [Johannes Gutenberg-Universität ↗

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

Here, we investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core–valence separated (FC-CVS) EOM-CCSD approach – a successful variant of EOM-CCSD tailored towards core excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Schiff Moments and CP Violation

This article reviews the calculation of nuclear Schiff moments, which one must know in order to interpret experiments that search for time-reversal-violating electric dipole moments in certain atoms and molecules. After briefly reviewing the connection between dipole moments and CP violation in and beyond the Standard Model of particle physics; Schiff's theorem, which concerns the screening of nuclear electric dipole moments by electrons; Schiff moments; and experiments to measure dipole moments in atoms and molecules, this review examines attempts to compute Schiff moments in nuclei such as 199 Hg and octupole-deformed isotopes such as 225 Ra, which are particularly useful in experiments. It then turns to ab initio nuclear-structure theory, describing ways in which both the in-medium similarity renormalization group and coupled-cluster theory can be used to compute important Schiff moments more accurately than the less controlled methods that have been applied so far.

CP violation↗

An open-shell restricted Hartree-Fock perturbation theory based on symmetric spin orbitals

A new open-shell perturbation theory is formulated in terms of symmetric spin orbitals. Only one set of spatial orbitals is required, thereby reducing the number of independent coefficients in the perturbed wavefunctions. For second order, the computational cost is shown to be similar to a closed-shell calculation. This formalism is therefore more efficient than the recently developed RMP, ROMP or RMP-MBPT theories. The perturbation theory described herein was designed to have a close correspondence with our recently proposed coupled-cluster theory based on symmetric spin orbitals. The first-order wavefunction contains contributions from only doubly excited determinants. Equilibrium structures and vibrational frequencies determined from second-order perturbation theory are presented for OH, NH, CH, 02, NH2 and CH2.

Lee, Timothy J.↗

Ultrafast infrared transient absorption spectroscopy of gas-phase Ni(CO) 4 photodissociation at 261 nm

Here, we employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin–orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Rise of Neural Networks for Materials and Chemical Dynamics

Machine learning (ML) is quickly becoming a premier tool for modeling chemical processes and materials. ML-based force fields, trained on large data sets of high-quality electron structure calculations, are particularly attractive due their unique combination of computational efficiency and physical accuracy. This Perspective summarizes some recent advances in the development of neural network-based interatomic potentials. Designing high-quality training data sets is crucial to overall model accuracy. One strategy is active learning, in which new data are automatically collected for atomic configurations that produce large ML uncertainties. Another strategy is to use the highest levels of quantum theory possible. Transfer learning allows training to a data set of mixed fidelity. A model initially trained to a large data set of density functional theory calculations can be significantly improved by retraining to a relatively small data set of expensive coupled cluster theory calculations. These advances are exemplified by applications to molecules and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic isotropic hyperfine properties for first row elements (B–F) revisited

Benchmark quality isotropic hyperfine properties have been obtained for first row elements (B–F) using a systematic composite approach consisting of a sequence of core/valence correlation consistent basis sets, up through aug-cc-pCV8Z, along with configuration interaction and coupled cluster theory methods. The best nonrelativistic final values (in MHz) are 10.64 (B), 20.22 (C), 10.59 (N), –31.74 (O), and 318.30 (F) and are in very good agreement with available experimental values for these difficult-to-describe properties. Agreement is especially close in the case of N, which has the most accurate experimental value. The spin densities derived from the best composite level of theory were found to closely follow a simple quadratic scaling with the atomic number, Z. Finally, observed convergence rates in the 1-particle and n-particle expansions obtained here may be useful in judging likely accuracy that can be expected in studies of molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗