Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

ab initio Sub-Mechanism Development for Cyclopentene Oxidation

To accurately predict low-temperature oxidation behavior, chemical kinetics mechanisms must contain complete reaction networks that include detailed consumption reactions of intermediates produced directly from hydroperoxyalkyl radicals, Q̇OOH, which undergo competing unimolecular reactions and bimolecular reactions with O2. Rates of chain-branching are governed by the flux between the two competing pathways, and inherently depend on temperature, pressure, and oxygen concentration. Neglect of consumption pathways for major oxidation intermediates leads to mechanism truncation error that is ameliorated by expanding the level of detail included in sub-mechanisms and employing ab initio methods for computing rates of elementary reactions and thermochemical properties of species involved. In the present work, an ab initio-derived sub-mechanism is developed using AutoMech to model the chemical kinetics of cyclopentene, a major product of cyclopentane oxidation. The ab initio sub-mechanism builds on a detailed mechanism developed using Reaction Mechanism Generator (RMG) for the specific purpose of determining the extent to which replacing cyclopentene-specific reactions and species with quantum chemical computations reduces model inaccuracies resulting from mechanism truncation error. In an effort to minimize interference from other reactions present during the formation of cyclopentene from cyclopentyl + O2, providing a narrower experimental scope, the model is compared against speciation measurements from jet-stirred reactor (JSR) experiments on cyclopentene oxidation. The experiments utilize vacuum ultraviolet-absorption spectroscopy and mass spectrometry for isomer-resolved speciation of intermediates at 835 Torr from 700 – 950 K. [O2]-dependent experiments were also conducted from 0.057 – 2.01 · 1018 molecules cm–3 at 825 K to examine the influence of oxygen on species profiles. Model predictions using the ab initio-revised mechanism yielded significant improvements in species profiles for both the temperature- and [O2]-dependent measurements, owing in part to increased rates of HOȮ and H2O2 production, which underscores the influence of theoretical calculations of reaction rates involving species produced from Ṙ + O2 such as cyclopentene.

AutoMech↗
DOE OSTI2019

Muon capture in nuclei: An ab initio approach based on Green's function Monte Carlo methods

An ab initio Green's function Monte Carlo (GFMC) method is introduced for calculating total rates of muon weak capture in light nuclei with mass number $A ≤ 12$. As a first application of the method, we perform a calculation of the rate in $^3$$\textbf{H}$ and $^4$$\textbf{He}$ in a dynamical framework based on realistic two- and three-nucleon interactions and realistic nuclear charge-changing weak currents. The currents include one- and two-body terms induced by $π$- and $ρ$-meson exchange, and $N$-to-$Δ$ excitation, and are constrained to reproduce the empirical value of the Gamow-Teller matrix element in tritium. We investigate the sensitivity of theoretical predictions to current parametrizations of the nucleon axial and induced pseudoscalar form factors as well as to two-body contributions in the weak currents. The large uncertainties in the measured $^4$$\textbf{He}$ rates obtained from bubble-chamber experiments (carried out over 50 years ago) prevent us from drawing any definite conclusions. No data exist for $^3$$\textbf{H}$, but results are compared to those of a recent Faddeev calculation as a validation of the present GFMC method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

pnnl/mol_dgnn

Molecular dynamics simulations often rely on parameterized atomistic potentials, which can be inaccurate compared to ab initio simulation methods. However, dynamic ab initio simulations come with high computational cost. Deep learning techniques offer a method to obtain both high accuracy and low cost predictions of dynamics. Our method gives a prediction of the structural dynamics of a molecule accounting for the prior dynamics, giving high-accuracy predictions at low computational cost

Ashby, Michael↗
DOE OSTI2023

Optical potentials for the rare-isotope beam era

Abstract We review recent progress and motivate the need for further developments in nuclear optical potentials that are widely used in the theoretical analysis of nucleon elastic scattering and reaction cross sections. In regions of the nuclear chart away from stability, which represent a frontier in nuclear science over the coming decade and which will be probed at new rare-isotope beam facilities worldwide, there is a targeted need to quantify and reduce theoretical reaction model uncertainties, especially with respect to nuclear optical potentials. We first describe the primary physics motivations for an improved description of nuclear reactions involving short-lived isotopes, focusing on its benefits for fundamental science discoveries and applications to medicine, energy, and security. We then outline the various methods in use today to build optical potentials starting from phenomenological, microscopic, and ab initio methods, highlighting in particular, the strengths and weaknesses of each approach. We then discuss publicly-available tools and resources facilitating the propagation of recent progresses in the field to practitioners. Finally, we provide a set of open challenges and recommendations for the field to advance the fundamental science goals of nuclear reaction studies in the rare-isotope beam era. This paper is the outcome of the Facility for Rare Isotope Beams Theory Alliance (FRIB-TA) topical program ‘Optical Potentials in Nuclear Physics’ held in March 2022 at FRIB. Its content is non-exhaustive, was chosen by the participants and reflects their efforts related to optical potentials.

