Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electronic structure theory”

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

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

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

Is the Matrix Completion of Reduced Density Matrices Unique?

Reduced density matrices are central to describing observables in many-body quantum systems. In electronic structure theory, the two-particle reduced density matrix (2-RDM) suffices to determine the energy and other key properties. Recent work has used matrix completion, leveraging the low-rank structure of RDMs and approximate theoretical models, to reconstruct the 2-RDM from partial data and thus reduce the computational cost. However, matrix completion is, in general, an under-determined problem. Revisiting Rosina’s theorem (Rosina, M. Queen’s Papers on Pure and Applied Mathematics , 1968, No. 11, 369), we here show that the matrix completion is unique under certain conditions, identifying the subset of 2-RDM elements that enables its exact reconstruction from incomplete information. Building on this, we introduce a hybrid quantum–stochastic algorithm that achieves exact matrix completion, demonstrated through applications to the Fermi–Hubbard model.

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

Discovering Ni/Cu Single-Atom Alloy as a Highly Active and Selective Catalyst for Direct Methane Conversion to Ethylene: A First-Principles Kinetic Study

Direct methane conversion to liquid fuels or value-added chemicals is a promising technology to utilize natural resources without resorting to further petroleum extraction. However, discovering efficient catalysts for this reaction is challenging due to either coke formation or unfavorable C-H bond activation. Herein, we design single-atom alloy (SAA) catalysts to simultaneously eliminate the above two bottlenecks based on mechanism-guided strategies: (1) the active single atom enables favorable C-H bond breaking and (2) the less reactive host metal facilitates C-C coupling and thus avoids strong binding of carbonaceous species. Employing electronic structure theory calculations, we screened the stability of multiple SAAs with 3d-5d transition metals atomically dispersed on a copper surface in terms of avoiding dopant aggregation and segregation. We then evaluated reactivities of the stable SAAs as catalysts for direct methane conversion to C2 products, including methane dehydrogenation and C-C coupling mechanisms. Combining selectivity analysis with kinetic modeling, we predicted that nickel dispersed on copper, i.e., Ni/Cu SAA, is a highly active and selective catalyst that can efficiently transform methane to ethylene. This work designs efficient SAA catalysts for direct methane activation and provides chemical insights into engineering compositions of SAAs to tune their catalytic performances.

Kothakonda, Manish↗
DOE OSTI2025

Unconventional low temperature decomposition of a saturated hydrocarbon over atomically-dispersed titanium-aluminum-boron catalyst

Sonochemically-synthesized atomically-dispersed titanium-aluminum-boron nanopowder (TiAlB NP) exhibits a remarkable low-temperature catalytic activation of aliphatic C-H bonds at 750 K followed by C-C bond activation thus emerging as a potent low-cost alternative to expensive platinum group metals. Here, the model saturated hydrocarbon, exo-tetrahydrodicyclopentadiene (C10H16), undergoes catalytic decomposition on TiAlB NPs in a chemical microreactor to produce 1,3-cyclopentadiene (c-C5H6), cyclopentene (c-C5H8), and molecular hydrogen (H2) as detected in situ via isomer-selective, single-photon ionization time-of-flight mass spectrometry. Extensive electronic structure theory calculations on model clusters of the catalyst decode a unique synergy among the atomic constituents of the catalyst and chemical bonding in this stepwise, retro Diels Alder reaction: Ti, although insensitive to C-H activation in its metallic state, initiates the catalysis via chemisorption of the hydrocarbon, adjacent B centers readily abstract hydrogen atoms and store them during the catalytic cycle, while Al stabilizes the catalyst structure yet providing space for critical docking sites for the departing hydrocarbons.

Biswas, Souvick↗
DOE OSTI2025

A foundation model for atomistic materials chemistry

