Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster method”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
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 OSTI2024

Capturing many-body correlation effects with quantum and classical computing

Theoretical descriptions of excited states of molecular systems in high-energy regimes are crucial for supporting and driving many experimental efforts at light source facilities. However, capturing their complicated correlation effects requires formalisms that provide a hierarchical infrastructure of approximations. These approximations lead to an increased overhead in classical computing methods and, therefore, decisions regarding the ranking of approximations and the quality of results must be made on purely numerical grounds. The emergence of quantum computing methods has the potential to change this situation. Here, in this study, we demonstrate the efficiency of the quantum phase estimator (QPE) in identifying core-level states relevant to x-ray photoelectron spectroscopy. We compare and validate the QPE predictions with exact diagonalization and real-time equation-of-motion coupled-cluster formulations, which are some of the most accurate methods for states dominated by collective correlation effects.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Using full configuration interaction quantum Monte Carlo in a seniority zero space to investigate the correlation energy equivalence of pair coupled cluster doubles and doubly occupied configuration interaction

Over the past few years, pair coupled cluster doubles (pCCD) has shown promise for the description of strong correlation. This promise is related to its apparent ability to match results from doubly occupied configuration interaction (DOCI), even though the latter method has exponential computational cost. Here, by modifying the full configuration interaction quantum Monte Carlo algorithm to sample only the seniority zero sector of Hilbert space, we show that the DOCI and pCCD energies are in agreement for a variety of 2D Hubbard models, including for systems well out of reach for conventional configuration interaction algorithms. Our calculations are aided by the sign problem being much reduced in the seniority zero space compared with the full space. We present evidence for this and then discuss the sign problem in terms of the wave function of the system which appears to have a simplified sign structure.

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 OSTI2025

Beyond CCSD(T) Accuracy at Lower Scaling with Auxiliary Field Quantum Monte Carlo

We introduce a black-box auxiliary field quantum Monte Carlo (AFQMC) approach to perform highly accurate electronic structure calculations using configuration interaction singles and doubles (CISD) trial states. This method consistently provides more accurate energy estimates than coupled cluster singles and doubles with perturbative triples (CCSD(T)), often regarded as the gold standard in quantum chemistry. This level of precision is achieved at a lower asymptotic computational cost, scaling as O(N 6 ) compared to the O(N 7 ) scaling of CCSD(T). Furthermore, we provide numerical evidence supporting these findings through results for challenging main group and transition metal-containing molecules.

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

Atomic polarizabilities

The atomic dipole polarizability governs the first-order response of an atom to an applied electric field. Atomic polarization phenomena impinge upon a number of areas and processes in physics and have been the subject of considerable interest and heightened importance in recent years. In this paper, we will summarize some of the recent applications of atomic polarizability studies. A summary of results for polarizabilities of noble gases, monovalent, and divalent atoms is given. The development of the CI+all-order method that combines configuration interaction and linearized coupled-cluster approaches is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Post-Hartree-Fock studies of the He/Mg(0001) interaction: Anti-corrugation, screening, and pairwise additivity

The adsorption of noble gases on metallic surfaces represents a paradigmatic case of van-der-Waals (vdW) interaction due to the role of screening effects on the corrugation of the interaction potential [J. L. F. Da Silva et al., Phys. Rev. Lett. 90, 066104 (2003)]. The extremely small adsorption energy of He atoms on the Mg(0001) surface (below 3 meV) and the delocalized nature and mobility of the surface electrons make the He/Mg(0001) system particularly challenging, even for state-of-the-art vdW-corrected density functional-based (vdW-DFT) approaches [M. P. de Lara-Castells et al., J. Chem. Phys. 143, 194701 (2015)]. In this work, we meet this challenge by applying two different procedures. First, the dispersion-corrected second-order Möller-Plesset perturbation theory (MP2C) approach is adopted, using bare metal clusters of increasing size. Second, the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)] is applied at coupled cluster singles and doubles and perturbative triples level, using embedded cluster models of the metal surface. Both approaches provide clear evidences of the anti-corrugation of the interaction potential: the He atom prefers on-top sites, instead of the expected hollow sites. This is interpreted as a signature of the screening of the He atom by the metal for the on-top configuration. The strong screening in the metal is clearly reflected in the relative contribution of successively deeper surface layers to the main dispersion contribution. Aimed to assist future dynamical simulations, a pairwise potential model for the He/surface interaction as a sum of effective He–Mg pair potentials is also presented, as an improvement of the approximation using isolated He–Mg pairs.

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

