Cosmic Ray Superdiffusion and Mirror Diffusion in Partially Ionized and Turbulent Medium

Abstract: Understanding cosmic ray (CR) diffusion in a partially ionized medium is both crucial and challenging. In this study, we investigate CR superdiffusion in turbulent, partially ionized media using high-resolution 3D two-fluid simulations that treat ions and neutrals separately. We examine the influence of neutral-ion decoupling and the associated damping of turbulence on CR propagation in both transonic and supersonic conditions. Our simulations demonstrate that neutral-ion decoupling significantly damps velocity, density, and magnetic field fluctuations at small scales, producing spectral slopes steeper than those of Kolmogorov and Burgers scaling. We also identify an intermediate coupling regime in which neutrals remain partially coupled with ions, leading to a steepening of kinetic energy spectra in both fluids. Shocks can enhance the neutral-ion coupling, thereby reducing the differences between neutral and ion density structures. Moreover, the damping of magnetic field fluctuations substantially decreases pitch angle scattering, which increases CR parallel mean free paths. As a result, CR perpendicular transport transitions between two distinct superdiffusive regimes: a scattering-dominated regime with perpendicular displacement proportional to $t^{3/4}$, and a scattering-free regime dominated by magnetic field line wandering, with perpendicular displacement scaling as $t^{3/2}$. When the pitch angle is large, the effects of magnetic mirroring, naturally arising in magnetohydrodynamic turbulence, become significant, enhancing the confinement (but not fully trapping) of CRs and inducing oscillations in their perpendicular displacement. These results highlight the necessity of incorporating two-fluid effects for accurately modeling CR transport in partially ionized environments such as molecular clouds and dense interstellar clumps.