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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Properties of NoO and YbO

Bond dissociation energies (BDEs) and spectroscopic parameters for the ground states of YbO and NoO were calculated at the Brueckner doubles (BD(T)) level and for NoO at the coupled cluster CCSD(T) level based on the closed-shell configurations Yb 2+ (4f 14 )O 2 – (2p 6 ) and No 2+ (5f 14 )O 2 – (2p 6 ). For YbO, the BD(T) BDE (3.93 eV) and vibrational frequency (686.3 cm –1 ) are consistent with the experiment; the BD(T) bond distance was 1.8492 Å, about 0.04 Å longer than the experiment. For NoO, the calculated BDEs are 3.07 eV (BD(T)) and 3.08 eV (CCSD(T)) and are expected to be good to ±3 kcal/mol (±0.13 eV). Potential energy curves (PECs) for the spin–orbit (Ω) states associated with the f 13 s 1 configurations of YbO and NoO were calculated at the SA-CASSCF/SO-CASPT2/aQ-DK level. PECs considering both the f 13 s 1 and f 14 configurations were calculated at this level, as the ground states are expected to contain both configurations. Due to a bias in the SA-CASSCF to the f 13 s 1 configurations, the f 14 configurations could not be properly described simultaneously with the f 13 s 1 ones. Finally, the natural population analysis based on the natural bond orbitals shows that YbO and NoO are quite ionic, with 5d orbitals being more important in YbO than the 6d orbitals in NoO.

basis sets↗
DOE OSTI2025

Energetic and Electronic Properties of UX +/0/– for X = Li and Be and Comparison of the Properties of the Uranium Atom Binding to 2nd Row Elements Li–F

The bonding and spectroscopic properties of ULi +/0/– and UBe +/0/– to complete the series for UX +/0/– for X = Li to F were investigated by high-level ab initio SO-CASPT2 and CCSD(T) electronic structure calculations. The low-lying spin–orbit states were obtained at the SA-CASPT2/aQ-PP level; bond dissociation energies (BDEs), ionization energies (IEs), adiabatic electronic affinities (AEAs), and vertical detachment energies (VDEs) were calculated at the Feller-Peterson-Dixon (FPD) level. A dense manifold of low-lying states was predicted for ULi +/0/– and UBe +/0/– . Here, the calculated BDEs for ULi (37.7 kJ/mol) and UBe (8.0 kJ/mol) show that UBe is weakly bound. For redox processes, the BDEs increased for ULi + (109.3 kJ/mol), ULi – (47.4 kJ/mol), UBe + (35.6 kJ/mol), and UBe – (72.3 kJ/mol). The IE(ULi) = 4.650 eV is lower than IE(Li); the IE(UBe) = 5.901 eV is close to the IE(U) and to the IEs of UB, UC, UN, UO, and UF. The AEAs of ULi (0.708 eV) and UBe (0.989 eV) are lower than those for UB, UC, UN, and UO but higher than that for EA(UF). Natural bond orbital (NBO) calculations show that ULi has the 5f 3 6d 1 7s 2 configuration for U and 2s 1 for Li, with a small partial negative charge slightly delocalized on U. UBe arises from the U(5f 3 6d 1 7s 2 ) and Be(2s 2 ) electron configurations with no charge separation. The same calculations were made for WX (X = Li, Be, C–F) to enable detailed comparisons of the properties for UX with WX (X = Li–F). For WX, BDE(WX) is higher than that for UX for X = Li to N and lower than BDE(UX) for X = O and F, mostly due to the higher IE of W than U as ionic character becomes more important going from Li to F.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Low-Energy Isomers of the Magic Number H + (H 2 O) 21 Cluster

Electronic structure calculations are used to characterize low-energy isomers of H + (H 2 O) 21 . Eleven different classes of isomers, based on the (H 2 O) 20 pentagonal dodecahedron with the excess proton localized on the surface (as a hydrated hydronium ion) and the “extra” water molecule located in the interior of the cluster, are characterized. In 10 of these classes, the internal water molecule is engaged in six 5-membered rings, but in the remaining class, which is predicted to start at only 0.6 kcal/mol above the global minimum, the internal water is engaged in a 4- membered ring, an additional 6-membered ring, and four 5- membered rings. In addition, isomers with two 4-membered rings and two 6-membered rings on the cluster surface are predicted to start at only ∼1.3 kcal/mol above the lowest-energy dodecahedralbased structure.

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

Chemist: A Domain-Specific Language by Chemists for Chemists

