We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, {\it viz.} the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in critically-balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field.

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

Recently, it has been realized that magnetic plasma turbulence and magnetic field reconnection are inherently related phenomena. Turbulent fluctuations generate regions of a sheared magnetic field that become unstable to the tearing instability and reconnection, thus modifying turbulence at the corresponding scales. In this contribution, we give a brief discussion of some recent results on tearing-mediated magnetic turbulence. We illustrate the main ideas of this rapidly developing field of study by concentrating on two important examples -- magnetohydrodynamic Alfv\'en turbulence and small-scale kinetic-Alfv\'en turbulence. Due to various potential applications of these phenomena in space physics and astrophysics, we specifically try not to overload the text by heavy analytical derivations, but rather present a qualitative discussion accessible to a non-expert in the theories of turbulence and reconnection.

Recent analytical works on strong magnetized plasma turbulence have hypothesized the existence of a range of scales where the tearing instability may govern the energy cascade. In this paper, we estimate the conditions under which such tearing may give rise to full nonlinear magnetic reconnection in the turbulent eddies, thereby enabling significant energy conversion and dissipation. When those conditions are met, a new turbulence regime is accessed where reconnection-driven energy dissipation becomes common, rather than the rare feature that it must be when they are not.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magnetic energy decays as $\tilde t^{-1}$, where $\tilde t$ is time normalized to the (appropriately defined) reconnection timescale; and the correlation length of the field grows as $\tilde t^{1/2}$. The magnetic energy spectrum is self-similar, and evolves as $\propto \tilde t ^{-3/2}k^{-2}$, where the $k$-dependence is imparted by the formation of thin current sheets.

We investigate analytically and numerically the semi-collisional regime of the plasmoid instability, defined by the inequality $\delta_{SP} \gg \rho_s \gg \delta_{in}$, where $\delta_{SP}$ is the width of a Sweet-Parker current sheet, $\rho_S$ is the ion sound Larmor radius, and $\delta_{in}$ is width of boundary layer that arises in the plasmoid instability analysis. Theoretically, this regime is predicted to exist if the Lundquist number $S$ and the length of the current sheet $L$ are such that $(L/\rho_S)^{14/9} < S <(L/\rho_S)^2$ (for a sinusoidal-like magnetic configuration; for a Harris-type sheet the lower bound is replaced with $(L/\rho_S)^{8/5}$). These bounds are validated numerically by means of simulations using a reduced gyrokinetic model (Zocco & Schekochihin, $\it{Physics~of~Plasmas}$, ${\bf 18}$, 2011) conducted with the code $\tt{Viriato}$. Importantly, this regime is conjectured to allow for plasmoid formation at relatively low, experimentally accessible, values of the Lundquist number. Our simulations obtain plasmoid instability at values of $S$ as low as $\sim 250$. The simulations do not prescribe a Sweet-Parker sheet; rather, one is formed self-consistenly during the nonlinear evolution of the initial tearing mode configuration. This proves that this regime of the plasmoid instability is realizable, at least at the relatively low values of the Lundquist number that are accessible to current dedicated experiments.

Recent studies have suggested that the tearing instability may play a significant role in magnetic turbulence. In this work, we review the theory of the magnetohydrodynamic tearing instability in the general case of an arbitrary tearing parameter, which is relevant for applications in turbulence. We discuss a detailed derivation of the results for the standard Harris profile and accompany it by the derivation of the results for a lesser known sine-shaped profile. We devote special attention to the exact solution of the inner equation, which is the central result in the theory of tearing instability. We also briefly discuss the influence of shear flows on tearing instability in magnetic structures. Our presentation is self-contained; we expect it to be accessible to researchers in plasma turbulence who are not experts in magnetic reconnection.

Alfv\'enic-type turbulence in strongly magnetized, low-beta pair plasmas is investigated. A coupled set of equations for the evolution of the magnetic and flow potentials are derived, covering both fluid and kinetic scales. In the fluid (MHD) range those equations are the same as for electron-ion plasmas, so turbulence at those scales is expected to be of the Alfv\'enic nature, exhibiting critical balance, dynamic alignment, and transition to a tearing mediated regime at small scales. The critical scale at which a transition to a tearing-mediated range occurs is derived, and the spectral slope in that range is predicted to be $k_\perp^{-8/3}$ (or $k_\perp^{-3}$, depending on details of the reconnecting configuration assumed). At scales below the electron (and positron) skin depth, it is argued that turbulence is dictated by a cascade of the inertial Alfv\'en wave, which we show to result in the magnetic energy spectrum $\propto k_\perp^{-11/3}$.