42 ENGINEERING↗
DOE OSTI2024

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Simulated XUV photoelectron spectra of THz-pumped liquid water

Highly intense, sub-picosecond terahertz (THz) pulses can be used to induce ultrafast temperature jumps (T-jumps) in liquid water. A supercritical state of gas-like water with liquid density is established, and the accompanying structural changes are expected to give rise to time-dependent chemical shifts. We investigate the possibility of using extreme ultraviolet photoelectron spectroscopy as a probe for ultrafast dynamics induced by sub-picosecond THz pulses of varying intensities and frequencies. To this end, we use ab initio methods to calculate photoionization cross sections and photoelectron energies of (H 2 O) 20 clusters embedded in an aqueous environment represented by point charges. The cluster geometries are sampled from ab initio molecular dynamics simulations modeling the THz-water interactions. Here we find that the peaks in the valence photoelectron spectrum are shifted by up to 0.4 eV after the pump pulse and that they are broadened with respect to unheated water. The shifts can be connected to structural changes caused by the heating, but due to saturation effects they are not sensitive enough to serve as a thermometer for T-jumped water.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗
DOE OSTI2020

Nature of Hyperfine Interactions in TbPc 2 Single-Molecule Magnets: Multiconfigurational Ab Initio Study

Lanthanide-based single-ion magnetic molecules can have large magnetic hyperfine interactions as well as large magnetic anisotropy. Recent experimental studies reported tunability of these properties by changes of chemical environments or by application of external stimuli for device applications. In order to provide insight onto the origin and mechanism of such tunability, here we investigate the magnetic hyperfine and nuclear quadrupole interactions for a 159 Tb nucleus in TbPc 2 (Pc = phthalocyanine) single-molecule magnets using multiconfigurational ab initio methods including spin–orbit interaction. Since the electronic ground and first-excited (quasi)doublets are well separated in energy, the microscopic Hamiltonian can be mapped onto an effective Hamiltonian with an electronic pseudospin S = 1/2. From the ab initio calculated parameters, we find that the magnetic hyperfine coupling is dominated by the interaction of the Tb nuclear spin with electronic orbital angular momentum. The asymmetric 4f-like electronic charge distribution leads to a strong nuclear quadrupole interaction with significant transverse terms for the molecule with low symmetry. The ab initio calculated electronic–nuclear spectrum including the magnetic hyperfine and quadrupole interactions is in excellent agreement with the experiment. We further find that the transverse quadrupole interactions significantly influence the avoided level crossings in magnetization dynamics and that the molecular distortions affect mostly the Fermi contact terms as well as the transverse quadrupole interactions. Synopsis Hyperfine interactions for 159 Tb nucleus in TbPc 2 single-molecule magnets are investigated from first-principles using multiconfigurational calculations. Strong nuclear quadrupole and magnetic hyperfine coupling are found, with the latter being dominated by the paramagnetic spin-orbital mechanism. Furthermore, we construct an ab initio pseudospin Hamiltonian and obtain the electronic-nuclear spectrum that is in excellent agreement with experiment. The Zeeman diagram is calculated, and the magnetization dynamics is discussed. The effects of molecular distortions are researched.

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

Deep learning of dynamically responsive chemical Hamiltonians with semiempirical quantum mechanics

Conventional machine-learning (ML) models in computational chemistry learn to directly predict molecular properties using quantum chemistry only for reference data. While these heuristic ML methods show quantum-level accuracy with speeds several orders of magnitude faster than traditional quantum chemistry methods, they suffer from poor extensibility and transferability; i.e., their accuracy degrades on large or new chemical systems. Incorporating quantum chemistry frameworks into the ML models directly solves this problem. Here we take the structure of semiempirical quantum mechanics (SEQM) methods to construct dynamically responsive Hamiltonians. SEQM methods use empirical parameters fitted to experimental properties to construct reduced-order Hamiltonians, facilitating much faster calculations than ab initio methods but with compromised accuracy. By replacing these static parameters with machine-learned dynamic values inferred from the local environment, we greatly improve the accuracy of the SEQM methods. Trained on molecular energies and atomic forces, these dynamically generated Hamiltonian parameters show a strong correlation with atomic hybridization and bonding. Trained with only about 60,000 small organic molecular conformers, the resulting model retains interpretability, extensibility, and transferability when testing on much larger chemical systems and predicting various molecular properties. Overall, this work demonstrates the virtues of incorporating physics-based descriptions with ML to develop models that are simultaneously accurate, transferable, and interpretable.

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

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2025

Python Library for Monte Carlo Simulations with Ab Initio and Machine-Learned Interatomic Potentials

