Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Importance of enforcing Hund’s rules in density functional theory calculations of rare earth magnetocrystalline anisotropy

Density functional theory (DFT) and its extensions, such as DFT+U and DFT+dynamical mean-field theory, are invaluable for studying magnetic properties in solids. However, rare-earth (R) materials remain challenging due to self-interaction errors and the lack of proper orbital polarization. We show how the orbital dependence of self-interaction error contradicts Hund’s rules and plagues magnetocrystalline anisotropy (MA) calculations, and how analyzing DFT states that respect Hund’s rules can mitigate this issue. We benchmark MA in RCo 5 , R 2 Fe 14 B, and RFe 12 , extending prior work on RMn 6 Sn 6 , achieving excellent agreement with experiments. Additionally, we illustrate a semi-analytical perturbation approach that treats crystal fields as a perturbation in the large spin-orbit coupling limit. Using Gd-4f crystal-field splitting, this method provides a microscopic understanding of MA and enables rapid screening of high-MA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Exploring Anomalous Photoelectron Angular Distributions in the Photoelectron Spectra of Gd 3 O 3 – : Study of Gd 3 O 2 – and Gd 3 O 3 – Using Photoelectron Spectroscopy and Density Functional Theory Calculations

Anion photoelectron (PE) spectra of lanthanide oxide clusters obtained previously have exhibited anomalous photoelectron angular distributions which were attributed to strong PE–valence electron (PEVE) interactions. Here, to further explore this effect, we have obtained the PE spectra of Gd 3 O 2 – and Gd 3 O 3 – , two clusters that have similarly complex electronic structures but contrasting symmetries. The spectra exhibit manifolds of detachment transitions at similar binding energies in a 0.5 eV window of energy. The electron affinity of Gd 3 O 2 is measured to be 1.29 ± 0.05 eV, and that of Gd 3 O 3 is 1.31 ± 0.05 eV. As seen in previous studies on lanthanide oxide cluster anions in lower than conventional oxidation states, transitions in spectra obtained lower photon energies are more congested than those obtained with higher photon energy, a signature of strong PEVE interactions. While the detachment transitions have predominantly parallel photoelectron angular distributions (PAD), the PAD varies across the manifold of transitions in the PE spectrum of Gd 3 O 3 – in a way that suggests four different subgroups of transitions. Results of calculations on Gd 3 O 2 – suggest kite or V-shape structures with antiferromagnetic coupling between one of the 4f 7 subshells with the two others. Calculations on Gd 3 O 3 – more definitively point to ring structures with a nearly isoenergetic ferromagnetically coupled high spin (24-tet) state and a dectet state in which one of the 4f 7 subshells is antiferromagnetically coupled with the other two. Taking these results as qualitative, we propose that strong mixing between the unperturbed states predicted computationally leads to overlapping transitions with different PADs.

anions↗
DOE OSTI2025

Unconventional Anomalous Hall Effect in Hexagonal Polar Magnet Y 3 Co 8 Sn 4

A rare realization of unconventional anomalous Hall effect (UAHE) both below and above the magnetic transition temperature (T C ) in a hexagonal noncentrosymmetric magnet Y 3 Co 8 Sn 4 , using a combined experimental and ab initio calculations. Occurrence of such UAHE is mainly attributed to the reciprocal ($\mathscr{KS}$) topology (i.e., the presence of topological Weyl points at/near the Fermi level), along with some contribution from the topological magnetic texture (at low temperatures), as inferred from the measured field-dependent ac susceptibility. The effect of UAHE on the measured transport behavior however evolves differently with temperature above and below T C = 53 K, suggesting different physical mechanism responsible in the two phases. A unique planar ferrimagnetic ordering is found to be the most stable state with ab-plane as the easy plane below T C . The simulated net magnetization and the moment per Co atom agrees fairly well with the experimentally measured values. A reasonably large anomalous Hall conductivity (AHC) is also observed in both the phases (above and below T C ) of the present compound, which is again not so ubiquitous. In conclusion, the results underscore the family of R 3 Co 8 Sn 4 (R = rare earth) polar magnets as a compelling backdrop for exploring the synergy of magnetism and topological non-trivial electronic states, pivotal for spintronic applications.

