Search NASA⌕ Search

SEARCH · Search NASA

Results for “Langley Working Papers”

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.

510 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Proteome-wide analysis of protein stability in Escherichia coli under acid stress

Knowledge of protein acid sensitivity remains sparse and is largely derived from low-throughput, enzyme-specific assays. We used a scalable framework to map acid stability across the Escherichia coli proteome to assess the acid stability of 1,675 unique proteins, estimating pH 50 values for over 90% of them. The parameter pH50 was defined as the pH value at which only 50% of the initial protein remains in solution following acid treatment. Proteome-wide pH 50 values ranged from 2.28 to 6.33 (median 5.11). Approximately 9% of detected proteins remained stable across all tested pH conditions. Our results align with published data and the assay of citrate synthase (GltA) performed here. Protein acid stability differed significantly by subcellular localization: periplasmic proteins were relatively more abundant in the acid-stable group, cytoplasmic proteins were abundant at pH 50 values 4.5–5.5, and inner membrane proteins at higher pH 50 between 5.5 and 6.0. Outer membrane proteins were too few to draw strong conclusions regarding enrichment within specific pH 50 groups. Notably, the periplasmic binding protein of the molybdate ABC transporter (ModA), was enriched after incubation at low pH. Estimated pH 50 values showed no correlation with protein isoelectric point and molecular weight. Together, this work provides the first proteome-wide map of protein acid stability and establishes a general framework for studying different chemical stressors.

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

Optical Vibrational Spectroscopic Investigation of Natural and Synthetic Analogs of Uranyl Oxyhydroxyhydrate Minerals

Uranyl oxy-hydroxy-hydrate minerals are common alteration products of uraninite (UO2+x) which is chemically and structurally analogous to uranium dioxide nuclear fuel. Therefore, structural and spectroscopic investigations of these alteration minerals and their analogous anthropogenic counterparts can provide insight into the environmental behavior of nuclear fuel cycle materials. Previously, we compiled available vibrational spectroscopic data for uranyl minerals in the Compendium of Uranium Raman and Infrared Experimental Spectra (CURIES) and found that only 37% of known uranyl oxy-hydroxy-hydrate minerals had spectra readily available in the literature and existing databases for inclusion therein. Furthermore, no available infrared spectra for this mineral group were included in CURIES. To expand our understanding of the spectroscopic features of uranyl oxy-hydroxy-hydrates and the structural origins thereof, we collected, and now include in CURIES, Raman and infrared spectra for an additional 12 uranyl hydroxide phases. To better understand the impact of structural and compositional variations of these phases on their spectroscopic features, we compare Raman spectra of different anion sheet topological groups and of analogous phases hosting different counter cations. We identify spectroscopic variations related to differences in equatorial bonding and structural changes as a result of cation substitution. We also prepare a uranyl hydroxide phase via hydrolysis of uranyl fluoride (UO2F2) as an analog of hydrolysis reactions that occur in nuclear fuel cycle materials; and we find that the alteration product of UO2F2, despite chemical and structural similarities to uranyl oxy-hydroxy-hydrate minerals, is readily distinguishable from related mineral phases using Raman spectroscopy. In this work, we provide new insights into the structural origins of spectroscopic features in uranyl oxy-hydroxy-hydrate minerals, improve the average Raman spectrum for this group of minerals, and thereby improve capabilities for identifying these mineral species and related anthropogenic phases using Raman spectroscopy.

Barth, Brodie [ORNL] (ORCID:0000000256142601)↗
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 OSTI2026

Scalable multiplexed machine learning gas sensor chips for food classification

