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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Scaling for Nuclear Reactor System: Overview

A scaling analysis is essential for predicting the behavior of a reactor system under various operating conditions, as well as for ensuring safe and reliable experimentation and adequate dataset generation for supporting the system code model development and assessment in accordance with USNRC guidelines 1.203. Testing and modeling can be carried out at various scales, from small test units to full-scale reactors, to validate the predictions produced by the analysis. However, due to the cost and challenges associated with full-scale testing, the nuclear industry uses scaled test facilities—integral effect test (IET) and separate effect test (SET)—to identify potential design issues and optimize the design to ensure regulatory compliance for safer operation and anticipated accident management. In general, IET facilities are non-nuclear thermal hydraulics test facilities. In contrast, SET facilities can be non-nuclear testing (e.g., steam generator, containment, pump system SETs) and nuclear testing (mostly fuel-material irradiation testing). IETs are used to study the system-level thermal-hydraulic behavior of the reactor coolant system (RCS) under different anticipated operational occurrences (AOOs) and design basis accidents (DBAs) such as loss-of-coolant-accident (LOCA) and loss-of-power (i.e., station blackout) and loss-of-flow (i.e., failure of RCS pump). Several LOCA scenarios could happen in reactor DBAs, such as main steam line break (MSLB), SG tube rapture (SGTR), break in direct vessel injection (DVI), and other small breaks LOCA. Scaling analysis determines the appropriate size and operating conditions for the test facility. This involves analyzing the system's thermal-hydraulic behavior at different scales and identifying the key parameters that affect the thermal-hydraulic behavior of the reactor. These key parameters can be identified so that the scaled facility can replicate the expected conditions in a prototypical setting. The scaled facility can simulate the same flow rates, pressure gradients, heat transfer coefficients, and fluid properties as the prototype. Scaling analysis is used to ensure that the results from these tests can be scaled up to predict the behavior of a full-scale prototype with an acceptable range of scaling distortions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Local atomic environment analysis of short and long-range solute-solute interactions in a symmetric tilt grain boundary

We develop an efficient molecular-statics algorithm for predicting local solute-solute interaction effects on grain-boundary segregation. This algorithm is applied to the segregation and interaction of Ni and Cu solutes in a Σ9 (221) symmetric tilt grain boundary in Ag. It is found that the strength of interaction between two solutes directly correlates with the change of average segregation energy in the local atomic environment of a fixed grain-boundary solute. It is shown that local interaction effects at grain boundaries can be either attractive or repulsive wherein long-range attractive effects increase the propensity for solute cluster nucleation. This study underscores the atomic nature of grain-boundary solute segregation and provides a tool that could yield meaningful insight into heterogeneous solute segregation phenomena in polycrystalline alloys beyond their dilute limit.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Imaging Wavepackets in Real-Space & Disentangling Ultrafast Solvation Dynamics with High-Energy Ultrafast X-ray Scattering

The microscopic information of solute and solute-solvent motions can be measured using ultrafast diffuse x-ray scattering, and compared to molecular dynamics simulations, as demonstrated in studies of solvated molecular systems done in X-ray free-electron lasers. However, for the typical photon energies used in these experiments, the photon momentum transfer range was limited to lower Q values where the signals arising from different parts of the studied system overlap significantly. Recently, this limitation was greatly alleviated by extending the photon energy at the Linac Coherent Light-Source (LCLS), significantly increasing the accessible momentum transfer range. This improvement is transformative for ultrafast diffuse x-ray scattering measurements, enabling for the first time to disentangle the solute scattering difference signal, that persists to higher Q values, from the bulk solvent and solute-solvent cross-terms difference signals. In addition, it enhances the prediction ability of modeling and simulation methods such as the hybrid QM/MM approach. The extended Q range also opens the way to resolve in real-space details regarding coherent wavepackets motions beyond their center-of-mass positions. Here, we present the first results on high-energy (18keV) diffuse x-ray scattering in solution, demonstrating high fidelity time-resolved scattering and analysis of the photoexcited model photocatalysts PtPOP ([Pt2(POP)4]4-), and IrDimen ([Ir2(dimen)4]2+) in several solvents. These complexes provide ideal systems for demonstrating the ability of high-Q ultrafast scattering and QM/MM simulations to decompose ultrafast solvation dynamics into specific changes in the solute-solvent pair distribution function with a particular focus on how electronically excited states change the interaction between the solvent and photo-catalytically active metal sites. We introduce an approach to further utilize the high-energy capability and demonstrate a single-shot ultra-wide-angle X-ray Scattering modality using two perpendicular detectors, spanning a scattering angle range of more than 100 degrees, allowing to extend the accessible momentum transfer range up to Q~14 Å-1. We develop a model-free method to invert the x-ray scattering signals and enable the recovery of multiple pair-density motions that happen simultaneously, allowing us to directly measure in real-space nuclear wavepacket motions. . [1] Natan, Adi. "Real-Space Inversion and Super-Resolution of Ultrafast X-ray Scattering using Natural Scattering Kernels." arXiv preprint arXiv:2107.05576 (2021)