We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $\beta$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the underlying mechanisms for the development of this turbulent state is the drift-kink instability (DKI) of the current sheets resulting from the nonlinear evolution of the channel modes. The onset of the DKI can only be observed for simulation domain sizes exceeding several linear MRI wavelengths. The DKI, together with ensuing magnetic reconnection, activate the turbulent motion of the plasma in the late stage of the nonlinear evolution of the MRI. At steady state, the magnetic energy has an MHD-like spectrum with a slope of $k^{-5/3}$ for $k\rho<1$ and $k^{-3}$ for sub-Larmor scale ($k\rho>1$). We also examine the role of the collisionless MRI and associated magnetic reconnection in the development of pressure anisotropy. We study the stability of the system due to this pressure anisotropy, observing the development of mirror instability during the early-stage of the MRI. We further discuss the importance of magnetic reconnection for particle acceleration during the turbulence regime. In particular, consistent with reconnection studies, we show that at late times the kinetic energy presents a characteristic slope of $\epsilon^{-2}$ in the high-energy region.

We present results from a three-dimensional particle-in-cell simulation of plasma turbulence, resembling the plasma conditions found at kinetic scales of the solar wind. The spectral properties of the turbulence in the subion range are consistent with theoretical expectations for kinetic Alfv\' en waves. Furthermore, we calculate the local anisotropy, defined by the relation $k_{\parallel}(k_{\perp})$, where $k_{\parallel}$ is a characteristic wave number along the local mean magnetic field at perpendicular scale $l_{\perp}\sim 1/k_{\perp}$. The subion range anisotropy is scale dependent with $k_{\parallel}<k_{\perp}$ and the ratio of linear to nonlinear time scales is of order unity, suggesting that the kinetic cascade is close to a state of critical balance. Our results compare favorably against a number of \emph{in situ} solar wind observations and demonstrate---from first principles---the feasibility of plasma turbulence models based on a critically balanced cascade of kinetic Alfv\' en waves.

Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfv\'enic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; Rutherford 1973) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal MHD, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. 1963; Coppi et al. 1976) can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time-dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

We report the results of a direct comparison between different kinetic models of collisionless plasma turbulence in two spatial dimensions. The models considered include a first principles fully-kinetic (FK) description, two widely used reduced models [gyrokinetic (GK) and hybrid-kinetic (HK) with fluid electrons], and a novel reduced gyrokinetic approach (KREHM). Two different ion beta ($\beta_i$) regimes are considered: 0.1 and 0.5. For $\beta_i=0.5$, good agreement between the GK and FK models is found at scales ranging from the ion to the electron gyroradius, thus providing firm evidence for a kinetic Alfv\'en cascade scenario. In the same range, the HK model produces shallower spectral slopes, presumably due to the lack of electron Landau damping. For $\beta_i=0.1$, a detailed analysis of spectral ratios reveals a slight disagreement between the GK and FK descriptions at kinetic scales, even though kinetic Alfv\'en fluctuations likely still play a significant role. The discrepancy can be traced back to scales above the ion gyroradius, where the FK and HK results seem to suggest the presence of fast magnetosonic and ion Bernstein modes in both plasma beta regimes, but with a more notable deviation from GK in the low-beta case. The identified practical limits and strengths of reduced-kinetic approximations, compared here against the fully-kinetic model on a case-by-case basis, may provide valuable insight into the main kinetic effects at play in turbulent collisionless plasmas, such as the solar wind.

It has recently been proposed (Loureiro & Boldyrev 2017; Mallet et al. 2017) that the inertial interval in magnetohydrodynamic (MHD) turbulence is terminated at small scales not by a Kolmogorov-like dissipation region, but rather by a new sub-inertial interval mediated by tearing instability. However, many astrophysical plasmas are nearly collisionless so that the MHD approximation is not applicable to turbulence at small scales. In this Letter, we propose the extension of the theory of reconnection-mediated turbulence to plasmas which are so weakly collisional that the reconnection occurring in the turbulent eddies is caused by electron inertia rather than by resistivity. We find that the transition scale to reconnection-mediated turbulence depends on the plasma beta and on the assumptions of the plasma turbulence model. However, in all cases analyzed, the energy spectra in the reconnection-mediated interval range from $E(k)dk\propto k^{-8/3}dk$ to $E(k)dk\propto k^{-3}dk$.

