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

Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfv\'enic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfv\'en velocity shear, whereas the median shear across 43 Alfv\'enic current sheets examined was 71% of the Alfv\'en velocity shear. This finding could suggest that large, albeit sub-Alfv\'enic, velocity shears suppress reconnection. An alternative interpretation is that the Alfv\'enic current sheets are isolated rotational discontinuities which do not undergo local reconnection.

We study spectral features of ion velocity and magnetic field correlations in the solar wind and in the magnetosheath using data from the Magnetospheric Multi-Scale (MMS) spacecraft. High resolution MMS observations enable the study of transition of these correlations between their magnetofluid character at larger scales into the sub-proton kinetic range, previously unstudied in spacecraft data. Cross-helicity, angular alignment and energy partitioning is examined over a suit- able range of scales, employing measurements based on the Taylor frozen-in approximation as well as direct two-spacecraft correlation measurements. The results demonstrate signatures of alignment at large scales. As kinetic scales are approached, the alignment between v and b is destroyed by demagnetization of protons.

We study the effects of plasma \b{eta} (ratio of plasma pressure to magnetic pressure) on the evolution of kinetic plasma turbulence using fully kinetic particle-in-cell simulations of decaying turbulence. We find that the plasma \b{eta} systematically affects spectra, measures of intermittency, decay rates of turbulence fluctuations, and partitioning over different channels of energy exchange More specifically, an increase in plasma \b{eta} leads to greater total heating, with proton heating preferentially more than electrons. Implications for achieving magnetosheath like temperature ratios are discussed.

Protons (ionized hydrogen) in the solar wind frequently exhibit distinct temperatures ($T_{\perp p}$ and $T_{\parallel p}$) perpendicular and parallel to the plasma's background magnetic-field. Numerous prior studies of the interplanetary solar-wind have shown that, as plasma beta ($\beta_{\parallel p}$) increases, a narrower range of temperature-anisotropy ($R_p\equiv T_{\perp p}\,/\,T_{\parallel p}$) values is observed. Conventionally, this effect has been ascribed to the actions of kinetic microinstabilities. This study is the first to use data from the Magnetospheric Multiscale Mission (MMS) to explore such $\beta_{\parallel p}$-dependent limits on $R_p$ in Earth's magnetosheath. The distribution of these data across the $(\beta_{\parallel p},R_p)$-plane reveals limits on both $R_p>1$ and $R_p<1$. Linear Vlasov theory is used to compute contours of constant growth-rate for the ion-cyclotron, mirror, parallel-firehose, and oblique-firehose instabilities. These instability thresholds closely align with the contours of the data distribution, which suggests a strong association of instabilities with extremes of ion temperature anisotropy in the magnetosheath. The potential for instabilities to regulate temperature anisotropy is discussed.

Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).

A prediction of the steady-state reconnection electric field in asymmetric reconnection is obtained by maximizing the reconnection rate as a function of the opening angle made by the upstream magnetic field on the weak magnetic field (magnetosheath) side. The prediction is within a factor of two of the widely examined asymmetric reconnection model [Cassak and Shay, Phys. Plasmas 14, 102114, 2007] in the collisionless limit, and they scale the same over a wide parameter regime. The previous model had the effective aspect ratio of the diffusion region as a free parameter, which simulations and observations suggest is on the order of 0.1, but the present model has no free parameters. In conjunction with the symmetric case [Liu et al., Phys. Rev. Lett. 118, 085101, 2017], this work further suggests that this nearly universal number 0.1, essentially the normalized fast reconnection rate, is a geometrical factor arising from maximizing the reconnection rate within magnetohydrodynamic (MHD)-scale constraints.

A long-standing problem in magnetic reconnection is to explain why it tends to proceed at or below a normalized rate of 0.1. This article gives a review of observational and numerical evidence for this rate and discusses recent theoretical work addressing this problem. Some remaining open questions are summarized.

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.

Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell (PIC) simulations. The degree of electron heating is well correlated with the inflowing Alfv\'en speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $\Delta T_e = 0.033 \,m_i\,c_{Ar}^2$, where $\Delta T_{e}$ is the increase in electron temperature. For the range of simulations performed, the heating shows almost no correlation with inflow total temperature $T_{tot} = T_i + T_e$ or plasma $\beta$. An out-of-plane (guide) magnetic field of similar magnitude to the reconnecting field does not affect the total heating, but it does quench perpendicular heating, with almost all heating being in the parallel direction. These results are qualitatively consistent with a recent statistical survey of electron heating in the dayside magnetopause, which also found that $\Delta T_e$ was proportional to the inflowing Alfv\'en speed. The net electron heating varies very little with distance downstream of the x-line. The simulations show at most a very weak dependence of electron heating on the ion to electron mass ratio. In the antiparallel reconnection case, the largely parallel heating is eventually isotropized downstream due a scattering mechanism such as stochastic particle motion or instabilities. The study highlights key properties that must be satisfied by an electron heating mechanism: (1) Preferential heating in the parallel direction; (2) Heating proportional to $m_i\,c_{Ar}^2$; (3) At most a weak dependence on electron mass; and (4) An exhaust electron temperature that varies little with distance from the x-line.

