SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Machine Learning-Accelerated First-Principles Molecular Dynamics Reveals C–C Coupling Mechanisms toward Ethylene on Cu(100)

Here, the Cu(100) termination has been identified as the most effective facet for converting CO and CO 2 into ethylene. To enhance both the activity and selectivity of ethylene production, we perform machine-learning-accelerated, first-principles molecular dynamics simulations at 298 K in an explicit solvent at pH 7 to elucidate the C–C coupling mechanism─the critical reaction step in forming C 2+ products. Among the six potential C–C coupling pathways, the most feasible are CO* dimerization and CO – CHO* and CHO* – CHO* couplings. Using the computational hydrogen electrode method, we demonstrate that all three pathways are equally accessible at −0.6 V vs RHE. At a potential below −1.0 V vs RHE, the thermodynamic barriers for the CO – CHO* and CHO* – CHO* pathways become negligible. Our computational findings explain the experimental observations, particularly the absence of C 2+ products above −0.4 V vs RHE and the peaks in ethylene production near −0.6 and −1.0 V vs RHE. Since CHO* acts as a key intermediate common to both C–C coupling and CH 4 formation, we propose that suppressing CHO* hydrogenation would inhibit CH 4 pathways, thereby maximizing ethylene selectivity.

CO2 reduction↗

DOE OSTI2026

Database and deep-learning scalability of anharmonic phonon properties by automated brute-force first-principles calculations

Understanding the anharmonic phonon properties of crystal compounds—such as phonon lifetimes and thermal conductivities—is essential for investigating and optimizing their thermal transport behaviors. These properties also impact optical, electronic, and magnetic characteristics through interactions between phonons and other quasiparticles and fields. In this study, we develop an automated first-principles workflow to calculate anharmonic phonon properties and build a comprehensive database encompassing more than 6500 inorganic compounds. Utilizing this dataset, we train a graph neural network model to predict thermal conductivity values and spectra from structural parameters, demonstrating a scaling law in which prediction accuracy improves with increasing training data size. High-throughput screening with the model enables the identification of materials exhibiting extreme thermal conductivities—both high and low. The resulting database offers valuable insights into the anharmonic behavior of phonons, thereby accelerating the design and development of advanced functional materials.

