Search NASA⌕ Search

SEARCH · Search NASA

Results for “98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL PROTECTION”

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

Modern insights into the mechanisms of neptunium oxalate decomposition

Neptunium oxalate (Np(C 2 O 4 ) 2 ·6H 2 O) is a historically relevant solid phase used in nuclear processing as a precursor for neptunium dioxide (NpO 2 ). Although Np oxalate has been synthesized and used for NpO 2 production for decades, the thermal decomposition mechanism of this phase remains poorly understood and has not been evaluated in over 30 years. Conflicting reports in historical literature suggest either a direct conversion from anhydrous oxalate to NpO 2 or a decomposition that includes the formation of Np carbonate or oxidized Np intermediate phases. In this work, we reexamine the decomposition pathway of Np(C 2 O 4 ) 2 ·6H 2 O using thermal analysis coupled with evolved gas analysis and temperature-dependent Raman spectroscopy to elucidate decomposition mechanisms and intermediate phases using modern analytical techniques. Thermal analysis revealed a three-stage decomposition process, including dehydration below 200 °C, oxalate breakdown between 170 and 370 °C, and NpO 2 formation by 500 °C. However, an unidentified plateau in the thermal data was observed during measurements. Raman spectroscopy confirmed the stages of decomposition, and in the analysis of potential intermediate phases, no carbonate phases or Np 2 O 5 were identified. Raman data suggest that residual oxalate or nonstoichiometric oxide are present during decomposition before pure NpO 2 is formed. These findings clarify aspects of the Np oxalate decomposition mechanism and address longstanding discrepancies in the literature, with a specific focus on Np-specific materials chemistry.

