Properties of plasmoid-dominated turbulent reconnection in a low-$\beta$ background plasma are investigated by resistive magnetohydrodynamic (MHD) simulations. In the $\beta_{\rm in} < 1$ regime, where $\beta_{\rm in}$ is plasma $\beta$ in the inflow region, the reconnection site is dominated by shocks and shock-related structures and plasma compression is significant. The effective reconnection rate increases from $0.01$ to $0.02$ as $\beta_{\rm in}$ decreases. We hypothesize that plasma compression allows faster reconnection rate, and then we estimate a speed-up factor, based on a compressible MHD theory. We validate our prediction by a series of MHD simulations. These results suggest that the plasmoid-dominated reconnection can be twice faster than expected in the $\beta \ll 1$ environment in a solar corona.

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

We construct Boris-type schemes for integrating the motion of charged particles in particle-in-cell (PIC) simulation. The new solvers virtually combine the 2-step Boris procedure arbitrary n times in the Lorentz-force part, and therefore we call them the multiple Boris solvers. Using Chebyshev polynomials, a one-step form of the new solvers is provided. The new solvers give n^2 times smaller errors, allow larger timesteps, and have a long-term stability. We present numerical tests of the new solvers, in comparison with other particle integrators.

We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the HLLD approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibited a convergent picture of the early nonlinear growth of the MRI in a global model, and indicated that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.

A simple form of the Boris solver in particle-in-cell (PIC) simulation is proposed. It employs an exact solution of the Lorentz-force part, and it is equivalent to the Boris solver with a gyrophase correction. As a favorable property for stable schemes, this form preserves a volume in the phase space. Numerical tests of the Boris solvers are conducted by test-particle simulations and by PIC simulations. The proposed form provides better accuracy than the popular form, while it only requires few additional computation time.

Basic properties of relativistic magnetic reconnection in electron-positron pair plasmas are investigated by using a particle-in-cell (PIC) simulation. We first revisit a problem by Hesse & Zenitani (2007), who examined the kinetic Ohm's law across the X line. We formulate a relativistic Ohm's law by decomposing the stress-energy tensor. Then, the role of the new term, called the heat-flow inertial term, is examined in the PIC simulation data. We further evaluate the energy balance in the reconnection system. These analyses demonstrate physically transparent ways to diagnose relativistic kinetic data.

Electron dynamics surrounding the X-line in magnetopause-type asymmetric reconnection is investigated using a two-dimensional particle-in-cell simulation. We study electron properties of three characteristic regions in the vicinity of the X-line. The fluid properties, velocity distribution functions (VDFs), and orbits are studied and cross-compared. On the magnetospheric side of the X-line, the normal electric field enhances the electron meandering motion from the magnetosheath side. The motion leads to a crescent-shaped component in the electron VDF, in agreement with recent studies. On the magnetosheath side of the X-line, the magnetic field line is so stretched in the third dimension that its curvature radius is comparable with typical electron Larmor radius. The electron motion becomes nonadiabatic, and therefore the electron idealness is no longer expected to hold. Around the middle of the outflow regions, the electron nonidealness is coincident with the region of the nonadiabatic motion. Finally, we introduce a finite-time mixing fraction (FTMF) to evaluate electron mixing. The FTMF marks the magnetospheric side of the X-line, where the nonideal energy dissipation occurs.

Particle dynamics in the electron current layer in collisionless magnetic reconnection is investigated by using a particle-in-cell simulation. Electron motion and velocity distribution functions are studied by tracking self-consistent trajectories. New classes of electron orbits are discovered: figure-eight-shaped regular orbits inside the electron jet, noncrossing regular orbits on the jet flanks, noncrossing Speiser orbits, and nongyrotropic electrons in the downstream of the jet termination region. Properties of a super-Alfv\'{e}nic outflow jet are attributed to an ensemble of electrons traveling through Speiser orbits. Noncrossing orbits are mediated by the polarization electric field near the electron current layer. The noncrossing electrons are found to be non-negligible in number density. The impact of these new orbits to electron mixing, spatial distribution of energetic electrons, and observational signatures, is presented.

Numerical algorithms to load relativistic Maxwell distributions in particle-in-cell (PIC) and Monte-Carlo simulations are presented. For stationary relativistic Maxwellian, the inverse transform method and the Sobol algorithm are reviewed. To boost particles to obtain relativistic shifted-Maxwellian, two rejection methods are proposed in a physically transparent manner. Their acceptance efficiencies are ${\approx}50\%$ for generic cases and $100\%$ for symmetric distributions. They can be combined with arbitrary base algorithms.

