Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis sets”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

NWChem: Recent and Ongoing Developments

In this paper we summarize developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0). Specifically, we focus on functionalities, along with input blocks, that are currently accessible in the current stable release (NWChem 7.2) and master branches, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem GitHub repository, and containerization of NWChem executable images. In conclusion, some of the ongoing developments that will be available in the near future are also discussed.

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

Robust Tensor Hypercontraction of the Particle–Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle–particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O(N 6 ) to O(N 5 ) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. Here, we find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5x.

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

RT-EOM-CCSD Calculations of Inner and Outer Valence Ionization Energies and Spectral Functions

Photoelectron spectroscopy (PES) is a standard experimental method for material characterization, but its interpretation can be hampered by its reliance on standard materials. To facilitate the study of unknown systems, theoretical methods are desirable. Here we present a real-time equation-of-motion coupled cluster (RT-EOM-CC) approach for valence PES, extending our core-level development. Here, we demonstrate that RT-EOM-CC yields ionization energies and spectral functions in good agreement with experiment and CI-based methods, even for some more correlated cases.

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

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

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

Capturing Correlation Effects in Positron Binding to Atoms and Molecules

A major challenge in contemporary electronic structure theory involves the development of methods to describe in a balanced manner the contribution of correlation effects to energy differences. This challenge can be even greater for multicomponent systems containing more than one type of quantum particle. In the present work, we describe a flexible code for carrying out self-consistent field and configuration interaction (CI) calculations on multicomponent systems and use it to generate trial wave functions for use in diffusion Monte Carlo (DMC) calculations of the positron affinity of Be, Be 2 , Be 4 , Mg, CS 2 , and benzene. The resulting positron affinities (PAs) are in good agreement with the best values from the literature.

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 OSTI2025

Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.

Basis sets↗
DOE OSTI2026

One-Step Relativistic Driven Similarity Renormalization Group Multireference Perturbation Theory

We present an efficient implementation of a one-step relativistic second-order multireference perturbation theory based on the multireference driven similarity renormalization group (MR-DSRG) using the exact two-component (X2C) Hamiltonian, which we denote X2C-DSRG-MRPT2. We show that the X2C-DSRG-MRPT2 method can accurately capture spin–orbit coupling (SOC) effects in the electronic structure of strongly correlated systems containing elements across the periodic table. We further demonstrate that the X2C-DSRG-MRPT2 method, through its variational treatment of SOC effects, can yield spin–orbit splittings with mean absolute percentage errors consistently below 7% with respect to experimental values for systems containing up to sixth row elements. With its modest computational scaling (fourth power in system size for the perturbative step) and high accuracy, X2C-DSRG-MRPT2 provides a promising avenue for the routine treatment of relativistic effects in strongly correlated molecular systems.

Hamiltonians↗
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↗
DOE OSTI2022

Molecular Properties of Thorium Hydrides: Electron Affinities and Thermochemistry

In this work, high-level electronic structure calculations of the ground and low-lying energy electronic states for ThH x and ThH x – for x = 2–5 are reported and compared to available anion photoelectron detachment experiments. The adiabatic electron affinities (EAs) are predicted to be 0.82, 0.88, 0.51, and 2.36 eV for x = 2 to 5, respectively, at the Feller– Peterson–Dixon (FPD) level. The vertical detachment energies (VDEs) are predicted to be 0.84, 0.88, 0.81, and 4.38 eV for x = 2–5, respectively. The corresponding experimental VDEs are 0.871 eV for x = 2, 0.88 eV for x = 3, and 4.09 eV for x = 5. As for ThH, there is a significant spin–orbit (SO) correction for the EA of ThH 2 , and this correction decreases substantially for x > 2. The observed ThH 2 – photoelectron spectrum has many transitions as predicted at the CASPT2-SO level. The FPD bond dissociation energies (BDEs) increase from 67 to 75 kcal/mol for x = 2 to x = 4 at the FPD level. The BDE for ThH 5 is much lower as it is a complex of H 2 with ThH 3 . The hydride affinities for x = 2 to 4 are all comparable and near 70 kcal/mol. A natural bond orbital analysis is consistent with a significant Th + –H– ionic contribution to the Th–H bonds. There is very little participation of the 5f orbitals in the bonding and the valence electrons on the Th are dominated by 7s and 6d for the neutrals and anions except for ThH 2 – where there is a significant contribution from the 7p.

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

Experimental and Computational Description of the Interaction of H and H – with U

In this work, the results of ab initio correlated molecular orbital theory electronic structure calculations for low-lying electronic states are presented for UH and UH – and compared to photoelectron spectroscopy measurements. The calculations were performed at the CCSD(T)/CBS and multireference CASPT2 including spin–orbit effects by the state interacting approach levels. The ground states of UH and UH – are predicted to be 4 Ι 9/2 and 5 Λ 6 , respectively. The spectroscopic parameters T e , r e , ω e , ω e x e , and Be were obtained, and potential energy curves were calculated for the low energy Ω states of UH. The calculated adiabatic electron affinity is 0.468 eV in excellent agreement with an experimental value of 0.462 ± 0.013 eV. The lowest vertical detachment energy was predicted to be 0.506 eV for the ground state, and the adiabatic ionization energy (IE) is predicted to be 6.116 eV. The bond dissociation energy (BDE) and heat of formation values of UH were obtained using the IE calculated at the Feller–Peterson–Dixon level. For UH, UH – , and UH + , the BDEs were predicted to be 225.5, 197.9, and 235.5 kJ/mol, respectively. The BDE for UH is predicted to be ~20% lower in energy than that for ThH. The analysis of the natural bond orbitals shows a significant U + H – ionic component in the bond of UH.

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

