Search NASA⌕ Search

SEARCH · Search NASA

Results for “mass and heat transfer”

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

Thermally Stable Co@C3N4 Single-Atom Catalysts for CO Oxidation: Atomic-Level Insights into Structure and Activity

Single Co atoms supported on C3N4 (Co@C3N4) have demonstrated high activity and selectivity in photocatalysis. However, the investigation of structure–function relationships and reaction mechanisms under photocatalytic conditions is very challenging due to the complex conditions of light absorption, charge transfer, and catalysis. In this study, we employed thermal CO oxidation as a prototypical probe reaction to benchmark the intrinsic catalytic performance and track the active-site evolution of Co@C3N4. Single Co atoms were identified and shown to be the catalytically active sites for CO oxidation based on control experiments and isotope-labeling experiments. The Co sites remained atomically dispersed before, during, and after the reaction with temperatures up to 400 °C, as established by in situ X-ray absorption fine structure (XAFS) combined with density functional theory (DFT), FDMNES simulations, and dynamic-time-warping (DTW)-assisted X-ray absorption near edge structure (XANES) matching. Together with theoretical calculations, the integrated analysis reveals a stable coordination environment under reaction conditions, which correlates with sustained activity, establishing Co@C3N4 single-atom catalysts as thermally stable CO oxidation catalysts. Beyond these findings, the current study provides a workflow for unambiguously assigning active sites in Co@C3N4 for thermal CO oxidation. This workflow will aid the understanding of their behavior in photocatalysis in the future, where light-driven dynamics obscure direct structure–function links. Notably, this study provides fundamental insights for the rational design of robust single-atom catalysts and a foundation for the broader application of Co@C3N4 catalysts in oxidation reactions.

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

Linking international technical specifications for acoustic characterization of marine energy converter sounds with environmental compliance criteria

As new ocean energy technologies emerge and are deployed for testing and operations, sound emissions are a potential concern for environmental effects to marine life. Consistent acoustic measurement and data analysis methods can help promote comparisons of technologies and transferability between project sites. In 2022, acoustic emissions from a prototype scale wave energy converter (WEC) were characterized for a range of environmental conditions and power generation states in the coastal waters off southern California using a set of international technical specifications. Results from the international technical specification analyses were applied to United States regulatory threshold criteria for acoustic impacts to marine mammals and examined in the context of European underwater noise monitoring guidelines. Weighted 24 hour cumulative sound exposure levels SEL24h calculated from the highest power generation state WEC sound pressure levels were often more than 20 dB below threshold criteria for temporary threshold shifts in five relevant marine mammal hearing groups. Following European Union recommendations for analyses and reporting, WEC sound characterization in third octave bands centered at 63 Hz and 125 Hz show clear spatial decay of WEC-generated noise, with more pronounced attenuation at 63 Hz, and a less marked but still detectable gradient at 125 Hz, collectively suggesting a relatively confined acoustic footprint under the observed conditions. The value of the international technical specification approach is highlighted by the isolation of WEC sounds from the surrounding soundscape. This allows for a robust characterization of acoustic emissions through a range of device power generation and sea states. Furthermore, in threshold-based regulatory contexts like the U.S., this facilitates direct evaluation of source contributions, while in broader monitoring frameworks used in the E.U. it provides a reproducible foundation for assessing the contribution of emerging ocean energy technologies to the underwater acoustic environment.

Haxel, Joseph H. (ORCID:0000000273864761)↗
DOE OSTI2023

Thermal Disorder‐Induced Strain and Carrier Localization Activate Reverse Halide Segregation

