Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent advances in particle acceleration and reconnection physics in the magnetically-dominated regime. More discussion is focused on the physics of particle acceleration, power-law formation as well as the reconnection rate problem. In addition, we provide an outlook for studying reconnection acceleration mechanisms and kinetic physics in the next step.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, since it opens up the possiblity of modeling the energetic spectra without resolving microscopic diffusion regions.

Within the resistive magnetohydrodynamic model, high-Lundquist number reconnection layers are unstable to the plasmoid instability, leading to a turbulent evolution where the reconnection rate can be independent of the underlying resistivity. However, the physical relevance of these results remains questionable for many applications. First, the reconnection electric field is often well above the runaway limit, implying that collisional resistivity is invalid. Furthermore, both theory and simulations suggest that plasmoid formation may rapidly induce a transition to kinetic scales, due to the formation of thin current sheets. Here, this problem is studied for the first time using a first-principles kinetic simulation with a Fokker-Planck collision operator in 3D. The low-$\beta$ reconnecting current layer thins rapidly due to Joule heating before onset of the oblique plasmoid instability. Linear growth rates for standard ($k_y = 0$) tearing modes agree with semi-collisional boundary layer theory, but the angular spectrum of oblique ($|k_y|>0$) modes is significantly narrower than predicted. In the non-linear regime, flux-ropes formed by the instability undergo complex interactions as they are advected and rotated by the reconnection outflow jets, leading to a turbulent state with stochastic magnetic field. In a manner similar to previous 2D results, super-Dreicer fields induce a transition to kinetic reconnection in thin current layers that form between flux-ropes. These results may be testable within new laboratory experiments.

Based on the rate of expansion of the solar wind, the plasma should cool rapidly as a function of distance to the Sun. Observations show this is not the case. In this work, a magnetic pumping model is developed as a possible explanation for the heating and the generation of power-law distribution functions observed in the solar wind plasma. Most previous studies in this area focus on the role that the dissipation of turbulent energy on microscopic kinetic scales plays in the overall heating of the plasma. However, with magnetic pumping particles are energized by the largest scale turbulent fluctuations, thus bypassing the energy cascade. In contrast to other models, we include the pressure anisotropy term, providing a channel for the large scale fluctuations to heat the plasma directly. In this work a complete set of coupled differential equations describing the evolution, and energization, of the distribution function are derived, as well as an approximate closed form solution. Numerical simulations using the VPIC kinetic code are applied to verify the model's analytical predictions. The results of the model for realistic solar wind scenario are computed, where thermal streaming of particles are important for generating a phase shift between the magnetic perturbations and the pressure anisotropy. In turn, averaged over a pump cycle, the phase shift permits mechanical work to be converted directly to heat in the plasma. The results of this scenario show that magnetic pumping may account for a significant portion of the solar wind energization.

Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients.

Simulations suggest collisionless steady-state magnetic reconnection of Harris-type current sheets proceeds with a rate of order 0.1, independent of dissipation mechanism. We argue this long-standing puzzle is a result of constraints at the magnetohydrodynamic (MHD) scale. We perform a scaling analysis of the reconnection rate as a function of the opening angle made by the upstream magnetic fields, finding a maximum reconnection rate close to 0.2. The predictions compare favorably to particle-in-cell simulations of relativistic electron-positron and non-relativistic electron-proton reconnection. The fact that simulated reconnection rates are close to the predicted maximum suggests reconnection proceeds near the most efficient state allowed at the MHD-scale. The rate near the maximum is relatively insensitive to the opening angle, potentially explaining why reconnection has a similar fast rate in differing models.

A model of global magnetic reconnection rate in relativistic collisionless plasmas is developed and validated by the fully kinetic simulation. Through considering the force balance at the upstream and downstream of the diffusion region, we show that the global rate is bounded by a value $\sim 0.3$ even when the local rate goes up to $\sim O(1)$ and the local inflow speed approaches the speed of light in strongly magnetized plasmas. The derived model is general and can be applied to magnetic reconnection under widely different circumstances.

Plasma energization through magnetic reconnection in the magnetically-dominated regime featured by low plasma beta ($\beta = 8 \pi nkT_0/B^2 \ll 1$) and/or high magnetization ($\sigma = B^2/(4 \pi nmc^2) \gg 1$) is important in a series of astrophysical systems such as solar flares, pulsar wind nebula, and relativistic jets from black holes, etc. In this paper, we review the recent progress on kinetic simulations of this process and further discuss plasma dynamics and particle acceleration in a low-$\beta$ reconnection layer that consists of electron-positron pairs. We also examine the effect of different initial thermal temperatures on the resulting particle energy spectra. While earlier papers have concluded that the spectral index is smaller for higher $\sigma$, our simulations show that the spectral index approaches $p=1$ for sufficiently low plasma $\beta$, even if $\sigma \sim 1$. Since this predicted spectral index in the idealized limit is harder than most observations, it is important to consider effects that can lead to a softer spectrum such as open boundary simulations. We also remark that the effects of 3D reconnection physics and turbulence on reconnection need to be addressed in the future.

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $\sigma_0=B^2/(4\pi n(m_i+m_e)c^2)$, the magnetization for each species is $\sigma_i \sim \sigma_0$ and $\sigma_e \sim (m_i/m_e) \sigma_0$, respectively. We have studied the magnetically dominated regime by varying $\sigma_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $\gamma_{be} \sim \sigma_e$ and $\gamma_{bi} \sim \sigma_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv\'{e}n waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of $k_{\perp}^{-1.3}$. The kinetic code shows a spectral slope of $k_{\perp}^{-1.5}$ for smaller simulation domain, and $k_{\perp}^{-1.3}$ for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.

