Search NASA⌕ Search

SEARCH · Search NASA

Results for “First principles calculations”

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

Alloying multiple halide perovskites on the same sublattice in search of stability and target band gaps

Single-component halide perovskites (HPs) rarely satisfy all the necessary criteria for optoelectronic applications, such as achieving an optimal band gap while maintaining high chemical and structural stability. Alloying halide perovskites has emerged as a promising strategy, not only to enhance stability but also to fine-tune their electronic and optical properties. In this work, we explore multiple degrees of freedom in alloy design, considering different substitution sublattices sites (A, B, or X in ABX3 perovskites), various chemical species (isovalent and hetero-valent elements), and multi-component compositions on a given sublattice. Using first-principles calculations based on density functional theory (DFT), we investigate how compositional variations influence the electronic (band gap) and structural properties (mixing enthalpy) of HP alloys. Our approach employs the polymorphous cell model, allowing full local relaxation which breaks local symmetry while preserving global cubic symmetry—an essential framework for accurately modeling HPs. Our results reveal that X-site mixing (halogen substitution) primarily affects the valence band maximum, allowing target band gap engineering. Additionally, variations in halogen radii introduce internal strain through octahedral distortions, influencing the mixing enthalpy. A-site substitution, while not directly contributing to the band edge states, modifies structural stability via volume effects, indirectly impacting the band gap. B-site alloying plays a dominant role in band gap modulation, leading to either positive or negative band gap bowing. Specifically, isovalent B-site mixing (Sn–Pb) induces strong positive bowing, where the alloy band gap is smaller than the average gap of parent compounds, whereas hetero-valent mixing (Cd–Pb) results in pronounced negative bowing. As an aside, we investigate the competition between the excess energy of disordered alloys vs. that of long-range ordered double perovskites of the same compositions, seeking examples of ordered phases emerging from disordered alloys. Furthermore, our findings provide fundamental insights into the electronic and structural behavior of HP alloys, offering valuable design principles for the development of stable and efficient materials for next-generation photovoltaic and optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2026

Harnessing ketones as hydrogen acceptors for atom-efficient upgrading of oxygenates to fuels over H/ZSM-5

Heterogeneous mixtures of bio-derived oxygenates are promising feedstocks for synthetic aviation fuel (SAF), but conversion strategies for one common component—short-chain (C 5–7 ) internal ketones—are lacking. Previous work has shown that cyclization of ketones over H/ZSM-5 is limited by its high productivity of light paraffins. We study 4-heptanone upgrading over H/ZSM-5 and show that aromatics and olefins can be formed at high carbon yield when operating at up to 90% conversion. The yield of desirable products is not impacted by the introduction of a recycle stream of the unconverted 4-heptanone and other products with similar boiling points. We hypothesize, based on first-principles calculations, that the higher olefin yield is driven by the ease of hydrogen transfer to unreacted ketones as opposed to hydrogenating olefin products. We demonstrate how this ease of hydrogen transfer to ketones can be leveraged to enhance olefin selectivity in the conversion of methanol to olefins as well by co-feeding ketones. Olefinic products of the cyclization reaction are then oligomerized to a SAF blendstock to demonstrate an end-to-end ketone-to-SAF process facilitated by upgrading ketones over H/ZSM-5 at partial conversion with a recycle stream. The results of this work demonstrate a strategy to improve the carbon yield from bio-derived acids to SAF to over 75%, representing a relative improvement of more than 50% compared to previously reported data.

08 HYDROGEN↗
DOE OSTI2025

Mechanistic insights into superionic thioarsenate argyrodite solid electrolytes via machine learning interatomic potentials

The lithium argyrodite sulfide solid electrolyte Li 6 PS 5 Cl has attracted considerable interest for all-solid-state batteries owing to its high ionic conductivity, which can be further enhanced through ionic substitution. Although a variety of substitutions have been investigated, thioarsenate argyrodites remain comparatively underexplored. Here, we systematically investigate the phase stability and Li-ion conduction mechanisms in superionic Br-incorporated thioarsenate argyrodites using first-principles calculations and molecular dynamics simulations based on machine learning interatomic potentials (MLIPs). Systematic variation of S/Br site inversion reveals that an optimal degree of anion disorder significantly enhances inter-cage connectivity and facilitates long-range Li-ion diffusion. Configurational entropy serves as an effective quantitative descriptor of anion disorder, exhibiting a strong correlation with ionic conductivity. While greater anion disorder induced by site inversion and higher Br content enhances ionic conductivity up to 50 mS cm −1 , it simultaneously reduces structural stability. This trade-off results in an optimal window in which a moderate level of disorder yields conductivities exceeding 20 mS cm −1 while maintaining synthetic feasibility. In conclusion, this work highlights the reliability and efficiency of MLIPs for elucidating ion-transport mechanisms and accelerating the design of novel superionic argyrodites.

