Search NASASearch

SEARCH · Search NASA

Results for “high pressure and high temperature”

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,689 records · Page 42

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID

Dirac Fermions and Flat Bands in Phosphorus Carbide Nanotubes: Structural and Quantum Phase Transitions in a Quasi-One-Dimensional Material

Chemically realistic quasi-one-dimensional (1D) materials in which Dirac Fermions and highly degenerate flat bands coexist intrinsically at the Fermi level are exceedingly rare, while representing a highly desirable platform for correlated and topological quantum phenomena. Here, in this work, using specialized symmetry-adapted first-principles calculations we predict a new class of nanomaterials─phosphorus carbide nanotubes (P 2 C 3 NTs)─obtained by rolling monolayer P 2 C 3 , a two-dimensional material shown in a previous letter to host “double Kagome bands”. Both armchair and zigzag P 2 C 3 NTs are stable at room temperature and feature the rare coexistence of Dirac crossings and multiple flat bands at the Fermi level inherited from the underlying honeycomb–Kagome lattice, with the flat bands resilient to elastic deformations. Under large strain, the structure transforms from honeycomb–Kagome to “brick-wall”, accompanied by multiple coupled structural and quantum phase transitions. We also uncover localized edge states, spin splitting from vacancies and dopants, and strain-tunable magnetism. Together, these results establish P 2 C 3 NTs as a chemically specific and mechanically tunable 1D material platform with potential applications in quantum hardware and spintronics.

carbon nanotubes

Acquisitive root exploration strategies help maintain higher peak sap flux rates during summer drought, but more root biomass does not

Roots are responsible for soil water uptake, yet little is known about how variation in fine-root traits relates to whole-tree water movement, particularly during periods of drought. By combining a 3-yr dataset monitoring sap flow rates with measures of fine-root biomass, length, and morphology across 10 tree species, we addressed hypotheses concerning the role of fine-root system size and morphology in determining tree responses to summer drought as well as potential changes in their relative importance under more severe droughts. Greater fine-root biomass and length did not enable trees to maintain high rates of sap flow during drought, whereas the morphological traits, specific root length and specific root area, were linked to sustained sap flow rates during drought. We found that all species, irrespective of root traits, progressively reduced their sap flow under more extreme drought conditions, although more acquisitive root morphology was still associated with smaller reductions. Our results run counter to long-standing assumptions that larger root systems are better able to access soil water and maintain photosynthetic activity during drought. Instead, we find evidence that root morphology at least partially determines the capacity for water uptake and movement as soil moisture declines.

drought

Considerations for Maritime Nuclear Technologies, Economic Viability and Public Acceptance

The Maritime Nuclear Application Group (MNAG) is a working group convened by the National Reactor Innovation Center at Idaho National Laboratory (INL), the American Bureau of Shipping, and Morgan, Lewis, and Bockius LLP. This report documents an MNAG examination of considerations relevant to implementing nuclear technology in commercial maritime applications. In general, two types of use case are examined: maritime nuclear power plants and nuclear reactors used on board shipping vessels for propulsion and other ship needs. The report finds that there may be economic benefits related to maritime nuclear technologies, including the flexible deployment of maritime nuclear reactors, which would allow them to complement land-based nuclear projects, and operational differences for nuclear cargo ships that may lead to an overall increase in revenue. High-level analyses in this report show that, based on general small modular reactor and microreactor cost estimates developed by INL, maritime nuclear reactors may be economically competitive for electricity production in remote regions and for use in the propulsion of large cargo ships. Besides economic viability, public acceptance will be key to implementing maritime nuclear technologies. The report discusses the public’s current perception of nuclear technologies. Engaging with the public will be important to improving this perception. The report discusses some key benefits and risks associated with maritime nuclear technologies. Benefits include the creation of jobs, the production of reliable energy, and the potential to improve air quality. Risks that concern the public are the potential for radioactive releases during operation and decommissioning, as well as those related to waste management. Communicating the benefits and the risks of maritime nuclear technologies will be essential to improving public perception.

22 - GENERAL STUDIES OF NUCLEAR REACTORS

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

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

Real-space visualization of a defect-mediated charge density wave transition

Here, we study the coupled charge density wave (CDW) and insulator-to-metal transitions in the 2D quantum material 1T-TaS 2 . By applying in situ cryogenic 4D scanning transmission electron microscopy with in situ electrical resistance measurements, we directly visualize the CDW transition and establish that the transition is mediated by basal dislocations (stacking solitons). We find that dislocations can both nucleate and pin the transition and locally alter the transition temperature T c by nearly ~75 K. This finding was enabled by the application of unsupervised machine learning to cluster five-dimensional, terabyte scale datasets, which demonstrate a one-to-one correlation between resistance—a global property—and local CDW domain-dislocation dynamics, thereby linking the material microstructure to device properties. This work represents a major step toward defect-engineering of quantum materials, which will become increasingly important as we aim to utilize such materials in real devices.