The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter, and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.

Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two- and three-dimensional kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron-positron plasmas starting with a magnetically dominated, force-free current sheet ($\sigma \equiv B^2/(4\pi n_e m_e c^2) \gg 1$). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra $f \propto (\gamma-1)^{-p}$ and approaches $p = 1$ for sufficiently large $\sigma$ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to non-relativistic regimes. In the three-dimensional simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications to a variety of high-energy astrophysical problems.

Using fully kinetic simulations, we demonstrate that magnetic reconnection in relativistic plasmas is highly efficient at accelerating particles through a first-order Fermi process resulting from the curvature drift of particles in the direction of the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra in parameter regimes where the energy density in the reconnecting field exceeds the rest mass energy density $\sigma \equiv B^2/(4 \pi n m_ec^2) > 1$ and when the system size is sufficiently large. In the limit $\sigma \gg 1$, the spectral index approaches $p=1$ and most of the available energy is converted into non-thermal particles. A simple analytic model is proposed which explains these key features and predicts a general condition under which hard power-law spectra will be generated from magnetic reconnection.

Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter $\sigma$ and approaches the light speed when $\sigma > O(100)$, leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains $\sim 0.1$ in both the non-relativistic and relativistic limits.

Results of the first validation of large guide field, $B_g / \delta B_0 \gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $\beta_i = 0.01$ and solar wind relevant $\beta_i = 1$ are presented, where $\delta B_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $\beta_i B_g / \delta B_0 \gtrsim 1$ and $B_g / \delta B_0 \gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.

Recent progress in understanding the physics of magnetic reconnection is conveniently summarized in terms of a phase diagram which organizes the essential dynamics for a wide variety of applications in heliophysics, laboratory and astrophysics. The two key dimensionless parameters are the Lundquist number and the macrosopic system size in units of the ion sound gyroradius. In addition to the conventional single X-line collisional and collisionless phases, multiple X-line reconnection phases arise due to the presence of the plasmoid instability either in collisional and collisionless current sheets. In particular, there exists a unique phase termed "multiple X-line hybrid phase" where a hierarchy of collisional islands or plasmoids is terminated by a collisionless current sheet, resulting in a rapid coupling between the macroscopic and kinetic scales and a mixture of collisional and collisionless dynamics. The new phases involving multiple X-lines and collisionless physics may be important for the emerging applications of magnetic reconnection to accelerate charged particles beyond their thermal speeds. A large number of heliophysical and astrophysical plasmas are surveyed and grouped in the phase diagram: Earth's magnetosphere, solar plasmas (chromosphere, corona, wind and tachocline), galactic plasmas (molecular clouds, interstellar media, accretion disks and their coronae, Crab nebula, Sgr A*, gamma ray bursts, magnetars), extragalactic plasmas (Active Galactic Nuclei disks and their coronae, galaxy clusters, radio lobes, and extragalactic jets). Significance of laboratory experiments, including a next generation reconnection experiment, is also discussed.

Detailed comparisons are reported between laboratory observations of electron-scale dissipation layers near a reconnecting X-line and direct two-dimensional full-particle simulations. Many experimental features of the electron layers, such as insensitivity to the ion mass, are reproduced by the simulations; the layer thickness, however, is about 3-5 times larger than the predictions. Consequently, the leading candidate 2D mechanism based on collisionless electron nongyrotropic pressure is insufficient to explain the observed reconnection rates. These results suggest that, in addition to the residual collisions, 3D effects play an important role in electron-scale dissipation during fast reconnection.

The results of kinetic simulations of magnetic reconnection in Harris current sheets are analyzed. A range of guide fields is considered to study reconnection in plasmas characterized by different $\beta$ values, $\beta > m_e/m_i$. Both an implicit Particle-in-Cell (PIC) simulation method and a parallel explicit PIC code are used. Simulations with mass ratios up to the physical value are performed. The simulations show that the reconnection rate decreases with the guide field and depends weakly on the mass ratio. The off-diagonal components of the electron pressure tensor break the frozen-in condition, even in low $\beta$ plasmas. In high $\beta$ plasmas, evidence is presented that whistler waves play a key role in the enhanced reconnection, while in low $\beta$ plasmas the kinetic Alfv\'en waves are important. The in-plane and the out-of-plane ion and electron motion are also considered, showing that they are influenced by the mass ratio and the plasma $\beta$.