Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “range prediction”

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.

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Symmetric U(1) and Z 2 spin liquids on the pyrochlore lattice

The geometrically frustrated three-dimensional pyrochlore lattice has been long predicted to host a quantum spin liquid, an intrinsic long-range entangled state with fractionalized excitations. To date, most proposals for pyrochlore materials have focused on quantum spin ice, a U(1) quantum spin liquid whose only low energy excitations are emergent photons of Maxwell type. Here, we explore the possibility of finding pyrochlore quantum spin liquids whose low energy theories go beyond this standard one. We give a complete classification of symmetric U ( 1 ) and Z 2 spin liquids on the pyrochlore lattice within the projective symmetry group framework for fermionic spinons. We find 18 U ( 1 ) spin liquids and 28 Z 2 spin liquids that preserve pyrochlore space-group symmetry while, upon further imposing time-reversal symmetry, the numbers of classes become 16 and 48, respectively. For each class, the most general symmetry-allowed spinon mean-field Hamiltonian is given. Interestingly, we find that several U(1) spin liquid classes possess an unusual gapless multi-nodal-line structure (“nodal star”) in the spinon bands, which is protected by the projective actions of the threefold rotation and screw symmetries of the pyrochlore space group. Through a simple model, we study the effect of gauge fluctuations on such a nodal star spin liquid and propose that the leading terms in the low temperature specific heat have the scaling form C / T ~ T + T / ln T , in contrast to the form C / T ~ T 2 of the standard U(1) pyrochlore spin liquid with gapped spinons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Measurement of the nucleon spin structure functions for 0.01 < 𝑄 2 < 1 GeV 2 using CLAS

The spin structure functions of the proton and the deuteron were measured during the EG4 experiment at Jefferson Lab in 2006. Data were collected for longitudinally polarized electron scattering off longitudinally polarized NH 3 and ND 3 targets, for 𝑄 2 values as small as 0.012 and 0.02 GeV 2 , respectively, using the CEBAF Large Acceptance Spectrometer. This is the archival paper of the EG4 experiment that summarizes the previously reported results of the polarized structure functions 𝑔 1 , 𝐴 1 ⁢𝐹 1 , and their moments $\bar{Γ}$ 1 , $\bar{𝛾}$ 0 , and $\bar{𝐼}$ TT , for both the proton and the deuteron. In addition, we report on new results on the neutron 𝑔 1 extracted by combining proton and deuteron data and correcting for Fermi smearing, and on the neutron moments $\bar{Γ}$ 1 , $\bar{𝛾}$ 0 , and $\bar{𝐼}$ TT formed directly from those of the proton and the deuteron. Our data are in good agreement with the Gerasimov-Drell-Hearn sum rule for the proton, deuteron, and neutron. Furthermore, the isovector combination was formed for 𝑔 1 and the Bjorken integral $\bar{Γ}^{𝑝−𝑛}_1$, and it was compared to available theoretical predictions. All of our results, to the best of our knowledge, provide for the first time extensive tests of spin observable predictions from chiral effective field theory (𝜒⁢EFT) in a 𝑄 2 range commensurate with the pion mass. Finally, they motivate further improvement in 𝜒⁢EFT calculations from other approaches such as the lattice gauge method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

High-Statistics Measurement of Neutrino Quasielasticlike Scattering at 6 GeV on a Hydrocarbon Target

We measure neutrino charged-current quasielasticlike scattering on hydrocarbon at high statistics using the wideband Neutrinos at the Main Injector beam with neutrino energy peaked at 6 GeV. The double-differential cross section is reported in terms of muon longitudinal ( p ∥ ) and transverse ( p ⊥ ) momentum. Cross section contours versus lepton momentum components are approximately described by a conventional generator-based simulation, however, discrepancies are observed for transverse momenta above 0.5 GeV / c for longitudinal momentum ranges 3–5 and 9 – 20 GeV / c . The single differential cross section versus momentum transfer squared ( d σ / d Q QE 2 ) is measured over a four-decade range of Q 2 that extends to 10 GeV 2 . The cross section turnover and falloff in the Q 2 range 0.3 – 10 GeV 2 is not fully reproduced by generator predictions that rely on dipole form factors. Our measurement probes the axial-vector content of the hadronic current and complements the electromagnetic form factor data obtained using electron-nucleon elastic scattering. These results help oscillation experiments because they probe the importance of various correlations and final-state interaction effects within the nucleus, which have different effects on the visible energy in detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

