SEARCH · Engineering Papers

Results for “first-principle 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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Pion and kaon PDFs from lattice QCD via large momentum effective theory and short-distance factorization

In this paper, we present a first-principles lattice-QCD calculation of the unpolarized quark PDF for the pion and the kaon. The lattice data rely on matrix elements calculated for boosted mesons coupled to non-local operators containing a Wilson line. The calculations on this lattice ensemble correspond to two degenerate light, a strange, and a charm quark (𝑁 𝑓 = 2 + 1 + 1), using maximally twisted mass fermions with a clover term. The lattice volume is 32 3 × 64, with a lattice spacing of 0.0934 fm, and a pion mass of 260 MeV. Matrix elements are calculated for hadron boosts of |𝑃 3 | = 0, 0.41, 0.83, 1.25, 1.66, and 2.07 GeV. To match lattice QCD results to their light-cone counterparts, we employ two complementary frameworks: the large-momentum effective theory (LaMET) and the short-distance factorization (SDF). Using these approaches in parallel, we also test the lattice data to identify methodology-driven systematics. Results are presented for the standard quark PDFs, as well as the valence sector. Beyond obtaining the PDFs, we also explore the possibility of extracting information on SU(3) flavor-symmetry-breaking effects. For LaMET, we also parametrize the momentum dependence to obtain the infinite-momentum PDFs. Since the present calculation is performed on a single ensemble at a pion mass of 260 MeV and fixed lattice spacing, the uncertainties reported are statistical only, and systematic uncertainties remain to be addressed in future multi-ensemble studies.

