Search NASASearch

SEARCH · Search NASA

Results for “21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS”

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.

3,200 records · Page 43

Hierarchical low Pt-loading “core-shell” electrocatalysts for the oxygen reduction reaction in fuel cells

The sluggish kinetics of the oxygen reduction reaction (ORR) hinder cost-effective polymer electrolyte fuel cells (PEFCs), which rely on scarce, expensive platinum-based electrocatalysts (ECs). Here, we present a novel synthesis method for ORR ECs achieving exceptional platinum utilization. The design features a hierarchical “multi-carbon” support comprising carbon nanoparticles interacting with graphene nanoplatelets as the “core”, encapsulated by a porous carbon nitride (CN) “shell”. This configuration promotes strong core/shell interactions and a bimodal active site distribution, consisting of chemically dispersed Pt and Ni single-atom complexes and PtNix alloy nanoclusters embedded in the CN shell. These advantages enable high activity and durability, achieving an ORR activity of 1.6 A mgPt−1 at 0.9 V vs. RHE-an order of magnitude higher than Pt/C (0.17 A mgPt−1). A proof-of-concept PEFC demonstrates a specific power of 12.0 kW gPt−1 at 0.60 V. This approach offers a significant step toward more efficient and sustainable PEFC technologies.

Pagot, Gioele [University of Padova, Italy]

Treatment of Lagoon Dairy Manure Wastewater via Iron Electrocoagulation, Microfiltration, and Adsorption

Dairy manure wastewater generated by flushing barn cow waste contains nutrients, pathogens, and organic and inorganic contaminants. This study utilized a process consisting of iron electrocoagulation (Fe-EC), microfiltration (MF), and activated carbon (AC) adsorption to treat farm wastewater and explore the reclamation of clean water for irrigation and livestock consumption. Significant removal (>99.9%) of chemical oxygen demand (COD), total organic carbon (TOC), phosphorus (P), turbidity, and microorganisms, as well as ions including magnesium, calcium, sulfur, and silica was achieved by the combined EC-MF-AC process. Specifically, a charge loading of ∼37,500 C/L in a continuous-flow EC configuration, followed by MF, achieved more than 95% removal of TOC and COD. Characterization of produced flocs and foam via scanning-electron microscopy with energy-dispersal spectroscopy and Fourier transform infrared spectroscopy confirmed the removal of ions, including calcium, sulfur, and silica. A key finding was the electrocatalytic conversion of nitrogen species to ammonia gas through the intermediate reduction of nitrate/nitrite, which led to ∼60% total nitrogen (TN) removal. AC treatment further improved TN removal to ∼70%. The Fe-EC process also eradicated >99.9% of bacteria. Preliminary process cost assessment, based on recycled materials for EC electrodes, showed significant cost savings (∼2 times) compared to commercial electrodes.

Dutta, Swapnamoy [ORNL]

From Rules to Reasoning: A Survey of Large Language Model-Based Approaches to Scientific Hypothesis and Idea Generation

Scientific hypothesis generation represents a fundamental challenge in contemporary research due to exponentially expanding literature volumes and increasing disciplinary specialization. Large language models (LLMs) have emerged as transformative tools for automated scientific discovery, moving beyond traditional rule-based and literature-mining approaches. Four paradigmatic approaches define current LLM-driven hypothesis generation: direct prompting and fine-tuning methods, knowledge-enhanced frameworks integrating retrieval-augmented generation (RAG), multi-agent collaborative systems simulating research teams, and reasoning-focused approaches implementing cognitive architectures. Domain-specific applications demonstrate statistical equivalence to human expert performance in social psychology, experimental validation in biomedical research, and near-expert quality in astronomy. Evaluation methodologies encompass human expert assessment, LLM-as-judge frameworks, and comprehensive benchmarking systems. Technical challenges include hallucination management, knowledge integration limitations, and balancing novelty with feasibility. Future directions emphasize hybrid neural-symbolic architectures and sophisticated human-AI collaboration models for responsible scientific discovery acceleration.

AI-driven discovery

Physics-Based Modeling and Simulation of Emerging Battery Technologies for Aerospace

