Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduce order modeling”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The MicroBooNE Single-Photon Low-Energy Excess Search (Public Note 1087)

MicroBooNE is a short baseline neutrino experiment at Fermilab designed to address the low energy excess observed by the MiniBooNE experiment. This note describes and presents preliminary results for the MicroBooNE analysis developed to address this excess as a single photon plus one or zero protons in the final state. The analysis assumes neutrino neutral current Δ resonance production followed by Δ radiative decay on argon (NC Δ → Nγ ) as the "signal model"; event reconstruction and selection have been developed and optimized in order to maximize efficiency and reduce cosmogenic and other beam-related backgrounds to the NC Δ → Nγ signal. We present the analysis methodology and validation checks performed on limited-statistics open data sets, corresponding to 5 x10 19 protons on target (POT), following a blind analysis, as well as the projected sensitivities for testing the Standard Model (SM) predicted rate for the NC Δ → Nγ process and for testing the interpretation of the previously observed MiniBooNE low energy excess as NC Δ → Nγ events, using the full anticipated MicroBooNE data set of 12.25 x 10 20 POT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The MicroBooNE Single-Photon Low-Energy Excess Search

MicroBooNE is a short baseline neutrino experiment at Fermilab designed to address the low energy excess observed by the MiniBooNE experiment. This note describes and presents preliminary results for the MicroBooNE analysis developed to address this excess as a single photon plus one or zero protons in the final state. The analysis assumes neutrino neutral current Δ resonance production followed by Δ radiative decay on argon (NC Δ → $N_γ$) as the “signal model”; event reconstruction and selection have been developed and optimized in order to maximize efficiency and reduce cosmogenic and other beam-related backgrounds to the NC Δ → $N_γ$ signal. We present the analysis methodology and validation checks performed on limited-statistics open data sets, corresponding to 5×10 19 protons on target (POT), following a blind analysis, as well as the projected sensitivities for testing the Standard Model (SM) predicted rate for the NC Δ → $N_γ$ process and for testing the interpretation of the previously observed MiniBooNE low energy excess as NC Δ → $N_γ$ events, using the full anticipated MicroBooNE data set of 12.25×10 20 POT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

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

Stacked-mosaic amplifier and diode delivery concept for kJ-class inertial fusion energy laser driver

We report on a laser amplifier architecture designed as a modular element for scaling to multi-megajoule laser facilities intended for inertial fusion energy (IFE) power plants. This kilojoule-class module features a stacked mosaic of gain media integrated with a diode delivery system for efficient optical pumping of the mosaic-structured gain medium, enabling high-repetition-rate operation at high wall plug efficiencies, necessary for an IFE driver. Details of a diode delivery system capable of pumping the mosaic architecture are presented. Comprehensive numerical modeling demonstrates that the stacked-mosaic approach reduces transverse gain by five orders of magnitude compared to a conventional full-aperture Yb:YAG based amplifier, substantially suppressing transverse amplified spontaneous emission (TASE) and enabling enhanced longitudinal energy extraction. Detailed analysis of thermal management and wavefront distortion in a 4 × 4 mosaic array indicates that temperature gradients and thermally induced aberrations are effectively controlled using gas cooling and commercially available phase-plate and deformable mirror technologies. We further discuss the applicability of a stacked-mosaic architecture to direct-drive IFE schemes, where broad spectral bandwidth is critical, and its compatibility with frequency conversion modules for up-conversion to blue wavelengths. Finally, an example point design for a 10 kJ, 10 Hz Yb:YAG module operating at 175 K with wall-plug efficiency exceeding 10 % is presented, underscoring the feasibility of this approach for next-generation high-energy lasers for IFE drivers. The results establish the stacked-mosaic amplifier as a scalable, robust platform not only for IFE but also for a broad range of advanced scientific and industrial laser applications.

Lasers↗
DOE OSTI2023

A variational encoder–decoder approach to precise spectroscopic age estimation for large Galactic surveys

