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

Vibrational frequencies for Be3 and Be4

The harmonic frequencies of Be3 and Be4, evaluated by multireference configuration-interaction (MRCI) and singles and doubles coupled-cluster (CCSD) methods, are reported, as are IR intensities evaluated by using the double harmonic approximation. A comparison of the results shows that the CCSD method augmented with a correction for connected triple excitations, CCSD(T), can be a viable alternative to MRCI for determining the ground-state properties of small Be clusters (other than the dimer). Results of normal mode analyses indicate that the anharmonic vibrational frequencies of Be3 and Be4 are reasonably consistent at the CCSD, CCSD(T), or MRCI levels of treatment.

Rendell, Alistair P.↗
NASA NTRS1990

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 contrast, good agreement with CFI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗
DOE OSTI2014

The electronic structure of vanadium monochloride cation (VCl{sup +}): Tackling the complexities of transition metal species

Six electronic states (X {sup 4}Σ{sup −}, A {sup 4}Π, B {sup 4}Δ, {sup 2}Φ, {sup 2}Δ, {sup 2}Σ{sup +}) of the vanadium monochloride cation (VCl{sup +}) are described using large basis set coupled cluster theory. For the two lowest quartet states (X {sup 4}Σ{sup −} and A {sup 4}Π), a focal point analysis (FPA) approach was used that conjoined a correlation-consistent family of basis sets up to aug-cc-pwCV5Z-DK with high-order coupled cluster theory through pentuple (CCSDTQP) excitations. FPA adiabatic excitation energies (T{sub 0}) and spectroscopic constants (r{sub e}, r{sub 0}, B{sub e}, B{sub 0}, D{sup ¯}{sub e}, H{sub e}, ω{sub e}, v{sub 0}, α{sub e}, ω{sub e}x{sub e}) were extrapolated to the valence complete basis set Douglas-Kroll (DK) aug-cc-pV∞Z-DK CCSDT level of theory, and additional treatments accounted for higher-order valence electron correlation, core correlation, and spin-orbit coupling. Due to the delicate interplay between dynamical and static electronic correlation, single reference coupled cluster theory is able to provide the correct ground electronic state (X {sup 4}Σ{sup −}), while multireference configuration interaction theory cannot. Perturbations from the first- and second-order spin orbit coupling of low-lying states with quartet spin multiplicity reveal an immensely complex rotational spectrum relative to the isovalent species VO, VS, and TiCl. Computational data on the doublet manifold suggest that the lowest-lying doublet state ({sup 2}Γ) has a T{sub e} of ∼11 200 cm{sup −1}. Overall, this study shows that laboratory and theoretical rotational spectroscopists must work more closely in tandem to better understand the bonding and structure of molecules containing transition metals.

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

Computational simulation and interpretation of the low-lying excited electronic states and electronic spectrum of thioanisole

Three singlet states, namely a closed-shell ground state and two excited states with 1 ππ* and 1 nσ* character, have been suggested to be responsible for the radiationless decay or photochemical reaction of photoexcited thioanisole. The correct interpretation of the electronic spectrum is critical for understanding the character of these low-lying excited states, but the experimental spectrum is yet to be fully interpreted. In the work reported here, we investigated the nature of those three states and a fourth singlet state of thioanisole using electronic structure calculations by multireference perturbation theory, by completely-renormalized equation-of-motion coupled cluster theory with single and double excitations and noniterative inclusion of connected triples (CR-EOM-CCSD(T)), and by linear-response time-dependent density functional theory (TDDFT). We clarified the assignment of the electronic spectrum by simulating it using a normal-mode sampling approach combined with TDDFT in the Tamm–Dancoff approximation (TDA). The understanding of the electronic states and of the accuracy of the electronic structure methods lays the foundation of our future work of constructing potential energy surfaces.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

An Ab Initio Electronic Structure Investigation of the Ground and Excited States of ScH + , YH + , and LaH +

Multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), coupled-cluster singles, doubles, and perturbative triples [CCSD(T)], and frozen-core full configuration interaction (fcFCI) calculations were carried out using large, correlation-consistent basis sets to investigate the excited states of the Sc atom and the spin–free and spin–orbit coupled potential energy profiles, energetics, spectroscopic constants, and electron populations of low-lying states of MH + (M = Sc, Y, La). The core electron correlation effects, complete basis set effects, and spin–orbit coupling effects were also evaluated. The first four electronic states of all MH + are 1 2 Δ, 1 2 Σ + , 1 2 Π, and 2 2 Σ + with 1σ 2 1δ 1 , 1σ 2 2σ 1 , 1σ 2 1π 1 , and 1σ 2 3σ 1 single-reference electron configurations, respectively. These states of MH + can be represented by the M 2+ H – ionic structure. The ground states of ScH + , YH + , and LaH + are 1 2 Δ 3/2 , 1 2 Σ + 1/2 , and 1 2 Δ 3/2 with 55.45, 60.54, and 62.34 kcal/mol bond energies, respectively. The core electron correlation was found to be vital for gaining accurate predictions on the ground and excited state properties of MH + . The spin–orbit coupling effects are minor for ScH + but become substantial moving to YH + and LaH + . Overall, the results of this work are in good agreement with the limited set of experimental findings of MH + available in the literature and will be of use for future investigations. Furthermore, the theoretical approaches, findings, and trends reported here are expected to aid studies of similar species.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1989

