- Award ID(s):
- 2109154
- Publication Date:
- NSF-PAR ID:
- 10351959
- Journal Name:
- The Astrophysical Journal
- Volume:
- 933
- Issue:
- 1
- Page Range or eLocation-ID:
- 73
- ISSN:
- 0004-637X
- Sponsoring Org:
- National Science Foundation
More Like this
-
We investigate particle acceleration in an MHD-scale system of multiple current sheets by performing 2D and 3D MHD simulations combined with a test particle simulation. The system is unstable for the tearing-mode instability, and magnetic islands are produced by magnetic reconnection. Due to the interaction of magnetic islands, the system relaxes to a turbulent state. The 2D (3D) case both yield −5/3 (− 11/3 and −7/3) power-law spectra for magnetic and velocity fluctuations. Particles are efficiently energized by the generated turbulence, and form a power-law tail with an index of −2.2 and −4.2 in the energy distribution function for the 2D and 3D case, respectively. We find more energetic particles outside magnetic islands than inside. We observe super-diffusion in the 2D (∼ t 2.27 ) and 3D (∼ t 1.2 ) case in the energy space of energetic particles.
-
Abstract Magnetic reconnection is invoked as one of the primary mechanisms to produce energetic particles. We employ large-scale 3D particle-in-cell simulations of reconnection in magnetically dominated ( σ = 10) pair plasmas to study the energization physics of high-energy particles. We identify an acceleration mechanism that only operates in 3D. For weak guide fields, 3D plasmoids/flux ropes extend along the z -direction of the electric current for a length comparable to their cross-sectional radius. Unlike in 2D simulations, where particles are buried in plasmoids, in 3D we find that a fraction of particles with γ ≳ 3 σ can escape from plasmoids by moving along z , and so they can experience the large-scale fields in the upstream region. These “free” particles preferentially move in z along Speiser-like orbits sampling both sides of the layer and are accelerated linearly in time—their Lorentz factor scales as γ ∝ t , in contrast to γ ∝ t in 2D. The energy gain rate approaches ∼ eE rec c , where E rec ≃ 0.1 B 0 is the reconnection electric field and B 0 the upstream magnetic field. The spectrum of free particles is hard, dN free / d γ ∝ γmore »
-
Abstract We perform a 2.5-dimensional particle-in-cell simulation of a quasi-parallel shock, using parameters for the Earth’s bow shock, to examine electron acceleration and heating due to magnetic reconnection. The shock transition region evolves from the ion-coupled reconnection dominant stage to the electron-only reconnection dominant stage, as time elapses. The electron temperature enhances locally in each reconnection site, and ion-scale magnetic islands generated by ion-coupled reconnection show the most significant enhancement of the electron temperature. The electron energy spectrum shows a power law, with a power-law index around 6. We perform electron trajectory tracing to understand how they are energized. Some electrons interact with multiple electron-only reconnection sties, and Fermi acceleration occurs during multiple reflections. Electrons trapped in ion-scale magnetic islands can be accelerated in another mechanism. Islands move in the shock transition region, and electrons can obtain larger energy from the in-plane electric field than the electric potential in those islands. These newly found energization mechanisms in magnetic islands in the shock can accelerate electrons to energies larger than the achievable energies by the conventional energization due to the parallel electric field and shock drift acceleration. This study based on the selected particle analysis indicates that the maximum energy inmore »
-
Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated
γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LATγ -ray bands as well as closely correlated optical andγ -ray flares. The swings in optical polarization angle are also accompanied byγ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fastγ -ray flares or orphanγ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very highmore » -
ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstreammore »