Search NASA⌕ Search

SEARCH · Search NASA

Results for “Vehicle-to-grid”

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.

1,494 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hamiltonian switching control of noisy bipartite qubit systems

Abstract We develop a Hamiltonian switching ansatz for bipartite control that is inspired by the quantum approximate optimization algorithm, to mitigate environmental noise on qubits. We demonstrate the control for a central spin coupled to bath spins via isotropic Heisenberg interactions, and then make physical applications to the protection of quantum gates performed on superconducting transmon qubits coupling to environmental two-level-systems (TLSs) through dipole-dipole interactions, as well as on such qubits coupled to both TLSs and a Lindblad bath. The control field is classical and acts only on the system qubits. We use reinforcement learning with policy gradient to optimize the Hamiltonian switching control protocols, using a fidelity objective for specific target quantum gates. We use this approach to demonstrate effective suppression of both coherent and dissipative noise, with numerical studies achieving target gate implementations with fidelities over 0.9999 (four nines) in the majority of our test cases and showing improvement beyond this to values of 0.999 999 999 (nine nines) upon a subsequent optimization by GRadient Ascent Pulse Engineering (GRAPE). We analyze how the control depth, total evolution time, number of environmental TLS, and choice of optimization method affect the fidelity achieved by the optimal protocols and reveal some critical behaviors of bipartite control of quantum gates.

Physics↗
DOE OSTI2024

Active Learning for Rapid Targeted Synthesis of Compositionally Complex Alloys

The next generation of advanced materials is tending toward increasingly complex compositions. Synthesizing precise composition is time-consuming and becomes exponentially demanding with increasing compositional complexity. An experienced human operator does significantly better than a novice but still struggles to consistently achieve precision when synthesis parameters are coupled. The time to optimize synthesis becomes a barrier to exploring scientifically and technologically exciting compositionally complex materials. This investigation demonstrates an active learning (AL) approach for optimizing physical vapor deposition synthesis of thin-film alloys with up to five principal elements. We compared AL-based on Gaussian process (GP) and random forest (RF) models. The best performing models were able to discover synthesis parameters for a target quinary alloy in 14 iterations. We also demonstrate the capability of these models to be used in transfer learning tasks. RF and GP models trained on lower dimensional systems (i.e., ternary, quarternary) show an immediate improvement in prediction accuracy compared to models trained only on quinary samples. Furthermore, samples that only share a few elements in common with the target composition can be used for model pre-training. We believe that such AL approaches can be widely adapted to significantly accelerate the exploration of compositionally complex materials.

Chemistry↗
DOE OSTI2026

Enter the AHU (36th Chamber of ASHRAE): A Multi-site Field Study of ASHRAE G36

Despite being recognized as the best practice for advanced building controls, ASHRAE Guideline 36 (G36) has seen slow adoption in retrofit cases. Decisionmakers lack credible field evidence to justify the time and person-power investment. Most prior analyses have relied on software simulations, which overlook implementation challenges and fail to persuade owners to move from models to real-world deployment. This paper presents a multi-site field study of G36 performance, drawing on measured results from 17 projects across diverse building types and climate zones. The analysis disaggregates outcomes by the most widely adopted air handling unit (AHU) based G36 strategies, including trim-and-respond approaches to supply air temperature (SAT) and duct static pressure (DSP) reset, and economizer controls. Results for controls re programming implementations are encouraging, with HVAC savings ranging from 2% - 49%, with a median of 18%. The range aligns with simulation study findings showing 1% - 46% savings through these three strategies. These findings show that even without capital investment in control infrastructure upgrades, existing building owners are reaping significant benefits from updating HVAC sequences of operation to industry best-practice solutions. The results give practitioners and decisionmakers a reference point on what to expect, which strategies deliver, and how field performance may compare to simulation.

Deshpande, Reva↗
DOE OSTI2026

Design of a Cryogenic Assembly for Comparing Qubit Infrared Packaging Inefficiencies