The reversal of halide ions is studied under various conditions. However, the underlying mechanism of heat-induced reversal remains unclear. This work finds that dynamic disorder-induced localization of self-trapped polarons and thermal disorder-induced strain (TDIS) can be co-acting drivers of reverse segregation. Localization of polarons results in an order of magnitude decrease in excess carrier density (polaron population), causing a reduced impact of the light-induced strain (LIS – responsible for segregation) on the perovskite framework. Meanwhile, exposing the lattice to TDIS exceeding the LIS can eliminate the photoexcitation-induced strain gradient, as thermal fluctuations of the lattice can mask the LIS strain. Under continuous 0.1 W cm -2 illumination (upon segregation), the strain disorder is estimated to be 0.14%, while at 80 °C under dark conditions, the strain is 0.23%. However, in situ heating of the segregated film to 80 °C under continuous illumination (upon reversal) increases the total strain disorder to 0.25%, where TDIS is likely to have a dominant contribution. Therefore, the contribution of entropy to the system's free energy is likely to dominate, respectively. Various temperature-dependent in situ measurements and simulations further support the results. These findings highlight the importance of strain homogenization for designing stable perovskites under real-world operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polycarbonate‐Based Solid‐Polymer Electrolytes for Solid‐State Sodium Batteries

Solid-polymer electrolytes comprised of polypropylene carbonate (PPC) and varied sodium bis(fluorosulfonyl)imide (NaFSI) salt concentrations are investigated for implementation as a conductive solid polymer electrolyte into solid-state cathode composites utilizing a sodium-layered oxide active material. The ionic conductivity generally increases with NaFSI salt content, reaching ≈1 mS cm −1 at 80 °C at the highest salt concentration (PPC:NaFSI = 0.5:1). Through an all-in-one slurry casting method, Na 2/3 Ni 1/3 Mn 2/3 O 2 cathode composites are fabricated in which the dispersed PPC electrolyte acts as the primary binder. Enabled by a bilayer polymer electrolyte system, cycling performance with the PPC cathode electrolyte is optimized with respect to salt concentration and anode material. The best cyclability is achieved with a moderate salt concentration electrolyte (PPC:NaFSI = 5:1), showcasing an initial capacity of 83 mA h g −1 with a remarkable 80% capacity retention after 150 cycles at C/5 rate and 60 °C. The superior performance of the lower salt concentration electrolyte is attributed to better electrochemical stability, as confirmed by linear sweep voltammetry and online electrochemical mass spectrometry measurements. In conclusion, these results underscore the potential of carbonate-based polymer electrolytes and the importance of balancing electrolyte conductivity and stability in cell design.

25 ENERGY STORAGE↗
DOE OSTI2025

Decoupling the capacity fade contributions in polymer electrolyte-based high-voltage solid-state batteries

Polymer electrolyte (PE)-based solid-state batteries (PE-SSBs) made with high-voltage cathodes are known to suffer from severe capacity fade, stemming primarily from the poor oxidative stability of most PEs under high-voltage cycling conditions. PEs also suffer from greater ion-transport limitations compared to liquid or solid electrolytes. However, often, these limitations are collectively stated to be responsible for the observed capacity fade, and it is challenging to decouple the contributions of different factors. Herein, a tunable cell fabrication platform was developed to systematically investigate and decouple the two primary capacity fade drivers (cell impedance growth and kinetic limitations), while keeping the other cell parameters constant. Three PE types with distinct transport characteristics were compared. By utilizing a voltage profile analysis method, the contribution of the cell's internal impedance growth was quantitatively decoupled from the kinetic limitations stemming from the high concentration gradient in the polymer catholyte and slow charge transfer reactions. We demonstrate that the high interfacial impedance did not necessarily correlate with the high capacity fade rate. Kinetic limitations that are not reflected by impedance measurements can play a dominant role in causing cumulative capacity decay.

Ock, Ji-young [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2023

Speciation of Transition Metal Dissolution in Electrolyte from Common Cathode Materials

Significant capacity loss has been observed across extended cycling of lithium-ion batteries cycled to high potential. One of the sources of capacity fade is transition metal dissolution from the cathode active material, ion migration through the electrolyte, and deposition on the solid-electrolyte interphase on the anode. While much research has been conducted on the oxidation state of the transition metal in the cathode active material or deposited on the anode, there have been limited investigations of the oxidation state of the transition metal ions dissolved in the electrolyte. Here, in this work, X-ray absorption spectroscopy has been performed on electrolytes extracted from cells built with four different cathode active materials (LiMn 2 O 4 (LMO), LiNi 0.5 Mn 1.5 O 4 (LNMO), LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), and (x Li 2 MnO 3 *(1-x) LiNi a Mn b Co c O 2 , with a+b+c=1) (LMRNMC)) that were cycled at either high or standard potentials to determine the oxidation state of Mn and Ni in solution. Inductively coupled plasma-mass spectrometry has been performed on the anodes from these cells to determine the concentration of deposited transition metal ions. While transition metal ions were found dissolved in all electrolytes, the oxidation state(s) of Mn and Ni were determined to be dependent on the cathode material and independent of cycling potential.

