Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonthermal particles”

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 OSTI2024

Plasma Dynamics and Nonthermal Particle Acceleration in 3D Nonrelativistic Magnetic Reconnection

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β ($\ll$1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission

Magnetic reconnection—a fundamental plasma physics process, where magnetic field lines of opposite polarity annihilate—is invoked in astrophysical plasmas as a powerful mechanism of nonthermal particle acceleration, able to explain fast-evolving, bright high-energy flares. Near black holes and neutron stars, reconnection occurs in the relativistic regime, in which the mean magnetic energy per particle exceeds the rest mass energy. This review reports recent advances in our understanding of the kinetic physics of relativistic reconnection (RR): ▪ Kinetic simulations have elucidated the physics of plasma heating and nonthermal particle acceleration in RR. ▪ The physics of radiative RR, with its self-consistent interplay between photons and reconnection-accelerated particles—a peculiarity of luminous, high-energy astrophysical sources—is the new frontier of research. ▪ RR plays a key role in global models of high-energy sources, in terms of both global-scale layers and reconnection sites generated as a by-product of local magnetohydrodynamic instabilities. We summarize themes of active investigation and future directions, emphasizing the role of upcoming observational capabilities, laboratory experiments, and new computational tools.

acceleration of particles↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interplay of Three-Dimensional Instabilities and Magnetic Reconnection in the Explosive Onset of Magnetospheric Substorms

Magnetospheric substorms are preceded by a slow growth phase of magnetic flux loading and current sheet thinning in the tail. Extensive data sets have provided evidence of the triggering of instabilities at substorm onset, including magnetic reconnection and ballooning instabilities. Using an exact kinetic magnetotail equilibrium we present particle-in-cell simulations which capture the explosive nature of substorms through a disruption of the dipolarization front by the ballooning instability. We use self-consistent particle tracking to determine the nonthermal particle acceleration mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Particle Acceleration in Magnetic Reconnection with Ad Hoc Pitch-angle Scattering

Particle acceleration during magnetic reconnection is a long-standing topic in space, solar, and astrophysical plasmas. Recent 3D particle-in-cell simulations of magnetic reconnection show that particles can leave flux ropes due to 3D field-line chaos, allowing particles to access additional acceleration sites, gain more energy through Fermi acceleration, and develop a power-law energy distribution. This 3D effect does not exist in traditional 2D simulations, where particles are artificially confined to magnetic islands due to their restricted motions across field lines. Full 3D simulations, however, are prohibitively expensive for most studies. Here, we attempt to reproduce 3D results in 2D simulations by introducing ad hoc pitch-angle scattering to a small fraction of the particles. We show that scattered particles are able to transport out of 2D islands and achieve more efficient Fermi acceleration, leading to a significant increase of energetic particle flux. We also study how the scattering frequency influences the nonthermal particle spectra. This study helps achieve a complete picture of particle acceleration in magnetic reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Particle Acceleration in Relativistic Alfvénic Turbulence

Abstract Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0 . If δ B 0 ≪ B 0 , the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0 , interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2024

Magnetic Reconnection and Associated Particle Acceleration in High-Energy Astrophysics

Abstract Magnetic reconnection occurs ubiquitously in the universe and is often invoked to explain fast energy release and particle acceleration in high-energy astrophysics. The study of relativistic magnetic reconnection in the magnetically dominated regime has surged over the past two decades, revealing the physics of fast magnetic reconnection and nonthermal particle acceleration. Here we review these recent progresses, including the magnetohydrodynamic and collisionless reconnection dynamics as well as particle energization. The insights in astrophysical reconnection strongly connect to the development of magnetic reconnection in other areas, and further communication is greatly desired. We also provide a summary and discussion of key physics processes and frontier problems, toward a better understanding of the roles of magnetic reconnection in high-energy astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Particle Acceleration and Pitch-angle Evolution in Relativistic Turbulence

Synchrotron radiation detected from relativistic astrophysical objects such as pulsar-wind nebulae and jets from active galactic nuclei depends on the magnetic fields and the distribution functions of energetic electrons in these systems. Relativistic magnetically dominated turbulence has been recognized as an efficient mechanism for structure formation and nonthermal particle acceleration in these environments. Recent numerical simulations of relativistic turbulence have provided insights into the energy distribution functions of accelerated electrons. Much less is currently understood about their pitch angle distributions, which are crucial for accurately interpreting the spectra of synchrotron radiation. We perform a detailed case study of the pitch angle distributions formed during the process of turbulent acceleration for B 0 /δB 0 = 10 and $\tilde{σ}$ 0 ~40, where B0 is the uniform component of the magnetic field, δB0 is the fluctuating component, and $\tilde{σ}$ 0 is the plasma magnetization based on the magnetic fluctuations. We find that even minimal numerical noise can cause substantial pitch angle scattering, but we demonstrate techniques for overcoming the numerical challenges associated with the evolution of very small pitch angles. Our numerical results are consistent with the phenomenological considerations found in C. Vega et al. (2024a, 2025).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

First-Principles Theory of the Relativistic Magnetic Reconnection Rate in Astrophysical Pair Plasmas

In this paper, we develop a first-principles model for the relativistic magnetic reconnection rate in strongly magnetized pair plasmas. By considering the energy budget and required current density near the x-line, we analytically show that in the magnetically dominated relativistic regime, the x-line thermal pressure is significantly lower than the upstream magnetic pressure due to the extreme energy needed to sustain the current density, consistent with kinetic simulations. This causes the upstream magnetic field lines to collapse in, producing the open outflow geometry which enables fast reconnection. The result is important for understanding a wide range of extreme astrophysical environments, where fast reconnection has been evoked to explain observations such as transient flares and nonthermal particle signatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

First-principles-integrated Study of Blazar Synchrotron Radiation and Polarization Signatures from Magnetic Turbulence

Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principles-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multiwavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability timescale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle swings with arbitrary amplitudes and duration, in either direction, that are not associated with changes in flux or PD. Surprisingly, these swings, which are stochastic in nature, can appear either bumpy or smooth, although large-amplitude swings (>180°) are very rare, as expected. Our radiation and polarization signatures from first-principles-integrated simulations are consistent with the TEMZ model, except that in the latter, there is a weak correlation, with zero lag, between flux and degree of polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Fully Kinetic Simulations of Proton-beam-driven Instabilities from Parker Solar Probe Observations

The expanding solar wind plasma ubiquitously exhibits anisotropic nonthermal particle velocity distributions. Typically, proton velocity distribution functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by the Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam population that drift with respect to each other. This configuration triggers a proton beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing perpendicular heating at the expense of the parallel temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

IceCube's convex all-sky neutrino spectrum consistent with the magnetically powered corona scenario for active galactic nuclei

High-energy multimessenger background analyses over the past decade have provided evidence for a population of hidden neutrino sources that are opaque to GeV-TeV gamma rays, a picture bolstered by recent observations of the nearby active galaxy NGC 1068. The coronal regions in the hearts of active galactic nuclei (AGNs) have been proposed as the most promising sites for such hidden nonthermal particle production, and NGC 1068 is expected to be the most neutrino-active galaxy for IceCube. We demonstrate that the latest all-sky neutrino spectrum, exhibiting a spectral bend around 3-30 TeV, is consistent with predictions of the magnetically powered corona scenario, and the models for the all-sky neutrino flux can simultaneously explain the multimessenger data from NGC 1068 within observational and modeling uncertainties. We further show, in a largely model-independent way, that the contribution from NGC 1068-like sources does not overshoot the observed medium-energy neutrino flux. Finally, we highlight the key role of the Eddington ratio, which can drive substantial variations in the predicted neutrino fluxes of nearby AGNs, and we encourage systematic multimessenger searches for the neutrino-brightest AGNs.

Murase, Kohta [Penn State U.; Penn State U., Unive↗
DOE OSTI2025

Cold freeze out of superheavy dark matter and Hubble tension

We present a unified framework, the "X miracle", in which dark matter consists of superheavy, nonthermal X particles whose relic abundance is determined not by the conventional weak-scale, semi-relativistic ("hot") freeze-out of WIMPs, but by annihilation or decay occurring within the smallest and earliest gravitationally bound objects. Unlike thermal WIMPs, which decouple at velocities of order 0.3c with relic abundance ρ∞ set by weak-scale interactions, X particles are produced nonthermally with an initial overabun dance ρ ini >> ρ ∞ . They become nonrelativistic extremely early, redshift to ultra-cold velocities, allowing collapse into compact bound structures characterized by a novel quantum gravitational scale, r X = 4$\hbar$ 2 $/Gm^3_X$ = 10 −13 m $\hbar$$/m_Xc$, much larger than the Compton wavelength. The framework predicts a particle mass of 10 12 GeV and an enhanced cross section of 10 −21 m 3 /s. Overlapping particle wavefunctions in these compact structures drive annihilation or decay into additional radiation, leading to a "cold" freeze-out that converts most of ρ ini into radiation while leaving a relic density ρ ∞ . Solutions to the Boltzmann equation indicate that an extreme ("big") depletion, with only one particle in a billion surviving, yields an additional radiation contribution $ΔN_{eff}$ ≈ 0.4, which could help alleviate the Hubble tension. For particles of 10 12 GeV, the scenario predicts a dark coupling constant α X = 0.09 that is responsible for an instanton-induced decay process, consistent with current UHECR bounds. Early collapse at 10 −6 s may release binding energy as high-frequency (100kHz) gravitational waves or ultralight GUT-scale axions (10 −9 eV). Superheavy sterile neutrinos provide a natural particle realization, linking dark matter to neutrino mass and baryogenesis. If gravitationally produced, this framework favors high-scale inflation and effi cient reheating. The "X miracle" thus demonstrates that dark matter need not be weak-scale: gravitational dynamics can control freeze-out and evolution, producing multi-messenger observational signatures in UHECRs, axions, gravitational waves, and small-scale structures.

Xu, Zhijie Jay [Pacific Northwest National Laborat↗
DOE OSTI2024

Proton Acceleration in Low- β Magnetic Reconnection with Energetic Particle Feedback

Magnetic reconnection regions in space and astrophysics are known as active particle acceleration sites. There is ample evidence showing that energetic particles can take a substantial amount of converted energy during magnetic reconnection. However, there has been a lack of studies understanding the backreaction of energetic particles at magnetohydrodynamical scales in magnetic reconnection. To address this, we have developed a new computational method to explore the feedback by nonthermal energetic particles. This approach considers the backreaction from these energetic particles by incorporating their pressure into magnetohydrodynamics (MHD) equations. The pressure of the energetic particles is evaluated from their distribution evolved through Parker's transport equation, solved using stochastic differential equations (SDEs), so we coin the name MHD-SDE. Applying this method to low-β magnetic reconnection simulations, we find that reconnection is capable of accelerating a large fraction of energetic particles that contain a substantial amount of energy. When the feedback from these particles is included, their pressure suppresses the compression structures generated by magnetic reconnection, thereby mediating particle energization. Consequently, the feedback from energetic particles results in a steeper power-law energy spectrum. These findings suggest that feedback from nonthermal energetic particles plays a crucial role in magnetic reconnection and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Collisionless relativistic magnetic reconnection driven by electron vortices in laser-plasma interaction

Magnetic reconnection (MR) is a fundamental process in space and laboratory plasmas. The appearance of high power lasers opens a new way to investigate MR under the relativistic condition. In this paper, relativistic collisionless MR driven by two ultra-intense lasers and a pair of asymmetric targets is studied numerically via the kinetic simulations. The static magnetic fields produced by the electron vortex structures with opposite magnetic polarities approach each other driven by the magnetic pressure and the density gradient. The antiparallel magnetic fields annihilate accompanied with the topological variation and the corresponding magnetic field energy is being dissipated to the kinetic energy of the nonthermal charged particles. Besides the outflows along the current sheet, a fast particle bunch is accelerated perpendicularly contributed by the displacement current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