X-Ray Diffraction of Ramp-Compressed Silicon to 390 GPa

Silicon (Si) exhibits a rich collection of phase transitions under ambient-temperature isothermal and shock compression. This report describes in-situ diffraction measurements of ramp compressed Si between 40 and 389 GPa. Angle-dispersive x-ray scattering reveals that Si assumes an hexagonal close-packed (hcp) structure between 40 and 93 GPa, and, at higher pressure, a face-centered cubic (fcc) structure which persists to at least 389 GPa, the highest pressure for which the crystal structure of Si has been investigated. Furthermore, the range of hcp stability extends to higher pressures and temperatures than predicted by theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Finite-Time Teleportation Phase Transition in Random Quantum Circuits

How long does it take to entangle two distant qubits in a quantum circuit evolved by generic unitary dynamics? Here, we show that if the time evolution is followed by measurements of all but two infinitely separated test qubits, then the entanglement between them can undergo a phase transition and become nonzero at a finite critical time t c . The fidelity of teleporting a quantum state from an input qubit to an infinitely distant output qubit shows the same critical onset. Specifically, these finite-time transitions occur in short-range interacting two-dimensional random unitary circuits and in sufficiently long-range interacting one-dimensional circuits. The phase transition is understood by mapping the random continuous-time evolution to a finite-temperature thermal state of an effective spin Hamiltonian, where the inverse temperature equals the evolution time in the circuit. In this framework, the entanglement between two distant qubits at times t > t c corresponds to the emergence of long-range ferromagnetic spin correlations below the critical temperature. We verify these predictions using numerical simulation of Clifford circuits and propose potential realizations in existing platforms for quantum simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tallo: A global tree allometry and crown architecture database

Abstract Data capturing multiple axes of tree size and shape, such as a tree's stem diameter, height and crown size, underpin a wide range of ecological research—from developing and testing theory on forest structure and dynamics, to estimating forest carbon stocks and their uncertainties, and integrating remote sensing imagery into forest monitoring programmes. However, these data can be surprisingly hard to come by, particularly for certain regions of the world and for specific taxonomic groups, posing a real barrier to progress in these fields. To overcome this challenge, we developed the Tallo database, a collection of 498,838 georeferenced and taxonomically standardized records of individual trees for which stem diameter, height and/or crown radius have been measured. These data were collected at 61,856 globally distributed sites, spanning all major forested and non‐forested biomes. The majority of trees in the database are identified to species (88%), and collectively Tallo includes data for 5163 species distributed across 1453 genera and 187 plant families. The database is publicly archived under a CC‐BY 4.0 licence and can be access from: https://doi.org/10.5281/zenodo.6637599 . To demonstrate its value, here we present three case studies that highlight how the Tallo database can be used to address a range of theoretical and applied questions in ecology—from testing the predictions of metabolic scaling theory, to exploring the limits of tree allometric plasticity along environmental gradients and modelling global variation in maximum attainable tree height. In doing so, we provide a key resource for field ecologists, remote sensing researchers and the modelling community working together to better understand the role that trees play in regulating the terrestrial carbon cycle.

Jucker, Tommaso↗
DOE OSTI2025

Simplified Model and Approach to Transform Infrared Surface Temperature to Film Effectiveness in a Conjugate Heat Transfer Experiment

In the pursuit of more efficient gas turbines, film cooling is a critical technology. This article describes a simplified engineering model based on a one-dimensional thermal resistance network. The model is used to relate film-cooling effectiveness and heat transfer augmentation to local overall cooling effectiveness in a conjugate flat plate experiment. Here, this article presents experimental proof-of-concept data to demonstrate the potential for this model. In contrast to previous approaches, neither the wall heat flux nor the adiabatic wall temperature is required to estimate the local film-cooling performance parameters. The model predicts surface temperatures that are within the experimental uncertainties over the range for which the model is trained and to within five percent when the model is extrapolated to higher coolant channel Reynolds numbers. This article is relevant to conjugate test rigs that can measure the hot-surface temperature distribution with and without film cooling. This information may also be relevant to designers as a method to approximate surface temperatures or used as an approximate heat transfer model for optimization studies.

