Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electronic properties of ThC − and ThC

The present study investigates the properties of low-lying states of ThC − and ThC using correlated wave function theories. To this end, we employed multireference calculations and various coupled-cluster approaches in combination with large correlation-consistent basis sets. These methods were applied to examine potential energy curves (PECs), electron configurations, energetics, spectroscopic constants, and spin–orbit coupling effects for 9 states of ThC − and 18 states of ThC. The ground states of ThC − and ThC were identified as single-reference 2 Σ + 1/2 (I) and $^3Σ^+_{0^-} (\textrm{I})$, respectively. Electron detachment from the 7s orbital of the Th center in ThC − [ 2 Σ + (I); 1σ 2 2σ1 3 σ2 1 π 4 ] yields ThC [ 3 Σ + (I); 1σ 2 2σ 1 3σ 1 1π 4 ]. The calculated adiabatic detachment energy (ADE) and vertical detachment energy (VDE) for this process are 1.591 and 1.604 eV, respectively. Furthermore, the dissociation energy (D 0 ) of ThC [$^3Σ^+_{0^-} (\textrm{I})$] is predicted to be 5.099 eV. The standard enthalpy of formation, $ΔH^{°}_{\textrm{f}}$ (298 K), of ThC is estimated to be 822.52 ± 6 kJ mol −1 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Exploring the exact limits of the real-time equation-of-motion coupled cluster cumulant Green’s functions

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green’s function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green’s function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green’s function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green’s function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green’s function approach and other Green’s function approaches in the numerical exact limit.

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

Quantum many-body linear algebra, Hamiltonian moments, and a coupled-cluster inspired framework

Here, we propose a general strategy to develop quantum many-body approximations of primitives in linear algebra algorithms. As a practical example, we introduce a coupled-cluster inspired framework to produce approximate Hamiltonian moments and demonstrate its application in various linear algebra algorithms for ground state estimation. Through numerical examples, we illustrate the difference between the ground-state energies arising from quantum many-body linear algebra and those from the analogous many-body perturbation theory. Our results support the general idea of designing quantum many-body approximations outside of perturbation theory, providing a route to new algorithms and approximations.

Algorithms and data structure↗
NASA NTRS1992

The calculation of the vibrational frequencies of CuCO+, NiCO and CuCH3

The vibrational frequencies of CuCO+, NiCO, and CuCH3 are computed at both the self-consistent-field (SCF) and coupled-cluster singles and doubles, with a perturbational estimate of the triples (CCSD(T)), levels of theory. Overall the SCF frequencies are in qualitative agreement with the CCSD(T) results. The calculations show that the frequencies for the electrostatically bound CuCO+ are quite different from the datively bound NiCO. The CCSD(T) calculations support a previous suggestion that the C-H stretches in CuCH3 have been misassigned.

Sodupe, M.↗
DOE OSTI2025

Repartitioning the Hamiltonian in many-body second-order Brillouin–Wigner perturbation theory: Uncovering new size-consistent models

Second-order Møller-Plesset perturbation theory is well-known as a computationally inexpensive approach to the electron correlation problem that is size-consistent with a size-consistent reference but fails to be regular. On the other hand, the less well-known many-body version of Brillouin-Wigner perturbation theory has the reverse properties: it is regular but fails to be size-consistent when used with the standard MP partitioning. Consequently, its widespread use remains limited. In this work, we analyze the ways in which it is possible to use alternative non-MP partitions of the Hamiltonian to yield variants of BW2 that are size-consistent as well as regular. We show that there is a vast space of such BW2 theories and also show that it is possible to define a repartitioned BW2 theory from the ground state density alone, which regenerates the exact correlation energy. We also provide a general recipe for deriving regular, size-consistent, and size-extensive partitions from physically meaningful components, and we apply the result to small model systems. The scope of these results appears to further set the stage for a revival of BW2 in quantum chemistry.

Ab initio perturbation↗
DOE OSTI2025

Benchmarking third-order cluster perturbation theory for electronically excited states

In this study, we investigate the reliability of cluster perturbation (CP) theory applied to the calculation of electronically excited states through a comprehensive benchmark. In CP theory, perturbative corrections are added to the properties of a parent excitation space, which converge toward the properties of a target excitation space. For the CPS(D-n) model, perturbative corrections through order n are added to the coupled cluster singles (CCS) excitation energies to target the coupled cluster singles and doubles (CCSD) excitation energies. Through a comparative analysis of excitation energy calculations across a diverse set of molecules and wavefunction methods, we present a comprehensive evaluation of the accuracy of the third-order CPS(D) model, CPS(D-3), in calculating excitation energies. Further, our findings demonstrate that CPS(D-3) is a reliable alternative to established methods, particularly CCSD, while systematically overestimating the excitation energies compared to high-level coupled cluster methods such as CC3. These results highlight the strengths and limitations of CPS(D-3), as well as the promising directions for its future development.

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