It is well known that everything emits infrared (IR) light in the form of thermal radiation. IR radiation on superconducting qubits carries enough energy to cause decoherence and so-called “quasiparticle poisoning.” Therefore, proper IR shielding is needed to keep the qubit in a high coherence state. However, trying to extract IR-specific information from a qubit is difficult, so a microwave kinetic inductance detector (MKID) can be used instead. We show that the IR shielding inefficiencies of MKID packages can only be compared with other packaging inefficiencies. We create a cryogenic assembly CAD design, consisting of magnetic shield cans, a copper mounting plate and mounts, and an MKID enclosure, which thermalizes the assembly to dilution refrigerator temperatures. Phase 1 of this design involves using a narrowband IR laser source to test the IR shielding of packaging at the laser’s wavelength. When the best IR shielding packaging is determined, the superconducting qubit can replace the MKID, allowing the qubit to avoid decoherence due to IR radiation and boost its performance and coherence lifetime.

Mailliard, Robert J. [Fermilab] (ORCID:00090000277↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (↗
DOE OSTI2026

Dielectric Spectroscopy Cable NDE for Unequal Aging Over Different Lengths

This Pacific Northwest National Laboratory milestone report assesses the effect of unequal aging over different lengths of cable on dielectric spectroscopy (DS) bulk impedance measurements from the cable end. DS measurements can indicate the remaining useful life of cables; however, most benchmark tests are based on accelerated aging of the entire cable. If only a portion of the cable is exposed to a harsh environment (as is frequently the case), the DS measurement will indicate a significantly less-aged cable than if the entire cable were exposed to a uniform environmental stress. The DS measurement is primarily related to the intrinsic insulation material relative permittivity—typically between 1.5 and 3.0. Insulation aging increases permittivity and, correspondingly, the cable’s overall capacitance measured from the cable end. To judge the health of the overall cable system, operators are primarily interested in the condition of the most severely aged section since that is where a failure is most likely to occur. This can be expressed as the dimensionless relative permittivity of the severely aged section divided by the permittivity of a pristine cable. This relative permittivity is an intrinsic material parameter independent of cable length and can serve as a quantitative indication of the insulation condition. If the relative percentages of cable exposed to a harsh and benign environment are known, the DS measurement of the entire cable length, including both pristine and aged sections, can be compensated to predict the relative permittivity (relative to pristine permittivity) of the most severly aged segment. A test was performed on three 50 ft (16 m) cables with 10%, 50%, and 90% of the cable inside a thermal aging oven. Predictably, the aging effect on the DS response was greatest in the 90% sample, followed by the 50% sample, and then the 10% sample. The more interesting result is the ratio of the aged insulation permittivity to the pristine insulation permittivity. Based on prior work, an end-of-life cable has a value of approximately 1.35. If the relative percentages of pristine and severely aged cable lengths are known, and the total cable-length capacitance of the initial pristine cable and the combined pristine plus aged cable are known, the relative aged permittivity (ratio of aged section permittivity / pristine section permittivity) can be estimated. The relative cable lengths exposed to a harsh environment may be known or estimated based on plant layout drawings or cable reflectometry tests that can locate the oven entry and exit points. The compensation was verified using a lumped-parameter model to predict the aged-segment permittivity e2/e1, where e2 is the aged-segment permittivity and e1 is the pristine-segment permittivity. For all three percentages of aged cable lengths, the ratio was quite similar, as would be expected since all cable segments were exposed to the same aging environment. This ratio of aged permittivity to pristine permittivity can be compared to other damage indicating tests, including elongation at break, to assess cable condition.

ARENA Test Bed↗
DOE OSTI2026

Resolving root causes of experiment discrepancies guided by machine learning

Abstract Scientists rely on accurate experimental data to explain nature and then harness this knowledge for applications addressing human needs. However, discrepancies between experiments of the same observable can impede scientific progress if one does not understand the underlying causes. Here, we developed a process that unravels data discrepancies by first using Bayesian machine learning to relate discrepancies to few of many, potentially biasing metadata features that encode experiment procedures. This machine learning output guides human experts to study discrepancy causes by simulating suspicious aspects of historical experiments or designing modern ones to address open questions. The study findings then lead to rejecting or correcting historical data on firm scientific bases. This process is demonstrated for the energy spectrum of neutrons emitted promptly (<1 ns) after fission of 252 Cf, a trusted nuclear physics Standard. It reduces the spread in experimental 252 Cf spectra by up to a factor of 6.

Neudecker, D. (ORCID:0000000339200627)↗
DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗
DOE OSTI2026

Dimensional Reduction for Sampled Priors and Application to Photometric Redshift Distributions

A typical Bayesian inference on the values of some parameters of interest q from some data D involves running a Markov Chain (MC) to sample from the posterior $p$($q$,$n$|$D$) $\propto$ $\mathcal{L}$($D$|$q$,$n$)$p$(q)$p$($n$), where n are some nuisance parameters with a separable prior. In some cases, the nuisance parameters are high-dimensional, and their prior p(n) is itself defined only by a set of samples that have been drawn from some other MC. The MC for the posterior will typically require evaluation of p(n) at arbitrary values of n, i.e., one needs to provide a density estimator over the full n space from the provided samples. But the high dimensionality of n hinders both the density estimation and the efficiency of the MC for the posterior. We describe a solution to this problem: a linear compression of the n space into a much lower-dimensional space u, which projects away directions in n space that cannot appreciably alter $\mathcal{L}$. The algorithm for doing so is a slight modification to principal components analysis, and is less restrictive on p(n) than other proposed solutions to this issue. We demonstrate this “mode projection” technique using the analysis of 2-point correlation functions of weak lensing fields and galaxy density in the Dark Energy Survey, where n is a binned representation of the redshift distribution n(z) of the galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
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

In-Situ Visualization of a Growing Brittle Crack in Aluminum Oxynitride Using Synchrotron X-Rays and the Double-Cleavage Drilled Compression Geometry

Brittle fracture is difficult to study in situ due to the speed of a growing crack and the often-catastrophic nature of failure in brittle materials. As a result, the influence of microstructural considerations, such as orientation, grain boundary locations, and strain field, on the crack path remains poorly understood. Presented in this study is a method addressing this knowledge gap, which utilizes the double-cleavage drilled compression geometry to achieve quasi-stable fracture in aluminum oxynitride (AlON). Synchrotron X-ray micro-computed tomography is used to characterize the crack shape and length, while high-energy diffraction microscopy provides information on the strains, orientations, and shapes of grains in the microstructure surrounding the crack tip. During testing, the crack grew in discrete and irregular jumps while the fracture toughness falls within reported ranges. The crack in AlON is found to have no greater tendency to crack intergranularly as compared to transgranularly, and grains which are cracked transgranularly do not display a trend in orientation or stress when compared to those around which the crack followed a grain boundary. The high resolution of the crack path and microstructural data provides a path forward for modeling and understanding 3D brittle fracture.

Gorske, Sara F.↗
DOE OSTI2026

Ultrafast Raman probe of the photoinduced superconducting to normal state transition in the cuprate Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝛿

Here, we report an ultrafast Time-Resolved Raman scattering study of the out-of-equilibrium photoinduced dynamics across the superconducting to normal state phase transition of the cuprate Bi 2 ⁢Sr 2⁢ CaCu 2 ⁢O 8+𝛿 . Using the polarization-resolved momentum space selectivity of Raman scattering, we track the superconducting condensate destruction and recovery dynamics with subpicoseconds time resolution in the antinodal region of the Fermi surface where the superconducting gap is maximum. Leveraging ultrafast Raman thermometry, we find a significant dichotomy between the superconducting condensate and the quasiparticle temperature dynamics near the antinodes, which cannot be framed in terms of a single effective electron temperature. The present work demonstrates the ability of Time-Resolved Raman scattering to selectively probe out-of-equilibrium pathways of different electronic subdegrees of freedom during a photoinduced phase transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mechanochemically responsive polymer enables shockwave visualization

Abstract Understanding the physical and chemical response of materials to impulsive deformation is crucial for applications ranging from soft robotic locomotion to space exploration to seismology. However, investigating material properties at extreme strain rates remains challenging due to temporal and spatial resolution limitations. Combining high-strain-rate testing with mechanochemistry encodes the molecular-level deformation within the material itself, thus enabling the direct quantification of the material response. Here, we demonstrate a mechanophore-functionalized block copolymer that self-reports energy dissipation mechanisms, such as bond rupture and acoustic wave dissipation, in response to high-strain-rate impacts. A microprojectile accelerated towards the polymer permanently deforms the material at a shallow depth. At intersonic velocities, the polymer reports significant subsurface energy absorption due to shockwave attenuation, a mechanism traditionally considered negligible compared to plasticity and not well explored in polymers. The acoustic wave velocity of the material is directly recovered from the mechanochemically-activated subsurface volume recorded in the material, which is validated by simulations, theory, and acoustic measurements. This integration of mechanochemistry with microballistic testing enables characterization of high-strain-rate mechanical properties and elucidates important insights applicable to nanomaterials, particle-reinforced composites, and biocompatible polymers.

Science & Technology - Other Topics↗
DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗
DOE OSTI2025

Electron Inversion and Tunneling at Silicon Thermal Oxide Interfaces for Solar-Driven Molecular Catalysis to Syngas

Semiconductor photoelectrodes are regularly coupled to solid-state heterogeneous catalysts to perform solar-driven reduction of CO 2 . Less frequently, molecular catalysts are employed to better control the reactivity toward desired products, yet the development of robust semiconductor/molecule interfaces has proven challenging. Here, we demonstrate that a 2–3 nm thermal oxide layer on Si exhibits stability in aqueous solution, high photovoltage, and a photocurrent density of ∼10 mA/cm 2 for the solar-driven photoelectrochemical reduction of a homogeneous molecular catalyst, producing syngas with an ∼2:1 H 2 to CO ratio. Because of a low defect density, the oxide interface forms an electron inversion layer with metal-like electron density at cathodic potentials. This inversion layer facilitates electron transfer to redox-active molecules via tunneling even if the molecule’s reduction potential is beyond the semiconductor’s conduction band edge. Using an electrolyte solution composed of a homogeneous cobalt bis(terpyridine) catalyst in a water/organic solvent mixture, stable photoelectrochemistry was observed under 1-sun illumination, exhibiting an ∼30% Faradaic efficiency for CO that was similar to a glassy carbon electrode under comparable conditions. Furthermore, the results demonstrate that an ultrathin thermal oxide interface is a robust platform for development of aqueous-stable, molecule-driven photoelectrocatalysis.

Catalysts↗
DOE OSTI2024

Benchmarking of ENDF/B-VIII.1 Thermal Scattering Library for Hydrogen

In correctly characterizing the energy and momentum transfer at thermal and cold energies between neutrons and its interacting medium, thermal scattering libraries, which details energy states due to the intra- and inter-molecular bond effects for the medium materials, are applied in place of free-gas cross section libraries in a particle transport simulation code. They are essential to the neutron performance of a neutron facility like SNS, where thermalized neutrons from 20 K liquid hydrogen and ambient (~300 K) water are transported to the beamlines for neutron scattering experiments in studying materials. Recently, a new version of thermal scattering library for parahydrogen and orthohydrogen at 14-20 K was developed and to be released in ENDF/B-VIII.1. It is, therefore, important to benchmark its impacts on the prediction of moderator performance due to the updates in the thermal scattering library. In this study, the recent ENDF/B-VIII.1 thermal scattering library was compared to the current ENDF/B-VII.1 one in the neutron performance calculations for the decoupled and coupled hydrogen moderators at SNS under theorized and real working conditions. In addition, the predictions using both thermal scattering libraries were benchmarked to the measurements of moderator performance. The consistency between the libraries was observed mostly for parahydrogen and the difference in orthohydrogen at cold neutron energies was noted.

Lu, Wei [Oak Ridge National Laboratory (ORNL), Oak↗
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.↗