Managing the complexity of quantum chemistry (QC) software is key to ensuring it remains accessible, maintainable, and reusable. Noticeably missing from the QC ecosystem are modules targeting bottleneck routines. Here we argue that this is likely due to the difficulty in defining interfaces for such modules. To that end, we introduce the open-source, publicly available Chemist library https://github.com/NWChemEx/Chemist. Chemist is a domain-specific language targeting the QC domain. Chemist has been developed focusing on performance and user-friendliness. Using Chemist, QC tasks are defined using familiar domain concepts such as molecules, wave functions, and operators. The domain objects are hierarchical to ensure a systematic encapsulation of information. Key features of Chemist include: extensibility, the ability to alias existing data, and the ability to succinctly define many common QC tasks. The usefulness of Chemist is demonstrated by discussing the interface of NWChemEx’s Fock build module and by showcasing a proof-of-concept self-consistent field algorithm containing uncertainty propagation.

Algorithms↗
DOE OSTI2025

Understanding the Existence of a Na 2 Dimer in a High-Spin State

The recent observation of a high-spin Na 2 dimer formed on the surface of liquid helium nanodroplets raises some fundamental questions, as the ground state of Na 2 is known to have zero spin. Is it protected against spontaneous dissociation? What is its binding energy and interatomic distance? Is it stable at a higher temperature? Using calculations based on density functional theory (with and without long-range interaction) and coupled cluster methods, CCSD(T), we show that the bonding in the high-spin Na 2 dimer is governed by van der Waals interaction with binding energy (bond length) varying between −0.030 eV (5.108 Å) and −0.192 eV (4.231 Å), depending on the computational method used. Thus, the experimental method used by Kresin and coworkers can be very useful to study larger metastable high-spin clusters such as Li 4 which was predicted in 1985 to have a tetrahedral structure carrying a magnetic moment of 2 μ B , while its ground state is planar and nonmagnetic.

Basis sets↗
DOE OSTI2025

The Singlet–Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC( P ; Q ) Study

An accurate determination of singlet−triplet gaps in biradicals, including cyclobutadiene in the automerization barrier region where one has to balance the substantial nondynamical many-electron correlation effects characterizing the singlet ground state with the predominantly dynamical correlations of the lowest-energy triplet, remains a challenge for many quantum chemistry methods. High-level coupled-cluster (CC) approaches, such as the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT), are often capable of providing reliable results, but routine application of such methods is hindered by their high computational costs. We have recently proposed a practical alternative to converging the CCSDT energetics at small fractions of the computational effort, even when electron correlations become stronger and connected triply excited clusters are larger and nonperturbative, by merging the CC(P;Q) moment expansions with the selected configuration interaction methodology abbreviated as CIPSI. We demonstrate that one can accurately approximate the highly accurate CCSDT potential surfaces characterizing the lowest singlet and triplet states of cyclobutadiene along the automerization coordinate and the gap between them using tiny fractions of triply excited cluster amplitudes identified with the help of relatively inexpensive CIPSI Hamiltonian diagonalizations.

Basis sets↗
DOE OSTI2023

Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory

The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density-matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies, relative to reference values from configuration-interaction-based calculations.

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

Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

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

Performance of Diffusion Monte Carlo Calculations for Predicting the Relative Energies of Quinoidal and Nonquinoidal Species

Coupled cluster singles and doubles with perturbative triples [CCSD(T)] and single determinant fixed-node diffusion Monte Carlo (SD-DMC) have emerged as two of the most useful methods for providing benchmark reaction and interaction energies of chemical systems without strong static correlation. The errors in DMC energies are dominated by an inexact description of the nodal surfaces for electron exchange. One of the main approaches to addressing the fixed-node error is to use multideterminant (MD) trial wave functions. We consider here the energy differences between pairs of related molecules with aromatic and quinoidal structures as well as between quinoidal isomers. Quinoidal systems tend to have some diradical character, leading one to anticipate that SD-DMC calculations may face challenges in accurately describing their energetics. The MD trial wave functions were generated from the complete active space calculations. A comparison is made with the predictions of well-converged CCSD(T) calculations.

basis sets↗
DOE OSTI2025

Tests of the DFT Ladder for the Fulminic Acid Challenge

Properties of the historically pivotal fulminic acid (HCNO) molecule have been computed with a panoply of 473 density functionals of all varieties, providing a snapshot of the performance of contemporary density functional theory (DFT) for a challenging chemical system. Exhaustive tabulations and statistical analyses have been carried out for geometric parameters, vibrational frequencies, barriers to linearity, and the HCN–O dissociation energy. As the DFT ladder is climbed, confusion rather than consensus ensues regarding the details of the distinctive, extremely flat H–C–N bending potential of fulminic acid and whether the equilibrium structure is linear or bent. While high-ranking DFT functionals produce the smallest errors for the HCN + O( 3 P) → HCNO reaction energy, lower rungs emerge as the best performers for many of the bond distances and harmonic vibrational frequencies. This research shows that the current DFT zoo of approximations does not constitute a transparent ladder of increasingly accurate methods that consistently converges on definitive predictions for various properties of HCNO. Additional analyses are performed on the side effects of popular dispersion corrections on the covalently bonded properties and thermochemistry of HCNO.

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

