Search NASA⌕ Search

SEARCH · Search NASA

Results for “RELATIVISTIC PLASMA”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Tearing-mediated Alfvénic Turbulence in a Relativistic Plasma

In astrophysical relativistic plasmas, Alfvénic turbulence exists across a wide range of scales. Similar to a nonrelativistic case, this turbulence can be influenced by the tearing instability. We argue that in an ultrarelativistic pair plasma the tearing effects become significant at the critical scale given by ${a}_{t}\sim {d}_{\mathrm{rel}}{\left({L}_{\perp }/{d}_{\mathrm{rel}}\right)}^{n/(4+5n)}$, where L ⊥ represents the outer scale of the turbulence and d rel the relativistic electron inertial scale. Here n = 1 corresponds to the Harris magnetic configurations, while n = 2 applies to localized (non-Harris) magnetic profiles of turbulence-generated current sheets. Below the critical scale, the energy spectrum of magnetic fluctuations steepens to W(k)dk ∝ k -(7n+2)/3n dk. However, this happens only when $\tilde{\sigma }\ll {\left({L}_{\perp }/{d}_{\mathrm{rel}}\right)}^{(2+4n)/(4+5n)}$, where $\tilde{\sigma }$ is the magnetization parameter based on the magnetic fluctuations. If the opposite inequality holds, the Alfvénic cascade is instead affected by charge starvation at scales below ${a}_{\mathrm{ch}}\sim {d}_{\mathrm{rel}}{\left({d}_{\mathrm{rel}}/{L}_{\perp }\right)}^{1/3}{\tilde{\sigma }}^{2/3}$, and the energy spectrum becomes steeper than W(k)dk ∝ k -3 dk. When d rel is replaced by its nonrelativistic equivalent, our results agree with the nonrelativistic analysis of N. F. Loureiro & S. Boldyrev (2018). In the Harris case n = 1, they agree with the recent relativistic studies by I. Demidov & Y. Lyubarsky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Synchrotron-driven Instabilities in Relativistic Plasmas of Arbitrary Opacity

Recent work has shown that synchrotron emission from relativistic plasmas leads the electron distribution to form an anisotropic ring in momentum space, which can be unstable to both kinetic and hydrodynamic instabilities. Fundamental to these works was the assumption that the plasma was optically thin, allowing all emitted radiation to escape. Here, we examine the behavior of these instabilities as the plasma becomes more optically thick. To do this, we extend a recently developed Fokker–Planck operator for synchrotron emission and absorption in mildly relativistic plasmas to ultrarelativistic plasmas. For a given set of plasma parameters, photons emitted by higher-energy electrons tend to be higher frequency, and thus more easily escape the plasma. As a result, the ratio of the photon emission rate (radiative drag) to absorption rate (radiative diffusion) for a given electron is extremely energy dependent. Given this behavior, we determine the critical parameters that control the opacity, and show how the plasma gradually transitions to become more isotropic and stable at higher opacity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Efficiency-optimized relativistic plasma harmonics for extreme fields

Bright harmonic radiation from relativistically oscillating laser plasmas offers a direct route for generating extreme electromagnetic fields. Theory predicts that under optimized conditions, the plasma medium can support strong spatiotemporal compression of laser energy in a coherent harmonic focus (CHF), delivering intensity boosts many orders of magnitude greater than the incident driving laser pulse. Although diffraction-limited performance (spatial compression) and attosecond phase locking (temporal compression) have been demonstrated experimentally, efficient coupling of relativistically intense laser pulse energy into the emitted harmonic cone has not been realized so far. Here we demonstrate that this highly nonlinear interaction can be tailored to deliver the maximum conversion efficiencies predicted from simulations. By fine-tuning the temporal profile of the driving laser on sub-picosecond (<10 −12 s) timescales, energies >9 mJ between the 12th and 47th harmonics are observed. These results are in agreement with the theoretically expected efficiency dependence on harmonic order, verifying that optimal conditions have been achieved in the generation process. This is the important final element required to achieve the expected intensity boosts from a CHF in experiments. Although obtaining spatiotemporal compression and optimal efficiency simultaneously remains challenging, the path to realizing extreme optical field strengths approaching the critical field of quantum electrodynamics (the Schwinger limit at >10 16 V cm −1 or >10 29 W cm −2 ) is now open, permitting all-optical studies of the quantum vacuum and new frontiers for intense attosecond science.

