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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

State-of-the-art ab initio potential energy curve for the krypton atom pair and thermophysical properties of dilute krypton gas

A new reference krypton-krypton interatomic potential energy curve was developed by means of quantum-chemical ab initio calculations for 36 interatomic separations. Highly accurate values for the interaction energies at the complete basis set limit were obtained using the coupled-cluster method with single, double, and perturbative triple excitations as well as t-aug-cc-pV5Z and t-aug-cc-pV6Z basis sets including mid-bond functions, with the 6Z basis set being newly constructed for this study. Higher orders of coupled-cluster terms were considered in a successive scheme up to full quadruple excitations. Core-core and core-valence correlation effects were included. Furthermore, relativistic effects were studied not only at a scalar relativistic level using second-order direct perturbation theory, but also utilizing full four-component and Gaunt-effect computations. An analytical pair potential function was fitted to the interaction energies, which is characterized by a depth of 200.88 K with an estimated standard uncertainty of 0.51 K. Thermophysical properties of low-density krypton were calculated for temperatures up to 5000 K. Second and third virial coefficients were obtained from statistical thermodynamics. Viscosity and thermal conductivity as well as the self-diffusion coefficient were computed using the kinetic theory of gases. The theoretical results are compared with experimental data and with results for other pair potential functions from the literature, especially with those calculated from the recently developed ab initio potential of Waldrop et al. [J. Chem. Phys. 142, 204307 (2015)]. Highly accurate experimental viscosity data indicate that both the present ab initio pair potential and the one of Waldrop et al. can be regarded as reference potentials, even though the quantum-chemical methods and basis sets differ. However, the uncertainties of the present potential and of the derived properties are estimated to be considerably lower.

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

Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application

Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A Holistic Approach for Elucidating Local Structure, Dynamics, and Speciation in Molten Salts with High Structural Disorder

To examine ion solvation, exchange, and speciation for minority components in molten salts (MS) typically found as corrosion products, we propose a multimodal approach combining extended X-ray absorption fine structure (EXAFS) spectroscopy, optical spectroscopy, ab initio molecular dynamics (AIMD) simulations, and rate theory of ion exchange. Going beyond conventional EXAFS analysis, our method can accurately quantify populations of different coordination states of ions with highly disordered coordination environments via linear combination fitting of the EXAFS spectra of these coordination states computed from AIMD to the experimental EXAFS spectrum. In a case study of dilute Ni(II) dissolved in the ZnCl 2 +KCl melts, our method reveals heterogeneous distributions of coordination states of Ni(II) that are sensitive to variations in temperature and melt composition. These results are fully explained by the difference in the chloride exchange dynamics at varied temperatures and melt compositions. Finally, this insight will enable a better understanding and control of ion solubility and transport in MS.

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

Development of interatomic potential for Al–Tb alloys using a deep neural network learning method

An interatomic potential for the Al–Tb alloy around the composition of Al 90 Tb 10 is developed using the deep neural network (DNN) learning method. The atomic configurations and the corresponding total potential energies and forces on each atom obtained from ab initio molecular dynamics (AIMD) simulations are collected to train a DNN model to construct the interatomic potential for the Al–Tb alloy. Here we show that the obtained DNN model can well reproduce the energies and forces calculated by AIMD simulations. Molecular dynamics (MD) simulations using the DNN interatomic potential also accurately describe the structural properties of the Al 90 Tb 10 liquid, such as partial pair correlation functions (PPCFs) and bond angle distributions, in comparison with the results from AIMD simulations. Furthermore, the developed DNN interatomic potential predicts the formation energies of the crystalline phases of the Al–Tb system with an accuracy comparable to ab initio calculations. The structure factors of the Al 90 Tb 10 metallic liquid and glass obtained by MD simulations using the developed DNN interatomic potential are also in good agreement with the experimental X-ray diffraction data. The development of short-range order (SRO) in the Al 90 Tb 10 liquid and the undercooled liquid is also analyzed and three dominant SROs, i.e., Al-centered distorted icosahedron (DISICO) and Tb-centered ‘3661’ and ‘15551’ clusters, respectively, are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Peierls transition, ferroelectricity, and spin-singlet formation in monolayer VOI 2

