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 253 records · Page 14

TCR-H: explainable machine learning prediction of T-cell receptor epitope binding on unseen datasets

Artificial-intelligence and machine-learning (AI/ML) approaches to predicting T-cell receptor (TCR)-epitope specificity achieve high performance metrics on test datasets which include sequences that are also part of the training set but fail to generalize to test sets consisting of epitopes and TCRs that are absent from the training set, i.e., are ‘unseen’ during training of the ML model. We present TCR-H, a supervised classification Support Vector Machines model using physicochemical features trained on the largest dataset available to date using only experimentally validated non-binders as negative datapoints. TCR-H exhibits an area under the curve of the receiver-operator characteristic (AUC of ROC) of 0.87 for epitope ‘hard splitting’ (i.e., on test sets with all epitopes unseen during ML training), 0.92 for TCR hard splitting and 0.89 for ‘strict splitting’ in which neither the epitopes nor the TCRs in the test set are seen in the training data. Furthermore, we employ the SHAP (Shapley additive explanations) eXplainable AI (XAI) method for post hoc interrogation to interpret the models trained with different hard splits, shedding light on the key physiochemical features driving model predictions. TCR-H thus represents a significant step towards general applicability and explainability of epitope:TCR specificity prediction.

60 APPLIED LIFE SCIENCES↗

Quantification of Cas9 binding and cleavage across diverse guide sequences maps landscapes of target engagement

The RNA-guided nuclease Cas9 has unlocked powerful methods for perturbing both the genome through targeted DNA cleavage and the regulome through targeted DNA binding, but limited biochemical data have hampered efforts to quantitatively model sequence perturbation of target binding and cleavage across diverse guide sequences. We present scalable, sequencing-based platforms for high-throughput filter binding and cleavage and then perform 62,444 quantitative binding and cleavage assays on 35,047 on- and off-target DNA sequences across 90 Cas9 ribonucleoproteins (RNPs) loaded with distinct guide RNAs. We observe that binding and cleavage efficacy, as well as specificity, vary substantially across RNPs; canonically studied guides often have atypically high specificity; sequence context surrounding the target modulates Cas9 on-rate; and Cas9 RNPs may sequester targets in nonproductive states that contribute to “proofreading” capability. Lastly, we distill our findings into an interpretable biophysical model that predicts changes in binding and cleavage for diverse target sequence perturbations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementation of a 9-point stencil in SOLPS-ITER and implications for Alcator C-Mod divertor plasma simulations

The SOLPS-ITER code suite is used worldwide for plasma edge modeling, the interpretation of experiments, as well as for the design of the ITER divertor. The numerical scheme of the plasma solver of the code, B2.5, is based on the assumption of perfectly field-aligned grids, while in practice grids are often strongly distorted to match divertor target shapes. Neglecting these grid distortion leads to qualitatively and quantitatively incorrect results for fluid neutral simulations, and may affect results in cold (detached) divertors even when using kinetic neutral simulations. In this contribution, we present the first results of a newly implemented 9-point stencil in B2.5 to properly handle misaligned grids. The new scheme is then applied to fluid neutral simulations of a well-diagnosed and previously modeled Alcator C-Mod discharge. Results are compared with the original 5-point scheme neglecting grid distortion effects, as well as with simulations including a full kinetic neutral model. We conclude that the 9-point stencil is essential to correctly model the transport of fluid neutrals on distorted grids, and to capture the effects of divertor closure on the fluid neutral behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrochemical Model Assessment of Strategies for Fast Charge

The United States Department of Energy (DOE) has identified extreme fast charging (XFC) as a critical challenge that must be overcome in order to achieve widespread adoption of electric vehicles. DOE’s eXtreme Fast Charge and Cell Evaluation (XCEL) program seeks to reduce charge time to 10-15 minutes for a 200+ Wh/kg Li-ion battery. This presentation for DOE’s Vehicle Technologies Program Annual Merit Review summarizes the XCEL team’s work developing electrochemical models to interpret data and applying them to identify best strategies to achieve fast charge. Strategies include advanced electrode architectures, elevated temperature charging, novel charge protocols, advanced electrolytes, thermal management and suppression of heterogeneities that cause early onset of the damaging Li plating side reaction.