Theoretical investigations of the structures and binding energies of Be(sub n) and Mg(sub n) (n = 3-5) clusters

Researchers determined the equilibrium geometries and binding energies of Be and Mg trimers, tetramers and pentamers using single and double excitation coupled cluster (CCSD) and complete active space self-consistent-field (CASSCF) multireference configuration interaction (MRCI) wave functions in conjunction with extended atomic basis sets. Best estimates of the cluster binding energies are 24, 83 and 110 kcal/mole for Be3, Be4 and Be5; and 9, 31 and 41 kcal/mole for Mg3, Mg4 and Mg5, respectively. A comparison of the MRCI and CCSD results shows that even the best single-reference approach (limited to single and double excitations) is not capable of quantitative accuracy in determining the binding energies of Be and Mg clusters.

Lee, Timothy J.↗
NASA NTRS1990

Theoretical investigations of the structures and binding energies of Be(n) and Mg(n) (n = 3-5) clusters

The equilibrium geometries and binding energies of Be and Mg trimers, tetramers and pentamers have been determined using single and double excitation coupled cluster (CCSD) and complete active space self-consistent-field (CASSCF) multireference configuration interaction (MRCI) wave functions in conjunction with extended atomic basis sets. The best estimates of the cluster binding energies are 24, 83, and 110 kcal/mol for Be3, Be4, and Be5; and 9, 31, and 41 kcal/mol for Mg3, Mg4, and Mg5, respectively. A comparison of the MRCI and CCSD results shows that even the best single-reference approach (limited to single and double excitations) is not capable of quantitative accuracy in determining the binding energies of Be and Mg clusters.

Lee, Timothy J.↗
DOE OSTI2024

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Accurate prediction of diradical chemistry from a single-reference density-matrix method: Model application to the bicyclobutane to gauche-1,3-butadiene isomerization

Multireference correlation in diradical molecules can be captured by a single-reference 2-electron reduced-density-matrix (2-RDM) calculation with only single and double excitations in the 2-RDM parametrization. The 2-RDM parametrization is determined by N-representability conditions that are non-perturbative in their treatment of the electron correlation. Conventional single-reference wave function methods cannot describe the entanglement within diradical molecules without employing triple- and potentially even higher-order excitations of the mean-field determinant. In the isomerization of bicyclobutane to gauche-1,3-butadiene the parametric 2-RDM (p2-RDM) method predicts that the diradical disrotatory transition state is 58.9 kcal/mol above bicyclobutane. This barrier is in agreement with previous multireference calculations as well as recent Monte Carlo and higher-order coupled cluster calculations. The p2-RDM method predicts the Nth natural-orbital occupation number of the transition state to be 0.635, revealing its diradical character. The optimized geometry from the p2-RDM method differs in important details from the complete-active-space self-consistent-field geometry used in many previous studies including the Monte Carlo calculation.

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

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
NASA NTRS1991

Quantum mechanical calculations to chemical accuracy

The accuracy of current molecular-structure calculations is illustrated with examples of quantum mechanical solutions for chemical problems. Two approaches are considered: (1) the coupled-cluster singles and doubles (CCSD) with a perturbational estimate of the contribution of connected triple excitations, or CCDS(T); and (2) the multireference configuration-interaction (MRCI) approach to the correlation problem. The MRCI approach gains greater applicability by means of size-extensive modifications such as the averaged-coupled pair functional approach. The examples of solutions to chemical problems include those for C-H bond energies, the vibrational frequencies of O3, identifying the ground state of Al2 and Si2, and the Lewis-Rayleigh afterglow and the Hermann IR system of N2. Accurate molecular-wave functions can be derived from a combination of basis-set saturation studies and full configuration-interaction calculations.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2016

Ab Initio Reaction Kinetics of CH 3 O$\dot{C}$(=O) and $\dot{C}$H 2 OC(=O)H Radicals

