SEARCH · Search NASA

Results for “RELATIVISTIC PARTICLE”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Bosenovae with quadratically-coupled scalars in quantum sensing experiments

Abstract Ultralight dark matter (ULDM) particles of massm ϕ ≲ 1 eV can form boson stars in DM halos. Collapse of boson stars leads to explosive bosenova emission of copious relativistic ULDM particles. In this work, we analyze the sensitivity of terrestrial and space-based experiments to detect such relativistic scalar ULDM particles interacting through quadratic couplings with Standard Model constituents, including electrons, photons, and gluons. We highlight key differences with searches for linear ULDM couplings. Screening of ULDM with quadratic couplings near the surface of the Earth can significantly impact observations in terrestrial experiments, motivating future space-based experiments. We demonstrate excellent ULDM discovery prospects, especially for quantum sensors, which can probe quadratic couplings orders below existing constraints by detecting bosenova events in the ULDM mass range 10 −23 eV ≲m ϕ ≲ 10 −5 eV. We also report updated constraints on quadratic couplings of ULDM in case it comprises cold DM.

Physics↗

DOE OSTI2023

Verification and benchmarking relativistic electron beam transport through a background gas

It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. Additionally, we also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Particle Acceleration in Collisionless Magnetically Arrested Disks

We present the first collisionless realization of two-dimensional axisymmetric black hole accretion consistent with a persistent magnetically arrested disk state. The accretion flow, consisting of an ion-electron disk plasma combined with magnetospheric pair creation effects, is simulated using first-principles general-relativistic particle-in-cell methods. The simulation is evolved over significant dynamical timescales during which a quasisteady accretion state is reached with several magnetic flux eruption cycles. We include a realistic treatment of inverse Compton scattering and pair production, which allows for studying the interaction between the collisionless accretion flow and pair-loaded jet. Our findings indicate that magnetic flux eruptions associated with equatorial magnetic reconnection within the black hole magnetosphere and the formation of spark gaps are locations of maximal particle acceleration. Flux eruptions, starting near the central black hole, can trigger Kelvin-Helmholtz-like vortices at the jet-disk interface that facilitate efficient mixing between disk and jet plasma in this region. Transient periods of increased pair production following magnetic flux eruptions and reconnection events are responsible for most of the highly accelerated particles.

Accretion disk & black-hole plasma↗

DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗

DOE OSTI2020

Cosmological Constant in Coherent Quantum Gravity

It is argued that quantum states of geometry, like those of particles, should be coherent on light cones of any size. An exact classical solution, the gravitational shock wave of a relativistic point particle, is used to estimate gravitational drag from coherent energy flows, and the expected gravitational effect of virtual transverse vacuum energy fluctuations on surfaces of causal diamonds. It is proposed that the appropriately spacetime-averaged gravitational effect of the Standard Model vacuum state leads to the observed small nonzero value of the cosmological constant, dominated by gravitational drag of virtual gluonic strings at the strong interaction scale.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Electron bunch dynamics and emission in particle-in-cell simulations of relativistic laser–solid interactions: On density artifacts, collisions, and numerical dispersion

Sub-optical-cycle dynamics of dense electron bunches in relativistic-intensity laser–solid interactions lead to the emission of high-order harmonics and attosecond light pulses. The capacity of particle-in-cell simulations to accurately model these dynamics is essential for the prediction of emission properties because the attosecond pulse intensity depends on the electron density distribution at the time of emission and on the temporal distribution of individual electron Lorentz-factors in an emitting electron bunch. Here, we show that in one-dimensional collisionless simulations, the peak density of the emitting electron bunch increases with the increase in the spatial resolution of the simulation grid. When collisions are added to the model, the peak electron density becomes independent of the spatial resolution. Collisions are shown to increase the spread of the peaks of Lorentz-factors of emitting electrons in time, especially in the regimes far from optimum generation conditions, thus leading to lower intensities of attosecond pulses as compared to those obtained in collisionless simulations.

Physics↗

DOE OSTI2022

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Magnetic Field Amplification by a Plasma Cavitation Instability in Relativistic Shock Precursors

Abstract Plasma streaming instabilities play an important role in magnetic field amplification and particle acceleration in relativistic shocks and their environments. However, in the far shock precursor region where accelerated particles constitute a highly relativistic and dilute beam, streaming instabilities typically become inefficient and operate at very small scales when compared to the gyroradii of the beam particles. We report on a plasma cavitation instability that is driven by dilute relativistic beams and can increase both the magnetic field strength and coherence scale by orders of magnitude to reach near-equipartition values with the beam energy density. This instability grows after the development of the Weibel instability and is associated with the asymmetric response of background leptons and ions to the beam current. The resulting net inductive electric field drives a strong energy asymmetry between positively and negatively charged beam species. Large-scale particle-in-cell simulations are used to verify analytical predictions for the growth and saturation level of the instability and indicate that it is robust over a wide range of conditions, including those associated with pair-loaded plasmas. These results can have important implications for the magnetization and structure of shocks in gamma-ray bursts, and more generally for magnetic field amplification and asymmetric scattering of relativistic charged particles in plasma astrophysical environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Effect of magnetic field inclination on black hole jet power and particle acceleration