The shock structure of a plasmoid in magnetic reconnection in low-beta plasmas is investigated by two-dimensional magnetohydrodynamic simulations. Using a high-accuracy code with unprecedented resolution, shocks, discontinuities, and their intersections are resolved and clarified. Contact discontinuities emanate from triple-shock intersection points, separating fluids of different origins. Shock-diamonds inside the plasmoid appear to decelerate a supersonic flow. New shock-diamonds and a slow expansion fan are found inside the Petschek outflow. A sufficient condition for the new shock-diamonds and the relevance to astrophysical jets are discussed.

It was recently proposed that the electron-frame dissipation measure, the energy transfer from the electromagnetic field to plasmas in the electron's rest frame, identifies the dissipation region of collisionless magnetic reconnection [Zenitani et al. Phys. Rev. Lett. 106, 195003 (2011)]. The measure is further applied to the electron-scale structures of antiparallel reconnection, by using two-dimensional particle-in-cell (PIC) simulations. The size of the central dissipation region is controlled by the electron-ion mass ratio, suggesting that electron physics is essential. A narrow electron jet extends along the outflow direction until it reaches an electron shock. The jet region appears to be anti-dissipative. At the shock, electron heating is relevant to a magnetic cavity signature. The results are summarized to a unified picture of the single dissipation region in a Hall magnetic geometry.

We investigate the one-dimensional interaction of a relativistic jet and an external medium. Relativistic magnetohydrodynamic simulations show an anomalous boost of the jet fluid in the boundary layer, as previously reported. We describe the boost mechanism using an ideal relativistic fluid and magnetohydrodynamic theory. The kinetic model is also examined for further understanding. Simple scaling laws for the maximum Lorentz factor are derived, and verified by the simulations.

Kinetic aspects of the ion current layer at the center of a reconnection outflow exhaust near the X-type region are investigated by a two-dimensional particle-in-cell (PIC) simulation. The layer consists of magnetized electrons and unmagnetized ions that carry a perpendicular electric current. The ion fluid appears to be nonideal, sub-Alfvenic, and nondissipative. The ion velocity distribution functions contain multiple populations such as global Speiser ions, local Speiser ions, and trapped ions. The particle motion of the local Speiser ions in an appropriately rotated coordinate system explains the ion fluid properties very well. The trapped ions are the first demonstration of the regular orbits in the chaotic particle dynamics [Chen and Palmadesso, J. Geophys. Res., 91, 1499 (1986)] in self-consistent PIC simulations. They would be observational signatures in the ion current layer near reconnection sites.

Plasmoid structures in fast reconnection in low-beta plasmas are investigated by two-dimensional magnetohydrodynamic simulations. A high-resolution shock-capturing code enables us to explore a variety of shock structures: vertical slow shocks behind the plasmoid, another slow shocks in the outer-region, and the shock-reflection in the front side. The Kelvin-Helmholtz-like turbulence is also found inside the plasmoid. It is concluded that these shocks are rigorous features in reconnection in low-beta plasmas, where the reconnection jet speed or the upstream Alfven speed exceeds the sound speed.

The structure of the diffusion regions in antiparallel magnetic reconnection is investigated by means of a theory and a Vlasov simulation. The magnetic diffusion is considered as relaxation to the frozen-in state, which depends on a reference velocity field. A field-aligned component of the frozen-in condition is proposed to evaluate a diffusion-like process. Diffusion signatures with respect to ion and electron bulk flows indicate the ion and electron diffusion regions near the reconnection site. The electron diffusion region resembles the energy dissipation region. These results are favorable to a previous expectation that an electron-scale dissipation region is surrounded by an ion-scale Hall-physics region.

Particle acceleration in the magnetic reconnection of electron-positron plasmas is studied by using a particle-in-cell simulation. It is found that a significantly large number of nonthermal particles are generated by the inductive electric fields around an X-type neutral line when the reconnection outflow velocity, which is known to be an Alfv\'{e}n velocity, is on the order of the speed of light. In such a relativistic reconnection regime, we also find that electrons and positrons form a power-law-like energy distribution through their drift along the reconnection electric field under the relativistic Speiser motion. A brief discussion of the relevance of these results to the current sheet structure, which has an antiparallel magnetic field in astrophysical sources of synchrotron radiation, is presented.

