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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electron transport in a tokamak scrape-off layer: impact of toroidal nonuniformities of divertor targets

An analysis of electron dynamics is carried out in the tokamak scrape-off layer (SOL) region. Small non-axisymmetric perturbations of the divertor target sheath potential affect electron drift orbits in the SOL and may lead to non-ambipolar electron radial transport. The resulting rates of electron convection, diffusion, and mobility are calculated analytically, and the analytic results are compared with direct numerical simulations of electron drift orbits. The proposed mechanism of non-ambipolar electron transport may be relevant to the sustainment of the SOL plasma quasi-neutrality in the ‘heuristic’ model of the SOL width (Eich et al 2011 Phys. Rev. Lett. 107 215001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗
DOE OSTI2025

Comparison of laser-produced plasma spatio-temporal electron density evolution measured using interferometry with simulation results

Due to the difficulties associated with experimental measurements of laser-produced plasma (LPP) properties during the earliest stages of plasma evolution, radiation hydrodynamic codes are often used. However, although these codes have been extensively validated in the higher intensity regimes, validation at low to moderate intensities has been limited. In this study, the spatio-temporal electron density evolution in an LPP generated at moderate laser intensities and at various laser wavelengths was validated against the FLASH code for times up to 20 ns. The LPP was generated by focusing the fundamental and various harmonics radiation (1064, 532, and 266 nm) from a 6 ns full width half maximum Nd:YAG laser, at a laser intensity of 10 GW cm -2 , onto a copper target. The spatio-temporal density evolution of the expanding plasma was analyzed using Nomarski interferometry. Experimental measurements were found to be consistent with FLASH simulations, and the dependence of electron densities on wavelength was found to be in agreement with analytical models, varying as n e ∝ λ -0.7 . However, slight differences were noted in the widths and shape functions of the experimental and simulated electron density profiles.

2D fast imaging↗
DOE OSTI2024

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

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

Experimental study of energy-dependent angular broadening of MeV electron beams for high-resolution imaging in thick samples

In scanning transmission electron microscopy (STEM), spatial resolution is primarily influenced by the projected size of the electron probe within the specimen. In thin samples, a large semi-convergence angle enables a tightly focused beam and sub-nanometer resolution. However, in thick specimens, resolution is fundamentally limited by transverse beam broadening from multiple large-angle scattering events—for example, a probe with 10 mrad angular divergence can broaden by ∼100 nm over a 10 μm path. Since this broadening scales inversely with beam energy, MeV-STEM offers a promising route for high-resolution imaging in thick materials. To quantitatively assess this effect, we performed high-precision measurements at UCLA’s PEGASUS beamline, characterizing beam divergence and intensity profiles for 3–8 MeV electrons transmitted through a wedged-silicon sample of varying thickness. Our results reconcile discrepancies among analytical models and validate Monte Carlo simulations. Here, we find that increasing beam energy from 3.0 to 5.8 MeV reduces angular broadening by a factor of 2.6, with diminishing returns observed at 7.6 MeV. These findings provide a quantitative framework for optimizing MeV-STEM parameters in high-resolution imaging of thick biological and microelectronic specimens, and for guiding beam energy selection in other advanced imaging modes beyond STEM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers

Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.

Tkachenko, Nikolay V. [University of California, B↗
DOE OSTI2024

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

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

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗
DOE OSTI2025

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Characterizing the Excited States and Electronic Absorption Spectra of Small Alkylperoxy (RO 2 •) and Hydroperoxy (•QOOH) Radicals

Organic peroxy (ROO•) and hydroperoxy (•QOOH) radicals are key reactive intermediates that are formed via the oxidation of volatile organic compounds during combustion or in the Earth’s atmosphere. Their primary fate is continued unimolecular decay or bimolecular chemistry, the relative branching for which is heavily structure- and temperature-dependent. This article outlines a combined single- and multi-reference quantum chemical study to characterize the near-UV accessible electronically excited states of the prototypical ROO• and •QOOH intermediates, tert-butyl peroxy and hydroperoxy-tert-butyl radicals—the ground-state chemistries of which have been well studied both experimentally and computationally. Additionally, we simulate the electronic absorption profiles of these ROO• and •QOOH intermediates with a variety of multi- and single-reference methods. The results show an interesting conformer dependence on the electronically excited-state character and electronic absorption maxima of •QOOH. The results show promise for electronic absorption spectroscopy to be used as a selected probe for determining •QOOH conformers. Additionally, electronic absorption may contribute to the daytime removal of long-lived •QOOH intermediates formed in the troposphere. We expect that our studies will motivate experiments on the electronic absorption spectra of experimentally achievable ROO• and •QOOH.

excited states↗
DOE OSTI2025

Enhanced quantum radiation with flying-focus laser pulses

The emission of a photon by an electron in an intense laser field is one of the most fundamental processes in electrodynamics and underlies the many applications that utilize high-energy photon beams. This process is typically studied for electrons colliding head-on with a stationary-focus laser pulse. Here, we show that the energy lost by electrons in the quantum regime and the yield of emitted photons can be substantially increased by replacing a stationary-focus pulse with an equal-energy flying-focus pulse whose focus co-propagates with the electrons. Furthermore, these advantages of the flying focus result from the energy loss and the photon yield scaling more favorably with the interaction time than the laser intensity in the quantum regime, with the latter also holding in the classical regime. Monte Carlo simulations of electrons colliding with equal-energy stationary and flying-focus laser pulses demonstrate these advantages.

Compton scattering↗
DOE OSTI2026

Band-gap reduction and band alignments of dilute bismide III–V alloys

Adding a few atomic percent of Bi to III–V semiconductors leads to significant changes in their electronic structure and optical properties. Bismuth substitution on the pnictogen site leads to a large increase in spin-orbit splitting Δ SO at the top of the valence band (Γ 8⁢𝑣 −Γ 7⁢𝑣 ) and a large reduction in the band gap, creating unique opportunities in semiconductor device applications. Quantifying these changes is key to the design and simulation of electronic and optoelectronic devices. Using hybrid functional calculations, we predict the band gap of III–Vs (III = Al, Ga, In and V = As, Sb) with low concentrations of Bi (3.125% and 6.25%), the effects of adding Bi on the valence- and conduction-band edges, and the band offset between these dilute alloys and their III–V parent compounds. As expected, adding Bi raises the valence-band maximum (VBM). However, contrary to previous assumptions, the conduction-band minimum (CBM) is also significantly lowered, and both effects contribute to the sizable band-gap reduction. Changes in band gap and Δ SO are notably larger in the arsenides than in the antimonides. In conclusion, we also predict cases of band-gap inversion (Γ 6⁢𝑐 below Γ 8⁢𝑣 ), and Δ SO larger than the band gap, which are key parameters for designing topological materials and for minimizing losses due to Auger recombination in infrared lasers.

Saboor, Abdul [Univ. of Delaware, Newark, DE (Unit↗
DOE OSTI2025

Polariton spectra under the collective coupling regime. II. 2D non-linear spectra

In our previous work [Mondal et al., J. Chem. Phys. 162, 014114 (2025)], we developed several efficient computational approaches to simulate exciton–polariton dynamics described by the Holstein–Tavis–Cummings (HTC) Hamiltonian under the collective coupling regime. Here, we incorporated these strategies into the previously developed Lindblad-partially linearized density matrix (⁠$\mathscr{L}$-PLDM) approach for simulating 2D electronic spectroscopy (2DES) of exciton–polariton under the collective coupling regime. In particular, we apply the efficient quantum dynamics propagation scheme developed in Paper I to both the forward and the backward propagations in the PLDM and develop an efficient importance sampling scheme and graphics processing unit vectorization scheme that allow us to reduce the computational costs from $\mathscr{O}$($\mathscr{K}$ 2 )$\mathscr{O}$(T 3 ) to $\mathscr{O}$($\mathscr{K}$)$\mathscr{O}$(T 0 ) for the 2DES simulation, where $\mathscr{K}$ is the number of states and T is the number of time steps of propagation. As a result, we further simulated the 2DES for an HTC Hamiltonian under the collective coupling regime and analyzed the signal from both rephasing and non-rephasing contributions of the ground state bleaching, excited state emission, and stimulated emission pathways.

2D non-linear spectra↗
DOE OSTI2026

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗
DOE OSTI2025

Chiral population analysis: a real space visualization of X-ray circular dichroism

The microscopic understanding of probing and controlling molecular chirality is of considerable interest. Numerous spectroscopic techniques are capable of monitoring molecular asymmetry and its consequences, ranging from the infrared to the X-ray regime. Resonant X-rays have long been used to investigate local atomic sites within molecules thanks to the localized nature of core electronic transitions. These techniques can be used to determine the extent to which chirality is a local versus a delocalized property. However, how to systematically partition dichroic contributions from the point of view of electronic structure simulations remains an open question. Here, we introduce the concept of chiral population analysis that connects chirality to the atomic orbital picture. In analogy with Mulliken population analysis, which assigns charges to atomic orbitals, chiral populations allow the dichroic response to be distributed among the participating atomic orbitals. This decomposition can be further visualized in real space by representing it in terms of isosurface plots, providing an intuitive way to connect the dichroic response to its origins. Thus chiral population analysis can be particularly useful to assess the extent to which a given electronic transition is sensitive to chirality as a local or global feature of the molecular geometry.

36 MATERIALS SCIENCE↗
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↗
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↗
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↗