Structure of odd-mass Ne, Na, and Mg nuclei

The island of inversion is a region of neutron-rich nuclei around 31 Na that are deformed in their ground states. In this region, less is known about the energy levels of odd-mass nuclei, how they evolve with increasing neutron numbers, and how they can be organized into rotational bands. We perform ab initio coupled-cluster calculations of spectra in odd-mass Ne, Na, and Mg nuclei based on an interaction of chiral effective field theory. In conclusion, our results confirm some tentative spin and parity assignments, predict the structure of nuclei near the neutron dripline, and inform us about rotational bands in this region of the nuclear table.

20 ≤ A ≤ 38↗
DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
DOE OSTI2025

Exploring the nexus of many-body theories through neural network techniques: the tangent model

Abstract In this paper, we present a physically informed neural network (NN) representation of the effective interactions associated with coupled-cluster downfolding models to describe chemical systems and processes. The NN representation not only allows us to evaluate the effective interactions efficiently for various geometrical configurations of chemical systems corresponding to various levels of complexity of the underlying wave functions, but also reveals that the bare and effective interactions are related by a tangent function of some latent variables. We refer to this characterization of the effective interaction as a tangent model. We discuss the connection between this tangent model for the effective interaction with the previously developed theoretical analysis that examines the difference between the bare and effective Hamiltonians in the corresponding active spaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

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

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
NASA NTRS1995

The Ozonide Anion: A Theoretical Study

Complete active space self-consistent field (CASSCF) and CASSCF second order perturbation theory (CASPT2) methods have been used for the geometry optimization and calculation of harmonic and fundamental frequencies of the ozonide ion O3(-) and the ozonide lithium complex Li(+)O3(-). For O3(-) harmonic frequencies have also been obtained using the coupled-cluster method, CCSD(T). Infrared intensities are computed from dipole moment derivatives at the CASSCF level. The predicted equilibrium geometry for O3(-) is: Roo = 1.361 Angstroms and angle ooo = 115.4 degrees, and the fundamental frequencies are: nu(sub 1) = 989 per centimeter, nu(sub 2) 556 per centimeter, nu(sub 3) = 870 per centimeter (experimental values are: Roo =1.36 plus or minus 0.02 Angstroms, angle (ooo) = 111.8 plus or minus 2.0 degrees, nu(sub 1) = 975(50) per centimeter, nu(sub 2) = 550(50) per centimeter, nu(sub 3) = 880(50) per centimeter). Corresponding data for the lithium ozonide complex have also been obtained. The presented data contradict the previous interpretation of the IR and Raman spectrum obtained after deposition of ozone in N2, argon, or neon matrices with atomic beams of alkali metals. The presence of the lithium cation raises the asymmetric stretch frequency to about 940 per centimeter, which is contradictory to assumptions made in the assignments of the matrix spectra. Calculations made in a dielectric medium strongly suggest that the effect of the matrix on the IR spectrum is small for O3(-) itself. The dissociation and atomization energies of O3(-) are found to be in agreement with experiment.

Borowski, Poitr↗
NASA NTRS1991

Analytic gradients for coupled-cluster energies that include noniterative connected triple excitations - Application to cis- and trans-HONO

An efficient formulation of the analytic energy gradient for the single and double excitation coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations is presented. The formulation has a small computational cost, and the algebraic manipulations may be applied generally to the analytic gradient of Moller-Plesset perturbation theory energies. The new formulation has been implemented in an efficient set of programs that utilize highly vectorized algorithms and has been used to investigate the equilibrium structures, harmonic vibrational frequencies, IR intensities, and energy separation of cis- and trans-HONO.

Lee, Timothy J.↗
NASA NTRS1994

A Doubles Correction to Electronic Excited States from Configuration Interaction in the Space of Single Substitutions

A perturbative correction to the method of configuration interaction with single substitutions (CIS) is presented. This CIS(D) correction approximately introduces the effect of double substitutions which are absent in CIS excited states. CIS(D) is a second-order perturbation expansion of the coupled-cluster excited state method, restricted to single and double substitutions, in a series in which CIS is zeroth order, and the first-order correction vanishes. CIS (D) excitation energies are size consistent and the calculational complexity scales with the fifth power of molecular size, akin to second-order Moller-Plesset theory for the ground state. Calculations on singlet excited states of ethylene, formaldehyde, acetaldehyde, butadiene and benzene show that CIS (D) is a uniform improvement over CIS. CIS(D) appears to be a promising method for examining excited states of large molecules, where more accurate methods are not feasible.

