SEARCH · Search NASA

Results for “Electron Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Probing cold nuclear matter with energy correlators

The future electron-ion collider (EIC) will produce the first-ever high energy collisions between electrons and a wide range of nuclei, opening a new era in the study of cold nuclear matter. Quarks and gluons produced in these collisions will propagate through the dense nuclear matter of nuclei, imprinting its structure into subtle correlations in the energy flux of final state hadrons. In this article, we apply recent developments from the field of jet substructure, namely the energy correlator observables, to decode these correlations and provide a new window into nuclear structure. The energy correlators provide a calibrated probe of the scale dependence of vacuum quantum chromodynamics (QCD), enabling medium modifications to be imaged and interpreted as a function of scale. Using the eHIJING parton shower to simulate electron-nucleus collisions, we demonstrate that the size of the nucleus is imprinted as an angular scale in the correlators, with a magnitude that is visible for realistic EIC kinematics. We can observe the size difference between the proposed EIC nuclear targets 3 He, 4 He, 12 C, 40 Ca, 64 Cu, 197 Au, and 238 U, showing that the energy correlators can image femtometer length scales using asymptotic energy flux. Our approach offers a unified view of jet substructure across collider experiments, and provides numerous new theoretical tools to unravel the complex dynamics of QCD in extreme environments, both hot and cold.

Electron-ion collisions↗

NASA NTRS2010

Composition Changes After the "Halloween" Solar Proton Event: The High-Energy Particle Precipitation in the Atmosphere (HEPPA) Model Versus MIPAS Data Intercomparison Study

We have compared composition changes of NO, NO2, H2O2,O3, N2O, HNO3 , N2O5, HNO4, ClO, HOCl, and ClONO2 as observed by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on Envisat in the aftermath of the "Halloween" solar proton event (SPE) in October/November 2003 at 25-0.01 hPa in the Northern hemisphere (40-90 N) and simulations performed by the following atmospheric models: the Bremen 2D model (B2dM) and Bremen 3D Chemical Transport Model (B3dCTM), the Central Aerological Observatory (CAO) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the ECHAM5/MESSY Atmospheric Chemistry (EMAC) model, the modeling tool for SO1ar Climate Ozone Links studies (SOCOL and SOCOLi), and the Whole Atmosphere Community Climate Model (WACCM4). The large number of participating models allowed for an evaluation of the overall ability of atmospheric models to reproduce observed atmospheric perturbations generated by SPEs, particularly with respect to NOS, and ozone changes. We have further assessed the meteorological conditions and their implications on the chemical response to the SPE in both the models and observations by comparing temperature and tracer (CH4 and CO) fields. Simulated SPE-induced ozone losses agree on average within 5% with the observations. Simulated NO(y) enhancements around 1 hPa, however, are typically 30% higher than indicated by the observations which can be partly attributed to an overestimation of simulated electron-induced ionization. The analysis of the observed and modeled NO(y) partitioning in the aftermath of the SPE has demonstrated the need to implement additional ion chemistry (HNO3 formation via ion-ion recombination and water cluster ions) into the chemical schemes. An overestimation of observed H2O2 enhancements by all models hints at an underestimation of the OH/HO2 ratio in the upper polar stratosphere during the SPE. The analysis of chlorine species perturbations has shown that the encountered differences between models and observations, particularly the underestimation of observed ClONO2 enhancements, are related to a smaller availability of ClO in the polar night region already before the SPE. In general, the intercomparison has demonstrated that differences in the meteorology and/or initial state of the atmosphere in the simulations causes a relevant variability of the model results, even on a short timescale of only a few days.

Funke, B.↗

DOE OSTI2025

Fast calculation of diffraction patterns from an ensemble of aligned molecules

We report an algorithm to calculate electron diffraction patterns for molecules with anisotropic angular distribution, which is significantly faster than existing methods. The algorithm uses a transform to convert the molecular orientation distribution, which is a function of three Euler angles, to the atom-pair distribution functions which depend on the polar and azimuthal angles. The diffraction signal can then be calculated from the atom-pair distributions. We demonstrate the computation method numerically by calculating electron diffraction patterns for a symmetric top molecule (trifluoroiodomethane) and an asymmetric top molecule (formaldehyde) and show that it reduces the calculation time by approximately two orders of magnitude compared to the standard brute-force method. Here, the method can also be applied to the calculation of x-ray diffraction patterns.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1975