The current understanding of MHD turbulence envisions turbulent eddies which are anisotropic in all three directions. In the plane perpendicular to the local mean magnetic field, this implies that such eddies become current-sheet-like structures at small scales. We analyze the role of magnetic reconnection in these structures and conclude that reconnection becomes important at a scale $\lambda\sim L S_L^{-4/7}$, where $S_L$ is the outer-scale ($L$) Lundquist number and $\lambda$ is the smallest of the field-perpendicular eddy dimensions. This scale is larger than the scale set by the resistive diffusion of eddies, therefore implying a fundamentally different route to energy dissipation than that predicted by the Kolmogorov-like phenomenology. In particular, our analysis predicts the existence of the sub-inertial, reconnection interval of MHD turbulence, with the Fourier energy spectrum $E(k_\perp)\propto k_\perp^{-5/2}$, where $k_\perp$ is the wave number perpendicular to the local mean magnetic field. The same calculation is also performed for high (perpendicular) magnetic Prandtl number plasmas ($Pm$), where the reconnection scale is found to be $\lambda/L\sim S_L^{-4/7}Pm^{-2/7}$.

The recent realization that Sweet-Parker current sheets are violently unstable to the secondary tearing (plasmoid) instability implies that such current sheets cannot occur in real systems. This suggests that, in order to understand the onset of magnetic reconnection, one needs to consider the growth of the tearing instability in a current layer as it is being formed. Such an analysis is performed here in the context of nonlinear resistive MHD for a generic time-dependent equilibrium representing a gradually forming current sheet. It is shown that two onset regimes, single-island and multi-island, are possible, depending on the rate of current sheet formation. A simple model is used to compute the criterion for transition between these two regimes, as well as the reconnection onset time and the current sheet parameters at that moment. For typical solar corona parameters this model yields results consistent with observations.

(abridged) Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergone a major paradigm shift. The steady-state model of reconnection described by the famous Sweet-Parker (SP) theory, which dominated the field for ~50 years, has been replaced with an essentially time-dependent, bursty picture of the reconnection layer, dominated by the continuous formation and ejection of multiple secondary islands (plasmoids). Whereas in the SP model reconnection was predicted to be slow, a major implication of this new paradigm is that reconnection in fluid systems is fast (i.e., independent of the Lundquist number), provided that the system is large enough. This conceptual shift hinges on the realization that SP-like current layers are violently unstable to the plasmoid instability - implying, therefore, that such current sheets are super-critically unstable and thus can never form in the first place. This suggests that the formation of a current sheet and the subsequent reconnection process cannot be decoupled, as is commonly assumed. This paper provides an introductory-level overview of the recent developments in reconnection theory and simulations that led to this essentially new framework. We briefly discuss the role played by the plasmoid instability in selected applications, and describe some of the outstanding challenges that remain at the frontier of this subject. Amongst these are the analytical and numerical extension of the plasmoid instability to (i) 3D and (ii) non-MHD regimes. New results are reported in both cases.

We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.

Magnetic reconnection in strongly magnetized (low-beta), weakly collisional plasmas is investigated using a novel fluid-kinetic model [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] which retains non-isothermal electron kinetics. It is shown that electron heating via Landau damping (linear phase mixing) is the dominant dissipation mechanism. In time, electron heating occurs after the peak of the reconnection rate; in space, it is concentrated along the separatrices of the magnetic island. For sufficiently large systems, the peak reconnection rate is $cE_{max}\approx 0.2v_AB_{y,0}$, where $v_A$ is the Alfv\'en speed based on the reconnecting field $B_{y,0}$. The island saturation width is the same as in MHD models except for small systems, when it becomes comparable to the kinetic scales.

