Results for “Compact model” — Page 25

DOE OSTI2025

Weak Form Scientific Machine Learning: Test Function Construction for System Identification

Weak form Scientific Machine Learning (WSciML) is a recently developed framework for data-driven modeling and scientific discovery. It leverages the weak form of equation error residuals to provide enhanced noise robustness in system identification via convolving model equations with test functions, reformulating the problem to avoid direct differentiation of data. The performance, however, relies on wisely choosing a set of compactly supported test functions. In this work, we mathematically motivate a novel data-driven method for constructing Single-scale-Local reference functions for creating the set of test functions. Our approach numerically approximates the integration error introduced by the quadrature and identifies the support size for which the error is minimal, without requiring access to the model parameter values. Through numerical experiments across various models, noise levels, and temporal resolutions, we demonstrate that the selected supports consistently align with regions of minimal parameter estimation error. We also compare the proposed method against the strategy for constructing Multi-scale-Global (and orthogonal) test functions introduced in our prior work, demonstrating the improved computational efficiency.

FOS: Computer and information sciences↗

DOE OSTI2022

A Square Pulse Thermoreflectance Technique for the Measurement of Thermal Properties

We report on a laser-based square pulse thermoreflectance (SPTR) technique for the measurement of thermal properties for a wide range of materials. SPTR adopts the pump-probe thermoreflectance principle to monitor the evolution of local temperature after square pulse excitation. The technique features a compact setup, high spatial resolution, and fast data collection. By comparing the acquired SPTR signals with a continuum heat transfer model, material thermal properties can be obtained. Taking advantage of various spot sizes and modulation frequencies, SPTR can measure both the thermal diffusivity and thermal conductivity of poorly to moderately conductive materials and the thermal conductivity of conductive materials with satisfactory accuracy, with potential to be applied to more conductive materials. The technique was validated on three materials: fused silica, single crystal CaF2 and single crystal nickel (with conductivities ranging from 1 W·m -1 ·K -1 to 100 W·m -1 ·K -1 ) with typical measurement errors of 5 % to 20 %. The leading sources of error have been identified by Monte Carlo simulations, and the primary limitations of SPTR are discussed. The compact, fiberized platform we describe here will allow instruments based on this methodology to be deployed in complex, multi-analytical environments for the type of high-throughput correlative analyses that are key to materials design and discovery.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Elucidating the Role of Electric Fields in Fe Oxidation via an Environmental Atom Probe

Abstract We quantify the effects of intensely applied electric fields on the Fe oxidation mechanism. The specimen are pristine Fe single crystals exposing a variety of surface structures identified by field ion microscopy. These crystals are simultaneously exposed to low pressures of pure oxygen gas, on the order of 10 −7 mbar, while applying intense electric fields on their surface of several tens of volts per nanometer. The local composition of the different surface structures is probed directly and in real time using an Environmental Atom Probe and successfully compared with first principles‐based models. We found that rough Fe{244} and Fe{112} facets are more reactive toward oxygen than compact Fe{024} and Fe{011} facets. Results demonstrate that the influence of an electric field on the oxidation kinetics depends on the timescales that are involved as the system evolves toward equilibrium. The initial oxidation kinetics show that strong increases in electric fields facilitate the formation of an oxide. However, as one approaches equilibrium, high field values mitigate this formation. Ultimately, this study elucidates how high externally applied electric fields can be used to dynamically exploit reaction dynamics at the nanoscale towards desired products in a catalytic reaction at mild reaction conditions.