Lawson, Kathryn M. [Oak Ridge National Laboratory ↗
DOE OSTI2025

Review on Preprocessing Strategies, Deactivation, Thermal Safety, and Future Perspectives in Lithium-Ion Battery Recycling

The rapid growth in the use of lithium-ion batteries (LIBs) in electric vehicles, consumer electronics, and renewable energy storage has made effective end-of-life management essential. Recycling LIBs is critical not only for resource recovery and environmental protection but also for ensuring safety and economic viability. This review focuses on the preprocessing technologies that precede typical recovery processes, including disassembly, sorting, discharging, electrolyte removal, dismantling, thermal treatment, separation, and flotation. These steps play a foundational role in determining the efficiency, safety, and environmental impact of LIB recycling. LIBs pose substantial fire and explosion risks due to residual charge, flammable electrolytes, and reactive materials. The conditions and successive progression of the exothermic reactions which lead to thermal runaway has been discussed. It also explores secure deactivation techniques such as external circuit discharge, saline immersion, and thermomechanical methods, alongside fire prevention strategies including the use of flame retardants, elimination of oxidants, and reduction of heat generation and accumulation. Challenges and future directions are outlined, highlighting the need for standardized designs, automation, and safer, more sustainable recycling infrastructure. Furthermore, this review is distinguished by its focused analysis of preprocessing and deactivation steps, with particular attention to the thermal safety engineering aspects of LIB recycling.

Battery deactivation↗
DOE OSTI2026

Dark Energy Survey year 6 results: Magnification modeling and its impact on galaxy clustering and galaxy-galaxy lensing cosmology

Gravitational lensing magnification alters the observed spatial distribution of galaxies and must be accounted for to prevent biases in cosmological probes of the large-scale structure. We investigate its effects on the Dark Energy Survey Year 6 galaxy clustering and galaxy-galaxy lensing analyses using the fiducial lens (position tracer) sample M ag L im++. Magnification bias is parameterized by a coefficient that describes the response of the number of selected objects per unlensed area element to a change in the lensing convergence. We quantify this coefficient using the BALROG synthetic source injection catalog to account for the complexity of the selection function, and compare these results with simplified estimates. The resulting values of the magnification coefficients for each redshift bin are [3.16 ± 0.08, 2.76 ± 0.21, 4.09 ± 0.15, 4.42 ± 0.16, 4.90 ± 0.29, 4.83 ± 0.25]. Relative to Year 3, this analysis provides more precise and accurate magnification bias estimates through a larger BALROG area and reweighting to better match the data properties. Here, the cosmological results are robust when tested against various magnification parameter prior choices and also when adding cross-clustering between lens redshift bins. Neglecting magnification, however, introduces significant systematic shifts: relative to the fiducial analysis with Gaussian priors centered on the BALROG -derived estimates, we observe shifts of 1.37σ in S 8 and -0.84σ in Ω m (with cosmic shear included: -0.61σ in S 8 and -0.71σ in Ω m ), in agreement with findings from simulated data, demonstrating that magnification must be modeled to avoid biases. Freeing the magnification bias in lens bin 2 leads to unphysical negative values, further justifying its exclusion from the fiducial Year 6 analysis.

Cosmological parameters↗
DOE OSTI2026

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗
DOE OSTI2025

Preliminary Evaluation of Joint Electricity-Hydrogen Concept of Operations

An initial thermal power dispatch (TPD) concept of operations was evaluated that couples a nuclear power plant to a nearby hydrogen production plant. GSE Systems’ generic pressurized water reactor full-scope simulator was modified with a TPD model comprised of a thermal power extraction and delivery system. A prototype human-system interface (HSI) was developed to interact with the TPD model and allow participants to execute the basic operating scenarios for normal operations. Four retired operators performed the evaluation, and due to COVID-19 travel restrictions, the original in-person experimental design was restructured to support a remote participator evaluation using a web meeting platform. Data from operator feedback, observations from the research team, and quantitative survey responses revealed that the initial TPD concept of operations is feasible. The operators were comfortable with the engineered system and HSI and could manage it without adverse impacts to reactor power, plant safety, or equipment. Findings are discussed in terms of both the TPD system design and HSI performance.

human factors↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
DOE OSTI2026

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

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

Dark Energy Survey Year 6 results: Redshift calibration of the MagLim++ lens sample

In this work, we derive and calibrate the redshift distribution of the MagLim++ lens galaxy sample used in the Dark Energy Survey Year 6 (DES Y6) 3 x 2pt cosmology analysis. The 3 x 2pt analysis combines galaxy clustering from the lens galaxy sample and weak gravitational lensing. The redshift distributions are inferred using the SOMPZ method - a Self-Organizing Map framework that combines deep-field multi-band photometry, wide-field data, and a synthetic source injection ( B alrog) catalog. Key improvements over the DES Year 3 (Y3) calibration include a noise-weighted SOM metric, an expanded Balrog catalogue, and an improved scheme for propagating systematic uncertainties, which allows us to generate O(10 8 ) redshift realizations that collectively span the dominant sources of uncertainty. These realizations are then combined with independent clustering-redshift measurements via importance sampling. The resulting calibration achieves typical uncertainties on the mean redshift of 1-2%, corresponding to a 20-30% average reduction relative to DES Y3. We compress the n(z) uncertainties into a small number of orthogonal modes for use in cosmological inference. Marginalizing over these modes leads to only a minor degradation in cosmological constraints. Here, this analysis establishes the MagLim++ sample as a robust lens sample for precision cosmology with DES Y6 and provides a scalable framework for future surveys.

dark energy↗
DOE OSTI2026

Advancing Protein Display on Bacterial Spores through an Extensive Survey of Coat Components

The profound stability of bacterial spores makes them a promising platform for biotechnological applications like biocatalysis, bioremediation, drug delivery, etc. However, though the Bacillus subtilis spore is composed of >40 types of proteins, only ∼12 have been explored as fusion carriers for protein display. Here, we assessed the suitability of 33 spore proteins (SPs) as enzyme display carriers by direct allele tagging at native genomic loci. Of the 33 SPs investigated, 26 formed functional fusions with β-glucuronidase (GUS)─a ∼272 kDa homotetramer. This almost triples the number of SPs assessed for enzyme display and doubles the number of functional fusions documented in the literature. We quantitatively assessed 1) SP promoter activation dynamics, 2) GUS activity on spores, 3) surface availability, and 4) protection from thermal and proteolytic degradation. Multicopy expression and pairwise coexpression of the most promising SP-GUS fusions highlighted the complexity of spore structure/assembly and the difficulty in predicting compatibility between different SP fusions. We also assessed the suitability of engineered spores to degrade PET (polyethylene terephthalate) films and found that surface-exposed SPs were most effective. Beyond the broad survey, a key outcome of our work was the identification of SscA (small spore coat assembly protein A) as an effective spore display carrier. SscA supported enzyme activity at least 4-fold higher than any other SP, including the well-established anchor, CotY. We attribute this to its promoter, which demonstrated early and sustained activation relative to other SPs and its small size (∼3 kDa), which likely minimally interferes with enzyme folding, oligomerization, and activity. Labeling and genetic studies, its hydrophobic nature, and low surface availability suggest that SscA assembles within the inner spore coat, which makes it stabilizing and suitable for many biocatalytic applications. Overall, this work serves as a knowledge base to advance the biotechnological utility of B. subtilis spores.