Recently there is a growing interest in the aviation sector to reduce air and noise pollution. Electrochemical energy storage devices such as batteries coupled with a distributed electric propulsion system can reduce noise concerns as well as emissions and allow the concept of Urban Air Mobility to come to fruition. The battery performance needs to improve considerably compared from current state-of-art Li-ion battery (about 200Wh/Kg) to realize all-electric passenger jet for short flights (up to 690 miles). NASA is exploring lithium-oxygen battery chemistry to power hybrid (battery powered electrical system) and all-electric aircraft for short distance and long-distance flights. Li-O2 is one of the advanced Li-ion technologies that promise to provide specific energy of more than 750Wh/Kg. For this presentation, we present our work on improving power density of Li-O2 batteries through the use of multiphysics simulations. Next, a path is outlined to port these model to simulate performance for a new battery chemistry for space application, Li-CO2. Li-CO2 uses carbon dioxide as the active material instead of oxygen. Although this technology is in its early development, the offers two benefits: it can be used as a CO2 scrubber, as oxygen is one of the by-products on charging, and as a backup or a standalone battery for various Mars or Venus missions, where the carbon dioxide content in the atmosphere is high and need battery to operate at higher temperatures.

Mehta, Mohit

The Origin of Palimpsests

Palimpsests are the most enigmatic impact features observed on icy worlds of the Solar System. Palimpsests are large impact features which are bright, old, mostly flat, and roughly circular. They range from >150 to ~350 km in diameter and sometimes exhibit recognizable associated radial secondary ejecta crater strings beyond their rims. Curiously, among all the Solar System's icy worlds, palimpsests are exclusively seen on the icy Galilean moons. Palimpsests are obviously different than similarly sized impact features seen on other icy satellites. For example, Evander, a 350 km diameter peak-ring impact basin on Saturn’s moon Dione, has a very low depth/diameter ratio of 0.008, indicating significant viscous relaxation. Nevertheless, it still retains short wavelength topography including a prominent central peak, peak ring, and a continuous rim, none of which are seen at palimpsests. Palimpsests display an approximately circular, and dominating, expanse of undulating plains, which is level with surrounding terrain and shows no recognizable rim. At the center of palimpsests are smooth central plains that are enclosed by a low scarp. Central plains material superposes the surrounding undulating plains and is littered with small blocks reminiscent of “small chaos” on Europa. Palimpsests have yet to be explored using modern impact simulation tools. We have begun to investigate two competing hypotheses: (A) Palimpsests form on targets which have a preexisting near-subsurface layer that is liquid or solid but very weak, and is excavated as proximal ejecta; or (B) They form in a solid ice target underlain by a steep thermal gradient that strongly influences the resulting impact feature, including the production of impact-induced melting at depth that is then a major component of proximal ejecta. Our results will have profound implications for the early thermal history of the icy Galilean satellites. We will test our two hypotheses for the formation of palimpsests on Ganymede (where they are both best observed and preserved) by means of 2D iSALE hydrodynamic simulations of impacts into ice targets that will be structured according to the scenarios for each hypothesis, with select runs in 3D. We have already conducted a pilot study using 2D iSALE that has yielded promising results, which we will report out at the meeting.

palimpsests

Scientific Highlights and Status of the Mars Science Laboratory Mission, Years 9 through 12

We present the scientific highlights and status from years 9-12 of NASA’s Mars Science Laboratory mission, with its Curiosity rover. A thermochemolysis-based sample analysis revealed a variety of organic molecules preserved in clay mineral-rich lake/lake margin deposits of the Glen Torridon region. The stratigraphically higher clay-sulfate mineral transition on Aeolis Mons was found to be associated with a lithological change from lacustrine to dry aeolian depositional environments. Water intermittently returned to the surface after the increase in aridity, including during the deposition of wave ripple strata which indicate a climate capable of sustaining an ice-free lake in the early Hesperian. The overlying Mg sulfate-bearing unit has both poly- and monohydrated Mg sulfates and is also enriched in iron carbonate of a purity consistent with chemical sedimentation and therefore atmospheric sequestration. A channel-incised canyon within Gediz Vallis and distal alluvial deposits recorded numerous fluvial and debris flow events after the deposition, lithification, and erosion of the clay- and sulfate-bearing units. Bright clasts found in Gediz Vallis canyon are composed of native sulfur. The investigation of the modern atmosphere and environment focused on the influence of the crater and mountain topography and captured high-altitude noctilucent and iridescent clouds. Radiation monitoring now spans an 11-year solar cycle. Rover systems and science instruments remain capable of addressing mission scientific objectives. All ten scientific instruments continue to return high-value data despite degradation of the Chemistry and Camera instrument’s high-voltage system, loss of multispectral imaging capability, and loss of active neutron spectroscopy. The mission has mitigated the loss of redundant brakes on two arm joints, additional wheel wear, and highly reduced capability of the redundant flight computer. Efficiency improvements in rover energy utilization and operations staffing have allowed the mission team to maintain a high level of productivity despite declining rover power generation and funding.

