Sweet-Parker Magnetic Reconnection

Sweet-Parker magnetic reconnection is the first quantitative, self-consistent magnetohydrodynamic (MHD) model of magnetic reconnection—the topological rearrangement of oppositely directed magnetic field lines in a conducting plasma that converts stored magnetic energy into plasma bulk flow kinetic energy, thermal heating, and non-thermal particle acceleration. Independently proposed by British astrophysicist Peter Alan Sweet at a 1956 International Astronomical Union symposium in Stockholm (published 1958) and formulated mathematically by American astrophysicist Eugene N. Parker in 1957, the Sweet-Parker model analyzes steady-state resistive reconnection across an elongated two-dimensional current sheet of macroscopic half-length L and microscopic half-thickness δ. By coupling mass conservation, magnetic tension acceleration, and resistive Ohm's law, Parker demonstrated that the dimensionless reconnection rate M_A = v_in / v_A scales inversely with the square root of the dimensionless Lundquist number (M_A = S^−1/2). Within astronomy and plasma physics, the Sweet-Parker scaling law established the benchmark of resistive MHD while exposing the fast reconnection problem: because S ≈ 10¹² to 10¹⁴ in the solar corona due to low classical Spitzer resistivity, pure Sweet-Parker reconnection predicts solar flare timescales of weeks to years rather than the observed 100 to 1,000 seconds, motivating the discovery of Petschek magnetic reconnection, collisionless Hall reconnection, and the tearing/plasmoid instability.

In ideal magnetohydrodynamics (where electrical resistivity η = 0), Alfvén's theorem of flux freezing (1942) dictates that magnetic field lines are permanently frozen into the conducting plasma fluid: two plasma parcels initially connected by a magnetic field line remain connected by that field line for all time, and distinct magnetic flux tubes can never interpenetrate or change their topological connectivity. If global magnetic diffusion acts across a macroscopic length scale L with magnetic diffusivity η_m = η / μ₀, the resistive diffusion timescale is τ_diff = L² / η_m. Compared to the ideal Alfvén transit time τ_A = L / v_A—where v_A = B₀ / (μ₀ ρ)^1/2 is the Alfvén wave speed—the ratio of these two timescales defines the dimensionless Lundquist number S:

S = τ_diff / τ_A = (L v_A) / η_m = (μ₀ L v_A) / η

In astrophysical and space plasmas, S is enormous: S ≈ 10⁸ in laboratory tokamaks, S ≈ 10¹¹ in Earth's magnetotail, S ≈ 10¹²–10¹⁴ in solar coronal loops, and S > 10¹⁸ in the interstellar medium.

Sweet (1958) and Parker (1957) realized that when two oppositely magnetized plasma flux systems (+B₀x̂ and −B₀x̂) are pressed together over a macroscopic system length 2L, the magnetic field gradient ∇ × B is compressed into a thin boundary layer (current sheet) of half-thickness δ ≪ L. Inside this thin diffusion region of thickness 2δ, the spatial second derivative ∇²B ~ B₀ / δ² becomes large enough that resistive diffusion η_m∇²B breaks Alfvén's frozen-in condition, allowing anti-parallel field lines to sever, cross-connect into sharp hairpin bends, and snap sideways out of the ends of the current sheet under magnetic tension.

Consider a steady-state (∂/∂t = 0), two-dimensional, incompressible plasma of mass density ρ entering a planar current sheet of half-length L (along x) and half-thickness δ (along y) from above and below with inflow velocity v_in and upstream magnetic field B₀, and exiting horizontally from the two narrow ends of the sheet with outflow velocity v_out. Parker (1957) derived the reconnection rate from three conservation laws:

1. Conservation of Mass (Continuity Equation)

In steady state, the mass flux entering the two wide faces of the diffusion layer (each of area 2L per unit depth) must equal the mass flux ejected through the two narrow ends (each of area 2δ per unit depth). For an approximately incompressible density ρ:

ρ v_in L = ρ v_out δ ⇒ v_in / v_out = δ / L

Because the sheet is thin (δ ≪ L) while the inflow must pass through the narrow throat of width δ to escape, this geometric continuity bottleneck restricts the inflow velocity v_in unless the outflow velocity v_out is large and δ/L is not too small.

2. Conservation of Momentum (Alfvénic Outflow Speed)

Inside the current sheet, newly reconnected field lines form sharp bends with radius of curvature ~δ that exert an outward magnetic tension force (J × B)_x ≈ B₀² / (μ₀L) (or equivalently, a magnetic pressure drop ΔP = B₀² / (2μ₀) between the stagnation center of the sheet and the external upstream field). Applying Bernoulli's equation (or the steady-state MHD momentum equation ρ(v · ∇)v = −∇P + J × B) along the outflow channel from the stagnation point (x = 0, v = 0) to the exit nozzle (x = ±L, v = v_out) yields:

