Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled cluster”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hilbert space multireference coupled cluster tailored by matrix product states

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants—the state universal one and the Brillouin–Wigner and Mukherjee’s state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO–LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Chemistry↗
DOE OSTI2020

Toward DMRG-tailored coupled cluster method in the 4c-relativistic domain

There are three essential problems in computational relativistic chemistry: Electrons moving at relativistic speeds, close lying states, and dynamical correlation. Currently available quantum-chemical methods are capable of solving systems with one or two of these issues. However, there is a significant class of molecules in which all the three effects are present. These are the heavier transition metal compounds, lanthanides, and actinides with open d or f shells. For such systems, sufficiently accurate numerical methods are not available, which hinders the application of theoretical chemistry in this field. In this paper, we combine two numerical methods in order to address this challenging class of molecules. These are the relativistic versions of coupled cluster methods and the density matrix renormalization group (DMRG) method. To the best of our knowledge, this is the first relativistic implementation of the coupled cluster method externally corrected by DMRG. The method brings a significant reduction of computational costs as we demonstrate on the system of TlH, AsH, and SbH.

Brandejs, Jan (ORCID:0000000221073095)↗
DOE OSTI2020

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

Coupled Cluster Study of the Interactions of AnO 2 , AnO 2 + , and AnO 2 2+ (An = U, Np) with N 2 and CO

We report thermochemical and spectroscopic properties for actinyl complexes involving UO 2 2+/1+/0 and NpO 2 2+/1+/0 with N 2 and CO, together with the UO 2 -O 2 , UO 2 + -O 2 , and UO 2 + -NO complexes, have been studied for the first time using an accurate composite coupled cluster approach. Two general bonding motifs were investigated, end-on (η 1 ) and side-on (η 2 ) relative to the metal center of the actinyls. For end-on CO complexes, both C-coordinated (An-C) and O-coordinated (An-O) structures were examined, with the former always being lower in energy. All of the η 1 complexes were calculated to be stable, with dissociation energies ranging from 2 to 36 kcal/mol, except for that of UO 2 + -O 2 (the η 1 orientation for UO 2 + -NO was not amenable to single reference coupled cluster). In agreement with a previous study, the η 2 structure for UO 2 + -O 2 was calculated to be relatively strongly bound, by 22.3 kcal/mol in this work. The closely related NO complex, however, had a calculated dissociation energy of just 4.0 kcal/mol. The binding energy of O 2 to neutral UO 2 in a η 2 orientation was calculated to be very strong, 75.4 kcal/mol, and strongly resembled a UO 2 + (O 2 – ) complex at equilibrium. The N–N and C–O bonds were found to be somewhat activated for all the side-on (η 2 ) neutral An(IV) complexes, with stretching frequencies of N 2 or CO being red-shifted by as much as 480 cm –1 with a 0.06 Å bond length elongation. Dissociation energies for the η 1 complexes are strongly correlated with the extent of electron transfer from ligand to actinyl. The nature of bonding in the actinyl complexes is examined using natural resonance theory (NRT). The correlation between bonding motif and small molecule activation is in agreement with experiments in condensed phases.

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

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

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

Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems

In this work, we discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas (UEG) at warm, dense conditions, and some simple materials. The performance of model systems at high temperatures is encouraging: for the one-dimensional Hubbard model, FT-CCSD provides a qualitatively accurate description of finite-temperature correlation effects even at U = 8, and it allows for the computation of systematically improvable exchange–correlation energies of the warm, dense UEG over a wide range of conditions. We highlight the obstacles that remain in using the method for realistic ab initio calculations on materials

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

Equation-of-motion internally contracted multireference unitary coupled-cluster theory

The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version

An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin–orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S o → S 1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. Here, a UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.

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

High-level coupled-cluster energetics by merging moment expansions with selected configuration interaction

Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics, we propose a novel form of the CC(P; Q) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the perturbative selection made iteratively (CIPSI) algorithm. The advantages of the resulting CIPSI-driven CC(P; Q) methodology are illustrated by a few molecular examples, including the dissociation of F 2 and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.

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

Spectral function for 4 He using the Chebyshev expansion in coupled-cluster theory

Here, we compute spectral function for 4 He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows us to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark our predictions with other calculations in the literature and with electron-scattering data in the quasi-elastic peak. The spectral function formalism allows one to extend ab initio lepton-nucleus cross sections into the relativistic regime. This makes it a promising tool for modeling this process at higher-energy transfers. The results we present open the door for studies of heavier nuclei, important for the neutrino oscillation programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Coupled-Cluster Calculations of Neutrinoless Double- β Decay in Ca 48

We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-$\beta$ decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-$\beta$ decay of $^{48}$Ca with a quenching factor deduced from two-body currents in recent ab initio calculation of the Ikeda sum rule in $^{48}$Ca [Gysbers et al., Nat. Phys. 15, 428 (2019)].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

On the functional dependence of transition-potential coupled cluster

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

Coupled-cluster methods↗
DOE OSTI2025

Relativistic core–valence-separated equation-of-motion coupled-cluster singles and doubles method: Efficient implementation and benchmark calculations

An efficient implementation for the relativistic exact two-component core–valence-separated equation-of-motion coupled-cluster singles and doubles (X2C-CVS-EOM-CCSD) method is reported. The explicit exclusion of pure valence excitations in the EOM-CCSD excited-state eigenvalue equations significantly improves the efficiency for calculations of core-excited states. Benchmark relativistic CVS-EOM-CC calculations with systematic inclusion of relativistic, correlation, and basis-set effects are shown to provide highly accurate results for core ionized and excited states involving heavy atoms.

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

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian

Here simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.

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

Analysis of Fourth-, Fifth-, and Infinite-Order Triple Excitations in Unitary Coupled Cluster Theory

Here, in this work, we introduce a perturbative correction to the unitary coupled cluster method with single and double excitations (UCCSD) that incorporates the effects of missing triple excitations through fifth-order in many-body perturbation theory (MBPT). Referred to as UCCSD[T-5], this method is benchmarked alongside the previously developed UCCSD[T] to lend insight into the behavior of perturbative triples corrections relative to UCCSDT, which inherently provides an infinite-order treatment of triple excitations, as well as full configuration interaction (FCI). Two key findings emerge from this analysis. First, UCCSD[T] consistently yields ground-state energies in closest agreement with FCI, outperforming both UCCSD[T-5] and UCCSDT. Second, UCCSD[T-5] largely emulates the behavior of UCCSDT, suggesting a close relationship between fifth- and infinite-order triple-excitation contributions. Given the growing interest in UCC ansätze for quantum computing and the limitations of current quantum hardware, these results highlight the potential of classically computed perturbative corrections within UCC theory to capture triple-excitation effects without the additional quantum resources required by the UCCSDT ansatz.

Windom, Zachary W. [Oak Ridge National Laboratory ↗