Two-Fluid Magnetohydrodynamic Simulations of Relativistic Magnetic Reconnection

We investigate the large scale evolution of a relativistic magnetic reconnection in an electron-positron pair plasma by a relativistic two-fluid magnetohydrodynamic (MHD) code. We introduce an inter-species friction force as an effective resistivity to dissipate magnetic fields. We demonstrate that magnetic reconnection successfully occurs in our two-fluid system, and that it involves Petschek-type bifurcated current layers in later stage. We further observe a quasi-steady evolution thanks to an open boundary condition, and find that the Petschek-type structure is stable over the long time period. Simulation results and theoretical analyses exhibit that the Petschek outflow channel becomes narrower when the reconnection inflow contains more magnetic energy, as previously claimed. Meanwhile, we find that the reconnection rate goes up to ~1 in extreme cases, which is faster than previously thought. The role of the resistivity, implications for reconnection models in the magnetically dominated limit, and relevance to kinetic reconnection works are discussed.

💡 Research Summary

The paper presents a comprehensive study of relativistic magnetic reconnection in an electron‑positron (pair) plasma using a relativistic two‑fluid magnetohydrodynamic (MHD) code. Unlike conventional single‑fluid relativistic MHD, which treats the plasma as a single conducting fluid and therefore cannot capture inter‑species effects, the authors model electrons and positrons as separate fluids and introduce an inter‑species friction term that acts as an effective resistivity. This friction term dissipates magnetic energy into heat within the current sheet, providing a physically motivated analogue of kinetic collision or charge‑separation effects that are otherwise absent in fluid models.

The numerical setup consists of a two‑dimensional Harris current sheet placed in the x‑z plane with anti‑parallel magnetic fields. Open boundary conditions are employed so that plasma and electromagnetic fields can freely leave the computational domain, preventing artificial reflections that could otherwise contaminate the reconnection dynamics. The simulation begins with a Sweet‑Parker‑type thin current layer and a modest reconnection rate (R≈0.05–0.1). As the system evolves, a Petschek‑type bifurcated current sheet emerges, characterized by two narrow outflow jets separated by slow‑mode shocks. This transition demonstrates that the two‑fluid framework can naturally produce the fast‑reconnection geometry originally proposed by Petschek, even in the relativistic regime.

A central parameter explored in the study is the magnetization σ = B²/(4πw), where w is the plasma enthalpy. By varying σ from ≈1 up to ≈20, the authors confirm the theoretical prediction that the Petschek outflow channel narrows as σ⁻¹/². In highly magnetized cases (σ≫1) the outflow speed approaches the speed of light, and the opening angle of the slow‑mode shocks becomes extremely small, leading to a very thin exhaust channel. This behavior matches analytic models (e.g., Lyubarsky 2005) and provides a clear illustration of how magnetic energy dominance reshapes the reconnection geometry.

Perhaps the most striking result concerns the reconnection rate itself. By adjusting the friction coefficient (i.e., the effective resistivity), the authors observe rates that exceed the conventional relativistic Petschek limit of ~0.3. In extreme configurations the normalized reconnection rate reaches values of order unity (R≈0.8–1.0), indicating that magnetic flux can be reconnected almost as fast as the Alfvén speed in the inflow region. The authors attribute this acceleration to the combination of a very thin current sheet (produced by strong friction) and the relativistic inertia of both species, which together allow the electric field to approach the ideal limit E≈B.

The paper also discusses the physical interpretation of the friction term. While it is introduced phenomenologically, it can be viewed as a coarse‑grained representation of kinetic processes such as particle scattering, anomalous resistivity, or micro‑instabilities that generate effective momentum exchange between electrons and positrons. Consequently, the two‑fluid model occupies a middle ground: it lacks the full kinetic detail of particle‑in‑cell (PIC) simulations (e.g., electron‑scale diffusion regions, non‑thermal particle acceleration) but captures large‑scale dynamics, energy conversion, and the formation of Petschek‑type structures with far lower computational cost.

In the discussion, the authors compare their findings with recent kinetic reconnection studies. They note that PIC simulations of relativistic pair plasmas also report fast reconnection and narrow exhausts, but kinetic effects such as guide‑field dependence and particle energization are not directly addressed in the fluid model. They suggest future work that could couple the two‑fluid MHD equations with a dynamic prescription for the friction coefficient derived from kinetic theory, thereby bridging the gap between fluid and kinetic descriptions.

Overall, the study demonstrates that relativistic two‑fluid MHD can successfully reproduce Petschek‑type reconnection, that the outflow geometry contracts with increasing magnetization, and that reconnection can proceed at rates approaching unity when the effective resistivity is sufficiently strong. These results have important implications for high‑energy astrophysical environments—such as pulsar wind nebulae, magnetar flares, and relativistic jets—where magnetic energy dominates the plasma dynamics and rapid magnetic energy release is required to explain observed high‑energy transients.