Simultaneous Optimization of Nuclear–Electronic Orbitals

Accurate modeling of important nuclear quantum effects, such as nuclear delocalization, zero-point energy, and tunneling, as well as non-Born-Oppenheimer effects, requires treatment of both nuclei and electrons quantum mechanically. The nuclear–electronic orbital (NEO) method provides an elegant framework to treat specified nuclei, typically protons, on the same level as the electrons. In conventional electronic structure theory, finding a converged ground state can be a computationally demanding task; converging NEO wavefunctions, due to their coupled electronic and nuclear nature, is even more demanding. Herein, we present an efficient simultaneous optimization method that uses the direct inversion in the iterative subspace method to simultaneously converge wavefunctions for both the electrons and quantum nuclei. In conclusion, benchmark studies show that the simultaneous optimization method can significantly reduce the computational cost compared to the conventional stepwise method for optimizing NEO wavefunctions for multicomponent systems.

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

Electronic Properties of UN and UN – from Photoelectron Spectroscopy and Correlated Molecular Orbital Theory

In this work, the results of calculations of the properties of the anion UN – including electron detachment are described, which further expand our knowledge of this diatomic molecule. High-level electronic structure calculations were conducted for the UN and UN – diatomic molecules and compared to photoelectron spectroscopy measurements. The low-lying Ω states were obtained using multireference CASPT2 including spin-orbit effects up to ~20,000 cm –1 . At the Feller–Peterson–Dixon (FPD) level, the adiabatic electron affinity (AEA) of UN is estimated to be 1.402 eV and the vertical detachment energy (VDE) is 1.423 eV. The assignment of the UN excited states shows good agreement with the experimental results with a VDE of 1.424 eV. An Ω = 4 ground state was obtained for UN – which is mainly associated with the 3 H ΛS state. Thermochemical calculations estimate a bond dissociation energy (BDE) for UN – (U – + N) of 665.9 kJ/mol, ~15% larger than that of UN and UN + . The NBO analysis reveals U–N triple bonds for the UN, UN – , and UN + species.

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

Theoretical and Experimental Study of the Spectroscopy and Thermochemistry of UC +/0/–

A combination of high-level ab initio calculations and anion photoelectron detachment (PD) measurements is reported for the UC, UC – , and UC + molecules. To better compare the theoretical values with the experimental photoelectron spectrum (PES), a value of 1.493 eV for the adiabatic electron affinity (AEA) of UC was calculated at the Feller–Peterson–Dixon (FPD) level. The lowest vertical detachment energy (VDE) is predicted to be 1.500 eV compared to the experimental value of 1.487 ± 0.035 eV. A shoulder to lower energy in the experimental PD spectrum with the 355 nm laser can be assigned to a combination of low-lying excited states of UC – and excited vibrational states. The VDEs calculated for the low-lying excited electronic states of UC at the SO-CASPT2 level are consistent with the observed additional electron binding energies at 1.990, 2.112, 2.316, and 3.760 eV. Potential energy curves for the Ω states and the associated spectroscopic properties are also reported. Compared to UN and UN + , the bond dissociation energy (BDE) of UC (411.3 kJ/mol) is predicted to be considerably lower. The natural bond orbitals (NBO) calculations show that the UC 0/+/– molecules have a bond order of 2.5 with their ground-state configuration arising from changes in the oxidation state of the U atom in terms of the 7s orbital occupation: UC (5f 2 7s 1 ), UC – (5f 2 7s 2 ), and UC + (5f 2 7s 0 ). Furthermore, the behavior of the UN and UC sequence of molecules and anions differs from the corresponding sequences for UO and UF.

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

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

Here, the thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller–Peterson–Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

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

Coupled Cluster Study of the Heats of Formation of UF 6 and the Uranium Oxyhalides, UO 2 X 2 (X = F, Cl, Br, I, and At)

The atomization enthalpies of the U(VI) species UF6 and the uranium oxyhalides UO 2 X 2 (X=F, Cl, Br, I, At) were calculated using a composite relativistic Feller-Peterson-Dixon (FPD) approach based on scalar relativistic DKH3-CCSD(T) with extrapolations to the CBS limit. The inherent multideterminantal nature of the U atom was mitigated by utilizing the singly charged atomic cation in all calculations with correction back to the neutral asymptote via the accurate ionization energy of the U atom. The effects of SO coupling were recovered using full 4-component CCSD(T) with contributions due to the Gaunt Hamiltonian calculated using Dirac-Hartree-Fock. The final atomization enthalpy for UF 6 (752.2 kcal/mol) was within 2.5 kcal/mol of the experimental value, but unfortunately the latter carries a ±2.4 kcal/mol uncertainty that is predominantly due to the experimental uncertainty in the formation enthalpy of U atom. The analogous value for UO 2 F 2 (607.6 kcal/mol) was in nearly exact agreement with experiment, but the latter has a stated experimental uncertainty of ±4.3 kcal/mol. The FPD atomization enthalpy for UO 2 Cl 2 (540.4 kcal/mol) was within the experimental error limits of ±5.5 kcal/mol. FPD atomization energies for the non-U-containing molecules (used for reaction enthalpies) H 2 O and HX (X=F, Cl, Br, I, At) were within at most 0.3 kcal/mol of their experimental values where available. The FPD atomization enthalpies, together with FPD reaction enthalpies for two different reactions, were used to determine heats of formation for all species of this work with estimated uncertainties of ±4 kcal/mol. The calculated heat of formation for UF 6 (-511.0 kcal/mol) is within 2.5 kcal/mol of the accurately-known (±0.45 kcal/mol) experimental value.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