Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vapor–liquid–solid”

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.

273 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Progress and Perspectives: Zirconium Electrodeposition from Different Electrolytes

Zirconium (Zr) possesses outstanding properties, including exceptional chemical resistance and a high melting point, and is therefore desirable for use in a wide variety of challenging environments, such as in nuclear reactors and the chemical processing industry. To minimize the amount of Zr required for any given application, developing methods for generating metallic Zr coatings is highly advantageous. Owing to Zr’s highly negative reduction potential, electrodeposition in traditional solvents at near ambient temperatures remains challenging and thus not well understood. Due to the extreme conditions required to deposit this metal, extensive work has been conducted in molten salt electrolytes, however the broad applicability of this methodology is limited due to its corrosivity. The primary focus of this article is to present an overview of Zr’s electrochemical behavior and to consolidate the efforts of researchers in exploring electrodeposition techniques for Zr involving aqueous, organic, ionic liquid, deep eutectic, and molten salt solvents. With this information, we highlight trends across solvent systems and opportunities for future research.

42 ENGINEERING↗
DOE OSTI2026

Unraveling interphase-driven failure pathways in LiMn0.6Fe0.4PO4/graphite pouch cells

LiMnxFe1−xPO4 (LMFP) is a promising high-voltage, thermally stable, and earth-abundant cathode material, yet its practical application is limited by interphase instability and Mn dissolution. In this work, we systematically evaluate LiMn0.6Fe0.4PO4/graphite pouch cells using three electrolyte formulations including control carbonate electrolyte, control + 2 wt% vinylene carbonate (VC), and control + 2 wt% VC + 1 wt% 1,3,2-dioxathiolane 2,2-dioxide (DTD), to establish how electrolyte composition governs interphase chemistry and long-term degradation. Electrochemical testing shows that both additives are preferentially reduced prior to ethylene carbonate (EC) during cell formation, generating robust cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI) layers that suppress gas evolution and raise the first-cycle coulombic efficiency to 89.3%. Additionally, the dual-additive electrolyte delivers the most stable performance, retaining over 85% capacity after 600 cycles while minimizing impedance growth under long-term cycling at C/3 and 40 °C. Soft X-ray absorption spectroscopy confirms that VC + DTD effectively suppresses electrolyte oxidation at the cathode surface, and micro-X-ray fluorescence shows substantially reduced Mn dissolution and deposition on the graphite anode. Density functional theory simulations further provided insights into the structural and energetic influences of alkoxide species on the cathode surface, proposing a Mn2+ extraction mechanism. The combined experimental and computational findings establish a mechanistic link between electrolyte composition and interphase evolution, highlighting the effectiveness of electrolyte engineering for extending the operational lifetime of LMFP-based lithium-ion batteries.

Chak, Chanmonirath Michael↗
DOE OSTI2026

Catalytic disproportionation on carbon superstructures enables long-life, high-loading Li–S batteries

Electrocatalysis has been widely explored as an effective strategy to accelerate polysulfide (PS) conversion and suppress the shuttle effect in lithium–sulfur (Li–S) batteries. However, the underlying mechanisms remain elusive, and electrocatalytic reactions are inactive during cell resting. In this work, we reveal and quantitatively analyze a previously unrecognized sulfur reduction route (SRR) driven by catalytic disproportionation at the carbon cathode surface—fundamentally distinct from conventional electrocatalysis. Unlike conventional stepwise pathways, this SRR enables high-order polysulfides (Sₓ²⁻, x = 5–8) to directly convert into S₈ and Li₂S₂, bypassing low-order intermediates. This sulfur-reduction shortcut is systematically elucidated through high-performance liquid chromatography, revealing the intrinsic catalytic contribution of carbon frameworks and the dynamic evolution of PS species. We demonstrate that carbon superstructures (CSS-0.5), assembled from nanosheet subunits with abundant N/O functionalities and interconnected charge-migration channels, synergistically promote this catalytic process. Benefiting from these features, CSS-0.5 delivers superior electrochemical performance under practical conditions, enabling high sulfur loading (6.0 mg cm⁻²) pouch cells with 80.5% capacity retention over 210 cycles. This study provides the first quantitative evidence of electrocatalytic disproportionation in Li–S batteries, offering mechanistic insights and design principles for advanced sulfur cathodes.

25 ENERGY STORAGE↗
DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2023

Dynamic Oscillation and Motion of Oil‐in‐Water Emulsion Droplets