Magnetic reconnection associated to the double tearing mode (DTM) is investigated by means of resistive relativistic magnetohydrodynamic (RRMHD) simulations. A linearly unstable double current sheet system in two dimensional cartesian geometry is considered. For initial perturbations of large enough longitudinal wavelengths, a fast reconnection event is triggered by a secondary instability that is structurally driven by the nonlinear evolution of the magnetic islands. The latter reconnection phase and time scale appear to weakly depend on the plasma resistivity and magnetization parameter. We discuss the possible role of such explosive reconnection dynamics to explain the MeV flares observed in the Crab pulsar nebula. Indeed the time scale and the critical minimum wavelength give constraints on the Lorentz factor of the striped wind and on the location of the emission region respectively.

Signatures of the dissipation region of collisionless magnetic reconnection are investigated by the Geotail spacecraft for the 15 May 2003 event. The energy dissipation in the rest frame of the electron's bulk flow is considered in an approximate form D*_e, which is validated by a particle-in-cell simulation. The dissipation measure is directly evaluated from the {plasma moments}, the electric field, and the magnetic field. Using D*_e, a compact dissipation region is successfully detected in the vicinity of the possible X-point in Geotail data. The dissipation rate is 45 pWm**{-3}. The length of the dissipation region is estimated to 1--2 local ion inertial length. The Lorentz work W, the work rate by Lorentz force to plasmas, is also introduced. It is positive over the reconnection region and it has a peak around the pileup region away from the X-point. These new measures D*_e and W provide useful information to understand the reconnection structure.

The fluid-scale evolution of relativistic magnetic reconnection is investigated by using two-fluid and magnetofluid simulation models. Relativistic two-fluid simulations demonstrate the meso-scale evolution beyond the kinetic scales, and exhibit quasi-steady Petschek-type reconnection. Resistive relativistic MHD simulations further show new shock structures in and around the downstream magnetic island (plasmoid). General discussions on these models are presented.

A new measure to identify a small-scale dissipation region in collisionless magnetic reconnection is proposed. The energy transfer from the electromagnetic field to plasmas in the electron's rest frame is formulated as a Lorentz-invariant scalar quantity. The measure is tested by two-dimensional particle-in-cell simulations in typical configurations: symmetric and asymmetric reconnection, with and without the guide field. The innermost region surrounding the reconnection site is accurately located in all cases. We further discuss implications for nonideal MHD dissipation.

Resistive relativistic magnetohydrodynamic (RRMHD) simulations are applied to investigate the system evolution of relativistic magnetic reconnection. A time-split Harten--Lan--van Leer method is employed. Under a localized resistivity, the system exhibits a fast reconnection jet with an Alfv\'{e}nic Lorentz factor inside a narrow Petschek-type exhaust. Various shock structures are resolved in and around the plasmoid such as the post-plasmoid vertical shocks and the "diamond-chain" structure due to multiple shock reflections. Under a uniform resistivity, Sweet--Parker-type reconnection slowly evolves. Under a current-dependent resistivity, plasmoids are repeatedly formed in an elongated current sheet. It is concluded that the resistivity model is of critical importance for RRMHD modeling of relativistic magnetic reconnection.

The nonlinear evolution of relativistic magnetic reconnection in sheared magnetic configuration (with a guide field) is investigated by using two-dimensional relativistic two-fluid simulations. Relativistic guide field reconnection features the charge separation and the guide field compression in and around the outflow channel. As the guide field increases, the composition of the outgoing energy changes from enthalpy-dominated to Poynting-dominated. The inertial effects of the two-fluid model play an important role to sustain magnetic reconnection. Implications for the single fluid magnetohydrodynamic approach and the physics models of relativistic reconnection are briefly addressed.

We investigate the large scale evolution of a relativistic magnetic reconnection in an electron-positron pair plasma by a relativistic two-fluid magnetohydrodynamic (MHD) code. We introduce an inter-species friction force as an effective resistivity to dissipate magnetic fields. We demonstrate that magnetic reconnection successfully occurs in our two-fluid system, and that it involves Petschek-type bifurcated current layers in later stage. We further observe a quasi-steady evolution thanks to an open boundary condition, and find that the Petschek-type structure is stable over the long time period. Simulation results and theoretical analyses exhibit that the Petschek outflow channel becomes narrower when the reconnection inflow contains more magnetic energy, as previously claimed. Meanwhile, we find that the reconnection rate goes up to ~1 in extreme cases, which is faster than previously thought. The role of the resistivity, implications for reconnection models in the magnetically dominated limit, and relevance to kinetic reconnection works are discussed.

