Plasmoid Cascade in Magnetic Reconnection

• Plasmoid cascade is the multiscale fragmentation of current sheets into magnetic islands and flux ropes, combining forward tearing with inverse coalescence.
• It is governed by critical parameters like the Lundquist number and aspect ratios, which set distinct onset thresholds across different plasma regimes.
• Nonlinear evolution in plasmoid cascades features chain formation, merger dynamics, and turbulent energy transfer, as evidenced in solar, laboratory, and astrophysical studies.

Plasmoid cascade denotes the multiscale disruption of an elongated reconnecting current sheet into magnetic islands or flux ropes and the inter-plasmoid current sheets between them. In the strongest usage, the term refers to recursive current-sheet fragmentation across successive scales; in a broader and common usage, it also includes finite plasmoid chains, coalescence-driven restructuring, and reconnection regimes in which plasmoids control the nonlinear evolution of the sheet. The literature is therefore not uniform: some studies explicitly show secondary fragmentation and tertiary sheets, whereas others show only a short plasmoid chain, or a merger-dominated state without a fully self-similar hierarchy (Comisso et al., 2014, Karlický et al., 2011, Ali et al., 2019, Markidis et al., 2013).

1. Conceptual scope and relation to current-sheet reconnection

In resistive MHD, the canonical precursor of a plasmoid cascade is a Sweet–Parker-like current sheet that becomes sufficiently long and thin to tear into multiple islands. In forced reconnection, this mechanism can replace slower nonlinear pathways. In the Taylor problem, for very small plasma resistivity and viscosity, the linear inertial phase is followed by a nonlinear Sweet–Parker evolution, which gives way to a faster reconnection regime characterized by a chain of plasmoids instead of a slower Rutherford phase (Comisso et al., 2014). The same article treats the plasmoid chain as a set of shorter, marginally stable interplasmoid current sheets, rather than a single global resistive bottleneck.

Solar-flare literature introduced an explicitly bidirectional picture. In a standard flare current sheet, stretching-tearing fragmentation produces plasmoids and thinner sub-sheets, while successive merging of plasmoids produces larger plasmoids; the sheet between two merging plasmoids can itself fragment into smaller and smaller plasmoids and current sheets (Karlický et al., 2011). This formulation is close to the older language of “fractal reconnection,” but it is more specific about the coupled roles of tearing and coalescence.

A central terminological distinction follows from these studies. “Plasmoid chain” is appropriate for a finite sequence of islands along a sheet; “plasmoid cascade” is most precise when secondary breakup of interplasmoid sheets is demonstrated. Several later papers adopt the broader usage because real systems often show mixed forward fragmentation and inverse coalescence rather than a clean scale-by-scale hierarchy (Markidis et al., 2012, Ali et al., 2019).

2. Onset criteria, thresholds, and control parameters

The standard macroscopic control parameter is the Lundquist number. In semi-collisional theory, the relevant ordering is

$\delta_{SP} \gg \rho_s \gg \delta_{in},$

with

$$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$

For a sinusoidal-like current sheet, the semi-collisional plasmoid window is

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$

while for a Harris-type sheet the lower bound becomes

$$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$

In reduced-gyrokinetic simulations, plasmoid instability was obtained at values of $S$ as low as $\sim 250$, and the study emphasized that the current sheet was formed self-consistently rather than prescribed (Bhat et al., 2018).

Tokamak sawtooth simulations show that the threshold depends strongly on geometry and drive. In a cylindrical full-resistive-MHD model of the resistive internal kink mode, the secondary current sheet becomes unstable to plasmoids for

$S \ge 1.6\times10^7,$

with a critical aspect ratio

$L/\delta \sim 70.$

For a representative case with

$S=2.3\times10^8,$

the sheet breaks into five small plasmoids and four tertiary current sheets (Ali et al., 2019). This is not yet a deeply recursive hierarchy, but it is unambiguous secondary fragmentation.

Solar-eruption simulations place the same physics in an eruptive current sheet under a rising flux rope. In the plasmoid-fed prominence-formation scenario, the plasmoid instability starts at $t\sim 40\,\mathrm{s}$ when the current-sheet Lundquist number reaches

$$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$0

and the paper explicitly frames this onset against the standard critical range near $$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$1–$$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$2 (Zhao et al., 2022).

A plausible implication is that “critical $$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$3” is not a universal scalar independent of regime. Semi-collisional sheets, cylindrical kink-driven sheets, and eruptive solar sheets all realize plasmoid onset, but with different controlling scales, aspect ratios, and global constraints (Bhat et al., 2018, Ali et al., 2019, Zhao et al., 2022).

3. Nonlinear evolution: chain formation, coalescence, and finite versus recursive fragmentation

