---
title: Tearing-Mode Instability in Plasmas
url: https://www.emergentmind.com/topics/tearing-mode-instability-6727778c-8cd5-4cbf-b489-44565721a347
type: topic
---

# Tearing-Mode Instability in Plasmas

Searching arXiv for recent and foundational papers on tearing-mode instability, emphasizing the papers provided and closely related work.
Tearing-mode instability is a reconnecting instability of thin current sheets in which a magnetic equilibrium with field reversal becomes unstable to perturbations localized about a resonant surface, so that non-ideal effects break flux freezing and allow the formation of magnetic islands or plasmoids. In the classical picture the instability is organized by the tearing stability parameter $\Delta'$, the Lundquist number $S$, and the current-sheet aspect ratio $a/L$; in later developments it has been reformulated for “ideal” tearing, Hall and electron-MHD limits, partially ionized and gyrotropic plasmas, gyrokinetic tokamak configurations, and fully three-dimensional current-sheet geometries [1603.04995] [1504.07036].

## 1. Classical resistive-MHD formulation

In resistive MHD, tearing is usually formulated for a one-dimensional current sheet of half-thickness $a$ embedded in a system of length $L$, with reversing field $B(y)=B_0 F(y/a)\,\hat x$ or, in Harris-sheet form, $B_x(y)=B_0\tanh(y/a)$. The standard control parameters are the Alfvén speed $V_A=B_0/\sqrt{\mu_0\rho}$, the Alfvén time $\tau_A=L/V_A$, the Lundquist number $S\equiv L V_A/\eta$, and the inverse aspect ratio $\epsilon\equiv a/L\ll1$. In the classical resistive theory, the fastest-growing mode satisfies $\gamma_{\max}\tau_A\sim S^{-1/2}\epsilon^{-3/2}$, while the constant-$\psi$ approximation for a Harris sheet uses $\Delta'(k)=2\bigl(1/(ka)-ka\bigr)$ as the tearing-stability index [2006.03957] [2007.00607].

The inner–outer asymptotic structure is central. Outside the resistive layer the perturbation obeys ideal-MHD equations, and inside the narrow reconnecting layer resistive diffusion regularizes the singularity at the resonant surface. Instability requires $\Delta'>0$; in the Harris-sheet normalization this suggests that sufficiently long-wavelength modes satisfy the basic tearing criterion. In local sheet variables, the classical large-$\Delta'$ estimate is $\gamma\tau_A^\ast\sim {S^\ast}^{-1/2}$, where $S^\ast=a c_A/\eta$ and $\tau_A^\ast=a/c_A$ [1603.04995].

This classical framework already distinguishes tearing from other current-sheet instabilities. It is driven by magnetic shear and reconnection physics rather than by pure hydrodynamic shear, and it produces O-points and X-points through reconnection of oppositely directed magnetic field lines. That structure remains recognizable across later Hall, kinetic, and toroidal generalizations, even when the inner-layer physics is substantially altered.

## 2. Ideal tearing and the disruption of thin current sheets

A major reformulation arises when the current-sheet aspect ratio is allowed to scale with $S$ as $a/L\sim S^{-\alpha}$. For the fastest tearing mode, the growth rate then scales as $\gamma\tau_A\sim S^{-1/2+3\alpha/2}$. Requiring the exponent to vanish in the $S\to\infty$ limit gives the critical exponent $\alpha=1/3$, hence the “ideal-tearing” threshold
$$
a/L\sim S^{-1/3}.
$$
At that threshold the growth rate becomes independent of $S$ to leading order, with $\gamma_m\tau_A\to \text{const}\simeq0.62$ and $\gamma_m\simeq0.6/\tau_A$; the fastest mode also satisfies $k_m a\simeq0.6\,S^{-1/6}$ and $k_mL\simeq0.6\,S^{1/6}$ [1504.07036].

Two-dimensional compressible resistive-MHD simulations for $S=10^5,10^6,10^7$ confirm these scalings. They recover the predicted dispersion curves and eigenmode structure, including the even magnetic perturbation $b_x(x)$ peaked inside the sheet, the odd transverse flow $v_x(x)$ vanishing at the origin, and the even parallel flow $v_y(x)$ in quadrature with $b_x$ along the ignorable direction. For the least diffusive case reported, $S=10^7$, the tearing instability reaches $\gamma\approx0.6\,\tau_A^{-1}$, i.e. macroscopic Alfvénic growth [1504.07036].

