Search NASA⌕ Search

SEARCH · Search NASA

Results for “First Principles”

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

Hydrogen defects in LaBi 2 O 4 X (X = Cl, Br, and I) Sillén oxyhalide phases and their impacts on ionic transport

Sillén oxyhalides have recently emerged as promising materials for both photocatalytic and ionic transport applications, yet the role of likely-ubiquitous hydrogen-related defects in these layered compounds remains largely unexplored. Here, we employ first-principles defect calculations to investigate incorporation energetics for hydrogen- and oxygen-related defects, as well as their migration barriers in LaBi 2 O 4 X (X = Cl, Br, I) phases. We find that hydrogen interstitials, particularly protonic species (H i + ), are readily accommodated within the open Bi–O layers. Protons compete with oxygen vacancy donors (V O 2+ ) and charge-compensate with oxygen interstitial acceptors (O i 2− ). By linking hydrogen defect formation to water- and oxygen-related redox equilibria, we reveal that V O 2+ facilitates H i + incorporation, while O i 2− promotes interstitial hydroxide formation, establishing a direct connection between proton and oxide-ion transport. Calculated migration barriers indicate that ionic diffusion is confined to Bi–O layers with low barriers of 0.20–0.25 eV for H i + and 0.14–0.25 eV for V O 2+ , suggesting that the materials contain intrinsic pathways for mixed ionic conduction. These results provide a microscopic picture of hydrogen behavior in Sillén oxyhalides and point to design strategies for integrating protonic and oxide-ion transport in layered oxyhalide electrolytes. Band-edge alignment analysis shows that LaBi 2 O 4 I provides the optimal combination of hydrogen solubility, oxygen defect stability, and mixed ionic conductivity, highlighting its potential for low-temperature electrochemical and energy-conversion applications. Overall, this work establishes the defect-driven origin of hydrogen transport in Sillén oxyhalides and expands their applicability beyond photocatalysis to mixed ionic conduction and hydrogen electrochemistry.

Energy - Conversion↗
DOE OSTI2026

Non-adiabatic quantum interference and complex formation in ultracold collisions of Rb with KRb

Ultracold elastic collisions of 87 Rb with 40 K 87 Rb in its ground vibrational and rotational state are investigated using a first principles based theoretical methodology. Full-dimensional ab initio computed potential energy surfaces are reported that include the two lowest-lying electronic states, their conical intersection, non-adiabatic couplings and an accurate long-range behavior. A numerically exact time-independent quantum dynamics method in hyperspherical coordinates is used to compute the elastic scattering cross sections, rate coefficients and collision lifetime spectrum. The quantum scattering calculations include all degrees of freedom and treat both electronic states and their non-adiabatic couplings using a two-state diabatic representation. The theoretically computed elastic rate coefficient is in good agreement with the recently reported experimental value. Significant non-adiabatic quantum interference effects are shown to originate from the unique properties of ultracold collisions and the geometric phase associated with the conical intersection. A high-resolution collision energy grid is used to investigate the origin of the experimentally reported long-lived 3-body collision complexes.

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

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗
DOE OSTI2026

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗
DOE OSTI2026

Electron-doping-induced destabilization of the dimerized insulating state in monolayer IrTe 2

The ability to tune the electronic phases of two-dimensional (2D) materials through external perturbations provides a powerful route to engineer functional nanoscale systems. In particular, the ground state of monolayer (ML) 1T-IrTe 2 is highly sensitive to the interplay between local chemical bonding and global electronic topology, leading to a unique 2 × 1 dimerized insulating phase. Here, we present an angle-resolved photoemission study on the evolution of the electronic structure in ML IrTe 2 induced by in situ Rb adsorption. We find that Rb adsorption suppresses the 2 × 1 dimerized phase in ML IrTe 2 , inducing a clear insulator-to-metal transition. This transition is characterized by a reconstruction of the band topology toward a bilayer-like metallic configuration. Combined with first-principles calculations, our results demonstrate that the collapse of the insulating state originates from Ir valence change and suppression of Fermi surface nesting-driven instabilities. Our findings establish ML IrTe 2 as a model system for investigating instability-driven phase control in 2D materials and highlight its potential for tunable electronic applications.

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

Computational discovery of a tetragonal MBene phase: diamond-shaped boron lattices, viable MAB precursors, and selective CO 2 reduction

MBenes are a class of two-dimensional transition-metal borides derived from layered MAB phases. Here, we use first-principles calculations to predict a new tetragonal MB (1 : 1) MBene phase featuring a diamond-shaped B–B lattice and systematically assess its phase stability, synthetic accessibility, and catalytic performance towards the CO 2 reduction reaction (CO 2 RR). Four tetragonal configurations are benchmarked against known hexagonal and orthorhombic phases, and a multi-tier screening identifies six robust members (CrB, FeB, MoB, WB, IrB, and PtB). Exploration of the I4/mmm MAB space reveals 11 viable precursors, including the experimentally synthesized Ir 2 ZnB 2 , which validates our screening approach. CO 2 RR free-energy diagrams show that CrB favors CH 3 OH while FeB, MoB, WB, IrB, and PtB preferentially yield HCOOH. The corresponding limiting potentials are competitive with, and in some cases superior to, those reported for state-of-the-art orthorhombic and hexagonal MBenes. These findings expand the MBene landscape by establishing a tetragonal phase with structural, synthetic, and catalytic promise.

