Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically-dominated ($\sigma_0\gg1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim \sigma_0 \gamma_{th0}$ (here, $\gamma_{th0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are roughly isotropic. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu F_\nu\propto \nu^s$ that is much harder than for isotropic particles. Remarkably, for $\sigma_0\gg1$ we find that the solid-angle-averaged spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta B_{\rm rms0}^2/B_0^2\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Our findings may help explaining the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of Pulsar Wind Nebulae.

We construct a turbulent model of the Crab Nebula's non-thermal emission. The present model resolves a number of long-standing problems of the Kennel-Coroniti (1984) model: (i) the sigma problem; (ii) the hard spectrum of radio electrons; (iii) the high peak energy of gamma-ray flares; (iv) and the spacial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I accelerated at the wind termination shock via Fermi-I mechanism, and Component-II accelerated in reconnecting turbulence in highly magnetized ($\sigma$ $\gg 1$) plasma in the central part of the Crab Nebula. The reconnecting turbulence Component-II extends from radio to gamma rays: it accelerate radio electrons with a hard spectrum, destroy the large scale magnetic flux (and thus resolves the sigma-problem), and occasionally produces gamma-ray flares (from the largest scale reconnection events). The model reproduces the broad-band spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as spatially resolved evolution of the spectral indices in IR and optical bands.

The ideal magnetohydrodynamic theorem on the conservation of the magnetic connections between plasma elements is extended to non-ideal relativistic plasmas in curved spacetime. The existence of generalized magnetofluid connections that are preserved by the plasma dynamics is formalized by means of a covariant connection equation that includes different non-ideal effects. These generalized connections are constituted by 2-dimensional hypersurfaces, which are linked to an antisymmetric tensor field that unifies the electromagnetic and fluid fields. They can be interpreted in terms of generalized magnetofluid vorticity field lines by considering a 3+1 foliation of spacetime and a time resetting projection that compensates for the loss of simultaneity in different reference frames between spatially separated events. The worldsheets of the generalized magnetofluid vorticity field lines play a fundamental role in the plasma dynamics by prohibiting evolutions that do not preserve the magnetofluid connectivity.

Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level of fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power-law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. The analytic scaling relations are validated with numerical solutions of the model.

Magnetized turbulence and magnetic reconnection are often invoked to explain the nonthermal emission observed from a wide variety of astrophysical sources. By means of fully-kinetic 2D and 3D PIC simulations, we investigate the interplay between turbulence and reconnection in generating nonthermal particles in magnetically-dominated pair plasmas. A generic by-product of the turbulence evolution is the generation of a nonthermal particle spectrum with a power-law energy range. The power-law slope $p$ is harder for larger magnetizations and stronger turbulence fluctuations, and it can be as hard as $p < 2$. The Larmor radius of particles at the high-energy cutoff is comparable to the size $l$ of the largest eddies. Plasmoid-mediated reconnection, which self-consistently occurs in the turbulent plasma, controls the physics of particle injection. Then, particles are further accelerated by stochastic scattering off turbulent fluctuations. The work done by parallel electric fields - naturally expected in reconnection layers - is responsible for most of the initial energy increase, and is proportional to the magnetization $\sigma$ of the system, while the subsequent energy gain, which dominates the overall energization of high-energy particles, is powered by the perpendicular electric fields of turbulent fluctuations. The two-stage acceleration process leaves an imprint in the particle pitch-angle distribution: low-energy particles are aligned with the field, while the highest energy particles move preferentially orthogonal to it. The energy diffusion coefficient of stochastic acceleration scales as $D_\gamma\sim 0.1\sigma(c/l)\gamma^2$, where $\gamma$ is the particle Lorentz factor. This results in fast acceleration timescales $t_{acc}\sim (3/\sigma)\,l/c$. Our findings have important implications for understanding the generation of nonthermal particles in high-energy astrophysical sources.

We outline a model of the Crab Pulsar Wind Nebula with two different populations of synchrotron emitting particles, arising from two different acceleration mechanisms: (i) Component-I due to Fermi-I acceleration at the equatorial portion of the termination shock, with particle spectral index $p_I \approx 2.2$ above the injection break corresponding to $\gamma_{wind} \sigma_{wind} \sim 10^5$, peaking in the UV ($\gamma_{wind} \sim 10^2$ is the bulk Lorentz factor of the wind, $\sigma _{wind} \sim 10^3$ is wind magnetization); (ii) Component-II due to acceleration at reconnection layers in the bulk of the turbulent Nebula, with particle index $p_{II} \approx 1.6$. The model requires relatively slow but highly magnetized wind. For both components the overall cooling break is in the infra-red at $\sim 0.01$ eV, so that the Component-I is in the fast cooling regime (cooling frequency below the peak frequency). In the optical band Component-I produces emission with the cooling spectral index of $\alpha_o \approx 0.5$, softening towards the edges due to radiative losses. Above the cooling break, in the optical, UV and X-rays, Component-I mostly overwhelms Component-II. We hypothesize that acceleration at large-scale current sheets in the turbulent nebula (Component-II) extends to the synchrotron burn-off limit of $\epsilon_s \approx 100$ MeV. Thus in our model acceleration in turbulent reconnection (Component-II) can produce both hard radio spectra and occasional gamma-ray flares. This model may be applicable to a broader class of high energy astrophysical objects, like AGNe and GRB jets, where often radio electrons form a different population from the high energy electrons.

We analytically explore the effects of the gravitational electromotive force on magnetic reconnection around Schwarzschild black holes through a generalized general-relativistic magnetohydrodynamic model that retains two-fluid effects. It is shown that the gravitational electromotive force can couple to collisionless two-fluid effects and drive magnetic reconnection. This is allowed by the departure from quasi-neutrality in curved spacetime, which is explicitly manifested as the emergence of an effective resistivity in Ohm's law. The departure from quasi-neutrality is owed to different gravitational pulls experienced by separate parts of the current layer. This produces an enhancement of the reconnecion rate due to purely gravitational effects.

The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At sufficiently large magnetic Reynolds number ($R_m=10^6$), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index $-2.2$ as proposed by recent analytic theories. We also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25.

Due to its ubiquitous presence, turbulence is often invoked to explain the origin of nonthermal particles in astrophysical sources of high-energy emission. With particle-in-cell simulations, we study decaying turbulence in magnetically-dominated (or equivalently, "relativistic") pair plasmas. We find that the generation of a power-law particle energy spectrum is a generic by-product of relativistic turbulence. The power-law slope is harder for higher magnetizations and stronger turbulence levels. In large systems, the slope attains an asymptotic, system-size-independent value, while the high-energy spectral cutoff increases linearly with system size; both the slope and the cutoff do not depend on the dimensionality of our domain. By following a large sample of particles, we show that particle injection happens at reconnecting current sheets; the injected particles are then further accelerated by stochastic interactions with turbulent fluctuations. Our results have important implications for the origin of non-thermal particles in high-energy astrophysical sources.

Magnetohydrodynamic turbulence and magnetic reconnection are ubiquitous in astrophysical environments. In most situations, these processes do not occur in isolation, but interact with each other. This renders a comprehensive theory of these processes highly challenging. Here, we propose a theory of magnetohydrodynamic turbulence driven at large scale that self-consistently accounts for the mutual interplay with magnetic reconnection occurring at smaller scales. Magnetic reconnection produces plasmoids that grow from turbulence-generated noise and eventually disrupt the sheet-like structures in which they are born. The disruption of these structures leads to a modification of the turbulent energy cascade, which, in turn, exerts a feedback effect on the plasmoid formation via the turbulence-generated noise. The energy spectrum in this plasmoid-mediated range steepens relative to the standard inertial range and does not follow a simple power law. As a result of the complex interplay between turbulence and reconnection, we also find that the length scale which marks the beginning of the plasmoid-mediated range and the dissipation length scale do not obey true power laws. The transitional magnetic Reynolds number above which the plasmoid formation becomes statistically significant enough to affect the turbulent cascade is fairly modest, implying that plasmoids are expected to modify the turbulent path to dissipation in many astrophysical systems.

Magnetic reconnection in curved spacetime is studied by adopting a general relativistic magnetohydrodynamic model that retains collisionless effects for both electron-ion and pair plasmas. A simple generalization of the standard Sweet-Parker model allows us to obtain the first order effects of the gravitational field of a rotating black hole. It is shown that the black hole rotation acts as to increase the length of azimuthal reconnection layers, per se leading to a decrease of the reconnection rate. However, when coupled to collisionless thermal-inertial effects, the net reconnection rate is enhanced with respect to what would happen in a purely collisional plasma due to a broadening of the reconnection layer. These findings identify an underlying interaction between gravity and collisionless magnetic reconnection in the vicinity of compact objects.

The Principle of Maximum Entropy, a powerful and general method for inferring the distribution function given a set of constraints, is applied to deduce the overall distribution of 3D plasmoids (flux ropes/tubes) for systems where resistive MHD is applicable and large numbers of plasmoids are produced. The analysis is undertaken for the 3D case, with mass, total flux and velocity serving as the variables of interest, on account of their physical and observational relevance. The distribution functions for the mass, width, total flux and helicity exhibit a power-law behavior with exponents of $-4/3$, $-2$, $-3$ and $-2$ respectively for small values, whilst all of them display an exponential falloff for large values. In contrast, the velocity distribution, as a function of $v = |{\bf v}|$, is shown to be flat for $v \rightarrow 0$, and becomes a power law with an exponent of $-7/3$ for $v \rightarrow \infty$. Most of these results are nearly independent of the free parameters involved in this specific problem. A preliminary comparison of our results with the observational evidence is presented, and some of the ensuing space and astrophysical implications are briefly discussed.

The ideal MHD theorem on the conservation of the magnetic connections between plasma elements is generalized to relativistic plasmas in curved spacetime. The connections between plasma elements, which are established by a covariant connection equation, display a particularly complex structure in curved spacetime. Nevertheless, it is shown that these connections can be interpreted in terms of magnetic field lines alone by adopting a $3+1$ foliation of spacetime.

The scaling of plasmoid instability maximum linear growth rate with respect to Lundquist number $S$ in a Sweet-Parker current sheet, $\gamma_{max}\sim S^{1/4}$, indicates that at high $S$, the current sheet will break apart before it approaches the Sweet-Parker width. Therefore, a proper description for the onset of the plasmoid instability must incorporate the evolving process of the current sheet. We carry out a series of two-dimensional simulations and develop diagnostics to separate fluctuations from an evolving background. It is found that the fluctuation amplitude starts to grow only when the linear growth rate is sufficiently large ($\gamma_{max}\tau_{A}>O(1)$) to overcome convective losses. The linear growth rate continues to rise until the sizes of plasmoids become comparable to the inner layer width of the tearing mode. At this point the current sheet is disrupted and the instability enters the early nonlinear regime. The growth rate suddenly decreases, but the fluctuation amplitude continues to grow until it reaches nonlinear saturation. We identify important time scales of the instability development, as well as scalings for the linear growth rate, current sheet width, and dominant wavenumber at current sheet disruption. These scalings depend on not only the Lundquist number, but also the initial noise amplitude. A phenomenological model that reproduces scalings from simulation results is proposed. The model incorporates the effect of reconnection outflow, which is crucial for yielding a critical Lundquist number $S_{c}$ below which disruption does not occur. The critical Lundquist number $S_{c}$ is not a constant value but has a weak dependence on the noise amplitude.

The plasmoid instability has revolutionized our understanding of magnetic reconnection in astrophysical environments. By preventing the formation of highly elongated reconnection layers, it is crucial in enabling the rapid energy conversion rates that are characteristic of many astrophysical phenomena. Most of the previous studies have focused on Sweet-Parker current sheets, which, however, are unattainable in typical astrophysical systems. Here, we derive a general set of scaling laws for the plasmoid instability in resistive and visco-resistive current sheets that evolve over time. Our method relies on a principle of least time that enables us to determine the properties of the reconnecting current sheet (aspect ratio and elapsed time) and the plasmoid instability (growth rate, wavenumber, inner layer width) at the end of the linear phase. After this phase the reconnecting current sheet is disrupted and fast reconnection can occur. The scaling laws of the plasmoid instability are \emph{not} simple power laws, and depend on the Lundquist number ($S$), the magnetic Prandtl number ($P_m$), the noise of the system ($\psi_0$), the characteristic rate of current sheet evolution ($1/\tau$), as well as the thinning process. We also demonstrate that previous scalings are inapplicable to the vast majority of the astrophysical systems. We explore the implications of the new scaling relations in astrophysical systems such as the solar corona and the interstellar medium. In both these systems, we show that our scaling laws yield values for the growth rate, wavenumber, and aspect ratio that are much smaller than the Sweet-Parker based scalings.

A new formulation of the plasma resistivity that arises from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The exact expression for the collisional momentum-transfer rate is determined, and is used to formulate the nonlinear anisotropic resistivity. The new formalism recovers the Spitzer resistivity, as well as the concept of thermal force if the heat flux is assumed to be proportional to a temperature gradient. Furthermore, if the pressure tensor is related to viscous stress, the latter enters the expression for the resistivity. The relative importance of the nonlinear term(s) with respect to the well-established electron inertia and Hall terms is also examined. The subtle implications of the nonlinear resistivity, and its dependence on the fluid variables, are discussed in the context of magnetized plasma environments and phenomena such as magnetic reconnection.

A general theory of the onset and development of the plasmoid instability is formulated by means of a principle of least time. The scaling relations for the final aspect ratio, transition time to rapid onset, growth rate, and number of plasmoids are derived, and shown to depend on the initial perturbation amplitude $\left({\hat w}_0\right)$, the characteristic rate of current sheet evolution $\left(1/\tau\right)$, and the Lundquist number $\left(S\right)$. They are not simple power laws, and are proportional to $S^{\alpha} \tau^{\beta} \left[\ln f(S,\tau,{\hat w}_0)\right]^\sigma$. The detailed dynamics of the instability is also elucidated, and shown to comprise of a period of quiescence followed by sudden growth over a short time scale.

The magnetic reconnection process is analyzed for relativistic magnetohydrodynamical plasmas around rotating black holes. A simple generalization of the Sweet-Parker model is used as a first approximation to the problem. The reconnection rate, as well as other important properties of the reconnection layer, have been calculated taking into account the effect of spacetime curvature. Azimuthal and radial current sheet configurations in the equatorial plane of the black hole have been studied, and the case of small black hole rotation rate has been analyzed. For the azimuthal configuration, it is found that the black hole rotation decreases the reconnection rate. On the other hand, in the radial configuration, it is the gravitational force created by the black hole mass that decreases the reconnection rate. These results establish a fundamental interaction between gravity and magnetic reconnection in astrophysical contexts.

Numerical simulations have consistently shown that the reconnection rate in certain collisionless regimes can be fast, on the order of $0.1 v_A B_u$, where $v_A$ and $B_u$ are the Alfv\'en speed and the reconnecting magnetic field upstream of the ion diffusion region. This particular value has been reported in myriad numerical simulations under disparate conditions. However, despite decades of research, the reasons underpinning this specific value remain mysterious. Here, we present an overview of this problem and discuss the conditions under which the "0.1 value" is attained. Furthermore, we explain why this problem should be interpreted in terms of the ion diffusion region length.

The plasmoid instability in visco-resistive current sheets is analyzed in both the linear and nonlinear regimes. The linear growth rate and the wavenumber are found to scale as $S^{1/4} {\left( {1 + {P_m}} \right)}^{-5/8}$ and $S^{3/8} {\left( {1 + {P_m}} \right)}^{-3/16}$ with respect to the Lundquist number $S$ and the magnetic Prandtl number $P_m$. Furthermore, the linear layer width is shown to scale as $S^{-1/8} {(1+P_m)}^{1/16}$. The growth of the plasmoids slows down from an exponential growth to an algebraic growth when they enter into the nonlinear regime. In particular, the time-scale of the nonlinear growth of the plasmoids is found to be $\tau_{NL} \sim S^{-3/16} {(1 + P_m)^{19/32}}{\tau_{A,L}}$. The nonlinear growth of the plasmoids is radically different from the linear one and it is shown to be essential to understand the global current sheet disruption. It is also discussed how the plasmoid instability enables fast magnetic reconnection in visco-resistive plasmas. In particular, it is shown that the recursive plasmoid formation can trigger a collisionless reconnection regime if $S \gtrsim L_{cs} {(\epsilon_c l_k)^{-1}} {(1 + {P_m})^{1/2}}$, where $L_{cs}$ is the half-length of the global current sheet and $l_k$ is the relevant kinetic length scale. On the other hand, if the current sheet remains in the collisional regime, the global (time-averaged) reconnection rate is shown to be $\left\langle {{{\left. {d\psi /dt} \right|}_X}} \right\rangle \approx \epsilon_c v_{A,u} B_{u} {(1 + {P_m})^{-1/2}}$, where $\epsilon_c$ is the critical inverse aspect ratio of the current sheet, while $v_{A,u}$ and $B_{u}$ are the Alfv\'en speed and the magnetic field upstream of the global reconnection layer.

We extend the magnetic connection theorem of ideal magnetohydrodynamics to nonideal relativistic pair plasmas. Adopting a generalized Ohm's law, we prove the existence of generalized magnetofluid connections that are preserved by the plasma dynamics. We show that these connections are related to a general antisymmetric tensor that unifies the electromagnetic and fluid fields. The generalized magnetofluid connections set important constraints on the plasma dynamics by forbidding transitions between configurations with different magnetofluid connectivity. An approximated solution is explicitly shown where the corrections due to current inertial effects are found.

Recent progress in the understanding of how externally driven magnetic reconnection evolves is organized in terms of parameter space diagrams. These diagrams are constructed using four pivotal dimensionless parameters: the Lundquist number $S$, the magnetic Prandtl number $P_m$, the amplitude of the boundary perturbation $\hat \Psi_0$, and the perturbation wave number $\hat k$. This new representation highlights the parameters regions of a given system in which the magnetic reconnection process is expected to be distinguished by a specific evolution. Contrary to previously proposed phase diagrams, the diagrams introduced here take into account the dynamical evolution of the reconnection process and are able to predict slow or fast reconnection regimes for the same values of $S$ and $P_m$, depending on the parameters that characterize the external drive, never considered so far. These features are important to understand the onset and evolution of magnetic reconnection in diverse physical systems

A fundamental problem of forced magnetic reconnection has been solved taking into account the plasmoid instability of thin reconnecting current sheets. In this problem, the reconnection is driven by a small amplitude boundary perturbation in a tearing-stable slab plasma equilibrium. It is shown that the evolution of the magnetic reconnection process depends on the external source perturbation and the microscopic plasma parameters. Small perturbations lead to a slow nonlinear Rutherford evolution, whereas larger perturbations can lead to either a stable Sweet-Parker-like phase or a plasmoid phase. An expression for the threshold perturbation amplitude required to trigger the plasmoid phase is derived, as well as an analytical expression for the reconnection rate in the plasmoid-dominated regime. Visco-resistive magnetohydrodynamic simulations complement the analytical calculations. The plasmoid formation plays a crucial role in allowing fast reconnection in a magnetohydrodynamical plasma, and the presented results suggest that it may occur and have profound consequences even if the plasma is tearing-stable.

The concept of magnetic connections is extended to non-ideal relativistic magnetohydrodynamical plasmas. Adopting a general set of equations for relativistic magnetohydrodynamics including thermal-inertial, thermal electromotive, Hall and current-inertia effects, we derive a new covariant connection equation showing the existence of generalized magnetofluid connections that are preserved during the dissipationless plasma dynamics. These connections are intimately linked to a general antisymmetric tensor that unifies the electromagnetic and fluid fields, allowing the extension of the magnetic connection notion to a much broader concept.

The magnetic reconnection process is studied in relativistic pair plasmas when the thermal and inertial properties of the magnetohydrodynamical fluid are included. We find that in both Sweet-Parker and Petschek relativistic scenarios there is an increase of the reconnection rate owing to the thermal-inertial effects, both satisfying causality. To characterize the new effects we define a thermal-inertial number which is independent of the relativistic Lundquist number, implying that reconnection can be achieved even for vanishing resistivity as a result of only thermal-inertial effects. The current model has fundamental importance for relativistic collisionless reconnection, as it constitutes the simplest way to get reconnection rates faster than those accessible with the sole resistivity.

The first detailed comparison between gyrokinetic and gyrofluid simulations of collisionless magnetic reconnection has been carried out. Both the linear and nonlinear evolution of the collisionless tearing mode have been analyzed. In the linear regime, we have found a good agreement between the two approaches over the whole spectrum of linearly unstable wave numbers, both in the drift kinetic limit and for finite ion temperature. Nonlinearly, focusing on the small-$\Delta '$ regime, with $\Delta '$ indicating the standard tearing stability parameter, we have compared relevant observables such as the evolution and saturation of the island width, as well as the island oscillation frequency in the saturated phase.The results are basically the same, with small discrepancies only in the value of the saturated island width for moderately high values of $\Delta '$. Therefore, in the regimes investigated here, the gyrofluid approach can describe the collisionless reconnection process as well as the more complete gyrokinetic model.

The linear and nonlinear evolution of magnetic reconnection in collisionless high-temperature plasmas with a strong guide field is analyzed on the basis of a two-dimensional gyrofluid model. The linear growth rate of the reconnecting instability is compared to analytical calculations over the whole spectrum of linearly unstable wave numbers. In the strongly unstable regime (large \Delta '), the nonlinear evolution of the reconnecting instability is found to undergo two distinctive acceleration phases separated by a stall phase in which the instantaneous growth rate decreases. The first acceleration phase is caused by the formation of strong electric fields close to the X-point due to ion gyration, while the second acceleration phase is driven by the development of an open Petschek-like configuration due to both ion and electron temperature effects. Furthermore, the maximum instantaneous growth rate is found to increase dramatically over its linear value for decreasing diffusion layers. This is a consequence of the fact that the peak instantaneous growth rate becomes weakly dependent on the microscopic plasma parameters if the diffusion region thickness is sufficiently smaller than the equilibrium magnetic field scale length. When this condition is satisfied, the peak reconnection rate asymptotes to a constant value.