Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled-cluster methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Coupled-cluster based basis sets for valence correlation calculations

Novel basis sets are generated that target the description of valence correlation in atoms H through Ar. The new contraction coefficients are obtained according to the Atomic Natural Orbital (ANO) procedure from CCSD(T) (coupled-cluster singles and doubles with perturbative triples correction) density matrices starting from the primitive functions of Dunning et al. [J. Chem. Phys. 90, 1007 (1989); ibid. 98, 1358 (1993); ibid. 100, 2975 (1993)] (correlation consistent polarized valence X-tuple zeta, cc-pVXZ). The exponents of the primitive Gaussian functions are subject to uniform scaling in order to ensure satisfaction of the virial theorem for the corresponding atoms. These new sets, named ANO-VT-XZ (Atomic Natural Orbital Virial Theorem X-tuple Zeta), have the same number of contracted functions as their cc-pVXZ counterparts in each subshell. The performance of these basis sets is assessed by the evaluation of the contraction errors in four distinct computations: correlation energies in atoms, probing the density in different regions of space via 〈r{sup n}〉 (−3 ≤ n ≤ 3) in atoms, correlation energies in diatomic molecules, and the quality of fitting potential energy curves as measured by spectroscopic constants. All energy calculations with ANO-VT-QZ have contraction errors within “chemical accuracy” of 1 kcal/mol, which is not true for cc-pVQZ, suggesting some improvement compared to the correlation consistent series of Dunning and co-workers.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

Theoretical studies of the global minima and polarizabilities of small lithium clusters

Lithium clusters Lin (n=1-20) have been investigated with density functional theory (DFT) and coupled—cluster (CC) methods. The global-minimum structures are located via an improved basin---hopping algorithm and the lowest energy Lin isomers are confirmed with DFT geometry optimizations, CCSD(T) energy calculations, and by comparing simulated and experimental polarizabilities. The tetrahedral Li4 structure is found to be the basic building block of lithium clusters Lin (n=6-20). Simulated polarizabilities, including thermal effects at room temperature, are in good agreement with measured isotropic polarizabilities.

clusters↗
DOE OSTI2018

Green’s Function Coupled-Cluster Approach: Simulating Photoelectron Spectra for Realistic Molecular Systems

In this paper we present an efficient implementation for the analytical energy- dependent coupled cluster Green’s function with singles and doubles (GFCCSD) approach with our first practice being computing spectral functions of real molecular systems. Due to its algebraic structure, the presented method is highly scalable, and is capable of computing spectral function for a given molecular system in any energy region. Several typical examples have been given to demonstrate its capability of computing spectral functions not only in the valence band, but also in the core-level energy region. Satellite peaks have been observed in the inner valence band and core-level energy region where many-body effect becomes significant and single particle picture of ionization often breaks down. The accuracy test has been carried out by extensively comparing our GFCCSD results with experimental results as well as the theoretical results from other methods. It turns out GFCCSD method is able to provide qualitative or semi-quantitative level of description of ionization processes in both the core and valence regimes. To significantly improve the GFCCSD results for the main states, larger basis set can usually be employed, whereas the improvement of the GFCCSD results for the satellite states needs higher order many-body terms to be included in the GFCC implementation.

electronic structure theory↗
DOE OSTI2018

Generalized Unitary Coupled Cluster Wave functions for Quantum Computation

We introduce a unitary coupled-cluster (UCC) ansatz termed k-UpCCGSD that is based on a family of sparse generalized doubles operators, which provides an affordable and systematically improvable unitary coupled-cluster wave function suitable for implementation on a near-term quantum computer. k-UpCCGSD employs k products of the exponential of pair coupled-cluster double excitation operators (pCCD), together with generalized single excitation operators. We compare its performance in both efficiency of implementation and accuracy with that of the generalized UCC ansatz employing the full generalized single and double excitation operators (UCCGSD), as well as with the standard ansatz employing only single and double excitations (UCCSD). k-UpCCGSD is found to show the best scaling for quantum computing applications, requiring a circuit depth of O(kN) , compared with O(N 3 ) for UCCGSD, and O((N-η) 2 η) for UCCSD, where N is the number of spin orbitals and η is the number of electrons. We analyzed the accuracy of these three ansätze by making classical benchmark calculations on the ground state and the first excited state of H 4 (STO-3G, 6-31G), H 2 O (STO-3G), and N 2 (STO-3G), making additional comparisons to conventional coupled cluster methods. The results for ground states show that k-UpCCGSD offers a good trade-off between accuracy and cost, achieving chemical accuracy for lower cost of implementation on quantum computers than both UCCGSD and UCCSD. UCCGSD is also found to be more accurate than UCCSD but at a greater cost for implementation. Excited states are calculated with an orthogonally constrained variational quantum eigensolver approach. This is seen to generally yield less accurate energies than for the corresponding ground states. Lastly, we demonstrate that using a specialized multideterminantal reference state constructed from classical linear response calculations allows these excited state energetics to be improved.

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