Rotating black holes are known to launch relativistic jets and accelerate particles if they accrete a magnetized plasma. It remains unclear, however, how the global magnetic field orientation affects the jet powering efficiency. We propose the first kinetic study of a collisionless plasma around a Kerr black hole embedded in a magnetic field that is inclined with respect to the black hole spin axis. Using three-dimensional general relativistic particle-in-cell simulations, we show that while oblique magnetic field configurations significantly reduce the jet power, particle acceleration still remains highly efficient. This suggests that black holes producing a weak jet might still be bright sources of nonthermal radiation and cosmic rays.

acceleration of particles↗

DOE OSTI2026

The Effect of Inverse Compton Losses on Particle Acceleration in Three-dimensional Relativistic Reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived “free phase” where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a “trapped phase,” where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection, including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor γrad exceeds the magnetization σ. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is dN free /dγ ∝ γ −1 , as in the uncooled case. The spectrum of trapped electrons above the cooling break γcool (in the range γ cool < γ < γ rad ) is dN/dγ ∝ γ −3 , steeper than the scaling dN/dγ ∝ γ −2 of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by Zhang et al. for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived "free phase" where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a "trapped phase", where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor $γ_{\rm rad}$ exceeds the magnetization $σ$. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is $dN_{\rm free}/dγ\propto γ^{-1}$, as in the uncooled case. The spectrum of trapped electrons above the cooling break $γ_{\rm cool}$ (in the range $γ_{\rm cool}<γ<γ_{\rm rad}$) is $dN/dγ\propto γ^{-3}$, steeper than the scaling $dN/dγ\propto γ^{-2}$ of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by arXiv:2302.12269 for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure, that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

FOS: Physical sciences↗

DOE OSTI2025

Particle Beam Acceleration Using 3 Petawatt Laser Pulses

The Zettawatt-Equivalent Ultrashort pulse laser System (ZEUS) is presently operational at the Gerard Mourou Center for Ultrafast Optical Science (CUOS) at the University of Michigan. ZEUS is a significant upgrade of the previous high power laser systems at CUOS and consists of two beamlines thatoperate in perfect synchronization. The 500 TW beamline became operational in 2023, 2 PW operation started in 2025 and full 3 PW power levels will be available in 2026. It is presently the highest power laser system in the US. In this grant the high field science group at CUOS has leveraged this unique high power laser facility to investigate laser wake field acceleration (LWFA) in ultra-high power laser plasma interactions and have shown how this can scale for future electron–positron colliders at high energy. The dual beam experimental configuration enables flexibility for many frontier experiments in laser-driven acceleration research, in particular, enabling extended channelling/acceleration experiments, positron generation/acceleration experiments and proof-of-principle transverse pumping “dephasingless” electron acceleration experiment and theory. LWFA may be able to miniaturize particle accelerators for high energy physics and also enable new sources of ultrafast, extreme brightness and precise x-rays for a wide variety of applications. In laser wake field acceleration, an electron bunch “surfs” on the electron plasma wave (the “wake field”) generated by the ponderomotive force of an intense laser. The plasma wave has a strong longitudinal electric field that stays in phase with the relativistic driver. A relativistic charged particle may, therefore, remain in phase with the accelerating field over long distances and gain ultra-relativistic energies. The accelerating electric field strength that the plasma wave can support can be many orders of magnitude higher than that of conventional accelerators, which makes laser wakefield acceleration an exciting prospect as an advanced accelerator concept. In this research project we have investigated the scaling of this mechanism to laser powers of 2 PW and have measured the x-ray emission and radio frequency emission resulting from the acceleration process. We have also performed theoretical investigation of mechanisms to scale laser driven accelerators to much higher energy using dephasingless acceleration processes.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Relativistic plasma physics in supercritical fields

Since the invention of chirped pulse amplification, which was recognized by a Nobel Prize in physics in 2018, there has been a continuing increase in available laser intensity. Combined with advances in our understanding of the kinetics of relativistic plasma, studies of laser-plasma interactions are entering a new regime where the physics of relativistic plasmas is strongly affected by strong-field quantum electrodynamics (QED) processes, including hard photon emission and electron-positron (e⁻-e⁺) pair production. Additionally, this coupling of quantum emission processes and relativistic collective particle dynamics can result in dramatically new plasma physics phenomena, such as the generation of dense e⁻-e⁺ pair plasma from near vacuum, complete laser energy absorption by QED processes, or the stopping of an ultra-relativistic electron beam, which could penetrate a cm of lead, by a hair's breadth of laser light. In addition to being of fundamental interest, it is crucial to study this new regime to understand the next generation of ultra-high intensity laser-matter experiments and their resulting applications, such as high energy ion, electron, positron, and photon sources for fundamental physics studies, medical radiotherapy, and next generation radiography for homeland security and industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