Natan, Adi↗
DOE OSTI2021

Artificial neural network prediction of self-diffusion in pure compounds over multiple phase regimes

Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for pure components in liquid, gas and super critical phases. The ANNs were tested on an experimental database of 6625 self-diffusion constants for 118 different chemical compounds. The presence of multiple phases results in a heavy skew in the distribution of diffusion constants and multiple approaches were used to address this challenge. First, an ANN was developed with the raw diffusion values to assess what the main drawbacks of this direct method were. The first approach for improving the predictions involved taking the log 10 of diffusion to provide a more uniform distribution and reduce the range of target output values used to develop the ANN. The second approach involved developing individual ANNs for each phase using the raw diffusion values. Results show that the log transformation leads to a model with the best self-diffusion constant predictions and an overall average absolute deviation (AAD) of 6.56%. The resultant ANN is a generalized model that can be used to predict diffusion across all three phases and over a diverse group of compounds. The importance of each input feature was ranked using a feature addition method revealing that the density of the compound has the largest impact on the ANN prediction of self-diffusion constants in pure compounds.

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

A quantitative model for a nanoscale switch accurately predicts thermal actuation behavior

Manipulation of temperature can be used to actuate DNA origami nano-hinges containing gold nanoparticles. We develop a physical model of this system that uses partition function analysis of the interaction between the nano-hinge and nanoparticle to predict the probability that the nano-hinge is open at a given temperature. The model agrees well with experimental data and predicts experimental conditions that allow the actuation temperature of the nano-hinge to be tuned over a range of temperatures from 30 °C to 45 °C. Additionally, the model identifies microscopic interactions that are important to the macroscopic behavior of the system, revealing surprising features of the system. This combination of physical insight and predictive potential is likely to inform future designs that integrate nanoparticles into dynamic DNA origami structures or use strand binding interactions to control dynamic DNA origami behavior. Furthermore, our modeling approach could be expanded to consider the incorporation, stability, and actuation of other types of functional elements or actuation mechanisms integrated into nucleic acid devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Global climate mode resonance due to rapidly intensifying El Niño-Southern Oscillation

The El Niño-Southern Oscillation (ENSO) influences climate variability globally, encompassing various other modes of variability, and thus represents a key predictable climate signal on seasonal timescales. Yet, its response to greenhouse warming remains uncertain, with models projecting a range of outcomes. Here, we demonstrate that in response to warming, a state-of-the-art high-resolution climate model simulates a rapid transition from a moderate-amplitude irregular regime, as observed in the current climate, to a highly regular oscillation with intensifying amplitude. This behaviour can be attributed to increasing air-sea feedbacks, which approach criticality in the second half of this century, and growing atmospheric noise. As ENSO intensifies in this model, it synchronizes with other prominent climate modes, such as the North Atlantic Oscillation and the Indian Ocean Dipole, thereby imprinting its regular, predictable variability on them. If realized, this global climate mode resonance would have wide-ranging whiplash impacts on regional hydroclimates.

Atmospheric dynamics↗
DOE OSTI2024

Evaluation of the EarthSHAB Stratospheric Solar Hot Air Balloon Flight Prediction Model Using Balloon Trajectory Data