47 OTHER INSTRUMENTATION↗

Applications of explainable artificial intelligence in renewable energy research

Researchers in renewable energy are applying deep learning (DL) to a variety of problems from diverse renewable energy domains, such as biofuels, wind, solar, power systems, buildings, vehicles, and transportation systems. Improvements in accuracy may be demonstrated using DL in laboratory settings. However, the lack of interpretability of DL models poses a practical limitation to their utility in advancing scientific knowledge and in the deployment of DL models in safety-critical energy systems. In this article, we discuss explainable artificial intelligence (XAI) as one pathway toward more interpretable DL models. We explore a brief timeline of U.S. national laboratory interest in XAI, an overview and taxonomy of methods in the field of XAI, and a selection of applications across renewable energy research domains. We conclude by highlighting pivotal areas where XAI can accelerate innovation in artificial intelligence for renewable energy research and other essential future directions.

97 MATHEMATICS AND COMPUTING↗

Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces

Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used to generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, including many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process using exact (as opposed to variational) analysis. We relate the theory to the existing continuous-state Gaussian diffusion as well as other approaches to discrete diffusion, and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce “Blackout Diffusion”, which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10, Binarized MNIST, and CelebA datasets confirm the feasibility of our approach. Generalizing from specific (Gaussian) forward processes to discrete-state processes without a variational approximation sheds light on how to interpret diffusion models, which we discuss.

Santos, Javier E.↗

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗

Projecting Large Fires in the Western US With an Interpretable and Accurate Hybrid Machine Learning Method

More frequent and widespread large fires are occurring in the western United States (US), yet reliable methods for predicting these fires, particularly with extended lead times and a high spatial resolution, remain challenging. In this study, we proposed an interpretable and accurate hybrid machine learning (ML) model, that explicitly represented the controls of fuel flammability, fuel availability, and human suppression effects on fires. The model demonstrated notable accuracy with a F 1 -score of 0.846 ± 0.012, surpassing process-driven fire danger indices and four commonly used ML models by up to 40% and 9%, respectively. More importantly, the ML model showed remarkably higher interpretability relative to other ML models. Specifically, by demystifying the “black box” of each ML model using the explainable AI techniques, we identified substantial structural differences across ML fire models, even among those with similar accuracy. The relationships between fires and their drivers, identified by our model, were aligned closer with established fire physical principles. The ML structural discrepancy led to diverse fire predictions and our model predictions exhibited greater consistency with actual fire occurrence. With the highly interpretable and accurate model, we revealed the strong compound effects from multiple climate variables related to evaporative demand, energy release component, temperature, and wind speed, on the dynamics of large fires and megafires in the western US. Our findings highlight the importance of assessing the structural integrity of models in addition to their accuracy. They also underscore the critical need to address the rise in compound climate extremes linked to large wildfires.

54 ENVIRONMENTAL SCIENCES↗

Versatile Sequential Casting Processing for Highly Efficient and Stable Binary Organic Photovoltaics

Forming an ideal bulk heterojunction (BHJ) morphology is a critical issue governing the photon to electron process in organic solar cells (OSCs). Complementary to the widely-used blend casting (BC) method for BHJ construction, sequential casting (SC) can also enable similar or even better morphology and device performance for OSCs. Here, BC and SC methods on three representative donor:acceptor (D:A) blends are utilized, that is, PM6:PC 71 BM, PM6:IT-4F and PM6:L8-BO. Higher power conversion efficiencies (PCEs) in all cases by taking advantage of beneficial morphology from SC processing are achieved, and a champion PCE of 18.86% (certified as 18.44%) based on the PM6:L8-BO blend is reached, representing the record value among binary OSCs. The observations on phase separation and vertical distribution inspire the proposal of the swelling–intercalation phase-separation model to interpret the morphology evolution during SC processing. Further, the vertical phase segregation is found to deliver an improvement of device performance via affecting the charge transport and collection processes, as evidenced by the D:A-ratio-dependent photovoltaic properties. Besides, OSCs based on SC processing show advantages on device photostability and upscale fabrication. Finally, this work demonstrates the versatility and efficacy of the SC method for BHJ-based OSCs.