Mars

In‐Plane Anisotropy in van der Waals NiTeSe Ternary Alloy

Abstract The anisotropic properties of materials profoundly influence their electronic, magnetic, optical, and mechanical behaviors and are critical for a wide range of applications. In this study, the anisotropic characteristics of Ni‐based van der Waals materials, specifically NiTe 2 and its alloy NiTeSe, utilizing a combination of comprehensive scanning tunneling microscopy (STM), angle‐resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations, are explored. Unlike 1T‐NiTe 2 , which exhibits trigonal in‐plane symmetry, the substitution of Te with Se in NiTe 2 (resulting in the NiTeSe alloy) induces a pronounced in‐plane anisotropy. This anisotropy is clear in the STM topographs, which reveal a distinct linear order of charge distribution. Corroborating these observations, ARPES measurements and DFT calculations reveal an anisotropic Fermi surface centered at the point, which is notably elongated along the k y direction, leading to directional variations in in‐plane carrier velocities. Consequently, the Fermi velocity is highest along the k x direction where the linear charge distribution aligns in real space and is lowest along the k y direction. These findings offer valuable insights into the tunability of anisotropic properties in ternary transition metal dichalcogenide systems, highlighting their potential applications in the development of anisotropic electronic and optoelectronic devices.

Chemistry

Disorder-induced spin-cluster magnetism in a doped kagome spin liquid candidate

The search for new quantum spin liquid materials relies on systems with strong frustration such as spins on an ideal kagome lattice. However, lattice imperfections can have substantial effects which are as yet not well understood. In recent work, the two-dimensional kagome system YCu 3 ⁢(OH) 6 ⁢[(Cl 𝑥 ⁢Br (1−𝑥) ) 3−𝑦 ⁢(OH) 𝑦 ] has emerged as a leading candidate hosting a Dirac spin liquid which appears to survive at least for 𝑥 < 0.4, associated with alternating-bond-hexagon (ABH) disorder. Here in magnetic samples with 𝑥 = 0.58, 𝑦 = 0.1 we report unusual in-plane ferromagnetic canting (FM) of the in-plane antiferromagnet (AFM), with an unusually wide regime of short-ranged order, and propose theoretical models to explain this behavior. First, we show that Kitaev-type exchanges naturally arise on the kagome lattice to second order in the known Dzyaloshinskii-Moriya exchanges, and that these interactions can produce the unusual in-plane FM canting from antichiral AFM. Second, we propose a phenomenological model of weakly FM-canted spin clusters to describe the short-ranged regime and analyze quantum fluctuations in an ABH toy model to show how ABH disorder can stabilize this regime. Here, the combination of experimental observation and theory suggests that kagome-Kitaev interactions and ABH disorder are necessary for describing the magnetic fluctuations in this family of materials, with potential implications for the proposed proximate spin liquid phase.