Bacillus subtilis↗
DOE OSTI2024

A novel digital lifecycle for Material‐Process‐Microstructure‐Performance relationships of thermoplastic olefins foams manufactured via supercritical fluid assisted foam injection molding

Abstract This research significantly enhances the applicability of thermoplastic olefins (TPOs) in the automotive industry using supercritical N 2 as a physical foaming agent, effectively addressing the limitations of traditional chemical agents. It merges experimental results with simulations to establish detailed material‐process‐microstructure‐performance (MP2) relationships, targeting 5–20% weight reductions. This innovative approach labeled digital lifecycle (DLC) helps accurately predict tensile, flexural, and impact properties based on the foam microstructure, along with experimentally demonstrating improved paintability. The study combines process simulations with finite element models to develop a comprehensive digital model for accurately predicting mechanical properties. Our findings demonstrate a strong correlation between simulated and experimental data, with about a 5% error across various weight reduction targets, marking significant improvements over existing analytical models. This research highlights the efficacy of physical foaming agents in TPO enhancement and emphasizes the importance of integrating experimental and simulation methods to capture the underlying foaming mechanism to establish material‐process‐microstructure‐performance (MP2) relationships. Highlights Establishes a material‐process‐microstructure‐performance (MP2) for TPO foams Sustainably produces TPO foams using supercritical (ScF) N 2 with 20% lightweighting Shows enhanced paintability for TPO foam improved surface aesthetics Digital lifecycle (DLC) that predicts both foam microstructure and properties DLC maps process effects & microstructure onto FEA mesh for precise prediction

Engineering↗
DOE OSTI2024

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↗
DOE OSTI2026

Ionic Precursors Transformed Into Vinyl Acetate Synthesis Catalyst via Reaction-Driven Restructuring

Conventional preparation of supported bimetallic catalysts relies on solution-mediated metal salt immobilization and pre-formation of alloy nanoparticles before reaction. Here, we report a fundamentally different synthesis strategy of using a physical mixture of salt precursors to generate an active catalyst during reaction. The catalytic structure is generated in situ from Pd3(OAc)6, Au(OH)3, and KOAc through H2 treatment and reaction-driven restructuring under vinyl acetate monomer (VAM) synthesis conditions. Ascertained from in situ X-ray diffraction and operando infrared spectroscopy analyses, reduction treatment produces segregated Pd and Au domains, and subsequent exposure to a VAM reaction mixture triggers dynamic extraction of Pd from the metal surface. This latter process, mediated by acetate-assisted redox cycles, facilitates Pd migration toward Au domains to form a near-surface localized Pd50Au50 alloy phase. Monometallic Pd domains serve as a reservoir of Pd to the alloy phase, leading to and sustaining a more Pd-enriched active surface and a higher population of accessible Pd sites, compared to a conventionally prepared K-PdAu/SiO2 catalyst. Consequently, this leads to a twofold increase in the VAM formation rate, demonstrating highly active bimetallic catalysts can be generated through the gas-phase treatment of physically mixed ionic precursors.

