2D or not 2D? Exploring 3D relativistic magnetic reconnection dynamics with highly accurate numerical simulations

Fast reconnection in magnetically dominated plasmas is widely invoked in models of dissipation in pulsar winds, gamma-ray flares in the Crab nebula, and to explain the radio nanoshots of pulsars. When current sheets evolve reaching a critical inverse aspect ratio, scaling as $S^{-1/3}$ with the plasma Lundquist number, the so-called \textit{ideal} tearing instability sets in, with modes growing, independently of $S$, extremely rapidly on timescales of only a few light-crossing times of the sheet length. We present the first set of fully 3D simulations of current-sheet disruption triggered by the ideal tearing instability within the resistive relativistic MHD approximation, as appropriate in situations where the Alfvén velocity approaches the speed of light. We compare 3D setups with different initial conditions with their 2D counterparts, and we assess the impact of dimensionality and of the magnetic field topology on the onset, evolution, and efficiency of reconnection. In force-free configurations, 3D runs develop ideal tearing, secondary instabilities, and a thick, turbulent current layer, sustaining dissipation of magnetic energy longer than in 2D. In pressure-balanced current sheets with a null guide field, 2D reference runs show the familiar reconnection dynamics, whereas in 3D tearing dynamics is quenched after the linear phase, as pressure-driven modes growing on forming plasmoids outcompete plasmoid coalescence and suppress fast dissipation of magnetic energy. Taken together, these results suggest that the evolution and efficiency of reconnection depend sensitively on the local plasma conditions and current-sheet configuration, and can be properly captured only in fully 3D simulations.

💡 Research Summary

This paper presents the first fully three‑dimensional (3D) simulations of relativistic magnetic reconnection triggered by the ideal tearing instability within the resistive relativistic magnetohydrodynamics (ResRMHD) framework. The authors focus on current sheets whose inverse aspect ratio a/L scales as S⁻¹ᐟ³, the critical condition at which the tearing mode growth rate becomes independent of the Lundquist number S and proceeds on ideal Alfvénic timescales (∼L/c_A). They adopt a high‑order (fourth‑order) finite‑volume scheme with implicit treatment of the stiff resistive term, implemented on GPUs, allowing them to reach resolutions of 1024×768×768 for 3D runs and 1024×768 for 2D reference cases while keeping numerical diffusion low.

Two families of equilibria are investigated: (i) force‑free sheets (parameter ζ = 1) where magnetic pressure balances itself, and (ii) pressure‑balanced sheets (ζ = 0) with a null guide field where plasma pressure balances the magnetic pressure. Both configurations share the same upstream magnetization σ₀ = 10, plasma beta β₀ = 0.1, Lundquist number S = 10⁶, and perturbation amplitude ε = 10⁻³, ensuring that the ideal‑tearing condition a/L = S⁻¹ᐟ³ = 10⁻² is satisfied. The computational domain spans several current‑sheet half‑widths in the normal direction and two sheet lengths in the transverse directions, with reflective boundaries in x and periodic in y and z.

In the force‑free 3D case (ff3d), the linear phase exhibits the expected S‑independent exponential growth of the tearing mode. Once nonlinear, secondary instabilities (e.g., kink‑type, Kelvin‑Helmholtz‑type) develop, thickening the current layer into a turbulent, volume‑filling structure. Magnetic energy is converted into plasma kinetic and thermal energy over many Alfvén crossing times, and the dissipation persists longer than in the corresponding 2D run (ff2d), where the system quickly settles into a chain of plasmoids that merge and exhaust the available magnetic energy.

Conversely, in the pressure‑balanced 3D case (pb3d) the linear tearing growth is also present, but the subsequent evolution is dominated by pressure‑driven modes that grow on the forming plasmoids. These modes disrupt plasmoid coalescence, preventing the development of a sustained plasmoid cascade. As a result, after the initial linear phase the reconnection rate drops sharply, and magnetic energy dissipation is strongly quenched compared with the 2D pressure‑balanced reference (pb2d), which displays the classic fast plasmoid‑mediated reconnection.

The contrasting outcomes demonstrate that the dimensionality alone does not dictate reconnection efficiency; rather, the magnetic topology (presence or absence of a guide field) and the pressure balance across the sheet critically control whether secondary 3D instabilities enhance or suppress reconnection. The study also highlights the necessity of sufficiently high physical resistivity (or equivalently, high Lundquist number) and high‑order numerical methods to resolve the thin diffusion layer where ideal tearing initiates.

Overall, the work provides compelling evidence that realistic astrophysical reconnection—such as that occurring in pulsar winds, Crab nebula flares, or magnetar magnetospheres—cannot be faithfully captured by 2D models. Depending on the local plasma conditions, reconnection may either proceed as a prolonged, turbulent energy‑release process (force‑free case) or be rapidly self‑limited by pressure‑driven instabilities (pressure‑balanced case). These findings have direct implications for interpreting the variability and energetics of high‑energy transients in relativistic astrophysical environments.