A model of supra-arcade downflows (SADs), dark low density regions also known as tadpoles that propagate sunward during solar flares, is presented. It is argued that the regions of low density are flow channels carved by sunward-directed outflow jets from reconnection. The solar corona is stratified, so the flare site is populated by a lower density plasma than that in the underlying arcade. As the jets penetrate the arcade, they carve out regions of depleted plasma density which appear as SADs. The present interpretation differs from previous models in that reconnection is localized in space but not in time. Reconnection is continuous in time to explain why SADs are not filled in from behind as they would if they were caused by isolated descending flux tubes or the wakes behind them due to temporally bursty reconnection. Reconnection is localized in space because outflow jets in standard two-dimensional reconnection models expand in the normal (inflow) direction with distance from the reconnection site, which would not produce thin SADs as seen in observations. On the contrary, outflow jets in spatially localized three-dimensional reconnection with an out-of-plane (guide) magnetic field expand primarily in the out-of-plane direction and remain collimated in the normal direction, which is consistent with observed SADs being thin. Two-dimensional proof-of-principle simulations of reconnection with an out-of-plane (guide) magnetic field confirm the creation of SAD-like depletion regions and the necessity of density stratification. Three-dimensional simulations confirm that localized reconnection remains collimated.

We use more than 4,500 microflares from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) microflare data set (Christe et al., 2008, Ap. J., 677, 1385) to estimate electron densities and volumetric filling factors of microflare loops using a cooling time analysis. We show that if the filling factor is assumed to be unity, the calculated conductive cooling times are much shorter than the observed flare decay times, which in turn are much shorter than the calculated radiative cooling times. This is likely unphysical, but the contradic- tion can be resolved by assuming the radiative and conductive cooling times are comparable, which is valid when the flare loop temperature is a maximum and when external heating can be ignored. We find that resultant radiative and con- ductive cooling times are comparable to observed decay times, which has been used as an assumption in some previous studies. The inferred electron densities have a mean value of 10^11.6 cm^-3 and filling factors have a mean of 10^-3.7. The filling factors are lower and densities are higher than previous estimates for large flares, but are similar to those found for two microflares by Moore et al. (Ap. J., 526, 505, 1999).

There is currently no explanation of why the corona has the temperature and density it has. We present a model which explains how the dynamics of magnetic reconnection regulates the conditions in the corona. A bifurcation in magnetic reconnection at a critical state enforces an upper bound on the coronal temperature for a given density. We present observational evidence from 107 flares in 37 sun-like stars that stellar coronae are near this critical state. The model may be important to self-organized criticality models of the solar corona.

We study the transition between anti-parallel and component collisionless magnetic reconnection with 2D particle-in-cell simulations. The primary finding is that a guide field \approx 0.1 times as strong as the asymptotic reconnecting field -- roughly the field strength at which the electron Larmor radius is comparable to the width of the electron current layer -- is sufficient to magnetize the electrons in the vicinity of the x-line, thus causing significant changes to the structure of the electron dissipation region. This implies that great care should be exercised before concluding that magnetospheric reconnection is antiparallel. We also find that even for such weak guide fields strong inward-flowing electron beams form in the vicinity of the magnetic separatrices and Buneman-unstable distribution functions arise at the x-line itself. As in the calculations of {\it Hesse et al.} [2002] and {\it Yin and Winske} [2003], the non-gyrotropic elements of the electron pressure tensor play the dominant role in decoupling the electrons from the magnetic field at the x-line, regardless of the magnitude of the guide field and the associated strong variations in the pressure tensor's spatial structure. Despite these changes, and consistent with previous work, the reconnection rate does not vary appreciably with the strength of the guide field as it changes between 0 and a value equal to the asymptotic reversed field.

The nature of collisionless reconnection in a three-species plasma composed of a heavy species, protons, and electrons is examined. Besides the usual two length scales present in two-species reconnection, there are two additional larger length scales in the system: one associated with a "heavy whistler" which produces a large scale quadrupolar out-of-plane magnetic field, and one associated with the "heavy Alfven" wave which can slow the outflow speed and thus the reconnection rate. The consequences for reconnection in the magnetotail with an O+ population present are discussed.

We present particle-in-cell simulations of collisionless magnetic reconnection in a system (like the magnetopause) with a large density asymmetry across the current layer. In the presence of an ambient component of the magnetic field perpendicular to the reconnection plane the gradient creates a diamagnetic drift that advects the X-line with the electron diamagnetic velocity. When the relative drift between the ions and electrons is of the order the Alfven speed the large scale outflows from the X-line necessary for fast reconnection cannot develop and the reconnection is suppressed. We discuss how these effects vary with both the plasma beta and the shear angle of the reconnecting field and discuss observational evidence for diamagnetic stabilization at the magnetopause.