Using ab initio density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI 2 with a 3d 1 electronic configuration. Our phonon calculations indicate that the orthorhombic Pmm2 FE-II phase is the most likely ground state, involving a ferroelectric (FE) distortion along the a axis and V-V dimerization along the b axis. Specifically, the “pseudo Jahn-Teller” effect caused by the coupling between empty V (d xz/yz and d 3z 2 -r 2 ) and O 2p states is proposed as the mechanism that stabilizes the FE distortion from the paraelectric phase. Moreover, the half-filled metallic d xy band displays a Peierls instability along the b axis, inducing a V-V dimerization. Additionally, we found very short-range antiferromagnetic coupling along the V-V chain due to the formation of nearly decoupled spin singlets in the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

NuLattice: Ab initio computations of atomic nuclei on lattices

Here, we introduce NuLattice, a Python software package for ab initio computations of atomic nuclei on lattices. The computational tools consist of Hartree Fock, the coupled cluster method, the in-medium similarity renormalization group, and full configuration interaction. At present, the employed interactions are from pion-less effective field theory at leading order and consist of two-body and three-body contacts. We present results for light nuclei 2 H, 3,4 He, 8 Be, 12 C, and 16 O. NuLattice algorithms exploit the sparsity and locality of lattice interactions, and as a result computations can be run on laptops.