Abstract The heliotrope is a solar balloon design which is constructed out of painter’s plastic, and the exterior is coated in charcoal powder. Darkening the plastic gives the balloon a high solar absorptance, which allows it to ascend into the lower stratosphere and float for hours at a time. The balloons have previously been used to lift scientific instruments into the stratosphere to study chemical explosions, earthquakes, and stratospheric aerosols. They have also been proposed as a platform for planetary exploration. Flight predictions are crucial to preflight planning to reduce safety risks and meet flight objectives. However, there exists a wide range of possible flight paths due to varying environmental conditions and solar balloon configurations. EarthSHAB is one such software that was designed to support flight planning using the weather forecasts and balloon properties to predict the flight path of a solar balloon. We compare EarthSHAB-simulated flight paths to a set of observed flight paths for the 3.5-m diameter heliotrope design called the “Cloudskimmer.” Using the criteria that the modeled paths must fall within 5% of the observations to be considered successful, we found that EarthSHAB successfully predicted the Cloudskimmer ascent rate and average float altitude 10% and 90% of the time, respectively. We also found that the average difference in the observed and predicted landing locations was 97 km and landing times were 54 ± 38 min. Significant deviations between the observed and predicted ascent rates and excursions at float were found to be associated with heavy payloads and convective cloud development, respectively. Significance Statement Solar balloons are used to lift scientific instruments into the lower stratosphere for hours at a time to study chemical explosions, earthquakes, stratospheric aerosols, and more. The flight paths of solar balloons can be difficult to predict due to variability in their design and surrounding environment. We evaluate the accuracy of EarthSHAB, a software that predicts the altitude profile and horizontal trajectory of a solar balloon using inputs such as the weather, balloon size, and balloon mass. Our results suggest that for the balloon design used in this study, EarthSHAB is best suited for modeling the behavior of balloons with lightweight payloads that do not fly within or directly above clouds.