25 ENERGY STORAGE↗
DOE OSTI2026

How Proton Incorporation Reshapes Lattice Dynamics In BaSnO 3 ‐Type Proton Conductors

Proton conduction in acceptor-doped perovskites is fundamentally a vibronic process: mobile H + and D + do not move independently, but dynamically co-vibrate with the surrounding oxygen–metal framework. Direct experimental evidence for this behavior is presented using in situ 119⁢ 𝑆⁢𝑛 nuclear resonance vibrational spectroscopy (NRVS) on hydrated, deuterated, and dry 𝐵⁢𝑎⁢𝑆⁢𝑛 1−𝑥⁢ 𝑌 𝑥 ⁢𝑂 3−𝛿 . Hydration induces systematic redistributions in the Sn-projected phonon density of states (PDOS), including an upshift of the first spectral moment by about 0.4 meV, indicating a stiffening of the extended Sn–O network. H/D isotopic substitution leaves the Sn-projected PDOS largely unchanged, with only subtle isotope-dependent spectral reweighting, demonstrating that protonic degrees of freedom are not localized oscillators but are embedded in collective lattice modes. These results are rationalized using a classical coupled proton–phonon oscillator model that links the observed PDOS variations to changes in effective force constants and vibrational mass terms. The model captures how H + and D + participate in cooperative lattice dynamics rather than forming isolated OH/OD entities. Overall, NRVS probes proton–lattice coupling in ceramic proton conductors and quantitatively describes how protonic defects modulate host lattice dynamics to enable phonon-assisted long-range proton transport.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The ABCs of phase retrieval: Connecting the acronyms of scanning transmission electron microscopy

High-resolution scanning transmission electron microscopy (S/TEM) is an indispensable tool for characterizing the structure and properties of materials down to the atomic scale. Conventional S/TEM imaging, however, is limited by the phase problem, whereby the phase of the electron exit wave is lost upon detection. Recent advances in diffractive imaging and 4D-STEM have enabled a range of phase-retrieval techniques that computationally reconstruct the missing information encoded in the phase of the transmission function. These approaches offer improved dose efficiency and enhanced sensitivity to weakly scattering signals, extending quantitative imaging to beam-sensitive materials composed of light elements. In this work, we introduce the phase problem in electron microscopy and survey the diverse landscape of phase-retrieval techniques used in the field. Despite their many acronyms and algorithmic variations, these techniques share a common physical and mathematical foundation. We present a unified framework that connects these seemingly distinct methods, from parallax imaging and tilt-corrected bright-field (tcBF-STEM), to aberration-corrected bright-field (acBF-STEM), optimum bright-field (OBF-STEM) and single-sideband (SSB) ptychography, as well as first-moment integrated center of mass techniques (iCOM) and iterative ptychographic algorithms. Based on these insights, we discuss the opportunities and practical limitations of applying these methods across different materials systems, detector designs, and microscope configurations.Graphical abstractRepresentative electron microscopy configurations used for phase retrieval and diffractive imaging in S/TEM: (a) Zernike phase-contrast transmission electron microscopy (TEM), (b) small-convergence-angle four-dimensional scanning transmission electron microscopy (4D-STEM) for nanobeam-based phase reconstruction methods, and (c) large-convergence-angle 4D-STEM for ptychographic and related diffractive imaging techniques reviewed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Elucidating the impact of stress states on grain boundary passivation in Ni-5Cr using the Rhines pack method