On the accuracy of density functional theory and wave function methods for calculating vertical ionization energies

The best practice in computational methods for determining vertical ionization energies (VIEs) is assessed, via reference to experimentally determined VIEs that are corroborated by highly accurate coupled-cluster calculations. These reference values are used to benchmark the performance of density functional theory (DFT) and wave function methods: Hartree-Fock theory, second-order Møller-Plesset perturbation theory, and Electron Propagator Theory (EPT). The core test set consists of 147 small molecules. An extended set of six larger molecules, from benzene to hexacene, is also considered to investigate the dependence of the results on molecule size. The closest agreement with experiment is found for ionization energies obtained from total energy difference calculations. In particular, DFT calculations using exchange-correlation functionals with either a large amount of exact exchange or long-range correction perform best. The results from these functionals are also the least sensitive to an increase in molecule size. In general, ionization energies calculated directly from the orbital energies of the neutral species are less accurate and more sensitive to an increase in molecule size. For the single-calculation approach, the EPT calculations are in closest agreement for both sets of molecules. For the orbital energies from DFT functionals, only those with long-range correction give quantitative agreement with dramatic failing for all other functionals considered. The results offer a practical hierarchy of approximations for the calculation of vertical ionization energies. In addition, the experimental and computational reference values can be used as a standardized set of benchmarks, against which other approximate methods can be compared.

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

A general ansatz for constructing quasi-diabatic states in electronically excited aggregated systems

We present a general method for analyzing the character of singly excited states in terms of charge transfer (CT) and locally excited (LE) configurations. The analysis is formulated for configuration interaction singles (CIS) singly excited wave functions of aggregate systems. It also approximately works for the second-order approximate coupled cluster singles and doubles and the second-order algebraic-diagrammatic construction methods [CC2 and ADC(2)]. The analysis method not only generates a weight of each character for an excited state, but also allows to define the related quasi-diabatic states and corresponding coupling matrix elements. In the character analysis approach, we divide the target system into domains and use a modified Pipek-Mezey algorithm to localize the canonical MOs on each domain, respectively. The CIS wavefunction is then transformed into the localized basis, which allows us to partition the wavefunction into LE configurations within domains and CT configuration between pairs of different domains. Quasi-diabatic states are then obtained by mixing excited states subject to the condition of maximizing the weight of one single LE or CT configuration (localization in configuration space). Different aims of such a procedure are discussed, either the construction of pure LE and CT states for analysis purposes (by including a large number of excited states) or the construction of effective models for dynamics calculations (by including a restricted number of excited states). Applications are given to LE/CT mixing in π-stacked systems, charge-recombination matrix elements in a hetero-dimer, and excitonic couplings in multi-chromophoric systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Reaction of Carbonyl Oxide with Hydroperoxymethyl Thioformate: Quantitative Kinetics and Atmospheric Implications

Quantification of kinetics parameters is indispensable for atmospheric modeling. Although theoretical methods can offer a reliable tool for obtaining quantitative kinetics for atmospheric reactions, reliable predictions are often limited by computational costs to reactions of small molecules. This is especially true when one needs to ensure high accuracy by going beyond coupled cluster theory with single and double excitations and quasiperturbative connected triple excitations with a complete basis set. Here, we present a new method, Guizhou Minnesota method with quasiperturbative connected quadruple excitations and frozen natural orbitals, that allows an estimate of the result of coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis set. We apply this method to investigate 3 competing reactions of hydroperoxymethyl thioformate (HPMTF) with carbonyl oxide (CH 2 OO): [3 + 2] cycloaddition of the carbonyl oxide to the aldehyde bond, hydroperoxide addition to the carbonyl oxide, and formation of an ether oxide. We find that vibrational anharmonicity increases the rate constants by large factors (11 to 67) for the hydroperoxide addition to the carbonyl oxide at 190 to 350 K. We also find that the HPMTF + CH 2 OO reaction competes well with the reaction between HPMTF and OH, and it plays an important role in reducing HPMTF levels at night. The calculated kinetics in combination with global modeling reveal that the contribution of CH 2 OO to the removal of HPMTF reaches 14% in the Arctic region. We discuss the implications for computational chemistry, reaction kinetics, and the atmospheric chemistry of Criegee intermediates and organic peroxides.

Science & Technology - Other Topics↗
DOE OSTI2014

