Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Surface Phase Stability of Fe 2 O 3 (001) in Hydrogen Reducing Environments: A DFT and XPS Analysis

Here, this study combines density functional theory (DFT) and ab initio thermodynamics calculations with X-ray photoelectron spectroscopy (XPS) investigations to identify the reduction properties of the Fe 2 O 3 (001) surface with implications for corrosion resistance, hydrogen transport, and energy safety. Ab initio thermodynamics modeling predicts fully hydroxylated surface stability across a broad range of pressures (1 × 10 –23 to 1 × 10 5 mbar) and temperatures below 700 K, consistent with previous experimental studies. Above 800 K, exposures to 1 × 10 –4 mbar H 2 , 1 × 10 –4 mbar O 2 , or 1 × 10 –4 mbar H 2 + 1 × 10 –4 mbar O 2 each yield unique XPS signals indicating a loss of −OH coverage, aligning with DFT predictions. Insight into the mechanism of reduction as a function of H 2 exposure is provided, as well as conditions that promote further reduction toward Fe 3 O 4 . Theoretical and experimental investigations indicate the ability to maintain the Fe 2 O 3 protective layer of iron oxides that have been exposed to H 2 environments by including trace amounts of aqueous O2.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Hot-Hole versus Hot-Electron Transport at Cu/GaN Heterojunction Interfaces

Among all plasmonic metals, copper (Cu) has the greatest potential for realizing optoelectronic and photochemical hot-carrier devices, thanks to its CMOS compatibility and outstanding catalytic properties. Yet, relative to gold (Au) or silver (Ag), Cu has rarely been studied and the fundamental properties of its photoexcited hot carriers are not well understood. Here, we demonstrate that Cu nanoantennas on p-type gallium nitride (p-GaN) enable hot-hole-driven photodetection across the visible spectrum. Importantly, we combine experimental measurements of the internal quantum efficiency (IQE) with ab initio theoretical modeling to clarify the competing roles of hot-carrier energy and mean-free path on the performance of hot-hole devices above and below the interband threshold of the metal. We also examine Cu-based plasmonic photodetectors on corresponding n-type GaN substrates that operate via the collection of hot electrons. By comparing hot hole and hot electron photodetectors that employ the same metal/semiconductor interface (Cu/GaN), we further elucidate the relative advantages and limitations of these complementary plasmonic systems. In particular, we find that harnessing hot holes with p-type semiconductors is a promising strategy for plasmon-driven photodetection across the visible and ultraviolet regimes. Given the technological relevance of Cu and the fundamental insights provided by our combined experimental and theoretical approach, we anticipate that our studies will have a broad impact on the design of hot-carrier optoelectronic devices and plasmon-driven photocatalytic systems.

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

Computational and Experimental Characterization of the Ligand Environment of a Ni-Oxo Catalyst Supported in the Metal–Organic Framework NU-1000

Heterogeneous catalysts exhibit significant changes in composition due to the influence of operating conditions, and these compositional changes can have dramatic effects on catalytic performance. For traditional bulk metal heterogeneous catalysts, relationships between composition and catalytic operating conditions are well documented. However, the influence of operating conditions on the compositions of single-site heterogeneous catalysts remains largely unresolved. To address this, we report a combined computational and experimental characterization of a Ni oxo catalyst under catalytic hydrogenation conditions. Specifically, pair distribution function (PDF) analysis is combined with ab initio thermodynamic modeling to investigate ligand environments present on a Ni oxo cluster supported in the metal–organic framework NU-1000. Comparisons of the experimentally observed and simulated Ni–O coordination numbers and Ni–O, Ni···Ni, and Ni···Zr distances provide insight into the Ni ligand environment under H 2 (g). These comparisons suggest significant OH and H 2 O content and, further, that different Ni ions within the cluster and/or NU-1000 structure may comprise subtly different numbers of these ligands. Further, the observation of significant H 2 O content under H 2 (g) suggests that the NU-1000 support supplies H 2 O to the cluster. Examples of ligand environments that could lead to the observed PDFs are provided. Furthermore, the combination of simulations and experiments provides new insights into the ligand environment for Ni-NU-1000 catalysts that will be useful for understanding the ligand environments of other single-site Ni catalysts as well.

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

Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry

Due to intense interest in the potential applications of quantum computing, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics. Through numerical studies of quantum state preparation and empirical complexity analysis (including the error scaling) of classical heuristics, in both ab initio and model Hamiltonian settings, we conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for ground-state quantum chemistry through polynomial speedups, it may be prudent to assume exponential speedups are not generically available for this problem.

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

Machine learning for deep elastic strain engineering of semiconductor electronic band structure and effective mass

Abstract The controlled introduction of elastic strains is an appealing strategy for modulating the physical properties of semiconductor materials. With the recent discovery of large elastic deformation in nanoscale specimens as diverse as silicon and diamond, employing this strategy to improve device performance necessitates first-principles computations of the fundamental electronic band structure and target figures-of-merit, through the design of an optimal straining pathway. Such simulations, however, call for approaches that combine deep learning algorithms and physics of deformation with band structure calculations to custom-design electronic and optical properties. Motivated by this challenge, we present here details of a machine learning framework involving convolutional neural networks to represent the topology and curvature of band structures in k -space. These calculations enable us to identify ways in which the physical properties can be altered through “deep” elastic strain engineering up to a large fraction of the ideal strain. Algorithms capable of active learning and informed by the underlying physics were presented here for predicting the bandgap and the band structure. By training a surrogate model with ab initio computational data, our method can identify the most efficient strain energy pathway to realize physical property changes. The power of this method is further demonstrated with results from the prediction of strain states that influence the effective electron mass. We illustrate the applications of the method with specific results for diamonds, although the general deep learning technique presented here is potentially useful for optimizing the physical properties of a wide variety of semiconductor materials.

Tsymbalov, Evgenii↗
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 OSTI2025

Effective field theories for collective excitations of atomic nuclei

Collective modes emerge as the relevant degrees of freedom that govern low-energy excitations of atomic nuclei. These modes—rotations, pairing rotations, and vibrations—are separated in energy from non-collective excitations, making it possible to describe them in the framework of effective field theory. Rotations and pairing rotations are the remnants of Nambu–Goldstone modes from the emergent breaking of rotational symmetry and phase symmetries in finite deformed and finite superfluid nuclei, respectively. The symmetry breaking severely constrains the structure of low-energy Lagrangians and thereby clarifies what is essential and simplifies the description. The approach via effective field theories exposes the essence of nuclear collective excitations and is defined with a breakdown scale in mind. This permits one to make systematic improvements and to estimate and quantify uncertainties. Effective field theories of collective excitations have been used to compute spectra, transition rates, and other matrix elements of interest. In particular, predictions of the nuclear matrix element for neutrinoless double beta decay then come with quantified uncertainties. This review summarizes these results and also compares the approach via effective field theories to well-known models and ab initio computations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Prediction of double-Weyl points in the iron-based superconductor CaKFe 4 As 4

Employing a combination of symmetry analysis, low-energy modeling, and ab initio simulations, we predict the presence of magnetic-field-induced Weyl points close to the Fermi level in CaKFe 4 As 4 . Depending on the relative strengths of the magnetic field and of the spin-orbit coupling, the Weyl fermions can carry a topological charge of ±1 or ±2, making CaKFe 4 As 4 a rare realization of a double-Weyl semimetal. Here, we further predict experimental manifestations of these Weyl points, both in bulk properties, such as the anomalous Hall effect, and in surface properties, such as the emergence of prominent Fermi arcs. Because CaKFe 4 As 4 displays unconventional fully gapped superconductivity below 30 K, our findings open a route to investigate the interplay between superconductivity and Weyl fermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unraveling the magnetic ground state in the alkali-metal lanthanide oxide Na 2 ⁢Pr⁡O 3

