SEARCH · Search NASA

Results for “String Theory”

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.

184 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Spectroscopic Demarcation of Emergent Photons and Spinons in a Dipolar-Octupolar Quantum Spin Liquid

The identification of fractionalized excitations in quantum spin liquids (QSLs) remains a central challenge in condensed matter physics. In dipolar-octupolar (DO) pyrochlores, such as Ce 2⁢ Zr 2 ⁢O 7 , the candidate 𝜋-flux quantum spin ice (QSI) state is predicted to host both gapless emergent photons and a continuum of spinons. However, resolving these modes at zero field is complicated by their spectral overlap and the presence of nonmagnetic scattering near zero energy. Here, in this study, we report neutron scattering experiments on Ce 2 ⁢Zr 2 ⁢O 7 under a magnetic field along the [1, 1, 1] direction. In contrast to previous unpolarized studies at zero field that relied on high-temperature subtraction, we use a same-temperature high-field subtraction protocol to isolate the photon and spinon modes. Leveraging the selective coupling of the magnetic field to the dipolar degrees of freedom, we demonstrate the spectroscopic demarcation of these excitations. We observe that weak fields (≈ 0.15 T) suppress the low-energy photon weight while leaving the high-energy spinon continuum robust, albeit hardened. Our results, supported by gauge mean-field theory and exact diagonalization calculations, provide strong evidence for the 𝜋-flux QSI state and introduce a powerful field-tuning protocol for investigating DO-QSLs.