Atomistic simulations of matter, especially those that leverage first-principles (ab initio) electronic structure theory, provide a microscopic view of the world, underpinning much of our understanding of chemistry and materials science. Over the last decade or so, machine-learned force fields have transformed atomistic modeling by enabling simulations of ab initio quality over unprecedented time and length scales. However, early machine-learning (ML) force fields have largely been limited by (i) the substantial computational and human effort required to develop and validate potentials for each particular system of interest and (ii) a general lack of transferability from one chemical system to the next. Here, we show that it is possible to create a general-purpose atomistic ML model, trained on a public dataset of moderate size, that is capable of running stable molecular dynamics for a wide range of molecules and materials. We demonstrate the power of the MACE-MP-0 model-and its qualitative and at times quantitative accuracy-on a diverse set of problems in the physical sciences, including properties of solids, liquids, gases, chemical reactions, interfaces, and even the dynamics of a small protein. The model can be applied out of the box as a starting or "foundation" model for any atomistic system of interest and, when desired, can be fine-tuned on just a handful of application-specific data points to reach ab initio accuracy. Establishing that a stable force-field model can cover almost all materials changes atomistic modeling in a fundamental way: experienced users obtain reliable results much faster, and beginners face a lower barrier to entry. Foundation models thus represent a step toward democratizing the revolution in atomic-scale modeling that has been brought about by ML force fields.

Batatia, Ilyes↗
DOE OSTI2025

Many-body expansion based machine learning models for octahedral transition metal complexes

Abstract Graph-based machine learning (ML) models for material properties show great potential to accelerate virtual high-throughput screening of large chemical spaces. However, in their simplest forms, graph-based models do not include any 3D information and are unable to distinguish stereoisomers such as those arising from different orderings of ligands around a metal center in coordination complexes. In this work we present a modification to revised autocorrelation descriptors, a molecular graph featurization method, for predicting spin state dependent properties of octahedral transition metal complexes (TMCs). Inspired by analytical semi-empirical models for TMCs, the new modeling strategy is based on the many-body expansion (MBE) and allows one to tune the captured stereoisomer information by changing the truncation order of the MBE. We present the necessary modifications to include this approach in two commonly used ML methods, kernel ridge regression and feed-forward neural networks. On a test set composed of all possible isomers of binary TMCs, the best MBE models achieve mean absolute errors (MAEs) of 2.75 kcal mol −1 on spin-splitting energies and 0.26 eV on frontier orbital energy gaps, a 30%–40% reduction in error compared to models based on our previous approach. We also observe improved generalization to previously unseen ligands where the best-performing models exhibit MAEs of 4.00 kcal mol −1 (i.e. a 0.73 kcal mol −1 reduction) on the spin-splitting energies and 0.53 eV (i.e. a 0.10 eV reduction) on the frontier orbital energy gaps. Because the new approach incorporates insights from electronic structure theory, such as ligand additivity relationships, these models exhibit systematic generalization from homoleptic to heteroleptic complexes, allowing for efficient screening of TMC search spaces.

Meyer, Ralf (ORCID:0000000322360261)↗
DOE OSTI2025

High-resolution spectroscopy of barium monofluoride: Odd isotopologues, hyperfine structure, and isotope shifts

Barium monofluoride (BaF) is a promising molecular species for precision tests of fundamental symmetries and interactions. We present a combined theoretical and experimental study of BaF spectra and isotope shifts, focusing in particular on the poorly understood odd isotopologues 137 BaF and 135 BaF. By comparing state-of-the-art ab initio calculations with high-resolution fluorescence and absorption spectroscopy data, we provide a benchmark for electronic structure theory and disentangle the hyperfine and rovibrational spectra of the five most abundant isotopologues, from 138 BaF to 134 BaF. The comprehensive knowledge gained enables a King plot analysis of the isotope shifts that reveals the odd-even staggering of the barium nuclear charge radii. It also paths the way for improved laser cooling of rare BaF isotopologues and crucially supports future measurements of nuclear anapole and Schiff moments.

