Plasmoid-Mediated 2D Magnetic Reconnection in Partially Ionized Plasmas

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=10^5$) regime. We present high-resolution ($16384\times 4096$ cells) two-dimensional two-fluid (ion $+$ neutral) simulations of Harris-sheet reconnection with upstream plasma beta $\beta =2$, comparing fully ionized and partially ionized (ionization fraction $\xi =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 ${\ell }_{\mathrm{dec}}$, ions concentrate into the plasmoids, reaching peak overdensities ${\rho }_i/{\rho }_{i,0}\approx 10$ ($\xi =10^{-1}$) and $3-5\times 10^3$ ($\xi =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 $\xi =10^{-2}$ case achieves $R_{\mathrm{rec}}\approx 0.01$, whereas the $\xi =10^{-1}$ case rises to a rate $\approx 0.02$ and further $0.035$ when apparent coalescence occurs. In the $\xi =10^{-2}$ case, the ambipolar drift drives a rapid ion inflow $\sim 0.5v_{A,0}$ into the layer at the same reconnection sites, far above the neutral inflow velocity $\sim 0.1v_{A,0}$. Here, $v_{A,0}$ is the upstream total Alfvén speed.

24 pages, 16 figures, submitted to ApJ