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

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

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

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗
DOE OSTI2020

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

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

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

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

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

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

Isotope Shifts of Radium Monofluoride Molecules

Isotope shifts of 223–226,228 Ra 19 F were measured for different vibrational levels in the electronic transition A 2 II 1/2 ←X 2 Σ + . The observed isotope shifts demonstrate the particularly high sensitivity of radium monofluoride to nuclear size effects, offering a stringent test of models describing the electronic density within the radium nucleus. Ab initio quantum chemical calculations are in excellent agreement with experimental observations. These results highlight some of the unique opportunities that short-lived molecules could offer in nuclear structure and in fundamental symmetry studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Exploiting 20 Ne Isotopes for Precision Characterizations of Collectivity in Small Systems

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174⁢(8) stat ⁢(31) syst for NLEFT and 1.139⁢(6) stat ⁢(39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Atomic and Molecular Physics in Strong Fields

In this project, we have investigated a broad range of atomic, molecular, and optical (AMO) processes of current significance in strong-field and ultrafast science and technology. The main goals are to develop new theoretical formalisms and accurate computational methods for ab initio and in-depth non-perturbative study of time-resolved multiphoton quantum dynamics and spectroscopy, very-high-order harmonic generation, optimal control of ultrafast attosecond pulse generation, and resonance-enhanced below threshold harmonics in intense ultrashort laser fields. We have worked on the following projects and we have published 20 top 15 % SCI journal articles during the last 3.5 years (Jan. 1, 2016 to June 30, 2019), including one article in top raking Physical Review Letters, and 3 articles in top ranking (Nature) Scientific Reports, and one Physical Review A article was chosen as Editors’ suggestion paper, etc. (a) Development of new time-dependent methods for accurate and efficient solution of multi-dimensional time-dependent Schrödinger equation and time-dependent density functional theory. (b) Non-perturbative treatment of very-high order (1000th – 5000th order) harmonic generation and broadband supercontinuum spectra from XUV to soft x-ray regime. (c) Development of synchrosqueezing transform (SST) for precision time-frequency analysis of electron multi-rescattering dynamics. (d) Exploration of the dynamical origin of resonance-enhanced near- and below-threshold harmonics. (e) Probing the high-order harmonics generation in disordered semiconductor. (f) Multiphoton ionization of one-electron relativistic diatomic quasi-molecules in strong laser fields. and, (g) New development of non-adiabatic density functional theory. The theoretical and simulation projects we have pursued are closely related to the ongoing experimental researches in several leading laboratories in the US and Europe, etc. It is expected that the outcome of this project will provide timely and significant new information and insight for the advancement of strong-field atomic and molecular physics and ultrafast science and technology.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Silica Surface States and their Wetting Characteristics

The wetting characteristics of silica (SiO 2 ) surfaces can be described by molecular dynamics (MD) and ab initio simulations, including comparison of silica surfaces (talc (001), siloxanated quartz, tridymite (001) and quartz (001)), some of which have not been considered previously in the literature. Classic MD and ab initio simulation methods have been used to determine the contact angle, interfacial water structure, hydroxylation reaction and hydration energy, the results of which are compared with experimental results reported in the literature. Additionally, wetting of silica surfaces depends on surface polarity and extent of hydroxylation. The non-polar siloxane surfaces are characterized by a contact angle of about 80°, an MD ‘water exclusion zone’ of about 3 Å, a relaxed interfacial water orientation, inertness to hydroxylation and minimal hydration energy. The polar silica surfaces can be wetted by water and have a more ordered interfacial water structure. Silanol groups form at the polar silica surface during hydroxylation reactions, and the calculated hydration energy of -1·2 and -1·6 eV matches the experimental heat of immersion measurements reported in the literature, which correspond to hydrogen (H) bonding with interfacial water. In conclusion, fundamental understanding of silica surfaces is important for understanding flotation phenomena and fluid flow in silica nanopores.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Highly charged interface trap states in PbS 1– x govern electro-thermal transport

This work describes our discovery of the dominant role of highly charged interfaces on the electrothermal transport properties of PbS, along with a method to reduce the barrier potential for charge carriers by an order of magnitude. High temperature thermoelectrics such as PbS are inevitably exposed to elevated temperatures during postsynthesis treatment as well as operation. However, we observed that as the material was heated, large concentrations of sulfur vacancy (V $\ddot{S}$ ) sites were formed at temperatures as low as 266 °C. This loss of sulfur doped the PbS n -type and increased the carrier concentration, where these excess electrons were trapped and immobilized at interfacial defect sites in polycrystalline PbS with an abundance of grain boundaries. Sulfur deficient PbS 0.81 exhibited a large barrier potential for charge carriers of 0.352 eV, whereas annealing the material under a sulfur-rich environment prevented V $\ddot{S}$ formation and lowered the barrier by an order of magnitude to 0.046 eV. Through ab initio calculations, the formation of V $\ddot{S}$ was found to be more favorable on the surface compared to the bulk of the material with a 1.72 times lower formation energy barrier. Furthermore, these observations underline the importance of controlling interface-vacancy effects in the preparation of bulk materials comprised of nanoscale constituents.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Laser-induced electron diffraction of the ultrafast umbrella motion in ammonia

Visualizing molecular transformations in real-time requires a structural retrieval method with Ångström spatial and femtosecond temporal atomic resolution. Imaging of hydrogen-containing molecules additionally requires an imaging method sensitive to the atomic positions of hydrogen nuclei, with most methods possessing relatively low sensitivity to hydrogen scattering. Laser-induced electron diffraction (LIED) is a table-top technique that can image ultrafast structural changes of gas-phase polyatomic molecules with sub-Ångström and femtosecond spatiotemporal resolution together with relatively high sensitivity to hydrogen scattering. Here, we image the umbrella motion of an isolated ammonia molecule (NH3) following its strong-field ionization. Upon ionization of a neutral ammonia molecule, the ammonia cation (NH3+) undergoes an ultrafast geometrical transformation from a pyramidal (ΦHNH=107°) to planar (ΦHNH=120°) structure in approximately 8 femtoseconds. Using LIED, we retrieve a near-planar (ΦHNH=117 ± 5°) field-dressed NH3+ molecular structure 7.8−9.8 femtoseconds after ionization. Our measured field-dressed NH3+ structure is in excellent agreement with our calculated equilibrium field-dressed structure using quantum chemical ab initio calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simulated field-modulated x-ray absorption in titania

Here, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin–orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO 2 ) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t 2g peaks in the Ti L-edge was observed to red shift and the e g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

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

Improved Cycling Performance of Li-Excess Cation-Disordered Cathode Materials upon Fluorine Substitution

We report the recent discovery of Li-excess cation-disordered rock salt cathodes has greatly enlarged the design space of Li-ion cathode materials. Evidence of facile lattice fluorine substitution for oxygen has further provided an important strategy to enhance the cycling performance of this class of materials. Here, a group of Mn 3+ –Nb 5+ -based cation-disordered oxyfluorides, Li 1.2 Mn 3+ 0.6+0.5x Nb 5+ 0.2-0.5x O 2-x F x (x = 0, 0.05, 0.1, 0.15, 0.2) is investigated and it is found that fluorination improves capacity retention in a very significant way. Combining spectroscopic methods and ab initio calculations, it is demonstrated that the increased transition-metal redox (Mn 3+ /Mn 4+ ) capacity that can be accommodated upon fluorination reduces reliance on oxygen redox and leads to less oxygen loss, as evidenced by differential electrochemical mass spectroscopy measurements. Furthermore, it is found that fluorine substitution also decreases the Mn 3+ -induced Jahn–Teller distortion, leading to an orbital rearrangement that further increases the contribution of Mn-redox capacity to the overall capacity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Automatisierte Generierung von Mikrokinetiken für heterogen katalysierte Reaktionen unter Berücksichtigung korrelierter Unsicherheiten**

