Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

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 883 records · Page 49

Core of 25 F in the rotational model

In a recent experiment, carried out at the Radioactive Isotope Beam Factory/RIKEN, the 25 F(p,2p) 24 O reaction was studied at 270 MeV/A in inverse kinematics. Derived spectroscopic factors suggest that the effective core of 25 F significantly differs from a free 24 O nucleus. We interpret these results within the particle-rotor model and show that the experimental level scheme of 25 F can be understood in the rotation-aligned coupling scheme with its 5/2$^{+}_{1}$ ground state as the bandhead of a decoupled band. Here, the excitation energies of the observed 1/2$^{+}_{1}$ and 9/2$^{+}_{1}$ states correlate strongly with the rotational energy of the effective core, seen by the odd proton, and allow us to estimate its 2 + energy at ≈3.2MeV and a moderate quadrupole deformation ϵ 2 ≈0.15. The measured fragmentation of the πd 5/2 single-particle strength is discussed, and some further experiments suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Geologic CO 2 Migration Pathways in Farnsworth Field, NW Anadarko Basin

This study reports on analyses of natural, geologic CO 2 migration paths in Farnsworth Oil Field, northern Texas, where CO 2 was injected into the Pennsylvanian Morrow B reservoir as part of enhanced oil recovery and carbon sequestration efforts. We interpret 2D and 3D seismic reflection datasets of the study site, which is located on the western flank of the Anadarko basin, and compare our seismic interpretations with results from a tracer study. Petroleum system models are developed to understand the petroleum system and petroleum- and CO 2 -migration pathways. We find no evidence of seismically resolvable faults in Farnsworth Field, but interpret a karst structure, erosional structures, and incised valleys. These interpretations are compared with results of a Morrow B well-to-well tracer study that suggests that inter-well flow is up-dip or lateral. Southeastward fluid flow is inhibited by dip direction, thinning, and draping of the Morrow B reservoir over a deeper, eroded formation. Petroleum system models predict a deep basin-ward increase in temperature and maturation of the source rocks. In the northwestern Anadarko Basin, petroleum migration was generally up-dip with local exceptions; the Morrow B sandstone was likely charged by formations both below and overlying the reservoir rock. Based on this analysis, we conclude that CO 2 escape in Farnsworth Field via geologic pathways such as tectonic faults is unlikely. Abandoned or aged wellbores remain a risk for CO 2 escape from the reservoir formation and deserve further monitoring and research.

58 GEOSCIENCES↗

Visual Analytics of Crosstalk in Quantum Hardware

Crosstalk remains a major obstacle to building scalable and fault-tolerant quantum computers. Conventional diagnostic techniques-often based on numerical simulation or statistical modeling-struggle to scale with hardware complexity and offer limited interpretability. In this work, we present a visual analytics framework for diagnosing qubit crosstalk using lightweight, circuit-based models integrated with an interactive user interface. Our approach quantifies correlations between active and idle qubits under parameterized single- and twoqubit operations, enabling detection of both spatial and gateinduced crosstalk. The system incorporates qubit topology and gate performance data to support sector-based exploration and correlation mapping. This tool assists users in identifying correlated error sources, informing qubit placement strategies, and guiding noise-aware circuit design.

Chae, Junghoon [ORNL] (ORCID:0000000206016746)↗

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗

Using AI Segmentation Models to Improve Foreign Body Detection and Triage from Ultrasound Images

Medical imaging can be a critical tool for triaging casualties in trauma situations. In remote or military medicine scenarios, triage is essential for identifying how to use limited resources or prioritize evacuation for the most serious cases. Ultrasound imaging, while portable and often available near the point of injury, can only be used for triage if images are properly acquired, interpreted, and objectively triage scored. Here, we detail how AI segmentation models can be used for improving image interpretation and objective triage evaluation for a medical application focused on foreign bodies embedded in tissues at variable distances from critical neurovascular features. Ultrasound images previously collected in a tissue phantom with or without neurovascular features were labeled with ground truth masks. These image sets were used to train two different segmentation AI frameworks: YOLOv7 and U-Net segmentation models. Overall, both approaches were successful in identifying shrapnel in the image set, with U-Net outperforming YOLOv7 for single-class segmentation. Both segmentation models were also evaluated with a more complex image set containing shrapnel, artery, vein, and nerve features. YOLOv7 obtained higher precision scores across multiple classes whereas U-Net achieved higher recall scores. Using each AI model, a triage distance metric was adapted to measure the proximity of shrapnel to the nearest neurovascular feature, with U-Net more closely mirroring the triage distances measured from ground truth labels. Overall, the segmentation AI models were successful in detecting shrapnel in ultrasound images and could allow for improved injury triage in emergency medicine scenarios.

47 OTHER INSTRUMENTATION↗

A robust estimator of mutual information for deep learning interpretability

