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 1,027 records · Page 57

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune “bond” lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol’s role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA’s conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Altogether, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Roles of External Forcing and Internal Variability in Precipitation Changes of a Sub-Region of the U.S. Mid-Atlantic During 1979–2019

Observations of the seasonal mean precipitation averaged over a sub-region of the U.S. mid-Atlantic show significant increases in summer but no changes in winter during 1979–2019. Clarifying the drivers of these observed changes is important for understanding and interpreting future projections of regional precipitation by climate models. Herein, this study examines the relative effects of external forcing and internal variability in the observed multi-decadal changes using the historical and RCP8.5/SSP585 simulations from 37 models in the Coupled Model Intercomparison Project phase 5 (CMIP5), 35 models in CMIP6, as well as large ensemble simulations from 6 models, each including more than 20 members. During 1979–2019, enhanced warming over land relative to ocean due to external forcing induces low pressure anomalies over land that contribute to the wetting summer. While similar wetting also occurs during winter, such effect is offset by the drying associated with the positive phase of the North Atlantic Oscillation (NAO) in recent decades. Internal variability produces large uncertainty in the multi-decadal precipitation changes simulated by climate models through changes in large-scale circulation associated with the NAO in winter and the North Atlantic Subtropical High in summer. These multi-decadal relationships between precipitation and large-scale circulation under internal variability also manifest on inter-annual timescale. These relationships highlight enhanced moisture transport from the Gulf of Mexico and North Atlantic in winter and near the mid-Atlantic coast in summer that modulates precipitation change in the study region.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear ionization dynamics of hot dense plasma observed in a laser-plasma amplifier

Understanding the behaviour of matter under conditions of extreme temperature, pressure, density and electromagnetic fields has profound effects on our understanding of cosmologic objects and the formation of the universe. Lacking direct access to such objects, our interpretation of observed data mainly relies on theoretical models. However, such models, which need to encompass nuclear physics, atomic physics and plasma physics over a huge dynamic range in the dimensions of energy and time, can only provide reliable information if we can benchmark them to experiments under well-defined laboratory conditions. Due to the plethora of effects occurring in this kind of highly excited matter, characterizing isolated dynamics or obtaining direct insight remains challenging. High-density plasmas are turbulent and opaque for radiation below the plasma frequency and allow only near-surface insight into ionization processes with visible wavelengths. Here, the output of a high-harmonic seeded laser-plasma amplifier using eight-fold ionized krypton as the gain medium operating at a 32.8 nm wavelength is ptychographically imaged. A complex-valued wavefront is observed in the extreme ultraviolet (XUV) beam with high resolution. Ab initio spatio-temporal Maxwell–Bloch simulations show excellent agreement with the experimental observations, revealing overionization of krypton in the plasma channel due to nonlinear laser-plasma interactions, successfully validating this four-dimensional multiscale model. This constitutes the first experimental observation of the laser ion abundance reshaping a laser-plasma amplifier. The presented approach shows the possibility of directly modelling light-plasma interactions in extreme conditions, such as those present during the early times of the universe, with direct experimental verification.

36 MATERIALS SCIENCE↗

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗

Ab initio thermodynamics reveals the nanocomposite structure of ferrihydrite

Ferrihydrite is a poorly crystalline iron oxyhydroxide nanomineral that serves a critical role as the most bioavailable form of ferric iron for living systems. However, its atomic structure and composition remain unclear due in part to ambiguities in interpretation of X-ray scattering results. Prevailing models so far have not considered the prospect that at the level of individual nanoparticles multiple X-ray indistinguishable phases could coexist. Using ab initio thermodynamics we show that ferrihydrite is likely a nanocomposite of distinct structure types whose distribution depends on particle size, temperature, and hydration. Nanoparticles of two contrasting single-phase ferrihydrite models of Michel and Manceau are here shown to be thermodynamically equivalent across a wide range of temperature and pressure conditions despite differences in their structural water content. Higher temperature and water pressure favor the formation of the former, while lower temperature and water pressure favor the latter. For aqueous suspensions at ambient conditions, their coexistence is maximal for particle sizes up to 12 nm. The predictions inform and help resolve different observations in various experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting 1 H NMR relaxation in Gd 3+ -aqua using molecular dynamics simulations

