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

Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.

We develop first-principles theory of kinetic plasma turbulence governed by the Vlasov-Maxwell-Landau equations in the limit of vanishing collision rates. Following an exact renormalization-group approach pioneered by Onsager, we demonstrate the existence of a "collisionless range" of scales (lengths and velocities) in 1-particle phase space where the ideal Vlasov-Maxwell equations are satisfied in a "coarse-grained sense". Entropy conservation may nevertheless be violated in that range by a "dissipative anomaly" due to nonlinear entropy cascade. We derive "4/5th-law" type expressions for the entropy flux, which allow us to characterize the singularities (structure-function scaling exponents) required for its non-vanishing. Conservation laws of mass, momentum and energy are not afflicted with anomalous transfers in the collisionless limit. In a subsequent limit of small gyroradii, however, anomalous contributions to inertial-range energy balance may appear due both to cascade of bulk energy and to turbulent redistribution of internal energy in phase space. In that same limit the "generalized Ohm's law" derived from the particle momentum balances reduces to an "ideal Ohm's law", but only in a coarse-grained sense that does not imply magnetic flux-freezing and that permits magnetic reconnection at all inertial-range scales. We compare our results with prior theory based on the gyrokinetic (high gyro-frequency) limit, with numerical simulations, and with spacecraft measurements of the solar wind and terrestrial magnetosphere.

We develop first-principles theory of relativistic fluid turbulence at high Reynolds and P\'eclet numbers. We follow an exact approach pioneered by Onsager, which we explain as a non-perturbative application of the principle of renormalization-group invariance. We obtain results very similar to those for non-relativistic turbulence, with hydrodynamic fields in the inertial-range described as distributional or "coarse-grained" solutions of the relativistic Euler equations. These solutions do not, however, satisfy the naive conservation-laws of smooth Euler solutions but are afflicted with dissipative anomalies in the balance equations of internal energy and entropy. The anomalies are shown to be possible by exactly two mechanisms, local cascade and pressure-work defect. We derive "4/5th-law"-type expressions for the anomalies, which allow us to characterize the singularities (structure-function scaling exponents) required for their non-vanishing. We also investigate the Lorentz covariance of the inertial-range fluxes, which we find is broken by our coarse-graining regularization but which is restored in the limit that the regularization is removed, similar to relativistic lattice quantum field theory. In the formal limit as speed of light goes to infinity, we recover the results of previous non-relativistic theory. In particular, anomalous heat input to relativistic internal energy coincides in that limit with anomalous dissipation of non-relativistic kinetic energy.

We investigate dissipative anomalies in a turbulent fluid governed by the compressible Navier-Stokes equation. We follow an exact approach pioneered by Onsager, which we explain as a non-perturbative application of the principle of renormalization-group invariance. In the limit of high Reynolds and P\'eclet numbers, the flow realizations are found to be described as distributional or "coarse-grained" solutions of the compressible Euler equations, with standard conservation laws broken by turbulent anomalies. The anomalous dissipation of kinetic energy is shown to be due not only to local cascade, but also to a distinct mechanism called pressure-work defect. Irreversible heating in stationary, planar shocks with an ideal-gas equation of state exemplifies the second mechanism. Entropy conservation anomalies are also found to occur by two mechanisms: an anomalous input of negative entropy (negentropy) by pressure-work and a cascade of negentropy to small scales. We derive "4/5th-law"-type expressions for the anomalies, which allow us to characterize the singularities (structure-function scaling exponents) required to sustain the cascades. We compare our approach with alternative theories and empirical evidence. It is argued that the "Big Power-Law in the Sky" observed in electron density scintillations in the interstellar medium is a manifestation of a forward negentropy cascade, or an inverse cascade of usual thermodynamic entropy.

In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the solar wind, however, with a hyperbolic flux-tube or "$X$-line" extending over integral length-scales.

Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can transit to the turbulence state, what is demonstrated by 3D high resolution numerical simulations. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (1999) reconnection model and present numerical evidence supporting the model and demonstrate that it is closely connected to the concept of Richardson diffusion and compatible with the Lagrangian dynamics of magnetized fluids. We point out that the Generalized Ohm's Law, that accounts for turbulent motion, predicts the subdominance of the microphysical plasma effects for a realistically turbulent media. We show that on of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a corner stone of most theories dealing with magnetized plasmas and therefore its change induces fundamental shifts in accepted paradigms like turbulent reconnection entailing the diffusion process that is essential for understanding star formation. We argue, that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares predicted by turbulent reconnection and relate them to solar flares and gamma ray bursts. We analyze solar observations, measurements in the solar wind or heliospheric current sheet, and show their correspondence with turbulent reconnection predictions. Finally, we discuss 1st Order Fermi acceleration as a natural consequence of the turbulent reconnection.