(1 / 2) ρ v_out² ≈ B₀² / (2 μ₀) ⇒ v_out ≈ B₀ / (μ₀ ρ)^1/2 = v_A

Thus, regardless of the resistivity η_m, the plasma is slung out of the two ends of the Sweet-Parker current sheet at the upstream Alfvén speed v_A, converting roughly half of the incoming magnetic energy flux into bulk flow kinetic energy of the Alfvénic outflow jets and the other half into Ohmic thermal heating inside the layer.

3. Resistive Ohm's Law (Advection-Diffusion Balance)

In steady-state 2D reconnection, Faraday's law ∇ × E = −∂B/∂t = 0 requires that the out-of-plane reconnection electric field E_z be spatially uniform across both the ideal upstream inflow zone (where J ≈ 0 and E_z = v_inB₀) and the center of the resistive current sheet (where v = 0 and E_z = ηJ_z ≈ (η/μ₀)(B₀/δ) = η_mB₀/δ):

v_in B₀ = η_m ( B₀ / δ ) ⇒ δ = η_m / v_in

4. The Sweet-Parker Reconnection Rate (M_A = S^−1/2)

Substituting v_out = v_A and δ = η_m / v_in into the mass continuity relation v_in / v_A = δ / L gives:

v_in / v_A = η_m / (L v_in) ⇒ ( v_in / v_A )² = η_m / (L v_A) = 1 / S

Taking the square root yields the celebrated Sweet-Parker scaling relations for the dimensionless reconnection rate (Alfvén Mach number M_A), current-sheet aspect ratio δ/L, and reconnection timescale τ_SP = L / v_in:

• Dimensionless Reconnection Rate: M_A = v_in / v_A = S^−1/2
• Current-Sheet Aspect Ratio: δ / L = S^−1/2
• Sweet-Parker Reconnection Time: τ_SP = L / v_in = S^1/2 τ_A = ( τ_A τ_diff )^1/2 (the geometric mean of the Alfvén time and the global resistive diffusion time).

Why Sweet-Parker Reconnection Is Too Slow for Solar Flares

The table above illustrates why Parker (1963) recognized a fundamental discrepancy when applying the S^−1/2 law to solar flares. In an active solar coronal loop (L ≈ 10⁷ m, B₀ ≈ 100 G = 0.01 T, n_e ≈ 10¹⁵ m⁻³, T_e ≈ 2 × 10⁶ K), the Alfvén speed is v_A ≈ 10⁶ m s⁻¹, so the Alfvén crossing time is τ_A = L / v_A ≈ 10 seconds, and the Lundquist number based on classical Spitzer resistivity is S ≈ 10¹² to 10¹⁴.

Under Sweet-Parker scaling, the current sheet must thin to δ = L S^−1/2 ≈ 1 to 10 meters while remaining 20,000 kilometers long, choking the inflow velocity to v_in = v_A S^−1/2 ≈ 0.1 to 1 m s⁻¹ (M_A ≈ 10⁻⁷ to 10⁻⁶). Exhausting the coronal magnetic flux would therefore require τ_SP = S^1/2 τ_A ≈ 10⁷ to 10⁸ seconds (4 months to 3 years). Yet X-ray, ultraviolet, and Hα observations show that impulsive solar flares release 10²⁴ to 10²⁵ Joules of magnetic energy in only 100 to 1,000 seconds (10 to 100 τ_A), corresponding to an empirical dimensionless reconnection rate of M_A ≈ 0.01 to 0.1—five orders of magnitude faster than Sweet-Parker theory predicts.

Three Physical Pathways Beyond the Sweet-Parker Bottleneck

The physical reason Sweet-Parker reconnection is slow is geometric: the model forces the diffusion layer length L to equal the macroscopic system size, requiring every incoming fluid parcel to squeeze out through two tiny nozzles of width δ = L S^−1/2. Modern plasma physics has identified three mechanisms that break this geometric bottleneck at high Lundquist numbers:

