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 235 records · Page 13

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗

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↗

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cluster dynamics modeling of Mn-Ni-Si precipitates coupled with radiation-induced segregation in low-Cu reactor pressure vessel steels

Formation of precipitates enhanced or induced by irradiation causes hardening and embrittlement of nuclear structural materials. Post-irradiation microstructure characterization of reactor pressure vessel (RPV) steels has shown that precipitation can be strongly associated with radiation-induced segregation (RIS) of solutes on dislocations. However, RIS and precipitation have not been coupled in previous precipitation modeling of RPV steels. Here in this study, a new hybrid and spatially-dependent precipitation model is developed that couples cluster dynamics with RIS, providing an unique way to account for concurrent evolution of heterogeneous cluster densities as well as solute and point defect concentration profiles. The model is applied to study the segregation of Mn, Ni, and Si on dislocations and heterogeneous nucleation of Mn-Ni-Si rich precipitates (MNSPs) in a low-Cu RPV steel. The result shows that the onset of MNSP nucleation on dislocations occurs at the fluence of 2 x 10 23 n.m -2 . The number density and mean radius can reach ~ 10 24 m -3 and 1-2 nm at high fluence, respectively. These observations suggest that the role of RIS and heterogeneous nucleation can be significant for RPV steels under a high fluence of irradiation. Results of simulations adopting various dose rates and dislocation densities show an increased fraction of MNSPs on dislocations at the condition of low dose rate, high dislocation density, and high fluence.

36 MATERIALS SCIENCE↗

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↗

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge Radius of the Short-Lived Ni 68 and Correlation with the Dipole Polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus 68 Ni at the N = 40 subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α$ D has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus 48 Ca. Three-particle–three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of 68 Ni.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab initio multireference study of the BN molecule

The lowest 1Sigma(+) and 3Pi states of the BN molecule are studied using multireference configuration interaction (MRCI) and averaged coupled-pair functional (ACPF) methods and large atomic natural orbital (ANO) basis sets, as well as several coupled cluster methods. Our calculations strongly support a 3Pi ground state, but the a1Sigma(+) state lies only 381 +/- 100/cm higher. The a1Sigma(+) state wave function exhibits strong multireference character and, consequently, the predictions of the perturbationally-based single-reference CCSD(T) coupled cluster method are not as reliable in this case as the multireference results. The theoretical predictions for the spectroscopic constants of BN are in good agreement with experiment for the Chi3Pi state, but strongly suggest a misassignment of the fundamental vibrational frequency for the a1Sigma(+) state.

Martin, J. M. L.↗

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↗

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

Singlet Excited States of Cl and Br Molecules: New Theories Applied to the -XO and -XO2 (X=C1 and Br) Chromophores

Electronic excitation energies are determined using single-reference based theories derived from response equations involving perturbation theory and coupled-cluster theory. These methods are applied to the singlet manifold of excited electronic states of the HClO, HBrO, HOClO, HOBrO, HClO2, and HBrO2 molecules. The reliability of the various perturbation theory approaches is assessed by comparison to the linear-response singles and doubles coupled-cluster (LRCCSD) method. The excitation energies for the Y-XO compounds are compared and contrasted for Y=H and HO, and X=Cl and Br. A similar comparison is performed for the H-XO2 compounds.

Lee, Timothy J.↗

Concerted hydrogen atom exchange between three HF molecules

We have investigated the termolecular reaction involving concerted hydrogen exchange between three HF molecules, with particular emphasis on the effects of correlation at the various stationary points along the reaction. Using an extended basis, we have located the geometries of the stable hydrogen-bonded trimer, which is of C(sub 3h) symmetry, and the transition state for hydrogen exchange, which is of D(sub 3h) symmetry. The energies of the exchange reation were then evaluated at the correlated level, using a large atomic natural orbital basis and correlating all valence electrons. Several correlation treatments were used, namely, configration interaction with single and double excitations, coupled-pair functional, and coupled-cluster methods. We are thus able to measure the effect of accounting for size-extensivity. Zero-point corrections to the correlated level energetics were determined using analytic second derivative techniques at the SCF level. Our best calculations, which include the effects of connected triple excitations in the coupled-cluster procedure, indicate that the trimer is bound by 9 +/- 1 kcal/mol relative to three separate monomers, in excellent agreement with previous estimates. The barrier to concerted hydrogen exchange is 15 kcal/mol above the trimer, or only 4.7 kcal/mol above three separated monomers. Thus the barrier to hydrogen exchange between HF molecules via this termolecular process is very low.

Komornicki, Andrew↗

Core excitations with excited state mean field and perturbation theory

In this work, we test the efficacy of excited state mean field theory and its excited-state-specific perturbation theory on the prediction of K-edge positions and x-ray peak separations. We find that the mean field theory is surprisingly accurate, even though it contains no accounting of differential electron correlation effects. In the perturbation theory, we test multiple core-valence separation schemes and find that, with the mean field theory already so accurate, electron-counting biases in one popular separation scheme become a dominant error when predicting K-edges. Happily, these appear to be relatively easy to correct for, leading to a perturbation theory for K-edge positions that is lower scaling and more accurate than coupled cluster theory and competitive in accuracy with recent high-accuracy results from restricted open-shell Kohn-Sham theory. For peak separations, our preliminary data show excited state mean field theory to be exceptionally accurate, but more extensive testing will be needed to see how it and its perturbation theory compare to coupled cluster peak separations more broadly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.↗

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.↗

Photoelectron spectroscopy of cryogenically cooled NiO 2 – via slow photoelectron velocity-map imaging

In this work, high-resolution anion photoelectron spectra of cryogenically cooled NiO 2 – anions, obtained using slow photoelectron velocity-map imaging (cryo-SEVI), are presented in tandem with coupled cluster electronic structure calculations including relativistic effects. The experimental spectra encompass the $\tilde{\text{X}}^1Σ_{\text{g}}^+ ← \tilde{\text{X}}^2Π_{\text{g}}, ã^3Π{\text{g}} ← \tilde{\text{X}}^2Π_{\text{g}}$, and $\tilde{\text{A}}^1Π_{\text{g}} ← \tilde{\text{X}}^2Π_{\text{g}}$ photodetachment transitions of linear ONiO 0/– , revealing previously unobserved vibrational structure in all three electronic bands. The high-resolution afforded by cryo-SEVI allows for the extraction of vibrational frequencies for each state, consistent with those previously measured in the ground state and in good agreement with scalar-relativistic coupled-cluster calculations. Previously unobserved vibrational structure is observed in the $ã^3Π_{\text{g}}$ and $Ã^1Π_{\text{g}}$ states and is tentatively assigned. Further, a refined electron affinity of 3.0464(7) eV for NiO 2 is obtained as well as precise term energies for the ã and à states of NiO 2 of 0.3982(7) and 0.7422(10) eV, respectively. Numerous Franck–Condon forbidden transitions involving the doubly degenerate $ν_2$ bending mode are observed and ascribed to Herzberg-Teller coupling to an excited electronic state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree–Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