Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

Mode Energy Partition in Partially Ionized Compressible MHD Turbulence

Yue Hu

TL;DR

This study addresses how neutral–ion damping in partially ionized interstellar medium (ISM) environments reshapes compressible MHD turbulence. Using 3D two-fluid simulations with varying ion–neutral coupling, the authors perform a Fourier-based mode decomposition into Alfvén, slow, and fast components to quantify spectral changes and energy partitioning. They find that in strong coupling, Alfvén and slow modes exhibit near-Kolmogorov spectra ($E(k) o k^{-5/3}$) while fast modes trend toward $k^{-2}$; as coupling weakens and damping dominates, all spectra steepen—Alfvén approaching $k^{-4}$—and slow-mode energy increases dramatically toward small scales, reaching up to ~70% near the damping scale, with fast-mode energy remaining roughly constant at ~10%. These results imply that partial ionization shifts turbulent energy toward the slow mode in the damping regime, with important consequences for ISM processes such as cosmic-ray transport and acceleration.

Abstract

We investigate how neutral-ion collisional damping modifies the spectral properties and energy partition of compressible MHD turbulence using a suite of 3D two-fluid simulations. By systematically varying the neutral-ion coupling strength and decomposing the turbulent velocity field into Alfvén, slow, and fast modes, we quantify how each mode responds to the transition from strong to weak coupling. In the strong-coupling regime, the Alfvén and slow modes follow nearly Kolmogorov $k^{-5/3}$ spectra and dominate the kinetic energy budget, while fast modes exhibit a steeper spectrum and contribute $\sim$10\% of the total energy. As the coupling weakens and neutral-ion damping becomes significant, all mode spectra steepen, with the Alfvén mode approaching a dissipation-dominated $k^{-4}$ spectrum. The slow-mode energy fraction increases substantially toward small scales, reaching up to $\sim$70\% near the damping scale, whereas the Alfvén-mode fraction decreases correspondingly. In contrast, the fast-mode energy fraction remains largely insensitive to coupling strength. These results demonstrate that partial ionization not only steepens the turbulent spectra but also reshapes the mode energy distribution, enhancing the relative importance of the slow mode while suppressing Alfvén mode in the damping regime. Our findings have important implications for turbulence-driven processes in the partially ionized interstellar medium, including cosmic-ray transport and acceleration.

Mode Energy Partition in Partially Ionized Compressible MHD Turbulence

TL;DR

This study addresses how neutral–ion damping in partially ionized interstellar medium (ISM) environments reshapes compressible MHD turbulence. Using 3D two-fluid simulations with varying ion–neutral coupling, the authors perform a Fourier-based mode decomposition into Alfvén, slow, and fast components to quantify spectral changes and energy partitioning. They find that in strong coupling, Alfvén and slow modes exhibit near-Kolmogorov spectra () while fast modes trend toward ; as coupling weakens and damping dominates, all spectra steepen—Alfvén approaching —and slow-mode energy increases dramatically toward small scales, reaching up to ~70% near the damping scale, with fast-mode energy remaining roughly constant at ~10%. These results imply that partial ionization shifts turbulent energy toward the slow mode in the damping regime, with important consequences for ISM processes such as cosmic-ray transport and acceleration.

Abstract

We investigate how neutral-ion collisional damping modifies the spectral properties and energy partition of compressible MHD turbulence using a suite of 3D two-fluid simulations. By systematically varying the neutral-ion coupling strength and decomposing the turbulent velocity field into Alfvén, slow, and fast modes, we quantify how each mode responds to the transition from strong to weak coupling. In the strong-coupling regime, the Alfvén and slow modes follow nearly Kolmogorov spectra and dominate the kinetic energy budget, while fast modes exhibit a steeper spectrum and contribute 10\% of the total energy. As the coupling weakens and neutral-ion damping becomes significant, all mode spectra steepen, with the Alfvén mode approaching a dissipation-dominated spectrum. The slow-mode energy fraction increases substantially toward small scales, reaching up to 70\% near the damping scale, whereas the Alfvén-mode fraction decreases correspondingly. In contrast, the fast-mode energy fraction remains largely insensitive to coupling strength. These results demonstrate that partial ionization not only steepens the turbulent spectra but also reshapes the mode energy distribution, enhancing the relative importance of the slow mode while suppressing Alfvén mode in the damping regime. Our findings have important implications for turbulence-driven processes in the partially ionized interstellar medium, including cosmic-ray transport and acceleration.

Paper Structure

This paper contains 13 sections, 1 equation, 5 figures.

Table of Contents

  1. Introduction
  2. Methodology
  3. Two-fluid simulations
  4. Initial conditions
  5. Ion--neutral coupling
  6. Mode decomposition
  7. Results
  8. Neutral-ion collisional damping of MHD turbulence
  9. Mode decomposition
  10. Kinetic energy spectra of Alfvén, slow, and fast modes
  11. Kinetic energy partition among Alfvén, slow, and fast modes
  12. Conclusion
  13. Acknowledgments

Figures (5)

  • Figure 1: Left and middle panels: 2D slices of the ion (left) and neutral (middle) velocity fields, normalized by their mean amplitudes. The viewing direction is perpendicular to the mean magnetic field, which is oriented vertically. Right panels: Turbulent kinetic energy spectra of ions and neutrals. The shaded regions indicate temporal variations over several snapshots after the turbulence has reached a statistically steady state. The dotted lines show a $-5/3$ slope for reference to the Kolmogorov inertial-range scaling. The dash–dotted lines mark the theoretically predicted neutral–ion parallel decoupling and damping wavenumbers for Alfvénic MHD turbulence. These wavenumbers are not shown when they exceed the numerical dissipation scale, $k \sim 50$. Modified from 2024MNRAS.527.3945H.
  • Figure 2: 2D slices of the kinetic energy fraction in the decomposed MHD modes: Alfvén (left), slow (middle), and fast (right). Panel (a) shows the strong-coupling case with $\gamma_{\rm d}=10^5$, whereas panel (b) corresponds to the weak-coupling case with $\gamma_{\rm d}=25$.
  • Figure 3: Turbulent kinetic energy spectra $E_i(k)$ of ions' decomposed MHD modes: Alfvén (top), slow (middle), and fast (bottom). The subscript $i = a, s, {\rm or}~f$ denotes the corresponding mode. Panel (a) shows the strong-coupling case with $\gamma_{\rm d}=10^5$, while panel (b) corresponds to the weak-coupling case with $\gamma_{\rm d}=25$. The dashed and dash–dotted lines indicate the Kolmogorov ($-5/3$) and Burgers ($-2$) scalings expected in the inertial range, while the dotted line with slope $-4$ illustrates the steepened spectrum in the dissipation-dominated regime caused by neutral-ion damping.
  • Figure 4: Fraction of turbulent kinetic energy in the Alfvén (red), slow (blue), and fast (green) modes as a function of wavenumber. Each fraction is defined as $E_i(k)/E(k)$, where $E(k)=E_a(k)+E_s(k)+E_f(k)$ and the subscript $i = a, s, {\rm or}~f$ denotes the corresponding mode.
  • Figure 5: Turbulent kinetic energy fraction in the Alfvén (red), slow (blue), and fast (green) modes as a function of the coupling strength, parameterized by the drag coefficient $\gamma_{\rm d}$. $\gamma_{\rm d}=1\times10^5$ represents the strong-coupling case and the decoupling and damping of MHD turbulence gradually happen towards smaller $\gamma_{\rm d}$. Each fraction is defined as the total kinetic energy of each mode normalized by the total kinetic energy of all three modes.