Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron localization function”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

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

Atomic ionization: sd energy imbalance and Perdew–Zunger self-interaction correction energy penalty in 3d atoms

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as E ( 3 d n − 1 4 s 1 ) − E ( 3 d n − 2 4 s 2 ) . Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the 4s 2 configuration. Here, we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from the local spin density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the regularized and restored Strongly Constrained and Appropriately Normed (r 2 SCAN) meta-GGA for first-row transition metal atoms. However, we find large (∼2 eV) ground-state sd energy imbalances when applying a Perdew–Zunger 1981 self-interaction correction. This is attributed to an “energy penalty” associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Science & Technology - Other Topics↗
DOE OSTI2026

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗
DOE OSTI2024

Localized Conduction Channels in Memristors

Since the early 2000s, the impending end of Moore’s scaling, as the physical limits to shrinking transistors have been approached, has fueled interest in improving the functionality and efficiency of integrated circuits by employing memristors or two-terminal resistive switches. Formation (or avoidance) of localized conducting channels in many memristors, often called “filaments”, has been established as the basis for their operation. While we understand some qualitative aspects of the physical and thermodynamic origins of conduction localization, there are not yet quantitative models that allow us to predict when they will form or how large they will be. Here, in this work, we compile observations and explanations of channel formation that have appeared in the literature since the 1930s, show how many of these seemingly unrelated pieces fit together, and outline what is needed to complete the puzzle. This understanding will be a necessary predictive component for the design and fabrication of post-Moore’s-era electronics.

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

Effects of iron carbide crystal phases and dopants on the conversions of CO 2 into ethylene

The density functional theory method was used to investigate the conversions of CO 2 to ethylene formation on two common iron carbide surfaces: Fe 3 C(0 1 0) and Fe 5 C 2 (1 1 1). Based on the structure relaxation of reaction intermediates and the elementary reaction transition states. We deduced the most competitive reaction pathways for ethylene production. The main CO 2 -to-ethylene routes and the competition of side products, CO and CH 4 , are discussed. Our analyses showed that CO 2 conversion is surface structure sensitive, whereas CH 4 and C2+ hydrocarbon formations depend on the reactivity of native C atoms in the carbides. To modify the intrinsic catalyst performance, mixing dopants in Fe catalysts is an effective strategy. Furthermore, we demonstrate that doping Zn and Zr can alter the local electronic structure and enhance CO 2 adsorption on the catalyst surface.

CO2 hydrogenation↗
DOE OSTI2026

How Does Metal Spin State Affect Electronic Communication in Mixed-Valence Dimers? Insights from Ultrafast Near-Infrared and Soft X-ray Transient Absorption Spectroscopy

Recent advancements in photocatalysis, photovoltaics, and quantum information science take advantage of electron spin, and determining how spin multiplicity affects electron transfer is key to understanding these phenomena. Here, in this study, we examine how metal spin state affects electronic communication in an organometallic mixed-valence dimer, ferrocenyl cobaltocenium ([Fe II Cp 2 Co III Cp 2 ] + ). This complex can be photoexcited from its low-spin singlet Fe II ground state to form intermediate-spin triplet Fe II and high-spin quintet Fe II excited states. Using femtosecond optical transient absorption (OTA) spectroscopy with visible (vis), near-infrared (NIR), and short-wave IR (SWIR) probes, supported by time-dependent density functional theory (TD-DFT) calculations, we measure Fe II Co III →Fe III Co II intervalence charge transfer (IVCT) bands in each of the Fe II spin states. Mulliken–Hush analysis of the excited-state IVCT bands was used to compute the electronic coupling between the metal centers in all three spin states, which increased as quintet < triplet < singlet. Meanwhile, the peak energy of the bands, and thus the ΔG of the IVCT transition, increased as triplet < quintet < singlet. Then, to directly probe the electronic structure at each metal center, we employed picosecond soft X-ray transient absorption (XTA) spectroscopy at the Fe and Co L 3 edges. Our results show that the low-spin and high-spin states of [Fe II Cp 2 Co III Cp 2 ] + are valence-localized, while the intermediate-spin state is partially delocalized. The differences in charge delocalization are attributed to differences in orbital occupation and geometry that affect the free energy and superexchange coupling.

