Self-consistent Simulations of Plasma-Neutral in a Partially Ionized Astrophysical Turbulent Plasma

A local turbulence model is developed to study energy cascades in the heliosheath and outer heliosphere (OH) based on self-consistent two-dimensional fluid simulations. The model describes a partially ionized magnetofluid OH that couples a neutral hydrogen fluid with a plasma primarily through charge-exchange interactions. Charge-exchange interactions are ubiquitous in warm heliospheric plasma, and the strength of the interaction depends largely on the relative speed between the plasma and the neutral fluid. Unlike small-length scale linear collisional dissipation in a single fluid, charge-exchange processes introduce channels that can be effective on a variety of length scales that depend on the neutral and plasma densities, temperature, relative velocities, charge-exchange cross section, and the characteristic length scales. We find, from scaling arguments and nonlinear coupled fluid simulations, that charge-exchange interactions modify spectral transfer associated with large-scale energy-containing eddies. Consequently, the turbulent cascade rate prolongs spectral transfer among inertial range turbulent modes. Turbulent spectra associated with the neutral and plasma fluids are therefore steeper than those predicted by Kolmogorov’s phenomenology. Our work is important in the context of the global heliospheric interaction, the energization and transport of cosmic rays, gamma-ray bursts, interstellar density spectra, etc. Furthermore, the plasma-neutral coupling is crucial in understanding the energy dissipation mechanism in molecular clouds and star formation processes.

💡 Research Summary

The paper presents a self‑consistent two‑dimensional fluid model that couples a partially ionized magnetofluid (plasma) with a neutral hydrogen fluid through charge‑exchange (CX) interactions, and investigates how these interactions modify turbulent energy cascades in the heliosheath and outer heliosphere. The authors formulate the plasma dynamics using the full set of magnetohydrodynamic (MHD) equations (mass, momentum, induction, and energy) and the neutral component with compressible hydrodynamic equations. CX is introduced via momentum (Q_M) and energy (Q_E) source terms that depend on the relative velocity, densities, temperatures, and the CX cross‑section. A characteristic CX wave number k_ce (or length scale ℓ_ce) emerges, representing the inverse of the CX mean free path; in the heliospheric plasma k_ce is typically much larger than the turbulent wave number k.

Numerically, the coupled equations are integrated with a pseudo‑spectral Fourier method for spatial discretization and a fourth‑order Runge‑Kutta scheme for time stepping. The code is parallelized with MPI, allowing high‑resolution simulations (e.g., 2048 × 2048 grid points). Initial conditions are isotropic, zero‑mean‑field fluctuations with a k⁻² spectrum and random phases, ensuring that any developed turbulence is self‑generated (free‑decay) rather than externally forced.

The key finding is that CX dramatically alters the nonlinear interaction time. In ordinary (uncoupled) turbulence the nonlinear time τ_nl scales as ℓ/v_ℓ ≈ (k v_k)⁻¹. With CX, the effective nonlinear time becomes τ_NL ≈ (k_ce/k) τ_nl, i.e., it is lengthened by the factor k_ce/k > 1 for scales larger than ℓ_ce. This prolonged interaction time slows the cascade, causing the inertial‑range energy spectra to steepen relative to the classic Kolmogorov (k⁻⁵ᐟ³) or Kraichnan (k⁻³ᐟ²) predictions. The simulations show that the plasma magnetic and kinetic energy spectra follow approximately k⁻²·³³ (spectral index ≈ 2.33), while the neutral fluid exhibits a steeper k⁻³·⁶⁷ (index ≈ 3.67). Analytic scaling arguments based on the modified τ_NL and the associated energy transfer rate ε lead to theoretical spectra E(k) ∝ ε^{2/3} k^{-7/3} for the plasma and E(k) ∝ ε^{2/3} k^{-11/3} for the neutrals, in good agreement with the numerical results.

The CX source terms are found to be most active at low wave numbers (large spatial scales), indicating that charge exchange couples large‑scale plasma motions directly to the neutral component, efficiently transferring momentum and energy. At higher k the CX coupling weakens, and conventional dissipative processes dominate. This scale‑dependent behavior explains why the cascade is altered primarily in the large‑scale part of the spectrum, while small‑scale turbulence remains essentially uncoupled.

The authors discuss several astrophysical implications. In the heliosphere, observed power spectra of magnetic fluctuations often deviate from Kolmogorov scaling; CX‑induced steepening offers a plausible mechanism. The prolonged nonlinear time may enhance Fermi‑type acceleration of cosmic rays and affect the transport of energetic particles. Moreover, similar partially ionized environments exist in molecular clouds and star‑forming regions, where CX could influence turbulence‑driven support against gravity and the rate of energy dissipation.

Limitations of the study are acknowledged: the simulations are two‑dimensional, free‑decaying, and omit external large‑scale forcing (e.g., supernova shocks). Extending the model to three dimensions, incorporating sustained driving, and exploring parameter regimes relevant to different astrophysical settings are suggested as future work.

In summary, the paper demonstrates that charge‑exchange interactions introduce a distinct length and time scale that lengthens nonlinear eddy turnover times, thereby steepening inertial‑range spectra in both plasma and neutral fluids. This mechanism provides a new perspective on energy transfer and dissipation in partially ionized astrophysical plasmas, with broad relevance to heliospheric physics, cosmic‑ray modulation, and the dynamics of interstellar and star‑forming media.