advanced turbines↗
DOE OSTI2020

Control of polymorphism during epitaxial growth of hyperferroelectric candidate LiZnSb on GaSb (111)B

A major challenge for ferroelectric devices is the depolarization field, which competes with and often destroys long-range polar order in the limit of ultrathin films. Recent theoretical predictions suggest a new class of materials, termed hyperferroelectics, that should be robust against the depolarization field and enable ferroelectricity down to the monolayer limit. Here, the authors demonstrate the epitaxial growth of hexagonal LiZnSb, one of the hyperferroelectric candidate materials, by molecular-beam epitaxy on GaSb (111)B substrates. Due to the high volatility of all three atomic species, they find that LiZnSb can be grown in an adsorption-controlled window, using an excess zinc flux. Within this window, the desired polar hexagonal phase is stabilized with respect to a competing cubic polymorph, as revealed by x-ray diffraction and transmission electron microscopy measurements. First-principles calculations show that for moderate amounts of epitaxial strain and moderate concentrations of Li vacancies, the cubic LiZnSb phase is lower in formation energy than the hexagonal phase, but only by a few millielectronvolts per formula unit. Therefore, they suggest that kinetics plays a role in stabilizing the desired hexagonal phase at low temperatures. Their results provide a path toward experimentally demonstrating ferroelectricity and hyperferroelectricity in a new class of ternary intermetallic compounds.

42 ENGINEERING↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗
DOE OSTI2021

Watching a hydroperoxyalkyl radical (•QOOH) dissociate

A prototypical hydroperoxyalkyl radical (•QOOH) intermediate, transiently formed in the oxidation of volatile organic compounds, was directly observed through its infrared fingerprint and energy-dependent unimolecular decay to hydroxyl radical and cyclic ether products. Direct time-domain measurements of •QOOH unimolecular dissociation rates over a wide range of energies were found to be in accord with those predicted theoretically using state-of-the-art electronic structure characterizations of the transition state barrier region. Unimolecular decay was enhanced by substantial heavy-atom tunneling involving O-O elongation and C-C-O angle contraction along the reaction pathway. Master equation modeling yielded a fully a priori prediction of the pressure-dependent thermal unimolecular dissociation rates for the •QOOH intermediate—again increased by heavy-atom tunneling—which are required for global models of atmospheric and combustion chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Highly diverse and unknown viruses may enhance Antarctic endoliths’ adaptability

Abstract Background Rock-dwelling microorganisms are key players in ecosystem functioning of Antarctic ice free-areas. Yet, little is known about their diversity and ecology, and further still, viruses in these communities have been largely unexplored despite important roles related to host metabolism and nutrient cycling. To begin to address this, we present a large-scale viral catalog from Antarctic rock microbial communities. Results We performed metagenomic analyses on rocks from across Antarctica representing a broad range of environmental and spatial conditions, and which resulted in a predicted viral catalog comprising > 75,000 viral operational taxonomic units (vOTUS). We found largely undescribed, highly diverse and spatially structured virus communities which had predicted auxiliary metabolic genes (AMGs) with functions indicating that they may be potentially influencing bacterial adaptation and biogeochemistry. Conclusion This catalog lays the foundation for expanding knowledge of virosphere diversity, function, spatial ecology, and dynamics in extreme environments. This work serves as a step towards exploring adaptability of microbial communities in the face of a changing climate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Microscopic subsurface characterization of layered magnetic materials using magnetic resonance force microscopy

