Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular-orbital-based machine learning for open-shell and multi-reference systems with kernel addition Gaussian process regression

We introduce a novel machine learning strategy, kernel addition Gaussian process regression (KA-GPR), in molecular-orbital-based machine learning (MOB-ML) to learn the total correlation energies of general electronic structure theories for closed- and open-shell systems by introducing a machine learning strategy. The learning efficiency of MOB-ML(KA-GPR) is the same as the original MOB-ML method for the smallest criegee molecule, which is a closed-shell molecule with multi-reference characters. In addition, the prediction accuracies of different small free radicals could reach the chemical accuracy of 1 kcal/mol by training on one example structure. Accurate potential energy surfaces for the H10 chain (closed-shell) and water OH bond dissociation (open-shell) could also be generated by MOB-ML(KA-GPR). To explore the breadth of chemical systems that KA-GPR can describe, we further apply MOB-ML to accurately predict the large benchmark datasets for closed- (QM9, QM7b-T, and GDB-13-T) and open-shell (QMSpin) molecules.

Chemistry↗
DOE OSTI2023

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

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

Fast uncertainty estimates in deep learning interatomic potentials

Deep learning has emerged as a promising paradigm to give access to highly accurate predictions of molecular and material properties. A common short-coming shared by current approaches, however, is that neural networks only give point estimates of their predictions and do not come with predictive uncertainties associated with these estimates. Existing uncertainty quantification efforts have primarily leveraged the standard deviation of predictions across an ensemble of independently trained neural networks. This incurs a large computational overhead in both training and prediction, resulting in order-of-magnitude more expensive predictions. Here, we propose a method to estimate the predictive uncertainty based on a single neural network without the need for an ensemble. This allows us to obtain uncertainty estimates with virtually no additional computational overhead over standard training and inference. We demonstrate that the quality of the uncertainty estimates matches those obtained from deep ensembles. We further examine the uncertainty estimates of our methods and deep ensembles across the configuration space of our test system and compare the uncertainties to the potential energy surface. Finally, we study the efficacy of the method in an active learning setting and find the results to match an ensemble-based strategy at order-of-magnitude reduced computational cost.

Chemistry↗
DOE OSTI2023

Generating accurate density matrices on the tangent space of a Grassmann manifold

Interpolating a density matrix from a set of known density matrices is not a trivial task. This is because a linear combination of density matrices does not necessarily correspond to another density matrix. In this Communication, density matrices are examined as objects of a Grassmann manifold. Although this manifold is not a vector space, its tangent space is a vector space. As a result, one can map the density matrices on this manifold to their corresponding vectors in the tangent space and then perform interpolations on that tangent space. The resulting interpolated vector can be mapped back to the Grassmann manifold, which can then be utilized (1) as an optimal initial guess for a self-consistent field (SCF) calculation or (2) to derive energy directly without time-consuming SCF iterations. Such a promising approach is denoted as Grassmann interpolation (G-Int). The hydrogen molecule has been used to illustrate that the described interpolated method in this work preserves the essential attributes of a density matrix. For phosphorus mononitride and ferrocene, it was demonstrated numerically that reference points for the definition of the corresponding tangent spaces can be chosen arbitrarily. In addition, the interpolated density matrices provide a superior and essentially converged initial guess for an SCF calculation to make the SCF procedure itself unnecessary. Finally, this accurate, efficient, robust, and systematically improved G-Int strategy has been used for the first time to generate highly accurate potential energy surfaces with fine details for the difficult case, ferrocene.

Chemistry↗
DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗
DOE OSTI2023

First principles simulation of reacting hypersonic flow over a blunt wedge

This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.

Mechanics↗
DOE OSTI2024

GradDFT. A software library for machine learning enhanced density functional theory

