Skip to content

Plasmoid-Mediated 2D Magnetic Reconnection in Partially Ionized Plasmas

Yue Hu, Siyao Xu, Grzegorz Kowal, James M. Stone, Alex Lazarian, Hui Li

astro-ph.HEarXiv:2608.08448

Abstract

Magnetic reconnection in partially ionized plasmas is an important channel for energy release. While the plasmoid instability is well characterized in 2D fully ionized plasmas, its behavior in the presence of neutral-dominated plasma remains poorly understood in the nonlinear, high-Lundquist-number (S = 105) regime. We present high-resolution (16384 × 4096 cells) two-dimensional two-fluid (ion + neutral) simulations of Harris-sheet reconnection with upstream plasma beta β= 2, comparing fully ionized and partially ionized (ionization fraction ξ= 10-1 and 10-2) regimes. Neutral-ion decoupling accelerates the linear tearing stage and alters the plasmoid hierarchy: the large-scale ``monster'' plasmoid that dominates the fully ionized case is suppressed, and the sheet instead fragments into a dense chain of sub-scale plasmoids. Below the neutral-ion decoupling scale dec, ions concentrate into the plasmoids, reaching peak overdensities ρi/ρi,0 ≈ 10 (ξ= 10-1) and 3-5×103 (ξ= 10-2), while the neutrals remain comparatively smooth. This local pile-up raises the ionization fraction and recouples the two fluids within the plasmoids. Measured from the out-of-plane electric field at the reconnection sites, the reconnection rate in the ξ= 10-2 case achieves R rec≈0.01, whereas the ξ= 10-1 case rises to a rate ≈0.02 and further 0.035 when apparent coalescence occurs. In the ξ= 10-2 case, the ambipolar drift drives a rapid ion inflow 0.5\,vA,0 into the layer at the same reconnection sites, far above the neutral inflow velocity 0.1\,vA,0. Here, vA,0 is the upstream total Alfvén speed.

Create a lesson