Once onset occurs, nonlinear evolution is rarely a pure forward cascade. In the tokamak sawtooth problem, the breakup of the secondary sheet into a chain of small plasmoids initially speeds up reconnection, but subsequent coalescence produces a monster plasmoid that slows the reconnection process in the final stage and leads to partial reconnection (Ali et al., 2019). In that system, the nonlinear fate of the chain is as important as its onset.

Collisionless PIC studies reinforce the importance of coalescence. In a 2D plasmoid chain with guide field, multi-X-point tearing forms several plasmoids that then grow mainly by coalescing; the system does not develop a self-similar hierarchy of smaller and smaller plasmoids, and the dominant secondary process is anti-reconnection at merger sites (Markidis et al., 2012). Each chain plasmoid exhibits a strong out-of-plane core magnetic field and an out-of-plane electron current, while the disappearance of X-points during mergers is associated with anti-reconnection and a Hall quadrupole signature (Markidis et al., 2012).

Three-dimensional kinetic simulations show that whether a hierarchy emerges can depend on guide field. Without a guide field, a main reconnection site dominates and smaller reconnection regions are included in larger ones, producing a hierarchical structure of the plasmoid-dominated current sheet. With a uniform guide field

$$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$4

plasmoids have approximately the same size and the hierarchical structure does not emerge; instead, twisted flux ropes develop a strong core field and electron holes appear near plasmoids (Markidis et al., 2013). This directly contradicts any assumption that a plasmoid cascade is structurally universal.

Viscosity also modifies the nonlinear chain. In a Harris-sheet tearing simulation with full MHD, the nonlinear sequence is primary island growth, X-point collapse, a secondary elongated current sheet, and then plasmoid formation. Two distinct viscous regimes were identified: in the low-viscosity regime ($$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$5), plasmoid width increases sharply with viscosity, whereas in the viscosity-dominant regime ($$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$6), plasmoid size gradually decreases with viscosity (Ahmad et al., 2021). That result quantifies a non-monotonic dependence of secondary breakup on transport coefficients.

These studies collectively show that plasmoid cascade is best understood as a family of nonlinear pathways rather than a single asymptotic morphology. Some systems produce tertiary sheets and local recursion; others produce short chains that rapidly merge; others are dominated by a central monster plasmoid or by guide-field-stabilized, same-sized flux ropes (Ali et al., 2019, Markidis et al., 2013, Ahmad et al., 2021).

4. Turbulence-mediated plasmoid regimes

A major extension of the concept replaces the isolated current sheet by a turbulent ensemble of intermittent sheets. In analytic theory, strong MHD turbulence generates anisotropic sheet-like structures with

$$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$7

and plasmoids become dynamically important only when tearing can amplify turbulence-generated seed noise to nonlinear amplitude within the sheet lifetime, so that

$$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$8

The onset scale of the plasmoid-mediated range is then not a pure power law in $$S=\frac{L V_A}{\eta}, \qquad \delta_{SP}\sim L S^{-1/2}.$$9, and the predicted perpendicular energy spectrum in that range is steeper than $$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$0 because of Lambert-$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$1 corrections (Comisso et al., 2018).

Direct numerical simulations support a plasmoid-mediated turbulent subrange. In 2D resistive MHD turbulence, copious plasmoid formation appears at

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$2

while at

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$3

plasmoids are barely observed. The magnetic spectrum steepens from a standard inertial-range slope near $$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$4 to a plasmoid-mediated slope close to $$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$5, and the paper argues that dynamic alignment alone is insufficient; intermittency is essential because average aspect ratios are too small to explain onset (Dong et al., 2018).

High-resolution 3D MHD simulations strengthen this interpretation. In a periodic box with

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$6

elongated current sheets break into chains of small magnetic flux ropes, and the energy-transfer rate in the reconnection-driven range is controlled by the plasmoid growth rate. The spectrum steepens to

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$7

the break occurs near

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$8

and shell-to-shell transfer becomes self-similar only under a tearing-based normalization (Dong et al., 2022). In this usage, “plasmoid cascade” and “reconnection-driven energy cascade” are nearly synonymous.

Multiphase turbulence adds thermodynamic structure to this picture. In 3D simulations of the thermally unstable ISM, large current sheets unstable to plasmoid-mediated reconnection form within cold clumps, at cold–warm interfaces, and within the warm phase. The paper argues that many of these sheets fragment into plasmoids, although it does not provide recursive-cascade statistics (Fielding et al., 2022).

5. Collisionless generalizations and the limits of the term