Lien, Jessica M. [Sandia National Laboratories, Al↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2023

Harnessing Uncertainty through Functional Data Analysis in Gas Breakthrough Data

Detecting subsurface explosions from radionuclide gas migration through rock fractures is an effective way to identify nuclear activity. Los Alamos National Laboratory (LANL) has developed simulation methods, based on data from the 1962 Hardhat underground nuclear test, to predict gas breakthrough times at the surface. However, these methods rely on an imperfect understanding of the relationship between rock damage and fracture permeability. Our clinic project studies methods for predicting breakthrough curves that characterize total mass produced as a function of time, as well as quantifying the uncertainty associated with these predictions. The model that is currently employed to relate damage to permeability uses an empirically motivated power-law expression, with a range of parameter values that are compatible with the experimental Hardhat data. We develop emulators, built from functional data analysis techniques and trained on simulation data, that rapidly predict the gas breakthrough curve given a damage field and given the parameter values of the power-law equation. Using Bayesian regression, we address the problem of uncertainty quantification in our emulators. Finally, in order to test the robustness of the model, we further validate it on a damage field representing different physical conditions.

58 GEOSCIENCES↗
DOE OSTI2023

Physics‐Informed Machine Learning for Inverse Design of Optical Metamaterials

Optical metamaterials manipulate light through various confinement and scattering processes, offering unique advantages like high performance, small form factor and easy integration with semiconductor devices. However, designing metasurfaces with suitable optical responses for complex metamaterial systems remains challenging due to the exponentially growing computation cost and the ill‐posed nature of inverse problems. To expedite the computation for the inverse design of metasurfaces, a physics‐informed deep learning (DL) framework is used. A tandem DL architecture with physics‐based learning is used to select designs that are scientifically consistent, have low error in design prediction, and accurate reconstruction of optical responses. The authors focus on the inverse design of a representative plasmonic device and consider the prediction of design for the optical response of a single wavelength incident or a spectrum of wavelength in the visible light range. The physics‐based constraint is derived from solving the electromagnetic wave equations for a simplified homogenized model. The model converges with an accuracy up to 97% for inverse design prediction with the optical response for the visible light spectrum as input, and up to 96% for optical response of single wavelength of light as input, with optical response reconstruction accuracy of 99%.

Sarkar, Sulagna↗
DOE OSTI2021

DeepGRN: prediction of transcription factor binding site across cell-types using attention-based deep neural networks

Abstract Background Due to the complexity of the biological systems, the prediction of the potential DNA binding sites for transcription factors remains a difficult problem in computational biology. Genomic DNA sequences and experimental results from parallel sequencing provide available information about the affinity and accessibility of genome and are commonly used features in binding sites prediction. The attention mechanism in deep learning has shown its capability to learn long-range dependencies from sequential data, such as sentences and voices. Until now, no study has applied this approach in binding site inference from massively parallel sequencing data. The successful applications of attention mechanism in similar input contexts motivate us to build and test new methods that can accurately determine the binding sites of transcription factors. Results In this study, we propose a novel tool (named DeepGRN) for transcription factors binding site prediction based on the combination of two components: single attention module and pairwise attention module. The performance of our methods is evaluated on the ENCODE-DREAM in vivo Transcription Factor Binding Site Prediction Challenge datasets. The results show that DeepGRN achieves higher unified scores in 6 of 13 targets than any of the top four methods in the DREAM challenge. We also demonstrate that the attention weights learned by the model are correlated with potential informative inputs, such as DNase-Seq coverage and motifs, which provide possible explanations for the predictive improvements in DeepGRN. Conclusions DeepGRN can automatically and effectively predict transcription factor binding sites from DNA sequences and DNase-Seq coverage. Furthermore, the visualization techniques we developed for the attention modules help to interpret how critical patterns from different types of input features are recognized by our model.

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

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Early season prediction of within-field crop yield variability by assimilating CubeSat data into a crop model

Accurate early season predictions of crop yield at the within-field scale can be used to address a range of crop production, management, and precision agricultural challenges. While the remote sensing of within-field insights has been a research goal for many years, it is only recently that observations with the required spatio-temporal resolutions, together with efficient assimilation methods to integrate these into modeling frameworks, have become available to advance yield prediction efforts. Here we explore a yield prediction approach that combines daily high-resolution CubeSat imagery with the APSIM crop model. The approach employs APSIM to train a linear regression that relates simulated yield to simulated leaf area index (LAI). That relationship is then used to identify the optimal regression date at which the LAI provides the best prediction of yield: in this case, approximately 14 weeks prior to harvest. Instead of applying the regression on satellite imagery that is coincident, or closest to, the regression date, our method implements a particle filter that integrates CubeSat-based LAI into APSIM to provide end-of-season high-resolution (3 m) yield maps weeks before the optimal regression date. The approach is demonstrated on a rainfed maize field located in Nebraska, USA, where suitable collections of both imagery and in-situ data were available for assessment. The procedure does not require in-field data to calibrate the regression model, with results showing that even with a single assimilation step, it is possible to provide yield estimates with good accuracy up to 21 days before the optimal regression date. Yield spatial variability was reproduced reasonably well, with a strong correlation to independently collected measurements (R 2 = 0.73 and rRMSE = 12%). When the field averaged yield was compared, our approach reduced yield prediction error from 1 Mg/ha (control case based on a calibrated APSIM model), to 0.5 Mg/ha (using satellite imagery alone), and then to 0.2 Mg/ha (results with assimilation up to three weeks prior to the optimal regression date). Such a capacity to provide spatially explicit yield predictions early in the season has considerable potential to enhance digital agricultural goals and improve end-of-season yield predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of High Fuel Loading and CO2 Dilution on Oxy-Methane Ignition Inside a Shock Tube at High Pressure

Abstract Ignition delay times were measured for methane/O2 mixtures in a high dilution environment of either CO2 or N2 using a shock tube facility. Experiments were performed between 1044 K and 1356 K at pressures near 16 ± 2 atm. Test mixtures had an equivalence ratio of 1.0 with 16.67% CH4, 33.33% O2, and 50% diluent. Ignition delay times were measured using OH* emission and pressure time-histories. Data were compared to the predictions of two literature kinetic mechanisms (ARAMCO MECH 2.0 and GRI Mech 3.0). Most experiments showed inhomogeneous (mild) ignition which was deduced from five time-of-arrival pressure transducers placed along the driven section of the shock tube. Further analysis included determination of blast wave velocities and locations away from the end wall of initial detonations. Blast velocities were 60–80% of CJ-Detonation calculations. A narrow high temperature region within the range was identified as showing homogenous (strong) ignition which showed generally good agreement with model predictions. Model comparisons with mild ignition cases should not be used to further refine kinetic mechanisms, though at these conditions, insight was gained into various ignition behavior. To the best of our knowledge, we present first shock tube data during ignition of high fuel loading CH4/O2 mixtures diluted with CO2 and N2.

Energy & Fuels↗
DOE OSTI2022

Predicting the Mechanical Response of Polyhydroxyalkanoate Biopolymers Using Molecular Dynamics Simulations

Polyhydroxyalkanoates (PHAs) have emerged as a promising class of biosynthesizable, biocompatible, and biodegradable polymers to replace petroleum-based plastics for addressing the global plastic pollution problem. Although PHAs offer a wide range of chemical diversity, the structure–property relationships in this class of polymers remain poorly established. In particular, the available experimental data on the mechanical properties is scarce. In this contribution, we have used molecular dynamics simulations employing a recently developed forcefield to predict chemical trends in mechanical properties of PHAs. Specifically, we make predictions for Young’s modulus, and yield stress for a wide range of PHAs that exhibit varying lengths of backbone and side chains as well as different side chain functional groups. Deformation simulations were performed at six different strain rates and six different temperatures to elucidate their influence on the mechanical properties. Our results indicate that Young’s modulus and yield stress decrease systematically with increase in the number of carbon atoms in the side chain as well as in the polymer backbone. In addition, we find that the mechanical properties were strongly correlated with the chemical nature of the functional group. The functional groups that enhance the interchain interactions lead to an enhancement in both the Young’s modulus and yield stress. Finally, we applied the developed methodology to study composition-dependence of the mechanical properties for a selected set of binary and ternary copolymers. Overall, our work not only provides insights into rational design rules for tailoring mechanical properties in PHAs, but also opens up avenues for future high throughput atomistic simulation studies geared towards identifying functional PHA polymer candidates for targeted applications.

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

Precipitation Probability and Its Future Changes From a Global Cloud-Resolving Model and CMIP6 Simulations

This study investigates precipitation probability and its future changes using simulations from the Non-hydrostatic Icosahedral Atmospheric Model (NICAM), a global cloud-resolving model, and Coupled Model Intercomparison Project Phase 6 (CMIP6) models. For the precipitation probability comparison between models and Tropical Rainfall Measuring Mission (TRMM) during 1998–2008, the difference is greatly reduced by regridding to a common coarser resolution, but the difference signs largely remain. Both NICAM and CMIP6 model ensemble mean show less global-mean nonprecipitation fractions than observations, but NICAM agrees better with observations. For global-mean light precipitation probability, NICAM (CMIP6 model ensemble mean) is smaller (larger) than TRMM. NICAM shows higher heavy precipitation probability than TRMM in native resolutions, but it is close to TRMM and falls in the CMIP6 model range after regridding. For the precipitation probability changes from 1979–2008 to 2075–2104, NICAM predicts an increase in nonprecipitation and extreme heavy (>100 mm/day) precipitation probability but a slight decrease in the probability for precipitation <100 mm/day. After regridding to 2.5°, the precipitation intensity where the probability starts to increase becomes ~20 mm/day and the relative increases in nonprecipitation and heavy precipitation probability become much larger. Future heavy precipitation probability relative increases from CMIP6 models are much larger than NICAM at their native resolutions, but NICAM falls in the CMIP6 model range after regridding, indicating robust future increase (decrease) in heavy (light) precipitation probability regardless of underlying convection treatments. Discrepancy, however, exists in the projected nonprecipitation fraction changes: NICAM shows that nonprecipitation probability will increase over most of low and middle latitudes where CMIP6 models show less increases and even decreases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uncertainty quantification for misspecified machine learned interatomic potentials

The use of high-dimensional regression techniques from machine learning has significantly improved the quantitative accuracy of interatomic potentials. Atomic simulations can now plausibly target quantitative predictions in a variety of settings, which has brought renewed interest in robust means to quantify uncertainties. In many practical settings where model complexity is constrained (e.g., due to performance considerations), misspecification — the inability of any one choice of model parameters to exactly match all training data — is a key contributor to errors that is often disregarded. Here, we employ a recent misspecification-aware regression technique to quantify parameter uncertainties, which is then propagated to a broad range of phase and defect properties in tungsten. The propagation is performed through both brute-force resampling and implicit Taylor expansion. The propagated misspecification uncertainties robustly quantify and bound errors on a broad range of material properties. We demonstrate application to recent foundational machine learning interatomic potentials, accurately predicting and bounding errors in MACE-MPA-0 energy predictions across the diverse materials project database.

36 MATERIALS SCIENCE↗