Reduced Effective Viscosity from Anisotropic Transport and Plasma Instabilities in the Sloshing Cores of Galaxy Clusters

Technical Summary: Reduced Effective Viscosity from Anisotropic Transport and Plasma Instabilities in the Sloshing Cores of Galaxy Clusters

Problem Statement
The intracluster medium (ICM) of galaxy clusters is a hot, weakly magnetized plasma where the Coulomb mean free path ($\lambda_{mfp}$) is significantly smaller than the system size but larger than the ion Larmor radius ($\rho_L$). While the ICM is often modeled as a fluid, its microphysical transport properties—specifically viscosity—are highly sensitive to collisionality and magnetic field geometry. Standard isotropic Spitzer viscosity predicts a Reynolds number ($Re$) that could be low enough to suppress turbulence, yet observational evidence from cold fronts (CFs) and velocity dispersions suggests the effective viscosity is much lower (1–20% of Spitzer).

Theoretical considerations indicate that in high-$\beta$ plasmas, the magnetic field induces anisotropic viscous stress (Braginskii viscosity). Furthermore, pressure anisotropy ($\Delta p = p_\perp - p_\parallel$) driven by adiabatic invariance can trigger kinetic plasma instabilities (firehose, mirror, ion-cyclotron). These instabilities scatter particles, effectively limiting the pressure anisotropy to marginally stable values and reducing the effective viscous stress. The interplay between Braginskii viscosity, magnetic field geometry, and these instability-driven limiters remains poorly constrained in the context of the turbulent, sloshing gas motions observed in cool-core clusters.

Methodology
The authors performed high-resolution, parallel magnetohydrodynamic (MHD) simulations using the FLASH 4.7.1 code to model a cool-core galaxy cluster ($M \sim 10^{15} M_\odot$, $T \sim 7$ keV) undergoing a merger with a subcluster. The simulation setup includes:

• Physics: The simulations solve the conservative MHD equations including the Braginskii viscous stress tensor. The pressure anisotropy is modeled via a "hard-wall" limiter approach, where $\Delta p$ is capped at marginal stability thresholds derived from kinetic theory.
• Limiting Prescriptions: Three scenarios were tested:
  1. Unlimited: No limits on $\Delta p$ (Braginskii viscosity acts freely, though the equations are ill-posed at kinetic scales).
  2. Limiters v1: $\Delta p$ is restricted by the oblique firehose ($\delta p \approx -1.4/\beta$) and mirror ($\delta p \approx 1/\beta$) instability thresholds.
  3. Limiters v2: $\Delta p$ is restricted by the parallel firehose ($\delta p \approx -2/\beta$) and ion-cyclotron ($\delta p \approx 0.5/\beta^{0.5}$) thresholds.
• Initial Conditions: Two sets of simulations were run with different initial plasma beta values ($\beta_{ini} = 100$ and $400$) to explore the effect of varying magnetic tension.
• Analysis: The study analyzes the suppression of Kelvin-Helmholtz instabilities (KHI) at cold fronts, the spatial distribution of pressure anisotropy, effective viscosity reduction factors, velocity amplitude spectra (3D and 2D projected), and energy transfer functions.

Key Contributions and Results

1. Suppression of Kelvin-Helmholtz Instabilities (KHI):
   Braginskii viscosity suppresses KHI at cold front surfaces, but the effect is modest compared to isotropic viscosity. The suppression depends heavily on the magnetic field strength and the limiting prescription. In simulations with stronger initial magnetic fields ($\beta_{ini}=100$), magnetic tension alone provides significant KHI suppression, making the addition of viscosity less visually distinct. In weaker field cases ($\beta_{ini}=400$), viscosity plays a larger role, though the "hard-wall" limiters reduce the effective viscosity, resulting in less suppression than the "Unlimited" case.

2. Effective Viscosity Reduction:
   The combined effect of magnetic field anisotropy and plasma instability limiters drastically reduces the effective viscosity.

   • Geometric Suppression: The tangled magnetic field geometry alone reduces the viscous flux relative to the isotropic case (factor $f_{geom}$).
   • Instability Suppression: The limiters further reduce the flux (factor $f_{limit}$).
   • Net Result: In the cluster core, the effective viscosity is reduced to $\sim 20\%$ of the Spitzer value in 57–75% of the volume, and to $\sim 50\%$ in 90–96% of the volume. This reduction is most pronounced in regions where the magnetic field is weaker and not yet amplified by shear.

3. Pressure Anisotropy and Instabilities:
   The sloshing motions generate significant pressure anisotropies ($|\delta p| \sim 0.03-0.07$ at $\beta \sim 100-1000$). In "Unlimited" simulations, a fraction of the volume (up to 30% for mirror/oblique firehose thresholds) exceeds instability limits. When limiters are applied, the pressure anisotropy is constrained to marginal stability, with a non-negligible fraction of the volume sitting exactly at the mirror or firehose thresholds. This has implications for cosmic-ray transport via pitch-angle scattering.

4. Velocity Spectra and Turbulence:

   • Spectral Slope: In inviscid simulations, the 3D velocity-amplitude spectrum follows a Kolmogorov-like slope ($k^{-1/3}$). The inclusion of Braginskii viscosity steepens the spectrum to $\sim k^{-2/3}$ to $k^{-1}$.
   • Interpretation: The authors suggest this steepening may result from the damping of turbulent motions by viscosity, or it may reflect the underlying large-scale sloshing motions that have not fully decayed into turbulence. The steepening is less pronounced when limiters are applied, as the effective viscosity is lower.
   • Projected Observables: Projected bulk velocity and velocity dispersion maps show that small-scale features are smoothed in viscous runs compared to inviscid ones, though large-scale distributions remain similar.

5. Magneto-Immutability:
   The simulations provide evidence for "magneto-immutable" dynamics. The probability distribution functions (PDFs) of the field-aligned rate of strain ($S$) are significantly narrower in Braginskii simulations than in inviscid MHD runs. This indicates a self-organization process where turbulent flows rearrange to minimize magnetic-field-aligned strain, further suppressing viscous dissipation beyond what geometry alone would predict. However, the authors note this may be a "partial" realization, as the energy transfer functions show non-negligible dissipation.

Significance and Claims
The paper claims that if viscous momentum transport in the ICM is controlled by pressure anisotropy and limited by plasma instabilities, its observational signatures are subtle.

• Cold Fronts: It is difficult to constrain the effective viscosity solely from the smoothness of cold fronts, as magnetic tension and anisotropic viscosity produce similar suppression of KHI.
• Velocity Fields: The most promising observational signature lies in the velocity field properties accessible via microcalorimeter instruments (e.g., XRISM). The steepening of the velocity-amplitude spectrum is a potential indicator of anisotropic viscosity, though the physical origin of the specific slope ($k^{-2/3}$) requires further investigation.
• Self-Organization: The results support the existence of magneto-immutable dynamics in cluster cores, which enhances viscosity suppression and may extend turbulent cascades to smaller scales than previously thought.

The authors conclude that while their simulations capture the essential physics of anisotropic transport and instability limiting, future work with higher resolution, driven turbulence in stratified boxes, and more advanced kinetic treatments (e.g., CGL-MHD) is necessary to fully resolve the origin of the spectral slopes and the extent of magneto-immutability.