Abstract We develop the use of mutual information (MI), a well-established metric in information theory, to interpret the inner workings of deep learning (DL) models. To accurately estimate MI from a finite number of samples, we present GMM-MI (pronounced ‘Jimmie’), an algorithm based on Gaussian mixture models that can be applied to both discrete and continuous settings. GMM-MI is computationally efficient, robust to the choice of hyperparameters and provides the uncertainty on the MI estimate due to the finite sample size. We extensively validate GMM-MI on toy data for which the ground truth MI is known, comparing its performance against established MI estimators. We then demonstrate the use of our MI estimator in the context of representation learning, working with synthetic data and physical datasets describing highly non-linear processes. We train DL models to encode high-dimensional data within a meaningful compressed (latent) representation, and use GMM-MI to quantify both the level of disentanglement between the latent variables, and their association with relevant physical quantities, thus unlocking the interpretability of the latent representation. We make GMM-MI publicly available in this GitHub repository.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Is the HEFT matching unique?

Physics beyond the Standard Model (BSM) can be described in a consistent and general way through the Higgs effective field theory (HEFT). Measurements of model-independent HEFT coefficients allow one to constrain the parameter space of BSM models via a matching procedure. In this work, we show that this procedure is not unique and depends on the scalings of the parameters of the Lagrangian. As examples, we consider three BSM models: the real singlet extension of the SM with a Z 2 symmetry, the complex singlet extension of the SM and the two Higgs doublet model. We discuss several physical observables, and show that different scalings of the model parameters with the UV scale in the matching to the HEFT can yield quite different results. This complicates the interpretation of HEFT measurements in terms of parameters of BSM models. Additionally, as a byproduct, we report the first matching of the complex singlet extension to the HEFT. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

Study of the agreement between binary neutron star ejecta models derived from numerical relativity simulations

Neutron star mergers have recently become a tool to study extreme gravity, nucleosynthesis, and the chemical composition of the Universe. To date, there has been one joint gravitational and electromagnetic observation of a binary neutron star merger, GW170817, as well as a solely gravitational observation, GW190425. In order to accurately identify and interpret electromagnetic signals of neutron star mergers, better models of the matter outflows generated by these mergers are required. We compare a series of ejecta models to see where they provide strong constraints on the amount of ejected mass expected from a system, and where systematic uncertainties in current models prevent us from reliably extracting information from observed events. We also examine 2396 neutron star equations of state compatible with GW170817 to see whether a given ejecta mass could be reasonably produced with a neutron star of said equation of state, and whether different ejecta models provide consistent predictions. We find that the difference between models is often comparable to or larger than the error generally assumed for these models, implying better constraints on the models are needed. In conclusion, we also note that the extrapolation of outflow models outside of their calibration window, while commonly needed to analyze gravitational wave events, is extremely unreliable and occasionally leads to completely unphysical results.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovery of Innovative Polymers for Next-Generation Gas-Separation Membranes using Interpretable Machine Learning

Polymer membranes perform innumerable separations with far-reaching environmental implications. Despite decades of research on membrane technologies, design of new membrane materials remains a largely Edisonian process. To address this shortcoming, we demonstrate a generalizable, accurate machine-learning (ML) implementation for the discovery of innovative polymers with ideal separation performance. Specifically, multitask ML models are trained on available experimental data to link polymer chemistry to gas permeabilities of He, H2, O2, N2, CO2, and CH4. Here, we interpret the ML models and extract chemical heuristics for membrane design, through Shapley Additive exPlanations (SHAP) analysis. We then screen over nine million hypothetical polymers through our models and identify thousands of candidates that lie well above current performance upper bounds. Notably, we discover hundreds of never-before-seen ultrapermeable polymer membranes with O2 and CO2 permeability greater than 104 and 105 Barrer, respectively. These hypothetical polymers are capable of overcoming undesirable trade-off relationship between permeability and selectivity, thus significantly expanding the currently limited library of polymer membranes for highly efficient gas separations. High-fidelity molecular dynamics simulations confirm the ML-predicted gas permeabilities of the promising candidates, which suggests that many can be translated to reality.

Yang, Jason↗

Standard Model in conformal geometry: Local vs gauged scale invariance

Here, we discuss comparatively local versus gauged Weyl symmetry beyond Standard Model (SM) and Einstein gravity and their geometric interpretation. The SM and Einstein gravity admit a natural embedding in Weyl integrable geometry which is a special limit of Weyl conformal (non-metric) geometry. The theory has a local Weyl scale symmetry but no associated gauge boson. Unlike previous models with such symmetry, this embedding is truly minimal i.e. with no additional fields beyond SM and underlying geometry. This theory is compared to a similar minimal embedding of SM and Einstein gravity in Weyl conformal geometry (SMW) which has a full gauged scale invariance, with an associated Weyl gauge boson. At large field values, both theories give realistic, Starobinsky–Higgs like inflation. The broken phase of the current model is the decoupling limit of the massive Weyl gauge boson of the broken phase of SMW, while the local scale symmetry of the current model is part of the larger gauged scale symmetry of SMW. Hence, the current theory has a gauge embedding in SMW. Unlike in the SMW, we note that in models with local scale symmetry the associated current is trivial, which is a concern for the physical meaning of this symmetry. Therefore, the SMW is a more fundamental UV completion of SM in a full gauge theory of scale invariance that generates Einstein gravity in the (spontaneously) broken phase, as an effective theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for charged Higgs bosons produced in vector boson fusion processes and decaying into vector boson pairs in proton–proton collisions at $\sqrt{s} = 13\,{\text {TeV}} $

