Search NASA⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Beam dynamics in a storage ring with a free-electron laser

In this paper, we study the dynamics of the electron beam in a diffraction limited storage ring with a high-gain free-electron laser. We present in detail the design, including lattice, dynamic aperture, intrabeam scattering, and coherent synchrotron radiation. The simulation of the electron storage ring combined with the free-electron laser is carried out until the electron beam reaches its equilibrium. We find that the initial exponential radiation power is reduced approximately to linear with respect to the undulator length due to the large increase of the energy spread from the free-electron laser.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Fall SULI Report

This fall, I interned under Fan Guo and Xiaocan Li in T-2 at Los Alamos National Laboratory through the SULI program. I investigated how shock thickness is impacted by collisions between particles in plasmas and how particle acceleration in shocks is impacted by the mass ratio for protons and electrons used in simulations. Both investigations were carried out using the VPIC simulation code. I found that collisions do not meaningfully impact shock width and that higher mass ratios lead to somewhat more electron acceleration in shocks. This write-up is in three parts, the first discussing a description of what I was up to personally throughout the fall, the second discussing the results of the collisional shock study, and the third discussing the results of the high mass ratio shock study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Study of Particle Loss in Synchrotron Phase Space Injection for ESR Using Weak-Strong Beam-Beam Simulation with Nonlinear Lattice

In this report, we use tracking simulations to investigate synchrotron phase space injection for electron accumulation in the electron storage ring of the Electron Ion Collider. Our simulation model accounts for both beam-beam interactions and lattice nonlinearities. Specifically, we examine how particle loss is influenced by various parameters. Additionally, we conduct a theoretical analysis and derive an analytical formula for the rapid evaluation of particle loss. Our results demon strate the feasibility of synchrotron phase-space injection for the electron storage ring and provide insights to guide parameter selection for the design of the injection line.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗
DOE OSTI2024

Photoinduced Electron–Nuclear Dynamics of Fullerene and Its Monolayer Networks in Solvated Environments

The recently synthesized monolayer fullerene network in a quasi-hexagonal phase (qHP-C 60 ) exhibits superior electron mobility and optoelectronic properties compared to molecular fullerene (C 60 ), making it highly promising for a variety of applications. However, the microscopic carrier dynamics of qHP-C 60 remain unclear, particularly in realistic environments, which are of significant importance for applications in optoelectronic devices. Unfortunately, traditional ab initio methods are prohibitive for capturing the real-time carrier dynamics of such large systems due to their high computational cost. In this work, we present the first real-time electron–nuclear dynamics study of qHP-C 60 using velocity-gauge density functional tight binding, which enables us to perform several picoseconds of excited-state electron–nuclear dynamics simulations for nanoscale systems with periodic boundary conditions. When applied to C 60 , qHP-C 60 , and their solvated counterparts, we demonstrate that water/moisture significantly increases the electron–hole recombination time in C 60 but has little impact on qHP-C 60 . Our excited-state electron–nuclear dynamics calculations show that qHP-C 60 is extremely unique and enable exploration of time-resolved dynamics for understanding excited-state processes of large systems in complex, solvated environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Wiggler-based electron beam microbunching amplifier study

The wiggler-based amplifier of the microbunching instability in the electron beam is considered and compared to a longitudinal space charge amplifier consisting of the drifts and chicanes. The gain in the magnitude of a seed perturbation of the electron beam density is evaluated analytically as a function of the period of the perturbation, the electron beam parameters, the wiggler parameters, and the amplifier lattice parameters. The theory is compared with the results of the numerical calculations of the gain based on simulations of the electron beam transport through the amplifier lattice. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Assessing Ground State Energy of Molecules and Energy Profile of the NH3 Capturing CO2 System Using the Quantum Computing Algorithms

Molecule size correlates with the number of electrons on electronic energies and strength of anharmonicity on vibrational properties, however, it is challenging to address using classical computing. In this study, variational quantum eigensolver (VQE) algorithm was implemented on a quantum simulator to quantify electronic and vibrational energies and reaction pathways of CO2 + NH3 = NH2COOH. The VQE-based Hartree-Fock-Embedding algorithm was adopted to benchmark electronic energies for a series of molecules (doi.org/10.1063/5.0188249) and quantify the reaction energy profile of the CO2 capture reaction (doi.org/10.1116/5.0137750). The generated reaction profile is in good agreement with the classical high-level Coupled-Cluster-Singles-and-Doubles (CCSD) results. The quantum computing algorithm also helps enhance the calculation of vibrational ground-state energies by considering the many-body coupling using the Vibrational Self-Consistent Field method, providing results for CO2 and NH3 molecules with accuracy comparable to the direct diagonalization method. Our approach indicates quantum computing can be applied to solve practical problems.

Lee, Yueh-Lin↗
DOE OSTI2024

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Roadmap for Molecular Benchmarks in Nonadiabatic Dynamics

Simulating the coupled electronic and nuclear response of a molecule to light excitation requires the application of nonadiabatic molecular dynamics. However, when faced with a specific photophysical or photochemical problem, selecting the most suitable theoretical approach from the wide array of available techniques is not a trivial task. The challenge is further complicated by the lack of systematic method comparisons and rigorous testing on realistic molecular systems. This absence of comprehensive molecular benchmarks remains a major obstacle to advances within the field of nonadiabatic molecular dynamics. A CECAM workshop, Standardizing Nonadiabatic Dynamics: Towards Common Benchmarks, was held in May 2024 to address this issue. This Perspective highlights the key challenges identified during the workshop in defining molecular benchmarks for nonadiabatic dynamics. Specifically, this work outlines some preliminary observations on essential components needed for simulations and proposes a roadmap aiming to establish, as an ultimate goal, a community-driven, standardized molecular benchmark set.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electron Influence on the Parallel Proton Firehose Instability in 10-moment, Multifluid Simulations

Instabilities driven by pressure anisotropy play a critical role in modulating the energy transfer in space and astrophysical plasmas. For the first time, we simulate the evolution and saturation of the parallel proton firehose instability using a multifluid model without adding artificial viscosity. These simulations are performed using a 10-moment, multifluid model with local and gradient relaxation heat-flux closures in high-β proton–electron plasmas. When these higher-order moments are included and pressure anisotropy is permitted to develop in all species, we find that the electrons have a significant impact on the saturation of the parallel proton firehose instability, modulating the proton pressure anisotropy as the instability saturates. Even for lower β's more relevant to heliospheric plasmas, we observe a pronounced electron energization in simulations using the gradient relaxation closure. Our results indicate that resolving the electron pressure anisotropy is important to correctly describe the behavior of multispecies plasma systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗
DOE OSTI2026

Ion-Electron Coupling-Driven Redox Behavior in Metal–Organic Frameworks

Redox-active metal–organic frameworks (MOFs) have long been proposed as electronic transport platforms, yet the microscopic origin of their conductivity remains debated. A theoretical demonstration reveals charge transport in a Zn(pyrazole–naphthalene diimide (NDI)) MOF arising not from delocalized band-like states but from redox hopping between discrete linker sites. Using ab initio molecular dynamics simulations combined with electronic structure analysis, we established a direct link among electron injection, structural reorganization, and transport. Electron accumulation proceeds sequentially and site-selectively from imide and carbonyl groups of the NDI core progressively involving pyrazole N atoms at higher reduction states, through a hierarchy of redox-active sites. In contrast, Zn nodes remain essentially redox-inactive, which confirms their structural role. Density-of-states analysis corroborates a transport regime dominated by linker-centered states with evolving p-character upon reduction, resulting in dynamically reconfigured conduction networks. Real-time trajectories reveal anisotropic linker-to-linker electron transfer modulated by counterion coordination. This cooperative ion–electron regime emerges from potential energy surface collapse into a single low-barrier transition (ΔG ‡ ≈ 45 meV), where ionic and electronic motions evolve adiabatically on the same free-energy landscape. Elucidating redox conductivity in Zn(pyrazole–NDI) MOFs provides a theoretical framework for use in neuromorphic computing and related technologies.

Charge transfer↗
DOE OSTI2026

Machine-learning force-field models for dynamical simulations of metallic magnets

We review recent advances in machine-learning (ML) force-field methods for Landau–Lifshitz–Gilbert simulations of itinerant electron magnets, focusing on their scalability and transferability. Built on the principle of locality, a deep neural-network model is developed to efficiently and accurately predict electron-mediated forces governing spin dynamics. Symmetry-aware descriptors constructed through a group-theoretical approach ensure rigorous incorporation of both lattice and spin-rotation symmetries. The framework is demonstrated using the prototypical s-d exchange model widely employed in spintronics. ML-enabled large-scale simulations reveal novel nonequilibrium phenomena, including anomalous coarsening of tetrahedral spin order on the triangular lattice and the freezing of phase-separation dynamics in lightly hole-doped, strong-coupling square-lattice systems. These results establish ML force-field frameworks as scalable, accurate, and versatile tools for modeling nonequilibrium spin dynamics in itinerant magnets.

Artificial neural networks↗
DOE OSTI2025

Electron-Proton Scattering Event Generation using Structured Tokenization

Recent work such as Omnijet-$\alpha$ has demonstrated that effective tokenization combined with transformer-based architectures can produce effective foundation models for jet physics. While tokenization may help models capture generalizable event characteristics, it also introduces discretization errors that may compromise the precision required for downstream physics analyses. As the number and complexity of the particle features grow, these errors are likely to grow proportionally. In this study, we investigate new tokenization strategies to improve the application of generative transformer models to \textsc{Pythia8} simulations of electron-proton scattering at the Electron-Ion Collider. Specifically, we propose a feature-based structured tokenization approach that utilizes multiple tokens per particle, improving expressivity, while reducing the total number of unique tokens needed. We evaluate this method against grid-based binning, K-means clustering, and vector-quantized variational auto-encoders on the event simulations. Our results show that feature-based structured tokenization reduces discretization error, leading to more accurate generative modeling of particle-level events.

Goldenberg, Steven [Thomas Jefferson National Acce↗
DOE OSTI2025

Electron-Proton Scattering Event Generation using Structured Tokenization

Recent work such as Omnijet-$\alpha$ has demonstrated that effective tokenization combined with transformer-based architectures can produce effective foundation models for jet physics. While tokenization may help models capture generalizable event characteristics, it also introduces discretization errors that may compromise the precision required for downstream physics analyses. As the number and complexity of the particle features grow, these errors are likely to grow proportionally. In this study, we investigate new tokenization strategies to improve the application of generative transformer models to \textsc{Pythia8} simulations of electron-proton scattering at the Electron-Ion Collider. Specifically, we propose a feature-based structured tokenization approach that utilizes multiple tokens per particle, improving expressivity, while reducing the total number of unique tokens needed. We evaluate this method against grid-based binning, K-means clustering, and vector-quantized variational auto-encoders on the event simulations. Our results show that feature-based structured tokenization reduces discretization error, leading to more accurate generative modeling of particle-level events.

Goldenberg, Steven [Thomas Jefferson National Acce↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