Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Explanability”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Critical behavior in a chiral molecular model

Understanding the condensed-phase behavior of chiral molecules is important in biology as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature of T = 0.985Tc, and our findings indicate that the free energy barrier for this transformation increases with system size as N2/3, where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field.

Chemistry↗
DOE OSTI2023

Calculated solvent reorganization entropies, free energies, and fluctuations of the energy gaps for intramolecular electron transfer and excitation of the solvatochromic dye B30

Intramolecular electron transfer between two biphenyl groups linked by an androstane spacer and excitation of the pyridinium-N-phenolate betaine dye B30 to the first excited singlet state are studied by quantum/classical molecular-dynamics simulations at temperatures between 150 and 300 K in solvents with a range of polarities. Temperature dependences of the solvent reorganization energies, free energies, entropies, and the inhomogeneous broadening of B30’s absorption band are examined. The variances of fluctuations of the energy gap between the reactant and product states do not have the direct proportionality to temperature that often is assumed to hold. An explanation for the observations is suggested.

Chemistry↗
DOE OSTI2023

Laboratory study of the PFRC-2's initial plasma densification stages

Initial plasma densification by odd-parity rotating magnetic fields (RMF o ) applied to the linear magnetized Princeton field-reversed configuration (PFRC-2) device with fill gases at pressures near 1 mTorr proceeds through two phases: a slow one, characterized by a rise time $τ_{slow}$ ~ 100 $μ$s, followed by a fast one, characterized by $τ_{fast}$ ~ 10 $μ$s. The transition from slow to fast occurs at a line-integral-averaged electron density, t n e , near 2$\times$ 10 11 cm –3 , independent of magnetic field. Here, over most of the range of experimental parameters investigated, as the PFRC-2 axial magnetic field strength was increased, RMF o power decreased, gas fill pressure lowered, or lower atomic mass unit (AMU) fill gas used, the duration of the slow phase lengthened from 50 $μ$s to longer than 10 ms after the RMF o power began. The post-fast-phase maximum n e increases with the fill-gas AMU, exceeding 5 × 10 13 cm –3 for Ar. The slow phase is consistent with atomic physics processes and field-parallel sound-speed losses. The fast phase may be explained by improved axial confinement, possibly augmented by radial or axial contraction of the plasma. Another possible explanation, a large increase in electron temperature, is inconsistent with x-ray emission. The n e behavior is discussed in relation to the E to H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i > a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρ i effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Long-term relaxation of orientational disorder and structural modifications in molecular nitrogen at high pressure

Up to 17 GPa, the crystalline phases of N 2 are characterized by pronounced orientational disorder, whereas the higher-pressure phases of molecular N 2 are ordered. This raises the question about long-term relaxation of orientational disorder within the low- to intermediate-pressure regime. Here, in this study, this question is addressed by comparing synthetic with natural, chemically pure, solid N 2 that resides as inclusions in diamonds at 300 K for about 10 8 years at pressures up to 11 GPa. It is shown that disorder prevails at 8.7 GPa, 300 K, where both synthetic and natural N 2 assume the same structure. However, at 10.8 GPa, natural solid N 2 exhibits monoclinic distortion and partial orientational ordering of the molecules, both of which are not observed in synthetic material. This difference is interpreted as the result of long-term structural relaxation. The ordering mechanism is examined and placed into the context of the δ- to ε-N 2 transition. We present explanations for the absence of complete ordering of δ-N 2 .

Tschauner, Oliver↗
DOE OSTI2024

Characterization of the response of radiochromic film to quasi-monoenergetic x rays through a cross-calibration with image plates

Radiochromic film (RCF) and image plates (IPs) are both commonly used detectors in diagnostics fielded at inertial confinement fusion (ICF) and high-energy-density physics (HEDP) research facilities. Due to the intense x-ray background in all ICF/HEDP experiments, accurately calibrating the optical density of RCF as a function of x-ray dose, and the photostimulated luminescence per photon of IPs as a function of x-ray energy, is necessary for interpreting experimental results. Various measurements of the sensitivity curve of different IPs to x rays have been performed [Izumi et al., Proc. SPIE 8850, 885006 (2013) and Rosenberg et al., Rev. Sci. Instrum. 90(1), 013506 (2019)]; however, calibrating RCF is a tedious process that depends on factors such as the orientation in which the RCF is scanned in the film scanner and the batch of RCF used. These issues can be mitigated by cross-calibrating RCF with IPs to enable the use of IPs for the determination of dose on the RCF without scanning the RCF. Here, the first cross-calibration of RCF with IPs to quasi-monoenergetic titanium, copper, and molybdenum K-line x rays is presented. It is found that the IP-inferred dose rates on the RCF for the Ti and Mo x rays agree well with the measured dose rates, while the IP-inferred dose rate for the Cu x rays is larger than the measured dose rate by ∼2×. Explanations for this discrepancy and plans for future work are discussed.

