No abstract

In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be the primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue in studies of relativistic reconnection is the consequence of three-dimensional dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of three-dimensional, fully kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron–electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long-wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. We find that plasma dynamics in RTMR is vastly different from its 2D counterpart in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted by the secondary kink instability. This turbulent evolution leads to superdiffusive behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and the timescale for energy release can be very fast and do not depend strongly on the turbulence amplitude. The main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the nonideal electric field acceleration plays a subdominant role. We also discuss possible observational implications of three-dimensional RTMR in high-energy astrophysics.

Magnetic reconnection in the relativistic regime has been proposed as an important process for the efficient production of nonthermal particles and high-energy emission. Using fully kinetic particle-in-cell simulations, we investigate how the guide-field strength and domain size affect the characteristic spectral features and acceleration processes. We study two stages of acceleration: energization up until the injection energy γ inj and further acceleration that generates a power-law spectrum. Stronger guide fields increase the power-law index and γ inj, which suppresses acceleration efficiency. These quantities seemingly converge with increasing domain size, suggesting that our findings can be extended to large-scale systems. We find that three distinct mechanisms contribute to acceleration during injection: particle streaming along the parallel electric field, Fermi reflection, and the pickup process. The Fermi and pickup processes, related to the electric field perpendicular to the magnetic field, govern the injection for weak guide fields and larger domains. Meanwhile, parallel electric fields are important for injection in the strong guide-field regime. In the post-injection stage, we find that perpendicular electric fields dominate particle acceleration in the weak guide-field regime, whereas parallel electric fields control acceleration for strong guide fields. These findings will help explain the nonthermal acceleration and emission in high-energy astrophysics, including black hole jets and pulsar wind nebulae.

Kinetic Riemann simulations have been completed to explore particle heating during guide field reconnection in the low-β environment of the inner heliosphere and the solar corona. The reconnection exhaust is bounded by two rotational discontinuities (RD) and two slow shocks (SS) form within the exhaust as in magnetohydrodynamic (MHD) models. At the RDs, ions are accelerated by the magnetic field tension to drive the reconnection outflow as well as flows in the out-of-plane direction. The out-of-plane flows stream toward the midplane and meet to drive the SSs. The SSs differ greatly from those in the MHD model. The turbulence at the shock fronts and both upstream and downstream is weak so the shocks are laminar and produce little dissipation. Downstream of the SSs the counterstreaming ion beams lead to higher density, which leads to a positive potential between the SSs that acts to confine the downstream electrons to maintain charge neutrality. The potential accelerates electrons from upstream of the SSs to downstream region and traps a small fraction but only modestly increases the downstream electron temperature above the upstream value. In the low-β limit the released magnetic energy is split between bulk flow and ion heating with little energy going to electrons. That the model does not produce strong electron heating nor an energetic electron component suggests that other mechanisms, such as multiple x-line reconnection, are required to explain energetic electron production in large flares. The model can be tested with the expected data from the Parker Solar Probe.1

No abstract