Constraints on the formation and evolution of the Milky Way Galaxy require multidimensional measurements of kinematics, abundances, and ages for a large population of stars. Ages for luminous giants, which can be seen to large distances, are an essential component of studies of the Milky Way, but they are traditionally very difficult to estimate precisely for a large data set and often require careful analysis on a star-by-star basis in asteroseismology. Because spectra are easier to obtain for large samples, being able to determine precise ages from spectra allows for large age samples to be constructed, but spectroscopic ages are often imprecise and contaminated by abundance correlations. Here we present an application of a variational encoder–decoder on cross-domain astronomical data to solve these issues. The model is trained on pairs of observations from APOGEE and Kepler of the same star in order to reduce the dimensionality of the APOGEE spectra in a latent space while removing abundance information. The low dimensional latent representation of these spectra can then be trained to predict age with just ∼1000 precise seismic ages. We demonstrate that this model produces more precise spectroscopic ages (∼ 22 per cent overall, ∼ 11 per cent for red-clump stars) than previous data-driven spectroscopic ages while being less contaminated by abundance information (in particular, our ages do not depend on [α/M]). We create a public age catalogue for the APOGEE DR17 data set and use it to map the age distribution and the age-[Fe/H]-[α/M] distribution across the radial range of the Galactic disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

Analysis of Waste Material Feedstocks Using Laser-Induced Breakdown Spectroscopy and Machine Learning

Predicting properties such as heating value, ash fusion temperature, and mineral ash composition from Laser-Induced Breakdown Spectroscopy (LIBS) data can make gasifiers more flexible to different feedstocks. Understanding these feedstock properties in-situ improves feedstock conversion modelling methods that allow for consistent operation, higher carbon conversion, and reduced fouling and erosion rates. The purpose of this study is to demonstrate methods for model creation that take LIBS data as predictor features and estimate higher order material properties as a function of feedstock material properties. Six samples were chosen to represent a mixture of abundant and carbon rich waste materials. LIBS measurements were performed on these samples for elemental wavelengths and intensity values. Laboratory analytical results were obtained for each sample’s heating value, proximate and ultimate analysis, mineral ash composition, ash fusion temperatures, and viscosity temperatures. Thermal conductivity was measured using a HotDisk TPS 2500S. LIBS measurements were processed and used as predictor features for machine learning (ML) models to predict the sample’s material properties. Predictor feature selection algorithms, particularly minimum redundancy maximum relevance (mRMR), reduced the dimensionality of ML models. Many modelling methods such as Gaussian process regression (GPR), regression tree, neural networks (NN), and support vector machines (SVM) were demonstrated to be effective at predicting higher order properties; however, mRMR with GPR stood out as a clear winning combination.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Implementation of a Model Predictive Control Strategy to Regulate Temperature Inside Plug-Flow Solar Reactor With Countercurrent Flow

Abstract Solar-driven thermochemical energy storage systems are proven to be promising energy carriers (solar fuels) to utilize solar energy by using reactive solid-state pellets. However, the production of solar fuel requires a quasi-steady-state process temperature, which represents the main challenge due to the transient nature of solar power. In this work, an adaptive model predictive controller (MPC) is presented to regulate the temperature inside a tubular solar reactor to produce solid-state solar fuel for long-term thermal storage systems. The solar reactor system consists of a vertical tube heated circumferentially over a segment of its length by concentrated solar power, and the reactive pellets (MgMn2O4) are fed from the top end and flow downwards through the heated tube. A countercurrent flowing gas supplied from the lower end interacts with flowing pellets to reduce it thermochemically at a temperature range of 1000—1500 °C. A low-order physical model was developed to simulate the dynamics of the solar reactor including the reaction kinetics, and the proposed model was validated numerically by using a 7-kW electric furnace. The numerical model then was utilized to design the MPC controller, where the control system consists of an MPC code linked to an adaptive system identification code that updates system parameters online to ensure system robustness against external disturbances (sudden change in the flow inside the reactor), model mismatches, and uncertainty. The MPC controller parameters are tuned to enhance the system performance with minimum steady-state error and overshoot. The controller is tested to track different temperature ranges between 500 °C and 1400 °C with different particles/gas mass flowrates and ramping temperature profiles. Results show that the MPC controller successfully regulated the reactor temperature within ± 1 °C of its setpoint and maintained robust performance with minimum input effort when subjected to sudden changes in the amount of flowing media and the presence of chemical reaction.

Engineering↗
DOE OSTI2020

Testing general relativity on cosmological scales at redshift z ∼ 1.5 with quasar and CMB lensing