Abstract Die Studie stellt eine ab initio basierte Methode für die automatische Konstruktion von detaillierten Mikrokinetiken vor, welche korrelierte Unsicherheiten in allen energetischen Parametern und Schätzroutinen berücksichtigt. Hierfür wurden innerhalb des Unsicherheitsraums des BEEF‐vdW Funktionals 2000 Mikrokinetiken für die Oxidationsreaktionen repräsentativer Abgasemissionen von stöchiometrischen Verbrennungsmotoren auf Pt(111) generiert und durch Multiskalenmodellierung mit Experimenten verglichen. Das Ensemble der Simulationsergebnisse unterstreicht die Bedeutung der Berücksichtigung von Unsicherheiten. Innerhalb dieses Ensembles ist es möglich, eine Mikrokinetik zu identifizieren, welche mit den experimentellen Daten übereinstimmt. Diese Mikrokinetik lässt sich auf ein einziges Austausch‐Korrelations‐Funktional zurückführen und legt nahe, dass Pt(111) das aktive Zentrum für die Oxidation leichter Kohlenwasserstoffe sein könnte. Die Studie bietet einen universellen Rahmen für die automatisierte Konstruktion von Reaktionsmechanismen unter Berücksichtigung korrelierter Unsicherheiten und ermöglicht eine DFT‐basierte Optimierung der Mikrokinetik auch für andere heterogen katalysierte Systeme.

Kreitz, Bjarne↗
DOE OSTI2022

Coupled effect of Cr and Al on interactions between a prismatic interstitial dislocation loop and an edge dislocation line in Fe-Cr-Al alloy

Segregation of alloying elements to a prismatic dislocation loop under irradiation is an important phenomenon in understanding the role of loops in radiation effects in alloys. In this work, various segregation sequences of Cr and Al atoms to a $\langle$100$\rangle$ and a 1/2 $\langle$111$\rangle$ prismatic dislocation loop in a Fe-Cr-Al alloy are explored using molecular dynamics method and ab initio energy calculation. Here, the results reveal that Cr can readily segregate to the loops while Al atoms cannot by themselves, but the segregated Cr atoms are able to promote Al segregation at the loops, defined as a coupled segregation of Cr and Al. After interacting with Cr and Al, a previously mobile 1/2 $\langle$111$\rangle$ loop becomes sessile. Al segregation through the coupled segregation process can further affect the pinning behavior of a dislocation loop to an edge dislocation motion. By varying the Al fraction in the segregated solutes, a nonlinear dependence of the pinning effect on dislocation motion is determined. At lower fractions, the pinning strength increases with the increasing Al fraction until Al fraction reaches a critical value of about 19%, above which the pinning strength decreases with further increase of Al fraction. All these results suggest that, in addition to the effect of Al on the surface resistance to oxidation, Al effect on the dislocation loop behavior is important to optimize the composition of Fe-Cr-Al alloys for their applications in accident tolerant fuel concepts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic structure of reduced CeO 2 (111) surfaces interacting with hydrogen as revealed through electron energy loss spectroscopy in comparison with theoretical investigations

Based on both, ab-initio cluster calculations as well as periodic density functional based loss function calculations, we have assigned the origins of the valence electron excitation regime in electron energy loss spectra of well-ordered ceria films in (111) orientation at various states of reduction as well as after exposing the reduced films to hydrogen from the gas phase (Li et al, Angew. Chem. Int. Ed., 58 (2019) 14686–14693). The explicit calculation of intensity distributions using the dipole approximation allow us to draw conclusions about the nature of oxygen vacancies, which occur upon reduction and how those interact with hydrogen to form hydride species. Herein, we also reports a brief discussion of vibrational excitations of hydrogen loaded ceria and corroborates previous interpretations, based on inelastic neutron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