The nonlinear regime is hierarchical. In multi-mode runs, initially formed islands elongate the X-point regions into new local current sheets of half-length $L^\ast$ and half-thickness $a^\ast$. These secondary sheets satisfy the same critical relation on local scales,
$$
a^\ast/L^\ast\sim {S^\ast}^{-1/3},\qquad S^\ast=(L^\ast/L)S,
$$
with measured aspect ratios $L^\ast/a^\ast\simeq 50$ at $S=10^5$, $80$ at $S=10^6$, and $200$ at $S=10^7$, fitted by $L^\ast/a^\ast\simeq k\,{S^\ast}^{1/3}$ with $k\simeq2.1$–$2.3$ [1504.07036]. Because $\tau_A^\ast=L^\ast/c_A\ll\tau_A$, each successive generation reconnects faster than the previous one. This produces a cascading, accelerating plasmoid regime.

The comparison with Sweet–Parker scaling is decisive. A Sweet–Parker sheet has $a/L\sim S^{-1/2}$, which is much thinner than the critical $S^{-1/3}$ threshold. The simulations and scaling analysis therefore support the conclusion that a sheet disrupts through tearing before reaching the fully collapsed Sweet–Parker aspect ratio; in that sense, the Sweet–Parker configuration is likely to be never realized in nature for sufficiently large $S$ [1603.04995].

## 3. Hall, partial-ionization, and gyrotropic extensions

In electron magnetohydrodynamics, the tearing-type Hall instability is driven by the combined action of Hall-induced whistler oscillations and shear in the electron current velocity. A local necessary condition is that the background field and the vorticity of the electron current velocity be antiparallel:
$$
B(x)\,\frac{d^2B(x)}{dx^2}<0,
$$
or equivalently,
$$
\mathbf B_0\cdot\nabla\times(\nabla\times\mathbf B_0)>0.
$$
For large Hall number $R_H$, asymptotic matching yields $\sigma\propto R_H^{2/3}$ for vacuum-wall tearing and $\sigma\propto R_H^{3/4}$ for periodic or broad profiles; in physical units, $\sigma_{\rm phys}\propto B_0^\alpha\eta^{1-\alpha}$ with $\alpha=2/3$ or $3/4$ respectively. Perfect-conductor walls strongly suppress or completely quench the mode by forbidding the essential helicoidal motion at the boundaries [2106.13438].

When a guide field and Hall physics are included in resistive Hall-MHD, the fastest linear instability remains the parallel mode with $k_z=0$; a strong guide field does not move the global maximum to oblique wavevectors. Instead, it shifts the resonant surface for oblique modes, suppresses their wave-like inner-layer structure, and makes the eigenfunctions asymmetric. In the Hall case, oblique modes acquire a complex frequency and propagate along the guide-field direction, whereas for $B_g=0$ the propagation is non-dispersive with $\omega(k_z)=d_i B_x' k_z$ [2007.00607].

Partial ionization modifies both resistivity and inertia. Electron–neutral collisions enhance the magnetic diffusivity, while ion–neutral collisions change the effective Alfvénic response. In the coupled regime the ideal-tearing threshold becomes
$$
\epsilon_c=S^{-1/3}(\rho/\rho_i)^{1/6},
$$
in the intermediate regime
$$
\epsilon_c=S^{-1/3}(\nu_{in}\tau_A)^{-1/6},
$$
and in the decoupled regime the fully ionized result $\epsilon_c=S^{-1/3}$ is recovered up to negligible corrections. The intermediate regime is stabilizing in the sense that ion–neutral coupling acts like an effective viscosity, although the electron–neutral contribution hidden in $\eta_m$ can reduce $S$ and thereby shift the net threshold in the opposite direction [2006.03957].

Weakly collisional gyrotropic plasmas introduce an additional layer of modification. For a force-free Harris sheet with prescribed equilibrium pressure anisotropy, the classical $S^{-1/2}$ scaling of the maximum growth rate is retained but acquires the prefactor $\gamma_{\max}\tau_A\sim\mathcal A^{1/2}\mathcal R_0^{-1/4}S^{-1/2}$, where $\mathcal A=1-\Delta\beta_0/2$ and $\mathcal R_0$ depends on $\beta_0$, $\Delta\beta_0$, and the gyrotropic closure. Positive $\Delta\beta_0$ suppresses tearing and broadens the inner layer; negative $\Delta\beta_0$ enhances growth and shifts the fastest mode to larger wavenumber [2606.22535]. A separate non-ideal CGL-MHD treatment with self-consistent pressure-anisotropy fluctuations finds a distinct high-$\beta$ scaling even for initially isotropic equilibria,
$$
\gamma_{\max}\tau_A\propto \beta^{-1/4},
$$
thereby breaking the standard MHD $\beta$-independence in the high-$\beta$ regime [2503.12702].

