SEARCH · Search NASA

Results for “CHARGE DISTRIBUTION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

562 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗

DOE OSTI2026

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗

NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗

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↗

NASA NTRS1995

Analysis of Gas Absorption to A Thin Liquid Film in the Presence of A Zero-Order Chemical Reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

S Rajagopalan↗

DOE OSTI2025

Direct observation of distinct bulk and edge nonequilibrium spin accumulation in ultrathin MoTe 2

Low-symmetry two-dimensional (2D) topological materials such as MoTe 2 host efficient charge-to-spin conversion (CSC) mechanisms that can be harnessed for novel electronic and spintronic devices. However, the nature of the various CSC mechanisms and their correlation with underlying crystal symmetries remain unsettled. In this work, we use local spin-sensitive electrochemical potential measurements to directly probe the spatially dependent nonequilibrium spin accumulation in MoTe 2 flakes down to four atomic layers. We are able to clearly disentangle contributions originating from the spin Hall and Rashba-Edelstein effects and uncover an abundance of unconventional spin polarizations that develop uniquely in the sample bulk and edges with decreasing thickness. Using ab-initio calculations, we construct a unified understanding of all the observed CSC components in relation to the material dimensionality and stacking arrangement. Our findings not only illuminate previous CSC results on MoTe 2 but also have important ramifications for future devices that can exploit the local and layer-dependent spin properties of this 2D topological material.

Science & Technology - Other Topics↗

DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

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↗

NASA NTRS1994

Bipolar Ag-Zn battery

The silver-zinc (AgZn) battery system has been unique in its ability to safely satisfy high power demand applications with low mass and volume. However, a new generation of defense, aerospace, and commercial applications will impose even higher power demands. These new power demands can be satisfied by the development of a bipolar battery design. In this configuration the power consuming, interelectrode current conductors are eliminated while the current is then conducted via the large cross-section electrode substrate. Negative and positive active materials are applied to opposite sides of a solid silver foil substrate. In addition to reducing the weight and volume required for a specified power level, the output voltage performance is also improved as follows. Reduced weight through: elimination of the plastic cell container; elimination of plate leads and intercell connector; and elimination of internal plate current collector. Increased voltage through: elimination of resistance of current collector; elimination of resistance of plate lead; and elimination of resistance of intercell connector. EPI worked previously on development of a secondary bipolar silver zinc battery. This development demonstrated the electrical capability of the system and manufacturing techniques. One difficulty with this development was mechanical problems with the seals. However, recent improvements in plastics and adhesives should eliminate the major problem of maintaining a seal around the periphery of the bipolar module. The seal problem is not as significant for a primary battery application or for a requirement for only a few discharge cycles. A second difficulty encountered was with activation (introducing electrolyte into the cell) and with venting gas from the cell without loss of electrolyte. During previous work, the following projections for energy density were made from test data for a high power system which demonstrated in excess of 50 discharge/charge cycles. Projected system power = 100 kilowatts; discharge time = 30 seconds; discharge current density = 1.75 amps/sq in.; system weight = 86 lbs (9.7 WH/lb); and system volume = 1071 cu. in. (.78 WH/cu. in.). EPI is currently working on a development program to produce a bipolar silver-zinc battery design for NASA. The potential application would be to power electromechanical actuators for space launch vehicles.

Giltner, L. John↗

NASA NTRSDate not supplied

BEAST: Expanding Sustainable Data Infrastructure for High-Enthalpy Facilities

Reproducible, data-driven thermal protection system (TPS) research requires that experimental records from high-enthalpy testing be consistently structured, traceable, and accessible across campaigns and institutions. In practice, however, arcjet and plasma facilities data remain largely fragmented: raw diagnostics are stored in ad hoc formats, material sample histories are disconnected from test conditions, and metadata standards are absent, precluding systematic cross-campaign analysis and long-term reuse. BEAST (Backend for Experiment Analysis, Storage, and Traceability) is an open-source, web-based platform that addresses these limitations by providing a unified, queryable infrastructure for high-enthalpy ground-test data [1]. First presented at the 15th Ablation Workshop [2], BEAST has since undergone significant development. The platform ingests and structures multi-channel time-series diagnostics, facility configurations, and material property records within a common provenance model, ensuring end-to-end traceability from raw sensor acquisition to reduced experimental quantities. A versioned material library links specimen identity and processing history to the specific runs in which each sample was tested. An integrated modeling workbench enables training and evaluation of regression models directly on archived experimental data, supporting condition interpolation and the construction of empirical material response databases. Beyond its original deployment at NASA Ames Research Center, BEAST has been designed to be facility-agnostic, with ongoing efforts to extend its adoption to other facilities. Its modular architecture accommodates heterogeneous diagnostic setups and facility types, and its future open-source distribution allows institutions to build on a common data standard rather than maintaining isolated, bespoke solutions. BEAST is further integrated within a broader ecosystem of companion tools: arcjetCV [3] extracts recession rates and shock standoff distances from high-speed video using computer vision, and miniSTARscan [4] provides sub-minute, portable photogrammetric surface reconstruction of test articles before and after exposure. All tools share a common data schema, enabling seamless ingestion of surface geometry, imagery, and time-series data into a single, coherent experimental record.

Database↗

NASA NTRSDate not supplied