1. Petschek Magnetic Reconnection (1964) and Anomalous Resistivity: Harry Petschek showed that if the diffusion region is localized to a microscopic X-line of length L_e ≪ L while four standing slow-mode MHD shock waves deflect and accelerate the bulk plasma outside the diffusion region, the reconnection rate rises to a logarithmic scaling M_A ≈ π / (8 ln S) ≈ 0.01–0.05 (stabilized when anomalous resistivity is spatially localized at the X-point).
2. The Plasmoid (Loureiro-Schekochihin-Cowley) Instability (S > S_c ≈ 10⁴): In 2007, Nuno Loureiro, Alexander Schekochihin, and Steven Cowley proved that an elongated Sweet-Parker current sheet with aspect ratio L / δ = S^1/2 becomes dynamically unstable to a super-Alfvénic secondary tearing instability (growth rate γ τ_A ∝ S^1/4) whenever the Lundquist number exceeds a threshold S_c ≈ 10⁴. Above S_c, a monolithic Sweet-Parker sheet cannot exist in nature: it spontaneously fragments into a dynamic, hierarchical chain of magnetic islands (plasmoids or flux ropes) separated by secondary current sheets of local threshold length, saturating the time-averaged resistive MHD reconnection rate at a plateau of M_A ≈ S_c^−1/2 ≈ 0.01 independent of S (Bhattacharjee et al., 2009; Uzdensky et al., 2010).
3. Collisionless Hall-MHD and Kinetic Reconnection (δ ≤ d_i): When the Sweet-Parker thickness δ = L S^−1/2 shrinks below the ion skin depth (ion inertial length) d_i = c / ω_pi = (m_i / (μ₀ n_i e²))^1/2, ions decouple from the magnetic field while electrons remain magnetized down to the electron skin depth d_e = c / ω_pe. Quadrupolar Hall magnetic fields and dispersive kinetic Alfvén / whistler waves open the outflow exhaust cone, yielding universal fast collisionless reconnection at M_A ≈ 0.1 (Birn et al., 2001).

What is the physical reason why Sweet-Parker reconnection becomes slower as electrical conductivity increases?

In a steady-state Sweet-Parker current sheet, the inflowing plasma carries magnetic flux into the layer at rate v_inB₀, which must match resistive magnetic diffusion η_mB₀/δ across the layer thickness δ. As plasma conductivity increases (η_m decreases, S increases), the current sheet must become proportionally thinner (δ = η_m/v_in) to maintain resistive dissipation. However, because the outflow speed at the ends of the sheet is capped at the Alfvén speed v_A, a thinner exit nozzle δ can vent less mass per second (v_inL = v_Aδ), throttling the inflow rate to v_in/v_A = S^−1/2.

Has Sweet-Parker reconnection been verified in laboratory plasma experiments?

Yes. Dedicated laboratory experiments such as the Magnetic Reconnection Experiment (MRX) at the Princeton Plasma Physics Laboratory (Yamada et al., 1997, 2010) and the Versatile Toroidal Facility (VTF) at MIT directly measured steady-state reconnection rates in collisional plasmas where the Lundquist number is moderate (S ≈ 50 to 1,000) and the Sweet-Parker thickness δ exceeds the ion skin depth d_i. In this collisional regime, the measured reconnection rate closely follows the classical M_A = S^−1/2 scaling law before transitioning to fast reconnection (M_A ≈ 0.1) when δ drops below d_i.

What happens to a Sweet-Parker current sheet when the Lundquist number exceeds S_c ≈ 10⁴?

Above S_c ≈ 10⁴, the extreme aspect ratio of the sheet (L / δ > 100) triggers the plasmoid instability: small magnetic perturbations grow faster than the plasma can be ejected from the ends of the sheet (γ τ_A ~ S^1/4 ≫ 1). The elongated sheet breaks up into a fractal chain of magnetic islands (plasmoids) observed as bright moving blobs in solar flare current sheets and Earth's magnetotail, maintaining a fast reconnection rate of M_A ≈ 0.01 even in purely resistive MHD.

1. Parker, E. N. (1957). "Sweet's mechanism for merging magnetic fields in conducting fluids." Journal of Geophysical Research, 62(4), 509–520. https://doi.org/10.1029/JZ062i004p00509
2. Sweet, P. A. (1958). "The neutral point theory of solar flares." In B. Lehnert (Ed.), Electromagnetic Phenomena in Cosmical Physics (IAU Symposium No. 6, pp. 123–134). Cambridge University Press. https://ui.adsabs.harvard.edu/abs/1958IAUS....6..123S
3. Loureiro, N. F., Schekochihin, A. A., & Cowley, S. C. (2007). "Instability of current sheets and formation of plasmoid chains." Physics of Plasmas, 14(10), 100703. https://doi.org/10.1063/1.2783986
4. Yamada, M., Kulsrud, R., & Ji, H. (2010). "Magnetic reconnection." Reviews of Modern Physics, 82(1), 603–664. https://doi.org/10.1103/RevModPhys.82.603