4D-STEM

Machine learning pipeline for denoising low signal-to-noise ratio and out-of-distribution transmission electron microscopy datasets

High-resolution transmission electron microscopy (HRTEM) is crucial for observing material’s structural and morphological evolution at Angstrom scales, but the electron beam can alter these processes. Devices such as CMOS-based direct-electron detectors operating in electron-counting mode can be utilized to substantially reduce the electron dosage. However, the resulting images often lead to a low signal-to-noise ratio, which requires frame integration that sacrifices temporal resolution. Several machine learning (ML) models have been recently developed to successfully denoise HRTEM images. Yet, these models are often computationally expensive, and their inference speeds on GPUs are outpaced by the imaging speed of advanced detectors, precluding in situ analysis. Furthermore, the performance of these denoising models on datasets with imaging conditions that deviate from the training datasets has not been evaluated. To mitigate these gaps, we propose a new self-supervised ML denoising pipeline specifically designed for time-series HRTEM images. This pipeline integrates a blind-spot convolution neural network with pre-processing and post-processing steps, including drift correction and low-pass filtering. Results demonstrate that our model outperforms various other ML and non-ML denoising methods in noise reduction and contrast enhancement, leading to improved visual clarity of atomic features. Additionally, the model is drastically faster than U-Net-based ML models and demonstrates excellent out-of-distribution generalization. The model’s computational inference speed is in the order of milliseconds per image, rendering it suitable for application in in-situ HRTEM experiments.

36 MATERIALS SCIENCE

Chronology of ungrouped achondrites Northwest Africa 11119 and 7325/8486: implications for early evolution of a heterogeneous crust on a differentiated planetesimal

The chronology of ungrouped achondrites provides key insights into the timeline of igneous activity in the early Solar System. Two ungrouped achondrites – Northwest Africa (NWA) 11119 (andesite-dacite) and NWA 8486 and its pair NWA 7325 (olivine gabbro) – are the focus of this study, where their chronologies are investigated using multiple high-resolution techniques. Here we report the lead-lead (207Pb-206Pb) and manganese-chromium (53Mn-53Cr) systematics of NWA 11119, as well as the 207Pb-206Pb systematics of NWA 8486 alongside aluminum-magnesium (26Al-26Mg) systematics for NWA 7325. The U-corrected 207Pb-206Pb ages of NWA 11119 and the combined ages of NWA 7325/8486 are 4566.4 ± 0.8 Ma and 4563.8 ± 1.9 Ma, respectively. Additionally, we report the 53Mn-53Cr age of NWA 11119 to be 4564.4 ± 2.5 Ma and the 26Al-26Mg age of NWA 7325/8486 to be 4563.1 ± 0.3 Ma. The formation of NWA 11119 requires partial melting of a (likely chondritic) source reservoir leading to eruption of Si rich, alkali depleted magmas, while NWA 7325/8486 likely formed from a chemically fractionated reservoir with superchondritic Al/Mg. The clear geochemical and isotopic similarities of these achondrites, combined with the chronology reported here, is suggestive of formation of these two ungrouped achondrites on a common parent body which likely formed in the inner Solar System and experienced early differentiation under reducing conditions. If these achondrites did share a parent body, it would suggest that primary asteroids in the early Solar System commonly produced mineralogically and geochemically heterogeneous crusts. While mineralogical and geochemical heterogeneity is known to exist in the (mostly mafic) crust of asteroid Vesta, our findings show that even more significant crustal heterogeneity (representing felsic and mafic compositions) may exist on other asteroids.

Cartwright, Julia [University of Leicester, UK]

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]

Photophysical and Time‐resolved Infrared Properties of Long‐Lived Rhenium(I) 4,5‐Diazafluorene Tricarbonyl Chromophores