## 4. Three-dimensionality, transverse fields, and non-modal onset

The tearing instability is not confined to strictly two-dimensional sheets. For a flux-tube-like equilibrium of the form $\mathbf B_0(x,y)=B_0 f(x)g(y)\,\hat{\mathbf z}$, linear analysis and direct numerical simulations show that a tearing-like mode persists in three dimensions even without a guide field. The principal effect of the modulation $g(y)$ is a reduction in the net growth rate by the factor
$$
\gamma_{3D}/\gamma_{2D}=\frac{\int g(y)^{1/2}\,dy}{\int dy},
$$
while the $k$- and $S$-scalings remain similar to the two-dimensional case. The local resistive layer becomes $y$-dependent, with $\delta(y)\sim\bigl(\eta\gamma/(kB_0 g(y))\bigr)^{1/4}$, and the simulations show synchronized growth across different $y$-slices rather than a mere stack of uncoupled two-dimensional sheets [2412.10065].

A transverse magnetic-field component can instead quench tearing. In viscoresistive incompressible MHD, introducing $\xi=B_z/B_0$ into an initial Harris sheet destroys the static equilibrium and generates a widening neutral layer with half-width
$$
w(t)=\xi V_A t.
$$
The effective sheet width $a_{\rm eff}=a+w$ reduces $\Delta'$ and depresses the entire $\gamma(k)$ curve. An approximate stabilization threshold is $\xi\gtrsim S^{-3/4}$, and direct simulations confirm that sufficiently strong transverse field produces a broad, stable current layer with strong shear flow instead of a plasmoid chain [2407.09996].

The onset problem is also non-modal. For the classical incompressible visco-resistive linear operator, pseudospectral analysis shows strong non-normality and transient amplification even before the asymptotic tearing eigenmode dominates. The maximum possible $L_2$-norm growth scales as $G_{\max}\sim S^{1/4}$ on times $t_{\max}\sim S^{1/4}$ for $S\gg1$, which is faster than the modal tearing timescale $O(S^{1/2})$. The optimal initial perturbations are localized wave packets concentrated at the current sheet, i.e. pseudomodes rather than eigenmodes [1808.10208].

These results complicate a common simplification in which tearing onset is identified solely with the dominant unstable eigenvalue. A plausible implication is that, at large $S$, the early evolution of thin sheets can be controlled by transiently amplified localized noise before the classical exponential tearing branch becomes asymptotically dominant.

## 5. Kinetic, gyrokinetic, and tokamak manifestations

In toroidal gyrokinetic simulations of classical current-gradient-driven tearing, the instability exhibits a collisionless plateau at low collisionality, a smooth transition to a semi-collisional regime, and a growth rate that follows the Fitzpatrick large-$\Delta'$ scaling $\gamma\propto\eta^{1/7}$ rather than the steeper Drake–Lee or classical resistive exponents. Even in the absence of pressure gradients the mode has a residual finite rotation frequency, attributed to toroidal finite-Larmor-radius effects; with pressure gradients it rotates at the electron diamagnetic frequency at low collisionality, while at high collisionality the rotation can reverse to the ion diamagnetic direction [1408.0112].

Microtearing extends tearing parity to ion-Larmor-radius scales. In the Mega Ampere Spherical Tokamak, linear gyrokinetic-Maxwell calculations show a tearing-parity instability with radial width $d\sim\rho_i$, peak growth at $k_y\rho_i\approx0.5$–$0.7$, and instability for $0.05\lesssim\bar\nu_e/\omega\lesssim2$, with maximum growth near $\bar\nu_e\sim\omega$. The free energy comes from the electron temperature gradient, but the dominant destabilization is not the classical thermal-force or banana-boundary-layer mechanism; rather, the principal drive is from magnetic drifts, with electrostatic potential coupling also essential [1110.3277].