Cha, Byeong Jun [Rice University]↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
DOE OSTI2026

Radiative divertor detachment and impurity transport with nitrogen and neon seeding in KSTAR H-mode plasmas

Achieving core-edge compatible divertor detachment is a critical requirement for stable operation in future fusion devices. This study compares nitrogen and neon seeding in KSTAR H-mode plasmas with carbon walls, combining experiments and SOLPS-ITER modelling to evaluate their radiative dissipation and core-edge compatibility. Experimentally, N seeding achieved stronger divertor detachment, with larger reductions in target particle and heat fluxes, a higher divertor radiation fraction, and a lower core radiation fraction than Ne. In contrast, Ne seeding triggered a significant rise in core radiation followed by H–L back transitions, limiting the maximum total radiated power fraction to roughly half that of N. SOLPS-ITER simulations reproduced the experimental trends and revealed that the better core-edge compatibility of N arises from its higher divertor retention in addition to its higher cooling factor. The relative positions of the stagnation points of impurity poloidal velocity and ionization sources did not explain the different divertor compression. Instead, in the present modelling, the higher impurity parallel particle flux, resulting from the higher parallel impurity velocity, explains the stronger nitrogen impurity compression in the divertor region. The parallel temperature distribution with N was more favourable for achieving higher impurity parallel velocity than with Ne, because the impurity velocity is governed by modifications of the main ion flow due to friction and thermal forces, both of which strongly depend on the temperature. Ultimately, this behaviour is attributed to the strongly divertor-localized radiation of N. These results demonstrate that N is more effective than Ne in achieving radiative divertor detachment while maintaining low core contamination in KSTAR, consistent with observations in other present tokamaks.

KSTAR↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

Boron Coordination in Multicomponent Glasses: Analytical Models and Machine Learning With Uncertainty

Borosilicate glasses are extensively used in a variety of applications from kitchenware to nuclear waste immobilization due to the strong network formed by the Si-O-B bond that makes it resistant to chemical corrosion and gives it a low thermal expansion. Boron, however, exists in both trigonal BO3 and tetrahedral BO4 bonds in glass systems, which impacts the chemical durability and thermal resistance of the glass, amongst other properties. Boron coordination (N4), or the ratio of the amount of BO4 to BO3 within a glass, may aid in predicting these properties but is difficult to derive without experimental data due to the complexity of impacts from varied glass compositions and processing factors. For this reason, compositional models have been developed to predict boron coordination, but the models typically include a limited number of glass components. To help fill this gap in the models, in this work, a diverse multicomponent glass dataset of 809 glasses is compiled from a literature search, and then a number of analytical and machine learning (ML) models are trained on the dataset. Previously developed modified Bernstein and modified Du Stebbins analytical models were fitted to update parameters with the new dataset. Then, partially Bayesian neural networks, Gaussian process regressor, and heteroskedastic deterministic neural networks were evaluated. The ML models examined all have different strategies to overcome the potential for overfitting as a result of a limited training dataset, and return results that account for model uncertainty, which can be valuable for understanding model reliability. For the first time, cooling rate is introduced as an input parameter for ML models, showing consistent improvements in performance and solidifying the importance of including parameters outside of composition alone for N4 prediction. The machine learning models examined here show promise in accurate predictions of boron coordination in borosilicate glasses, all achieving R2 values of 0.91.

boron coordination↗