“Best” Iterative Coupled-Cluster Triples Model? More Evidence for 3CC

To follow up on the unexpectedly good performance of several coupled-cluster models with approximate inclusion of 3-body clusters we performed a more complete assessment of the 3CC method for accurate computational thermochemistry in the standard HEAT framework. New spin-integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-ζ basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean absolute errors in cc-pVDZ {3CC, CCSDT, and CCSD(T)} electronic (per valence electron) and atomization energies relative to the CCSDTQ reference for the HEAT data set, were {24, 70, 122} μE h /e and {0.46, 2.00, 2.58} kJ/mol, respectively. The mean absolute errors in the complete-basis-set limit {3CC, CCSDT, and CCSD(T)} atomization energies relative to the HEAT model reference, were {0.52, 2.00, and 1.07} kJ/mol, The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post- CCSD(T) thermochemistry applications, as well as the preferred alternative to CCSDT in general.

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

Probing electronic wave functions of sodium-doped clusters: Dyson orbitals, anisotropy parameters, and ionization cross-sections

We apply high-level ab initio methods to describe the electronic structure of small clusters of ammonia and dimethylether (DME) doped with sodium, which provide a model for solvated electrons. We investigate the effect of the solvent and cluster size on the electronic states. We consider both energies and properties, with a focus on the shape of the electronic wave function and the related experimental observables such as photoelectron angular distributions. The central quantity in modeling photoionization experiments is the Dyson orbital, which describes the difference between the initial N-electron and final (N-1)-electron states of a system. Dyson orbitals enter the expression of the photoelectron matrix element, which determines total and partial photoionization cross-sections. We compute Dyson orbitals for the Na(NH3)n and Na(DME)m clusters using correlated wave functions (obtained with equation-of-motion coupled-cluster model for electron attachment with single and double substitutions) and compare them with more approximate Hartree-Fock and Kohn-Sham orbitals. As a result, we also analyze the effect of correlation and basis sets on the shapes of Dyson orbitals and the experimental observables.

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

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↗
DOE OSTI2016

Coupled cluster Green function: Model involving single and double excitations

In this paper, we report on the development of a parallel implementation of the coupled-cluster (CC) Green function formulation (GFCC) employing single and double excitations in the cluster operator (GFCCSD). A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω). The detailed description of the underlying algorithm is provided, including approximations used that preserve the pole structure of the full GFCCSD method, thereby reducing the computational costs while maintaining an accurate character of methodology. Furthermore, for systems with strong local correlation, our formulation reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally. Here, to demonstrate the accuracy of our approach, several examples including calculations of ionization potentials for benchmark systems are presented and compared against experiment.

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

Equation of motion coupled-cluster cumulant approach for intrinsic losses in x-ray spectra

We present an equation of motion coupled cluster approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS). The method is based on a factorization of the transition amplitude in the time-domain, which leads to a convolution of an effective one-body spectrum and the core-hole spectral function. The spectral function characterizes these losses in terms of shake-up excitations and satellites, and is calculated using a cumulant representation of the core-hole Green’s function that includes non-linear corrections. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.

Rehr, John J.↗
DOE OSTI2022

The many-body expansion for metals. I. The alkaline earth metals Be, Mg, and Ca

We examine the many-body expansion (MBE) for alkaline earth metal clusters, Be n , Mg n , Ca n ( n = 4, 5, 6), at the Møller–Plesset second order perturbation theory, coupled-cluster singles and doubles with perturbative triples, multi-reference perturbation theory, and multi-reference configuration interaction levels of theory. The magnitude of each term in the MBE is evaluated for several geometrical configurations. We find that the behavior of the MBE for these clusters depends strongly on the geometrical arrangement and, to a lesser extent, on the level of theory used. Another factor that affects the MBE is the in situ (ground or excited) electronic state of the individual atoms in the cluster. For most geometries, the three-body term is the largest, followed by a steady decrease in absolute energy for subsequent terms. Though these systems exhibit non-negligible multi-reference effects, there was little qualitative difference in the MBE when employing single vs multi-reference methods. Useful insights into the connectivity and stability of these clusters have been drawn from the respective potential energy surfaces and quasi-atomic orbitals for the various dimers, trimers, and tetramers. Through these analyses, we investigate the similarities and differences in the binding energies of different-sized clusters for these metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2019