Atomistic molecular dynamics simulations are used to predict 1 H NMR T 1 relaxation of water from paramagnetic Gd 3+ ions in solution at 25 °C. Simulations of the T 1 relaxivity dispersion function r 1 computed from the Gd 3+ – 1 H dipole–dipole autocorrelation function agree within ≃8% of measurements in the range f 0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. Here the simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f 0 ≲ 5 MHz, the simulation overestimates r 1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r 1 at high frequencies in chelated Gd 3+ contrast-agents used for clinical MRI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The cosmic web of X-ray active galactic nuclei seen through the eROSITA Final Equatorial Depth Survey (eFEDS)

Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA’s contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited (FX > 6.5 × 10 -15 erg s -1 cm -2 ) AGN sample at low redshift, 0.05 < z < 0.55. Two summary statistics, the clustering using spectra from SDSS-V and galaxy-galaxy lensing with imaging from HSC, are measured and interpreted with halo occupation distribution and abundance matching models. Both models successfully account for the observations. We obtain an exceptionally complete view of the AGN halo occupation distribution. The population of AGNs is broadly distributed among halos with a mean mass of 3.9 -2.4 +2.0 × 10 12 M ⊙ . This corresponds to a large-scale halo bias of b(z = 0.34) = 0.99 -0.10 +0.08 . The central occupation has a large transition parameter, σ log 10 (M) = 1.28 ± 0.2. The satellite occupation distribution is characterized by a shallow slope, α sat = 0.73 ± 0.38. We find that AGNs in satellites are rare, with f sat < 20%. Most soft X-ray-selected AGNs are hosted by central galaxies in their dark matter halo. A weak correlation between soft X-ray luminosity and large-scale halo bias is confirmed (3.3σ). We discuss the implications of environmental-dependent AGN triggering. This study paves the way toward fully charting, in the coming decade, the coevolution of X-ray AGNs, their host galaxies, and dark matter halos by combining eROSITA with SDSS-V, 4MOST, DESI, LSST, and Euclid data.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic simulation of a 50 mTorr capacitively coupled argon discharge over a range of frequencies and comparison to experiments

The effect of driving frequency in the range of 13.56–73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct “temperatures” for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2–10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25% of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model, and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coherent fluctuations in a partially magnetized small magnetron discharge operating on molecular nitrogen

We report experimental data on rotating plasma disturbances in a quasi-steady direct-current magnetron discharge operating with molecular nitrogen at pressures of 27, 40, 67, and 133 Pa. Experiments with a segmented anode measuring the local current flow reveal up to 12 coherent structures (“spokes”) propagating in the +E×B and −E×B azimuthal directions. The current–voltage characteristics are studied and found to follow a negative resistance profile at the lower discharge currents. The existence of coherent structures and the number of spokes observed depend on the discharge current and operating pressure. Lower pressures are more favorable to the retrograde −E×B spoke rotation, which is observed in the 27, 40, and 67 Pa discharges. The results are interpreted with a gradient-driven drift wave model, and using the dispersion relation, we show that the propagation direction depends on the orientation of the local electric field within the plasma. Using an electron fluid model, we show how experimental conditions can lead to field reversals that generate −E×B spoke rotation in the lower pressure discharges.

Physics↗

Effects of drifts on scrape-off layer transport in W7-X

Abstract Drifts affect particle, momentum, and energy transport in the scrape-off layer (SOL) of tokamaks and stellarators, altering plasma flows and creating asymmetries between divertors. To understand how drifts affect SOL transport in the W7-X island divertor, an experiment was performed to compare plasmas with matched core parameters but opposite magnetic field directions, and therefore opposite drift transport directions. Parallel flow measurements made with coherence imaging spectroscopy are interpreted with the aid of a diagnostic forward model and a 1D simple SOL model that includes the E × B drift. In low-density plasmas ( n ‾ e < 2 × 10 19 m −3 ), the poloidal E × B drift induces a large poloidal density asymmetry within the island SOL, as measured by divertor Langmuir probes. This in turn causes the parallel flow stagnation point to shift from the position halfway between targets to the X-point in the drift direction, leading to near-unidirectional flow throughout the SOL. As density increases, the effects of the poloidal E × B drift decrease substantially, resulting in a smaller density asymmetry and the development of a counter-streaming flow pattern. For the entire density range probed in this experiment ( n ‾ e = 1.5 − 6 × 10 19 m −3 ), the experimental observations are more consistent with the effects of the poloidal E × B drift than the radial E × B drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temperature-dependent magnetic-field-induced ab -plane Gd spin reorientations in antiferromagnetic domains in a single crystal of the tetragonal A-type antiferromagnet GdRh 2 Si 2