The interplay of interfacial tensions on droplets results in a range of self‐powered motions that mimic those of living systems and serve as a tunable model to understand their complex non‐equilibrium behavior. Spontaneous shape deformations and oscillations are crucial features observed in nature but difficult to incorporate in synthetic artificial systems. Here, we report sessile oil‐in‐water emulsions that exhibit rapid oscillating behavior. The oscillations depend on the nature and concentration of the surfactant, the chemical composition of the oil, and the wettability of the solid substrate. The rapid changes in the contact angle per oscillation are observed using side‐view optical microscopy. We propose that the changes in the interfacial tension of the oil droplets is due to the partitioning of the surfactant into the oil phase and the movement of self‐emulsified oil out of the parent droplets giving rise to the rhythmic variation in droplet contact‐line. The ability to control and understand droplet oscillation can help model similar oscillations in out‐of‐equilibrium systems in nature and reproduce biomimetic behavior in artificial systems for various applications, such as microfluidic lab‐on‐a‐chip and adaptive materials.

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

Adaptive Scalpel Scanning Probe Microscopy for Enhanced Volumetric Sensing in Tomographic Analysis

Controlling nanoscale tip‐induced material removal is crucial for achieving atomic‐level precision in tomographic sensing with atomic force microscopy (AFM). While advances have enabled volumetric probing of conductive features with nanometer accuracy in solid‐state devices, materials, and photovoltaics, limitations in spatial resolution and volumetric sensitivity persist. This work identifies and addresses in‐plane and vertical tip‐sample junction leakage as sources of parasitic contrast in tomographic AFM, hindering real‐space 3D reconstructions. Novel strategies are proposed to overcome these limitations. First, the contrast mechanisms analyzing nanosized conductive features are explored when confining current collection purely to in‐plane transport, thus allowing reconstruction with a reduction in the overestimation of the lateral dimensions. Furthermore, an adaptive tip‐sample biasing scheme is demonstrated for the mitigation of a class of artefacts induced by the high electric field inside the thin oxide when volumetrically reduced. This significantly enhances vertical sensitivity by approaching the intrinsic limits set by quantum tunneling processes, allowing detailed depth analysis in thin dielectrics. The effectiveness of these methods is showcased in tomographic reconstructions of conductive filaments in valence change memory, highlighting the potential for application in nanoelectronics devices and bulk materials and unlocking new limits for tomographic AFM.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2026

Mitigation of polysulfide shuttle effect in Li-S batteries through catalytic disproportionation reaction

Polysulfides are poorly retained within porous cathodes and readily diffuse into the electrolyte over time, leading to the well-known shuttle effect that undermines the reversibility of Li-S batteries. Here, in this study, we demonstrate that catalytic disproportionation of polysulfides provides an effective pathway to suppress this process by rapidly converting dissolved species into solid sulfur and sulfides, thereby preventing their migration into the electrolyte. Fundamentally, the sluggish kinetics of sulfur redox reactions are responsible for the accumulation and redistribution of soluble polysulfides in the bulk electrolyte. By accelerating these kinetics, catalyzed disproportionation not only confines sulfur within the conductive cathode matrix but also promotes the homogeneous precipitation of Li₂S₂/Li₂S, which enhances electrochemical reversibility and cycling stability. Using nitrogen-doped carbon (NC800) as a model catalyst, we reveal its ability to drive a pseudo-16-electron reduction pathway, leading to a single dominant Li₂S product and uniform deposition within the porous framework. In contrast, a non-catalytic carbon (KB) yields multiple polysulfide intermediates and heterogeneous deposition. The mechanistic insights provided here highlight the pivotal role of catalytic disproportionation in reshaping sulfur redox pathways and offer a rational strategy for mitigating polysulfide shuttling in practical Li-S pouch cells.

25 ENERGY STORAGE↗
DOE OSTI2026

Jumpstart Opportunities to Unleash Leadership in Energy Storage (JOULES)

Current-generation Li-ion batteries with cobalt- and nickel-containing cathodes and graphite anodes are approaching performance and cost limits. In this program, 24M Technologies, Inc. (24M) is teaming with the Massachusetts Institute of Technology (MIT) and University of Michigan (UM) to develop low cost and fast charging sodium metal batteries with good low-temperature performance and high energy density, building upon previous work performed under ARPA-E programs. Key achievements include optimization of solid electrolyte and anode current collector, optimized cathode active materials, development of high-performance electrolyte formulations, and integration of these components into full cells. The cell design incorporates (1) an ultra-thick cathode (>9 mAh/cm 2 ) comprising advanced cobalt-free, sodium cathode active material, (2) advanced fast-charging electrolyte (up to 12 mS/cm) developed using machine learning and automated high-throughput screening technology by UM, and (3) ceramic modified separator that enable smooth Na transport and deposition, developed at MIT, enabling a high-energy density anode-free configuration and maximizing the energy density of sodium batteries. The team has successfully combined these approaches to sodium chemistry and paved the way to meeting the fast-charging, high-energy density, and low-cost requirements of next-generation drone, electric vertical take-off and -landing, and electric vehicle batteries. Performance for anode-free sodium cells developed under this program is more powerful than the commercial Li-ion batteries. The final deliverable cell design has achieved over 300 Wh/kg and volumetric energy density above 800 Wh/L (Table 1). Additionally, the team has achieved over (1) a lifetime of 340 cycles, (2) 80% capacity retention at -20 °C (compared 25 °C), and (3) the ability to fast charge to 80% SOC in 20 minutes.