A 2D linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of Reduced MHD. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, Phys. Plasmas {\bf 14}, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a large-wave-number chain of plasmoids ($k_{\rm max}\Lsheet \sim S^{3/8}$, where $k_{\rm max}$ is the wave-number of fastest growing mode, $S=\Lsheet V_A/\eta$ is the Lundquist number, $\Lsheet$ is the length of the sheet, $V_A$ is the Alfv\'en speed and $\eta$ is the plasma resistivity), which grows super-Alfv\'enically fast ($\gmax\tau_A\sim S^{1/4}$, where $\gmax$ is the maximum growth rate, and $\tau_A=\Lsheet/V_A$). For typical background profiles, the growth rate and the wave-number are found to {\it increase} in the outflow direction. This is due to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at the periphery of the layer, where the outflow velocity exceeds the Alfv\'en speed associated with the upstream magnetic field. The KH instability grows even faster than the plasmoid instability, $\gmax \tau_A \sim k_{\rm max} \Lsheet\sim S^{1/2}$. The effect of viscosity ($\nu$) on the plasmoid instability is also addressed. In the limit of large magnetic Prandtl numbers, $Pm=\nu/\eta$, it is found that $\gmax\sim S^{1/4}Pm^{-5/8}$ and $k_{\rm max} \Lsheet\sim S^{3/8}Pm^{-3/16}$, leading to the prediction that the critical Lundquist number for plasmoid instability in the $Pm\gg1$ regime is $\Scrit\sim 10^4Pm^{1/2}$. These results are verified via direct numerical simulation of the linearized equations, using a new, analytical 2D SP equilibrium solution.

A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {\it et al.} [Phys. Rev. Lett. {\bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is $\normEeff\sim 0.02$ independently of $S$ for $S\gg10^4$. The plasmoid flux ($\Psi$) and half-width ($w_x$) distribution functions scale as $f(\Psi)\sim \Psi^{-2}$ and $f(w_x)\sim w_x^{-2}$. The joint distribution of $\Psi$ and $w_x$ shows that plasmoids populate a triangular region $w_x\gtrsim\Psi/B_0$, where $B_0$ is the reconnecting field. It is argued that this feature is due to plasmoid coalescence. Macroscopic "monster" plasmoids with $w_x\sim 10$% of the system size are shown to emerge in just a few Alfv\'en times, independently of $S$, suggesting that large disruptive events are an inevitable feature of large-$S$ reconnection.

A conceptual model of resistive magnetic reconnection via a stochastic plasmoid chain is proposed. The global reconnection rate is shown to be independent of the Lundquist number. The distribution of fluxes in the plasmoids is shown to be an inverse square law. It is argued that there is a finite probability of emergence of abnormally large plasmoids, which can disrupt the chain (and may be responsible for observable large abrupt events in solar flares and sawtooth crashes). A criterion for the transition from magnetohydrodynamic to collisionless regime is provided.

Two-dimensional numerical simulations of the effect of background turbulence on 2D resistive magnetic reconnection are presented. For sufficiently small values of the resistivity ($\eta$) and moderate values of the turbulent power ($\epsilon$), the reconnection rate is found to have a much weaker dependence on $\eta$ than the Sweet-Parker scaling of $\eta^{1/2}$ and is even consistent with an $\eta-$independent value. For a given value of $\eta$, the dependence of the reconnection rate on the turbulent power exhibits a critical threshold in $\epsilon$ above which the reconnection rate is significantly enhanced.

A detailed numerical study of magnetic reconnection in resistive MHD for very large, previously inaccessible, Lundquist numbers ($10^4\le S\le 10^8$) is reported. Large-aspect-ratio Sweet-Parker current sheets are shown to be unstable to super-Alfv\'enically fast formation of plasmoid (magnetic-island) chains. The plasmoid number scales as $S^{3/8}$ and the instability growth rate in the linear stage as $S^{1/4}$, in agreement with the theory by Loureiro et al. [Phys. Plasmas {\bf 14}, 100703 (2007)]. In the nonlinear regime, plasmoids continue to grow faster than they are ejected and completely disrupt the reconnection layer. These results suggest that high-Lundquist-number reconnection is inherently time-dependent and hence call for a substantial revision of the standard Sweet-Parker quasi-stationary picture for $S>10^4$.

Current sheets formed in magnetic reconnection events are found to be unstable to high-wavenumber perturbations. The instability is very fast: its maximum growth rate scales as S^{1/4} v_A/L, where L is the length of the sheet, v_A the Alfven speed and S the Lundquist number. As a result, a chain of plasmoids (secondary islands) is formed, whose number scales as S^{3/8}.

We study the nonlinear evolution of the resistive tearing mode in slab geometry in two dimensions. We show that, in the strongly driven regime (large Delta'), a collapse of the X-point occurs once the island width exceeds a certain critical value ~1/Delta'. A current sheet is formed and the reconnection is exponential in time with a growth rate ~eta^1/2, where eta is the resistivity. If the aspect ratio of the current sheet is sufficiently large, the sheet can itself become tearing-mode unstable, giving rise to secondary islands, which then coalesce with the original island. The saturated state depends on the value of Delta'. For small Delta', the saturation amplitude is ~Delta' and quantitatively agrees with the theoretical prediction. If Delta' is large enough for the X-point collapse to have occured, the saturation amplitude increases noticeably and becomes independent of Delta'.