Effervescent heating: constraints from nearby cooling flow clusters observed with XMM-Newton

Rocco Piffaretti,Jelle Kaastra

DOI: https://doi.org/10.1051/0004-6361%3A20054365

2006-05-16

Abstract:We have used deprojected radial density and temperature profiles of a sample of 16 nearby CF clusters observed with XMM-Newton to test whether the effervescent heating model can satisfactorily explain the dynamics of CF clusters. For each cluster we derived the required extra heating as a function of cluster-centric distance for various values of the unknown parameters $\dot M$ (mass deposition rate) and $f_c$ (conduction efficiency). We fitted the extra heating curve using the AGN effervescent heating function and derived the AGN parameters $L$ (the time-averaged luminosity) and $r_0$ (the scale radius where the bubbles start rising in the ICM). While we do not find any solution with the effervescent heating model for only one object, we do show that AGN and conduction heating are not cooperating effectively for half of the objects in our sample. For most of the clusters we find that, when a comparison is possible, the derived AGN scale radius $r_0$ and the observed AGN jet extension have the same order of magnitude. The AGN luminosities required to balance radiative losses are substantially lowered if the fact that the AGN deposits energy within a finite volume is taken into account. For the Virgo cluster, we find that the AGN power derived from the effervescent heating model is in good agreement with the observed jet power.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to verify whether the **AGN (Active Galactic Nucleus) bubble heating model** can explain the dynamic phenomena of cooling - flow galaxy clusters (CF clusters). Specifically, the researchers used the de - projected radial density and temperature distributions of 16 nearby cooling - flow galaxy clusters observed by XMM - Newton to test whether the AGN bubble heating model can satisfactorily explain the dynamic behaviors of these galaxy clusters.

### Research Background

In the core regions of many galaxy clusters, due to the short radiation cooling time, the plasma can release energy equal to its thermal energy in less than one billion years, causing the gas to gradually condense and eventually form a cooling flow (CF). However, the early - estimated condensation mass deposition rate (about \(10^2 - 10^3 M_\odot \text{yr}^{-1}\)) indicates that the cores of cooling - flow galaxy clusters should contain a large amount of cold gas, but actual observations have not found the existence of these cold gases.

With the observations of XMM - Newton and Chandra satellites, our understanding of cooling - flow galaxy clusters has been significantly improved. In particular, the X - ray spectrum failed to detect gas emission lines below 2 keV, which refutes the standard cooling - flow model. In addition, the observations also found that the mass deposition rate is much lower than previous estimates, and the temperature in the central region drops to about one - third of the average temperature, indicating that the gas cannot cool to a lower temperature.

### Research Questions

Based on the above background, the researchers proposed the following questions:
1. **Can the AGN bubble heating model explain the dynamic behaviors of cooling - flow galaxy clusters?**
2. **Can AGN heating and thermal conduction effectively cooperate under certain conditions to balance radiation losses?**

To answer these questions, the researchers derived AGN parameters (such as the time - averaged luminosity \(L\) and the scale radius \(r_0\) at which the bubble begins to rise) by fitting an additional heating curve using the AGN bubble heating function, and analyzed the influence of different parameters (such as the mass deposition rate \(\dot{M}\) and the thermal conduction efficiency \(f_c\)) on the results.

### Key Formulas

1. **AGN Heating Function**:
\[
H_{\text{AGN}} = -h(r) \left(\frac{p}{p_0}\right)^{\frac{\gamma_b - 1}{\gamma_b}} \frac{1}{r} \frac{d \ln p}{d \ln r}
\]
where
\[
h(r) = \frac{L}{4\pi r^2} \left(1 - e^{-\frac{r}{r_0}}\right)^{q - 1}
\]
and
\[
q = \int_0^{+\infty} \left(\frac{p}{p_0}\right)^{\frac{\gamma_b - 1}{\gamma_b}} \frac{1}{r} \frac{d \ln p}{d \ln r} \left(1 - e^{-\frac{r}{r_0}}\right) dr
\]

2. **Thermal Conduction Heating**:
\[
\varepsilon_{\text{cond}} = \frac{1}{r^2} \frac{d}{dr} \left(r^2 \kappa \frac{dT}{dr}\right)
\]
where \(\kappa\) is the thermal conductivity, and for ionized plasma, \(\kappa\) is the Spitzer thermal conductivity:
\[
\kappa = \kappa_S = 1.84\times 10^{-5} T^{5/2} \ln \Lambda \, \text{erg cm}^{-1} \text{s}^{-1} \text{K}^{-1}
\]

Through these formulas, the researchers calculated the required additional heating and fitted it with the observed data, thereby evaluating the AGN bubble heating model.