ABSTRACT We test general relativity (GR) at the effective redshift $\bar{z} \sim 1.5$ by estimating the statistic EG, a probe of gravity, on cosmological scales $19 - 190\, h^{-1}{\rm Mpc}$. This is the highest redshift and largest scale estimation of EG so far. We use the quasar sample with redshifts 0.8 < z < 2.2 from Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 16 as the large-scale structure (LSS) tracer, for which the angular power spectrum $C_\ell ^{qq}$ and the redshift-space distortion parameter β are estimated. By cross-correlating with the Planck 2018 cosmic microwave background (CMB) lensing map, we detect the angular cross-power spectrum $C_\ell ^{\kappa q}$ signal at $12\, \sigma$ significance. Both jackknife resampling and simulations are used to estimate the covariance matrix (CM) of EG at five bins covering different scales, with the later preferred for its better constraints on the covariances. We find EG estimates agree with the GR prediction at $1\, \sigma$ level over all these scales. With the CM estimated with 300 simulations, we report a best-fitting scale-averaged estimate of $E_G(\bar{z})=0.30\pm 0.05$, which is in line with the GR prediction $E_G^{\rm GR}(\bar{z})=0.33$ with Planck 2018 CMB + BAO matter density fraction Ωm = 0.31. The statistical errors of EG with future LSS surveys at similar redshifts will be reduced by an order of magnitude, which makes it possible to constrain modified gravity models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Sensitivity Analysis and Effective Parametrization of PEM Fuel Cell Models

The cost of proton-exchange-membrane fuel cells (PEMFCs) remains a major hurdle in large-scale commercialization of this technology. To improve their performance and reduce cost, novel materials and electrode designs are continuously envisioned, e.g., non-PGM catalyst layers, ultra-thin Pt/Pt-Ni based catalyst layers, structured ionomer arrays or NSTF catalyst layers.1 Understanding the impact of these improvement strategies can be extremely time and cost intensive due to complex physical phenomena and large design space. We have previously developed a PEMFC modeling framework2 which has been a time and cost effective tool for understanding and optimizing the complex multi-physics phenomena within PEMFCs; however, several of the cell parameters used in the modeling have large spread in measured data.3 Furthermore, several transport parameters such as water adsorption kinetics have not been accurately measured and the approximations are spread over several orders of magnitude. These uncertainties cause problems in ascertaining accuracy of the modeling approach and they reduce the predictive power of the numerical models. The aim of this work is to identify the sensitivity of PEFC numerical model outputs to various input parameters. The previously in-house developed MEA modeling framework2 is used for PEMFC modeling. The sensitivity of the model outputs with respect to inputs parameters is obtained by analyzing the condition numbers for different output-input pairs at varying operating conditions. An example of the sensitivity analysis is shown in Figure 1. The condition numbers are obtained for the entire possible range of input parameters at varying operating conditions to identify the most crucial parameters of the PEMFC model. Based on our preliminary analysis, parameters related to kinetics (exchange current density and ECSA) and heat/water management in electrodes and ionomer (ionomer fraction, thermal conductivity) are most crucial. One of the major goals of this work is to identify the most crucial set of parameters towards which the model shows maximum sensitivity. This will guide future experimentalists to measure these properties with higher accuracy. Furthermore, the sensitivity analysis will also enable us to optimize the PEMFC performance by selectively targeting the most sensitive parameters and thereby making the largest impact. Acknowledgements The work is funded under the Fuel Cell Performance and Durability Consortium (FC-PAD), by the Fuel Cell Technologies Office (FCTO), Office of Energy Efficiency and Renewable Energy (EERE), of the U.S. Department of Energy under contract number DE-AC02-05CH11231. The authors would like to thank Nathan Craig at Robert Bosch LLC for his valuable input in designing the sensitivity analysis. The authors would also like to thank Giovanna Bucci and Matthias Hanauer at Robert Bosch for their valuable inputs and discussion. References P. K. Sinha, W. Gu, A. Kongkanand and E. Thompson, J. Electrochem. Soc., 158, B831 (2011). L. M. Pant, M. R. Gerhardt, N. Macauley, R. Mukundan, R. L. Borup and A. Z. Weber, Electrochim. Acta, 326, 134963 (2019). R. Vetter and J. O. Schumacher, ArXiv181110091 Phys. (2018). Figure 1