This report investigates the synthesis, structural characterization, fundamental molecular photophysics, electrochemistry, UV-Vis spectroelectrochemistry, and time-resolved infrared spectroscopic properties of eight [fac-Re(dafR)(CO) 3 L] 0/+ complexes, where R=ethyl [(dedaf); 1, 3, 5, 7] or H [(dafH); 2, 4, 6, 8] and L=Cl− (1, 2), imidazole [(Im); 3, 4], 4-ethylpyridine [(4-Etpy); 5, 6], or pyridine [(py); 7, 8]. Universally, 1–8 yield higher energy photoluminescence (PL) emission bands and higher PL quantum yields (up to 53 %) than the classic 2,2’-bipyridine (bpy) and 1,10-phenanthroline (phen) ligated Re(I) tricarbonyl complexes. The excited state lifetimes of 1–8 lie between those corresponding to the bpy and phen derivatives, ranging from 120 and 1300 ns at room temperature. Combinations of reductive UV-Vis spectroelectrochemistry, transient absorption spectroscopy, and time-resolved infrared spectroscopy consistently assigned the lowest excited states in 1–8 being of metal-to-ligand charge transfer (MLCT) character. These new ReI MLCT chromophores follow classic energy gap law behavior and possess the characteristics necessary for serving as valuable photosensitizers suitable to energize excited state electron and energy transfer photochemistry.

Charge transfer

Cross-scale interaction between microturbulence and meso-scale reversed shear Alfvén eigenmodes in DIII-D plasmas

Abstract This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.

Physics

On-Chip Tuning of Superconductivity in Fullerides via Current-Driven Rb + Intercalation

An in-operando electro-intercalation method for the on-chip synthesis of alkali-metal-intercalated materials and their Raman spectroscopic and transport characterization in ultrahigh vacuum (UHV) is developed. We apply this method to synthesize fulleride superconductors via Rb + intercalation into a C 60 film. During the intercalation, we monitor the stoichiometry via UHV-Raman spectroscopy and probe superconductivity via transport measurements. An increase of the superconducting transition temperature from 7.0 K to 14.5 K is observed when the stoichiometry is tuned from Rb 2.7 C 60 to Rb 3 C 60 . In our experiment, an ionic Rb+ flux into the host material is induced by an applied electronic current via a Butler–Volmer-type mechanism. Electro-intercalation captivates through improved stoichiometric precision, the ability to smoothly vary stoichiometry via duration of current application, and the absence of a lower limit of the volume of the host material. It represents a powerful concept for the on-chip synthesis of intercalated materials, battery research, and beyond.

Raman

Four Years of Type Ia Supernovae Observed by TESS: Early-time Light-curve Shapes and Constraints on Companion Interaction Models

Abstract We present 307 type Ia supernova (SN) light curves from the first 4 yr of the Transiting Exoplanet Survey Satellite mission. We use this sample to characterize the shapes of the early-time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise <10% of the peak flux to avoid biases in the early-time light-curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power-law index t β 1 of the early-time light curves isβ 1 = 1.93 ± 0.57, and the mean rise time to peak is 15.7 ± 3.5 days. The underlying population distribution forβ 1 may instead consist of a Gaussian component with mean 2.29, width 0.34, and a long tail extending to values less than 1.0. We find that the data can rarely distinguish between models with and without companion interaction models. Nevertheless, we find three high-quality light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence for the companion interactions is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power-law model. Taking the 74 SNe together, we calculate 3σupper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii >31R ⊙ (separations >5.7 × 10 12 cm, 99.9% confidence level) and disfavor companions having Roche lobe radii >10R ⊙ (separations >1.9 × 10 12 cm, 95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Astronomy & Astrophysics

Design and Realization of Ohmic and Schottky Interfaces for Oxide Electronics

Understanding band alignment and charge transfer at complex oxide interfaces is critical to tailoring and utilizing their diverse functionality. Toward this goal, both Ohmic- and Schottky-like charge transfers at oxide/oxide semiconductor/metal interfaces are designed and experimentally validated. A method for predicting band alignment and charge transfer in ABO 3 perovskites is utilized, where previously established rules for simple semiconductors fail. The prototypical systems chosen are the rare class of oxide metals, SrBO 3 with B = V–Ta, when interfaced with the multifaceted semiconducting oxide, SrTiO 3 . For B = Nb and Ta, it is confirmed that a large accumulation of charge occurs in SrTiO 3 due to the higher energy Nb and Ta states relative to Ti. Furthermore, this gives rise to a high mobility metallic interface, which is an ideal epitaxial oxide/oxide Ohmic contact. On the contrary, for B = V, there is no charge transfer into the SrTiO 3 interface, which serves as a highly conductive epitaxial gate metal. Going beyond these specific cases, this work opens the door to integrating the vast phenomena of ABO 3 perovskites into a wide range of practical devices.

36 MATERIALS SCIENCE

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

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

Manard, Benjamin [ORNL] (ORCID:0000000207400627)

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.