BEAST: Expanding Sustainable Data Infrastructure for High-Enthalpy Facilities

Database↗

NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

DOE OSTI2026

Reduction of baseplate distortion during directed energy deposition using compliant features

Distortion in additive manufacturing (AM) remains a barrier to its adoption in precision industries. Baseplate warpage is one such issue which compromises the feasibility of post-process precision machining. The restriction to thermal contraction of the deposited part, imposed by the baseplate, generates bending moments, that in turn causes warpage. A novel distortion mitigation strategy using baseplates with integrated compliant features, which enables thermal contraction of the build is presented. Six unique design concepts are evaluated, including a solid reference, using two deposition geometries. A laser, hot-wire, directed energy deposition (DED) process is used to deposit a symmetric cylindrical part (C-part) and a T-shaped asymmetric part (T-part). Flatness deviation of baseplates is measured using structured light 3D scanning. Measurements reveal a reduction in net flatness deviation of 58.8% for the C-part and 40.9% for the T-part, compared to the solid reference. While most designs yielded reductions exceeding 30% and 20% for the C- and T-parts, respectively, one configuration resulted in increased deviation. Finite element (FE) simulations are performed to elucidate the underlying mechanisms affecting distortion of compliant baseplates during DED. Despite variations between predictions and measurements, agreement in the general trend is observed. It also revealed that initial flatness errors in baseplates significantly affect its deviation during deposition. Predictions indicate that compliant features significantly affect the thermal distribution as well as the evolution of flatness deviation in the baseplate during deposition. Notably, one design exhibited a reduction in distortion during cooling, following its initial increase during deposition. FE predictions show a maximum reduction of 60.9% and 38.8% in net flatness deviation for the C- and T-parts, respectively. The performance of compliant baseplates is found to be governed by both its the thermal and mechanical characteristics, which are crucial factors to be considered during design.

Mathews, Ritin [ORNL] (ORCID:0000000301440828)↗

NASA NTRS1997

Analysis of Impact Induced Damage and its Effect on Structural Integrity of Space Flight Composite Overwrapped Pressure Vessels

The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.

Michael P Wnuk↗

DOE OSTI2023

The brighter-fatter effect in the JWST MIRI Si:As IBC detectors: I. Observations, impact on science, and modeling

Context.The Mid-Infrared Instrument (MIRI) on board theJames WebbSpace Telescope (JWST) uses three Si:As impurity band conduction (IBC) detector arrays. The output voltage level of each MIRI detector pixel is digitally recorded by sampling up the ramp. For uniform or low-contrast illumination, the pixel ramps become nonlinear in a predictable way, but in areas of high contrast, the nonlinearity curve becomes much more complex. The origin of the effect is poorly understood and currently not calibrated out of the data. Aims.We provide observational evidence of the brighter-fatter effect (BFE) in MIRI conventional and high-contrast coronagraphic imaging, low-resolution spectroscopy, and medium-resolution spectroscopy data, and we investigate the physical mechanism that gives rise to the effect on the MIRI detector pixel raw voltage integration ramps. Methods.We used public data from the JWST/MIRI commissioning and Cycle 1 phase. We also developed a numerical electrostatic model of the MIRI detectors using a modified version of the publicPoisson_CCDcode. Results.We find that the physical mechanism behind the BFE manifesting in MIRI data is fundamentally different to that of charge-coupled devices and photodiode arrays such as the Hawaii-XRG near-infrared detectors used by the NIRISS, NIRCam, and NIRSpec instruments on board JWST. Observationally, the BFE makes the JWST MIRI data yield 10–25% larger point sources and spectral line profiles as a function of the relative level of de-biasing of neighboring detector pixels. This broadening impacts the MIRI absolute flux calibration, time-series observations of faint companions, and point spread function modeling and subtraction. We also find that the intra-pixel 2D profile of the shrinking Si:As IBC detector depletion region directly impacts the accuracy of the pixel ramp nonlinearity calibration model.

Astronomy & Astrophysics↗

DOE OSTI2023

Spatiotemporal Studies of Soluble Inorganic Nanostructures with X‐rays and Neutrons

This Review addresses the use of X-ray and neutron scattering as well as X-ray absorption to describe how inorganic nanostructured materials assemble, evolve, and function in solution. We first provide an overview of techniques and instrumentation (both large user facilities and benchtop). We review recent studies of soluble inorganic nanostructure assembly, covering the disciplines of materials synthesis, processes in nature, nuclear materials, and the widely applicable fundamental processes of hydrophobic interactions and ion pairing. Reviewed studies cover size regimes and length scales ranging from sub-Ångström (coordination chemistry and ion pairing) to several nanometers (molecular clusters, i.e. polyoxometalates, polyoxocations, and metal-organic polyhedra), to the mesoscale (supramolecular assembly processes). Reviewed studies predominantly exploit 1) SAXS/WAXS/SANS (small- and wide-angle X-ray or neutron scattering), 2) PDF (pair-distribution function analysis of X-ray total scattering), and 3) XANES and EXAFS (X-ray absorption near-edge structure and extended X-ray absorption fine structure, respectively). While the scattering techniques provide structural information, X-ray absorption yields the oxidation state in addition to the local coordination. Our goal for this Review is to provide information and inspiration for the inorganic/materials science communities that may benefit from elucidating the role of solution speciation in natural and synthetic processes.

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.↗