Lambeets, Sten V. [Physical and Computational Scie↗

DOE OSTI2025

Elucidating the Role of Electric Fields in Fe Oxidation via an Environmental Atom Probe

We quantify the effects of intensely applied electric fields on the Fe oxidation mechanism. The specimen are pristine Fe single crystals exposing a variety of surface structures identified by field ion microscopy. These crystals are simultaneously exposed to low pressures of pure oxygen gas, on the order of 10 −7 mbar, while applying intense electric fields on their surface of several tens of volts per nanometer. The local composition of the different surface structures is probed directly and in real time using an Environmental Atom Probe and successfully compared with first principles-based models. We found that rough Fe{244} and Fe{112} facets are more reactive toward oxygen than compact Fe{024} and Fe{011} facets. Results demonstrate that the influence of an electric field on the oxidation kinetics depends on the timescales that are involved as the system evolves toward equilibrium. The initial oxidation kinetics show that strong increases in electric fields facilitate the formation of an oxide. However, as one approaches equilibrium, high field values mitigate this formation. Ultimately, this study elucidates how high externally applied electric fields can be used to dynamically exploit reaction dynamics at the nanoscale towards desired products in a catalytic reaction at mild reaction conditions.

09 BIOMASS FUELS↗

DOE OSTI2021

Gravitational waves from a black hole orbiting in a wormhole geometry

Current ground-based gravitational wave detectors are tuned to observe compact object mergers from stellar mass black holes and neutron stars; ~50 such events have been published to date. More exotic compact objects may exist, collisions of which would also generate copious gravitational waves. In this work, we model a stellar mass black hole inspiral into a stable, nonspinning, traversable wormhole, and find a characteristic waveform—an antichirp and/or burst—as the black hole outspirals into our region of the Universe. This is a characteristic signature which can be useful in wormhole searches in gravitational wave data or used to constrain wormhole geometries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

DOE OSTI2020

AGR-5/6/7 Post-Irradiation Examination Plan

This plan describes the PIE activities to be carried out on AGR-5/6/7 fuel and the non-fuel components of the AGR-5/6/7 irradiation test train. PIE of AGR-5/6/7 fuel and test train components will provide data on fuel performance over a range of irradiation conditions (e.g., burnup, neutron fluence, and temperature), test fuel under postulated accident conditions, and support development of fuel performance and fission product transport models. Examination of the non-fuel test train components will give information on TRISO fuel and compact performance and information on releases of radioactive fission products from the fuel compacts. PIE of the fuel compacts is used for many purposes including, but not limited to: updating the irradiation thermal calculations, measuring burnup, determining retention/release of fission products from the compacts and the TRISO particles, determining TRISO layer failure rates, determining fuel performance at elevated accident temperatures, determining fuel performance under oxidizing atmospheres, determining the distribution of fission products within the compact and fuel particles, and observing compact and particle morphology to explain mechanisms of TRISO layer failure and/or fission product release/retention.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Statistically-informed deep learning for gravitational wave parameter estimation

We introduce deep learning models to estimate the masses of the binary components of black hole mergers, $(m_1,m_2)$, and three astrophysical properties of the post-merger compact remnant, namely, the final spin, $a_\mathrm f$, and the frequency and damping time of the ringdown oscillations of the fundamental $\ell = m = 2$ bar mode, $(\omega_\mathrm R, \omega_\mathrm I)$. Our neural networks combine a modified WaveNet architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters $(m_1,m_2, a_\mathrm f, \omega_\mathrm R, \omega_\mathrm I)$ of five binary black holes: GW150914, GW170104, GW170814, GW190521 and GW190630. We use PyCBC Inference to directly compare traditional Bayesian methodologies for parameter estimation with our deep learning based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90% confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 NVIDIA GPU to produce median values and posterior distributions within two milliseconds for each event. Furthermore, this neural network, and a tutorial for its use, are available at the Data and Learning Hub for Science.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Graph neural networks predict energetic and mechanical properties for models of solid solution metal alloy phases

Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The Muon Collider

Muons offer a unique opportunity to build a compact high-energy electroweak collider at the 10 TeV scale. A Muon Collider enables direct access to the underlying simplicity of the Standard Model and unparalleled reach beyond it. It will be a paradigm-shifting tool for particle physics representing the first collider to combine the high-energy reach of a proton collider and the high precision of an electron-positron collider, yielding a physics potential significantly greater than the sum of its individual parts. A high-energy muon collider is the natural next step in the exploration of fundamental physics after the HL-LHC and a natural complement to a future low-energy Higgs factory. Such a facility would significantly broaden the scope of particle colliders, engaging the many frontiers of the high energy community. The last European Strategy for Particle Physics Update and later the Particle Physics Project Prioritisation Panel in the US requested a study of the muon collider, which is being carried on by the International Muon Collider Collaboration. In this comprehensive document we present the physics case, the state of the work on accelerator design and technology, and propose an R\&D project that can make the muon collider a reality.

Accelerator Physics (physics.acc-ph)↗

DOE OSTI2022

Prediction of strain, inter-layer interaction and critical current in CORC ® wires under axial strain by T-A modeling

Superconducting conductors on round core (CORC ® ) cables and wires can meet the needs of large high-field magnets, such as particle accelerators and compact nuclear fusion machines, due to their simple cabling process, high current-carrying capacity and reliable operation under high mechanical stresses. Many high-field magnets require CORC ® cables to carry a current of thousands of amperes in a background magnetic field exceeding 20 T. As a result, the large electromagnetic forces will deform the cable in the axial direction due to hoop stress and in the transverse direction by compressive stress. Therefore, it is essential to determine the irreversible deformation limit of the CORC ® cable under axial tensile load and optimize the cabling parameters to potentially extend this limit. Analytical and numerical methods are developed to assess the performance degradation of CORC ® wires under axial tensile load. The strain level, interlayer contact pressure and friction and their impact on the critical current are calculated by combining the mechanical response and the T-A method. Analyzing the results shows that the winding angle of the tape and the Poisson’s ratio of the inner core are key factors affecting the irreversible tensile strain limit of CORC ® wires. The smaller the winding angle and the higher the Poisson’s ratio of the inner core, the higher the irreversible tensile strain limit. For multi-layer CORC ® wires, the initial contact pressure caused by the cabling process must also be considered. The inter-layer interaction is coupled with the tape strain of each layer. The results of this research can serve as a basis for optimizing and designing CORC ® wires with extended irreversible strain limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Testing and Characterization to Develop a Mechanistic Explanation for Unsaturated Drift of Fiber Optic Sensors during High-Dose Irradiation

The primary limitation for any optical fiber-based sensor for nuclear reactor applications is radiation-induced attenuation (RIA) of the transmitted and/or reflected signals. Based on several recent studies, RIA is tolerable for some fused silica optical fibers with the proper choice of sensing wavelength and fiber dopants. For extreme temperature applications (> 1000°C), sapphire optical fibers have been proposed; however, recent optical transmission measurements performed on bulk sapphire samples showed prohibitively large RIA. For some sensors, radiation-induced dimensional changes in the fiber materials can also cause significant drift. Moreover, the drift that was observed in numerous experiments performed in the High Flux Isotope Reactor (HFIR), the Advanced Test Reactor, the Massachusetts Institute of Technology Reactor, and other international facilities far exceeded what would be expected based on compaction of fused silica glass. Clearly, additional work is needed to better understand the origins of both RIA and radiation-induced drift in both silica and sapphire optical fiber-based sensors before these sensors can be reliably deployed for nuclear applications. This work evaluated the underlying mechanisms that may be responsible for RIA and drift in silica and sapphire materials. First, detailed characterization was performed on bulk fused silica glass samples that were previously irradiated to different neutron fluences at different temperatures to better understand the structural changes that drive radiation-induced drift in the absence of coating effects that are discussed later. Results show that the non-monotonic compaction that occurs with increasing neutron fluence continues up to fast neutron fluences approaching 10 22 n/cm 2 , which has important implications for physics-based models that may be used to compensate for the sensor drift. Initial Raman spectroscopy and synchrotron x-ray diffraction provide insights into the nature of the structural changes. Next, detailed characterizations were performed on silica fibers with various coatings that were subjected to several different thermal treatments. The hypothesis is that the coatings convert to carbon-rich materials that compact under irradiation, putting a large compressive strain on the fiber. Out-of-pile testing confirms that both polyimide and acrylate fiber coatings convert to glassy carbon (GC) materials when heated under inert conditions, and the degree of order (i.e., graphitization) increases with increasing temperature. The results provide increasingly strong evidence that the combination of polymeric coatings and inert (or vacuum) conditions render fiber optic sensors susceptible to significant radiation-induced drift that would not otherwise exist in uncoated fibers. Finally, transmission electron microscopy was performed on bulk sapphire samples that were irradiated to two neutron fluences at different temperatures to gain insights into the potential mechanisms driving the prohibitive RIA at higher neutron fluences and temperatures. Contrary to previous hypotheses, results show that scattering from radiation-induced voids cannot explain the observed RIA. Similarly, models for scattering losses from dislocation loops also do not agree with the experimental results. Instead, fitting to the experimental data shows that increased absorption from aluminum vacancy centers is the most likely explanation for the the prohibitively large RIA that was observed at high irradiation temperature and dose. In addition, the voids that formed in these single-crystal samples were found to align along the basal plane (a-axis) as opposed to that seen in previous observations of c-axis alignment in polycrystalline samples, which could have important implications for anisotropic swelling and other phenomena that could affect sensor performance at high neutron fluence.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Testing and Characterization to Develop a Mechanistic Explanation for Unsaturated Drift of Fiber Optic Sensors during High-Dose Irradiation

The primary limitation for any optical fiber-based sensor for nuclear reactor applications is radiation-induced attenuation (RIA) of the transmitted and/or reflected signals. Based on several recent studies, RIA is tolerable for some fused silica optical fibers with the proper choice of sensing wavelength and fiber dopants. For extreme temperature applications (> 1000°C), sapphire optical fibers have been proposed; however, recent optical transmission measurements performed on bulk sapphire samples showed prohibitively large RIA. For some sensors, radiation-induced dimensional changes in the fiber materials can also cause significant drift. Moreover, the drift that was observed in numerous experiments performed in the High Flux Isotope Reactor (HFIR), the Advanced Test Reactor, the Massachusetts Institute of Technology Reactor, and other international facilities far exceeded what would be expected based on compaction of fused silica glass. Clearly, additional work is needed to better understand the origins of both RIA and radiation-induced drift in both silica and sapphire optical fiber-based sensors before these sensors can be reliably deployed for nuclear applications. This work evaluated the underlying mechanisms that may be responsible for RIA and drift in silica and sapphire materials. First, detailed characterization was performed on bulk fused silica glass samples that were previously irradiated to different neutron fluences at different temperatures to better understand the structural changes that drive radiation-induced drift in the absence of coating effects that are discussed later. Results show that the non-monotonic compaction that occurs with increasing neutron fluence continues up to fast neutron fluences approaching 10 22 n/cm 2 , which has important implications for physics-based models that may be used to compensate for the sensor drift. Initial Raman spectroscopy and synchrotron x-ray diffraction provide insights into the nature of the structural changes. Next, detailed characterizations were performed on silica fibers with various coatings that were subjected to several different thermal treatments. The hypothesis is that the coatings convert to carbon-rich materials that compact under irradiation, putting a large compressive strain on the fiber. Out-of-pile testing confirms that both polyimide and acrylate fiber coatings convert to glassy carbon (GC) materials when heated under inert conditions, and the degree of order (i.e., graphitization) increases with increasing temperature. The results provide increasingly strong evidence that the combination of polymeric coatings and inert (or vacuum) conditions render fiber optic sensors susceptible to significant radiation-induced drift that would not otherwise exist in uncoated fibers. Finally, transmission electron microscopy was performed on bulk sapphire samples that were irradiated to two neutron fluences at different temperatures to gain insights into the potential mechanisms driving the prohibitive RIA at higher neutron fluences and temperatures. Contrary to previous hypotheses, results show that scattering from radiation-induced voids cannot explain the observed RIA. Similarly, models for scattering losses from dislocation loops also do not agree with the experimental results. Instead, fitting to the experimental data shows that increased absorption from aluminum vacancy centers is the most likely explanation for the the prohibitively large RIA that was observed at high irradiation temperature and dose. In addition, the voids that formed in these single-crystal samples were found to align along the basal plane (a-axis) as opposed to that seen in previous observations of c-axis alignment in polycrystalline samples, which could have important implications for anisotropic swelling and other phenomena that could affect sensor performance at high neutron fluence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Quantum simulation of bosons with the contracted quantum eigensolver

Abstract Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where N quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.

Physics↗

DOE OSTI2024

Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation

Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.

58 GEOSCIENCES↗

DOE OSTI2026

Modeling Multi-View Impedance-Based Cross-Geometry SOH Estimator for Li-ion Batteries

Abstract: Accurately estimating battery’s State of Health (SOH) remains challenging when models must generalize across cell designs and operating conditions. Most Electrochemical Impedance Spectroscopy (EIS)-based approaches either (i) hand-engineer a few Nyquist-plot features for shallow models—fast but does not generalize across geometries—or (ii) learn directly from Nyquist plots with deep networks, which removes manual feature extraction, yet still limited to a single plot type. As a result, cross-geometry robustness and deployability on constrained Internet of Things (IoT) devices remain open problems. We propose a compact Convolutional Neural Network (CNN) (∼ 10k parameters) that takes multi-representation EIS inputs—Nyquist (real/imaginary) and phase–magnitude (|Z|/ϕ) stacked as four channels, so the model can learn complementary degradation signatures while remaining small enough for fast inference. We build a dataset from cyclic aging of two geometries (LG INR18650MJ1 cylindrical cells and LIR2032 coin cells), acquire EIS every ten cycles from 10 kHz to 10 mHz (10 points/decade), and evaluate with leave-one-cell-out testing strategy. We further study fusion vs. single-representation inputs and assess feasibility for on-device deployment (e.g., NVIDIA Jetson device). The results show that training on multiple EIS representations improves SOH estimation accuracy and cross-geometry generalization compared to single-representation models, which uses only Nyquist or phase–magnitude plots. This design targets accurate, generalizable SOH prediction without manual feature engineering while enabling practical real-time use.

Bakr, Ahmed [The University of Alabama (UA)]↗

DOE OSTI2021

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