Kogel, Felix [Univ. of Stuttgart (Germany)] (ORCID↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2025

The XPS of pyridine: A combined theoretical and experimental analysis

A detailed analysis of the N(1s) and C(1s) X-Ray Photoelectron Spectroscopy (XPS) is made, where the measured XPS is compared with theoretical Sudden Approximation (SA) intensities and theoretical XPS Binding Energies (BEs). There is remarkably good agreement between the theoretical predictions and the measured XPS; in particular, the different full width at half maximum values for the C(1s) and N(1s) BEs are explained in terms of unresolved C(1s) BEs for the different C atoms in pyridine. Here, this work demonstrates that the combination of theory and XPS measurements can extract analysis of the XPS relevant to the molecular electronic structure. The theory used is based on fully relativistic self-consistent field solutions of the Dirac–Coulomb Hamiltonian, and the SA is used to determine relative XPS intensities.

Binding Energies (BEs)↗
DOE OSTI2025

Effects of chemical disorder and spin-orbit coupling on the electronic-structure and Fermi-surface topology of YbSb-based monopnictides

Here, in this work, we study the influence of disorder on the electronic structure of YbSb—a rare-earth monopnictide featuring a simple rocksalt (B1) crystal structure and a well-defined Fermi-surface topology—by employing first-principles density-functional theory. We focus on chemical disorder introduced through Te and Al doping, selected based on their thermodynamic stability in alloyed configurations, to understand how such perturbations modify the electronic states of YbSb. Our results indicate that Te doping predominantly introduces electronlike states at the 𝑋 and 𝐿 points, while Al doping leads to a suppression of holelike states at Γ, effectively driving the system from a semimetallic state to one characterized by very narrow-gap behavior at Γ. This modulation of the Fermi surface, particularly the reduction of central hole pockets at Γ, plays a central role in altering interpocket scattering—a mechanism critical for tuning quantum transport properties, including superconductivity. This disorder-driven modulation of the Fermi surface, particularly the suppression of central hole pockets at Γ, controls interpocket scattering, which is essential for optimizing quantum transport properties, including superconductivity. Our results show that disorder can be effectively used for engineering band topology, thereby tuning quantum related response through a tailored electronic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2025

First principles free energy model with dynamic magnetism for δ -plutonium

We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.

dynamic↗
DOE OSTI2026

Structural Propensities in Cs2MBiX6 (M=Na, Ag; X=Cl, Br) Bismuth Halide Double Perovskites

A previously unreported low-temperature phase transition in the bismuth halide double perovskite Cs2AgBiCl6 is reported, thereby establishing trends in the structural ground state across Cs2NaBiCl6, Cs2AgBiCl6, and Cs2AgBiBr6. Using the combined toolkit of variable-temperature synchrotron X-ray and neutron powder diffraction, Raman spectroscopy, and density-functional theory–based electronic structure modeling, we demonstrate a cubic Fm¯3m → tetragonal I4/m transition upon cooling with distinct onset temperatures. Neutron powder diffraction refinements permit the unambiguously assignment of the low-temperature phase of Cs2NaBiCl6 to I4/m, correcting prior reports of an I4/mmm ground state. Cs2AgBiCl6 is also found to transforms to a structure crystallizing in the I4/m space group at low temperatures. Temperaturedependent Raman data and density-functional theory-based modeling capture the softening and freezing of out-of-phase octahedral-tilt modes and quantify relative instabilities. Solid-state nuclear magnetic resonance spectroscopy at room temperature completes the characterization and helps underpin the subtle differences in covalency across the compounds. Trends in the phase transition temperature Ts and tilt magnitudes emerge from coupled effects of halide identity, M(I)–site bonding character, and a mismatch between interatomic distances. These results establish the structure– dynamics–bonding framework for tuning tilt-driven instabilities in halide double perovskites.

Tian, Haowen↗
DOE OSTI2025

Interpretable machine learning-guided design of Fe-based soft magnetic alloys

Here, we present a machine learning (ML) guided approach to predict saturation magnetization (𝑀 S ) and coercivity (𝐻 C ) in Fe-rich soft magnetic alloys, particularly Fe-Si-B systems. ML models trained on experimental data reveal that increasing Si and B content reduces 𝑀 S from 1.81 T (DFT ≈ 2.04 T) to ≈1.54 T (DFT ≈ 1.56T) in Fe-Si-B, which is attributed to decreased magnetic density and structural modifications. Experimental validation of ML predicted magnetic saturation on Fe-1Si-1B (2.09 T), Fe-5Si-5B (2.01 T), and Fe-10Si-10B (1.54 T) alloy compositions further supports our findings. These trends are consistent with density functional theory predictions, which link increased electronic disorder and band broadening to lower 𝑀 S values. Experimental validation on selected alloys confirms the predictive accuracy of the ML model, with good agreement across compositions. Beyond predictive accuracy, detailed uncertainty quantification and model interpretability including through feature importance and partial dependence analysis reveal that 𝑀 S is governed by a nonlinear interplay between Fe content and early transition metal ratios, while 𝐻 C is more sensitive to processing conditions such as ribbon thickness and thermal treatment windows. The ML framework was further applied to Fe-Si-B/Cr/Cu/Zr/Nb alloys in a pseudoquaternary compositional space, which shows comparable magnetic properties to NANOMET (Fe 84.8 ⁢Si 0.5 ⁢B 9.4 ⁢Cu 0.8⁢ P 3.5 ⁢C 1 ), FINEMET (Fe 73.5 ⁢Si 13.5 ⁢B 9 Cu 1 ⁢Nb 3 ), NANOPERM (Fe 88 ⁢Zr 7⁢ B 4 ⁢Cu 1 ), and HITPERM (Fe 44 ⁢Co 44 ⁢Zr 7⁢ B 4 ⁢Cu 1 . Our findings demonstrate the potential of the ML framework for accelerated search of high-performance soft magnetic materials.

density functional theory↗
DOE OSTI2025

Relationship between atomic and electronic structure in Ln-bearing oxides

Metastable states of matter are of great interest as they offer the promise of novel functionality. They are often a natural consequence of exposure to nonequilibrium environments. In oxides, metastability can take the form of new polymorphs, chemical disorder, and even amorphization. While significant attention has been given to the impact those changes have on the atomic properties of the material, the corresponding changes in the electronic structure have received less attention. Here, using density functional theory, we consider how the electronic structure varies with potential metastable structures in two classes of lanthanide-bearing oxides—pyrochlores and interlanthanide sesquioxides. We find that the changes depend strongly on both the crystal structure and crystal chemistry of the compound with, for example, disordering and amorphization either increasing or decreasing the bandgap depending on the chemistry. For the 𝐴 2 ⁢𝐵 2 ⁢O 7 pyrochlores, we find different dependencies of the bandgap on the 𝐴 = 𝐿⁢𝑛 cations as the 𝐵 cation is changed, which we relate to the nature of the density of states at the conduction band minimum for different 𝐵 chemistries. Our calculations are validated by electron energy loss spectroscopy measurements for two pyrochlore compounds in which amorphization does reduce the bandgap, consistent with our calculations on these two compounds. In conclusion, our results highlight the relationship between atomic and electronic structure and how radiation can be used to modify and potentially control the electronic properties of oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nernst power factor and figure of merit in the compensated semimetal ScSb

Recently, topological semimetals have emerged as strong candidates for solid-state thermomagnetic refrigerators due to their enhanced Nernst effect. This enhancement arises from the combined contributions of the Berry-curvature-induced anomalous Nernst coefficient associated with topological bands and the normal Nernst effect resulting from synergistic electron-hole compensation. Generally, these two effects are intertwined in topological semimetals, making it challenging to evaluate them independently. Here, we report the observation of a high Nernst effect in the electron-hole compensated semimetal ScSb with topologically trivial electronic band structures. Remarkably, we find a high maximum Nernst power factor of 𝑃⁢𝐹 𝑁 ∼ 35 × 10 −4 W m −1 K −2 in ScSb. The Nernst thermopower (𝑆 𝑥⁢𝑦 ) exhibits a peak of ∼ 47 µ⁢V/K at 12 K and 14 T, yielding a Nernst figure of merit (𝑧 𝑁 ) of ∼ 28 × 10 −4 K −1 . Notably, despite its trivial electronic band structure, both the 𝑃⁢𝐹 𝑁 and 𝑧 𝑁 values of ScSb are comparable to those observed in topological semimetals with Dirac band dispersions. In conclusion, the origin of the large Nernst signal in ScSb is explained by compensated electron and hole carriers, through Hall resistivity measurements, angle-resolved photoemission spectroscopy, and density functional theory calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