Miller, Joshua [Temple University, Philadelphia, P↗

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 OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗

DOE OSTI2026

Understanding and Predicting Pedestal Scenarios on NSTX-U (Final Technical Report)

This project investigated gyrokinetic instabilities in the NSTX pedestal, identified the major transport mechanisms (some of them novel) in the NSTX pedestal, and developed and validated a predictive modeling capability for pedestal transport in spherical tokamaks, with particular emphasis on NSTX discharges. The work combined first-principles gyrokinetic simulations, reduced transport models, integrated transport calculations with ASTRA, and exploratory machine learning tools. The central outcome is a practical modeling capability for pedestal temperature profiles based on reduced models informed and constrained by gyrokinetic physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Synergistic Doping and Stabilization of Magnetically Tunable LnTi3(Sb,Sn)4 (Ln: Ce–Gd) Kagome Metals

Here, we present our synthesis and characterization of the LnTi3(Sb,Sn)4 (Ln: Ce, Pr, Nd, Sm, Gd) family of cleavable kagome metals. While these materials are isostructural to the LnTi3Bi4 family, they only form as (Sb,Sn) solid-solutions with no corresponding LnTi3Sb4 or LnTi3Sn4 phases. We use a combination of first-principles density functional theory (DFT) and crystal orbital Hamilton population (COHP) calculations to show that (Sb,Sn) alloying has a stabilizing effect on the structure by adjusting the Fermi level, filling bonding states, depopulating antibonding states, and adjusting the density of states (DOS) toward local minima, an effect we call “synergistic doping.” Through a detailed characterization of the SmTi3(Sb,Sn)4 series, we further demonstrate that the tunable Fermi level has a profound effect on the physical properties. We observe multiple magnetic phases that stem from a competition between antiferromagnetic (AFM) and ferromagnetic-like (FM) ground states. Furthermore, the (Sb,Sn) ratio allows us to tune from direct AFM-FM competition toward a complex admixture with properties reminiscent of the recently discovered TbTi3Bi4. Ultimately, our work demonstrates how the concept of synergistic doping provides a means to stabilize new structures while developing systems with intrinsically tunable chemical, magnetic, and electronic properties.

Ortiz, Brenden [ORNL] (ORCID:0000000213337003)↗

DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗

DOE OSTI2026

Er Al :Al 2 ⁢O 3 for telecom-band photonics: Electronic structure and optical properties

Er-doped Al 2 ⁢O 3 is a promising host for telecom-band integrated photonics. Here, in this study, we combine ab initio calculations with a symmetry-resolved analysis to elucidate substitutional Er on the Al site (Er Al ) in 𝛼−Al 2 ⁢O 3 . First-principles relaxations confirm the structural stability of Er Al . We then use the local trigonal crystal-field symmetry to classify the Er-derived impurity levels by irreducible representations and to derive polarization-resolved electric-dipole selection rules, explicitly identifying the symmetry-allowed 𝑓−𝑑 hybridization channels. Kubo-Greenwood absorption spectra computed from Kohn-Sham states quantitatively corroborate these symmetry predictions. Furthermore, we connect the calculated intra-4⁢𝑓 line strengths to Judd-Ofelt theory, clarifying the role of 4⁢𝑓−5⁢𝑑 admixture in enabling optical activity. Notably, we predict a characteristic absorption near 1.47 µ⁢m (telecom band), relevant for on-chip amplification and emission. To our knowledge, a symmetry-resolved first-principles treatment of Er:Al 2 ⁢O 3 with an explicit Judd-Ofelt interpretation has not been reported, providing a transferable framework for tailoring rare-earth dopants in wide-band-gap oxides for integrated photonics. Our results for the optical spectra are in good agreement with experimental data. The resulting symmetry-based selection rules translate directly to polarization-dependent coupling in Al 2 ⁢O 3 integrated photonic waveguides and resonators, enabling device-level design of TE/TM-mode interaction with Er emitters in the 1.5-µ⁢m telecom band.

Er-doped Al2O3↗

DOE OSTI2026

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗

DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗

DOE OSTI2026

Charged Defects in UO 2 Bulk and Surface: A First-Principles Study

Uranium dioxide (UO 2 ) is the primary fuel used in nuclear reactors. Under the extreme heat and radiation inside a reactor, this material inevitably develops defects in its crystal structure. To investigate the nature and behavior of these defects, DFT+U calculations were employed to investigate charged point defects in both bulk UO 2 and its most stable surface, the (111) plane. The formation of defects and their impact on the electronic structure were systematically examined. The results reveal that these defects introduce localized electronic states, alter magnetic behavior, and modify the structural properties. In general, such defects act as deep traps capable of capturing and retaining charge carriers. The stability of these defects depends strongly on the chemical environment and the position of the Fermi level. Surface defect calculations reveal that oxygen vacancies form more readily at the surface than in the bulk over a wide range of electron chemical potential, with subsurface oxygen vacancies being more stable than those in the top layer. Overall, the findings demonstrate how charged defects influence magnetism, transport, and stability in UO 2 , providing insights that may guide improvements in the safety and efficiency of nuclear fuel.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Coupled-channel approach to isotensor π π π scattering from lattice QCD

The quest to understand three-body dynamics from first-principle QCD includes the study of nonresonant and resonant systems. The isospin I = 2 system is of particular interest having no three-body resonance but featuring a resonance in a subchannel, while also being a coupled-channel problem. In this study, we calculate the finite-volume spectrum from lattice QCD at two different pion masses, map the amplitude to the infinite volume through a generalized Finite-Volume Unitarity three-body quantization condition, investigate the limit of a narrow ρ , and compare with an effective Lagrangian prediction at leading order. Chiral extrapolations between different pion masses are performed.

Feng, Yuchuan [The George Washington University] (↗

DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

DOE OSTI2026

First-Principles Insights into Proton-Coupled Electron Transfer versus Hydrogen Evolution Reaction Selectivity from a Base-Appended Cobaltocene Mediator

Performing selective proton-coupled electron transfer (PCET) to substrates such as N 2 , CO 2 , and unsaturated organic molecules under electrochemical conditions requires the suppression of the competing hydrogen evolution reaction (HER). To address this challenge, our laboratory previously demonstrated a PCET mediator strategy using a dimethylaniline-appended cobaltocene complex, [(CpCoCp NMe2 )H] + , which performs selective reductive chemistry while suppressing the HER. However, the origin of the suppressed, yet still observable, HER has not been thoroughly established. In this work, we perform density functional theory (DFT) calculations to elucidate the HER mechanism involving this redox mediator and to provide atomistic insights into the bifurcation between the PCET and HER pathways. We find that protonation of the aniline moiety to form [CpCoCp NMe2H ] + is more favorable, both kinetically and thermodynamically, than formation of the ring-protonated species [(CpCo(Cp-H) NMe2 )] + . Furthermore, PCET to acetophenone is energetically more favorable via [CpCoCp NMe2H ] + than via [(CpCo(Cp-H) NMe2 )] +1/0 . In contrast, the most favorable HER pathway involves the ring-protonated Co(I) species. These results offer mechanistic insights into HER versus PCET bifurcation and establish guiding principles for designing PCET mediators for selective electroreductive transformations.

evolution reactions↗

DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗

DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗

DOE OSTI2026

First-Principles Discovery of Novel LiInP 2 S 6 Polymorphs with Promising Optoelectronic Responses

Recent advances in two-dimensional (2D) van der Waals (vdW) metal thiophosphates have attracted considerable attention due to their promising ionic conductivity, optical characteristics, and tunable physical properties. Within this material family, LiInP 2 S 6 has emerged as an intriguing candidate, not only because of its sensitivity to air and moisture but also due to its suitable band gap within the UV−vis range, enabling potential optoelectronic and photocatalytic applications. In this study, through comprehensive first-principles investigations, we unveil two previously unreported polymorphs of LiInP 2 S 6 in the monoclinic C2/c and trigonal P3̅1c (in-gap) space groups, in addition to examining the experimentally synthesized P3̅1c (in-layer) phase. Our studies identify the C2/c structure as the ground state, lying 9 meV per unit cell lower in energy than the experimentally realized trigonal P3̅1c (in-layer) phase. Further, we systematically examine the elastic, mechanical, thermodynamical, dynamical, electronic, and optical properties of all three polymorphs, confirming their mechanical, thermal, and dynamical stability. Notably, the P3̅1c (in-gap) phase exhibits enhanced stiffness, while the calculated indirect band gaps and strong photon absorption in the UV−vis range (∼3 eV) highlight the potential of the studied LiInP 2 S 6 phases for iontronic devices and optoelectronic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