Seth, Arnab [Georgia Institute of Technology, Atla

Cryogenic Spray Quenching of A Simulated Propellant Storage Tank Wall With Heat Transfer Enhancement By A Thin-Film Coating and Flow Pulsing in Microgravity

Human space exploration to the Moon, Mars, and possibly asteroids is NASA’s biggest challenge for the new millennium. One of the critical elements to this mission is the effective, sufficient, and reliable supply of cryogenic propellant fluids. Future lower-earth-orbiting (LEO) propellant fuel depots and human-carrying orbital transfer spacecraft flying to the moon and Mars will have to utilize the high thrust and high efficiency of liquid cryogenic chemical propulsion or nuclear thermal propulsion. Efficient in-space tank-to-tank propellant transfer (propellant fuel depot to orbital transfer spacecraft) of cryogenic propellants is an enabling technology for the planned Crewed Mars Surface Mission. The transfer of cryogenic propellants in space, however, has yet to be accomplished, solely due to the unavailability of cryogenic quenching heat transfer data during chilldown (quenching) and filling of the propellant receiver tank in reduced gravity and microgravity as liquid propellant cannot be stored in a required liquid state until the tank is quenched down to the liquid temperature. Therefore, highly energy efficient thermal-fluid management breakthrough concepts to conserve and minimize the cryogen consumption during propellant transfer have become the focus of research and engineering development, especially for the deep-space mission to Mars. In this paper, we introduce such concepts and demonstrate their feasibility for cryogenic storage tank chilldown in parabolic flights under a simulated space microgravity condition. In order to maximize the storage tank chilldown efficiency for the least amount of cryogen consumption, the technology adopted included cryogenic spray cooling, Teflon thin-film coating of the simulated tank surface, and spray flow pulsing. The completed flight experiments successfully demonstrated that spray cooling is the most efficient cooling method for the tank chilldown in microgravity. In microgravity, Teflon coating alone can improve the efficiency up to 72% and the efficiency can be improved up to 59% by flow pulsing alone. However, Teflon coating together with flow pulsing was found to substantially enhance the chilldown efficiency in microgravity for up to 113%.

spray

Revealing the evolution of order in materials microstructures using multi-modal computer vision

The development of high-performance materials for microelectronics, energy storage, and extreme environments depends on our ability to describe and direct property-defining microstructural order. Our present understanding is typically derived from laborious manual analysis of imaging and spectroscopy data, which is difficult to scale, challenging to reproduce, and lacks the ability to reveal latent associations needed for mechanistic models. Here, we demonstrate a multi-modal machine learning (ML) approach to describe order from electron microscopy analysis of the complex oxide La 1−x Sr x FeO 3 . We construct a hybrid pipeline based on fully and semi-supervised classification, allowing us to evaluate both the characteristics of each data modality and the value each modality adds to the ensemble. We observe distinct differences in the performance of uni- and multi-modal models, from which we draw general lessons in describing crystal order using computer vision.

36 MATERIALS SCIENCE

Future Projection of Tropical Upper-Tropospheric Troughs and Implications for Tropical Cyclone Activity

Summertime tropical upper-tropospheric troughs (TUTTs) provide a unified framework to better understand how extratropical and tropical forcings jointly modulate basin-scale tropical cyclone (TC) activity. In this study, we examine future changes in TUTTs and their implications for TC activity. Multimodel ensemble-mean projections from 45 Coupled Model Intercomparison Project phase 6 (CMIP6) models suggest a contraction of the Pacific TUTT and an expansion of the Atlantic TUTT as the climate warms. Consistently, future changes in environment-based TC indices indicate that the large-scale conditions will become more favorable for TC genesis and intensification over the central North Pacific but less favorable over the tropical North Atlantic and Gulf of Mexico. Utilizing a TC-permitting large-ensemble dataset [i.e., the Database for Policy Decision-Making for Future Climate Change (d4PDF)] that adequately captures the observed interannual TUTT–TC relationships, we further confirm the impacts of projected TUTT changes on the TC activity in a warmer climate. In contrast, the TUTT–TC relationship is poorly represented in most CMIP6 High-Resolution Model Intercomparison Project (HighResMIP) models; such deficiencies call for caution when assessing future TC risk based on explicitly tracked TCs in these models. Additionally, CMIP6 projections show large intermodel spread in TUTT changes, implying uncertainty in projected TC activity, especially over the central-to-eastern Pacific and the North Atlantic. This intermodel spread is associated with interhemispheric sea surface temperature warming asymmetry, which leads to a meridional shift of the intertropical convergence zone (ITCZ) and the simultaneous weakening or strengthening of TUTTs in the North Pacific and North Atlantic. The potential contributions of anthropogenic aerosol forcing and oceanic circulation to this interhemispheric warming asymmetry are briefly examined.

Climate Change

Qualifying LaBr3:Ce+Sr detector performance for the Mu2e experiment at Fermilab Using the ELBE Accelerator

A LaBr3:Ce+Sr detector will be used to measure the stopped muon captures at the Mu2e experiment at Fermilab. It has been benchmarked in a test beam experiment performed at the ELBE electron accelerator located at the Helmholtz-Zentrum Dresden-Rossendorf, Germany. ELBE’s pulsed bremsstrahlung beam line was set to deliver an average γ -ray energy of between 4–5 MeV. The detector response was mapped to match Mu2e beam conditions, including rates up to 1 Mcps, energy flux, and time structure. A radioactive calibration source was used to mimic the characteristic 1808.7 keV γ -ray, emitted during the atomic muon nuclear capture in the Mu2e aluminum stopping target. The detector energy resolution was measured as a function of the average energy flux: up to 1 TeV/s for 0.34 s, the steady-operation beam-on time and up to 4 TeV/s for 5 ms to get a conservative estimate of the effect of high intensity Mu2e beam fluctuations. The PMT gain variation as a function of the beam spill length and average intensity has been parametrized and corrected for. When a PMT gain correction corresponding to the average beam-spill intensity is applied, the residual effect of beam intensity fluctuations around the average degrades the energy resolution, σ E γ /E γ , at 1808.7 keV from 0.66% to 0.83%.

Huang, Shihua [Purdue U., West Lafayette]

Revealing the Hidden Third Dimension of Point Defects in Two-Dimensional MXenes

Point defects govern many important functional properties of two-dimensional (2D) materials. However, resolving the three-dimensional (3D) arrangement of these defects in multi-layer 2D materials remains a fundamental challenge, hindering rational defect engineering. Here, we overcome this limitation using an artificial intelligence-guided electron microscopy workflow to map the 3D topology and clustering of atomic vacancies in Ti3C2TX MXene. Our approach reconstructs the 3D coordinates of vacancies across hundreds of thousands of lattice sites, generating robust statistical insight into their distribution that can be correlated with specific synthesis pathways. This large-scale data enables us to classify a hierarchy of defect structures-from isolated vacancies to nanopores-revealing their preferred formation and interaction mechanisms, as corroborated by molecular dynamics simulations. This work provides a generalizable framework for understanding and ultimately controlling point defects across large volumes, paving the way for the rational design of defect-engineered functional 2D materials.

2D materials

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path

Stationary Testing of a Pulse Tube Cryocooler Designed to Rotate at 6800 RPM with the Superconducting Rotor of an Electric Machine

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Most of these electric machines employ superconductors on their rotor, thereby requiring cryogenic cooling of the rotating system. An attractive solution is to conductively cool the rotor using a cryocooler that rotates with the rotor. NASA has been developing a 50 K pulse tube cryocooler for its High Efficiency Megawatt Motor (HEMM) that can operate while rotating at up to 6,800 rpm. The first non-rotating tests of the fully assembled HEMM cryocooler have been completed, and initial testing has yielded a no-load temperature of 112 K. A small helium leak into the vacuum system has been identified as the major bottleneck to performance, and based on measurements with the motor turned off there is an observed vacuum-side heat leak of approximately 30 W at 115 K. Immediate next steps for testing are to locate and mitigate the apparent helium leaks into the vacuum chamber. Thus far, compressor performance and the apparent heat lift after correcting for the heat leak are on par with expectations.

Electrified Aircraft Propulsion

Accessing bands with extended quantum metric in kagome Cs 2 Ni 3 S 4 through soft chemical processing

Flat bands that do not merely arise from weak interactions can produce exotic physical properties, such as superconductivity or correlated many-body effects. The quantum metric can differentiate whether flat bands will result in correlated physics or are merely dangling bonds. A potential avenue for achieving correlated flat bands involves leveraging geometrical constraints within specific lattice structures, such as the kagome lattice; however, materials are often more complex. In these cases, quantum geometry becomes a powerful indicator of the nature of bands with small dispersions. We present a simple, soft-chemical processing route to access a flat band with an extended quantum metric below the Fermi level. By oxidizing Ni-kagome material Cs 2 Ni 3 S 4 to CsNi 3 S 4 , we see a two orders of magnitude drop in the room temperature resistance. However, CsNi 3 S 4 is still insulating, with no evidence of a phase transition. Using experimental data, density functional theory calculations, and symmetry analysis, our results suggest the emergence of a correlated insulating state of unknown origin.

Science & Technology - Other Topics

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.