Search NASA⌕ Search

SEARCH · Search NASA

Results for “likelihood ratio test”

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.

300 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development of Electrolytes under Lean Condition in Lithium–Sulfur Batteries

Lithium–sulfur (Li–S) batteries stand out as one of the promising candidates for next-generation electrochemical energy storage technologies. A key requirement to realize high-specific-energy Li–S batteries is to implement low amount of electrolyte, often characterized by the electrolyte/sulfur (E/S) ratio. Low E/S ratio aggravates the known challenges for Li–S batteries and introduces new ones originated from the high concentration of polysulfides in limited electrolyte reservoir. Here, in this review, the connections between the fundamental properties of electrolytes and the electrochemical/chemical reactions in Li–S batteries under lean electrolyte condition are elucidated. The emphasis is on how the solvating properties of the electrolyte affect the fate of polysulfides. Built upon the mechanistic analysis, different strategies to design lean electrolytes to improve the overall process of Li–S reactions and Li anode protection are discussed.

25 ENERGY STORAGE↗
DOE OSTI2026

Single particle – MC-ICP-MS for isotopic analysis of uranium particles

Single particle – multi-collector – inductively coupled plasma – mass spectrometry (SP-MC-ICP-MS) was employed to measure a suspension of 1 µm U3O8 particles (∼1.3 pg total U/particle) to determine their individual isotopic compositions of 234U/238U, 235U/238U, and 236U/238U. The effects of different detector combinations for 235U and 238U, including secondary electron multipliers (SEM) and Faraday detectors (1011 and 1013 Ω amplifiers), were explored for accuracy and precision optimization on the observed 235U/238U. The minor isotopic ratios (i.e., 234U/238U and 236U/238U) were analyzed such that the 234U and 236U were monitored on SEM detectors and the 238U was monitored on a Faraday (1011 Ω) detector. Various integration times (5, 10, 25, and 50 ms) were investigated in all detector configurations to gain a better understanding of their impact on sensitivity, accuracy, and precision. For 235U/238U ratios of 1 µm U3O8 particles, a dual Faraday detector measurement with 1011 Ω was the optimal choice; measurement of 1021 particles yielded an average 235U/238U ratio of 0.00170 (14), a −1.8% relative difference (% RD) from the reference value. The minor isotopic compositions were determined to be 0.0000070 (14) and 0.0000758 (48) for the 234U/238U and 236U/238U, respectively. These measurements correspond to <8% and <1% RD from their reference value for the 234U/238U and 236U/238U, respectively. SP-MC-ICP-MS was also able to provide insight into measurement sensitivity. In these individual particles, merely 15 and 165 atto-grams (ag) of 234U and 236U were present (calculated). Initial limits of detection for SP-MC-ICP-MS were determined to be ∼1.0 ag (when measured via SEM detectors). This valuable approach is applicable to areas including nuclear forensics, nuclear safeguards, and geochemical analysis, which require high-precision measurements of uranium within micron-sized particles.

Manard, Benjamin [ORNL] (ORCID:0000000207400627)↗
DOE OSTI2025

Effect of Preparation Conditions of Fe@SiO2 Catalyst on Its Structure Using High-Pressure Activity Studies in a 3D-Printed SS Microreactor

Fischer–Tropsch synthesis (FTS) in a 3D-printed stainless steel (SS) microchannel microreactor was investigated using Fe@SiO2 catalysts. The catalysts were prepared by two different techniques: one pot (OP) and autoclave (AC). The mesoporous structure of the two catalysts, Fe@SiO2 (OP) and Fe@SiO2 (AC), ensured a large contact area between the reactants and the catalyst. They were characterized by N2 physisorption, H2 temperature-programmed reduction (H2-TPR), scanning electron microscopy (SEM), transmission electron microscopy (TEM), X-ray diffraction (XRD), Fourier-transform infrared spectroscopy (FTIR), X-ray photoelectron microscopy (XPS), and thermogravimetric analysis–differential scanning calorimetry (TGA-DSC) techniques. The AC catalyst had a clear core–shell structure and showed a much greater surface area than that prepared by the OP method. The activities of the catalysts in terms of FTS were studied in the 200–350 °C temperature range at 20-bar pressure with a H2/CO molar ratio of 2:1. The Fe@SiO2 (AC) catalyst showed higher selectivity and higher CO conversion to olefins than Fe@SiO2 (OP). Stability studies of both catalysts were carried out for 30 h at 320 °C at 20 bar with a feed gas molar ratio of 2:1. The Fe@SiO2 (AC) catalyst showed higher stability and yielded consistent CO conversion compared to the Fe@SiO2 (OP) catalyst.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Electron temperature relations and the direct N, O, Ne, S, and Ar abundances of 49 959 star-forming galaxies in DESI data release 2