Density functional theory (DFT) stands as a cornerstone method in computational quantum chemistry and materials science due to its remarkable versatility and scalability. Yet, it suffers from limitations in accuracy, particularly when dealing with strongly correlated systems. To address these shortcomings, recent work has begun to explore how machine learning can expand the capabilities of DFT: an endeavor with many open questions and technical challenges. In this work, we present GradDFT a fully differentiable JAX-based DFT library, enabling quick prototyping and experimentation with machine learning-enhanced exchange–correlation energy functionals. GradDFT employs a pioneering parametrization of exchange–correlation functionals constructed using a weighted sum of energy densities, where the weights are determined using neural networks. Moreover, GradDFT encompasses a comprehensive suite of auxiliary functions, notably featuring a just-in-time compilable and fully differentiable self-consistent iterative procedure. To support training and benchmarking efforts, we additionally compile a curated dataset of experimental dissociation energies of dimers, half of which contain transition metal atoms characterized by strong electronic correlations. The software library is tested against experimental results to study the generalization capabilities of a neural functional across potential energy surfaces and atomic species, as well as the effect of training data noise on the resulting model accuracy.

Chemistry↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2024

The photoionization of methoxymethanol: Fingerprinting a reactive C 2 oxygenate in a complex reactive mixture

Methoxymethanol (CH 3 OCH 2 OH) is a reactive C 2 ether-alcohol that is formed by coupling events in both heterogeneous and homogeneous systems. It is found in complex reactive environments—for example those associated with catalytic reactors, combustion systems, and liquid-phase mixtures of oxygenates. Using tunable synchrotron-generated vacuum-ultraviolet photons between 10.0 and 11.5 eV, we report on the photoionization spectroscopy of methoxymethanol. We determine that the lowest-energy photoionization process is the dissociative ionization of methoxymethanol via H-atom loss to produce [C 2 H 5 O 2 ] + , a fragment cation with a mass-to-charge ratio (m/z) = 61.029. Here, we measure the appearance energy of this fragment ion to be 10.24 ± 0.05 eV. The parent cation is not detected in the energy range examined. To elucidate the origin of the m/z = 61.029 (C 2 H 5 O 2 ) fragment, we used automated electronic structure calculations to identify key stationary points on the cation potential energy surface and compute conformer-specific microcanonical rate coefficients for the important unimolecular processes. The calculated H-atom dissociation pathway results in a [C 2 H 5 O 2 ] + fragment appearance at 10.21 eV, in excellent agreement with experimental results.

Hansen, Niko A.↗
DOE OSTI2024

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1.

Chemistry↗
DOE OSTI2024

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. Here, we address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

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

Molecular dynamics studies of knotted polymers

Molecular dynamics calculations have been used to explore the influence of knots on the strength of a polymer strand. In particular, the mechanism of breaking 31, 41, 51, and 52 prime knots has been studied using two very different models to represent the polymer: (1) the generic coarse-grained (CG) bead model of polymer physics and (2) a state-of-the-art machine learned atomistic neural network (NN) potential for polyethylene derived from electronic structure calculations. While there is a broad overall agreement between the results on the influence of the pulling rate on chain rupture based on the CG and atomistic NN models, for the simple 31 and 41 knots, significant differences are found for the more complex 51 and 52 knots. Notably, in the latter case, the NN model more frequently predicts that these knots can break not only at the crossings at the entrance/exit but also at one of the central crossing points. The relative smoothness of the CG potential energy surface also leads to stabilization of tighter knots compared to the more realistic NN model.

DelloStritto, Mark (ORCID:0000000206785860)↗
DOE OSTI2025

Extending Badger's rule. I. The relationship between energy and structure in hydrogen bonds