Lu, Linguo [University of Puerto Rico, San Juan, P↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2024

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

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

NiGa 2 O 4 interfacial layers in NiO/Ga 2 O 3 heterojunction diodes at high temperature

NiO/Ga 2 O 3 heterojunction diodes have attracted attention for high-power applications, but their high temperature performance and reliability remain underexplored. Here, we report the time evolution of the electrical properties in the widely studied p-NiO/n-Ga 2 O 3 heterojunction diodes and formation of NiGa 2 O 4 interfacial layers at high temperatures. Results of our thermal cycling experiment show an initial leakage current increase which stabilizes after sustained thermal load, due to reactions at the NiO–Ga 2 O 3 interface. High-resolution TEM microstructure analysis of the devices after thermal cycling indicates that the NiO–Ga 2 O 3 interface forms a ternary compound at high temperatures, and thermodynamic calculations suggest the formation of the spinel NiGa 2 O 4 layer between NiO and Ga 2 O 3 . First-principles defect calculations find that NiGa 2 O 4 shows low p-type intrinsic doping and hence can serve to limit electric field crowding at the interface. Vertical NiO/ Ga 2 O 3 diodes with intentionally grown 5 nm thin spinel-type NiGa 2 O 4 interfacial layers show an excellent device ON/OFF ratio of >10 10 (± 3 V), V ON of ~1.9 V, and increased breakdown voltage of ~1.2 kV for an initial unoptimized 300 lm diameter device. These p–n heterojunction diodes are promising for high-voltage, high temperature applications.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Localized phase contrast imaging at the Wendelstein 7-X stellarator

In its basic form, phase contrast imaging (PCI) provides line-integrated measurements of electron density fluctuations in plasmas. As turbulent fluctuations in magnetically confined plasmas have wave vectors almost perpendicular to the background magnetic field, the signals scattered by fluctuations from different parts of the PCI line-of-sight (LoS) are spatially separated in focal planes of the plasma. This allows localized PCI measurements by placing a mask in such a plane, to only permit signals from specific parts of the LoS to reach the PCI detectors. The present paper describes modeling and design of localization masks for the PCI system at the Wendelstein 7-X (W7-X) stellarator as well as the first results obtained using the masks in the recent long-pulse W7-X experimental campaign. During this project, we have extended the theory describing the mask response within the Fraunhofer diffraction model. As a novel development, we show from first principles that the mask response is determined by the fraction of power of the scattered beam spots that passes the mask. These insights have been used to select the W7-X mask design, consisting of a circular cutout, allowing the unscattered beam spot to pass the mask, with wedges covering a fixed angular range outside the central cutout. In the recent W7-X experimental campaign, the masks have verified the location of the main turbulence features observed by the PCI system and provided new information about the location of short-wavelength magnetohydrodynamic modes.

Physics↗
DOE OSTI2024

Review of the second charged-particle transport coefficient code comparison workshop

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24–27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron–ion properties such as the electron–ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient exact exchange using Wannier functions and other related developments in planewave-pseudopotential implementation of RT-TDDFT

The plane-wave pseudopotential (PW-PP) formalism is widely used for the first-principles electronic structure calculation of extended periodic systems. The PW-PP approach has also been adapted for real-time time-dependent density functional theory (RT-TDDFT) to investigate time-dependent electronic dynamical phenomena. In this work, we detail recent advances in the PW-PP formalism for RT-TDDFT, particularly how maximally localized Wannier functions (MLWFs) are used to accelerate simulations using the exact exchange. We also discuss several related developments, including an anti-Hermitian correction for the time-dependent MLWFs (TD-MLWFs) when a time-dependent electric field is applied, the refinement procedure for TD-MLWFs, comparison of the velocity and length gauge approaches for applying an electric field, and elimination of long-range electrostatic interaction, as well as usage of a complex absorbing potential for modeling isolated systems when using the PW-PP formalism.

Chemistry↗
DOE OSTI2024

Selective incorporation of antimony into gallium nitride

Dilute concentrations of antimony (Sb) incorporation into GaN induce strong bandgap bowing and tunable room-temperature photoluminescence from the UV to the green spectral regions. However, the atomistic details of the incorporation of Sb into the GaN host remain unclear. In this work, we use first-principles calculations to understand the thermodynamics of Sb substitution into GaN and its effect on the optical and Raman spectra. Although it is empirically considered that Sb is preferentially incorporated as an anion (Sb3−) into the N sublattice, we demonstrate that Sb can also be incorporated as a cation (Sb3+, Sb5+) into the metal sublattice. Our thermodynamic analysis demonstrates that SbN0, SbGa2+, and SbGa0 can co-exist under Ga-rich conditions in n-type samples. We further confirm the dual incorporation of Sb by calculating the vibrational frequencies of different anionic and cationic substitutions to explain the origins of experimentally observed additional Raman peaks of Sb-doped GaN. Moreover, the calculated band structures of different Sb substitutions into GaN explain the experimental photoluminescence and optical absorption spectra. Overall, our analysis suggests that the coexistence of Sb3−, Sb3+, and Sb5+ substitutions into GaN explains the totality of experimental measurements. Our results demonstrate that the selective incorporation of Sb into GaN (and potentially other group-V elements such as As, P, or Bi) by tuning the growth conditions can drastically modify the electronic properties, for applications in visible light emitters and photocatalysis.

Physics↗
DOE OSTI2024

A high-throughput and data-driven computational framework for novel quantum materials

Two-dimensional layered materials, such as transition metal dichalcogenides (TMDs), possess an intrinsic van der Waals gap at the layer interface, allowing for remarkable tunability of the optoelectronic features via external intercalation of foreign guests such as atoms, ions, or molecules. Herein, we introduce a high-throughput, data-driven computational framework for the design of novel quantum materials derived from intercalating planar conjugated organic molecules into bilayer transition metal dichalcogenides and dioxides. By combining first-principles methods, material informatics, and machine learning, we characterize the energetic and mechanical stability of this new class of materials and identify the fifty (50) most stable hybrid materials from a vast configurational space comprising ∼105 materials, employing intercalation energy as the screening criterion.

Kastuar, Srihari M. (ORCID:0000000279001561)↗
DOE OSTI2024

High-pressure characterization of Ag 3 AuTe 2 : Implications for strain-induced band tuning

Recent band structure calculations have suggested the potential for band tuning in the chiral semiconductor Ag 3 AuTe 2 to zero upon application of negative strain. In this study, we report on the synthesis of polycrystalline Ag 3 AuTe 2 and investigate its transport and optical properties and mechanical compressibility. Transport measurements reveal the semiconducting behavior of Ag 3 AuTe 2 with high resistivity and an activation energy E a of 0.2 eV. The optical bandgap determined by diffuse reflectance measurements is about three times wider than the experimental E a ⁠. Despite the difference, both experimental gaps fall within the range of predicted bandgaps by our first-principles density functional theory (DFT) calculations employing the Perdew–Burke–Ernzerhof and modified Becke–Johnson methods. Furthermore, our DFT simulations predict a progressive narrowing of the bandgap under compressive strain, with a full closure expected at a strain of –4% relative to the lattice parameter. To evaluate the feasibility of gap tunability at such substantial strain, the high-pressure behavior of Ag 3 AuTe 2 was investigated by in situ high-pressure x-ray diffraction up to 47 GPa. Mechanical compression beyond 4% resulted in a pressure-induced structural transformation, indicating the possibility of substantial gap modulation under extreme compression conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An energy-modified quantum defect method for the analysis of Rydberg spectra: Application to 2-butyne

The high resolution Rydberg absorption spectrum of 2-butyne C 4 H 6 recorded previously at the SOLEIL synchrotron facility has been interpreted using multichannel quantum defect theory (MQDT). The calculations are based on the continuum scattering calculations of Xu et al., J. Chem. Phys. 136, 154303 (2012) and of Jacovella et al., J. Phys. Chem. A 119, 12339 (2015) pertaining to the dipole-allowed excited state symmetries in absorption from the ground state. In contrast to the traditional approach of calculating low-lying electronic states first and then attempting to extend the calculations to ever higher energy, here the analysis proceeds through the extension of these detailed calculations of the electronic continuum scattering down into the discrete region of the spectrum. The continuum reaction matrices and dipole transition moments are adapted to the discrete Rydberg region via the use of an energy-modified formulation of MQDT theory and associated energy dependences of the quantum defects. The analysis reproduces more than 40 Rydberg states from n ≈ 10 down to the 3 d and 4 s levels with an rms error of better than 20 cm −1 . Further, these belong to five Rydberg series with three different molecular symmetries. While the approach predicts many additional series, most of these are calculated and observed to carry only little oscillator strength. The analysis shows that the Rydberg spectrum is dominated by the excitation of an e″ symmetry electron of fδ and gπ type, in line with what previous studies of the above-threshold shape resonance of 2-butyne have shown. The present study is intended to serve as an example showing how first principles continuum calculations may be useful for the interpretation of highly bound discrete states in a range that poses problems for the standard ab initio techniques. The quantitative treatment of the dipole absorption cross sections is deferred to a future paper.

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

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