Multiplexed gas sensor arrays combined with machine learning have unlocked previously inaccessible applications for scent-based sensing. Current platforms are limited by overlapping sensing materials with similar compositions, leading to highly correlated responses, or multistep deposition processes that hinder scalability. In this work, we developed a 16-element monolithic chip with fully distinct sensing layers, enabling a truly heterogeneous array. The system consists of highly sensitive carbon nanotube field effect transistors that are functionalized through a single-step microdispensing method compatible with automated pipetting systems. The resulting chip produces characteristic signal patterns in response to object-specific scent profiles and, when combined with machine learning algorithms, can perform automated object identification. We demonstrate the classification of 16 different objects, including food spoilage and nut allergens, with a 92.6% overall prediction accuracy.

Bassil, Carla [University of California, Berkeley,↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Structure of Self-Generated Magnetic Fields in Laser-Solid Interaction from Proton Tomography

Self-generated magnetic fields in laser-solid interactions are experimentally characterized to reveal the 3D location and local field strength, rather than path-integrated quantities, using multi-view proton radiography and tomographic inversion. We infer magnetic fields that extend several millimeters off the target into the hot, rarefied corona, sufficient to strongly magnetize the plasma (Ω e τ e ≫ 1). The data are compared to MHD simulations incorporating recent improvements in modeling magnetic field generation and transport; the volume-averaged coronal field strength and magnetic flux agree to within 25% using a model with magnetic re-localization of transport, although the near-target morphology is not reproduced. This work demonstrates tomographic proton radiography as a valuable tool for investigating magnetic fields in laser-produced plasmas.

High-energy-density plasmas↗
DOE OSTI2025

Surface deposition of nanometer scale carbon structures on Bi 2 Sr 2 CaCu 2 O 8+δ using a low-cost scanning electron microscope

In this work, we investigate optically active nanostructures on Bi 2 Sr 2 CaCu 2 O 8+δ (BSCCO). The nanostructures are constructed from carbon nanocrystals formed using the electron beam of a scanning electron microscope to locally decompose residual hydrocarbons at the BSCCO surface. Atomic force microscopy is used to characterize the dimensions of these structures as a function of beam exposure time. This technique has been used to induce localized intercalation to form carbon nanoparticles up to tens of nanometers into a variety of 2D materials. Cross-sectional scanning transmission electron microscopy in BSCCO shows the presence of these carbon deposits on the surface, but there is no evidence they penetrate into the BSCCO substrate. Therefore, this technique is not suitable for localized intercalation in BSCCO, and it does have promise as a cost-effective approach to pattern nanometer scale etch stops.

2D materials↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗
DOE OSTI2024

Establishing accretion flares from supermassive black holes as a source of high-energy neutrinos

ABSTRACT The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

Astronomy & Astrophysics↗
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↗
DOE OSTI2026

SILICON CARBIDE THERMOMETRY USING RAMAN SPECTROSCOPY ON IRRADIATED TRISO PARTICLES

Silicon carbide (SiC) passive thermometry has emerged as a promising post-irradiation examination (PIE) technique for estimating irradiation temperature near the end of irradiation. While dilatometry techniques have been traditionally used to analyze prismatic samples after irradiation, Raman spectroscopy has recently been shown to provide comparable results by analyzing Raman-active phonon modes. In this work, Raman spectroscopy has been applied to the SiC layer of cross-sectioned irradiation tristructural isotropic (TRISO) particles from the AGR-5/6/7 experiment to evaluate the feasibility of particle-scale passive thermometry. Two-dimensional Raman mapping was used to measure the position of the SiC longitudinal optical (LO) phonon, which was then converted to an apparent irradiation temperature using a previously established empirical correlation. Particles from AGR-5/6/7 Compacts 2-2-1 and 5-1-3 were selected as their calculated time-averaged, volume-averaged (TAVA) temperatures (828°C and 706°C, respectively) fall within the range of the sensitivity of the experimental approach. The use of SiC thermometry was anticipated to confirm or highlight potential deviations from calculated end-of-life TAVA temperature across compacts. Four particles from Compact 2-2-1 and two particles from Compact 5-1-3 were selected based on 110mAg inventory (either some measurable activity or below the minimum detection limit) which is commonly used as an indicator of in-pile temperature variation of particles within the same compact. Across every particle selected it was determined that the average LO peak position was located around 968 cm-1 to 969 cm-1. Using the previously determined empirical correlation, this corresponds to an irradiation temperature around 950°C to 966°C, which does not align with the reported TAVA temperature values. This discrepancy in apparent irradiation temperatures likely reflects a combination of uncertainties in the calculated particle temperatures and differences in irradiation history between the present specimens and those used to establish the empirical Raman calibration such as neutron flux (damage rate) and SiC microstructure (as fabricated and irradiated).

Vawdrey, Josh [ORNL]↗
DOE OSTI2026

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗
DOE OSTI2026

Tailored Solvent Treatment for Optimized Production of Upcycled Anodes from End-Of-Life Li-Ion Batteries

Recycling processes for lithium-ion batteries typically overlook graphite because of its lower market value relative to that of transition-metal-containing cathode materials. However, graphite recovered from cycled lithium-ion batteries holds additional engineered value associated with the solid-electrolyte interphase (SEI). The SEI contributes critical electronic passivation of the graphite surface but becomes highly resistive with extended cycling, yielding poor cell performance. In this work, we apply tailored solvent treatment to end-of-life (EOL) graphite anodes to selectively remove adverse SEI components while retaining beneficially passivating species. We evaluate a series of polar protic solvents to achieve targeted removal of SEI components and control selectivity through rational variation in solvent properties. The physiochemical properties of treatment solvents correlate with both the retained SEI composition and the corresponding electrochemical performance of solvent-treated “upcycled” graphite anodes. Within the initial set of solvents evaluated, top-performing candidates show capacity and Coulombic efficiency nearly equivalent to those of an analogous pristine anode, as well as promising electrochemical performance enhancement with regard to irreversible capacity-loss metrics. This study establishes critical design principles for an optimized anode upcycling method that enhances the value of recycled graphite by retaining and upgrading the SEI.

25 ENERGY STORAGE↗
DOE OSTI2025

Multi-contrast machine learning improves schistosomiasis diagnostic performance

Schistosomiasis currently affects over 250 million people and remains a public health burden despite ongoing global control efforts. Conventional microscopy is a practical tool for diagnosis and screening ofSchistosoma haematobium, but identification of eggs requires a skilled microscopist. Here we present a machine learning (ML)-based strategy for automated detection ofS. haematobiumthat combines two imaging contrasts, brightfield (BF) and darkfield (DF), to improve diagnostic performance. We collected BF and DF images of urine samples, many of them containingS. haematobiumeggs, during two different field studies in Côte d’Ivoire using a mobile phone-based microscope, the SchistoScope. We then trained separate egg-detection ML models and compared the patient-level performance of BF and DF models alone to combinations of BF and DF models, using annotations from trained microscopists as the gold standard. We found that models trained on DF images, and almost all BF and DF combinations, performed significantly better than models trained on BF images only. When models were trained on images from the first field study (n = 349 patients, 748 images of each contrast), patient-level classification performance on patient images from the second study (n = 375 patients, 752 images of each contrast) met the WHO Diagnostic Target Product Profile (TPP) sensitivity and specificity for the monitoring and evaluation use case (sensitivity for all models and combinations was >75% when evaluated at a confidence score threshold that resulted in specificity >96.5%). When we used images from both field studies for the training set, performance of the models was improved. Overall, this work shows that the use of DF and BF increases the performance of ML models on images from devices with low-cost optics, while retaining the portability, power, and time-to-results of the WHO’s diagnostic TPP. DF requires no additional sample preparation and does not increase the complexity of the imaging system. It thus offers a practical means to improve performance of automated diagnostics forS. haematobiumas well as other microscopy-based diagnostics.

Infectious Diseases↗
DOE OSTI2026

Solid neon as a noise-resilient host for electron qubits above 100 mK

Solid neon can be used as a solid host for single-electron qubits. At temperatures of around 10 mK, electron-on-solid-neon charge qubits exhibit long coherence times and high operation fidelities. However, a systematic characterization of the noise features of such systems is needed for the development of scalable quantum information architectures. Here, in this work, we show that solid neon can be used as a noise-resilient host for electron qubits above 100 mK. We examine the resilience of solid neon against charge and thermal noise when electron-on-solid-neon charge qubits are operated away from the charge-insensitive sweet spot and at elevated temperatures. We show that the extracted high-frequency charge noise density of electron-on-solid-neon qubits, projected as voltage fluctuations on nearby electrodes, is between 10 −4 μV 2 Hz −1 and 10 −6 μV 2 Hz −1 at 0.01 MHz to 1 MHz, which is comparable to common semiconductor hosts. We also show that the electron-on-solid-neon charge qubits operating at frequencies of around 5 GHz can maintain echo coherence times of over 1 μs at temperatures up to 400 mK.

42 ENGINEERING↗
DOE OSTI2025

Monitoring Depolymerization in Mesopores Using Dynamic Properties of Polymeric Melt Accessed via Dielectric Spectroscopy

Traditional design principles for heterogeneous catalysis guide the use of catalytic particles with mesosized (∼2–50 nm) pores to increase the number of surface-active sites by way of an increased surface area. However, the entry of long-chain polymers into such pores may be significantly limited by the size and entanglement of polymers in the melt state, thereby decreasing the number of accessible sites. Assessment of catalyst performance from traditional reactor-based studies averages over intrapore reaction events as well as reactions on the surface of a particle, resulting in an inability to distinguish between differences in site accessibility and activity. Techniques that assess the intrapore performance can inform the design of future heterogeneous catalysts for polymer upcycling. In this work, we demonstrate the use of broadband dielectric spectroscopy to monitor depolymerization of a polymer melt within mesopores via changes in the segmental relaxation time scale of amorphous polymer chains. In particular, we highlight the use of an anodic aluminum oxide (AAO) membrane as a readily available model for catalyst pores with a well-characterized pore morphology. The decrease in the segmental relaxation (α-relaxation) time of the melt with increasing chain scission emerges as a measure of the extent of polymer deconstruction inside mesopores. To demonstrate the utility of this technique, we demonstrate the decomposition of two commercial poly(propylene carbonate) polymers with different decomposition rates within mesopores. As the polymers depolymerize, their segmental relaxation time decreases as the molecular weight decreases (as predicted by the Fox–Flory equation). The BDS-measured change in segmental relaxation time mirrors the expected trend based on change in molecular weight measured by size exclusion chromatography.

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

Design and Realization of Ohmic and Schottky Interfaces for Oxide Electronics

Understanding band alignment and charge transfer at complex oxide interfaces is critical to tailoring and utilizing their diverse functionality. Toward this goal, both Ohmic- and Schottky-like charge transfers at oxide/oxide semiconductor/metal interfaces are designed and experimentally validated. A method for predicting band alignment and charge transfer in ABO 3 perovskites is utilized, where previously established rules for simple semiconductors fail. The prototypical systems chosen are the rare class of oxide metals, SrBO 3 with B = V–Ta, when interfaced with the multifaceted semiconducting oxide, SrTiO 3 . For B = Nb and Ta, it is confirmed that a large accumulation of charge occurs in SrTiO 3 due to the higher energy Nb and Ta states relative to Ti. Furthermore, this gives rise to a high mobility metallic interface, which is an ideal epitaxial oxide/oxide Ohmic contact. On the contrary, for B = V, there is no charge transfer into the SrTiO 3 interface, which serves as a highly conductive epitaxial gate metal. Going beyond these specific cases, this work opens the door to integrating the vast phenomena of ABO 3 perovskites into a wide range of practical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