14 SOLAR ENERGY↗

Search for higgsinos decaying to two Higgs bosons and missing transverse momentum in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

Results are presented from a search for physics beyond the standard model in proton-proton collisions at $ \sqrt{s} $ = 13 TeV in channels with two Higgs bosons, each decaying via the process H → b$ \overline{\mathrm{b}} $, and large missing transverse momentum. The search uses a data sample corresponding to an integrated luminosity of 137 fb$^{−1}$ collected by the CMS experiment at the CERN LHC. The search is motivated by models of supersymmetry that predict the production of neutralinos, the neutral partners of the electroweak gauge and Higgs bosons. The observed event yields in the signal regions are found to be consistent with the standard model background expectations. The results are interpreted using simplified models of supersymmetry. For the electroweak production of nearly mass-degenerate higgsinos, each of whose decay chains yields a neutralino $ \left({\overset{\sim }{\upchi}}_1^0\right) $ that in turn decays to a massless goldstino and a Higgs boson, $ \left({\overset{\sim }{\upchi}}_1^0\right) $ masses in the range 175 to 1025 GeV are excluded at 95% confidence level. For the strong production of gluino pairs decaying via a slightly lighter $ \left({\overset{\sim }{\upchi}}_2^0\right) $ to H and a light $ \left({\overset{\sim }{\upchi}}_1^0\right) $, gluino masses below 2330 GeV are excluded.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for the production of a Higgs boson in association with a single top quark in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search for the production of a Higgs boson in association with a single top quark, tH, is presented. The analysis uses proton-proton collision data corresponding to an integrated luminosity of 140 fb −1 at a centre-of-mass energy of 13 TeV, collected by the ATLAS detector at the LHC. The search targets Higgs-boson decays into $b\bar{b}$, WW * , ZZ * , and ττ, accompanied by an isolated lepton (electron or muon) from the top-quark decay. Multivariate techniques are employed to enhance the separation between signal and background processes. The observed signal strength, μ tH , defined as the ratio between the measured cross-section and the predicted Standard Model value, is μ tH = 8.1 ± 2.6 (stat.) ± 2.0 (syst.). The significance of the observed (expected) signal above the background-only expectation is 2.8 (0.4) standard deviations. The corresponding observed (expected) upper limit at the 95% confidence level on the tH cross-section is found to be 13.9 (6.1) times the value predicted by the Standard Model. An interpretation with an inverted sign of the top-quark Yukawa coupling is performed, and the signal strength and corresponding limit are reported.

Hadron-Hadron Scattering↗

Machine learning insights into microstructural origins of transport and mechanical properties in porous microstructures

Multifunctional porous materials are increasingly needed across various fields, but their complex microstructures create significant challenges due to the intricate microstructure-property relationships. This complexity, combined with limitations of traditional analysis methods, hinders efforts to understand and optimize microstructure–property relationships. Here, to address this, we integrate physics-based mesoscale modeling with interpretable machine learning (ML) to uncover how microstructural features govern effective diffusivity and elastic modulus. At constant porosity, we show diffusivity varies by over 150 × and modulus by ∼50 ×, highlighting the power of microstructure engineering. Statistical analysis reveals bimodal behavior in diffusivity and unimodal in modulus. ML identifies connectivity as the dominant factor, while modulus is also sensitive to domain size and feature interactions. Controlled simulations further highlight domain shape as a critical feature for modulus. This framework enables efficient exploration of microstructure-property correlations, offering new insights to guide the design of advanced porous materials.

Bicontinuous microstructure↗

Dating fractures using luminescence