This project resulted in the following accomplishments: 1) The invention of scanned probe ferromagnetic resonance microscopy. Upon implementation and refinement this approach was used to measure static and dynamic magnetic resonance properties of ferromagnets with an unprecedented combination of spectral and spatial resolution. 2) Demonstration and study of spin pumping in a pristine ferromagnet across a boundary defined entirely by a magnetic field gradient, and so free of materials discontinuities. 3) First demonstrated highly efficient dynamic spin injection from Y 3 Fe 5 O 12 (YIG) into NiO, an antiferromagnetic (AF) insulator and robust spin propagation in AF NiO over distances up to 100-nm thickness; both spin injection and spin transport are mediated by the AF spin correlations. 4) Demonstration of exceptionally high efficiency conversion of spin current to charge current achieved in spin pumping from the dynamic magnetization of a ferromagnetic insulator into a variety of adjacent metal films. 5) Seminal study of the nature of spin-orbit coupling in d-electron metals that revealed the dominant sensitivity of the d-orbital moment to the d-electron count and revealing the important roles of both d-electron configuration and atomic number in the physics of the spin Hall effect. 6) Demonstration of observation of FMR by means of optically-detected NV diamond relaxometry that provided a highly sensitive, spatially resolved FMR detection method and revealed spin pumping across 100-nm length scales by means of dipolar coupling. 7) Using magnetic field localization, developed for FMR imaging, demonstrated a novel approach to discretizing spin wave modes to improve the characteristics of auto-oscillators driven by pure spin currents produced by the spin Hall effect. 8) We demonstrated large spin pumping signals from propagating high wavevector spinwaves in thin film YIG find that the spin pumping efficiency was insensitive to wavevector thus emphasizing the promise of thin YIG films for spinwave-based spintronics by allowing for predictable down-scaling of device dimensions and enabling a broad range of operational frequencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hadron spectroscopy of pseudoscalar and vector meson photoproduction using linearly polarised photons at CLAS

The field of hadron spectroscopy is informed by experiments which study excited states of the nucleon by detecting the decay products following excitation with a hadronic, leptonic or electromagnetic probe. Experiments using a polarised beam of photons provide an important contribution to the world data. This thesis describes the analysis of data from one such experiment, the CEBAF Large Acceptance Spectrometer at Jefferson Laboratory, Virginia, USA. In the experiment, an incident beam of electrons passes through a diamond radiator producing a linearly polarised photon beam of energy from 1.1 to 2.1 GeV on a liquid hydrogen target. The analysis method and results for two reaction channels are presented. In the first analysis, the polarisation observables {?, P, T, Ox, Oz} for the reaction ??? p ? K0 S ?+ are extracted simultaneously using likelihood sampling with Markov-Chain Monte-Carlo. The results are presented as angular distributions (cos ?K0S) for bins in photon energy E? for a total of 21 bins. The values extracted for T , Ox and Oz are a first measurement of these quantities for this reaction and the beam asymmetry (?) measurement supplements the one previous measurement for this quantity and extends the energy range over which it has been extracted. The measurement of the recoil polarisation, P, is consistent with previous measurements. In a second analysis, the spin density matrix elements (SDMEs) for ? p ? p ? are extracted using similar analysis techniques. The results are presented as angular distributions (cos ?? ) for bins in photon energy E? for a total of 24 bins. This work is the first to extract all nine SDMEs for the full angular range of the ? meson production. A comparison is made to the predictions of the Vector Meson Dominance (VMD) model and the results provide evidence for contributing processes other than VMD, particularly at more backward angles of the ? meson production. The data also indicate non-helicity conserving process in both the s-channel and the t-channel. The results from the first analysis have been passed to the Juelich-Bonn theory group and the effect of including the new data within their dynamical coupled channel model is described. Several nucleon resonances are affected and the impact on the pole positions, widths and photon decay amplitudes of the affected resonances is described.

Clark, Louise↗
DOE OSTI2025

Simplified Model and Approach to Transform Infrared Surface Temperature to Film Effectiveness in a Conjugate Heat Transfer Experiment

This paper describes a simplified engineering model based on a one-dimensional thermal resistance network. The model is used to develop a new method to relate film cooling effectiveness and heat transfer augmentation to local overall cooling effectiveness in a conjugate flat plate experiment. This paper presents experimental proof-of-concept data to demonstrate the potential for this model. In contrast to previous approaches, neither the wall heat flux nor the adiabatic wall temperature is required to estimate the local film cooling performance parameters. The model predicts surface temperatures that are within the experimental uncertainties over the range for which the model is trained, and to within five percent when the model is extrapolated to higher coolant channel Reynolds numbers. This paper is relevant to conjugate test rigs that can measure the hot surface temperature distribution with and without film cooling. This information may also be relevant to designers as a method to approximate surface temperatures or used as an approximate heat transfer model for optimization studies.

