SEARCH · Search NASA

Results for “Order parameters”

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.

332 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Adsorptive behavior of poly (vinylidene fluoride) membranes for the recovery of lignin-derived hydrophobic deep eutectic solvents

Abstract Recently, membrane technology has gained significant traction as an energy-efficient alternative to traditional thermal processes for solvent recovery. Deep eutectic solvents (DESs) have emerged as sustainable alternatives to conventional organic solvents, yet a systematic methodology for selecting compatible membrane materials for their recovery remains underdeveloped. This study established a predictive framework for membrane material selection in hydrophobic DES applications using Hansen Solubility Parameters (HSP) with inverted criteria targeting materials with relative energy difference (RED) values greater than 1.0. Flat sheet membranes were fabricated via the non-solvent induced phase separation (NIPS) technique. Four NIPS fabricated polymer membranes were evaluated: polysulfone, cellulose acetate, polyvinylidene fluoride (PVDF) fabricated with polyethylene glycol (PEG) as a pore-forming agent, and polybenzimidazole (PBI). The HSP approach successfully predicted membrane-solvent compatibility, with polysulfone (RED = 0.6) and cellulose acetate (RED = 0.9) dissolving completely within 24 h, while PVDF (RED = 1.9) and PBI (RED = 1.1) maintained structural integrity throughout a 7-day exposure period. Furthermore, PVDF demonstrated superior performance with minimal weight gain (3.0%), hydrophobic surface characteristics (122° water contact angle), and enhanced mechanical properties following DES exposure. Comprehensive chemical and morphological characterization confirmed PVDF’s chemical stability and revealed a surface-selective interaction mechanism involving simultaneous PEG (pore-forming agent) extraction and DES component adsorption. Adsorption kinetics followed pseudo-first-order behavior with reversible characteristics, best described by the Temkin isotherm model (R² = 0.9987). PVDF membranes-maintained separation functionality with average lignin rejection (75.2 ± 7.69%) and demonstrated filtration permeability of 2.0 ± 0.34 LMH/bar. This methodology provides a rational approach for membrane selection in emerging solvent systems, contributing to the advancement of sustainable separation technologies for DES-based biomass processing applications.

Science & Technology - Other Topics↗

DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗

DOE OSTI2026

Growth of deuterium supersaturated surface layer with increasing ion flux and fluence in plasma-exposed tungsten

Deuterium supersaturated surface layer (DSSL) in tungsten, a few nm thick layer exhibiting extremely high deuterium content (> 5%), has been studied as a function of deuterium plasma ion flux and fluence. Tungsten samples were exposed at 400 K and deuterium ion energy of ∼ 60 eV. The deuterium ion flux spanned over an order of magnitude, being 7.3 × 10 20 , 4.2 × 10 21 , and 3.8 × 10 22 D/m 2 . The fluence ranged from 3.4 × 1024 to 3.4 × 1025 D/m2. The samples and their deuterium content were analyzed by nuclear reaction analysis (NRA), scanning transmission electron microscopy (STEM), and thermal desorption spectroscopy (TDS). The largest thickness of DSSL was found to be around 10.5 nm and was observed in the case of the highest exposure flux and fluence, whereas the layer was only about 3.8 nm thick in the case of the lowest value of the two parameters. The thickness of the DSSL was found to monotonically increase with both deuterium flux and fluence. Finally, similar to the thickness, the estimated D concentration seems to follow the same trend, increasing from the lowest value of around 3 at.% to the highest value of around 5 at.%.

Deuterium↗

DOE OSTI2023

Design and Synthesis of Cubic K 3−2 x Ba x SbSe 4 Solid Electrolytes for K–O 2 Batteries