A search for charged Higgs bosons produced in vector boson fusion processes and decaying into vector bosons, using proton–proton collisions at √s = 13 TeV at the LHC, is reported. The data sample corresponds to an integrated luminosity of 137 fb -1 collected with the CMS detector. Events are selected by requiring two or three electrons or muons, moderate missing transverse momentum, and two jets with a large rapidity separation and a large dijet mass. No excess of events with respect to the standard model background predictions is observed. Model independent upper limits at 95% confidence level are reported on the product of the cross section and branching fraction for vector boson fusion production of charged Higgs bosons as a function of mass, from 200 to 3000GeV. The results are interpreted in the context of the Georgi–Machacek model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unstructured clinical notes within the 24 hours since admission predict short, mid & long-term mortality in adult ICU patients

Mortality prediction for intensive care unit (ICU) patients is crucial for improving outcomes and efficient utilization of resources. Accessibility of electronic health records (EHR) has enabled data-driven predictive modeling using machine learning. However, very few studies rely solely on unstructured clinical notes from the EHR for mortality prediction. In this work, we propose a framework to predict short, mid, and long-term mortality in adult ICU patients using unstructured clinical notes from the MIMIC III database, natural language processing (NLP), and machine learning (ML) models. Depending on the statistical description of the patients’ length of stay, we define the short-term as 48-hour and 4-day period, the mid-term as 7-day and 10-day period, and the long-term as 15-day and 30-day period after admission. We found that by only using clinical notes within the 24 hours of admission, our framework can achieve a high area under the receiver operating characteristics (AU-ROC) score for short, mid and long-term mortality prediction tasks. The test AU-ROC scores are 0.87, 0.83, 0.83, 0.82, 0.82, and 0.82 for 48-hour, 4-day, 7-day, 10-day, 15-day, and 30-day period mortality prediction, respectively. We also provide a comparative study among three types of feature extraction techniques from NLP: frequency-based technique, fixed embedding-based technique, and dynamic embedding-based technique. Lastly, we provide an interpretation of the NLP-based predictive models using feature-importance scores.

60 APPLIED LIFE SCIENCES↗

Structure of 43 P and 42 Si in a two-level shape-coexistence model

Exclusive cross sections for the 43 P(-1p) 42 Si reaction to the lowest 0 + and 2 + states, measured at NSCL with GRETINA and the S800, are interpreted in terms of a two-level mixing (collective) model of oblate and prolate co-existing shapes. Here, using the formalism developed for deformed nuclei we calculate the spectroscopic amplitudes and exclusive cross-sections in the strong coupling limit, where for 43 P the schematic wavefunction includes the coupling of the Nilsson [211] 1/2 proton orbit. Good agreement with the experimental data is obtained when the amplitude of the oblate configuration is larger or equal than 80%, suggesting that both nuclei are predominantly oblate, in line with theoretical expectations. Model predictions for excitation energies and electromagnetic observables, that could motivate further studies of both nuclei, are also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An analysis of Bayesian estimates for missing higher orders in perturbative calculations

With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% credibility intervals in the Bayesian approach using the new public code MiHO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale mol­ecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

Integrated Operations for Nuclear Business Operation Model Analysis and Industry Validation

The purpose of this report is to refine and analyze five work reduction opportunities first presented in INL/EXT-21-64134, Process for Significant Nuclear Work Function Innovation Based on Integrated Operations Concepts. This report seeks to further refine and analyze five work reduction opportunities first presented in the original report. Researchers selected five work reduction opportunities from the full Integrated Operations for Nuclear (ION) suite. A selected group of utilities then verified details and inputs from the original report. Categories for verification included capital cost, technology requirements, and savings. Researchers then modeled the data points and data ranges using probabilistic analysis which predicts the likelihood of positive or negative net present value. Research results show four out of the five work reduction opportunities have a greater than fifty percent chance of a positive net present value outcome when analyzed independently. When the five work reduction opportunities are grouped and analyzed together the model indicates a sixty percent chance that the outcome of all five taken together will be positive. The nuclear industry should interpret these results as encouraging. In line with the ION model, positive financial analysis supports the investment of capital dollars into existing nuclear power plants along the ION model. Implementation of the five work reduction opportunities in this report is likely to result in substantive long-term savings for the owners and operators of domestic nuclear power plants.

99 GENERAL AND MISCELLANEOUS↗