The dissociation and isomerization kinetics of the methyl ester combustion intermediates methoxycarbonyl radical (CH3O$\dot{C}$(=O)) and (formyloxy)methyl radical ($\dot{C}$H2OC(=O)H) are investigated theoretically using high-level ab initio methods and Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Geometries obtained at the hybrid density functional theory (DFT) and coupled cluster singles and doubles with perturbative triples correction (CCSD(T)) levels of theory are found to be similar. We employ high-level ab initio wave function methods to refine the potential energy surface: CCSD(T), multireference singles and doubles configuration interaction (MRSDCI) with the Davidson–Silver (DS) correction, and multireference averaged coupled-pair functional (MRACPF2) theory. MRSDCI+DS and MRACPF2 capture the multiconfigurational character of transition states (TSs) and predict lower barrier heights than CCSD(T). The temperature- and pressure-dependent rate coefficients are computed using RRKM/ME theory in the temperature range 300–2500 K and a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. Dissociation of CH3O$\dot{C}$(=O) to $\dot{C}$H3 and CO2 is predicted to be much faster than dissociating to CH3$\dot{O}$ and CO, consistent with its greater exothermicity. Isomerization between CH3O$\dot{C}$(=O) and $\dot{C}$H2OC(=O)H is predicted to be the slowest among the studied reactions and rarely happens even at high temperature and high pressure, suggesting the decomposition pathways of the two radicals are not strongly coupled. The predicted rate coefficients and branching fractions at finite pressures differ significantly from the corresponding high-pressure-limit results, especially at relatively high temperatures. Finally, because it is one of the most important CH3$\dot{O}$ removal mechanisms under atmospheric conditions, the reaction kinetics of CH3$\dot{O}$ + CO was also studied along the PES of CH3O$\dot{C}$(=O); the resulting kinetics predictions are in remarkable agreement with experiments.

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

Interaction of Th with H 0/–/+ : Combined Experimental and Theoretical Thermodynamic Properties

In this work, high-level electronic structure calculations of the lowlying energy electronic states for ThH, ThH – , and ThH + are reported and compared to experimental measurements. The inclusion of spin–orbit coupling is critical to predict the ground-state ordering as inclusion of spin–orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH + . At the multireference spin–orbit SO-CASPT2 level, the ground states of ThH, ThH – , and ThH + are predicted to be the 2 Δ 3/2 , 3 Φ 2 , and 3 Δ 1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH – photoelectron spectrum has many transitions, which approximately correlate with excitations of Th + and/or Th. The adiabatic ionization energy of ThH including spin–orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th + H – ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH – , and ThH + using the Feller–Peterson–Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

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

A diagnostic for determining the quality of single-reference electron correlation methods

It was recently proposed that the Euclidian norm of the t sub 1 vector of the coupled cluster wave function (normalized by the number of electrons included in the correlation procedure) could be used to determine whether a single-reference-based electron correlation procedure is appropriate. This diagnostic, T sub 1, is defined for use with self consistent field molecular orbitals and is invariant to the same orbital rotations as the coupled cluster energy. T sub 1 is investigated for several different chemical systems which exhibit a range of multireference behavior, and is shown to be an excellent measure of the importance of nondynamical electron correlation and is far superior to C sub 0 from a singles and doubles configuration interaction wave function. It is further suggested that when the aim is to recover a large fraction of the dynamical electron correlation energy, a large T sub 1 (i.e., greater than 0.02) probably indicates the need for a multireference electron correlation procedure.

Lee, Timothy J.↗
NASA NTRS1989

A diagnostic for determining the quality of single-reference electron correlation methods

It was recently proposed that the Euclidian norm of the t(sub 1) vector of the coupled cluster wave function (normalized by the number of electrons included in the correlation procedure) could be used to determine whether a single-reference-based electron correlation procedure is appopriate. This diagnostic, T(sub 1) is defined for use with self-consistent-field molecular orbitals and is invariant to the same orbital rotations as the coupled cluster energy. T(sub 1) is investigated for several different chemical systems which exhibit a range of multireference behavior, and is shown to be an excellent measure of the importance of non-dynamical electron correlation and is far superior to C(sub 0) from a singles and doubles configuration interaction wave function. It is further suggested that when the aim is to recover a large fraction of the dynamical electron correlation energy, a large T(sub 1) (i.e., greater than 0.02) probably indicates the need for a multireference electron correlation procedure.

Lee, Timothy J.↗
DOE OSTI2024

The Localized Active Space Method with Unitary Selective Coupled Cluster

Here, we introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H 2 ) 2 , (H 2 ) 4 , and trans-butadiene, and the magnetic coupling constant for a bimetallic compound [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ . For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.

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

Size-extensivity-corrected multireference configuration interaction schemes to accurately predict bond dissociation energies of oxygenated hydrocarbons

Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this scheme's use on larger molecules and for more complex regions of combustion PESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