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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

In search of molecular ions for optical cycling: a difficult road

Optical cycling, a continuous photon scattering off atoms or molecules, plays a central role in the quantum information science. While optical cycling has been experimentally achieved for many neutral species, few molecular ions have been investigated. We present a systematic theoretical search for diatomic molecular ions suitable for optical cycling using equation-of-motion coupled-cluster methods. Inspired by the electronic structure patterns of laser-cooled neutral molecules, we establish the design principles for molecular ions and explore various possible cationic molecular frameworks. The results show that finding a perfect molecular ion for optical cycling is challenging, yet possible. Among various possible diatomic molecules we suggest several candidates, which require further attention from both theory and experiment: YF + , SiO + , PN + , SiBr + , and BO + .

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

Diffusion Monte Carlo evaluation of disiloxane linearisation barrier

The disiloxane molecule is a prime example of silicate compounds containing the Si-O-Si bridge. The molecule is of significant interest within the field of quantum chemistry, owing to the difficulty in theoretically predicting its properties. Herein, the linearisation barrier of disiloxane is investigated using a fixed-node diffusion Monte Carlo (FNDMC) approach, which is one of the most reliable ab initio methods in accounting for the electronic correlation. Calculations utilizing the density functional theory (DFT) and the coupled cluster method with single and double substitutions, including noniterative triples (CCSD(T)) are carried out alongside FNDMC for comparison. It is concluded that FNDMC successfully predicts the disiloxane linearisation barrier and does not depend on the completeness of the basis-set as much as DFT or CCSD(T), thus establishing its suitability.

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

The CH(X 2 Π) + H 2 O reaction: two transition state kinetics

he reaction of ground state methylidyne (CH) with water vapor (H 2 O) is theoretically re-investigated using high-level coupled cluster computations in combination with semi-classical transition state theory (SCTST) and two-dimensional master equation simulations. Insertion of CH into a H–O bond of H 2 O over a submerged barrier via a well-skipping mechanism yielding solely H and CH 2 O is characterized. The reaction kinetics is effectively determined by the formation of a pre-reaction van der Waals complex (PRC, HC—OH 2 ) and its subsequent isomerization to activated CH 2 OH in competition with PRC re-dissociation. The tunneling effects are found to be minor, while variational effects in the PRC → CH 2 OH step are negligible. The calculated rate coefficient k ( T ) is nearly pressure-independent, but strongly depends on temperature with pronounced down-up behavior: a high value of 2 x 10 -10 cm 3 s -1 at 50 K, followed by a fairly steep decrease down to 8 × 10 -12 cm 3 s -1 at 900 K, but increasing again to 5 × 10 -11 cm 3 s -1 at 3500 K. Over the T-range of this work, k ( T ) can be expressed as: k ( T , P = 0) = 2.31 × 10 -11 ( T /300 K) -1.615 exp(-38.45/ T ) cm 3 s -1 for T = 50-400 K k ( T , P = 0) = 1.15 × 10 -12 ( T /300 K) 0.8637 exp(892.6/ T ) cm 3 s -1 for T = 400-1000 K k ( T , P = 0) = 4.57 × 10 -15 ( T /300 K) 3.375 exp(3477.4/ T ) cm 3 s -1 for T = 1000–3500 K.

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

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

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

Mechanism, thermochemistry, and kinetics of the reversible reactions: C 2 H 3 + H 2 ⇌ C 2 H 4 + H ⇌ C 2 H 5

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations, was used in a study of the title reactions, which play an important role in the combustion of hydrocarbons. In the set of radical/radical reactions leading to soot formation in flames, the addition of H-atoms to alkenes is likely a common reaction, triggering the isomerization of complex hydrocarbons to aromatics. The heats of formation of C 2 H 3 , C 2 H 4 , and C 2 H 5 are established to be 301.26 ± 0.30 at 0 K (297.22 ± 0.30 at 298 K), 60.89 ± 0.11 (52.38 ± 0.11), and 131.38 ± 0.22 (120.63 ± 0.22) kJ mol -1 , respectively. The calculated rate constants from first principles agree well with experiments where they are available. Under conditions typical of high temperature combustion – where experimental work is very challenging with a consequent dearth of accurate data –here we provide high-level theoretical results for kinetic modeling.

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

High-accuracy first-principles-based rate coefficients for the reaction of OH and CH 3 OOH

Here, the ˙OH-initiated oxidation of methyl hydroperoxide, which plays an important role in the atmospheric chemistry of methane, was theoretically characterized using high-accuracy composite amHEAT-345(Q) coupled-cluster calculations followed by a two-dimensional E,J resolved master equation analysis. The reaction is found to proceed through two distinct hydrogen-bonded pre-reactive complexes leading to two product channels, in accord with the experimental observations: (i) ˙OH + CH 3 OOH → CH 3 OO˙ + H 2 O with a yield of 0.8 ± 0.1, and (ii) ˙OH + CH 3 OOH → HCHO + ˙OH + H 2 O with a yield of 0.2 ± 0.1. The calculated reaction enthalpies are within 0.2 kcal mol -1 of the benchmark ATcT values. Overall thermal rate coefficients obtained from first principles are found to be in the low-pressure limit at atmospheric pressure; the total rate coefficient can be expressed over the T = 200–450 K range as k(T) = 5.0 × 10 -12 × T -0.152 × exp(287/T) cm 3 s -1 , strongly supporting the experimental results of Vaghjiani and Ravishankara, with which this expression agrees within ca. 15%. The current results show that (i) is the principal reaction channel and support the view that, due to its inherently fast transformations, CH 3 OOH is an important redistribution species for HO x ˙ radicals in the Earth's atmosphere.

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