Instruments & Instrumentation↗
DOE OSTI2024

Proton rings from late-forming ballistic sheath fields

Many laser-driven ion experiments have seen ring-like patterns in the proton angular distribution across a wide array of laser and target parameters. These rings can impede measurement due to the small acceptance angle of detectors and often inhibit potential applications. A myriad of explanations for their formation have been proposed, yet most studies attribute them to some aspect of the laser–plasma interaction. Using 3D particle-in-cell simulations, we show that late-forming strong radial electric fields can arise due to charge separation while the beam is in flight, long after the laser–plasma interaction. These fields can accelerate ions to significant divergences (≈ 10°) as they propagate away from the target. We compare our results to a recent experiment where a high intensity, short pulse laser (I 0 ≈ 10 21 W/cm 2 , τ ≈ 30 fs) was incident upon thin (≈ 1 μm) liquid crystal targets. Our simulations capture all the main features of the experimental results—namely, robust ring formation and larger rings for higher energy protons. In addition, we show that rings do not form for sufficiently short preplasma scale lengths. Finally, we develop a phenomenological model to describe the spatiotemporal structure of the radial electric field and use this to explain the proton rings' energy and preplasma dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗
DOE OSTI2024

Advances in laser-based bremsstrahlung x-ray sources. II. Laser pulse propagation and guiding in nonuniform plasma media in the presence of self-focusing

An analytic Wentzel–Kramers–Brillouin model is presented of Gaussian laser pulse propagation through plasma with a quadratic transverse density profile and an arbitrarily varying, longitudinal density gradient under conditions of nonlinear self-focusing. From these solutions, it is shown that in the absence of nonlinear self-focusing and transverse nonuniformity, for exponential pre-plasma density profiles, the use of a low density coating of the laser target with electron density n0∼11 ncr (e.g., a CH foam of density 35 mg/cm3 for 1-micron laser light) maximizes laser intensity at best focus. Also, under laser and plasma conditions relevant to recent experiments on high-power laser systems, conditions are obtained for a Gaussian laser pulse to propagate stably through the pre-plasma medium. Such conditions would be expected to enhance the production of relativistic electrons from laser-target coupling, providing a possible explanation for the observed increase in MeV photon dose and enabling applications such as laser-based MeV X-ray radiography.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Detailed simulations of the first deuterium-tritium-filled double shell implosions on the National Ignition Facility

The first indirectly driven, liquid DT-filled double shell inertial confinement fusion (ICF) implosions have recently been successfully performed on the National Ignition Facility (NIF). Double shells are a class of alternative designs that use a low-Z outer shell to compress a foam cushion that accelerates a high-Z inner shell to efficiently compress a liquid DT core. Double shells are challenging to fabricate, field, and model. Important engineering features enabling double shell fabrication include a fill-tube penetrating all shells and a carefully designed and very narrow (few μm) step-joint in the ablator. Due to the higher density materials involved, high Atwood number instabilities are also important at many material interfaces. In this paper, numerical simulations of double shell implosions using the Los Alamos National Laboratory multi-physics radiation-hydrodynamics code xRAGE will be discussed. An extensive effort has been under way for several years to develop the code capabilities for ICF simulations in a common modeling framework to allow ease of simulation setup and standardization of the computational methodology. This paper will present a wide range of simulation results capturing, quantifying, and comparing the impact of all these degradation mechanisms on implosion performance. Brief comparisons with recent experimental results and suggestions for future improvements will also be discussed. Our results suggest that capsule surface roughness and the step-joint gap have the largest impact on implosion performance. Initial experimental data may suggest that the sensitivity to the step-joint gap could provide the dominant explanation for DT-filled double shell experiments that have been fielded on NIF thus far.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Characterization of electron spin polarization from positive electron affinity GaAs photocathodes

Negative Electron Affinity (NEA) GaAs photocathodes are widely used to generate spin-polarized electron beams, typically achieving Electron Spin Polarizations (ESPs) in the range of 35%–40%. However, when operated in a Positive Electron Affinity (PEA) state, where a potential barrier inhibits low-energy electrons from escaping into vacuum, measured ESPs can exceed 50%. This effect can occur naturally during photocathode operation, as the NEA surface activation layers can degrade easily over time, increasing the electron affinity. In this work, we investigate and characterize the behavior of enhanced ESPs under PEA conditions. We present experimental measurements of ESP using a retarding-field Mott polarimeter on GaAs photocathodes with controlled and varying electron affinities. These results are complemented by theoretical explanations considering the material band structure, the light excitation profile, and spin depolarization mechanisms.

Electronic band structure↗
DOE OSTI2026

Hot electron preheat effects in two-photon polymerization plastic lattice laser-driven shock-tube experiments on the OMEGA laser

We present simulations of Omega shock tube experiments designed to investigate hot electron preheat effects in 3D-printed, two-photon polymerization (2PP) plastic lattices. Preheat is inferred in the experiments from the expansion of a plastic witness disk embedded in the lattice. Using the Eulerian radiation-hydrodynamics code xRAGE, we model shock propagation and preheat from both radiative and hot electron energy sources to evaluate their relative impact. To simulate the transport of laser-generated hot electrons, the nonlocal electron heat transport model proposed by Schurtz, Nicolaï, and Busquet (SNB) is extended with a hot electron source term and an energy cascade algorithm. We explore how variations in ablator, lattice geometry, and laser drive affect the shock velocity and witness disk expansion. Simulations show that the inclusion of a 5 μm gold layer reduces shock pressure by 60% and shock speeds by 30%–40% but does not significantly reduce the hot electron preheat, and that different lattice geometries lead to enhanced shock velocities—up to 40% faster than in homogeneous foams. However, radiative and conductive preheat from classical mechanisms alone fail to match experiment. By including a hot electron source term, we reproduce experimental observables such as disk expansion rates and spatial radiographic features. We find that a hot electron population corresponding to 4%–8% of the incident laser energy with T hot = 50 keV produces expansion which agrees with the experimental data, suggesting hot electron preheat is the most plausible explanation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Analytic rate theory of polariton relaxation that explains long polariton lifetime

Hybridization of a molecular exciton with a quantized photon creates a polariton. Despite extensive experimental investigations, the apparent lifetime of the exciton–polariton is not well-understood. For this work, we examined the steady-state population dynamics for a Holstein–Tavis–Cumming Hamiltonian to illuminate the long-term polaritonic dynamics and lifetime of the exciton–polariton in an optical cavity. For a realistic description of polariton relaxation, cavity loss and various exciton decay channels are included in the model. We found that in the presence of weak but finite exciton loss, the apparent lifetime of the lower polariton coincides with the out-of-cavity exciton lifetime and is independent of cavity-matter detuning. This is a simple explanation for the experimentally observed lifetimes for exciton polaritons and theoretically justifies the dark state reservoir hypothesis. Furthermore, if the upper polariton is initially populated, the system reaches the steady state very quickly, leading to single-exponential polariton relaxation. Starting from the lower polariton leads to a longer pre-steady-state time period, leading to double-exponential relaxation. Finally, we considered the effect of site orientational disorders and the exciton frequency disorderers. Under the collective limit, the effects of this disorder can be included in Fermi’s golden rule population dynamics without explicit sampling. For the exciton energy disorders, numerical calculations are needed. Our theoretical framework is applicable to interpret exciton–polariton experiments, especially related to the measured apparent lifetime of polaritons.

Chemical dynamics↗
DOE OSTI2026

Attention-based explainability for structure–property relationships

Machine learning methods are emerging as a universal paradigm for constructing correlative structure–property relationships in materials science based on multimodal characterization. However, this necessitates the development of methods for the physical interpretability of the resulting correlative models. Here, we demonstrate the potential of attention-based neural networks for revealing structure–property relationships and the underlying physical mechanisms, using the ferroelectric properties of PbTiO3 thin films as a case study. Through the analysis of attention scores, we disentangle the influence of distinct domain patterns on the polarization switching process. The attention-based Transformer model is explored both as a direct interpretability tool and as a surrogate for explaining representations learned via unsupervised machine learning, enabling the identification of physically grounded correlations. We compare attention-derived interpretability scores with classical SHapley Additive exPlanations analysis and show that, in contrast to applications in natural language processing, attention mechanisms in materials science exhibit high efficiency in highlighting meaningful structural features.

Slautin, Boris [Independent Researcher]↗
DOE OSTI2026

Chaos as the cause of randomness in tearing mode onset times in DIII-D ITER baseline scenario plasmas

The first evidence of widespread chaotic plasma dynamics relevant to the onset of disruptive tearing modes is found in DIII-D tokamak discharges. This evidence is obtained by calculating positive maximal Lyapunov exponents from Mirnov signals using the Rosenstein algorithm, implying the chaotic growth and/or rotation of magnetohydrodynamic modes in DIII-D [Luxon, Nucl. Fusion 42, 614 (2002)] plasmas. Such chaotic dynamics offer an explanation for the random tearing mode onset times previously observed in ITER baseline scenario discharges. Further, the Lyapunov timescales associated with this chaos could inhibit reliable prediction of tearing mode onset beyond the angular momentum confinement timescale; this constraint should inform the design of plasma control systems on ITER.

Chaotic dynamics↗
DOE OSTI2026

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