Here, we derive a new expression for the strength of a hydrogen bond (V HB ) in terms of the elongation of the covalent bond of the donor fragment participating in the hydrogen bond (Δr HB ) and the intermolecular coordinates R (separation between the heavy atoms) and θ (deviation of the hydrogen bond from linearity). The expression includes components describing the covalent D–H bond of the hydrogen bond donor via a Morse potential, the Pauli repulsion, and electrostatic interactions between the constituent fragments using a linear expansion of their dipole moment and a quadratic expansion of their polarizability tensor. We fitted the parameters of the model using ab initio electronic structure results for six hydrogen bonded dimers, namely, NH 3 –NH 3 , H 2 O–H 2 O, HF–HF, H 2 O–NH 3 , HF–H 2 O, and HF–NH 3 , and validated its performance for extended parts of their potential energy surfaces, resulting in a mean absolute error ranging from 0.07 to 0.31 kcal/mol. The derived expression describes the energy–structure relationship in terms of a single structural parameter, namely, the elongation of the donor’s covalent bond (Δr HB ), and suggests the novel relationship of 8.0 kcal/mol pm −1 (or 0.8 kcal/mol per 0.001 Å elongation). This structural parameter is easily obtained from theory and can serve as the single descriptor of the strength of individual hydrogen bonds.

Santis, Garrett D. [Univ. of Washington, Seattle, ↗
DOE OSTI2026

Fragmentation dynamics of CS 2 dications and trications following S 2p ionization

Here, we present the results from a detailed study of the fragmentation dynamics of CS$_2^{2+}$ and CS$_2^{3+}$, formed in intense femtosecond soft x-ray pulses above the sulfur 2p edge, primarily through single core photoionization from the S 2p site, and subsequent Auger–Meitner decay(s). By combining three-dimensional velocity map imaging with covariance analysis, we determine the relative momenta of the ions produced in each two- and three-body fragmentation channel, at significantly higher ion count rates than conventional coincidence measurements. We shed new light on the wide range of fragmentation channels observed from the CS 2 dication and trication, including channels that involve ionization-induced bond formation and fragmentations producing undetected neutral cofragments. In the latter case, a “native frames” approach is used to isolate contributions from concerted and sequential fragmentations and extract dynamical information about each step of a concerted fragmentation process. While dications often fragment sequentially, the trication is dominated by concerted fragmentation. The main trication fragmentation channel into S + + C + + S + can be well-approximated by classical Coulombic simulations of the ground-state geometry distribution, reflecting both the nature of the trication potential energy surface and the rapid multiple ionization prior to substantial structural dynamics. This study demonstrates ways in which fundamental insights into the fragmentation dynamics of polycations following x-ray ionization may be extracted, which will be beneficial to future studies that employ time-resolved x-ray Coulomb explosion imaging to study ultrafast photochemistry.

Allum, Felix [Deutsches Elektronen-Synchrotron (DE↗
DOE OSTI2026

Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry

Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.

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

New reactions of diazene and related species for modelling combustion of amine fuels

Potential energy surfaces for reactions involving N 2 H 2 isomers of diazene (diimide) have been explored using density functional theory, with energies based on coupled-cluster theory. A focus is on processes that create or consume these species, and isomerisation between the E (trans) and Z (cis) forms of HNNH. Furthermore, these include isomerisation and dissociation pathways for HNNH, addition of H atoms to form N 2 H 3 , abstraction by H atoms yielding short-lived NNH, and abstraction reactions of H with N 2 H 3 . Transition state and capture theories are applied for high-pressure-limiting behaviour, while low-pressure and falloff regions are characterised via the methods of Troe and coworkers. Rate constants and thermochemistry are provided to improve models of diamine chemistry, relevant to the combustion of NH 3 especially at high concentrations, high pressures or under reducing conditions. Results indicate that amine radical recombination mainly yields the E HNNH isomer, while H-abstraction from N 2 H 3 results in E HNNH and H 2 NN. However, at elevated temperature E → Z isomerisation becomes competitive, and Z HNNH, being more reactive, acts to enhance the diazene consumption rate.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Observation of a new isomer in {sup 185}Au

Delayed γ-ray spectroscopy of 185Au was studied at the Argonne Gas-Filled Analyzer. A new isomer at an excitation energy of 1504.2(4) keV with a half-life of 630(80) ns was identified via γ-γ coincidence analysis, decaying via a 294.8(3) keV transition. Based on Weisskopf estimates, the multipolarity of the 295 keV transition is assigned to be E1, M1, E2, or M2. Possible configurations for this new isomer are discussed based on configuration-constrained potential energy surface calculations.

Fu, X. Y.↗