Rock fracturing (cracking) is a universal process that drives and limits chemical degradation, sediment production and erosion, and deterioration of infrastructure. Despite extensive research gains in rock mechanics on one hand and geochronology on the other, there remains a glaring gap in our ability to understand the long term evolution of natural, in situ fractures. Here we develop a novel fracture exposure dating technique, grounded in modern advances in luminescence geochronology. We apply our new dating method to a granitic boulder from a glacial outwash terrace in California, US. We conclude that the longest, clast-splitting E-W fracture appeared shortly after the boulder's deposit, whereas the secondary N-S fracture appeared 5 ka after the deposition, approximately correlating with the Last Glacial Maximum and Younger Dryas periods of the region, respectively. However, dating of the third fracture (<< 50 µm width) which does not fully split the rock, is ambiguous due to negligible daylight penetration and poor determination of fracture width. The fracture dating method presented herein brings with it the potential to decipher relationships that are crucial for the interpretation and modeling of, for example, long-term landscape and atmospheric evolution relating rock weathering to climate change and erosion.

58 GEOSCIENCES↗

Physics-guided logistic classification for tool life modeling and process parameter optimization in machining

This paper describes a physics-guided logistic classification method for tool life modeling and process parameter optimization in machining. Tool life is modeled using a classification method since the exact tool life cannot be measured in a typical production environment where tool wear can only be directly measured when the tool is replaced. Here, in this study, laboratory tool wear experiments are used to simulate tool wear data normally collected during part production. Two states are defined: tool not worn (class 0) and tool worn (class 1). The non-linear reduction in tool life with cutting speed is modeled by applying a logarithmic transformation to the inputs for the logistic classification model. A method for interpretability of the logistic model coefficients is provided by comparison with the empirical Taylor tool life model. The method is validated using tool wear experiments for milling. Results show that the physics-guided logistic classification method can predict tool life using limited datasets. A method for pre-process optimization of machining parameters using a probabilistic machining cost model is presented. The proposed method offers a robust and practical approach to tool life modeling and process parameter optimization in a production environment.

Machine learning↗

Machine learning-guided design of direct methanol fuel cells with a platinum group metal-free cathode

Direct methanol fuel cells (DMFCs) offer a promising solution for clean electricity generation, particularly in small electronics and remote auxiliary power units. However, optimizing their efficiency and performance is challenging due to the complex interactions between various factors. Here, we present a novel approach that integrates experiments with machine learning to model and predict the performance of these fuel cells using atomically dispersed platinum group metal (PGM)-free catalysts at the cathode. Further, our machine learning models, trained on diverse input parameters, allow for the comprehensive optimization of DMFC performance prior to fabrication and testing. Through extensive experimental validation, we demonstrate that this data-driven approach accurately predicts key performance metrics, such as maximum power output and polarization curves. By combining our models with interpretable game-theory methods, we provide deep insights into the factors governing fuel cell performance, ultimately paving the way for the design of scalable and efficient DMFC technologies.

25 ENERGY STORAGE↗

The exit-wave power-cepstrum transform for scanning nanobeam electron diffraction: robust strain mapping at subnanometer resolution and subpicometer precision

Scanning nanobeam electron diffraction (NBED) with fast pixelated detectors is a valuable technique for rapid, spatially resolved mapping of lattice structure over a wide range of length scales. However, intensity variations caused by dynamical diffraction and sample mistilts can hinder the measurement of diffracted disk centers as necessary for quantification. Robust data processing techniques are needed to provide accurate and precise measurements for complex samples and non-ideal conditions. Here we present an approach to address these challenges using a transform, called the exit wave power cepstrum (EWPC), inspired by cepstral analysis in audio signal processing. The EWPC transforms NBED patterns into real-space patterns with sharp peaks corresponding to inter-atomic spacings. We describe a simple analytical model for interpretation of these patterns that cleanly decouples lattice information from the intensity variations in NBED patterns caused by tilt and thickness. With tracking the inter-atomic spacing peaks in EWPC patterns, strain mapping is demonstrated for two practical applications: mapping of ferroelectric domains in epitaxially strained PbTiO 3 films and mapping of strain profiles in arbitrarily oriented core-shell Pt-Co nanoparticle fuel-cell catalysts. The EWPC transform enables lattice structure measurement at sub-pm precision and sub-nm resolution that is robust to small sample mistilts and random orientations.

42 ENGINEERING↗

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