Burke, John H. [Univ. of Illinois at Urbana-Champa↗
DOE OSTI2026

Transfer Learning Meets Embedded Correlated Wavefunction Theory for Chemically Accurate Molecular Simulations: Application to Calcium Carbonate Ion Pairing

Achieving chemical accuracy for molecular simulations remains a central challenge in computational chemistry. Here, we present an embedded correlated wavefunction transfer learning (ECW-TL) framework for accurately simulating molecular dynamics in the condensed phase. ECW-TL incorporates high-level electron exchange and correlation effects in ECW theory while preserving the training and computational efficiency of machine-learned interatomic potentials. We demonstrate the framework on Ca 2+ –CO 3 2– ion pairing in aqueous solution, a key process underlying CO 2 mineralization in seawater. As proof of principle, we first show that fine-tuning a DFT-revPBE-D3(BJ) baseline model with embedded-DFT-SCAN data reproduces the DFT-SCAN free-energy surface within 1 kcal/mol across all solvation states. Extending the framework to embedded MP2 and localized natural-orbital CCSD(T) further refines the free-energy profile, revealing the crucial role of exact electron exchange and correlation in determining ion-pair stability and structure. The computed ion-pair association free energy is in quantitative agreement with experimental measurements, further validating the accuracy of the ECW-TL framework. ECW-TL thus provides a general, data-efficient route for transferring CW accuracy to efficient simulations of complex aqueous and interfacial chemical processes.

cluster chemistry↗
DOE OSTI2026

Observation of Iso-Symmetric Structural and Lifshitz Transitions in Quasi-One-Dimensional CrNbSe 5

Chalcogenides-rich transition metal compounds host a rich landscape of emergent quantum phenomena that are intimately governed by their quasi-onedimensional chemical-bonding frameworks and their response to external perturbations such as pressure. Here, we report a pressure-induced iso-symmetric structural transition in the quasi-one-dimensional compound CrNbSe 5 , in which the electronic ground state is controlled not by symmetry breaking but by a continuous reorganization of local bonding interactions. Applied pressure reversibly tunes CrNbSe 5 between semiconducting and semimetallic states, enabling access to low- and high-carrier electronic regimes through direct modulation of metal−chalcogen bonding. High-pressure singlecrystal X-ray diffraction directly resolves the evolution of Cr−Se and Nb−Se bond distances, coordination polyhedra, and connectivity, revealing a fully reversible semimetal−semiconductor−semimetal transition driven by gradual yet cooperative bond rearrangements within a preserved crystallographic symmetry. In contrast to chemical substitution, which irreversibly alters composition and introduces disorder, pressure acts as a clean, continuous control parameter that reshapes the bonding landscape without disrupting structural symmetry. These results establish CrNbSe 5 as a model system for electronically driven phase switching via tunable chemical bonding, highlighting iso-symmetric bond reorganization as a powerful design principle for pressure-controlled electronic and spintronic functionalities.

Compression↗
DOE OSTI2025

Long-Range Allosteric Communication Modulated by Active Site Mn(II) Coordination Drives Catalysis in Xanthobacter autotrophicus Acetone Carboxylase

Acetone carboxylase (AC) from Xanthobacter autotrophicus is a 360 KDa α2β2γ2 heterohexamer that catalyzes the ATP-dependent formation of phosphorylated acetone and bicarbonate intermediates that react at Mn(II) metal active sites to form acetoacetate. Structural models of X. autotrophicus AC (XaAC) with and without nucleotides reveal that the binding and phosphorylation of the two substrates occurs ~40 Å from the Mn(II) active sites where acetoacetate is formed. Based on the crystal structures, a significant conformational change was proposed to open and close a tunnel that facilitates the passage of reaction intermediates between the sites for nucleotide binding and phosphorylation of substrates and Mn(II) sites of acetoacetate formation. We have employed electron paramagnetic resonance (EPR), kinetic assays, and hydrogen/deuterium exchange mass spectrometry (HDX-MS) of poised ligand-bound states and site-specific amino acid variants to complete an in-depth analysis of Mn(II) coordination and allosteric communication throughout the catalytic cycle. In contrast with the established paradigms for carboxylation, our analyses of XaAC suggested a carboxylate shift that couples both local and long-range structural transitions. Shifts in the coordination mode of a single carboxylic acid residue (αE89) mediate both catalysis proximal to a Mn(II) center and communication with an ATP active site in a separate subunit of a 180 kDa α2β2γ2 complex at a distance of 40 Å. This work demonstrates the power of combining structural models from X-ray crystallography with solution-phase spectroscopy and biophysical techniques to elucidate functional aspects of a multi-subunit enzyme.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Evolution of Structure and Magnetism in FeCl 2 and FeCl 3 : From Clusters to Monolayers

In this work, we address one of the most fundamental questions in cluster science─how do the structure and properties evolve from clusters to crystals? Using density functional theory (DFT), we focus our study on the evolution of structure and magnetism in iron-chloride systems, from clusters to monolayers. The choice of this system is motivated by the recent experimental confirmation of one of the author’s earlier theoretical prediction that the FeCl 2 cluster is magnetic with a spin magnetic moment of 4 μ B localized at the Fe site, while its dimer, Fe 2 Cl 4 , is antiferromagnetic. Similarly, FeCl 3 cluster is magnetic with a total spin magnetic moment of 5 μ B , with 4 μ B localized at the Fe site and 1 μ B distributed over the Cl sites. The dimer clusters Fe 2 Cl 4 and Fe 2 Cl 6 have an antiferromagnetic ground state, and upon Li-functionalization, both can be magnetically transformed from antiferromagnetic to ferromagnetic states. In contrast, FeCl 2 and FeCl 3 monolayers exhibit different magnetic ground states in their periodic forms: FeCl 2 is ferromagnetic (FM), but in FeCl 3 , the antiferromagnetic (AFM) and FM states are energetically nearly degenerate. Such a difference arises due to the different chemical coordination of the Fe atoms with the Cl atoms, caused by their different oxidation states, which is +2 in FeCl 2 and +3 in FeCl 3 , respectively. Interestingly, Li-functionalization allows both FeCl 3 and FeCl 2 monolayers to be ferromagnetic. Our study highlights that several, but not all, electronic and magnetic characteristics of isolated clusters are preserved in the extended periodic structures. This systematic investigation of iron-halide clusters is expected to inspire further experimental and theoretical exploration into the magnetism of other transition metal halides.

cluster chemistry↗
DOE OSTI2026

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗
DOE OSTI2025

Multiple Slater Determinants and Strong Spin-Fluctuations as Key Ingredients of the Electronic Structure of Electron- and Hole-Doped Pb 10− x Cu x (PO 4 ) 6 O

LK-99, with chemical formula Pb 10−x Cu x (PO 4 ) 6 O, was recently reported to be a room-temperature superconductor. While this claim has met with little support in a flurry of ensuing work, a variety of calculations (mostly based on density-functional theory) have demonstrated that the system possesses some unusual characteristics in the electronic structure, in particular flat bands. We have established previously that within DFT, the system is insulating with many characteristics resembling the classic cuprates, provided the structure is not constrained to the P3(143) symmetry nominally assigned to it. Here we describe the basic electronic structure of LK-99 within self-consistent many-body perturbative approach, quasiparticle self-consistent GW (QSGW) approximation and their diagrammatic extensions. QSGW predicts that pristine LK-99 is indeed a Mott/charge transfer insulator, with a bandgap gap in excess of 3 eV, whether or not constrained to the P3(143) symmetry. When Pb 9 Cu(PO 4 ) 6 O is hole-doped, the valence bands modify only slightly, and a hole pocket appears. However, two solutions emerge: a high-moment solution with the Cu local moment aligned parallel to neighbors, and a low-moment solution with Cu aligned antiparallel to its environment. In the electron-doped case the conduction band structure changes significantly: states of mostly Pb character merge with the formerly dispersionless Cu d state, and high-spin and low spin solutions once again appear. Thus we conclude that with suitable doping, the ground state of the system is not adequately described by a band picture, and that strong correlations are likely. Irrespective of whether this system class hosts superconductivity or not, the transition of Pb 10 (PO 4 ) 6 O from being a band insulator to Pb 9 Cu(PO 4 ) 6 O, a Mott insulator, and multi-determinantal nature of doped Mott physics make this an extremely interesting case-study for strongly correlated many-body physics.

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

Tip-Induced 3D Printing on the Nanoscale with Field Emission Scanning Probes

3D printing down to the nanoscale remains a significant challenge. In this paper, the study explores the use of scanning probes that emit low-energy electrons (<100 eV) coupled with the localized injection and electron-induced decomposition of precursor molecules, for the precise localized deposition of 3D nanostructures. The experiments are performed inside the chamber of a scanning electron microscope (SEM), enabling the use of the in-built gas injector system (GIS) with gaseous naphthalene precursor for carbon deposition, as well as immediate inspection of the deposits by SEM. Substrate materials are planar fused silica with thin conductive coatings and non-planar copper wedges. After investigation of the deposition process parameters, various 2D and 3D carbon deposits are grown. Vertical nanowires several microns in length with a diameter <100 nm are achieved and 3D deposits with a high degree of nanoscale branching are also obtained, presumably due to a charging effect. High aspect ratio carbon nanostructures such as those demonstrated here can be employed as miniaturized electrodes or field emitters. The tip-based approach presented thus paves the way toward 3D nanoscale printing of various materials and functional nanostructures.

3D printing↗
DOE OSTI2025

Competing conduction mechanisms in high performance carbon nanotube fibers

The performance of carbon nanotube (CNT) cables, a contender for copper-wire replacement, is tied to its metallic and semi-conducting-like conductivity responses with temperature; the origin of the semi-conducting-like response however is an underappreciated incongruity in literature. With controlled aspect-ratio and doping-degree, over 61 unique cryogenic experiments including anisotropy and Hall measurements, CNT cable performance is explored at extreme temperatures (65 mK) and magnetic-fields (60 T). A semi-conducting-like conductivity response with temperature becomes temperature-independent approaching absolute-zero, uniquely demonstrating for the first time the necessity of heterogeneous fluctuation induced tunneling; complete de-doping leads to localized hopping, contrasting graphite's pure metallic-like response. High-field magneto-resistance (including novel +22 % longitudinal magneto-resistance near room-temperature) is analyzed with hopping and classical two-band models, both yielding a similar parameter useful for conductor development. Varying field-orientation angle uncovers significant two- and four-fold symmetries that are shown to be from Aharonov-Bohm-like corrections to the curvature-induced bandgap, a first for macroscale CNT fibers. Tight-binding calculations using Green's Function formalism model the largest, coherent transport to-date in commensurate CNT bundles in magnetic-field, revealing non-uniform transmission across bundle cross-sections with doping restoring uniformity; independent of doping, transport in bundle-junction-bundle systems are predominantly from CNTs adjacent to the other bundle—demonstrating that smaller bundles are more efficient for electronic transport. The final impact is predicting the ultimate conductivity of heterogeneous CNT cables using temperature and field-dependent transport, surpassing conductivity of traditional metals.

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

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