Theoretical study of the CO 2 –O 2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO 2 –O 2 ) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l –3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO 2 –O 2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ~ –243 cm –1 . Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ~ –218 cm –1 . The transition structure connecting the two minima is found at V ~ –211 cm –1 . We also performed comparisons with some CO 2 –X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO 2 and O 2 . Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO 2 –O 2 mixtures.

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

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S 1 and S 2 . Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S 1 –S 6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S 0 –S 6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. In conclusion, we compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Rishi, Varun↗
DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

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

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗
DOE OSTI2024

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
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 OSTI2026

Gaunt and Breit two-electron contributions to mean-field transformations and fine structure splitting

Materials utilized by novel energy systems are often studied using weakly correlated mean-field theories. However, if these systems incorporate heavy elements, relativistic effects must be included. Therefore, a Kramers unrestricted coupled cluster with singles and doubles excitation formalism within a molecular mean-field exact two-component framework (X2C mmf ) using a four-component Dirac–Hartree–Fock (DHF) reference state is presented. The exact X2C mmf transformed normal-order Hamiltonian incorporates all one-electron and two-electron (2e) contributions from the Coulomb, Gaunt, and Breit operators and is used with the equation of motion method to calculate the excitation energies of the alkali group of elements. Using this framework, the effects of 2e Gaunt and Breit integrals are studied. Results demonstrate growing contributions from these integrals to the generated X2C mmf mean-fields and electronic fine structure calculations with increasing atomic number. Overall, this paper outlines the method, its effect within the X2C mmf approach, and lays the foundation for future theoretical development of relativistic calculations within this framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The unusual symmetry of hexafluoro-o-xylene—A microwave spectroscopy and computational study

The rotational constants and quartic centrifugal distortion coefficients of hexafluoro-o-xylene were precisely derived from the 8 GHz to 18 GHz gas phase microwave spectrum. In addition, the rotational constants of all singly substituted 13 C isotopologues were determined. Instead of the intuitively expected symmetry of C 2v , as in o-xylene, calculations with a variety of methods (B3LYP, CAM-B3LYP, ωB97XD, MP2, and coupled-cluster singles, doubles, and perturbative triples) predict a C 2 symmetry structure in which the two CF 3 groups rotate in opposite directions by about 16°. While the experimental results in this study are not capable of proving this unusual symmetry, they can support future microwave, gas phase electron diffraction or nuclear magnetic resonance studies. From the presented data, a preliminary r 0 structure was determined, reproducing the experimental rotational constants with deviations of no more than 15 kHz. Analysis of the interactions between the two CF 3 groups using an effective fragment potential approach identifies two major contributions to their interaction, due to exchange-repulsion and electrostatic repulsion, with electrostatic repulsion responsible for the barrier at the C 2v geometry.

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

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

Spectroscopic and theoretical studies of UN and UN +

The low-energy electronic states of UN and UN+ have been examined using high-level electronic structure calculations and two-color photoionization techniques. The experimental measurements provided an accurate ionization energy for UN (IE = 50 802 +/- 5 cm -1 ). Spectra for UN+ yielded ro-vibrational constants and established that the ground state has the electronic angular momentum projection omega = 4. Ab initio calculations were carried out using the spin-orbit state interacting approach with the complete active space second-order perturbation theory method. A series of correlation consistent basis sets were used in conjunction with small-core relativistic pseudopotentials on U to extrapolate to the complete basis set limits. The results for UN correctly obtained an omega = 3.5 ground state and demonstrated a high density of configurationally related excited states with closely similar ro-vibrational constants. Similar results were obtained for UN + , with reduced complexity owing to the smaller number of outer-shell electrons. The calculated IE for UN was in excellent agreement with the measured value. Improved values for the dissociation energies of UN and UN + , as well as their heats of formation, were obtained using the Feller-Peterson-Dixon composite thermochemistry method, including corrections up through coupled cluster singles, doubles, triples and quadruples. An analysis of the ab initio results from the perspective of the ligand field theory shows that the patterns of electronic states for both UN and UN + can be understood in terms of the underlying energy level structure of the atomic metal ion.

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

Three-dimensional potential energy surfaces of ArNO ($\tilde{X}$ 2 Π)

Until now, the potential energy surfaces (PESs) of the ArNO complex found in the literature were two-dimensional, with the NO interatomic distance being fixed. In this work, we present the first accurate three-dimensional ground state X 2 Π PESs (both A ' and A ") of ArNO computed at the CCSD(T)/CBS level of theory. The equilibrium geometries and the well depths ( D e ) are compared to several other electronic structure methods. We found that using the multireference method, MRCI-F12 makes the surfaces much shallower (by 25%) and the depth of the surfaces does not agree with experimental data. The explicitly correlated coupled-cluster method underestimates the well depth as well. Analytic representations for both A ' and A " surfaces were fit to 4380 ab initio points to within 2.71 cm -1 .A three-dimensional Numerov propagator method in Delves coordinates is used to compute the bound state spectrum up to J tot = 6.5. The recommended dissociation energies are D 0 = 97.2 cm -1 for the adiabatic ground state and D e = 133.7 (128.1) cm -1 for A ' ( A ")

36 MATERIALS SCIENCE↗