The compound GdRh 2 Si 2 crystallizes in the body-centered tetragonal ThCr 2 Si 2 structure and exhibits A-type antiferromagnetic order below its Néel temperature T N = 107 K where the Gd moments are ferromagnetically aligned in the ab-plane with the moments in adjacent planes aligned antiferromagnetically. Previously, the nonlinear low-field magnetization versus field M ab ( H ∥ [1, 0, 0]) measurements on single crystals were interpreted in terms of a mean-field model. Here, in this work, we reinterpret this nonlinearity more simply in terms of the response of the ordered moments in four noninteracting tetragonal domains to the applied field and find good agreement with our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Search in diphoton and dielectron final states for displaced production of Higgs or Z bosons with the ATLAS detector in s = 13 TeV p p collisions

A search is presented for displaced production of Higgs bosons or Z bosons, originating from the decay of a neutral long-lived particle (LLP) and reconstructed in the decay modes H → γ γ and Z → ee . The analysis uses the full Run 2 dataset of proton-proton collisions delivered by the LHC at an energy of $\sqrt{s}$ = 13 TeV between 2015 and 2018 and recorded by the ATLAS detector, corresponding to an integrated luminosity of 139 fb -1 . Exploiting the capabilities of the ATLAS liquid argon calorimeter to precisely measure the arrival times and trajectories of electromagnetic objects, the analysis searches for the signature of pairs of photons or electrons which arise from a common displaced vertex and which arrive after some delay at the calorimeter. The results are interpreted in a gauge-mediated supersymmetry breaking model with pair-produced Higgsinos that decay to LLPs, and each LLP subsequently decays into either a Higgs boson or a Z boson. The final state includes at least two particles that escape direct detection, giving rise to missing transverse momentum. No significant excess is observed above the background expectation. The results are used to set upper limits on the cross section for Higgsino pair production, up to a $\bar{χ}^0_1$ mass of 369 (704) GeV for decays with 100% branching ratio of $\bar{χ}^0_1$ to Higgs (Z) bosons for a $\bar{χ}^0_1$ lifetime of 2 ns. A model-independent limit is also set on the production of pairs of photons or electrons with a significant delay in arrival at the calorimeter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coincident learning for beam-based rf station fault identification using phase information at the SLAC linac coherent light source

Anomalies in radio-frequency (rf) stations can result in unplanned downtime and performance degradation in linear accelerators such as SLAC’s Linac Coherent Light Source (LCLS). Detecting these anomalies is challenging due to the complexity of accelerator systems, high data volume, and scarcity of labeled fault data. Prior work identified faults using beam-based detection, combining rf amplitude and beam position monitor data. Due to the simplicity of the rf amplitude data, classical methods are sufficient to identify faults, but the recall is constrained by the low-frequency and asynchronous characteristics of the data. In this work, we leverage high-frequency, time-synchronous rf phase data to enhance anomaly detection in the LCLS accelerator. Due to the complexity of phase data, classical methods fail, and we instead train deep neural networks within the Coincident Anomaly Detection (CoAD) framework. We find that applying CoAD to phase data detects nearly 3 times as many anomalies as when applied to amplitude data, while achieving broader coverage across rf stations. Furthermore, the rich structure of phase data enables us to cluster anomalies into distinct physical categories. Through the integration of auxiliary system status bits, we link clusters to specific fault signatures, providing additional granularity for uncovering the root cause of faults. We also investigate interpretability via Shapley values, confirming that the learned models focus on the most informative regions of the data and providing insight for cases where the model makes mistakes. This work demonstrates that phase-based anomaly detection for rf stations improves both diagnostic coverage and root cause analysis in accelerator systems and that deep neural networks are essential for effective analysis.

Accelerator Physics (physics.acc-ph)↗

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

54 ENVIRONMENTAL SCIENCES↗

Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation

The extensive heterogeneity of biological data poses challenges to analysis and interpretation. Construction of a large-scale mechanistic model of Escherichia coli enabled us to integrate and cross-evaluate a massive, heterogeneous dataset based on measurements reported by various groups over decades. We identified inconsistencies with functional consequences across the data, including that the total output of the ribosomes and RNA polymerases described by data are not sufficient for a cell to reproduce measured doubling times, that measured metabolic parameters are neither fully compatible with each other nor with overall growth, and that essential proteins are absent during the cell cycle—and the cell is robust to this absence. Finally, considering these data as a whole leads to successful predictions of new experimental outcomes, in this case protein half-lives.

59 BASIC BIOLOGICAL SCIENCES↗

INTERDUNE DEPOSITS AND STRATIGRAPHIC FRAMEWORK OF THE NAVAJO SANDSTONE IN THE LAKE POWELL AREA, UTAH, USA

