Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multireference 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Coupled-cluster theory for strong entanglement in nuclei

Atomic nuclei can exhibit shape coexistence and multireference physics that enters in their ground states, and accurately capturing the ensuing correlations and entanglement is challenging. Here, we address this problem by applying single-reference coupled-cluster theory based on spherical and deformed reference states and the tailored coupled-cluster method. The latter combines configuration interactions to capture static correlations with coupled-cluster theory for dynamic correlations. We compute the atomic nuclei 12 C, 28 Si, and 56 Ni and find that the tailored coupled-cluster method and the single-reference approach based on a deformed Hartree-Fock state yield the most accurate results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Ab Initio Kinetics of Hydrogen Abstraction from Methyl Acetate by Hydrogen, Methyl, Oxygen, Hydroxyl, and Hydroperoxy Radicals

The kinetics of hydrogen abstraction by five radicals (H, O( 3 P), OH, CH 3 , and HO 2 ) from methyl acetate (MA) is investigated theoretically in order to gain further understanding of certain aspects of the combustion chemistry of biodiesels, such as the effect of the ester moiety. Here, we employ ab initio quantum chemistry methods, coupled cluster singles and doubles with perturbative triples correction (CCSD(T)) and multireference averaged coupled pair functional theory (MRACPF2), to predict chemically accurate reaction energetics. Overall, MRACPF2 predicts slightly higher barrier heights than CCSD(T) for MA + H/CH 3 /O/OH, but slightly lower barrier heights for hydrogen abstraction by HO 2 . Based on the obtained reaction energies, we also report high-pressure-limit rate constants using transition state theory (TST) in conjunction with the separable-hindered-rotor approximation, the variable reaction coordinate TST, and the multi-structure all-structure approach. The fitted modified Arrhenius expressions are provided over a temperature range of 250 to 2000 K. The predictions are in good agreement with available experimental results. Abstractions from both of the methyl groups in MA are expected to contribute to consumption of the fuel as they exhibit similar rate coefficients. Moreover, the reactions involving the OH radical are predicted to have the highest rates among the five abstracting radicals, while those initiated by HO 2 are expected to be the lowest.

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

An ab initio investigation of the ground and low-lying singlet and triplet electronic states of XNO{sub 2} and XONO (X = Cl, Br, and I)

A systematic ab initio treatment of the nitryl halides (XNO{sub 2}) and the cis- and trans- conformers of the halide nitrites (XONO), where X = Cl, Br, and I, have been carried out using highly correlated methods with sequences of correlation consistent basis sets. Equilibrium geometries and harmonic frequencies have been accurately calculated in all cases at the explicitly correlated CCSD(T)-F12b level of theory, including the effects of core-valence correlation for the former. Where experimental values are available for the equilibrium structures (ClNO{sub 2} and BrNO{sub 2}), the present calculations are in excellent agreement; however, the X-O distances are slightly too long by about 0.01 Å due to missing multireference effects. Accurate predictions for the iodine species are made for the first time. The vertical electronic excitation spectra have been calculated using equation-of-motion coupled cluster methods for the low-lying singlet states and multireference configuration interaction for both singlet and triplet states. The latter also included the effects of spin-orbit coupling to provide oscillator strengths for the ground state singlet to excited triplet transitions. While for ClNO{sub 2} the transitions to excited singlet states all occur at wavelengths shorter than 310 nm, there is one longer wavelength singlet transition in BrNO{sub 2} and two in the case of INO{sub 2}. The long wavelength tail in the XNO{sub 2} species is predicted to be dominated by transitions to triplet states. In addition to red-shifting from X = Cl to I, the triplet transitions also increase in oscillator strength, becoming comparable to many of the singlet transitions in the case of INO{sub 2}. Hence in particular, the latter species should be very photolabile. Similar trends are observed and reported for the halogen nitrites, many of which for the first time.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium–Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D e and bond length r e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15–0.19. For the single-reference methods, Kohn–Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D e and r e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core–valence terms; this results in highly accurate dissociation energies and good bond lengths. In conclusion, averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

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

Ground and excited state properties of ThBe and AcBe

In this work, the ground and excited states of ThBe and AcBe were investigated by performing high-level multireference and single-reference coupled-cluster quantum chemical calculations with large correlation consistent basis sets. Full potential energy curves (PECs), chemical bonding patterns, energetics, spectroscopic parameters (T e , r e , ω e , and ω e x e ), and spin–orbit effects of 13 and 8 electronic states of ThBe and AcBe, respectively, are reported. The ground electronic states of ThBe and AcBe are single-reference 1 3 Σ − (1σ 2 2σ 2 1π 2 ) and 1 2 Π (1σ 2 2σ 2 1π 1 ), respectively, and originate from their corresponding ground state fragments. The chemical bonding of ThBe (1 3 Σ − ) and AcBe (1 2 Π) are π-dative in character and are formed by d-electron transfers from Th/Ac to the empty 2p x and 2p y of the Be atom. The electron populations of the f-orbitals of both ThBe (1 3 Σ − ) and AcBe (1 2 Π) are minor which exhibit their “transition-metal-like” nature. The estimated bond energies of the spin–orbit ground states of ThBe (1 3 Σ − 0+ ) and AcBe (1 2 Π 1/2 ) are 12.79 and 11.02 kcal mol −1 , respectively. Finally, the bond energy of ThBe was used to estimate its heat of formation ΔH 0 f (298 K) of 869.61 ± 6 kJ mol −1 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