Developing K-ion conducting solid-state electrolytes (SSEs) plays a critical role in the safe implementation of potassium batteries. In this work, a chalcogenide-based potassium ion SSE is reported, K 3 SbSe 4 , which adopts a trigonal structure at room temperature. Single-crystal structural analysis reveals a trigonal-to-cubic phase transition at the low temperature of 50 °C, which is the lowest among similar compounds and thus provides easy access to the cubic phase. The substitution of barium for potassium in K 3 SbSe 4 leads to the creation of potassium vacancies, expansion of lattice parameters, and a transformation from a trigonal phase to a cubic phase. As a result, the maximum conductivity of K 3−2x Ba x SbSe 4 reaches around 0.1 mS cm −1 at 40 °C for K 2.2 Ba 0.4 SbSe 4 , which is over two orders of magnitude higher than that of undoped K 3 SbSe 4 . This novel SSE is successfully employed in a K–O 2 battery operating at room temperature where a polymer-laminated K 2.2 Ba 0.4 SbSe 4 pellet serves as a separator between the oxygen cathode and the potassium metal anode. Effective protection of the K metal anode against corrosion caused by O 2 is demonstrated.

25 ENERGY STORAGE↗

DOE OSTI2026

Anion-dependent phase behavior of methylimidazolium-based ionic liquids mixed with water: Correlation between local molecular structure and mesoscale behaviors

The mesoscopic phase behavior of decylmethylimidazolium (C 10 mim) ionic liquids (ILs) bearing three monovalent anions—thiocyanate (SCN − ), nitrate (NO 3 − ), and chloride (Cl − )—mixed with water at relatively high IL contents (50–95 wt%) was investigated. Small-angle and wide-angle X-ray scattering (SAXS/WAXS) were employed to follow the evolution of both local and mesoscale structures across this composition range. In the absence of water, C 10 mimSCN and C 10 mimCl behaved as disordered liquids, whereas C 10 mimNO 3 spontaneously formed a hexagonally ordered cylindrical mesophase and displayed a sticky-solid macroscopic appearance. Upon addition of water, C 10 mimSCN remained a viscous liquid and only weakly ordered lamellar domains were observed. This limited ordering is attributed to the lack of hydrogen bonding and weak interaction energy between SCN − anions and C 10 mim + cations. In contrast, the trigonal-planar NO 3 − and point-like Cl − anions promoted the formation of well-defined hexagonal mesophases up to 35–45 wt% water. The formation of hydrogen bonding of the two anions with imidazolium ring protons likely enabled the creation of compact ion clusters that effectively exclude water molecules from the immediate vicinity of the IL aggregates. These findings demonstrate that the interaction energy between ion pairs dictate IL-water interactions and therefore control the transition from disordered liquids to ordered mesophases in IL/water mixtures with high IL contents. In conclusion, the combined SAXS/WAXS analysis reveals a correlation between local intermolecular structure and the emergence of mesoscopic order, providing a systematic framework for tailoring mesoscale structures in alkylimidazolium-based IL/water systems.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Coexistence of Synchronization and Stochasticity in Thermally Coupled Mott Oscillators

Synchronization is conventionally regarded as a mechanism for suppressing variability and enforcing order in coupled systems, from pendula and lasers to neurons and electronic oscillators. Here, we show that synchronization can also embed stochasticity at finer scales. We observe this phenomenon in thermally coupled VO 2 neuristors, where robust in-phase synchronization at the microsecond scale coexists with spike onset fluctuations at the nanosecond scale, with no fixed leader. The coexistence of order and disorder originates from stochastic domain-level physics of the insulator–metal and metal–insulator transitions, where local variations in transition temperature drive cycle-to-cycle randomness in nucleation, percolation, and relaxation. A stochastic domain model reproduces this effect by generating synchronized spike trains with random lead–lag jitter, and experimental interspike interval statistics confirm the persistence of fine-scale variability despite macroscopic phase locking. These findings establish that synchronization and stochasticity can coexist within the same physical platform, revealing hidden disorder within collective order. Furthermore, this insight reframes synchronization as not purely deterministic, but as a universal context where microscopic variability can persist, with implications for electronics, cryptography, and the fundamental physics of order–disorder coexistence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Decorated Honeycomb Lattice and Successive Magnetic Transitions of the Delafossite Derivative K 4 Ni 5 Te 3 O 16