Pant, Lalit↗
DOE OSTI2021

Feasibility of full-core pin resolved CFD simulations of small modular reactor with momentum sources

Complex flow structure interactions and heat transfer processes take place in nuclear reactor cores. Given the extreme pressure/temperature and radioactive conditions inside the core, numerical simulations offer an attractive and sometimes more feasible approach to study the related flow and heat transfer phenomena in addition to the experiments. Under the Exascale Computing Project, the full-core simulation of a small modular reactor (SMR) has been pursued coupling Computational Fluid Dynamics (CFD) and neutronics. A key aspect of the modeling of SMR fuel assemblies is the presence of spacer grids and the mixing promoted by mixing vanes or the equivalent. A reduced order methodology is adopted based on momentum sources to mimic the mixing of the vanes. The momentum sources have been carefully calibrated with detailed Large Eddy Simulations (LES) of spacer grids performed with Nek5000. Modeling the spacer grid and mixing vanes (SGMV) effect without body-fitted computational grid avoids the excessive costs in resolving the local geometric details, and thus supports the simulation to be scaled up to the full core. Besides the progress on momentum source modeling, this paper also features the first full-core pin resolved CFD simulation ever performed to the authors' knowledge. This represents a significant advancement in capability for the CFD of nuclear reactors, which will hopefully serve as an inspiration for further integrating high-fidelity numerical simulations in actual engineering designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Water structure and electric fields at the interface of oil droplets

Interfacial water exhibits rich and complex behaviour, playing an important part in chemistry, biology, geology and engineering. However, there is still much debate on the fundamental properties of water at hydrophobic interfaces, such as orientational ordering, the concentration of hydronium and hydroxide, improper hydrogen bonds and the presence of large electric fields. This controversy arises from the challenges in measuring interfacial systems, even with the most advanced experimental techniques and theoretical approaches available. Here we report on an in-solution, interface-selective Raman spectroscopy method using multivariate curve resolution to probe hexadecane-in-water emulsions, aided by a monomer-field theoretical model for Raman spectroscopy. Our results indicate that oil-water emulsion interfaces can exhibit reduced tetrahedral order and weaker hydrogen bonding, along with a substantial population of free hydroxyl groups that experience about 95 cm -1 redshift in their stretching mode compared with planar oil-water interfaces. Given the known electrostatic zeta potential characteristic of oil droplets, we propose the existence of a strong electric field (about 50-90 MV cm -1 ) emanating from the oil phase. This field is inferred indirectly but supported by control experiments and theoretical estimates. These observations are either absent or opposite in the molecular hydrophobic interface formed by small solutes or at planar oil-water interfaces. Instead, water structural disorder and enhanced electric fields emerge as unique features of the mesoscale interface in oil-water emulsions, potentially contributing to the accelerated chemical reactivity observed at hydrophobic-water interfaces.

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

Language models for the prediction of SARS-CoV-2 inhibitors

The COVID-19 pandemic highlights the need for computational tools to automate and accelerate drug design for novel protein targets. We leverage deep learning language models to generate and score drug candidates based on predicted protein binding affinity. We pre-trained a deep learning language model (BERT) on ∼9.6 billion molecules and achieved peak performance of 603 petaflops in mixed precision. Our work reduces pre-training time from days to hours, compared to previous efforts with this architecture, while also increasing the dataset size by nearly an order of magnitude. For scoring, we fine-tuned the language model using an assembled set of thousands of protein targets with binding affinity data and searched for inhibitors of specific protein targets, SARS-CoV-2 Mpro and PLpro. We utilized a genetic algorithm approach for finding optimal candidates using the generation and scoring capabilities of the language model. Our generalizable models accelerate the identification of inhibitors for emerging therapeutic targets.

Blanchard, Andrew E.↗
DOE OSTI2024

Paradigm for approaching the forbidden phase transition in the one-dimensional Ising model at fixed finite temperature: Single chain in a magnetic field