High-harmonic generation↗
DOE OSTI2024

Electron tail suppression and effective collisionality due to synchrotron emission and absorption in mildly relativistic plasmas

Synchrotron radiation losses are a significant cause of concern for high-temperature aneutronic fusion reactions such as proton–Boron 11. The fact that radiation losses occur primarily in the high-energy tail, where the radiation itself has a substantial impact on the electron distribution, necessitates a self-consistent approach to modeling the diffusion and drag induced by synchrotron absorption and emission. Furthermore, an accurate model must account for the fact that the radiation emission spectrum is momentum-dependent, and the plasma opacity is frequency-dependent. Here, we present a simple Fokker–Planck operator, built on a newly solved-for blackbody synchrotron diffusion operator, which captures all relevant features of the synchrotron radiation. Focusing on magnetic mirror fusion plasmas, we show that significant suppression of the electron distribution occurs for relativistic values of the perpendicular electron momentum, which therefore emit much less radiation than predicted under the assumption of a Maxwell–Jüttner distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Relativistic plasma effects and ion acceleration using Scarlet and ELI-NP

High energy density science (HEDS) explores the nature of matter under extreme conditions of temperature and pressure. It has wide ranging fundamental importance from understanding the physical structure of our planet to extreme phenomena in the cosmos. It also has many applications such as facilitating imaging with ions, neutrons, x-rays, gamma rays and more with new applications being developed, including materials processing and medical therapies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Physics of Ultra-High Energy Density Relativistic Plasmas from Ordered Nanostructures

This award funded a four-year program at Colorado State University on the physics of ultra-high energy density plasmas generated by focusing femtosecond, >1021 Wcm-2 laser pulses onto ordered nanowire arrays. We used the ALEPH petawatt-class Ti:Sa laser with its high contrast 400 nm second harmonic beamline to irradiate arrays of Ni and TiO2 nanowires made by template assisted electrodeposition and atomic layer deposition over a wide range of average densities. We met two of the three original project goals: radial ion acceleration via doppler shifted K shell emission, and time and spectrally resolved K shell measurements. The main experimental advance was integrating the Lawrence Livermore National Laboratory (LLNL) TREX ultrafast x-ray streak camera with a suite of spherically bent quartz Bragg crystal spectrometers built using Princeton Plasma Physics Laboratory (PPPL) methodology. This is the first TREX-on-bent-quartz setup fielded at a femtosecond, >1021 Wcm-2 facility. It is now a permanent diagnostic at ALEPH and is offered to the LaserNetUS user community. The program produced three peer reviewed journal articles that acknowledge this award. The platform will continue at the upgraded CSU ATLAS facility, where the ALEPH laser is being upgraded to 2 PW.

Hollinger, Reed [Colorado State University]↗
DOE OSTI2024

Laboratory realization of relativistic pair-plasma beams

Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black-hole and neutron-star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with electron-positron pairs. Their role in the dynamics of such environments is in many cases believed to be fundamental, but their behavior differs significantly from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies, which are rather limited. We present the first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN’s Super Proton Synchrotron (SPS) accelerator. Monte Carlo simulations agree well with the experimental data and show that the characteristic scales necessary for collective plasma behavior, such as the Debye length and the collisionless skin depth, are exceeded by the measured size of the produced pair beams. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Thermal Weibel instability induced magnetic fields co-exist with linear wakes in laser-ionized plasmas