A common misconception is that any multi-island collisionless reconnection automatically constitutes plasmoid-mediated reconnection in the resistive-MHD sense. The recent collisionless literature is more restrictive. In 2.5D antiparallel PIC simulations, secondary plasmoids can form inside the electron diffusion region, but the paper argues that they are generated by collisionless electron tearing rather than by the resistive plasmoid instability. It further finds that they do not enhance the reconnection rate and may disappear for realistic ion-electron mass ratio, with a threshold estimate

$$\left(\frac{L}{\rho_s}\right)^{14/9} \ll S \ll \left(\frac{L}{\rho_s}\right)^2,$$9

for secondary-plasmoid formation in the specific setup studied (Akutagawa et al., 16 Jul 2025). On that view, a true collisionless plasmoid cascade is not generically established in 2D antiparallel systems.

That caution does not eliminate kinetic cascade-like behavior; it changes its form. In 3D kinetic simulations without guide field, embedded reconnection regions produce a partially hierarchical chain, whereas guide field suppresses that hierarchy (Markidis et al., 2013). In a different kinetic route, lower-hybrid drift instability can grow at sheet edges, nonlinearly merge into larger structures, and produce magnetic islands through an inverse cascade from kinetic to fluid scales; thicker sheets or weaker density gradients instead favor Kelvin–Helmholtz vortices that later generate plasmoids secondarily (Thatikonda et al., 18 Jan 2026). This suggests that in collisionless plasmas the phrase “plasmoid cascade” may refer either to forward tearing of current sheets or to microinstability-driven upscale transfer that seeds plasmoids.

The most defensible generalization is therefore regime-qualified. In resistive and semi-collisional MHD, plasmoid cascade usually means hierarchical current-sheet fragmentation. In kinetic systems, one must distinguish chain formation, merger-dominated evolution, Hall-mediated flux-rope dynamics, and inverse-cascade seeding by microinstabilities (Akutagawa et al., 16 Jul 2025, Thatikonda et al., 18 Jan 2026).

6. Manifestations in solar, laboratory, fusion, and astrophysical systems

In solar eruptive current sheets, plasmoids are not only reconnection byproducts but also transport agents. In the PF$$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$0 scenario, multiple magnetic islands formed in the current sheet beneath an erupting flux rope transfer cool, dense chromospheric material into the rope; about half of the final prominence mass, $$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$1, is attributed to plasmoid-fed transport, with the other half supplied by later condensation (Zhao et al., 2022). Earlier flare-current-sheet work had already proposed that fragmentation between merging plasmoids could continue down to scales of order $$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$2 and could be linked to above-the-loop-top hard X-ray sources and narrowband dm-spikes (Karlický et al., 2011).

In fusion plasmas, plasmoid chains alter global relaxation dynamics rather than merely local reconnection geometry. In the Taylor problem, plasmoid chains replace the slower Rutherford phase under very small resistivity and viscosity (Comisso et al., 2014). In tokamak sawteeth, a secondary sheet driven by the nonlinear $$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$3 internal kink becomes plasmoid-unstable above

$$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$4

and the resulting monster plasmoid provides a mechanism for partial reconnection (Ali et al., 2019).

In laboratory streamer analogs, plasmoid production can be triggered by a pressure-curvature-driven loss of equilibrium rather than by a pre-existing supercritical sheet. In the Big Red Ball / WiPAL Parker-spiral current sheet, laminar quasi-periodic plasmoids appear at modest drive, while higher drive produces turbulent plasmoid ejection with plasmoids of many sizes; the heuristic ejection frequency is

$$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$5

and the transition to multiple simultaneous plasmoids is estimated near

$$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$6

in the experiment-specific Hall-MHD regime (Peterson et al., 2021).

In global astrophysical flows, current-sheet fragmentation can become recurrent rather than isolated. Axisymmetric GRMHD simulations of accretion onto a Kerr black hole with alternating-polarity magnetic loops show repeated current-sheet formation and plasmoid chains near the hole and funnel. Plasmoids are tracked by thresholds

$$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$7

and a representative tracked plasmoid grows while moving from roughly $$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$8 to $$\left(\frac{L}{\rho_s}\right)^{8/5} \ll S \ll \left(\frac{L}{\rho_s}\right)^2.$$9 (Nathanail et al., 2020). In the multiphase ISM, large current sheets unstable to plasmoid-mediated reconnection form regularly throughout the volume and are associated with high magnetic curvature over a broad temperature range (Fielding et al., 2022).

Taken together, these results suggest that plasmoid cascade is not a niche instability confined to idealized sheets. It is a recurrent organizing principle of reconnecting structures across driven, turbulent, radiatively cooling, semi-collisional, and kinetic plasmas, albeit with system-dependent morphology and with important caveats whenever a fully recursive hierarchy is inferred rather than directly demonstrated (Comisso et al., 2014, Dong et al., 2022, Fielding et al., 2022).