The oxidation behavior of a Ni–5Cr (at.%) alloy was evaluated at 420 ° C using the Rhines pack method and simultaneous tensile and compressive stress states via a miniature four-point bending fixture. At this moderate temperature, grain boundaries dominate mass transport and the resulting oxidation response. Oxidation produced approximately 1μ⁢m-wide protective Cr 2 O 3 films capping some grain boundaries and penetrative, intergranular Cr-rich oxides at other grain boundaries. Externally applied tensile and compressive stress during oxidation increased the prevalence of Cr 2 O 3 cap formation compared to no applied stress, with tensile stress resulting in more Cr 2 O 3 caps than compressive stress. The observed oxide cap morphology was similar across all test conditions. Regions under compressive stress showed an order of magnitude greater Cr depletion depth along the grain boundary. Chromium nitrides (CrN), likely from N contamination of the Rhines pack cell, were observed both at the oxide–metal interface and intergranularly within the alloy. Collectively these results demonstrate that applied stress promotes localized protective oxide cap formation over grain boundaries, with compressive stress additionally promoting deeper Cr depletion. Furthermore the experimental approach helped separate out the effects of stress on local oxide formation and grain boundary passivation.

4 point bend↗
DOE OSTI2026

Limited proteolysis of pea protein to promote aggregation/gelation: role of enzyme concentration and molecular characteristics

Limited proteolysis promotes aggregation/gelation of plant proteins, but its relationship with the molecular characteristics of hydrolysates remains poorly understood. In here, response surface methodology was employed to evaluate the effects of protein concentration, Alcalase level, and hydrolysis time on the storage modulus (G′) of pea protein hydrolysates (PPH). 7.5% protein, 1.65% Alcalase, and a 6-min reaction time resulted in the highest relative increase in G′ upon heating PPH, whereas both higher (3.55%) and lower (0.3%) Alcalase levels led to less increment. Increasing the Alcalase level increased the degree of hydrolysis while decreasing the molecular weight, surface hydrophobicity, and total sulfhydryl and disulfide bond contents of PPH. Confocal laser scanning microscopy and X-ray scattering analyses found PPH produced with 1.65% Alcalase formed interconnected aggregates with dual structural domains, but not at other Alcalase conditions. These results demonstrate the role of controlled proteolysis in governing the aggregation/gelation behavior of PPH.

X-ray scattering↗
DOE OSTI2026

Photovoltage behaviour of p-Sb 2 S 3 photocathodes for hydrogen evolution: effect of n-In 2 S 3 passivation layers

The 1.76 eV band gap of antimony(iii) sulphide (Sb 2 S 3 ) makes this semiconductor material a promising light absorber for photoelectrochemical water splitting, but scalable fabrication approaches to efficient devices are still lacking. Here we show that compact Sb 2 S 3 films on FTO can be obtained by electrochemical growth from aqueous colloidal sulphur and antimony trichloride solutions, followed by mild annealing. These films can be converted into hydrogen evolution photocathodes after coating with In 2 S 3 passivation layers and the addition of Pt proton reduction co-catalysts. For the first time, vibrating Kelvin probe surface photovoltage (VKP-SPV) spectroscopy is used to observe the carrier dynamics in such photoelectrodes. While the bare Sb 2 S 3 films suffer from high surface recombination rates and poor electron extraction, the In 2 S 3 overlayer is found to raise the photovoltage and cathodic photocurrent density, due to passivation of surface defects and formation of a p–n heterojunction. In thick In 2 S 3 films, these benefits are offset by shading and slow electron transfer. Also, we find that O 2 strongly affects the band bending in the Sb 2 S 3 –air and In 2 S 3 –air junctions and their photovoltage. The optimised devices evolve H 2 at 77.5% Faradaic efficiency and with 0.084% applied bias photon-to-current efficiency (ABPE) at 0.12 V vs. RHE. The low ABPE value is attributed to Sb 2 S 3 sub-bandgap defects visible in SPV spectra, the random orientation of Sb 2 S 3 crystallites in the films, which inhibits charge transport, the absence of crystal facets of Sb 2 S 3 , and a detrimental Schottky junction at the FTO|Sb 2 S 3 interface.