25 ENERGY STORAGE↗
DOE OSTI2023

A Self‐Healing System for Polydicyclopentadiene Thermosets

Self-healing offers promise for addressing structural failures, increasing lifespan, and improving durability in polymeric materials. Implementing self-healing in thermoset polymers faces significant manufacturing challenges, especially due to the elevated temperature requirements of thermoset processing. To introduce self-healing into structural thermosets, the self-healing system must be thermally stable and compatible with the thermoset chemistry. This article demonstrates a self-healing microcapsule-based system stable to frontal polymerization (FP), a rapid and energy-efficient manufacturing process with a self-propagating exothermic reaction (≈200 °C). A thermally latent Grubbs-type complex bearing two N-heterocyclic carbene ligands addresses limitations in conventional G2-based self-healing approaches. Under FP's elevated temperatures, the catalyst remains dormant until activated by a Cu(I) co-reagent, ensuring efficient polymerization of the dicyclopentadiene (DCPD) upon damage to the polyDCPD matrix. The two-part microcapsule system consists of one capsule containing the thermally latent Grubbs-type catalyst dissolved in the solvent, and another capsule containing a Cu(I) coagent blended with liquid DCPD monomer. Using the same chemistry for both matrix fabrication and healing results in strong interfaces as demonstrated by lap-shear tests. In an optimized system, the self-healing system restores the mechanical properties of the tough polyDCPD thermoset. Self-healing efficiencies greater than 90% via tapered double cantilever beam tests are observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The impact of metastability on the high-pressure behavior of cerium