When a moderately intense, few-picosecond-long laser pulse ionizes gas to produce an underdense plasma column, a linear relativistic plasma wave or wake can be excited by the self-modulation instability that may prove useful for multi-bunch acceleration of externally injected electrons or positrons to high energies in a short distance. At the same time, due to the anisotropic temperature distributions of the ionized plasma electrons, the Weibel instability can self-generate magnetic fields throughout such a plasma on a few picoseconds timescale that can persist even longer than the lifetime of the wake. In the present paper, we first show using simulations that both these effects do indeed co-exist in space and time in the plasma. Using our simulations, we make preliminary estimates of the contribution to the transverse emittance growth of an externally injected beam due to the Weibel magnetic fields in a few-millimeter-long plasma. We then present the results of an experiment that has allowed us to measure the spatiotemporal evolution of the magnetic fields using an ultrashort relativistic electron probe beam. Both the topology and the lifetime of the Weibel instability induced magnetic fields in the experiment are in reasonable agreement with the fields induced by the Weibel instability in the simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Relativistic Laser Plasma Interactions At The Highest Intensities

High energy density science (HEDS) explores the nature of matter under extreme conditions of temperature and pressure. It is of fundamental importance and has many applications such as facilitating imaging with ions, neutrons, x-rays, and gamma rays with new applications being developed, including materials processing and medical therapies. In this project, we used high power, ultrashort pulse lasers to reach HEDS conditions. We have shown that low-cost, liquid crystal film based, double plasma mirror systems can be used to greatly improve laser pulse contrast while still maintaining high power and excellent spatial mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The role of electric dominance for particle injection in relativistic reconnection

ABSTRACT Magnetic reconnection in relativistic plasmas – where the magnetization $\sigma \gg 1$ – is regarded as an efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. We employ two-dimensional (2D) particle-in-cell simulations of relativistic pair-plasma reconnection with vanishing guide field and outflow boundaries to quantify the impact of the energy gain occurring in regions of electric dominance ($E\gt B$) for the early stages of particle acceleration (i.e. the ‘injection’ stage). Given an injection threshold energy $\epsilon ^\ast =\sigma /4$ for the particles that eventually attain energy $\epsilon _{\rm T}$ by time T, we calculate the mean fractional contribution $\zeta (\epsilon ^\ast ,\epsilon _{\rm T})$ by $E\gt B$ fields to particle energization at the time when the threshold $\epsilon ^\ast$ is reached. We find that $\zeta$ monotonically increases with $\sigma$ and $\epsilon _{\rm T}$; for $\sigma \gtrsim 50$ and $\epsilon _{\rm T}/\sigma \gtrsim 8$, we find that $\gtrsim 80~{{\ \rm per\ cent}}$ of the energy gain obtained before reaching $\epsilon ^\ast =\sigma /4$ occurs in $E\gt B$ regions. We find that $\zeta$ is independent of simulation box size $L_x$, as long as $\epsilon _{\rm T}$ is normalized to the maximum particle energy, which scales as $\epsilon _{\rm max}\propto L_{\rm x}^{1/2}$ in 2D. The distribution of energy gains $\epsilon _{\chi }$ acquired in $E\gt B$ regions can be modelled as $\mathrm{ d}N/\mathrm{ d}\epsilon _{\chi }\propto \epsilon _{\chi }^{-0.35}\exp [-(\epsilon _{\chi }/0.06\, \sigma)^{0.5}]$. Our results help assess the role of electric dominance in relativistic reconnection with vanishing guide fields, which is realized in the magnetospheres of black holes and neutron stars.

Gupta, Sanya (ORCID:0000000151944384)↗
DOE OSTI2025

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗
DOE OSTI2024

Force-free Wave Interaction in Magnetar Magnetospheres: Computational Modeling in Axisymmetry

Abstract Crustal quakes of highly magnetized neutron stars can disrupt their magnetospheres, triggering energetic phenomena like X-ray and fast radio bursts. Understanding plasma wave dynamics in these extreme environments is vital for predicting energy transport across scales to the radiation length. This study models relativistic plasma wave interaction in magnetar magnetospheres with force-free electrodynamics simulations. For propagation along curved magnetic field lines, we observe the continuous conversion of Alfvén waves to fast magnetosonic (FMS) waves. The conversion efficiency can be up to three times higher when counter-propagating Alfvén waves interact in the equatorial region. Alfvén waves generate FMS waves of twice their frequency during their first crossing of the magnetosphere. After the initial transient burst of FMS waves, Alfvén waves convert to FMS waves periodically, generating variations on timescales of the magnetospheric Alfvén wave crossing time. This decaying FMS wave tail carries a significant portion (half) of the total energy emitted. Plastic damping of “bouncing” Alfvén waves by the magnetar crust has minimal impact on the FMS efficiency. We discuss the implications of the identified wave phenomena for magnetar observations. Outgoing FMS waves can develop electric zones, potential sources of coherent radiation. Long wavelength FMS waves could generate FRBs through reconnection beyond the light cylinder.