Application of the coupled-cluster CC( P ; Q ) approaches to the magnesium dimer

Several coupled-cluster approaches with an approximate treatment of connected triply or quadruply excited clusters, including the CCSD(T), CCSD(2) T , and CR-CC(2,3) corrections to CCSD, the active-space CCSDt and CCSDTq methods, and the CC(t;3) and CC(q;4) corrections to CCSDt and CCSDTq derived using the CC(P;Q) framework, and their CCSDT and CCSDTQ parents are applied to the ground-state potential curve and vibrational term values of the magnesium dimer. The correlation-consistent aug-cc-pV(n+d)Z and aug-cc-pwCVnZ (n = T and Q) basis sets are used. Among the noniterative triples corrections to CCSD, the CR-CC(2,3) approach performs the best, but, in analogy to the previously studied beryllium dimer [I. Magoulas et al., J. Phys. Chem. A 122, 1350 (2018)], the CC(t;3) and CC(q;4) methods outperform other employed approaches in reproducing the CCSDT and CCSDTQ data. Composite calculations combining the nearly all-electron CCSDT and valence CC(q;4) and CCSDTQ computations reproduce the experimentally resolved part of the vibrational spectrum of Mg2 to within ~1 cm -1 when the vibrational spacings are examined, while predicting five extra vibrational states near the dissociation threshold, which have not been experimentally resolved. The corresponding binding energies and equilibrium bond lengths are within a few cm -1 and 0.01 A from their experimentally derived values.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
DOE OSTI2025

Direct reactions with the AT-TPC

Direct reactions are crucial tools for accessing properties of the atomic nucleus. Fundamental and exotic phenomena such as collective modes, pairing, weakbinding effects and evolution of single-particles energies can be investigated in peripheral collisions between a heavy nucleus and a light target. The necessity of using inverse kinematics to reveal how these structural properties change with isospin imbalance renders direct reactions a challenging technique when using the missing mass method. In this scenario, Active Target Time Projection Chambers (AT-TPC) have demonstrated an outstanding performance in enabling these types of reactions even under conditions of very low beam intensities. The AT-TPC of the Facility for Rare Isotope Beams (FRIB) is a next generation multipurpose Active Target. When operated inside a solenoidal magnet, direct reactions benefit from the measurement of the magnetic rigidity that enables particle identification and the determination of the excitation energy with high resolution without the need of auxiliary detectors. Additionally, the AT-TPC can be coupled to a magnetic spectrometer improving even further its spectroscopic investigation capability. In this contribution, we discuss inelastic scattering and transfer reaction data obtained via the AT-TPC and compare them to theory. In particular, we present the results for the 14 C(p,p′) and 12 Be (p,d) 11 Be reactions. For 14 C, we compare the experimental excitation energy of the first 1 – excited state with coupled-cluster calculationsbased on nuclear interactions from chiral effective field theory and with available shell-model predictions. For 12 Be, we determine the theoretical spectroscopic factors of the 12 Be (p,d) 11 Be transfer reaction in the shell modeland compare them to the experimental excitation spectrum from a qualitative standpoint.

active target↗
NASA NTRS1990

Vibrations in small Mg clusters

The equilibrium geometries, binding energies, and harmonic frequencies of Mg3 and Mg4 have been determined using large atomic natural orbital basis sets in conjunction with high levels of electron correlation. The correlation treatments comprise multireference configuration interaction (MRCI), singles and doubles coupled-cluster (CCSD) theory and the CCSD(T) extension that includes a perturbational estimate of connected triple excitations. As with previous studies of small Be clusters, the CCSD(T) method is found to reproduce the MRCI results with a remarkable degree of accuracy and at a fraction of the computational cost. Using the CCSD(T) method, full quartic force fields for Mg3 and Mg4 have been determined and anharmonic analyses have been performed using second-order perturbation theory. Vibration-rotation interaction constants and centrifugal distortion constants have also been determined.

Lee, Timothy J.↗