Synchronization in electric power networks with inherent heterogeneity up to 100% inverter-based renewable generation

The synchronized operation of power generators is the foundation of electric power network stability and a key to the prevention of undesired power outages and blackouts. Here, we derive the conditions that guarantee synchronization in power networks with inherent generator heterogeneity when subjected to small perturbations, and perform a parametric sensitivity analysis to understand synchronization with varied types of generators. As inverter-based resources, which are the primary interfacing technology for many renewable sources of energy, have supplanted synchronous generators in ever growing numbers, the center of attention on associated integration challenges have resided primarily on the role of declining system inertia. Our results instead highlight the critical role of generator damping in achieving a stable state of synchronization. Additionally, we report the feasibility of operating interconnected electric grids with up to 100% power contribution from inverter-based renewable generation technologies. Our study has important implications as it sets the basis for the development of advanced control architectures and grid optimization methods that ensure synchronization and further pave the path towards the decarbonization of the electric power sector.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Lossy compression of statistical data using quantum annealer

Abstract We present a new lossy compression algorithm for statistical floating-point data through a representation learning with binary variables. The algorithm finds a set of basis vectors and their binary coefficients that precisely reconstruct the original data. The optimization for the basis vectors is performed classically, while binary coefficients are retrieved through both simulated and quantum annealing for comparison. A bias correction procedure is also presented to estimate and eliminate the error and bias introduced from the inexact reconstruction of the lossy compression for statistical data analyses. The compression algorithm is demonstrated on two different datasets of lattice quantum chromodynamics simulations. The results obtained using simulated annealing show 3–3.5 times better compression performance than the algorithm based on neural-network autoencoder. Calculations using quantum annealing also show promising results, but performance is limited by the integrated control error of the quantum processing unit, which yields large uncertainties in the biases and coupling parameters. Hardware comparison is further studied between the previous generation D-Wave 2000Q and the current D-Wave Advantage system. Our study shows that the Advantage system is more likely to obtain low-energy solutions for the problems than the 2000Q.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Enhanced descriptor identification and mechanism understanding for catalytic activity using a data-driven framework: revealing the importance of interactions between elementary steps

We report accurate identification of descriptors for catalytic activities has long been essential to the in-depth understanding of catalysis and recently to set the basis for catalyst screening. However, commonly used methods suffer from low accuracy in predictability. This study reports an enhanced approach to accurately identify the descriptors from a kinetic dataset using a machine learning (ML) surrogate model. CO hydrogenation to methanol over Cu-based catalysts was taken as a case study. Our model captures not only the contribution from individual elementary steps but also the interaction between relevant steps within a reaction network, which was found to be essential for high accuracy. As a result, six effective descriptors are identified, which are accurate enough to ensure the trained gradient boosted regression (GBR) model for good prediction of the methanol turnover frequency (TOF) over metal (M)-doped Cu(111) model surfaces (M = Au, Cu, Pd, Pt, Ni). More importantly, going beyond the purely mathematical ML model, the catalytic role of each identified descriptor can be revealed by using model-agnostic interpretation tools, which enhances the insight into the promoting effect of alloying. The trained GBR model outperforms the conventional derivative-based methods in terms of both the predictability and the mechanism understanding. It opens alternative possibilities toward accurate descriptor-based rational catalyst optimization.

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

An analytic asymmetric-piston model for the impact of mode-1 shell asymmetry on ICF implosions

For many years, low mode asymmetry in inertially confined fusion (ICF) implosions has been recognized as a potential performance limiting factor, but analysis has been limited to using simulations and searching for data correlations. In this paper, an analytically solvable model based upon the simple picture of an asymmetric piston is presented. Asymmetry of the shell driving the implosion, as opposed to asymmetry in the hot-spot, is key to the model. The model provides a unifying framework for the action of mode-1 shell asymmetry and the resulting connections between various diagnostic signatures. A key variable in the model is the shell asymmetry fraction, f, which is related to the areal density variation of the shell surrounding the hot-spot. It is shown that f is simply related to the observed hot-spot mode-1 velocity and to the concept of residual energy in an implosion. The model presented in this paper yields explicit expressions for the hot-spot diameter, stagnation pressure, hot-spot energy, inertial confinement-time, Lawson parameter, hot-spot temperature, and fusion yield under the action of mode-1 asymmetry. Agreement is found between the theory scalings when compared to ICF implosion data from the National Ignition Facility and to large ensembles of detailed simulations, making the theory a useful tool for interpreting data. The theory provides a basis for setting tolerable limits on asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