Head-Gordon, Martin↗
NASA NTRS1994

An Accurate Quartic Force Field and Vibrational Frequencies for HNO and DNO

An accurate ab initio quartic force field for HNO has been determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with the correlation consistent polarized valence triple zeta (cc-pVTZ) basis set. Improved harmonic frequencies were determined with the cc-pVQZ basis set. Fundamental vibrational frequencies were determined using a second-order perturbation theory analysis and also using variational calculations. The N-0 stretch and bending fundamentals are determined well from both vibrational analyses. The H-N stretch, however, is shown to have an unusually large anharmonic correction, and is not well determined using second-order perturbation theory. The H-N fundamental is well determined from the variational calculations, demonstrating the quality of the ab initio quartic force field. The zero-point energy of HNO that should be used in isodesmic reactions is also discussed.

Dateo, Christopher E.↗
NASA NTRS1994

Symmetry breaking in O4(+): An application of the Brueckner coupled-cluster method

A recent calculation of the antisymmetric stretch frequency for the rectangular structure of quartet O4(+) using the singles and doubles quadratic configuration interaction method with a perturbational estimate of connected triple excitations (QCISD(T)) method gave a value of 3710 cm(exp -1). This anomalous frequency is shown to be a consequence of symmetry breaking effects, which occur even though the QCISD(T) solution derived from a delocalized SCF reference function lies energetically well below the two localized (symmetry-broken) solutions at the equilibrium geometry. The symmetry breaking is almost eliminated at the CCSD level of theory, but the small remaining symmetry breaking effects are magnified at the CCSD(T) level of theory so that the antisymmetric stretch frequency is still significantly in error. The use of the Brueckner coupled cluster method, however, leads to a symmetrical solution which is free of symmetry breaking effects, with an antisymmetric stretch frequency of 1322 cm(exp -1), in good agreement with our earlier calculations using the complete active space self consistent field/complete active space state interaction (CASSCF/CASSI) method.

Barnes, Leslie A.↗
NASA NTRS1995

An Accurate ab initio Quartic Force Field and Vibrational Frequencies for CH4 and Isotopomers

A very accurate ab initio quartic force field for CH4 and its isotopomers is presented. The quartic force field was determined with the singles and doubles coupled-cluster procedure that includes a quasiperturbative estimate of the effects of connected triple excitations, CCSD(T), using the correlation consistent polarized valence triple zeta, cc-pVTZ, basis set. Improved quadratic force constants were evaluated with the correlation consistent polarized valence quadruple zeta, cc-pVQZ, basis set. Fundamental vibrational frequencies are determined using second-order perturbation theory anharmonic analyses. All fundamentals of CH4 and isotopomers for which accurate experimental values exist and for which there is not a large Fermi resonance, are predicted to within +/- 6 cm(exp -1). It is thus concluded that our predictions for the harmonic frequencies and the anharmonic constants are the most accurate estimates available. It is also shown that using cubic and quartic force constants determined with the correlation consistent polarized double zeta, cc-pVDZ, basis set in conjunction with the cc-pVQZ quadratic force constants and equilibrium geometry leads to accurate predictions for the fundamental vibrational frequencies of methane, suggesting that this approach may be a viable alternative for larger molecules. Using CCSD(T), core correlation is found to reduce the CH4 r(e), by 0.0015 A. Our best estimate for r, is 1.0862 +/- 0.0005 A.

Lee, Timothy J.↗
NASA NTRS1990

The vibrational frequencies of ozone

The singles and doubles coupled-cluster (CCSD) method is used, together with the recently proposed CCSD(T) and CCSDT-1 methods, to determine the equilibrium structure and harmonic frequencies of ozone at various levels, on the basis of large one-particle basis sets. It is demonstrated that f-type functions are required in order to obtain quantitative accuracy for the equilibrium structure and harmonic frequencies. The CCSD(T) level of theory, in conjunction with a (4s 3p 2d 1f) atomic natural orbital set, yields an equilibrium geometry and harmonic frequencies that are in excellent agreement with experimental results; an effort has accordingly been made to determine the cubic and quartic force constants, and evaluate the anharmonic corrections to the fundamental vibrational frequencies.

Lee, Timothy J.↗