Jang, Myeongcho [Korea Institute of Science and Te↗
DOE OSTI2025

Charting the chemical space of Zintl phases with graph neural networks and bonding insights

A large number of Zintl phases have been discovered by solid-state chemists driven by empirical knowledge, chemical intuition and in some cases, through serendipitous accidents. These discoveries have only scratched the surface, given the vast compositional and structural diversity that Zintl phases can accommodate. The large chemical space of Zintl phases, as well as intermetallic compounds in general, remain under-explored. Here, we use graph neural networks and the upper bound energy minimization approach to efficiently scan a large chemical space of >90 000 hypothetical Zintl phases and accurately discover 1810 new thermodynamically stable phases with 90% precision, as validated with first-principles calculations. We show that our approach is more than 2× more accurate in predicting DFT stability than M3GNet (40% precision) on the same dataset. Using a random forest model and SHAP analysis, we demonstrate the critical role of ionic bonding in the thermodynamic stability of Zintl phases. Our results not only expand the known chemical landscape of Zintl phases but also highlight the efficacy of machine learning frameworks combined with domain knowledge in uncovering chemically meaningful insights across complex intermetallics.

36 MATERIALS SCIENCE↗
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 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

NEXAFS spectroscopy of alkylated benzothienobenzothiophene thin films at the carbon and sulfur K-edges

Alkylated benzothienobenzothiophenes are an important class of organic semiconductors that exhibit high performance in solution-processed organic field-effect transistors. In this work, we study the near-edge x-ray absorption fine-structure (NEXAFS) spectra of 2,7-didecyl[1]benzothieno[3,2-b][1]benzothiophene (C 10 -BTBT) at both the carbon and sulfur K-edges. Angle-resolved experiments of thin films are performed to characterize the dichroism associated with molecular orientation. First-principles calculations using the density functional theory-based many-body x-ray absorption spectroscopy (MBXAS) method are also performed to correlate the peaks observed and their dichroism with transitions to specific antibonding molecular orbitals. Interestingly, the dichroism of the dominant, lowest energy peak is opposite at the carbon and sulfur K-edges. While the low-energy peak at the carbon K-edge is assigned to carbon 1s → π* transitions with transition dipole moment (TDM) perpendicular to the planar BTBT core, the dominant low energy peak at the sulfur K-edge is assigned to sulfur 1s → σ* transitions with TDM oriented along the long axis of the BTBT core. These differences at the sulfur and carbon K-edges are understood through the MBAXS simulations that find a reordering of the energy of the lowest energy π* and σ* transitions at the sulfur K-edge due to the strong localization of the σ* orbital over the sulfur atom. This work highlights differences in the NEXAFS spectra of organic semiconductors at carbon and sulfur K-edges and provides new insights into peak assignment and x-ray dichroism relevant for studying the molecular orientation of organic semiconductor films.

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

Modeling of the metal–insulator transition temperature in alio-valently doped VO 2 through symbolic regression

The correlated semiconductor vanadium dioxide (VO 2 ) exhibits an insulator–metal transition (IMT) near room temperature, which is of interest in various device applications. Precise IMT temperature control is crucial to determine the use cases across technologies such as thermochromic windows, actuators for robots or neuronal oscillators. Doping the cation or anion sites can modulate the IMT by several tens of degrees and control hysteresis. However, modeling the effects of control parameters (e.g., doping concentration, type of dopants) is challenging due to complex experimental procedures and limited data, hindering the use of traditional data-driven machine learning approaches. Symbolic regression (SR) can bridge this gap by identifying nonlinear expressions connecting key input parameters to target properties, even with small data sets. In this work, we develop SR models to capture the IMT trends in VO 2 influenced by different dopant parameters. Using experimental data from the literature, our study reveals a dual nature of the IMT temperature with varying tungsten (W) doping concentrations. The symbolic model captures data trends and accounts for experimental variability, providing a complementary approach to first-principles calculations. Our feature-driven analysis across a broader class of dopants informs selectivity and provides qualitative insights into tuning phase transition properties valuable for neuromorphic computing and thermochromic windows.

36 MATERIALS SCIENCE↗
DOE OSTI2026

First-principles insights into Si substitution effects in Sm 2 (Fe,Si) 17 C x magnet

The partial substitution of Fe by Si enhances the phase stability of Sm 2 Fe 17 C x magnets with x > 1.0. Here, we elucidate the Si-substitution scheme and its impact on phase stability and magnetic properties in Sm 2 (Fe,Si) 17 C 3 from first-principles calculations and chemical bonding analysis. The calculated substitution energies for Si at various Fe sites are negative, indicating improved phase stability. Si preferentially substitutes Fe atoms at the 9d site in Sm 2 (Fe,Si) 17 C 3 while it tends to enter the Fe 18h site in Sm 2 (Fe,Si) 17 . This difference in site preference is attributed to the distinct chemical environments surrounding the Fe (Si) sites in the two compounds. Si substitution favors the formation of Sm–Si bonds while minimizing the Si–C and Si–Si interactions. Crystal orbital Hamilton populations and crystal orbital bond index calculations indicate that the partial replacement of Fe with Si strengthens the chemical bonding of Sm–Fe 3 (18f) and Sm–Fe 4 (18h) and improves overall phase stability in Sm 2 (Fe,Si) 17 C 3 . Beyond the dilution effect, Si substitution also reduces the magnetic moments of neighboring Fe atoms, a phenomenon linked to the strong Fe–Si bonding. These findings highlight the dual role of Si in modifying both the structural and magnetic characteristics of Sm 2 Fe 17 -based magnetic compounds.

Chemical bonding↗
DOE OSTI2026

Solvent effects on triplet yields in BODIPY-based photosensitizers

Here, we employ molecular dynamics simulations and quantum rate theories to elucidate the complex condensed-phase dynamics underpinning triplet-state formation in organic photosensitizers. Using models informed by first-principles calculations complete with a molecular representation of solvents of different polarities, we elucidate the interplay of the internal and environmental interactions underlying triplet yield. We find that triplet yields depend sensitively on the dielectric stabilization of the charge transfer intermediate that facilitates a transition into the triplet manifold. Our results illustrate the importance of molecularly detailed models in understanding the excited-state internal charge-transfer dynamics of photochemically relevant organic molecules.

Coello Escalante, Leonardo [University of Californ↗
DOE OSTI2025

Morphogenesis of spin cycloids in a noncollinear antiferromagnet

Pattern formation in spin systems with continuous-rotational symmetry (CRS) provides a powerful platform to study emergent complex magnetic phases and topological defects in condensed-matter physics. However, its understanding and correlation with unconventional magnetic order along with high-resolution nanoscale imaging are challenging. Here, we employ scanning nitrogen vacancy (NV) magnetometry to unveil the morphogenesis of spin cycloids at both the local and global scales within a single ferroelectric domain of (111)-oriented BiFeO 3 , which is a noncollinear antiferromagnet, resulting in formation of a glassy labyrinthine pattern. We find that the domains of locally oriented cycloids are interconnected by an array of topological defects and exhibit isotropic energy landscape predicted by first-principles calculations. We propose that the CRS of spin-cycloid propagation directions within the (111) drives the formation of the labyrinthine pattern and the associated topological defects such as antiferromagnetic skyrmions. Unexpectedly, reversing the as-grown ferroelectric polarization from [$\overline{1}$ $\overline{1}$ $\overline{1}$] to [111] produces a noncycloidal NV image contrast which ¯ could be attributed to either the emergence of a uniformly magnetized state or a reversal of the cycloid polarity. These findings highlight that (111)-oriented BiFeO 3 is not only important for studying the fascinating subject of pattern formation but could also be utilized as an ideal platform for integrating novel topological defects in the field of antiferromagnetic spintronics.

magnetoelectric↗
DOE OSTI2025

Deciphering decomposition pathways of high explosives with cryogenic X-ray Raman spectroscopy

We employed cryogenic X-ray Raman spectroscopy to investigate the early-stage decomposition of the high explosive molecule hexanitrohexaazaisowurtzitane (CL-20). By systematically varying the radiation dose under cryogenic conditions, we induced the decomposition of the molecule using ionizing radiation and observed the evolution of spectral features at the carbon, nitrogen, and oxygen K edges. Through extensive first-principles calculations, we identified key intermediates in the early stages of the decomposition process, resulting from C–C and C–N bond cleavage which leads to the opening of the internal cage structure. A detailed analysis of spectral trends and fingerprints provided evidence supporting N–NO 2 homolytic cleavage as the primary initial decomposition pathway. The combination of advanced core-level spectroscopy methods and state-of-the-art theoretical calculations enabled a comprehensive characterization of the molecular changes induced by controlled radiation dose exposures. In conclusion, our findings establish a benchmark for understanding the decomposition chemistry of high-explosive materials, offering important insights into their stability and reactivity under extreme conditions.

X-ray Raman spectroscopy↗
DOE OSTI2024

Lattice dynamics of 119 Sn impurity in a bcc-Cr crystal

The chromium crystal doped with 119 Sn isotope was studied using the nuclear resonance inelastic x-ray scattering and first principles calculations. The Sn partial phonon density of states (PDOS) was obtained for three temperatures that correspond to different magnetic states of Cr. At all temperatures, the energy spectrum consists of a broad band around 18 meV and a narrow peak at 43 meV. The additional peak around 39 meV is observed only in the magnetically ordered phases, indicating the influence of magnetic order in chromium on lattice dynamics. The partial PDOS calculated with the antiferromagnetic order on Cr atoms show a very good agreement with the experimental data. It is revealed that the high-energy peak is lying above the phonon spectra of the pure bcc-Cr crystal. These are the local modes with the increased energies due to a strongly reduced distance between Sn and the nearest-neighbor Cr atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Theoretical investigation of decoherence channels in athermal phonon sensors

The creation and evolution of nonequilibrium phonons is central in applications ranging from cosmological particle searches to decoherence processes in qubits. However, the fundamental understanding of decoherence pathways for athermal phonon distributions in solid-state systems remains an open question. Using first-principles calculations, we investigate the primary decay channels of athermal phonons in two technologically relevant semiconductors—Si and GaAs. We quantify the contributions of anharmonic, isotopic, and interfacial scattering in these materials. From this, we construct a model to estimate the thermal power in a readout scheme as a function of time. We discuss the implication of our results on noise limitations in current phonon sensor designs and strategies for improving coherence in next-generation phonon sensors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Large magnetoresistance and quantum oscillations observed in superconducting β -IrSn 4

Diverse quantum phenomena have been observed in TSn 4 transition metal (T) stannides, including superconductivity, nontrivial topology, and large magnetoresistance (MR) at low temperatures. Here, we report the experimental and theoretical investigation of tetragonal β-IrSn 4 (space group I 4 1 /acd) with the lattice parameters a = 6.362(3) Å and c = 22.723(0) Å. The temperature dependence of the electrical resistivity indicates that β-IrSn 4 is a good metal with Fermi-liquid behavior above the superconducting transition at T c ∼ 0.7 K. Magnetic susceptibility measurements indicate it is non-magnetic. However, positive transverse MR with linear field dependence is observed. At T= 2 K and Μ 0 H = 14 T, the MR reaches 700% without sign of saturation. Shubnikov–de Haas oscillations are observed from proximity detector oscillator measurements up to 60 T. We experimentally extract several frequencies, and through comparison to first-principles calculations, identify corresponding electronic bands. For H // c, two strong oscillations, F 1 ≈ 109 and F 2 ≈ 302 T, allow us to determine effective masses 0.266m 0 and 0.300m 0 (m 0 is the free-electron mass), respectively. The topology of these bands is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Role of chalcogen vacancies and hydrogen in the optical and electrical properties of bulk transition-metal dichalcogenides

Abstract Like in any other semiconductor, point defects in transition-metal dichalcogenides (TMDs) are expected to strongly impact their electronic and optical properties. However, identifying defects in these layered two-dimensional materials has been quite challenging with controversial conclusions despite the extensive literature in the past decade. Using first-principles calculations, we revisit the role of chalcogen vacancies and hydrogen impurity in bulk TMDs, reporting formation energies and thermodynamic and optical transition levels. We show that the S vacancy can explain recently observed cathodoluminescence spectra of MoS 2 flakes and predict similar optical levels in the other TMDs. In the case of the H impurity, we find it more stable sitting on an interstitial site in the Mo plane, acting as a shallow donor, and possibly explaining the often observed n-type conductivity in some TMDs. We also predict the frequencies of the local vibration modes for the H impurity, aiding its identification through Raman or infrared spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural mode coupling in perovskite oxides using hypothesis-driven active learning

Abstract Finding the ground-state structure with minimum energy is paramount to designing any material. In ABO 3 -type perovskite oxides with Pnma symmetry, the lowest energy phase is driven by an inherent trilinear coupling between the two primary order parameters such as rotation and tilt with antiferroelectric displacement of the A-site cations as established via hybrid improper ferroelectric mechanism. Conventionally, finding the relevant mode coupling driving phase transition requires performing first-principles calculations which is computationally time-consuming as well as expensive. It involves following an intuitive iterative hit and trial method of (a) adding two or multiple mode vectors, followed by (b) evaluating which combination would lead to the ground-state energy. In this study, we show how a hypothesis-driven active learning framework can identify suitable mode couplings within the Landau free energy expansion with minimal information on amplitudes of modes for a series of double perovskite oxides with A-site layered, columnar and rocksalt ordering. This scheme is expected to be applicable universally for understanding atomistic mechanisms derived from various structural mode couplings behind functionalities, for e.g. polarization, magnetization and metal–insulator transitions.

36 MATERIALS SCIENCE↗