An ab initio study of the electronic structure of the boron oxide neutral (BO), cationic (BO{sup +}), and anionic (BO{sup −}) species

The BO neutral, cationic, and anionic molecular species have been painstakingly studied through multireference configuration interaction and single reference coupled cluster methods employing basis sets of quintuple cardinality. Potential energy curves have been constructed for 38 (BO), 37 (BO{sup +}), and 12 (BO{sup −}) states and the usual molecular parameters have been extracted most of which are in very good agreement with the scarce experimental data. Numerous avoided crossings appear on more or less all of the studied states of the neutral and cationic species challenging the validity of the Born Oppenheimer approximation. Finally, all excited states of the anionic system lie above the ground state of the neutral BO system and are therefore resonances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

On the Reaction of FNO2 with CH3, t-butyl, and C13H21

Theoretical studies are reported for the reaction of FN02 with the radicals CH3, t-butyl, and C13H21, which are templates for the radical site of a hydrogenated diamond surface. All structures axe fully optimized using density functional theory (DFT) based on the B3LYP functional. Calibration calculations axe performed for CH3 + FNO2 using the coupled cluster approach, the the internally contracted multireference configuration. interaction method, and second order perturbation theory based upon complete active space SCF (CASSCF) reference wave function. These calibration calculations support the B3LYP approach for the calculation of bond energies, but show the B3LYP barrier is too low. Combining the calibration calculations with the larger clusters yields our best estimate of a barrier of about 10 kcal/mol. for the reaction of FNO2 with a radical site on hydrogenated diamond.

Thuemmel, H. T.↗
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.↗
DOE OSTI2017

Merging symmetry projection methods with coupled cluster theory: Lessons from the Lipkin model Hamiltonian

Coupled cluster and symmetry projected Hartree-Fock are two central paradigms in electronic structure theory. However, they are very different. Single reference coupled cluster is highly successful for treating weakly correlated systems but fails under strong correlation unless one sacrifices good quantum numbers and works with broken-symmetry wave functions, which is unphysical for finite systems. Symmetry projection is effective for the treatment of strong correlation at the mean-field level through multireference non-orthogonal configuration interaction wavefunctions, but unlike coupled cluster, it is neither size extensive nor ideal for treating dynamic correlation. We here examine different scenarios for merging these two dissimilar theories. We carry out this exercise over the integrable Lipkin model Hamiltonian, which despite its simplicity, encompasses non-trivial physics for degenerate systems and can be solved via diagonalization for a very large number of particles. We show how symmetry projection and coupled cluster doubles individually fail in different correlation limits, whereas models that merge these two theories are highly successful over the entire phase diagram. In conclusion, despite the simplicity of the Lipkin Hamiltonian, the lessons learned in this work will be useful for building an ab initio symmetry projected coupled cluster theory that we expect to be accurate in the weakly and strongly correlated limits, as well as the recoupling regime.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Merging symmetry projection methods with coupled cluster theory: Lessons from the Lipkin model Hamiltonian

Coupled cluster and symmetry projected Hartree-Fock are two central paradigms in electronic structure theory. However, they are very different. Single reference coupled cluster is highly successful for treating weakly correlated systems, but fails under strong correlation unless one sacrifices good quantum numbers and works with broken-symmetry wave functions, which is unphysical for finite systems. Symmetry projection is effective for the treatment of strong correlation at the mean-field level through multireference non-orthogonal configuration interaction wavefunctions, but unlike coupled cluster, it is neither size extensive nor ideal for treating dynamic correlation. We here examine different scenarios for merging these two dissimilar theories. We carry out this exercise over the integrable Lipkin model Hamiltonian, which despite its simplicity, encompasses non-trivial physics for degenerate systems and can be solved via diagonalization for a very large number of particles. We show how symmetry projection and coupled cluster doubles individually fail in different correlation limits, whereas models that merge these two theories are highly successful over the entire phase diagram. Despite the simplicity of the Lipkin Hamiltonian, the lessons learned in this work will be useful for building an ab initio symmetry projected coupled cluster theory that we expect to be accurate in the weakly and strongly correlated limits, as well as the recoupling regime.

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

Reinterpretation of the electronic absorption spectrum of the methylene amidogen radical (H 2 CN)