Astronomy & Astrophysics↗
DOE OSTI2024

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Anisotropic Particle Acceleration in Alfvénic Turbulence

Alfvénic turbulence is an effective mechanism for particle acceleration in strongly magnetized, relativistic plasma. In this study, we investigate a scenario where turbulent plasma is influenced by a strong guide magnetic field, resulting in highly anisotropic turbulent fluctuations. In such cases, the magnetic moments of particles are conserved, which means that acceleration can only occur along the direction of the magnetic field. Consistent with previous analytic studies, we find through particle-in-cell simulations of magnetically dominated pair plasma that the momenta of accelerated particles are closely aligned with the magnetic field lines. Notably, the alignment angle decreases as particle energy increases, potentially limited only by the inherent curvature and gradients of the turbulent magnetic fluctuations. This finding has significant implications for interpreting the synchrotron radiation emitted by highly accelerated particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. II. General Relativistic Module

Black hole (BH) environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae, where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general relativistic module for Entity—the first open-source, coordinate-agnostic, performance-portable, particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around BHs on any modern high-performance computing architecture (both GPUs and CPUs).

Galishnikova, Alisa [Flatiron Institute, New York,↗
DOE OSTI2025

Minimizing Selection Bias in Inclusive Jets in Heavy-Ion Collisions with Energy Correlators

The first-ever measurement of energy correlators within inclusive jets produced in heavy-ion collisions, revealed by the CMS Collaboration, shows a clear enhancement at large angles relative to the proton-proton (p-p) baseline. However, interpreting this enhancement is complicated due to selection bias from energy loss, which also distorts the energy correlator heavy-ion to p-p ratio in the hadronization region, hindering our understanding of parton/hadron dynamics in a colored medium. In this Letter, we introduce a new ratio of energy correlator observables that removes the leading effects of selection bias from the two-point energy correlator spectrum (E2C). PYTHIA and HERWIG simulations show that the impact of selection bias in the E2C is reduced by an order of magnitude, while sensitivity to any other medium modifications is retained. This quantity can be obtained directly from the experimental measurements presented by CMS, as illustrated in an accompanying note [C. Andres and J. Holguin, arXiv:2409.07526].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Precision Measurement of Net-Proton-Number Fluctuations in Au + Au Collisions at RHIC

We report precision measurements on cumulants (𝐶 𝑛 ) and factorial cumulants (𝜅 𝑛 ) of (net) proton number distributions up to fourth order in Au + Au collisions over center-of-mass energies $\sqrt{s_{NN}}$ = 7.7–27 GeV from phase II of the Beam Energy Scan program at RHIC. (Anti)protons are selected at midrapidity (|𝑦| < 0.5) within a transverse momentum range of 0.4 < 𝑝 𝑇 < 2.0 GeV/𝑐. Relative to various noncritical-point model calculations and peripheral collision 70%–80% data, the net proton 𝐶 4 /𝐶 2 measurement in 0%–5% collisions shows a minimum around 19.6 GeV for significance of deviation at ∼2–5⁢𝜎. A minimum in 𝐶 4 /𝐶 2 with respect to a noncritical baseline is expected to be a characteristic feature of the signature associated with a critical point in the QCD phase diagram. In addition, deviations from noncritical baselines around the same collision energy region are also seen in proton factorial cumulant ratios, especially in 𝜅 2 /𝜅 1 and 𝜅 3 /𝜅 1 . As a result, dynamical model calculations including a critical point are called for in order to understand these precision measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