Gao, Bin [Rice University, Houston, TX (United Sta↗

DOE OSTI2024

Circumventing data imbalance in magnetic ground state data for magnetic moment predictions

Abstract Magnetic materials play a crucial role in the transition to more sustainable forms of energy and electric vehicles. There is an anticipated shortage in magnetic materials in the future, and as a result there is an urgent need to discover and design new magnetic materials. Computational magnetic material design using density functional theory is daunting because of the challenge in identifying magnetic ground states from a combinatorially large set of possibilities. Machine learning offers a path forward by enabling efficient surrogate models that can more readily enumerate these states, but there is a dearth of training data available, and what is available tends to be imbalanced with too much non-magnetic data. In this work we show that the discrete and previously tackled data imbalance that exists at the level of the magnetic ordering leads to an imbalanced continuous distribution with many zeros when the data is unraveled at the atomic magnetic moment level, which subsequently leads to models with low accuracy for magnetic properties. We mitigate this by using a two-part model framework. Our scheme is able to classify atoms into magnetic and non-magnetic with an F1 score and Matthew’s correlation coefficient (MCC) of ~91% and then to provide an implicit embedding representation that maps directly onto the magnitude of the magnetic moment with a mean absolute error of 0.1 μ B . Beyond screening for new magnetic materials, we demonstrate an additional practical use case of our scheme: the provision of good initial guesses for magnetic moments in first-principles electronic relaxations. Such initialization is shown to lead to faster convergence to configurations that lie closer to the ground state.

Computer Science↗

DOE OSTI2022

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO 2 C 2 H 7 + and PO 2 CH 4 + , oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO 2 C 2 H 7 + and PO 2 CH 4 + . The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Effect of the Nature of Both Cation and Anion Substitution on the Structural Symmetry of Li‐Rich 3 d ‐Metal Chalcogenide Electrodes

Abstract Li‐rich layered chalcogenides have recently led to better understanding of the anionic redox process and its associated high capacity while providing ways to overcome its practical limitations of voltage fade and irreversibility. This study reports on the feasibility of triggering anionic activity in Li 2 TiS 3 , through anionic substitution (Se for S) or cationic substitution (Fe for Ti). Herein, the chalcogenide chemical space is further explored to prepare mono‐substituted Li 1.7 Ti 0.85 Mn 0.45 Ch 3 (Ch = S/Se) and doubly substituted cationic and anionic phases (Li 1.7 Ti 0.85 Fe 0.45 S 3‐z Se z ) which crystallize either in the O3‐ or O1‐type structures depending upon substituents. All series show a bell‐shape capacity variation as function of the transition metal (TM) substitution degree with values up to 240 mAh g −1 . For specific compositions, a structural O3 to O1 phase transition is observed upon Li removal, which is not reversible upon Li re‐insertion due to kinetic limitations and negatively affects long‐term cycling performance. Density functional theory (DFT) calculations confirm the O3/O1 relative stability along the different series and point subtle electronic differences in the TM‐doping, rationalizing the structural and electrochemical behaviors of these phases upon cycling. These findings provide further insights into the link between structural and electronic stability, which is of key importance for designing chalcogenide‐based anionic redox compounds.

Chemistry↗

DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗

DOE OSTI2025

Field-tunable BKT and quantum phase transitions in spin-$\frac{1}{2}$ triangular lattice antiferromagnet

Quantum magnetism is one of the most active fields for exploring exotic phases and phase transitions. The recently synthesized Na 2 BaCo(PO 4 ) 2 (NBCP) is an ideal material incarnation of the spin-$\frac{1}{2}$ easy-axis triangular lattice antiferromagnet (TLAF). Experimental evidence shows that NBCP hosts the spin supersolid state with a giant magnetocaloric effect. Theory further predicts that magnetic fields can drive NBCP through Berezinskii-Kosterlitz-Thouless (BKT) and other richer quantum phase transitions. However, detecting these transitions is challenging, as they onset at ultralow temperatures near 60 mK and require high magnetization sensitivity. Using a newly developed gradient force magnetometer in a dilution refrigerator, we mapped the magnetic susceptibility phase diagram down to 30 mK. Our results provide a more comprehensive and accurate understanding of BKT melting of spin supersolidity and several field-tunable quantum phase transitions, which establish NBCP as a model platform for frustrated magnetism and highlight potential applications of its giant magnetocaloric effects.

BKT transition↗

DOE OSTI2022

Efficient electrocatalytic valorization of chlorinated organic water pollutant to ethylene

Electrochemistry can provide an efficient and sustainable way to treat environmental waters polluted by chlorinated organic compounds. However, the electrochemical valorization of 1,2-dichloroethane (DCA) is currently challenged by the lack of a catalyst that can selectively convert DCA in aqueous solutions into ethylene. Here, in this work, we report a catalyst comprising cobalt phthalocyanine molecules assembled on multiwalled carbon nanotubes that can electrochemically decompose aqueous DCA with high current and energy efficiencies. Ethylene is produced at high rates with unprecedented ~100% Faradaic efficiency across wide electrode potential and reactant concentration ranges. Kinetic studies and density functional theory calculations reveal that the rate-determining step is the first C–Cl bond breaking, which does not involve protons—a key mechanistic feature that enables cobalt phthalocyanine/carbon nanotube to efficiently catalyse DCA dechlorination and suppress the hydrogen evolution reaction. The nanotubular structure of the catalyst enables us to shape it into a flow-through electrified membrane, which we have used to demonstrate >95% DCA removal from simulated water samples with environmentally relevant DCA and electrolyte concentrations.

Choi, Chungseok [Yale Univ., New Haven, CT (United↗

DOE OSTI2026

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗

DOE OSTI2026

Pinning ångström-size solid ionic channels for rare-earth element separation

High-purity rare-earth elements are essential for modern technologies, yet current solvent extraction processes are energy-intensive and environmentally harmful because of inadequate selectivity and ligand toxicity. Although combining size exclusion and binding affinity can improve lanthanide separation, the role of long-range confinement remains underexplored. Here we report lanthanide separation in aqueous systems using extremely confined manganese oxide solid ionic channels with optimized layer spacing. Different lanthanides induce distinct solid-state phase transformations in manganese oxide, creating a strong driving force for separation. Two lanthanide groups, differing by ~1.4 Å in spacing, were identified and confirmed to be stable by density functional theory. The narrower confinement of heavier Group II lanthanides improves cross-group separation by increasing the dehydration barrier for lighter Group I lanthanides without inducing strong binding. Here, we further developed a strategy to pin the confinement dimensions and enhance same-group separation, increasing enrichment factors for La–Nd and La–Pr pairs from 1.6 ± 0.1 and 1.5 ± 0.1 to 5.4 ± 0.1 and 4.2 ± 0.1, respectively.

Chemical engineering↗

DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗

DOE OSTI2024

Why Perovskite Thermal Stress is Unaffected by Thin Contact Layers

Metal halide perovskite photovoltaics have emerged as a high efficiency, low-cost alternative that can potentially rival or enhance conventional silicon technology. Despite exceptional initial power conversion efficiencies, achieving compliance with international standards and widespread adoption requires further enhancements to their operational stability. Notably, addressing mechanical strain and stress in brittle perovskites has emerged as a pivotal approach to mitigate chemical degradation and improve reliability during thermal cycling. Here, in this study, a popularized strain engineering strategy is investigated in which a high coefficient of thermal expansion (CTE) hole transport layer (i.e., PDCBT) is cast onto inorganic perovskite (CsPbI 2 Br) at 100 °C. Contrary to previously published results, the X-ray diffraction (XRD):Sin 2 ψ and substrate curvature measurement techniques show that the hole transport layer has no discernible impact on perovskite strain. The accuracy of the XRD:Sin 2 ψ method for measuring strain is highlighted in contrast to an analysis based on shifts of single XRD peaks which can be influenced by multiple artifacts. The findings in this study are in accordance with mechanics theory: thin layers are unable to induce significant strain changes in perovskite thin films as the force they apply is negligible compared to that applied by a thick and stiff substrate.

14 SOLAR ENERGY↗

NASA NTRSDate not supplied

Fluid-Thermal-Structural Interactions Induced by an Asymmetric Shock-Wave/Boundary-Layer Interaction in a Mach-6 Compression Corner

An experimental study is conducted of the fluid-thermal-structural interaction of a clamped compliant panel exposed to a three dimensional shock-wave/boundary-layer interaction (SWBLI) induced by a Mach-6 compression ramp with a spanwise nonuniform incoming boundary layer. The nonuniform boundary layer was produced by placing trips on one side of the upstream flat plate, resulting in largely turbulent flow on the tripped side and transitional flow on the untripped side. Measurements of the flowfield confirmed that the tripped boundary layer contained elevated levels of unsteadiness, and the SWBLI was observed to vary from attached to fully separated as the ramp angle was increased from 10◦ to 38◦; the separation region on the tripped side of the panel was noticeably smaller, showing the elevated turbulence levels of the tripped-side flow to remain relatively localized rather than diffusing across the whole model. Full-field, time-resolved panel deformations were measured using high-speed photogrammetry and the vibrational response at each compression angle was characterized. Although the measured modes conformed largely to those from classical clamped-plate theory, some skewing of the mode shapes was observed. IR thermography highlighted regions of the compliant region where elevated temperatures were likely to promote thermal softening effects to the transient panel response. The quasi-static deformation and stress field was used to characterize the internal stress factor of each mode and showed a meaningful relationship between transient panel response and stress contained within each mode: modes with antinodes lying in high-stress areas of the plate tended to exhibit increases in vibrational frequency and decreases in vibrational power, whereas the opposite was true for modes with antinodes in low-stress areas.

Spectral Proper Orthogonal Decomposition↗

DOE OSTI2025

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗

DOE OSTI2026

Physical Modeling and Design of a Nonvolatile Optically Gated High‐Power Diamond Transistor

In this work, we present the theory and modeling framework of a diamond optically gated junction field‐effect transistor (DOGFET). The device utilizes nitrogen substitutional centers in type‐1b diamond to optically modulate a p‐ boron doped diamond channel. Using sub‐gap lasers with intensities as low as 100 W/cm 2 , electrons are optically excited from substitutional nitrogen sites to the conduction band of the diamond substrate, thus enabling the optical gate to exercise control on modulating the space‐charge region at the junction and therefore the channel conductivity. We show that the device can deliver a current of 7 μA/μm, or equivalently 1750 A/cm 2 , while switching at a frequency greater than 100 kHz, in a form factor of 5 μm 2 . The breakdown voltage is found to be greater than 1850 V, with a breakdown field strength of ~13 MV/cm. Moreover, the device supports nonvolatile operation with a “memory effect” enabling single transistor state retention. The presented simulation framework provides a physically grounded insight into the limits and opportunities of optoelectronic diamond systems.

Engineering - Electronic and electrical engineerin↗

NASA NTRS1992

On the isotopic signature of recent solar-wind nitrogen

One of the most intriguing discoveries yielded by the Apollo samples was evidence pointing towards a significant long-term change in the composition of the sun. Such a change, of the size inferred from the lunar sample data, is inconsistent with present theories of solar evolution. Consequently, there is much interest in exploring this phenomenon as closely as possible, to determine exactly what compositional changes have taken place and whether those changes really did take place in the sun, or whether the cause lies elsewhere. The reason why we can use the moon to analyze the elements in the sun is that the sun emits a stream of ions, known as the solar wind, whose composition, on average, is believed to be the same as that in the surface regions of the sun. When the solar-wind ions hit the surface of the moon, many of them penetrate a short distance into the dust grains lying on the lunar surface. Thus, after a grain has sat on the lunar surface for a while, it has a rim of material that is partly lunar and partly solar in composition. For most chemical elements, the difference between lunar and solar composition is so sufficiently small that the solar elements cannot be detected, but for a handful of elements that are missing from the moon, their solar 'signature' can be observed in samples of lunar soil brought back by the astronauts. Among those elements is nitrogen, the most common element in the air we breathe, but very rare indeed on the moon. Our analytical techniques are not sophisticated enough yet to enable us to analyze individual lunar soil grains for nitrogen, much less to zero in on just the nitrogen in the surface of such a grain. Consequently we are forced to analyze samples consisting of many different grains, each of which could have experienced its own individual history. This makes it difficult to identify the nitrogen implanted in grain surfaces, and also to define the age of a sample.

Kim, Y.↗

DOE OSTI2026

Studies of laser stimulated photodetachment from nanoparticles for particle charge measurements

Determining nanoparticle charge is more challenging than that for microparticles due to change in the particle size during the synthesis substantial plasma property variations, and difficulties in visualizing individual particles, rendering conventional microparticle charge diagnostics ineffective in dusty plasma. In this work, we utilized laser-stimulated photodetachment (LSPD) to deduce the mean charge of nanoparticles. Nanoparticles were grown in an Ar/C 2 H 2 mixture using a capacitively coupled RF discharge and the LSPD induced changes in the electron current monitored by a cylindrical Langmuir probe. LSPD signals were obtained and analyzed across different dust growth phases. The prolonged decay of electron current pulses was attributed to the presence of residual negative ions, caused by the effective electrostatic trapping of these ions and the potential post—LSPD re-formation of new ones. The particle charge was estimated by combining the laser-stimulated photodetachment signal from the probe with the dust density obtained from laser-light extinction using the measured nanoparticle size distribution. For a nanoparticles size range of approximately 100–250 nm and mean diameter of $d$ p ∼154.37 nm, the effective mean charge was estimated to be $\langle$$Q$$\rangle$ d $≈$ 37 elementary charge units. The measured charge values are lower than those predicted by orbital motion limited theory, which may be attributed to significant electron depletion in the nanodusty plasma. LSPD results in Ar/C 2 H 2 nano-dusty plasma confirm the applicability of this method for estimating individual nanoparticle charges. However, it has also been demonstrated that electron detachment from residual background negative ions can influence the detachment current decay and must be carefully considered.

dust particle charge↗

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