Here, a comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na 2 ⁢PrO 3 . μSR reveals a Néel antiferromagnetic order below T N ~ 4.9K, with a small static magnetic moment m static ≤ 0.22 μ B /Pr collinearly aligned along the c axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr 4+ ground-state wave function deviates significantly from the Γ 7 limit that is relevant to the Kitaev model. Single- and two-magnon excitations are observed in the ordered state below T N =4.6K and are well described by nonlinear spin wave theory from the Néel state using a magnetic Hamiltonian with Heisenberg exchange J=1 meV and symmetric anisotropic exchange Γ/J=0.1, corresponding to an XY model. Intense two magnon excitations are accounted for by g-factor anisotropy g z /g ± = 1.29. A fluctuating moment δm 2 = 0.57 ⁢(22) ⁢μ$^2_B$/Pr extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local J = 1/2 moment with average g factor g avg ≈ 1.1. Together, the results demonstrate that the small moment in Na 2 ⁢PrO 3 arises from crystal field and covalency effects and the material does not exhibit significant quantum fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coulomb excitation of the | T z | = 1 2 , A = 23 mirror pair

Background: Electric-quadrupole ( E 2 ) strengths relate to the underlying quadrupole deformation of a nucleus and present a challenge for many nuclear theories. Additionally, mirror nuclei in the vicinity of the line of N = Z represent a convenient laboratory for testing deficiencies in such models, making use of the isospin symmetry of the systems. Purpose: Uncertainties associated with literature E 2 strengths in Mg 23 are some of the largest in T z = | 1 2 | nuclei in the s d shell. The purpose of the present paper is to improve the precision with which these values are known, to enable better comparison with theoretical models. Methods: Coulomb-excitation measurements of Mg 23 and Na 23 were performed at the TRIUMF-ISAC facility using the TIGRESS spectrometer. They were used to determine the E 2 matrix elements of mixed E 2 / M 1 transitions. Results: Reduced E 2 transition strengths, B ( E 2 ) , were extracted for Mg 23 and Na 23 . Their precision was improved by factors of approximately 6 for both isotopes, while agreeing within uncertainties with previous measurements. Conclusions: A comparison was made with both shell-model and ab initio valence-space in-medium similarity renormalization group calculations. Valence-space in-medium similarity renormalization group calculations were found to underpredict the absolute E 2 strength, in agreement with previous studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment

36 MATERIALS SCIENCE↗
DOE OSTI2022

Controlling Structural, Electronic, and Energy Flow Dynamics of Catalytic Processes through Tailored Nanostructures

