Are X-ray Atmospheres Heated by Turbulent Dissipation? XRISM Constraints

1. Problem Statement

The central problem addressed is the mechanism by which radio jets from active galactic nuclei (AGN) prevent the radiative cooling of hot intracluster medium (ICM) in galaxy clusters. While it is established that jet power roughly matches the X-ray cooling luminosity, the specific mechanism of energy dissipation remains debated.

• The Hypothesis: The paper investigates whether the dissipation of turbulence injected by rising bubbles and jets is sufficient to offset radiative cooling flows.
• The Challenge: Previous observations (e.g., XMM-Newton, Hitomi) provided upper limits on velocity dispersions ($\sigma_v$), but the new XRISM observatory offers high-resolution spectroscopy (5 eV near the 6.7 keV iron line) to measure $\sigma_v$ with precision ($\sim$ tens of km/s). The authors aim to determine if the measured turbulent velocities are physically capable of heating the entire cooling volume, or if they are merely bulk motions that fail to thermalize efficiently.

2. Methodology

The authors employ a multi-faceted approach combining observational data analysis with phenomenological and physically motivated modeling.

• Data Sample: The study utilizes XRISM data for a sample of 10 galaxy clusters (including Perseus, Virgo/M87, Hydra A, Abell 2029, Cygnus A, etc.). They analyze central velocity dispersions ($\sigma_v$) and radial profiles where available.
• Correlation Analysis: The authors test for trends between jet power ($P_{jet}$) and:
  • Central atmospheric velocity dispersion ($\sigma_v$).
  • The ratio of kinetic to thermal energy ($E_k/E_{th}$).
  • Specific kinetic energy (energy per unit gas mass).
• Modeling Frameworks:
  1. Phenomenological Models: They test idealized scenarios where turbulent power ($P_{turb} \propto M(r)\sigma_v^3 / l$) balances X-ray cooling luminosity ($L_X$). They vary the injection scale ($l$) and velocity dispersion ($\sigma_v$) to see if equilibrium is possible.
  2. HAK19 Model: A model by Zhoolideh Haghighi et al. (2019) where turbulent heating timescales are anchored to the sound speed ($c_s$) via a viscosity parameter $\alpha$, assuming a radially rising injection scale ($l \propto R$).
  3. Buoyancy-Driven Model: A new physically motivated model where bubbles rise at a terminal speed ($v_{term} = \epsilon c_s$). This model assumes the injection scale grows with altitude (as bubbles expand) and derives the required turbulent power to balance cooling.
• Timescale Analysis: The authors compare the turbulent dissipation timescale ($\tau_{turb} \sim l/\sigma_v$) against the radiative cooling time ($\tau_{rad}$) and the jet duty cycle. They also consider the constraints imposed by a stably stratified atmosphere (Brunt-Väisälä frequency) on the size of radial eddies.

3. Key Results

A. Lack of Scaling Trends

• No Correlation: Over nearly four decades of jet power ($10^{42}$ to $10^{46}$ erg s$^{-1}$), there is no trend between jet power and central velocity dispersion, kinetic-to-thermal energy ratio, or specific kinetic energy.
• Implication: Jets appear to disperse energy "gently," imparting a roughly constant amount of kinetic energy per gram of gas regardless of the total jet power. High $\sigma_v$ in systems like Cygnus A and M87 likely reflects unresolved bulk motions rather than fully developed turbulence.

B. Radial Profiles and Injection Scales

• Flat/Declining Profiles: In most systems (Hydra A, Perseus, M87), $\sigma_v$ is flat or declines with radius, contrary to models requiring rising velocity dispersions to maintain heating.
• Injection Scale Constraints: To balance cooling, models require injection scales ($l$) that are often much smaller than the bubble sizes (a few kpc) and dissipation timescales ($\tau_{turb}$) that are extremely short ($\lesssim 10^7$ yr) in inner regions.
• Diffusion Limit: Turbulence generated centrally cannot diffuse fast enough to fill the cooling volume. The diffusion coefficient for jet-driven turbulence is too low; turbulence must be injected locally throughout the volume, which is difficult for centrally located jets.

C. Model Performance

• Abell 2029: The buoyancy model (bubbles rising at $\sim 0.5 c_s$) successfully balances cooling. This is aided by the system's high sound speed and a single, flat $\sigma_v$ measurement across the cooling volume.
• Hydra A: The buoyancy model works well if bubbles rise at $0.2 - 0.5 c_s$. The flat $\sigma_v$ profile supports this. However, the HAK19 model (requiring rising $\sigma_v$) is inconsistent with observations.
• Perseus & Virgo (M87): The buoyancy model fails to balance cooling in the outer regions ($>20$ kpc).
  • The observed $\sigma_v$ declines with radius.
  • To maintain equilibrium, the model requires bubble speeds exceeding the sound speed ($\epsilon > 1$) or injection scales of only a few kpc, which are physically implausible for subsonic bubbles.
  • The HAK19 model also fails because it predicts rising $\sigma_v$ profiles not seen in data.

D. Energy Budget

• Only a fraction (roughly 20–50% at most) of the jet enthalpy can be converted into turbulent kinetic energy.
• In Perseus and Virgo, turbulent dissipation likely struggles to offset cooling except perhaps in the innermost regions.

4. Significance and Conclusions

• Turbulence is Likely Insufficient: The paper concludes that while turbulent dissipation is a plausible heating mechanism in specific, high-sound-speed systems (like Abell 2029) or inner regions, it struggles and probably fails to balance cooling in the vast majority of cluster atmospheres (specifically Perseus and Virgo).
• Physical Constraints: The combination of low velocity dispersions, anisotropic injection (jets/bubbles), long turbulent diffusion timescales, and the requirement for short duty cycles presents severe challenges for jet-driven turbulence heating.
• Alternative Mechanisms: The failure of pure turbulent dissipation models suggests that other mechanisms (shocks, sound waves, cosmic rays, or mixing) or a combination thereof are likely required to achieve thermal balance.
• Future Work: The authors emphasize that a larger sample of spatially resolved cluster atmospheres is needed to definitively confirm these findings. The current results highlight the power of XRISM to constrain the physics of AGN feedback by ruling out simple turbulent heating scenarios.

Summary of Key Equations Used

• Turbulent Power: $P_{turb} \simeq \frac{3}{2} \frac{M(r)\sigma_v^3}{l}$
• Buoyancy Model Power: $P_{buoy} \simeq \frac{3}{2} \frac{M(r)\sigma_v^2 \epsilon c_s}{R}$ (where $\epsilon = v_{term}/c_s$)
• Dissipation Timescale Limit (Stratified): $t_{turb} \gtrsim \sqrt{\frac{\gamma}{\eta}} \frac{R}{v_K}$

The paper represents a critical step in moving from "jet power equals cooling power" to a detailed understanding of how that energy is distributed and dissipated, ultimately casting doubt on turbulence as the sole or primary heating agent in many major clusters.