We present the largest direct-method abundance catalogue of galaxies to date, containing measurements of 49 959 star-forming galaxies at z<0.96 from DESI (Dark Energy Spectroscopic Instrument) data release 2. By directly measuring electron temperatures across multiple ionization zones, we provide constraints on a number of electron temperature relations. Using the temperature measurements, we derive reliable abundances for N, O, Ne, S, and Ar, and measure the evolution of abundances and abundance ratios of as a function of metallicity and other galaxy properties. Our measurements include direct oxygen abundances for 49 507 galaxies, leading to the discovery of the two most metal-poor galaxies in the nearby Universe, with oxygen abundances of 12+log⁡(O/H)=6.77−0.03+0.03 dex (1.2 per cent Z⊙⁠) and 12+log⁡(O/H)=6.81−0.04+0.04 dex (1.3 per cent Z⊙⁠). We identify a rare outlier population of 24 galaxies with high-N/O ratios at low metallicity, reminiscent of galaxy abundances observed in the early Universe. We find the Ne/O ratio is constant at low metallicity but increases gradually at 12+log(O/H)>8.105±0.004 dex. We show that the S/O and Ar/O abundance ratios are strongly correlated, consistent with the expected additional Type Ia enrichment channel for S and Ar. In this work, we present an initial survey of the key properties of the sample, with this data set serving as a foundation for extensive future work on galaxy abundances at low redshift.