Plasma flows with an MHD-like turbulent inertial range, such as the solar wind, require a generalization of General Magnetic Reconnection (GMR) theory. We introduce the slip-velocity source vector, which gives the rate of development of slip velocity per unit arc length of field line. The slip source vector is the ratio of the curl of the non ideal electric field in the Generalized Ohm's Law and the magnetic field strength. It diverges at magnetic nulls, unifying GMR with magnetic null-point reconnection. Only under restrictive assumptions is the slip velocity related to the gradient of the quasi potential (integral of parallel electric field along field lines). In a turbulent inertial range the curl becomes extremely large while the parallel component is tiny, so that line slippage occurs even while ideal MHD becomes accurate. The resolution of this paradox is that ideal MHD is valid for a turbulent inertial-range only in a weak sense which does not imply magnetic line freezing. The notion of weak solution is explained in terms of spatial coarse-graining and renormalization group (RG) theory. We give an argument for the weak validity of the ideal Ohm's law in the inertial range, via rigorous estimates of the terms in the Generalized Ohm's Law for an electron-ion plasma. All of the nonideal terms (from collisional resistivity, Hall field, electron pressure anisotropy, and electron inertia) are shown to be irrelevant in the RG sense and large-scale reconnection is thus governed solely by ideal dynamics. We briefly discuss some implications for heliospheric reconnection, in particular for deviations from the Parker spiral model of interplanetary magnetic field. Solar wind observations show that reconnection in a turbulence broadened heliospheric current sheet, consistent with the Lazarian-Vishniac theory, leads to slip velocities that cause field lines to lag relative to the spiral model.

Magnetic reconnection is a fundamental process of magnetic field topology change. We analyze the connection of this process with turbulence which is ubiquitous in astrophysical environments. We show how Lazarian & Vishniac (1999) model of turbulent reconnection is connected to the experimentally proven concept of Richardson diffusion and discuss how turbulence violates the generally accepted notion of magnetic flux freezing. We note that in environments that are laminar initially turbulence can develop as a result of magnetic reconnection and this can result in flares of magnetic reconnection in magnetically dominated media. In particular, magnetic reconnection can initially develop through tearing, but the transition to the turbulent state is expected for astrophysical systems. We show that turbulent reconnection predictions corresponds to the Solar and solar wind data.

As an alternative way of describing the cosmological velocity field, we discuss the evolution of rotational invariants constructed from the velocity gradient tensor. Compared with the traditional divergence-vorticity decomposition, these invariants, defined as coefficients of characteristic equation of the velocity gradient tensor, enable a complete classification of all possible flow patterns in the dark-matter comoving frame, including both potential and vortical flows. We show that this tool, first introduced in turbulence two decades ago, proves to be very useful in understanding the evolution of the cosmic web structure, and in classifying its morphology. Before shell-crossing, different categories of potential flow are highly associated with cosmic web structure, because of the coherent evolution of density and velocity. This correspondence is even preserved at some level when vorticity is generated after shell-crossing. The evolution from the potential to vortical flow can be traced continuously by these invariants. With the help of this tool, we show that the vorticity is generated in a particular way that is highly correlated with the large-scale structure. This includes a distinct spatial distribution and different types of alignment between cosmic web and vorticity direction for various vortical flows. Incorporating shell-crossing into closed dynamical systems is highly non-trivial, but we propose a possible statistical explanation for some of these phenomena relating to the internal structure of the three-dimensional invariants space.

Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Traditionally, magnetic reconnection was associated mostly with solar flares. In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen in field lines and magnetic reconnection. We consider magnetic reconnection in realistic 3D geometry in the presence of turbulence. This turbulence in most astrophysical settings is of pre-existing nature, but it also can be induced by magnetic reconnection itself. In this situation turbulent magnetic field wandering opens up reconnection outflow regions, making reconnection fast. We discuss Lazarian \& Vishniac (1999) model of turbulent reconnection, its numerical and observational testings, as well as its connection to the modern understanding of the Lagrangian properties of turbulent fluids. We show that the predicted dependences of the reconnection rates on the level of MHD turbulence make the generally accepted Goldreich \& Sridhar (1995) model of turbulence self-consistent. Similarly, we argue that the well-known Alfv\'en theorem on flux freezing is not valid for the turbulent fluids and therefore magnetic fields diffuse within turbulent volumes. This is an element of magnetic field dynamics that was not accounted by earlier theories. For instance, the theory of star formation that was developing assuming that it is only the drift of neutrals that can violate the otherwise perfect flux freezing, is affected and we discuss the consequences of the turbulent diffusion of magnetic fields mediated by reconnection.