The width of electron-photon cascades in air

Results of a Monte-Carlo simulation of electron-photon cascades are presented with particular emphasis on cascade lateral spread. The effects of single and multiple Coulomb scattering are represented by giving the particles suitable displacements in position and angle at each interaction point. Some of the more oblique particle tracks are first divided into a number of steps between interaction points and the scattering procedure is applied to each step separately. To represent the particle defections due the geomagnetic field, a field of 0.3 gauss is supposed present and the appropriate helical path is imposed on each particle between interaction points. The results are in good agreement with the calculations of Nordheim et al., with simulations using the Berger-Seltzer electron photon transport program and with experimental measurements in water.

Allan, H. R.↗

NASA NTRS2010

Estimating the Probability of Electrical Short Circuits from Tin Whiskers

To comply with lead-free legislation, many manufacturers have converted from tin-lead to pure tin finishes of electronic components. However, pure tin finishes have a greater propensity to grow tin whiskers than tin-lead finishes. Since tin whiskers present an electrical short circuit hazard in electronic components, simulations have been developed to quantify the risk of said short circuits occurring. Existing risk simulations make the assumption that when a free tin whisker has bridged two adjacent exposed electrical conductors, the result is an electrical short circuit. This conservative assumption is made because shorting is a random event that had an unknown probability associated with it. Note however that due to contact resistance electrical shorts may not occur at lower voltage levels. In our first article we developed an empirical probability model for tin whisker shorting. In this paper, we develop a more comprehensive empirical model using a refined experiment with a larger sample size, in which we studied the effect of varying voltage on the breakdown of the contact resistance which leads to a short circuit. From the resulting data we estimated the probability distribution of an electrical short, as a function of voltage. In addition, the unexpected polycrystalline structure seen in the focused ion beam (FIB) cross section in the first experiment was confirmed in this experiment using transmission electron microscopy (TEM). The FIB was also used to cross section two card guides to facilitate the measurement of the grain size of each card guide's tin plating to determine its finish .

Courey, Karim J.↗

NASA NTRS2009

Developing an Empirical Model for Estimating the Probability of Electrical Short Circuits from Tin Whiskers

Courey, Karim J.↗

DOE OSTI2025

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations↗

NASA NTRS2006

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves, Precipitating Ring Current Ions, and Thermal Electron Heating - 2