Scholte, D. [Edinburgh U., Inst. Astron.] (ORCID:0↗
DOE OSTI2024

Efficient Leaching of Metal Ions from Spent Li-Ion Battery Combined Electrode Coatings Using Hydroxy Acid Mixtures and Regeneration of Lithium Nickel Manganese Cobalt Oxide

Extensive use of Li-ion batteries in electric vehicles, electronics, and other energy storage applications has resulted in a need to recycle valuable metals Li, Mn, Ni, and Co in these devices. In this work, an aqueous mixture of glycolic and lactic acid is shown as an excellent leaching agent to recover these critical metals from spent Li-ion laptop batteries combined with cathode and anode coatings without adding hydrogen peroxide or other reducing agents. An aqueous acid mixture of 0.15 M in glycolic and 0.35 M in lactic acid showed the highest leaching efficiencies of 100, 100, 100, and 89% for Li, Ni, Mn, and Co, respectively, in an experiment at 120 °C for 6 h. Subsequently, the chelate solution was evaporated to give a mixed metal-hydroxy acid chelate gel. Pyrolysis of the dried chelate gel at 800 °C for 15 h could be used to burn off hydroxy acids, regenerating lithium nickel manganese cobalt oxide, and the novel method presented to avoid the precipitation of metals as hydroxide or carbonates. The Li, Ni, Mn, and Co ratio of regenerated lithium nickel manganese cobalt oxide is comparable to this metal ratio in pyrolyzed electrode coating and showed similar powder X-ray diffractograms, suggesting the suitability of α-hydroxy carboxylic acid mixtures as leaching agents and ligands in regeneration of mixed metal oxide via pyrolysis of the dried chelate gel.

Electrochemistry↗
DOE OSTI2026

Strong gradient neoclassical transport in the plateau regime

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

fusion plasma↗
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2025

Novel Zwitterionic Polyurethane-in-Salt Electrolytes with High Ion Conductivity, Elasticity, and Adhesion for High-Performance Solid-State Lithium Metal Batteries

This study presents a novel polymer-in-salt (PIS) zwitterionic polyurethane-based solid polymer electrolyte (zPU-SPE) that offers high ionic conductivity, strong interaction with electrodes, and excellent mechanical and electrochemical stabilities, making it promising for high-performance all solid-state lithium batteries (ASSLBs). The zPU-SPE exhibits remarkable lithium-ion (Li+) conductivity (3.7 × 10⁻⁴ S cm−1 at 25 °C), enabled by exceptionally high salt loading of up to 90 wt.% (12.6 molar ratio of Li salt to polymer unit) without phase separation. It addresses the limitations of conventional SPEs by combining high ionic conductivity with a Li+ transference number of 0.44, achieved through the incorporation of zwitterionic groups that enhance ion dissociation and transport. The high surface energy (338.4 J m−2) and elasticity ensure excellent adhesion to Li anodes, reducing interfacial resistance and ensuring uniform Li+ flux. When tested in Li||zPU||LiFePO₄ and Li||zPU||S/C cells, the zPU-SPE demonstrated remarkable cycling stability, retaining 76% capacity after 2000 cycles with the LiFePO4 cathode, and achieving 84% capacity retention after 300 cycles with the S/C cathode. Molecular simulations and a range of experimental characterizations confirm the superior structural organization of the zPU matrix, contributing to its outstanding electrochemical performance. The findings strongly suggest that zPU-SPE is a promising candidate for next-generation ASSLBs.

Wang, Kun↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2026

Advanced blade-shaped thermal energy storage device: Development and application

Thermal energy storage (TES) using phase change materials (PCMs) is a promising approach for capturing and reusing excess thermal energy, yet widespread adoption is limited by low thermal conductivity, bulky configurations, and inadequate scalability. Here, this study presents a modular, blade-shaped TES prototype designed to address these challenges. The device integrates a lightweight aluminum shell, an embedded serpentine coil for active or passive heat exchange, and a cost-effective corrugated metal mesh for enhanced PCM thermal conductivity. With thickness-to-length and thickness-to-width ratios of 0.03 and 0.08, respectively, the blade-shaped TES achieves a compact, modular form factor suitable for space-constrained applications. Experimental testing demonstrated the efficient charge and discharge behavior of blade-shaped TES, capturing PCM superheating, phase-change transitions, and subcooling dynamics, with charging and discharging efficiencies of 94.9% and 94.6%, respectively. Also, the system can potentially achieve higher energy density than that of conventional TES designs. When integrated into a household refrigerator during the study, three blade-shaped TES modules successfully shifted 100% of peak-time compressor operation to off-peak hours, reducing energy consumption while maintaining more stable compartment temperatures. The blade-shaped TES's thin geometry, modularity, and enhanced thermal performance support scalable deployment across residential, commercial, and industrial applications, providing a versatile, cost-effective solution for high-efficiency, demand-flexible thermal energy management.

Blade-shaped↗
DOE OSTI2026

Nuclear Structure and Decay Data for A=35 Isobars

Here, this work presents a comprehensive and critical evaluation of experimental nuclear spectroscopic data from reactions and decays for all 11 known nuclides with mass number 35 (Ne, Na, Mg, Al, Si, P, S, Cl, Ar, K, Ca). Recommended values are produced for level energies, spins and parities, half-lives, and radiation properties including energies, branching ratios, and multipolarities of γ rays, as well as characteristics of β radiation decays, based on a rigorous assessment of all available experimental data. Discrepancies among existing results are carefully addressed. This work supersedes earlier full evaluations of A=35 published by 2011Ch48, 1990En08 (also 1998En04 update) and 1978En02.

Sun, Lijie [Michigan State University, East Lansin↗
DOE OSTI2026

Working Fluid Characterization and Performance Assessment of Subcritical Organic Rankine Cycles Based on the Lee–Kesler Approach for Energy Recovery

Here, a generalized model using the Lee–Kesler approach based on the corresponding states principle is developed to assess the performance of subcritical Organic Rankine Cycles operating with different working fluids. Each fluid is characterized by five parameters: the acentric factor, critical temperature, critical pressure, molar mass, and the ideal-gas ratio of specific heats at the critical temperature. The model was developed using the compressibility factor modified version of the Benedict–Webb–Rubin equation proposed by Lee and Kesler and the enthalpy and entropy functions to calculate thermodynamic state properties. The model was validated by comparing the results calculated with the model and working fluid thermodynamic properties obtained with the CoolProp database. This comparison was conducted for 91 working fluids, obtaining a relative error below 5% for 88 out of the 91 fluids (∼97%). A generalized parametric study was conducted to determine the influence of the pinch point and each fluid parameter on the performance of Organic Rankine Cycle (ORC) systems. It was found that efficiency increases with critical temperature, ideal-gas ratio of specific heats at the critical temperature, and acentric factor, reaching up to 13%. The developed model enables the evaluation of ORC system performance for existing working fluids. It also allows the formulation and evaluation of new fluids to enhance the performance of the ORC while retrieving energy from any kind of source; and likewise, the methodology can be applied to other power generation cycles.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2026

Novel PtNi single-atom–nanocluster (SA–NC) ensembles promote Tafel kinetics and ampere-class AEM hydrogen evolution

The development of efficient, durable, and low-PGM electrocatalysts for the hydrogen evolution reaction (HER) in alkaline media is critical for next-generation electrolysis technologies. We report a facile two-step synthesis of highly dispersed PtNi and PtNi-nitride nanoclusters (NCs) (~2.3 nm) with ultralow Pt content (0.5 at.%) anchored on N-doped Vulcan carbon. Structural and compositional characterization via XAS, XPS, HAADF-STEM, HRTEM, and EDS mapping established key structure–activity relationships across varying Pt/Ni ratios and pyrolysis temperatures. The Pt0.5Ni0.5/C-750 catalyst, an ensemble of PtNi M-N-C type single-atom (SA) moieties with neighboring PtNi nanoclusters (NC), exhibited superior HER performance in alkaline media, achieving overpotentials of 30, 115, and 210 mV at 10, 100, and 500 mA cm−2, respectively. Despite at a lower Pt content, this novel SA–NC ensemble outperformed commercial Pt/C by ~36%. A standardized literature comparison with contemporary Pt- and Ru-doped analogues reveals the as-prepared Pt0.5Ni0.5/C-750 to sit at the apex of Tafel-limited kinetics and low overpotential at 100 mA cm−2. Tafel-limited Tafel slopes in both alkaline and acidic regimes confirm favorable proton recombination kinetics. Mass activities at 200 mV reached 13.8 and 18.84 A mgPt−1 in alkaline and acidic media, respectively. However, excessive nitridation (e.g., at 650 °C) adversely altered Pt electronic structure and HER kinetics. While Ni enhanced alkaline HER, acidic HER favored Ni-free analogues. Pt0.5Ni0.5/C-750 also demonstrated robust temperature responsiveness and 300-h operational stability at high current densities (0.5–1.0 A cm−2) in MEA tests. This work presents a scalable strategy for designing thermally responsive, durable, and compositionally tunable NC catalysts with neighboring SA moieties for alkaline electrolysis.

AEMWE↗
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Dark Energy Survey Year 6 results: Clustering redshifts and importance sampling of self-organized-maps 𝑛⁡(𝑧) realizations for 3 × 2 ⁢pt samples

This work is part of a series establishing the redshift framework for the 3 × 2 ⁢pt analysis of the Dark Energy Survey Year 6 (DES Y6). For DES Y6, photometric redshift distributions are estimated using self-organizing maps (SOMs), calibrated with spectroscopic and many-band photometric data. To overcome limitations from color-redshift degeneracies and incomplete spectroscopic coverage, we enhance this approach by incorporating clustering-based redshift constraints (clustering-z, or WZ) from angular cross-correlations with BOSS and eBOSS galaxies and eBOSS quasar samples. We define a WZ likelihood and apply importance sampling to a large ensemble of SOM-derived 𝑛⁡(𝑧) realizations, selecting those consistent with the clustering measurements to produce a posterior sample for each lens and source bin. The analysis uses angular scales corresponding to 1.5–5 Mpc to optimize signal-to-noise ratio while mitigating modeling uncertainties and marginalizes over redshift-dependent galaxy bias and other systematics informed by the N-body simulation CARDINAL . While a sparser spectroscopic reference sample limits WZ constraining power at 𝑧 >1.1, particularly for source bins, we demonstrate that combining SOM with WZ improves redshift accuracy and enhances the overall cosmological constraining power of DES Y6. As a result, we estimate an improvement in 𝑆 8 of approximately 10% for cosmic shear and 3 ×2⁢pt analysis, primarily due to the WZ calibration of the source samples.

Cosmological parameters↗
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.↗
DOE OSTI2025

Influence of gradation on failure mode of saturated sand at particle scale

Extensive research has been reported in the literature to characterize the failure mode of dry sand using various experimental techniques such as surface optical imaging, photo-elastic materials, three-dimensional (3D) computed tomography (CT), and 3D synchrotron micro-computed tomography (SMT). However, there is a limited literature about the behavior of saturated sand. This paper presents the results of axisymmetric triaxial compression (ATC) experiments that were conducted on saturated sand specimens. The behavior of specimens composed of a uniform sand with grain size between US sieves #40 and #50 is compared to specimens conducted on the same sand that has a wider gradation. 3D SMT technique was used to acquire 3D scans while shearing the specimens to probe localized events that are completely missed or misinterpreted when analyzing ATC measurements based on global standard measurements. The results show a higher effective principal stress ratio (EPSR) for the non-uniform specimen and a thicker shear band when compared to uniform specimen.

Elnur, Mohammed [University of Tennessee]↗