In a previous paper, the forbidden spontaneous phase transition in the one-dimensional Ising model was found to be approachable arbitrarily closely in decorated ladders by ultra-narrow phase crossover (UNPC) at a given finite temperature T 0 with the crossover width 2δT reduced exponentially, which resemble a genuine first-order transition with large latent heat. Here, I reveal that the forbidden phase transition can be approached at fixed T 0 as well in decorated single-chain Ising models in the presence of a magnetic field, in which T 0 is determined by the interactions involving only the decorated parts and the magnetic field, while 2δT is independently, exponentially reduced (δT = 0 means a genuine transition) by restoring the ferromagnetic interaction between the ordinary spins on the chain backbone-which was neglected in the previous studies of pseudo-transition-thus manifesting that this asymptoticity to the forbidden transition is essentially the buildup of coherence in preformed crossover of local states. Furthermore, I show that the UNPC can be realized even in the absence of the conventional geometric frustration because the magnetic field itself can induce previously unnoticed hidden spin frustration. In conclusion, these findings make the doors wide open to the engineering and utilization of UNPC as a new paradigm for exploring exotic phenomena and 1D device applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Spin wave Hamiltonian and anomalous scattering in NiPS 3

Here, we report a comprehensive spin wave analysis of the semiconducting honeycomb van der Waal antiferromagnet NiPS 3 . Using single-crystal inelastic neutron scattering, we map out the full Brillouin zone and fit the observed modes to a spin wave model with rigorously defined uncertainty. We find that the third-nearest-neighbor exchange J 3 dominates the Hamiltonian, a feature which we fully account for by ab initio density functional theory calculations. We also quantify the degree to which the threefold rotation symmetry is broken and account for the Q = 0 excitations observed in other measurements, yielding a spin exchange model which is consistent across multiple experimental probes. We also identify a strongly reduced static ordered moment and reduced low-energy intensity relative to the linear spin wave calculations, signaling unexplained features in the magnetism which requires going beyond the linear spin wave approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

EMI Mitigation of a &Cacute;uk-Based Power-Electronic System Using Switching-Sequence-Based Control

Switching-sequence-based control (SBC) laws when designed based on topological switching behavior can have positive effects on slow- and fast-scale dynamics of a power-electronic system (PES). The slow-scale control can encompass fast PES state regulation and tracking, based on predefined objective, while fastscale control can address differential-mode (DM) and commonmode (CM) spectral-peak energy associated with PES switching operation. Such control laws may offer enhanced programmability to conventional PES design where bulky electromagnetic interference (EMI) filters have been traditionally used to reduce EMI of switching power converters to meet EMI regulatory standards. An EMI filter is always a less programmable solution since it is usually designed for the worst-case EMI mitigation and usually overkill for a PES operating under reduced load condition. The control scheme outlined in this article offers EMI mitigation across wide operating regions without compromising PES regulation. Moreover, it does so by use of switching sequences that guarantee the reachability of the PES dynamics using an advanced Lyapunov-function-based approach. SBC is a powerful tool to generate control actions for a PES based on multivariate PES state constraints. Hence, contemporary EMI regulatory standards are used as constraints in the SBC formulation to operate the PES under wide operating regime while autonomously mitigating the EMI levels. The work may be of paramount importance for operating the ultra-fast-transition recent wide-bandgap semiconductor devices like GaN–FET and SiC MOSFET under higher power with increasing switching frequencies, which is usually desirable for increased power density and reduced switching losses. Here, a hardware Cuk–PES operated ´ with GaN–FETs is fabricated and is used for case illustration. It is shown by experimental results how SBC mitigate DM and CM EMI noise of the PES while maintaining regulation even for the higher order nonminimum phase PES, while reducing sensor requirements using state observer derived from the switching model of the PES.

42 ENGINEERING↗
DOE OSTI2020

Constraining Properties of the Next Nearby Core-collapse Supernova with Multimessenger Signals

With the advent of modern neutrino and gravitational wave (GW) detectors, the promise of multimessenger detections of the next galactic core-collapse supernova (CCSN) has become very real. Such detections will give insight into the CCSN mechanism and the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible CCSN events, and their multimessenger signals, is needed. We present detailed predictions of neutrino and GW signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from 9 to 120 M ⊙ . In order to achieve explosions in 1D, we use the Supernova Turbulence In Reduced-dimensionality model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the GW emission from the 1D simulations using an astroseismology analysis of the protoneutron star. We find that the neutrino and GW signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and GW emission from a galactic CCSN may be able to provide constraints on stellar evolution independent of preexplosion imaging and the mass of the compact object remnant prior to fallback accretion.

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↗