Effects of Viscosity on Sloshing Cold Fronts in Galaxy Clusters

The viscous properties of the intracluster medium (ICM) remain poorly constrained. Cold fronts-sharp discontinuities formed during cluster mergers-offer a potential avenue to probe the effective viscosity of the ICM. Velocity shear across these fronts should generate Kelvin-Helmholtz instabilities (KHI), unless viscosity or magnetic tension suppresses them. We perform cluster merger simulations incorporating four ICM viscosity models: (A) inviscid, (B) isotropic Spitzer viscosity, (C) anisotropic Braginskii viscosity, and (D) Braginskii viscosity limited by microinstabilities. The isotropic Spitzer viscosity (case B) strongly suppresses KHI, producing smooth cold front surfaces, while the inviscid (A) and microinstability-limited (D) cases show prominent ripples. The Braginskii case (C) yields intermediate suppression. We also vary the plasma $β$ parameter ($β\approx$ 100 and 1600) to examine how a changing magnetic field strength affects the results. Stronger magnetic fields further suppress KHI, leading to smoother fronts and reduced differences between different viscosity models, while also widening the range of permitted pressure anisotropies when microinstability-based limiters are present. These results indicate that both viscosity and magnetic fields play crucial roles in stabilising sloshing cold fronts in galaxy clusters.

💡 Research Summary

This paper investigates how the effective viscosity of the intracluster medium (ICM) influences the morphology and stability of sloshing cold fronts in galaxy clusters. Cold fronts are sharp contact discontinuities that arise when a sub‑cluster perturbs the core of a larger, relaxed cluster, causing the low‑entropy gas to “slosh” in the gravitational potential well. The resulting shear flow across the front can trigger Kelvin‑Helmholtz instabilities (KHI), which would appear as ripples or eddies in high‑resolution X‑ray images. Whether such ripples are present, and how smooth the fronts appear, provides a diagnostic of the ICM’s microphysical properties, in particular its viscosity and magnetic field strength.

The authors perform a suite of three‑dimensional magnetohydrodynamic (MHD) simulations using the FLASH code with adaptive mesh refinement. The merger setup follows ZuHone et al. (2011): a primary cool‑core cluster of mass 1.25 × 10¹⁵ M⊙ and a sub‑cluster of 2.5 × 10¹⁴ M⊙ start 3 Mpc apart with an impact parameter of 0.5 Mpc and a relative velocity of 1466 km s⁻¹. The computational domain is a cube 8 Mpc on a side, refined down to 3.9 kpc in the central 150 kpc where sloshing occurs. The dark‑matter potential follows a Hernquist profile, while the gas is initialized with a cool core and a β‑model atmosphere.

Four viscosity prescriptions are explored: (A) inviscid, (B) isotropic Spitzer viscosity, (C) anisotropic Braginskii viscosity, and (D) Braginskii viscosity limited by kinetic micro‑instabilities (firehose and mirror). The isotropic Spitzer viscosity uses the classic μ ∝ T⁵⁄² / ln Λ scaling, yielding kinematic viscosities ν≈10²⁸–10³¹ cm² s⁻¹ in the core. The Braginskii model makes the viscous stress tensor depend on the magnetic‑field direction b̂, with parallel and perpendicular components differing by a factor of three. In model D the pressure anisotropy Δp = p⊥ − p∥ is constrained to lie within the firehose–mirror stability bounds (−2β ≤ Δp ≤ β), effectively capping the anisotropic viscosity when the plasma becomes unstable.

Magnetic fields are initialized as tangled, power‑law spectra ranging from 500 kpc down to 43 kpc, with two choices of plasma β (β = 100 and β = 1600) representing strong (≈5 μG) and weak (≈0.4 μG) fields, respectively. The β parameter is used throughout to quantify magnetic energy relative to thermal pressure.

The simulations are analyzed at two epochs after core passage: 2.2 Gyr (early sloshing) and 4 Gyr (more developed spiral patterns). The main diagnostics are the visual appearance of the cold‑front surface, the thickness of the density/temperature jump, and the presence of KHI‑like ripples with characteristic wavelengths 10–20 kpc.

Key findings:

1. Isotropic Spitzer viscosity (Model B) strongly suppresses KHI. The high ν reduces the Reynolds number dramatically, so the shear layer remains laminar. Cold‑front surfaces are smooth, with thicknesses ≤ 5 kpc and no visible ripples at either epoch, regardless of β. This reproduces the “perfectly sharp” fronts seen in early Chandra observations.

2. Inviscid (A) and micro‑instability‑limited Braginskii (D) produce prominent ripples. With essentially no effective viscosity (or a viscosity limited by the instability ceiling), the shear flow develops KHI modes with wavenumbers k≈0.02–0.05 kpc⁻¹. The fronts display undulations of 10–20 kpc amplitude, and the jump thickness expands to ≈10–15 kpc.

3. Anisotropic Braginskii viscosity (C) yields intermediate suppression. Because the viscous stress acts mainly along magnetic‑field lines, the degree of KHI damping varies locally with the angle between the shear flow and the field. In regions where the field is roughly parallel to the front, KHI is only weakly damped, producing modest ripples; where the field is perpendicular, suppression is comparable to the isotropic case. The net result is a front that is smoother than A/D but not as immaculate as B.

4. Magnetic field strength (β) modulates all effects. For β = 100 (strong field), magnetic tension alone can suppress KHI, reducing the differences between viscosity models. Moreover, the firehose/mirror stability bounds widen (since they scale with β), allowing larger pressure anisotropies before the micro‑instability limiter kicks in. Consequently, Model D behaves more like Model C, and even Model A shows fewer ripples. For β = 1600 (weak field), magnetic tension is negligible; KHI suppression relies almost entirely on viscosity, so the hierarchy A/D > C > B is most pronounced.

5. Temporal evolution. By 2.2 Gyr the first cold front forms; by 4 Gyr the spiral pattern is well developed. In Model B the front remains smooth throughout; in Model C the early ripples fade as the field becomes more aligned with the front; in Models A and D the ripples persist, indicating that once KHI is triggered it can survive for several gigayears in the absence of strong damping.

The authors compare these synthetic fronts with observed clusters. Some observed fronts (e.g., in the Perseus cluster) are extremely smooth, suggesting an effective viscosity of order 5–10 % of the Spitzer value, possibly combined with a moderately strong magnetic field (β ≈ 100–300). Other clusters show clear KHI eddies, implying lower viscosity or weaker fields. The study therefore constrains the ICM viscosity to be sub‑Spitzer but not negligible, and highlights that magnetic tension can mask the viscous signature, making it essential to measure both magnetic field strength (e.g., via Faraday rotation or radio polarization) and front morphology.

In conclusion, the paper demonstrates that (i) isotropic Spitzer viscosity can alone produce the smooth cold fronts seen in some clusters, (ii) anisotropic Braginskii viscosity leads to a direction‑dependent suppression of KHI, (iii) kinetic micro‑instabilities limit the effective anisotropic viscosity, allowing KHI to develop, and (iv) strong magnetic fields both directly suppress KHI and broaden the allowed pressure‑anisotropy range, reducing the observable differences between viscosity models. These results provide a roadmap for interpreting high‑resolution X‑ray observations and for future work that couples kinetic plasma physics with cluster‑scale MHD simulations.