A distinct collisionless fine-scale tearing-parity mode exists even at zero collision frequency in slab gyrokinetics. It is also driven by electron temperature gradient, requires finite electron FLR effects, has $k_x\rho_e\sim1$, persists in the electrostatic limit, and is actually stabilized by electromagnetic feedback. The analytic reduction identifies it with a tearing-parity slab ETG mode that can be more unstable than the more commonly studied twisting-parity ETG branch [2004.10073].

Self-consistent collisionless tokamak tearing in global ORB5 simulations adds nonlinear structure. For flat profiles, the linear growth rate obeys $\gamma_{TM}\propto 1/\beta_e$ for $\beta<0.4\%$ and $\propto m_e^{1/2}$ in the mass-ratio scan; the initial nonlinear island width satisfies $\delta_{\rm sat}\propto \gamma_{TM}^{1/2}$. Large islands then drive zonal current redistribution and island-induced $\mathbf E\times\mathbf B$ shear flows at the separatrix, which destabilize Kelvin–Helmholtz turbulence and cause strong island decay. At sufficiently high $\beta_e$, the classical tearing mode is suppressed and a twisting-parity mode with the same helicity is destabilized instead [2410.11498].

In neoclassical tokamak physics, tearing also appears as a seeded nonlinear branch. NIMROD simulations of impurity radiation cooling show a two-step mechanism in which a local helical pressure perturbation drives a diamagnetic-current seed island, after which a heuristic neoclassical electron-viscosity closure generates a perturbed bootstrap current that drives neoclassical tearing-mode growth. In those simulations the nonlinear growth rate is proportional to the electron neoclassical viscosity, while the pure-resistive seed behavior depends on impurity radiation power [2410.20787].

## 6. Observational evidence and astrophysical significance

Direct in-situ evidence for tearing-mode onset prior to reconnection has been reported in the terrestrial magnetotail. Using Cluster electron pitch-angle distributions and neural-network outlier detection, fifteen candidate tearing events were identified, all satisfying the collisionless threshold
$$
\frac{B_n}{B_0}\le \sqrt{\frac{m_e T_e}{m_i T_i}}
$$
for the longest-wavelength mode. The characteristic sequence is a field-aligned counter-streaming beam phase followed by isotropization and heating, with reconnection signatures occurring minutes later; in the statistical sample, reconnection follows ten of the fifteen events with delays between $5$ min $32$ s and $8$ min $20$ s, with mean $\simeq7$ min [2202.12389].

The ideal-tearing framework has immediate astrophysical implications. In the solar corona, the collisional Lundquist number is quoted as $S\sim10^{12}$, so the critical threshold is $a/L\sim S^{-1/3}\sim10^{-4}$ rather than the much thinner Sweet–Parker value. Once a coronal current sheet reaches that critical aspect ratio, tearing and plasmoid formation proceed on Alfvénic times $\tau_A\sim O(1$–$10\,{\rm s})$, providing a route to impulsive flaring. The same hierarchical threshold has been proposed for the magnetotail, pulsar wind nebulae, relativistic pair plasmas, and turbulent current-sheet ensembles in which sheets of many lengths thin until they reach the local $S^{-1/3}$ condition [1504.07036].

In relativistic collisionless pair plasmas, the growth rate depends sensitively on how the current sheet is supported. For pressure-supported sheets with a uniform guide field, hardening the particle spectrum suppresses tearing; for force-free sheets the maximum growth rate becomes independent of spectral hardness. Applied to relativistic MHD turbulence, the same theory yields a tearing-mediated disruption scale for turbulent eddies and a critical magnetization above which charge starvation prevents the tearing instability [2410.16167].

Taken together, these developments place tearing-mode instability at the center of contemporary reconnection theory, but not as a monolithic mechanism. In some regimes it is a classical resistive eigenmode; in others it is Hall-mediated, semi-collisional, microtearing, or collisionless and gyrotropic. In two dimensions it readily produces plasmoid chains, whereas in three-dimensional stochastic reconnection it may be important only at the initial stage, with Kelvin–Helmholtz instability later dominating turbulence generation [1909.09179]. The unifying feature is the conversion of current-sheet free energy into reconnecting magnetic topology once the relevant inner-layer physics allows the ideal constraint to fail.

Source: https://www.emergentmind.com/topics/tearing-mode-instability-6727778c-8cd5-4cbf-b489-44565721a347