The linear and non-linear evolution of a relativistic current sheet of pair ($e^{\pm}$) plasmas is investigated by three-dimensional particle-in-cell simulations. In a Harris configuration, it is obtained that the magnetic energy is fast dissipated by the relativistic drift kink instability (RDKI). However, when a current-aligned magnetic field (the so-called "guide field") is introduced, the RDKI is stabilized by the magnetic tension force and it separates into two obliquely-propagating modes, which we call the relativistic drift-kink-tearing instability (RDKTI). These two waves deform the current sheet so that they trigger relativistic magnetic reconnection at a crossover thinning point. Since relativistic reconnection produces a lot of non-thermal particles, the guide field is of critical importance to study the energetics of a relativistic current sheet.

We investigate properties of the reconnecting current layer in relativistic pair plasma reconnection. We found that the current layer self-regulates its thickness when the current layer runs out current carriers, and so relativistic reconnection retains a fast reconnection rate. Constructing a steady state Sweet-Parker model, we discuss conditions for the current sheet expansion. Based on the energy argument, we conclude that the incompressible assumption is invalid in relativistic Sweet-Parker reconnection. The guide field cases are more incompressible than the anti-parallel cases, and we find a more significant current sheet expansion.

The role of the Weibel instability is investigated for the first time in the context of the large-scale magnetic reconnection problem. A late-time evolution of magnetic reconnection in relativistic pair plasmas is demonstrated by particle-in-cell (PIC) simulations. In the outflow regions, powerful reconnection jet piles up the magnetic fields and then a tangential discontinuity appears there. Further downstream, it is found that the two-dimensional extension of the relativistic Weibel instability generates electro-magnetic fields, which are comparable to the anti-parallel or piled-up fields. In a microscopic viewpoint, the instability allows plasma's multiple interactions with the discontinuity. In a macroscopic viewpoint, the instability leads to rapid expansion of the current sheet and then the reconnection jet front further propagates into the downstream. Possible application to the three-dimensional case is briefly discussed.

We study the role of the guide field in relativistic magnetic reconnection in a Harris current sheet of pair ($e^{\pm}$) plasmas, using linear theories and particle-in-cell (PIC) simulations. Two-dimensional PIC simulations exhibit the guide field dependence to the linear instabilities; the tearing or reconnection modes are relatively insensitive, while the relativistic drift-kink instability (RDKI), the fastest mode in a relativistic current sheet, is stabilized by the guide field. Particle acceleration in the nonlinear stage is also investigated. A three-dimensional PIC simulation demonstrates that the current sheet is unstable to the RDKI, although small reconnection occurs in the deformed current sheet. Another three-dimensional PIC simulation with a guide field demonstrates a completely different scenario. Secondary magnetic reconnection is triggered by nonlinear coupling of oblique instabilities, which we call the relativistic drift-sausage tearing instability. Therefore, particle acceleration by relativistic guide field reconnection occurs in three-dimensional configuration. Based on the plasma theories, we discuss an important role of the guide field: to enable non-thermal particle acceleration by magnetic reconnection.

We study linear and nonlinear development of relativistic and ultrarelativistic current sheets of pair plasmas with antiparallel magnetic fields. Two types of two-dimensional problems are investigated by particle-in-cell simulations. First, we present the development of relativistic magnetic reconnection, whose outflow speed is an order of the light speed c. It is demonstrated that particles are strongly accelerated in and around the reconnection region, and that most of magnetic energy is converted into "nonthermal" part of plasma kinetic energy. Second, we present another two-dimensional problem of a current sheet in a cross-field plane. In this case, the relativistic drift kink instability (RDKI) occurs. Particle acceleration also takes place, but the RDKI fast dissipates the magnetic energy into plasma heat. We discuss the mechanism of particle acceleration and the theory of the RDKI in detail. It is important that properties of these two processes are similar in the relativistic regime of T > mc^2, as long as we consider the kinetics. Comparison of the two processes indicates that magnetic dissipation by the RDKI is more favorable process in the relativistic current sheet. Therefore the striped pulsar wind scenario should be reconsidered by the RDKI.

Two-dimensional particle simulations of a relativistic Harris current sheet of pair plasmashave demonstrated that the system is unstable to the relativistic drift kink instability (RDKI) and that a new kind of acceleration process takes place in the deformed current sheet. This process contributes to the generation of non-thermal particles and contributes to the fast magnetic dissipation in the current sheet structure. The acceleration mechanism and a brief comparison with relativistic magnetic reconnection are presented.