This paper is dedicated to further presentations and discussions of the results from our new global self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2006; here referred to as Paper 1]. In order to adequately take into account the wave propagation and refraction in a multi-ion plasmasphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation [for details see Paper 1]. To demonstrate the effects of the EMIC wave propagation and refraction on the RC proton precipitations and heating of the thermal plasmaspheric electrons we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. Firstly, the wave induced precipitations have a quite fine structure, and are highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 (raised dot) 10(exp 6) [(cm (raised dot) s (raised dot) sr)(sup -l)] are observed during the main and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not simply connected to the most intense EMIC waves. The character of the EMIC wave power spectral density distribution over the equatorial wave normal angle is an extremely crucial for the effectiveness of the RC ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from other ring current model [Kozyra et al., 1997] reveals that although we observe a qualitative agreement between localizations of the wave induced fluxes in the models, there is no quantitative agreement between the magnitudes of these fluxes. These differences are mainly due to a qualitative difference between the characters of the EMIC wave power spectral density distributions over the equatorial wave normal angle. Finally, the two energy sources to the plasmaspheric electrons are considered; (i) the heat fluxes caused by the EMIC wave energy absorption due to Landau resonance, and (ii) the heat fluxes due to Coulomb energy degradation of the RC o(+) ions. The heat fluxes caused by the EMIC wave energy absorption due to Landau resonance are observed in the postnoon-premidnight MLT sector, and maximize at the magnitude of 10l1 (eV/(cm(sup 2)(raised dot) s) at L=3.25, MLT=22 at 3400 UT after 1 May, 0000 UT. The greatest Coulomb energy deposition rates are about 2 (raised dot) 10(sup 10)(eV/(cm(sup 2)(raised dot) s) and observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the RC with plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, the RC, and the topside ionosphere [Gurgiolo et al., 20051.

Khazanov, G. V.↗

NASA NTRS2007

Self-consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves. 2. Wave Induced Ring Current Precipitation and Thermal Electron Heating

This paper continues presentation and discussion of the results from our new global self-consistent theoretical model of interacting ring current ions and propagating electromagnetic ion cyclotron waves [Khazanov et al., 2006]. To study the effects of electromagnetic ion cyclotron wave propagation and refraction on the wave induced ring current precipitation and heating of the thermal plasmaspheric electrons, we simulate the May 1998 storm. The main findings after a simulation can be summarized as follows. Firstly, the wave induced ring current precipitation exhibits quite a lot of fine structure, and is highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 x 10(exp 6) (cm(raised dot) s(raised dot) sr(raised dot) (sup -1)) are observed during the maill and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not connected to the most intense waves in simple manner. The characteristics of the wave power spectral density distribution over the wave normal angle are extremely crucial for the effectiveness of the ring current ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from RAM [Kozyra et al., 1997a] reveals that although we observe a qualitative agreement between the localizations of the wave induced precipitations in the models, there is no quantitative agreement between the magnitudes of the fluxes. The quantitative differences are mainly due to a qualitative difference between the characteristics of the wave power spectral density distributions over the wave normal angle in RAM and in our model. Thirdly, the heat fluxes to plasmaspheric electrons caused by Landau resonate energy absorption from electromagnetic ion cyclotron waves are observed in the postnoon-premidnight MLT sector, and can reach the magnitude of 10(exp 11) eV/(cm(sup 2)(raised dot)s). The Coulomb energy degradation of the RC H(+) and O(+) ions maximizes at about 10(exp 11) (eV/(cm(sup 2) (raised dot) s), and typically leads to electron energy deposition rates of about 2(raised dot) 10(exp 10) (eV/(cm(sup 2)(raised dot)s) which are observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the ring current with the plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, ring current, and topside ionosphere [Gurgiolo et al., 2005]. Finally, the wave induced intense electron heating has a structure of the spot-like patches along the most enhanced density gradients in the plasmasphere boundary layer and can be a possible driver to the observed but still not explained small-scale structures of enhanced emissions in the stable auroral red arcs.

Khazanov, G. V.↗

DOE OSTI2025

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗

DOE OSTI2025

Computational Study of Indium Oxide Photoelectrodes

Using a combination of first-principles molecular dynamics simulations and electronic structure calculations, we characterize the atomistic structure and vibrational properties of a photocatalytic surface of In 2 O 3 , a promising photoelectrode for the production of hydrogen peroxide. We then investigate the surface in contact with water and show that the electronic states of In 2 O 3 are appropriately positioned in energy to facilitate the two-electron water oxidation reaction (WOR) over the competing four-electron oxygen evolution reaction. In conclusion, we further propose that the use of strained thin films interfaced with water is beneficial in decreasing the optical gap of In 2 O 3 and thus in utilizing a wider portion of the solar spectrum for the WOR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Vibronic coupling-driven symmetry breaking and solvation in the photoexcited dynamics of quadrupolar dyes

Quadrupolar dyes, such as acceptor–donor–acceptor molecules, are highly relevant for applications in nonlinear optics and photovoltaics. They are also versatile models for exploring photoinduced charge-transfer dynamics. The interplay between electronic and vibronic couplings in these molecules may break excited-state symmetry, resulting in intramolecular charge separation and pronounced solvatochromism. Experimentally, distinguishing the roles of intramolecular vibronic coupling and solvent reorganization for the initial charge-transfer dynamics has been challenging so far. Here we investigate a prototypical quadrupolar dye in polar and non-polar solvents using ultrafast pump–probe and two-dimensional electronic spectroscopy. Our results reveal that vibronic couplings initiate excited-state symmetry breaking during the first ~50 fs of the photoinduced charge transfer, whereas solvent-induced charge localization sets in at later times. Quantum dynamics and electronic structure simulations support our experimental findings. Our results reveal the details of solvation dynamics in photoexcited molecules and suggest strategies for their manipulation through vibronic couplings.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Magnon spectroscopy in the electron microscope

Abstract The miniaturization of transistors is approaching its limits owing to challenges in heat management and information transfer speed 1 . To overcome these obstacles, emerging technologies such as spintronics 2 are being developed, which make use of the electron’s spin as well as its charge. Local phenomena at interfaces or structural defects will greatly influence the efficiency of spin-based devices, making the ability to study spin-wave propagation at the nanoscale and atomic scale a key challenge 3,4 . The development of high-spatial-resolution tools to investigate spin waves, also called magnons, at relevant length scales is thus essential to understand how their properties are affected by local features. Here we detect bulk THz magnons at the nanoscale using scanning transmission electron microscopy (STEM). By using high-resolution electron energy-loss spectroscopy with hybrid-pixel electron detectors, we overcome the challenges posed by weak signals to map THz magnon excitations in a thin NiO nanocrystal. Advanced inelastic electron scattering simulations corroborate our findings. These results open new avenues for detecting magnons and exploring their dispersions and their modifications arising from nanoscale structural or chemical defects. This marks a milestone in magnonics and presents exciting opportunities for the development of spintronic devices.

Science & Technology - Other Topics↗

DOE OSTI2024

Physics-constrained machine learning for electrodynamics without gauge ambiguity based on Fourier transformed Maxwell’s equations

We utilize a Fourier transformation-based representation of Maxwell’s equations to develop physics-constrained neural networks for electrodynamics without gauge ambiguity, which we label the Fourier–Helmholtz–Maxwell neural operator method. In this approach, both of Gauss’s laws and Faraday’s law are built in as hard constraints, as well as the longitudinal component of Ampère–Maxwell in Fourier space, assuming the continuity equation. An encoder–decoder network acts as a solution operator for the transverse components of the Fourier transformed vector potential, $\hat{A}_⟂(k,t)$, whose two degrees of freedom are used to predict the electromagnetic fields. This method was tested on two electron beam simulations. Among the models investigated, it was found that a U-Net architecture exhibited the best performance as it trained quicker, was more accurate and generalized better than the other architectures examined. We demonstrate that our approach is useful for solving Maxwell’s equations for the electromagnetic fields generated by intense relativistic charged particle beams and that it generalizes well to unseen test data, while being orders of magnitude quicker than conventional simulations. We show that the model can be re-trained to make highly accurate predictions in as few as 20 epochs on a previously unseen data set.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Real-space Kohn–Sham density functional theory for complex energy applications

Real-space Kohn-Sham density functional theory (real-space KS-DFT) enables large-scale electronic structure simulations that is particularly well-suited for the modern high-performance computing (HPC) architectures. This feature article reviews its theoretical foundations, highlights the algorithmic advances and recent developments, and showcases applications in complex nano systems. We aim to provide a perspective on the trajectory of real-space KS-DFT as an emerging tool for computational chemistry and materials science in the exascale era.

Zhang, Zeyi↗

DOE OSTI2024

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗

DOE OSTI2025

Unified architecture for quantum lookup tables

Quantum access to arbitrary classical data encoded in unitary black-box oracles underlies interesting data-intensive quantum algorithms, such as machine learning or electronic structure simulation. The feasibility of these applications depends crucially on gate-efficient implementations of these oracles, which are commonly some reversible versions of the Boolean circuit for a classical lookup table. Here, we present a general parametrized architecture for quantum circuits implementing a lookup table that encompasses all prior work in realizing a continuum of optimal trade-offs between qubits, non-Clifford gates, and error resilience, up to logarithmic factors. Our architecture assumes only local 2D connectivity, yet recovers results, with the appropriate parameters, polylogarithmic error scaling. We also identify regimes, such as simultaneous sublinear scaling, in all parameters. These results enable tailoring implementations of the commonly used lookup table primitive to any given quantum device with constrained resources.

quantum circuits↗

DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