The Western Interior of the USA was dominated by desert environments for much of the Jurassic. Of these desert environments, the most extensive by far was the Navajo Sandstone erg, which, in conjunction with its correlative Nugget and Aztec Sandstones, is thought to have covered well over 500,000 km2, and is one of the largest ergs preserved in the rock record. The Navajo Sandstone is also one of a few geologic units which show promise as a potential CO2 sequestration reservoir in Utah. It outcrops along the eastern and southern margin of the state, and subcrops throughout most of the central area. The Navajo Sandstone consists of stacked eolian deposits, composed of well sorted, clean sandstone, deposited in commonly tabular, high angle cross bedded paleodunes, which are on average 15 m thick. Interspersed between these paleodunes are interdune deposits. Interdune deposits consist of sandy to calcareous, predominantly planar bedded, lensoid units, ranging from <1 m up to 7 m thick. In some cases, these interdune deposits occur along more extensive stratigraphic bounding surfaces. Dating within the Navajo Sandstone has historically been very difficult, as a result of its predominantly clean, aeolian sandstone composition. In addition, complexity of internal stratigraphic surfaces coupled with a low structural dip and generally flat topography has resulted in a poor understanding of the internal stratigraphy of the Navajo Sandstone. Exceptional three-dimensional exposures of the Navajo Sandstone in the upper portions of Lake Powell permit us to relate interdune carbonates to a hierarchical bounding surfaces model. Here we use this combined approach to interpret the internal stratigraphy of the Navajo Sandstone in this region. This work forms part of a larger study of the Navajo Sandstone across Utah and its potential as a CO2 reservoir, by developing a more comprehensive stratigraphic framework, and increasing understanding of stratigraphic complexity within one of the most significant eolian systems in the world.

Mahon, Elizabeth↗

Effect of heating/cooling dynamics in the hysteresis loop and tunable IR emissivity of VO 2 thin films

We experimentally investigate the semiconductor-to-metal transition (SMT) in vanadium dioxide thin films using an infrared thermographic technique. During the semiconductor to metal phase change process, VO 2 optical properties dynamically change and infrared emission undergoes a hysteresis loop due to differences between heating and cooling stages. The shape of the hysteresis loop was accurately monitored under different dynamic heating/cooling rates. In order to quantify and understand the effects of different rates, we used a numerical modelling approach in which a VO 2 thin layer was modeled as metamaterial. The main experimental findings are interpreted assuming that both the rate of formation and shape of metallic inclusions are tuned with the heating/cooling rate. The structural transition from monoclinic to tetragonal phases is the main mechanism for controlling the global properties of the phase transition. However, our experimental results reveal that the dynamics of the heating/cooling process can become a useful parameter for further tuning options and lays out a macroscopic optical sensing scheme for the microscopic phase change dynamics of VO 2 . Our study sheds light on phase-transition dynamics and their effect on the infrared emission spectra of VO 2 thin films, therefore enabling the heating/cooling rate to be an additional parameter to control infrared emission characteristics of thermal emitters. The hysteresis loop represents the phase coexistence region, thus being of fundamental importance for several applications, such as the operation of radiative thermal logic elements based on phase transition materials. For such applications, the phase transition region is shifted for heating and cooling processes. We also show that, depending on the way the phase change elements are heated, the temperature operation range will be slightly modified.

42 ENGINEERING↗

TX$^2$: Transformer eXplainability and eXploration

The Transformer eXplainability and eXploration (Martindale & Stewart, 2021), or TX 2 software package, is a library designed for artificial intelligence researchers to better understand the performance of transformer models (Vaswani et al., 2017) used for sequence classification. The tool is capable of integrating with a trained transformer model and a dataset split into training and testing populations to produce an ipywidget (Project Jupyter Contributors, 2021) dashboard with a number of visualizations to understand model performance with an emphasis on explainability and interpretability. The TX 2 package is primarily intended to integrate into a workflow centered around Jupyter Notebooks (Kluyver et al., 2016), and currently assumes the use of PyTorch (Paszke et al., 2019) and Hugging Face transformers library (Wolf et al., 2020). The dashboard includes visualization and data exploration features to aid researchers, including an interactive UMAP embedding graph (McInnes et al., 2018) to understand classification clusters, a word salience map that can be updated as researchers alter textual entries in near real time, a set of tools to understand word frequency and importance based on the clusters in the UMAP embedding graph, and a set of traditional confusion matrix analysis tools.

97 MATHEMATICS AND COMPUTING↗