Density functional theory calculations↗
DOE OSTI2026

Modeling an S N 2 Reaction Mechanism for Hydrolysis of Siloxane Linkages with Density Functional Theory under Basic Conditions and Implications for Dissolution of Quartz

An extended silicate molecular cluster was used to perform density functional theory calculations to determine if an S N 2 mechanism, where OH – attacks a Q 1 Si coupled with Na + charge balancing an adjacent siloxane bridging O atom, could explain the observed energy of activation of quartz dissolution under basic conditions.

chemical calculations↗
DOE OSTI2025

d -Electron Heavy Fermion Behavior in a Near-Room-Temperature Polar Metallic Ferrimagnet: A Case of Mn 5 SiC

Polycrystalline Mn 5 SiC was synthesized by using a high-temperature solid-state method. Mn 5 SiC adopts a polar space group (Cmc2 1 ) with six crystallographic Mn sites confirmed by X-ray and neutron diffraction, transmission electron microscopy, and second harmonic generation experiments. The complex crystal structure features edge-sharing trigonal prisms and icosahedra, as well as face/edge-sharing pentagonal prisms. Magnetic measurements indicate ferrimagnetic ordering with a transition temperature of 284 K. The ferrimagnetic structure (magnetic space group Cm’c’2 1 ) was further identified by powder neutron diffraction, where collinear Mn spins align along the crystallographic c-axis. The refined magnetic moment for each crystallographic Mn site at 4 K is 1.8(2), −2.42(9), −1.72(8), 0.51(6), 0.50(4), and 1.7(2) μ B . Density functional theory calculations confirm both the metallic behavior and the ferrimagnetic structure observed experimentally and further provide insight into the observed Mn moment dependence across crystallographic sites. The resistivity and specific heat measurements and density functional theory calculations reveal a substantially large Kadowaki–Woods ratio of 5 × 10 –5 μΩ·cm/(mJ/mol) 2 and a many-body renormalization factor of 5.5, indicating the unusual heavy Fermion behavior in such an itinerant magnetic metal.