Ohnishi, Masato [University of Tokyo (Japan); Inst↗

DOE OSTI2025

Tunable energy landscape of screw dislocation cores by compositional fluctuations in bcc high-entropy alloys from first-principles calculations

The energy landscape of screw dislocation cores plays a central role in dislocation-mediated deformation mechanisms in body-centered cubic (bcc) metals. In bcc high-entropy alloys (HEAs), this energy landscape is modulated by local compositional fluctuations, which has important implications for deformation processes in these materials. Through first-principles calculations, this study investigates high-symmetry screw dislocation core structures in NbTaMoW and NbTaTiHf bcc HEAs. The results show that alloying group IV transition metals lead to large local lattice distortions at dislocation cores, which is demonstrated to be an important factor governing fluctuations in core configurations along a dislocation line. Importantly, group IV elements near the core induce features in the energy landscape that are exclusive for HEAs, specifically lowering the energy of core configurations that are unstable in elemental bcc metals. A combined influence of these chemical effects with crystallographic details enables the activation of glide planes, a feature that has been linked to ductility improvements in bcc HEAs. These findings provide new insights into the atomic-scale mechanisms underlying dislocation mobility in bcc HEAs, offering a pathway for designing materials with tailored mechanical properties.

Borges, Pedro P P O↗

DOE OSTI2025

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗

DOE OSTI2024

First-principles study of the magneto-Raman effect in van der Waals layered magnets

Magneto-Raman spectroscopy has been used to study spin-phonon coupling in two-dimensional (2D) magnets. Raman spectra of CrI 3 show a strong dependence on the magnetic order within a layer and between the layers. Here we carry out the first systematic theoretical investigation of the magneto-Raman effect in 2D magnets by performing density functional theory calculations and developing a generalized polarizability model. Our first-principles simulations well reproduce experimental Raman spectra of CrI 3 with different magnetic states. The model reveals how the change of spin orientation in each layer is coupled to the layer’s vibration to induce or eliminate the spin-dependent anti-symmetric off-diagonal terms in the Raman tensor for altering the selection rules. We also uncover that the correlation between phonon modes and magnetic orders is a universal phenomenon, which should exist in other phonon modes and 2D magnets. Our predictive simulations and modeling are expected to guide the research in 2D magnets.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Effects of interstitial oxygen on ω transformations and twin formation in bcc NbTaTiHf multi-principal element alloy from first-principles

Transformation- and twinning-induced plasticity (TRIP and TWIP) have been reported to contribute to the low-temperature deformation of some body-centered cubic (bcc) multi-principal element alloys (MPEAs) containing large fractions of group IV transition metals. The influence of interstitial solutes on the mechanisms underlying these forms of plasticity, however, remains unclear. Using first-principles calculations, we study the effects of interstitial O atoms on the relative stability of bcc and ω phases and on unstable and twin boundary stacking fault energy profiles in a representative bcc MPEA with high group-IV elemental fraction: NbTaTiHf. We find that O additions generally promote the relaxation of ω configurations back to their parent bcc structure, therefore inhibiting ω transformation. Calculations of the Rice parameter for bulk bcc and phases, as well as bcc-ω interfaces, further show that ω formation is a potent embrittlement factor, an effect that is enhanced by O additions, suggesting that the formation of bcc-ω interfaces is energetically preferred over the formation of the bulk ω phase. By contrast, the Rice parameter for twin boundaries indicates that these interfaces do not embrittle the material, even with O atoms at twin boundaries, providing a more favorable pathway for plastic deformation compared to ω transformation.

Density functional theory↗

DOE OSTI2025

First-principles investigation of the resistive switching energetics in monolayer MoS 2 : insights into metal diffusion and adsorption

A deeper understanding of resistive switching (RS) in 2D materials is essential for advancing neuromorphic computing. The Dissociation-Diffusion-Adsorption (DDA) model offers a useful framework for probing RS mechanisms in non-volatile memory (NVM) and in-memory computing. We have employed first-principles density functional theory (DFT) to explore dissociation, diffusion, and adsorption phenomena within the DDA model, focusing on the interactions between exemplary metal atoms (Au, Ag, Cu) and monolayer MoS 2 . Nudged elastic band (NEB) calculations evaluated diffusion barriers in pristine and sulfur-vacancy MoS 2 . Charged systems were modeled to assess the impact of applied bias on migration pathways. We also examined metal dissociation from bulk electrodes and adsorption at S vacancies. Ag/MoS 2 shows the lowest dissociation barrier (~0.034 eV), while Au and Cu exhibit similar values (~0.32 eV). These insights highlight Ag as a promising candidate for low-energy RS applications and provide guidance for optimizing switching efficiency in 2D memory devices.

Atomistic models↗

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 OSTI2025

First-principles computations of the Stark shift of a defect-bound exciton: The case of the T center in silicon

The T center in silicon has recently drawn a lot of attention for its potential in quantum information science. The sensitivity of the zero-phonon line (ZPL) to electrical field was recently investigated by a combination of different experimental methods but there are still few first principles studies on the Stark shift of the T center. Dealing with the defect-bound exciton nature of the excited state is particularly challenging using density functional theory because of the large spatial delocalization associated with the wavefunction. Here, in this work, we tackle this issue by performing a convergence study over the supercell size. We obtain an exciton binding energy of 28.5 meV, in good agreement with experimental results. We then calculate the Stark shift through the dipole moment change of the ZPL transition of the T center using the modern theory of polarization formalism and find a modest linear coefficient of Δ⁢𝜇=0.79⁢D along X and Δ⁢𝜇=0.03⁢𝐷 along Y. We discuss our results in light of the recent experimental measurements of the Stark shift. Our analysis suggests that bound-exciton defects could be particularly sensitive to local field effect as a result of their large spatial extent.

color centers↗

DOE OSTI2026

Chemical bonding, phase stability and magnetic property in Sm 2 Fe 17 X 3 (X=H, C, N): A first-principles perspective

As a promising alternative to Nd–Fe–B magnets, the critical rare earth free Sm 2 Fe 17 X 3 (X = C, N) exhibits potential for high-performance magnets. However, their poor phase stability remains a major obstacle to developing bulk magnets. We investigated the phase stability and intrinsic magnetic properties of Sm 2 Fe 17 X 3 (X = H, C, N) using first-principles calculations and chemical bond analysis. The formation energies are negative, while the decomposition energies are −1.53, 0.348, and −0.74 eV per formula unit for X = H, C, and N, respectively, which is responsible for the weak thermal stability. Our chemical bond analysis reveals that the bonding asymmetry between Sm–X and Fe–X interactions creates local structural distortions and degrades the phase stability of Sm 2 Fe 17 X 3 . The project Crystal Orbital Hamilton Population (-pCOHP) analysis indicates that the Sm–X bonding remains positive up to the Fermi level, indicating stable bonding interactions. Here, in contrast, the Fe–X bonding becomes negative near the Fermi level, signifying anti-bonding contributions that reduce structural stability. Interstitial atoms X expand the lattice and enhance Fe magnetic moments, but Fe–X bonding suppresses neighboring Fe moments. Electron transfer from Sm to X modifies the valence state of Sm and the crystal field at the site, contributing to enhanced magnetocrystalline anisotropy in Sm 2 Fe 17 X 3 . Among the interstitial elements, carbon and nitrogen—with their larger atomic radius and higher electronegativity—induce greater lattice expansion and form stronger bonds with neighboring Sm and Fe atoms compared to hydrogen. Consequently, Sm 2 Fe 17 X 3 (X = C and N) exhibits better phase stability and significant improvement in magnetic properties.

Chemical bonding↗

DOE OSTI2025

Exciton thermalization dynamics in monolayer MoS2: A first-principles Boltzmann equation study

Understanding exciton thermalization is critical for optimizing optoelectronic and photocatalytic processes in many materials. However, it is hard to access the dynamics of such processes experimentally, especially on systems such as monolayer transition metal dichalcogenides, where various low-energy excitations pathways can compete for exciton thermalization. Here, we study exciton dynamics due to exciton-phonon scattering in monolayer MoS2 from a first-principles, interacting Green's function approach, to obtain the relaxation and thermalization of low-energy excitons following different initial excitations at different temperatures. We find that the thermalization occurs on a picosecond time scale at 300 K but can increase by an order of magnitude at 100 K. The long total thermalization time, owing to the nature of its excitonic band structure, is dominated by slow spin-flip scattering processes in monolayer MoS2. In contrast, thermalization of excitons in individual spin-aligned and spin-anti-aligned channels can be achieved within a few hundred fs when exciting higher-energy excitons. We further simulate the intensity spectrum of time-resolved angle-resolved photoemission spectroscopy experiments and anticipate that such calculations may serve as a map to correlate spectroscopic signatures with microscopic exciton dynamics.

Chan, Yang-hao↗

DOE OSTI2025

First-Principles study of interface capture of Ni during high-temperature oxidation of NiCr alloys

The formation of precipitated phases of less reactive metals within an oxide scale during alloy oxidation can disrupt its structural integrity, significantly compromising the protective function of the oxide. To better understand this phenomenon, it is crucial to investigate the atomic-level diffusion mechanisms that drive precipitate formation. In this study, we employ first-principles calculations to examine the diffusion and aggregation of Ni atoms across the NiCr/Cr 2 O 3 interface during the high-temperature oxidation of NiCr alloys. Our results reveal that, in the absence of Cr vacancies, Ni atoms predominantly remain within the NiCr substrate, with minimal migration toward the Cr 2 O 3 layer. However, when Cr vacancies are present at both the NiCr/Cr 2 O 3 interface and within the Cr 2 O 3 bulk, Ni atoms exhibit a strong propensity to migrate from the NiCr substrate into the Cr 2 O 3 bulk, where they aggregate into Ni-rich regions with a corundum structure. This study provides valuable atomic-level insights into the diffusion and aggregation of less-reactive metals at the alloy/oxide interface, enhancing our understanding of solute capture and its contribution to the degradation of oxidation resistance in alloys.

Chen, Jianmin [Xiangtan University (China)] (ORCID↗

DOE OSTI2025

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗

DOE OSTI2025

First-principles elucidation of defect-mediated Li transport in hexagonal boron nitride

Hexagonal boron nitride (hBN) is a promising candidate as a protective membrane or separator in Li-ion and Li–S batteries, given its excellent chemical stability, mechanical robustness, and high thermal conductivity. In addition, hBN can be functionalized by introducing defects and dopants, or be directly integrated into other active components of batteries, which further augments its appeal to the field. Here, we use first-principles simulations to evaluate the role of atomic defects in hBN in regulating the Li-ion diffusion mechanism and associated kinetics. Specifically, the following four distinct types of vacancy defects are considered: isolated single B and N vacancies, a B–N vacancy pair, and a B 3 N vacancy cluster. It is found that these defect sites generally favor Li intercalation and out-of-plane diffusion but slow down in-plane Li-ion diffusion due to a strong Li trapping effect at the defect sites. Such a trapping effect is, however, highly local such that it does not necessarily affect the overall Li-ion conductivity in defected hBN layers. The present systematic evaluation of the impact of atomic defects on Li ion migration and accompanied charge analysis of hBN lattice in response to Li-ion diffusion provide a mechanistic understanding of Li-ion transport behavior in defected hBN and highlight the potential of defect engineering to achieve optimal material performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Chemical trends favoring interstitial cluster formation in bcc high-entropy alloys from first-principles calculations

Achieving high strength and ductility is a common goal in structural alloy design. Body-centered cubic high-entropy alloys (HEAs) commonly highlight the conflict between these properties, with stronger alloys being brittle and vice versa. Recent reports suggest interstitial solutes can be used to overcome this trade-off, in some cases providing both strength and ductility enhancements. This effect has been correlated with interstitial cluster formation, although the conditions favoring their formation remain incompletely understood. Using first-principles calculations of solution energies and diffusivities, we provide insights into thermodynamic and kinetic factors favoring interstitial solute clusters. Among C, N and O solutes, O interstitials display most desirable diffusion kinetics. Further, the results highlight the importance of local composition fluctuations in the HEAs to enable the formation of clusters of appreciable size. The results are explained in terms of bonding and distortion trends across solutes and HEA compositions to provide guidelines for alloy design.

Borges, Pedro P P O↗

DOE OSTI2026

First-principles investigation of structure-property relationships in stable and metastable MXenes

Understanding the structure–property relationships in layered transition-metal carbides or nitrides, known as MXenes, is of critical importance for their rational design, synthesis, and application. However, the vast chemical and structural diversity of MXenes, stemming from their wide range of M and X elements, surface terminations, and different atomic coordination environments, makes it challenging to clearly understand these structure–property relationships. Here, in this work, we perform first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to comprehensively investigate the stability and a variety of physical properties of MXenes with different coordination environments. Using Ti- and Mo-based carbide MXenes as model systems, energetic calculations reveal that Ti-based MXenes are most stable in octahedral coordination, whereas Mo-based MXenes preferentially adopt prismatic coordination. This fundamental difference in preferred atomic coordination gives rise to markedly distinct properties between these two systems as a function of the fraction of octahedral and prismatic sites. For instance, the in-plane stiffness of Ti-based MXenes increases as octahedral coordination becomes dominant, but it decreases in the Mo-based MXenes under the same conditions. Additional stability analyses based on mechanical, lattice-dynamical, and temperature-dependent thermodynamic properties demonstrate that many metastable MXenes not only satisfy the strict stability criteria but can also undergo phase transitions among different structures and even become stabilized at elevated temperatures. Although surface terminations, such as F and O atoms, do not alter the energetic ordering or the overall stiffness trends among stable and metastable MXenes, they influence other material properties. For instance, O termination can induce semiconducting behavior in both stable and metastable Ti 2 ⁢CO 2 MXenes. This study significantly advances the fundamental understanding of structure–property relationships in MXenes and provides valuable guidance for developing coordination-based design principles to precisely engineer MXenes with improved properties.

Oyeniran, Noah [University of Alabama, Tuscaloosa,↗

DOE OSTI2025

First-principles electron-phonon interactions and polarons in the parent cuprate La 2 CuO 4

Understanding electronic interactions in high-temperature superconductors is an outstanding challenge. In the widely studied cuprate materials, experimental evidence points to strong electron-phonon ( e -ph) coupling and broad photoemission spectra. Yet, the microscopic origin of this behavior is not fully understood. Here, we study e -ph interactions and polarons in a prototypical parent (undoped) cuprate, La 2 CuO 4 (LCO), by means of first-principles calculations. Leveraging parameter-free Hubbard-corrected density functional theory, we obtain a ground state with the band gap and Cu magnetic moment in nearly exact agreement with experiments. This enables a quantitative characterization of e -ph interactions. Our calculations reveal two classes of longitudinal optical (LO) phonons with strong e -ph coupling to hole states. These modes consist of bond stretching and bond bending in the Cu-O plane as well as vibrations of apical O atoms. The hole spectral functions, obtained with a cumulant method that can capture strong e -ph coupling, exhibit broad quasiparticle peaks with a small spectral weight ( Z ≈ 0.25 ) and pronounced LO-phonon sidebands characteristic of polaron effects. Our calculations predict features observed in photoemission spectra, including a 40-meV peak in the e -ph coupling distribution function not explained by existing models. These results show that the universal strong e -ph coupling found experimentally in doped lanthanum cuprates is also present in the parent compound, and elucidate its microscopic origin. Published by the American Physical Society 2025

Chang, Benjamin K. (ORCID:0000000313049324)↗

NASA NTRSDate not supplied

First-Principles Investigation of Phase Transitions in RuNb, RuTa, and Ru-Nb-Ta Ultrahigh Temperature Shape Memory Alloys

Ultrahigh temperature shape memory alloys (UHT-SMAs) have transition temperatures above 600 C, and they have found applications for sensing and actuating devices in aerospace industry. Among very few such UHT-SMAs currently known are Ru-based alloys such as RuNb and RuTa, whose martensite structures and phase transitions are totally different from those of NiTi-based SMAs and were poorly understood. In this work, we carried out a systematical study of RuNb using first-principles total energy calculations and molecular dynamics (MD) simulations. The transition paths and mechanisms in cubic → tetragonal →monoclinic transitions are revealed. and the transition sequence and martensitic transition temperatures are determined (MTTs) by evaluating the Gibbs free energies using thermodynamic integration. The calculated MTTs are in very good agreement with the experimental data. We found that the monoclinic phase at the second transition has the 𝑃21/m symmetry instead of experimentally identified 𝑃2/m. Our calculations demonstrate that RuTa has very similar phase transitions to those of RuTa. Furthermore, we studied the Ru0.5Nb0.25Ta0.25 ternary. Our results can explain the measured significant drop of MTT (~ 200 C) in the second transition for Ru0.5Nb0.25Ta0.25 compared with those of RuNb and RuTa, while in the first transition its MTT is between MTTs of RuNb and RuTa. The insights gained by this study and the verified ab initio methods for accurate MTT calculations can be applied to fast screen and quantitatively design novel UHT-SMAs having similar properties with desirable MTTs and much reduced cost.

Zhigang Wu↗