Rothman, Maxwell [Univ. of Tennessee, Knoxville, T↗
DOE OSTI2017

Theoretical Investigation of Kinetic Processes in Small Radicals of Importance in Combustion

Our group studies inelastic and reactive collisions of small molecules, focusing on radicals important in combustion environments. The goal is the better understanding of kinetic processes that may be difficult to access experimentally. An essential component is the accurate determination and fitting of potential energy surfaces (PESs). After fitting the ab initio points to obtain global PESs, we treat the dynamics using time-independent (close-coupling) methods. Cross sections and rate constants for collisions of are determined with our Hibridon program suite . We have studied energy transfer (rotationally, vibrationally, and/or electronically inelastic) in small hydrocarbon radicals (CH 2 and CH 3 ) and the CN radical. We have made a comparison with experimental measurements of relevant rate constants for collisions of these radicals. Also, we have calculated accurate transport properties using state-of-the-art PESs and to investigate the sensitivity to these parameters in 1-dimensional flame simulations. Of particular interest are collision pairs involving the light H atom.

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

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light–matter hybrid states, so-called polaritons, are formed due to quantum light–matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light–matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light–matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule–cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

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

Dynamical mean field theory simulations with the adaptive sampling configuration interaction method

In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero-temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. Here, we present the ASCI-DMFT method and example calculations on the one- and two-dimensional Hubbard models that exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground-state DMFT simulations, which allows us to study the bath discretization error in simulations with small clusters, as well as to address cluster sizes beyond the current state of the art. Finally, our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Site-averaged kinetics for catalysts on amorphous supports: an importance learning algorithm

Ab initio calculations have greatly advanced our understanding of homogeneous catalysts and crystalline heterogeneous catalysts. In contrast, amorphous heterogeneous catalysts remain poorly understood. The principal difficulties include (i) the nature of the disorder is quenched and unknown; (ii) each active site has a different local environment and activity; (iii) active sites are rare, often less than ~20% of potential sites, depending on the catalyst and its preparation method. Few (if any) studies of amorphous heterogeneous catalysts have ever attempted to compute site-averaged kinetics, because the exponential dependence on variable activation energy requires an intractable number of ab initio calculations to converge. We present a new algorithm using machine learning techniques (metric learning kernel regression) and importance sampling to efficiently learn the distribution of activation energies. We demonstrate the algorithm by computing the site-averaged activity for a model amorphous catalyst with quenched disorder.

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

An Ab Initio -Derived Force Field for Amorphous Silica Interfaces for Use in Molecular Dynamics Simulations

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

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

Quantum chemical ab initio prediction of proton exchange barriers between CH{sub 4} and different H-zeolites

A hybrid MP2:DFT (second-order Møller–Plesset perturbation theory–density functional theory) method that combines MP2 calculations for cluster models with DFT calculations for the full periodic structure is used to localize minima and transition structures for proton jumps at different Brønsted sites in different frameworks (chabazite, faujasite, ferrierite, and ZSM-5) and at different crystallographic positions of a given framework. The MP2 limit for the periodic structures is obtained by extrapolating the results of a series of cluster models of increasing size. A coupled-cluster (CCSD(T)) correction to MP2 energies is calculated for cluster models consisting of three tetrahedra. For the adsorption energies, this difference is small, between 0.1 and 0.9 kJ/mol, but for the intrinsic proton exchange barriers, this difference makes a significant (10.85 ± 0.25 kJ/mol) and almost constant contribution across different systems. The total values of the adsorption energies vary between 22 and 34 kJ/mol, whereas the total proton exchange energy barriers fall in the narrow range of 152–156 kJ/mol. After adding nuclear motion contributions (harmonic approximation, 298 K), intrinsic enthalpy barriers between 134 and 141 kJ/mol and apparent energy barriers between 105 and 118 kJ/mol are predicted for the different sites examined for the different frameworks. These predictions are consistent with experimental results available for faujasite, ferrierite, and ZSM-5.

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

Ab-initio study of thermoelectric properties of Mg{sub 2}Ge

In this paper we investigate the thermoelectric properties of Mg{sub 2}Ge material using first principles pseudo potential method based on density functional theory and Boltzmann transport equations. The calculations show n-type conduction, indicating that the electrical conduction is due to electron. The electrical conductivity decrease with temperature; the negative value of Seebeck Coefficient also shows that the conduction is due to electron. In this paper we have calculated Seebeck coefficient, electrical conductivity and thermal conductivity. The thermoelectric properties of system have been calculated in this temperature range 100 K-1200 K. The value of Figure of Merit (ZT) of this material is also increasing with temperature and maximum value is 4.58×10{sup −5}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Research in Lattice Gauge Theory and in the Phenomenology of Neutrinos and Dark Matter

Research in theoretical elementary particle physics was performed by the PI Yannick Meurice and co-PI Mary Hall Reno. New techniques designed for precision calculations of strong interaction physics were developed using the tensor renormalization group method. Large-scale Monte Carlo simulations with dynamical quarks were performed for candidate models for Higgs compositeness. Ab-initio lattice gauge theory calculations of semileptonic decays of B-mesons observed in collider experiments and relevant to test the validity of the standard model were performed with the Fermilab/MILC collaboration. The phenomenology of strong interaction physics was applied to new predictions for physics processes in accelerator physics experiments and to cosmic ray production and interactions. A research focus has been on heavy quark production and their decays to neutrinos. The heavy quark contributions to atmospheric neutrino and muon fluxes have been evaluated, as have the neutrino fluxes from accelerator beams incident on heavy targets. Results are applicable to current and future particle physics experiments and to astrophysical neutrino detectors such as the IceCube Neutrino Observatory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

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

Identification of Auger mechanisms responsible for low energy electron emission from graphene on copper using Auger-gamma coincidence spectroscopy

We have applied positron annihilation induced Auger-gamma coincidence spectroscopy to identify important mechanisms responsible for the emission of low energy electrons following the sudden creation of holes in bilayer graphene on copper substrate. The novel surface spectroscopic method measures the energy of the Doppler shifted annihilation gamma photon in coincidence with the Auger electron emitted following the relaxation of the hole created by the annihilation of a surface electron with a surface trapped positron. By extracting and theoretically modelling the annihilation gamma spectra coincident with low energy electrons, we associate majority of the intensity in the low energy (7 eV–25 eV) region of the Auger spectrum to electron emission following Auger decays of 2s holes in adsorbed oxygen and deep valence holes in graphene. We provide additional support to this conclusion by showing that most of the Auger electrons with energy less than ~ 25 eV are coincident with gamma photons with small Doppler shift (511 keV–512 keV) indicating that the primary electron whose annihilation resulted in low energy Auger electron emission had small momentum parallel to the gamma emission direction. On the other hand, majority of the Auger electrons measured in coincidence with the photons having maximum Doppler shift (515 keV–521 keV) were emitted following the decay of core holes (C 1s and O 1s). Furthermore, by selecting annihilation gamma photons coincident only with the O KVV (from surface adsorbed oxygen), C KVV (graphene), and Cu MVV (copper substrate) Auger electrons, we have also experimentally derived the energy spectrum of Doppler shifted gamma photons resulting from the annihilation of positrons with 1s electrons of O, 1s electrons of C and 3p electrons of Cu respectively. Model Doppler broadened gamma spectra produced using ab-initio calculations agree well with the experimentally derived line shapes. This demonstrates the ability of Auger-gamma coincidence method to experimentally resolve the Doppler broadened annihilation gamma spectrum from surfaces into its veiled electronic level constituents which has, heretofore, relied solely on theoretical analysis. Our results also demonstrate that the line shape of the Doppler broadened annihilation gamma peak reflect the chemical composition of the topmost atomic layer and can thus be used to characterize both external and inaccessible surfaces of porous materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