de Araújo, Moisés A. [University of California, Da↗
DOE OSTI2026

Induction Melter Processing Alternatives for High Level Radioactive Wastes – 26198

This work investigates the potential to vitrify nuclear fuel directly, as well as the vitrification potential associated with experimental dissolver solutions. Greater understanding of the exothermicity is needed to quantify processing risk, especially with respect to potential phase changes and associated explosion hazards present in some systems. Thermal analysis was carried out on various simulants to elucidate the exothermic reaction potential from solid metal dissolution in glass and from the drying of alternative process dissolver solutions. Results from experimental testing of simulants to demonstrate vitrification potential and compatibility indicate that alternative dissolver flowsheets suppress the heat released during processing and that common silicate- and phosphate- based glass systems are potential candidates for direct vitrification. Direct vitrification (conversion) of fuel simulants was assessed using laboratory scale glass melts to obtain qualitative information on dissolution rates and waste loadings. Initial tests were successful to dissolve and incorporate metal directly into glass, although the kinetics and limits of dissolution and incorporation into glass are not fully understood.

Amoroso, Jake [Savannah River National Laboratory ↗
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

Vapor-phase pillarization of MXenes for engineering hierarchical interlayer porosity

MXenes, a family of two-dimensional (2D) multilamellar materials, possess excellent thermal and electronic properties for a range of applications. Their use in heterogeneous catalysis, however, is limited by the low surface area resulting from stacked layers. Pillarization with inorganic oxides can create more open, mesoporous MXene structures, improving accessibility for guest species to diffuse, reside or react in the space between 2D layers. A previous liquid-phase pillarization method, however, involves excessive use of solvent-based precursors and multiple processing steps. Here, we report a vapor-phase pillarization (VPP) strategy to introduce pillars, exemplified by silica pillars, with high pillar precursor usage efficiency and a simplified processing workflow. The resulting silica-pillared mesoporous MXene exhibits significantly increased surface area and porosity. These textural properties can be easily tuned by the VPP synthesis conditions. When applied as a ruthenium (Ru) catalyst support for the hydrogenolysis of low-density polyethylene (LDPE), the silica-pillared MXene enabled high Ru dispersion and catalytic activity. This study highlights the potential of the VPP method for engineering mesoporous, 2D MXene materials and demonstrates the effectiveness of mesoporous MXene as a catalyst support in overcoming mass transport and active-site accessibility challenges in heterogeneous catalysis involving bulky substances, such as plastics upcycling.

Luo, Song [University of Delaware, Newark, DE (Uni↗
DOE OSTI2023

Oxidation Chemistry of Bicarbonate and Peroxybicarbonate: Implications for Carbonate Management in Energy Storage

Carbonate formation presents a major challenge to energy storage applications based on low-temperature CO 2 electrolysis and recyclable metal–air batteries. While direct electrochemical oxidation of (bi)carbonate represents a straightforward route for carbonate management, knowledge of the feasibility and mechanisms of direct oxidation is presently lacking. Herein, we report the isolation and characterization of the bis(triphenylphosphine)iminium salts of bicarbonate and peroxybicarbonate, thus enabling the examination of their oxidation chemistry. Infrared spectroelectrochemistry combined with time-resolved infrared spectroscopy reveals that the photoinduced oxidation of HCO 3 – by an Ir(III) photoreagent results in the generation of the short-lived bicarbonate radical in less than 50 ns. The highly acidic bicarbonate radical undergoes proton transfer with HCO 3 – to furnish the carbonate radical anion and H 2 CO 3 , leading to the eventual release of CO 2 and H 2 O, thus accounting for the appearance of H 2 O and CO 2 in both electrochemical and photochemical oxidation experiments. Here, the back reaction of the carbonate radical subsequently oxidizes the Ir(II) photoreagent, leading to carbonate. In the absence of this back reaction, dimerization of the carbonate radical provides entry into peroxybicarbonate, which we show undergoes facile oxidation to O 2 and CO 2 . Together, the results reported identify tangible pathways for the design of catalysts for the management of carbonate in energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