advanced gas turbines↗
DOE OSTI2025

LDRD conclusion poster - Synthesizing Heterometallic Uranium Single Crystals to Understand the Influence of the Secondary Metals on Uranyl Axial Bond Strength

Understanding how transition metals influenced the chemistry of lanthanide and actinide (f-element) materials is critical for advancing separation technologies, materials design, and coordination chemistry. This project examined how incorporating first-row transition metals affected the structural and spectroscopic properties of f-element coordination polymers. In uranium(VI)-based systems synthesized with 2,6-pyridinedicarboxylic acid (PDC) ligands, single-crystal X-ray diffraction and Raman spectroscopy revealed that the presence of transition metals shortened the uranyl axial bond and induced a blue shift in its symmetric stretching vibration—evidence of increased bond strength. Electronic structure analysis, including Density of States (DOS) calculations using density functional theory (DFT), revealed altered orbital overlaps and highlighted the role of transition metal d-orbitals in modulating bonding. Raman modes were modeled using truncated structural fragments in collaboration with the University of Notre Dame, and although the predicted frequencies were lower than experimental values, they remained within expected ranges. In parallel, similar experiments with cerium (Ce) in the presence of cobalt (Co) and PDC demonstrated multi-step single-crystal-to-single-crystal transformations—behavior not observed in the uranium systems. Initial products included light yellow, orange, and polycrystalline materials. Single-crystal X-ray diffraction studies, conducted in collaboration with the Colorado School of Mines, identified the yellow phase as monometallic Ce(PDC)2(H2O)2·4H2O and the orange phase as heterometallic Ce2Co(PDC)4(H2O)6. After standing in solution for one week, both phases fully transformed into a dark yellow crystalline phase, [Ce3(PDC)5(H2O)8].6(H2O). Remarkably, this transformation was reversible—disturbing the equilibrium by removing some crystals caused reversion to the initial Ce(PDC)2(H2O)2·4H2O phase, highlighting dynamic behavior. All three structures were previously unreported. Solid-state UV-visible and Raman spectroscopy further distinguished these phases, revealing ligand-to-metal charge transfer involving Ce and characteristic d–d transitions from Co(II). The precise mechanism driving these transformations remained unclear; however, pH-dependent experiments confirmed that the transformation did not occur when the pH decreased. Overall, the project demonstrated that transition metals could be employed to tune bonding interactions, structural dimensionality, and optical properties in f-element materials, establishing new pathways for designing functional heterometallic systems. The work resulted in several novel structural discoveries and fostered productive collaborations with the University of Notre Dame and the Colorado School of Mines.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2022

Further Evidence of Modified Spin-down in Sun-like Stars: Pileups in the Temperature–Period Distribution

We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature–period distribution for main-sequence stars with temperatures exceeding ~5500 K. The long-period pileup is well described by a curve of constant Rossby number, with a critical value of Ro crit ≲ Ro ⊙ . The long-period pileup was predicted by van Saders et al. as a consequence of weakened magnetic braking, in which wind-driven angular momentum losses cease once stars reach a critical Rossby number. Stars in the long-period pileup are found to have a wide range of ages (~2–6 Gyr), meaning that, along the pileup, rotation period is strongly predictive of a star's surface temperature but weakly predictive of its age. The short-period pileup, which is also well described by a curve of constant Rossby number, is not a prediction of the weakened magnetic braking hypothesis but may instead be related to a phase of slowed surface spin-down due to core-envelope coupling. The same mechanism was proposed by Curtis et al. to explain the overlapping rotation sequences of low-mass members of differently aged open clusters. The relative dearth of stars with intermediate rotation periods between the short- and long-period pileups is also well described by a curve of constant Rossby number, which aligns with the period gap initially discovered by McQuillan et al. in M-type stars. These observations provide further support for the hypothesis that the period gap is due to stellar astrophysics, rather than a nonuniform star formation history in the Kepler field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