Here, we report the synthesis and characterization of the delafossite derivative K 4 Ni 5 Te 3 O 16 , which hosts a decorated honeycomb Ni 2+ lattice. Single-crystal X-ray diffraction, electron microscopy, and neutron diffraction establish a noncentrosymmetric Pmm2 crystal structure that combines regular honeycomb units with elongated motifs formed by Ni 2+ trimers along extended edges. Magnetic susceptibility reveals strong antiferromagnetic interactions, while heat capacity measurements identify two successive magnetic transitions at 30 and 12 K. Neutron diffraction shows that these transitions correspond to a progression from partial to full long-range magnetic order on distinct Ni sites, reflecting competition among the various exchange interactions. The resulting magnetic point group, mm2.1′, inherits the lattice polarity and permits high-order magnetoelectric coupling. Consistent with this symmetry, a quadratic magnetic-field-induced polarization is experimentally observed below the magnetic ordering temperature, likely arising from exchange-driven magnetostriction coupled to the polar lattice. These results establish K 4 Ni 5 Te 3 O 16 as an interesting platform for engineering B-site ordering in delafossite derivatives to realize complex quantum magnetism and functional responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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 OSTI2025

A topological superconductor tuned by electronic correlations

A topological superconductor, characterized by either a chiral order parameter or a topological surface state in proximity to bulk superconductivity, is foundational to topological quantum computing. A key open challenge is whether electron-electron interactions can tune such topological superconducting phases. Here, we provide experimental signatures of a unique topological superconducting phase in competition with electronic correlations in 10-unit-cell thick FeTe x Se 1-x films grown on SrTiO 3 substrates. When the Te content x exceeds 0.7, we observe a topological transition marked by the emergence of a superconducting surface state. Near the FeTe limit, the system undergoes another transition where the surface state disappears, and superconductivity is suppressed. Theory suggests that electron-electron interactions in the odd-parity xy− band drives this second topological transition. The flattening and eventual decoherence of d xy -derived bands track the superconducting dome, linking correlation effects directly to superconducting coherent transport. Our work establishes many-body electronic correlations as a sensitive knob for tuning topology and superconductivity, offering a pathway to engineer new topological phases in correlated materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics↗

DOE OSTI2023

Host galaxies of ultra-strong Mg ii absorbers at z ∼ 0.7

ABSTRACT We report spectroscopic identification of the host galaxies of 18 ultra-strong Mg ii systems (USMg ii) at 0.6 ≤ z ≤ 0.8. We created the largest sample by merging these with 20 host galaxies from our previous survey within 0.4 ≤ z ≤ 0.6. Using this sample, we confirm that the measured impact parameters ($\rm 6.3\leqslant D[kpc] \leqslant 120$ with a median of 19 kpc) are much larger than expected, and the USMg ii host galaxies do not follow the canonical $\rm {\it W}_{2796}-{\it D}$ anticorrelation. We show that the presence and significance of this anticorrelation may depend on the sample selection. The $\rm {\it W}_{2796}-{\it D}$ anticorrelation seen for the general Mg ii absorbers show a mild evolution at low $\rm W_{2796}$ end over the redshift range 0.4 ≤ z ≤ 1.5 with an increase of the impact parameters. Compared to the host galaxies of normal Mg ii absorbers, USMg ii host galaxies are brighter and more massive for a given impact parameter. While the USMg ii systems preferentially pick star-forming galaxies, they exhibit slightly lower ongoing star-forming rates compared to main sequence galaxies with the same stellar mass, suggesting a transition from star-forming to quiescent states. For a limiting magnitude of mr < 23.6, at least 29 per cent of the USMg ii host galaxies are isolated, and the width of the Mg ii absorption in these cases may originate from gas flows (infall/outflow) in isolated haloes of massive star forming but not starbursting galaxies. We associate more than one galaxy with the absorber in $\ge 21~{{\ \rm per\ cent}}$ cases, where interactions may cause wide velocity spread.