The odd electronic absorption spectrum of H 2 CN has been of great interest to experiment. Herein, this system is reported on extensively by applying theoretical methods to the ground and low-lying excited electronic states. Employing a large breadth of high-level ab initio computations, including coupled cluster [CCSD(T) and CCSDT(Q)] and multireference configuration interaction [MRCISD+Q] methods, we comprehensively demonstrate that the most recent experimental and theoretical interpretations of the electronic spectrum of H 2 CN are in error. The previous assignments of the two broad features in the spectrum as the origin 0$^{0}_{0}$ (~35 050 cm –1 ) and 4$^{2}_{0}$ (~35 600 cm –1 )$\widetilde{B}$ 2 A 1 ←$\widetilde{X}$ 2 B 2 transitions are both found to be incorrect. The presently reported transition energies suggest that the higher energy band near 35 600 cm –1 is the true origin band. Moreover, from the computed anharmonic vibrational frequencies of the $\widetilde{X}$ 2 B 2 and $\widetilde{B}$ 2 A 1 states, we demonstrate that this ~550 cm –1 band spacing cannot be attributed to a simple vibronic transition, as claimed by the 4$^{2}_{0}$ assignment. Possible alternative explanations for the appearance of the lower intensity band near 35 050 cm –1 are discussed.

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

Anion Photoelectron Spectroscopy and Ab Initio Studies of the UF – Anion

A synergistic anion photoelectron spectroscopic and ab initio computational study of photodetachment of UF – is reported. The measurement determined a vertical detachment energy of 0.63(03) eV, which is consistent with a spinor-based relativistic coupled-cluster CCSD(T) value of 0.61 eV. The complex spectral features due to excited electronic states and vibrational progressions of UF are analyzed and assigned with the help of spin–orbit-coupled multireference perturbation theory and spinor-based relativistic coupled-cluster calculations. UF and UF – are confirmed to be dominated by ionic bonding. Furthermore, the usefulness of the spinor CCSD(T) approach is demonstrated.

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

Modeling Oxidative Dehydrogenation of Propane with Supported Vanadia Catalysts Using Multireference Methods

The oxidative dehydrogenation of propane over supported vanadium oxide catalysts poses significant computational challenges due to complex electronic structure changes along the reaction coordinate, driven primarily by changes in the oxidation states of vanadium. To address these challenges, we systematically test quantum chemical methods, including multireference (MR) approaches, domain-based local pair natural orbital coupled cluster theory (DLPNO-CCSD(T)), and density functional theory (DFT). The initial C–H bond-breaking transition state requires MR treatment due to its multireference character, while subsequent steps permit efficient single-reference calculations. For the rate-limiting C–H activation step mediated by the vanadyl moiety, complete 1 active space second-order perturbation theory (CASPT2) yields an apparent activation barrier (E app 600K ) of 138 kJ/mol, consistent with experimental values (134 ± 4 kJ/mol; Gruene et al. Catal. Today 2010, 157, 137). In contrast, DLPNO-CCSD(T) overestimates this barrier (198 kJ/mol), whereas DFT predictions span 125–150 kJ/mol, depending on the functional. Our multireference investigation of this transition metal oxide-catalyzed process demonstrates that an active space that incorporates the C–H σ and V=O σ/π bonding orbitals, oxygen lone pairs, and their antibonding counterparts adequately captures electronic structure changes along the chemical transformation. Furthermore, these findings provide a general strategy for active space selection in transition metal oxide-catalyzed C/O–H bond activation reactions. The reference dataset from this work, which includes MR calculations with manually selected active spaces for all intermediates and transition states in the propane ODH reaction network, will serve as a benchmark for automating active space selection in similar systems.

Catalysts↗
DOE OSTI2025

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

Clean and Convenient Tessellations for Number Counting Jastrow Factors

We introduce a basis of counting functions that, by cleanly tessellating three-dimensional space, allows real space number counting Jastrow factors to be straightforwardly applied to general molecular situations. Here, by exerting direct control over electron populations in local regions of space and encoding pairwise correlations between these populations, these Jastrow factors allow even very simple reference wave functions to adopt nodal surfaces well suited to many strongly correlated settings. Being trivially compatible with traditional Jastrow factors and diffusion Monte Carlo and having the same cubic per-sample cost scaling as a single determinant trial function, these Jastrow factors thus offer a powerful new route to the simultaneous capture of weak and strong electron correlation effects in a wide variety of molecular and materials settings. In multiple strongly correlated molecular examples, we show that even when paired with the simplest possible single determinant reference, these Jastrow factors allow quantum Monte Carlo to out-perform coupled cluster theory and approach the accuracy of traditional multireference methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

The structures, binding energies and vibrational frequencies of Ca3 and Ca4: An application of the CCSD(T) method

The Ca3 and Ca4 metallic clusters have been investigated using state-of-the-art ab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca3 and Ca4 have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both the e' mode of Ca3 and the t2 mode of Ca4 are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.

Lee, Timothy J.↗
NASA NTRS1991

Triple and quadruple excitation contributions to the binding in Be clusters: Calibration calculations on Be3

The contribution of connected triple and quadruple excitations to the binding in Be3 is investigated by comparing various coupled-cluster (CC) and truncated configuration interaction (CI) treatments with multireference CI (MRCI) and full CI(FCI) calculations. The CC method with single and double excitations (CCSD) produces results that differ substantially from more elaborate treatments, but most extensions to CCSD that account approximately for connected triple excitations perform very well. In constrast, good agreement with FCI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗
DOE OSTI2024

State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space–unitary coupled cluster (LAS–UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. Here, in this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C–C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