The structures adopted by solids under pressure are often assumed to reflect thermodynamic equilibrium, yet in many materials phase selection is strongly influenced by kinetic pathways and microstructural inheritance. Elemental cerium (Ce) exemplifies this challenge, with decades of conflicting reports describing different high-pressure crystal structures emerging under nominally identical conditions. Here we use neutron diffraction from large ( ~ 60 mm 3 ) sample volumes to follow the structural evolution in ultra-high-purity Ce during controlled pressure-temperature cycling between 85 and 295 K and up to 8 GPa. We find that the crystal structure formed at high pressure depends on the compression pathway: slow compression ( ~ 0.25 GPa hr −1 ) at room temperature favors an orthorhombic phase (α'), whereas slow ( ~ 0.25 GPa hr −1 ) and also moderately faster ( ~ 0.5 GPa hr −1 ) compression at low temperature stabilizes a pure monoclinic phase (α"). The low-temperature phase persists metastably over a wide temperature range but transforms irreversibly upon heating above ~ 280 K, or modest pressure cycling. Remarkably, the lower-pressure γ phase remains trapped far beyond the expected stability field, persisting to the highest pressures studied. These observations show that phase selection in Ce is governed by kinetics and microstructural memory rather than equilibrium thermodynamics or sample morphology alone, establishing path dependence as a defining feature of its high-pressure behavior.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2026

Modern insights into the mechanisms of neptunium oxalate decomposition

Neptunium oxalate (Np(C 2 O 4 ) 2 ·6H 2 O) is a historically relevant solid phase used in nuclear processing as a precursor for neptunium dioxide (NpO 2 ). Although Np oxalate has been synthesized and used for NpO 2 production for decades, the thermal decomposition mechanism of this phase remains poorly understood and has not been evaluated in over 30 years. Conflicting reports in historical literature suggest either a direct conversion from anhydrous oxalate to NpO 2 or a decomposition that includes the formation of Np carbonate or oxidized Np intermediate phases. In this work, we reexamine the decomposition pathway of Np(C 2 O 4 ) 2 ·6H 2 O using thermal analysis coupled with evolved gas analysis and temperature-dependent Raman spectroscopy to elucidate decomposition mechanisms and intermediate phases using modern analytical techniques. Thermal analysis revealed a three-stage decomposition process, including dehydration below 200 °C, oxalate breakdown between 170 and 370 °C, and NpO 2 formation by 500 °C. However, an unidentified plateau in the thermal data was observed during measurements. Raman spectroscopy confirmed the stages of decomposition, and in the analysis of potential intermediate phases, no carbonate phases or Np 2 O 5 were identified. Raman data suggest that residual oxalate or nonstoichiometric oxide are present during decomposition before pure NpO 2 is formed. These findings clarify aspects of the Np oxalate decomposition mechanism and address longstanding discrepancies in the literature, with a specific focus on Np-specific materials chemistry.

Lawson, Kathryn M. [Oak Ridge National Laboratory ↗
DOE OSTI2026

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)↗
DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
DOE OSTI2026

Self-aligned heterogeneous quantum photonic integration

Integrated quantum photonics holds significant promise for scalable photonic quantum information processing, quantum repeaters, and quantum networks, but its development is hindered by the mismatch between materials hosting high-quality quantum emitters and those compatible with mature photonic technologies. Heterogeneous integration offers a potential solution to this challenge, yet practical implementations have been limited by inevitable insertion losses at material interfaces. Here, we present a self-aligned heterogeneous quantum photonic integration approach that enables near-unity coupling efficiency at the interface. To showcase our approach, we demonstrate Purcell enhancement of a silicon vacancy (SiV) center in diamond induced by a heterogeneous photonic crystal cavity defined by titanium dioxide (TiO 2 ), as well as optical spin control and readout via a TiO 2 photonic circuit. We further show that, when combined with inverse photonic design, our approach enables efficient and broadband collection of single photons from a color center into a heterogeneous waveguide. Our approach is not restricted to SiV centers or TiO 2 ; it has the potential to be broadly applied to integrate diverse solid-state quantum emitters with thin-film photonic devices where conformal deposition is possible. Together, these results establish a practical route to scalable quantum photonic integrated circuits that combine high-quality quantum emitters with technologically mature photonic platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Maximizing dynamic range and performance of anatase TiO 2 ECRAM through structure and programming

Here, in this study, we investigate the structure-dependent modulation characteristics of all-solid-state three-terminal electrochemical random-access memory (ECRAM) based on an anatase Li x TiO 2 channel. By directly comparing “asymmetric” and “symmetric” ECRAM device architectures, we reveal significant insight into the impact of a non-zero gate-drain open-circuit voltage and its influence on voltage vs. current-controlled gating. We also explore the impact of potentiation/depression write parameters on the symmetry, linearity, and dynamic range of the device response. Together, initial results from optimizing structure and programming approaches yielded unprecedented G max /G min ratios of >1,000 for ECRAM and hundreds of tunable memory states with excellent linearity and symmetry. Simulations based on these ECRAM devices further illustrate the promise of this analog memory technology, achieving near 2% classification error in the MNIST digit recognition benchmark for a range of training parameters compared to a theoretical best of 1.66% and outperforming other device models extracted from the literature.

AIHWKit↗
DOE OSTI2026

Analysis of an irradiated uranium sample for source attribution without chemical separation using microplasma ionization and ultrahigh resolution mass spectrometry

The use of element isotope ratios has great potential in not only determining the reactor type used to produce plutonium (Pu) but also in determining the burnup and the time since irradiation. While a powerful nuclear forensic technique, determining element isotope ratios is complicated by severe isobaric interferences when performed on typical inductively coupled plasma mass spectrometers. Such analyses require extensive chemical separations prior to analysis to alleviate the inter-elemental isobars. Ultrahigh mass resolution spectrometry provides a potential alternative, greatly reducing the complexity of sample preparation and turnaround times for these critical measurements. To demonstrate the power of the approach, a sample of irradiated, depleted uranium was analyzed with the liquid sampling—atmospheric pressure glow discharge ion source coupled to an Orbitrap mass spectrometer. The Orbitrap is augmented with an external data acquisition system, Spectroswiss’s FTMS-Booster X2T, allowing collection of extended ion transients, providing higher mass resolution. In using this approach, the 150 Sm/ 149 Sm and 152 Sm/ 149 Sm isotope ratios were found to be within 20% of predicted values without any chemical separations and without mass bias corrections. In addition, the 240 Pu/ 239 Pu isotope ratio was determined, free from the 238 UH + interferences common to the ICP-MS platforms, while at the same time allowing for the determination of U isotopic signatures. While these demonstrative results are from a single sample, the advantages of the microplasma/ultrahigh mass resolution approach to intra-element isotope ratio determinations are clear.

Fuel burnup↗
DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