Astronomy & Astrophysics↗

DOE OSTI2025

The role of fast and slow dynamics in nonlinear resonant ultrasound spectroscopy of consolidated granular materials

Abstract Elastic nonlinearity observed in consolidated granular media can be attributed to the combination of slow and fast effects, which give rise to hysteresis and relaxation of both modulus and damping after the sample is perturbed. A consequence is a high level of complexity in the measurements of the sample linear and nonlinear elastic parameters. The results of experiments are dependent on the experimental protocol that is adopted to measure the relevant quantities and it is hard to quantify parameters with accuracy and repeatability. Here we focus on examining Nonlinear Resonant Ultrasound Spectroscopy, showing experimentally the role of slow dynamics in the process and quantifying/discussing its influence on the quantification of nonlinearity. We also propose a model to describe the process, which shows that different contributions to nonlinearity (e.g., classical and hysteretic) could be due to physical features (defects) relaxing with different relaxation times.

Science & Technology - Other Topics↗

DOE OSTI2026

Open quantum system approach to inclusive jet production in heavy-ion collisions

We derive a factorization formula for inclusive jet production in heavy-ion collisions using the tools of Effective Field Theory (EFT). We show how physics at widely separated scales in this process can be systematically separated by matching to EFTs at successively lower virtualities. Owing to a strong scale separation, we recover a vacuum-like DGLAP evolution above the jet scale, while the additional low-energy scales induced by the medium effectively probe the internal structure of the jet. As a result, the cross section can be written as a series with an increasing number of subjets characterized by perturbative matching coefficients each of which is convolved with a distinct function. These functions encode broadening, medium-induced radiations as well as quantum interference such as the Landau-Pomeranchuk-Migdal effect and color coherence dynamics to all orders in perturbation theory. As a first application of this EFT framework, we investigate the case of an unresolved jet and show how the cross section can be factorized and fully separate the jet dynamics from the universal physics of the medium. To compare to the existing literature, we explicitly compute the medium jet function at next-to-leading order in the coupling and leading order in medium opacity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Manipulation of Localized Excitons in CrPS4 by Temperature and Magnetic Field

Layered van der Waals magnetic semiconductors provide a versatile platform for exploring excitonic phenomena intertwined with spin and lattice degrees of freedom, enabling excitons to act as sensitive probes of magnetic order. CrPS4 is a layered antiferromagnetic semiconductor that hosts rich excitonic features whose microscopic origin and connection to magnetic ordering remain incompletely understood. Here, we investigate the electronic and excitonic properties of bulk CrPS4 using a combination of many-body perturbation theory, dynamical mean-field theory, and photoluminescence-based experiments. Our calculations establish CrPS4 as a direct-gap semiconductor with a bandgap of 2.48 eV in the antiferromagnetic phase. Several subbandgap excitonic transitions are predicted by theory, comprising multiple spin-allowed excitons and an additional spin-flip excitation, predominantly localized on the Cr3+ ions. Temperature- and magnetic-field-dependent optical measurements reveal thermally driven exciton redistribution among localized states and identify characteristic energy shifts that provide clear optical signatures of magnetic phase transitions in CrPS4. These results provide insights into the excitonic transitions of antiferromagnets and suggest potential routes for all-optical sensing and light-driven control of their magnetic order.

2D materials↗

DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗

DOE OSTI2026

Assembly bias and local Primordial non-Gaussianity from DESI DR1 quasars

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by f NL . The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on f NL , however, are fully degenerate with the local PNG bias parameter b ϕ , dependent on the assembly bias parameter p. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for p, and thus for b ϕ , building on the findings of Fondi et al. 2024. We find a distribution with mean p̅ ≃ 1.4 with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining f NL = -3.3±9.2.

cosmological parameters from LSS↗

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↗