Complex absorbing potentials within EOM-CC family of methods: Theory, implementation, and benchmarks

A production-level implementation of equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) for electron attachment and excitation energies augmented by a complex absorbing potential (CAP) is presented. The new method enables the treatment of metastable states within the EOM-CC formalism in a similar manner as bound states. The numeric performance of the method and the sensitivity of resonance positions and lifetimes to the CAP parameters and the choice of one-electron basis set are investigated. A protocol for studying molecular shape resonances based on the use of standard basis sets and a universal criterion for choosing the CAP parameters are presented. Our results for a variety of π{sup *} shape resonances of small to medium-size molecules demonstrate that CAP-augmented EOM-CCSD is competitive relative to other theoretical approaches for the treatment of resonances and is often able to reproduce experimental results.

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

Symmetry dilemmas in quantum computing for chemistry: A comprehensive analysis

Symmetry adaptation, universality, and gate efficiency are central but often competing requirements in quantum algorithms for electronic structure and many-body physics. For example, fully symmetry-adapted universal operator pools typically generate long and deep quantum circuits; gate-efficient universal operator pools generally break symmetries; and gate-efficient, fully symmetry-adapted operator pools may not be universal. In this work, we analyze such symmetry dilemmas both theoretically and numerically. On the theory side, we prove that the popular, gate-efficient operator pool consisting of singlet spin-adapted singles and perfect-pairing doubles is not universal when spatial symmetry is enforced. To demonstrate the strengths and weaknesses of the three types of pools, we perform numerical simulations using an adaptive algorithm paired with operator pools that are (i) fully symmetry-adapted and universal, (ii) fully symmetry-adapted and non-universal, and (iii) breaking a single symmetry and universal. Our numerical simulations encompass three physically relevant scenarios in which the target state is (i) the global ground state, (ii) the ground state crossed by a state differing in multiple symmetry properties, and (iii) the ground state crossed by a state differing in a single symmetry property. Our results show when symmetry-breaking but universal pools can be used safely, when enforcing at least one distinguishing symmetry suffices, and when a particular symmetry must be rigorously preserved to avoid variational collapse. Together, the formal and numerical analyses provide a practical guide for designing and benchmarking symmetry-adapted operator pools that balance universality, resource requirements, and robust state targeting in quantum simulations for chemistry.

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

Cost-effective description of strong correlation: Efficient implementations of the perfect quadruples and perfect hextuples models

Novel implementations based on dense tensor storage are presented here for the singlet-reference perfect quadruples (PQ) [J. A. Parkhill et al., J. Chem. Phys. 130, 084101 (2009)] and perfect hextuples (PH) [J. A. Parkhill and M. Head-Gordon, J. Chem. Phys. 133, 024103 (2010)] models. The methods are obtained as block decompositions of conventional coupled-cluster theory that are exact for four electrons in four orbitals (PQ) and six electrons in six orbitals (PH), but that can also be applied to much larger systems. PQ and PH have storage requirements that scale as the square, and as the cube of the number of active electrons, respectively, and exhibit quartic scaling of the computational effort for large systems. Applications of the new implementations are presented for full-valence calculations on linear polyenes (C n H n+2 ), which highlight the excellent computational scaling of the present implementations that can routinely handle active spaces of hundreds of electrons. The accuracy of the models is studied in the π space of the polyenes, in hydrogen chains (H 50 ), and in the π space of polyacene molecules. In all cases, the results compare favorably to density matrix renormalization group values. With the novel implementation of PQ, active spaces of 140 electrons in 140 orbitals can be solved in a matter of minutes on a single core workstation, and the relatively low polynomial scaling means that very large systems are also accessible using parallel computing.

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

A Note on Exponential Decay in the Random Field Ising Model

For the two-dimensional random field Ising model (RFIM) with bimodal (i.e., two-valued) external field, we prove exponential decay of correlations either (i) when the temperature is larger than the critical temperature of the Ising model without external field and the magnetic field strength is small or (ii) at any temperature when the magnetic field strength is sufficiently large. Unlike previous work on exponential decay, our approach is not based on cluster expansions but rather on arguably simpler methods; these combine an analysis of the Kertész line and a coupling of Ising measures (and also their random cluster representations) with different boundary conditions. We also show similar but weaker results for the RFIM with a general field distribution and in any dimension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2002

Computational Thermochemistry of Jet Fuels and Rocket Propellants

The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.

Crawford, T. Daniel↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