There is a growing need in the simulation community for software that provides a transparent, reproducible, usable, and extensible (TRUE) Monte Carlo (MC) simulation framework employing energies from ab initio methods and machine-learning interatomic potentials (MLIPs). We introduce a Python library (ASE-MC) that adds Monte Carlo functionality to the Atomic Simulation Environment (ASE) package. Now, we can combine the powerful tools used to build systems and perform ab initio and MLIP in ASE with MC simulation algorithms to sample the configurational space with a concise Python script. After presenting the design philosophy, we demonstrate the flexibility of our approach using selected examples. These example simulations include liquid water described with a message-passing MLIP in the canonical and isothermal–isobaric ensembles, sampling the characteristic dihedral angle of biphenyl and comparing an MLIP to first-principles calculations, and a grand canonical Monte Carlo simulation of ammonia adsorption on Pt(111). These examples showcase the main features of the software, which include flexibility in the choice of ab initio or MLIP engine, ab initio or MLIP grand canonical MC with cavity bias insertions and deletions, the ability to add custom MC moves to the move set, and how users can condense complex MC workflows into a single Python script. Finally, this library serves as a framework for reproducible Monte Carlo simulations, facilitating easy reproduction of the work and application to new systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Ab Initio Computation of Charge Densities for Sn and Xe Isotopes

We present the first ab initio calculations for open-shell nuclei past the tin isotopic line, focusing on Xe isotopes as well as doubly magic Sn isotopes. We show that, even for moderately hard interactions, it is possible to obtain meaningful predictions and that the NNLO sat chiral interaction predicts radii and charge density distributions close to the experiment. We then make a new prediction for 100 Sn. This paves the way for ab initio studies of exotic charge density distributions at the limit of the present ab initio mass domain, where experimental data is becoming available. The present study closes the gap between the largest isotopes reachable by ab initio methods and the smallest exotic nuclei accessible to electron scattering experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

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

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

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

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Benchmarking the performance of uncertainty quantification methods for neural network-based interatomic potentials

Machine-learned interatomic potentials (ML-IAPs) continue to gain popularity as accurate, computationally efficient replacements for traditional, physics-based interatomic potentials and expensive ab initio methods. Uncertainty quantification (UQ) of ML-IAPs is a growing area of research as UQ is critical in many applications of IAPs, such as developing curated datasets, active learning-based data augmentation, self-improving models, and estimating the uncertainty of molecular dynamics simulations. In this paper, we construct and benchmark a series of different neural network potentials (NNPs) with varying network architectures to determine the performance of these models with respect to both the mean and uncertainty calibration error. Each NNP method is specifically designed to predict either epistemic or aleatoric uncertainty with particular focus on the differences in behavior between the epistemic and aleatoric uncertainty estimates. We benchmark these methods using multiple datasets common in the ML-IAP literature. The results show that the aleatoric uncertainty from single-shot model architectures is a competitive alternative to ensemble-based epistemic uncertainty predictions in regions of sufficient data-density. However, in regions where the representative data is sparse, aleatoric uncertainty models tend to overpredict and epistemic methods tend to underpredict the actual model error. We conclude that the type of UQ is crucial when discussing performance of probabilistic model results as different methods have different performance characteristics depending on the regime in which they are evaluated. Therefore, the type of UQ method should be carefully evaluated against both the data characteristics and requirements for the intended application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Be 7 and Li 7 nuclei within the no-core shell model with continuum

The production of $^7Be$ and $^7Li$ nuclei plays an important role in primordial nucleosynthesis, nuclear astrophysics, and fusion energy generation. The $^3He(α,γ) ^7Be$ and $^3H(α,γ) 7Li$ radiative-capture processes are important to determine the $^7Li$ abundance in the early universe and to predict the correct fraction of $pp$ -chain branches resulting in $^7Be$ versus $^8B$ neutrinos. The $^6Li(p,γ) ^7Be$ has been investigated recently hinting at a possible cross section enhancement near the threshold. The $^6Li (n,^3H) ^4He$ process can be utilized for tritium breeding in machines dedicated to fusion energy generation through the deuteron-tritium reaction, and is a neutron cross section standard used in the measurement and evaluation of fission cross sections. In this work we study the properties of $^7Be$ and $^7Li$ within the no-core shell model with continuum (NCSMC) method, using chiral nucleon-nucleon interactions as the only input, and analyze all the binary mass partitions involved in the formation of these systems. The NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states and thus is well suited to investigate systems with many resonances and pronounced clustering like $^7Be$ and $^7Li$. Our calculations reproduce all the experimentally known states of the two systems and provide predictions for several new resonances of both parities. Some of these new possible resonances are built on the ground states of $^6Li$ and $^6He$, and thus represent a robust prediction. We do not find any resonance in the $p + ^6Li$ mass partition near the threshold. On the other hand, in the $p + ^6He$ mass partition of $^7Li$ we observe an $S$ -wave resonance near the threshold producing a very pronounced peak in the calculated $S$ factor of the $^6He (p,γ)$ $^7Li$ radiative-capture reaction. While we do not find a resonance near the threshold in the $p + ^6Li$ channel, in the case of $^6He + p$ reaction a resonant S -wave state is predicted at a very low energy above the reaction threshold, which could be relevant for astrophysics and its implications should be investigated. We note though that this state lies above the three-body breakup threshold not included in our method and may be influenced by three-body continuum correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