Astrophysical fluids are generically turbulent and this must be taken into account for most transport processes. We discuss how the preexisting turbulence modifies magnetic reconnection and how magnetic reconnection affects the MHD turbulent cascade. We show the intrinsic interdependence and interrelation of magnetic turbulence and magnetic reconnection, in particular, that strong magnetic turbulence in 3D requires reconnection and 3D magnetic turbulence entails fast reconnection. We follow the approach in Eyink, Lazarian & Vishniac 2011 to show that the expressions of fast magnetic reconnection in Lazarian & Vishniac 1999 can be recovered if Richardson diffusion of turbulent flows is used instead of ordinary Ohmic diffusion. This does not revive, however, the concept of magnetic turbulent diffusion which assumes that magnetic fields can be mixed up in a passive way down to a very small dissipation scales. On the contrary, we are dealing the reconnection of dynamically important magnetic field bundles which strongly resist bending and have well defined mean direction weakly perturbed by turbulence. We argue that in the presence of turbulence the very concept of flux-freezing requires modification. The diffusion that arises from magnetic turbulence can be called reconnection diffusion as it based on reconnection of magnetic field lines. The reconnection diffusion has important implications for the continuous transport processes in magnetized plasmas and for star formation. In addition, fast magnetic reconnection in turbulent media induces the First order Fermi acceleration of energetic particles, can explain solar flares and gamma ray bursts. However, the most dramatic consequence of these developments is the fact that the standard flux freezing concept must be radically modified in the presence of turbulence.

Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfv\'en flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy spectra. Lagrangian particle trajectories in such environments become "spontaneously stochastic", so that infinitely-many magnetic field-lines are advected to each point and must be averaged to obtain the resultant magnetic field. The relative distance between initial magnetic field lines which arrive to the same final point depends upon the properties of two-particle turbulent dispersion. We develop predictions based on the phenomenological Goldreich & Sridhar theory of strong MHD turbulence and on weak MHD turbulence theory. We recover the predictions of the Lazarian & Vishniac theory for the reconnection rate of large-scale magnetic structures. Lazarian & Vishniac also invoked "spontaneous stochasticity", but of the field-lines rather than of the Lagrangian trajectories. More recent theories of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in most astrophysical environments. For example, the predictions of the Lazarian-Vishniac theory are unchanged in Hall MHD turbulence with an extended inertial range, whenever the ion skin depth $\delta_i$ is much smaller than the turbulent integral length or injection-scale $L_i.$

We study the Lagrangian mechanism of the fluctuation dynamo at zero Prandtl number and infinite magnetic Reynolds number, in the Kazantsev-Kraichnan model of white-noise advection. With a rough velocity field corresponding to a turbulent inertial-range, flux-freezing holds only in a stochastic sense. We show that field-lines arriving to the same point which were initially separated by many resistive lengths are important to the dynamo. Magnetic vectors of the seed field that point parallel to the initial separation vector arrive anti-correlated and produce an "anti-dynamo" effect. We also study the problem of "magnetic induction" of a spatially uniform seed field. We find no essential distinction between this process and fluctuation dynamo, both producing the same growth-rates and small-scale magnetic correlations. In the regime of very rough velocity fields where fluctuation dynamo fails, we obtain the induced magnetic energy spectra. We use these results to evaluate theories proposed for magnetic spectra in laboratory experiments of turbulent induction

We investigate the scale-locality of cascades of conserved invariants at high kinetic and magnetic Reynolds numbers in the ``inertial-inductive range'' of magnetohydrodynamic (MHD) turbulence, where velocity and magnetic field increments exhibit suitable power-law scaling. We prove that fluxes of total energy and cross-helicity--or, equivalently, fluxes of Els\"asser energies--are dominated by the contributions of local triads. Corresponding spectral transfers are also scale-local when defined using octave wavenumber bands. Flux and transfer of magnetic helicity may be dominated by non-local triads. The magnetic stretching term also may be dominated by non-local triads but we prove that it can convert energy only between velocity and magnetic modes at comparable scales. We explain the disagreement with numerical studies that have claimed conversion nonlocally between disparate scales. We present supporting data from a $1024^3$ simulation of forced MHD turbulence.

This paper presents both rigorous results and physical theory on the breakdown of magnetic flux conservation for ideal plasmas, by nonlinear effects. Our analysis is based upon an effective equation for magnetohydrodynamic (MHD) modes at length-scales $>\ell,$ with smaller scales eliminated, as in renormalization-group methodology. We prove that flux-conservation can be violated for an arbitrarily small length-scale $\ell,$ and in the absence of any non-ideality, but only if singular current sheets and vortex sheets both exist and intersect in sets of large enough dimension. This result gives analytical support to and rigorous constraints on theories of fast turbulent reconnection. Mathematically, our theorem is analogous to Onsager's result on energy dissipation anomaly in hydrodynamic turbulence. As a physical phenomenon, the breakdown of magnetic-flux conservation in ideal MHD is similar to the decay of magnetic flux through a narrow superconducting ring, by phase-slip of quantized flux lines. The effect should be observable both in numerical MHD simulations and in laboratory plasma experiments at moderately high magnetic Reynolds numbers.