MoS 2 (molybdenum disulfide) is a highly-versatile catalyst material for support of numerous reactions from hydrodesulfurization and denitrogenation to the focus of this renewal proposal: hydrogenation of CO/CO 2 towards (higher) alcohols. At the same time, MoS 2 is a non-toxic, environmentally-benign and rather inert material- which under ambient conditions for some time has served as a lubricant and, more recently, as next-generation electronic material. The apparent contrast between inertness and stability in ambient, and catalytic activity under reactive conditions is puzzling and calls for a synergistic theoretical and experimental investigation with the long-term objective of enabling the rational design of MoS 2 -based catalysts for alcohol-formation reactions by providing a microscopic understanding of the environmental factors that determines site activity and selectivity. Our research project seeks answers to the questions (a) what conformation does MoS 2 adopt under reaction conditions (as opposed to that under ultrahigh vacuum and low temperatures)?; (b) what reaction pathways exist on such a material?; (c) how can the local environment of the active sites be manipulated so as to make MoS 2 an efficient catalyst for production of higher alcohol from syngas? In particular, research strategies will explore how the basal plane composed of sulfur atoms can be activated so as to exhibit a reactivity of its own, by addressing three research targets and building on extensive preliminary and enabling work: (1) vacancies and vacancy aggregates on the basal plane; (2) non-local catalyst transformation through alkali doping, hydrogenation and phase transition; (3) fabrication of a metal-nanoparticle-activated MoS 2 system, in which particle anchoring, reactive sites and pathways as well as selectivity are controlled by design, as an example of predictive development of a catalyst material,. All strategies are directed to improve the efficacy of the key reactive sites and selectivity of chemical pathways by design, to replace the inefficient methodology of trial and error in catalyst development. This research project represents a synergistic combination of computational guidance, foundational surface-science-based experiments and validation under reactive conditions that aims at transformative new insights into the working of MoS 2 -based hydrogenation-catalysts, a topic squarely at the center of the interest of DOE BES. Guided and led by Talat Rahman, a computational physicist, this project will apply density functional theory to understand structure, reaction pathways and chemical potentials associated with MoS 2 -based CO/CO 2 hydrogenation, augmented by kinetic Monte Carlo methods for reaction rates and prefactors as well as ab-initio molecular dynamics for evaluation of thermal stability. Complementary experimental input and validation will originate from co-PI Ludwig Bartels, a physical chemist and materials scientist, whose group focuses on local imaging and preparation of MoS 2 materials, from co-PI Peter Dowben, an experimental physicist, whose group is expert in the spectroscopy of occupied and unoccupied electronic states, and from senior collaborator Michael White of Brookhaven National Laboratory, whose group generates high-resolution electronic and activity information on size-selected well-defined metal chalcogenide clusters. This collaborative effort will enable a comprehensive understanding of the correlation of structural integrity and catalytic activity of MoS 2 in forms ranging from extended films to individual particles with known geometries and binding sites. Alcohol formation from syngas is a rapidly emerging application that has great potential through facile, economic, and decentralized biomass gasification. CO 2 activation is one of the most pressing concerns of our time: increasing CO 2 levels in the atmosphere change the climate and expose the globe to environmental transformations with the potential for enormous economic and societal impact. We will investigate CO/CO 2 hydrogenation via an interdisciplinary research collaboration with established synergy – one that, in accordance with the mission of the DOE, involves accredited Hispanic-Serving Institutions and that, through student exchange with international collaborators and National Labs, will directly benefit a broad spectrum of communities and generate human resources in sciences essential to their future.

2D materials↗
DOE OSTI2022

Mutations of Omicron Variant at the Interface of the Receptor Domain Motif and Human Angiotensin-Converting Enzyme-2

The most recent Omicron variant of SARS-CoV-2 has caused global concern and anxiety. The only thing certain about this strain, with a large number of mutations in the spike protein, is that it spreads quickly, seems to evade immune defense, and mitigates the benefits of existing vaccines. Based on the ultra-large-scale ab initio computational modeling of the receptor binding motif (RBM) and the human angiotensin-converting enzyme-2 (ACE2) interface, we provide the details of the effect of Omicron mutations at the fundamental atomic scale level. In-depth analysis anchored in the novel concept of amino acid-amino acid bond pair units (AABPU) indicates that mutations in the Omicron variant are connected with (i) significant changes in the shape and structure of AABPU components, together with (ii) significant increase in the positive partial charge, which facilitates the interaction with ACE2. We have identified changes in bonding due to mutations in the RBM. The calculated bond order, based on AABPU, reveals that the Omicron mutations increase the binding strength of RBM to ACE2. Our findings correlate with and are instrumental to explain the current observations and can contribute to the prediction of next potential new variant of concern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1977

The molecular composition of dense interstellar clouds

Presented in this paper is an ab initio chemical model for dense interstellar clouds that incorporates 598 grain surface reactions, with small grains providing the reaction area. Gas-phase molecules are depleted through collisions with grains. The abundances of 372 chemical species are calculated as a function of time and are found to be of sufficient magnitude to explain most observations. Peak abundances are achieved on time scales of the order of 100,000 to 1 million years, depending on cloud density and kinetic temperature. The reaction rates for ion-molecule chemistry are approximately the same, indicating that surface and gas-phase chemistry may be coupled in certain regions. The composition of grain mantles is shown to be a function of grain radius. In certain grain-size ranges, large molecules containing two or more heavy atoms are more predominant than lighter 'ices' - H2O, NH3, and CH4. It is possible that absorption due to these large molecules in the mantle may contribute to the observed 3-micron band in astronomical spectra.

Allen, M.↗
NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