Messegee, Zachary T. [George Mason University, Fai↗
DOE OSTI2025

Constraints on magnetism and correlations in RuO 2 from lattice dynamics and Mössbauer spectroscopy

We provide experimental evidence for the absence of a magnetic moment in bulk RuO 2 , a candidate altermagnetic material, by using a combination of Mössbauer spectroscopy, nuclear forward scattering, inelastic X-ray and neutron scattering, and density functional theory calculations. Using complementary Mössbauer and nuclear forward scattering, we determine the Ru magnetic hyperfine splitting to be negligible. Inelastic X-ray and neutron scattering-derived lattice dynamics of RuO 2 are compared to density functional theory calculations of varying flavors. Comparisons among theory with experiments indicate that electronic correlations, rather than magnetic order, are key in describing the lattice dynamics.

Mössbauer spectroscopy↗
DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗
DOE OSTI2026

Confinement of quasi-atomic structures in Ti 2 N and Ti 3 N 2 MXene electrides

Metal carbides, nitrides, or carbonitrides of early transition metals, better known as MXenes, possess notable structural, electrical, and magnetic properties. Analyzing electronic structures by calculating structural stability, band structure, density of states, Bader charge transfer, and work functions utilizing first principle calculations, we revealed that titanium nitride MXenes, namely Ti 2 N and Ti 3 N 2 , have excess anionic electrons in their lattice voids, making them MXene electrides. Bulk Ti 3 N 2 has competing antiferromagnetic (AFM) and ferromagnetic (FM) configurations with slightly more stable AFM configuration, while the Ti 2 N MXene is nonmagnetic. Although Ti 3 N 2 favors AFM configuration with hexagonal crystal systems having 6/ mmm point group symmetry, Ti 3 N 2 does not support altermagnetism. The monolayer of the Ti 3 N 2 MXene is a ferromagnetic electride. These unique properties of having non-nuclear interstitial anionic electrons in the electronic structure of titanium nitride MXene have not yet been reported in the literature. Density functional theory calculations show TiN is neither an electride, MXene, or magnetic.

Anionic electrons↗
DOE OSTI2026

Ultrafast One-Dimensional Peierls-Distortion Dynamics in 1⁢T′−Re⁢S 2 Revealed by 4D Electron Microscopy

Rhenium disulfide (ReS 2 ), a prototypical 2D semiconductor with an anisotropic 1⁢T′ structure due to the pronounced Peierls distortion, has demonstrated great potential for polarization-dependent optoelectronics and photonics. Here, we report an ultrafast phase transition occurring within 1 ps, accompanied by a one-dimensional Peierls-distortion relaxation in 1⁢T′−ReS 2 revealed with combined 4D electron microscopy and time-dependent density functional theory calculations. Upon femtosecond laser pulse excitation, the Re-Re dimerization morphology rapidly transforms from a diamond cluster to zigzag chains and persists for several nanoseconds. Here, this ultrafast Peierls-distortion relaxation is further verified by transient changes in optical anisotropy via polarization-dependent transient absorption spectroscopy. Time-dependent density functional theory calculations attribute this novel phase transition to strong correlation between Peierls distortion and impulsive photoexcited carrier doping, predicting a transient band-gap collapse. This Letter opens up an exciting avenue for ultrafast control of Peierls-distortion-induced anisotropy and the metal-insulator transition in 1⁢T′−ReS 2 .

1T’-ReS2↗
DOE OSTI2024

Design and Synthesis of PtPdNiCoMn High‐Entropy Alloy Electrocatalyst for Enhanced Alkaline Hydrogen Evolution Reaction: A Theoretically Supported Predictive Design Approach

Electrocatalytic hydrogen generation requires a multifunctional electrocatalyst with abundant active sites to drive multielectron transfer reactions. High entropy alloys (HEA) are five or more-elements with high configurational entropy are considered unique materials for next-generation electrocatalysts. Here, in this work, based on new screening guidelines for catalyst selections that combine density-functional theory calculated Gibbs formation-enthalpy with bond length and electronegativity variance, a novel HEA electrocatalyst consisting of five elements, namely, Pt, Pd, Ni, Co, and Mn has been designed. By simple room temperature electrodeposition, the designed catalyst is prepared and its hydrogen evolution reaction (HER) is explored and validated through experimental and theoretical approaches. The HEA demonstrated a superior HER activity with an overpotential of 22.6 mV at -10 mA cm -2 which outperforms Pt/C commercial catalyst. No evident degradation of the material is detected even after 100 hours of continuous operation under high current density. Moreover, the HEA has shown exceptional performance in harsh electrolyte conditions such as in simulated seawater and actual seawater. Remarkably, the density-functional theory calculated Gibbs formation-enthalpy is small (≈0 eV) compared to Pt/C placing the new HEA near the apex of Trasatti's model of Volcano plot, which is also suggestive of superior HER activity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Oxidative Deboronation of Boronic Acids by Hydrogen Peroxide in Planta Generates Borate for Cross-Linking of Rhamnogalacturonan II

Vascular plants require boron to cross-link the rhamnogalacturonan-II (RG-II) domain of pectin to form functional cell walls. Boronic acids, which form reversible esters with cis-diols like borate, have been proposed to influence RG-II cross-linking, though the mechanism remains unclear. We used suspension-cultured rose cells adapted to grow without boron to investigate the effect of boronic acids on RG-II dimerization. When grown with phenylboronic acid (PBA) as the sole boron source, nearly all RG-II was crosslinked, whereas methylboronic acid (MBA) only partially restored cross-linking. In contrast, in vitro assays showed that homogeneous RG-II monomers did not dimerize with alkyl or aryl boronic acids unless supplemented with hydrogen peroxide (H2O2), which oxidatively converts boronic acids to boric acid. Real-time NMR spectroscopy and density functional theory calculations provided insight into the reaction mechanism and energetics of oxidation respectively. Together, our data show that exogenous boronic acids are a source of boric acid for plants, and that the deboronation reaction generates aryl or alkyl alcohol byproducts that can undergo further chemical modification in planta. The fate and potential roles of these byproducts in planta remain to be determined.

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

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