ABSTRACT
We present the results of a Chandra study of the Hydra A galaxy cluster, where a powerful active galactic nucleus (AGN) outburst created a large-scale cocoon shock. We investigated possible azimuthal variations in shock strength and shape, finding indications for a weak shock with a Mach number in the range ∼1.2–1.3. We measured the temperature change across the shock front. However, the detection of a temperature rise in the regions immediately inside of the front is complicated by the underlying temperature profile of the cluster atmosphere. We measured the global temperature profile of the cluster up to 700 kpc, which represents the farthest measurement obtained with Chandra for this cluster. A "plateau" in the temperature profile in the range ∼70–150 kpc indicates the presence of cool gas, which is likely the result of uplift of material by the AGN outburst. After masking the cool filaments visible in the hardness ratio map, the plateau disappears and the temperature profile recovers a typical shape with a peak around 190 kpc, just inside the shock front. However, it is unlikely that such a temperature feature is produced by the shock as it is consistent with the general shape of the temperature profiles observed for relaxed galaxy clusters. We studied the spectral properties of the cool filaments finding evidence that ∼1011M☉ of low-entropy material has been dredged up by the rising lobes from the central 30 kpc to the observed current position of 75–150 kpc. The energy required to lift the cool gas is ≳ 2.2 × 1060 erg, which is comparable to the work required to inflate the cavities and is ∼25% of the total energy of the large-scale shock. Our results show that the AGN feedback in Hydra A is acting not only by directly heating the gas, but also by removing a substantial amount of potential fuel for the supermassive black hole.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
Observational and theoretical evidence has been growing in the past decade in favor of the existence in galaxy clusters of a feedback mechanism that prevents cool cores from establishing "cooling flows" at the rates predicted by earlier X-ray observations (e.g., see Peterson & Fabian 2006 for a review). The dominant cD galaxies, which are present at the cluster center in all cool core clusters, host the most massive black holes in the local universe and usually show nuclear activity. They accordingly provide a natural feedback mechanism for the regulation of the cooling process. Feedback is also required to suppress the overproduction of massive galaxies predicted by dark-matter-only simulations and to break the self-similarity of clusters (e.g., Benson et al. 2003). The nature of this feedback, vital to our understanding of galaxy and structure evolution, is one of the most important unresolved questions in extragalactic astronomy. Based on the detection of cavities and active galactic nucleus (AGN)-driven shocks, the primary source of feedback in cluster has been identified as radio galaxies acting through outbursts and accompanying energy injection, likely intermittent, from the central AGN (e.g., for a review, see McNamara & Nulsen 2007 and references therein). In galaxy clusters, where cooling rates should be highest, the current generation of X-ray observatories Chandra and XMM-Newton have shown that there is not a significant amount of gas cooling below about one-third of its virial temperature (Peterson et al. 2003; Kaastra et al. 2004). Images from these telescopes also reveal highly disturbed structures in the cores of many clusters, including shocks, cavities, and sharp density discontinuities. At radio wavelengths, it is clear that AGN jets are the cause of many of these disturbances. The incidence and variety of bubbles, cavities, shocks, and ripples observed both in the radio and in X-rays in galaxy clusters provide direct evidence of the widespread presence of AGN-driven phenomena (see, e.g., Fabian et al. 2003 for a discussion of the Perseus cluster properties in terms of AGN-generated viscously damped sound waves; Nulsen et al. 2005 for shock heating in Hydra A; Gitti et al. 2007 for a study of giant cavities created by the most powerful AGN outburst currently known; Bîrzan et al. 2008 and Diehl et al. 2008 for a survey of cavities and the implied cavity heating rates). Such AGN feedback has a wide range of impacts, from the formation of galaxies, to the explanation of the observed Mbh–σ relation (which indicates a causal connection or feedback mechanism between the formation of bulges and their central black holes, e.g., Magorrian et al. 1998), to the regulation of cool cores. In most cases, the energy introduced by the AGN is more than sufficient to counteract putative cooling flows (Bîrzan et al. 2004, 2008; Rafferty et al. 2006). However, the details of how the feedback loop operates are still unknown. Only by studying striking examples of interaction between the central radio galaxy and the intracluster medium (ICM) can we understand why cooling and star formation still proceeds at a reduced rate, and this is likely to reveal the coupling between AGN feedback and the ICM.
The galaxy cluster Hydra A has a well-known, large-scale system of X-ray cavities embedded in a "cocoon" shock surrounding the central, powerful radio source (McNamara et al. 2000; Nulsen et al. 2005). It is considered one of the prototypes of cool core clusters with cavities, which has served as an early test of the AGN feedback paradigm, and it has been extensively studied both in the radio and X-rays (Taylor et al. 1990; McNamara et al. 2000; David et al. 2001; Nulsen et al. 2002, 2005; Lane et al. 2004; Wise et al. 2007; Simionescu et al. 2009a, 2009b; Kirkpatrick et al. 2009). By analyzing the archival ∼200 ks Chandra exposure in this paper, we study the azimuthal properties of the large-scale shock and attempt to measure the temperature jump associated with the shock in different directions. We also perform a detailed spectral analysis of the cool X-ray filaments extending out to 150 kpc, finding evidence for extensive mass dredge-up from the central 30 kpc, which affects the global temperature profile of the cluster.
With H0 = 70kms−1Mpc−1 and ΩM = 1 − ΩΛ = 0.3, the luminosity distance to Hydra A (z = 0.0538) is 240 Mpc and 1'' corresponds to 1.05 kpc in the rest frame of the cluster.
2. CHANDRA OBSERVATIONS AND DATA REDUCTION
Hydra A has been imaged four times by Chandra ACIS for a total exposure of 240 ks. The two shorter exposures (ObsIDs 575 and 576, which collectively comprise only ∼17% of the total exposure) were taken early in the Chandra mission when the ACIS detector was operated at the higher focal plane temperature of −110 °C. Due to the higher quality calibration at −120 °C, the spectral analysis presented below was performed on the two longer, more recent exposures: ∼97 ks collected on 2004 January 13 (ObsID 4969) and ∼99 ks collected on 2004 October 22 (ObsID 4970) with ACIS-S. We use data from the S3 and S2 CCDs to study the central part of the cluster emission, where the cavity system and the radio source are located, and from the I2 and I3 CCDs to measure the temperature and surface brightness in the cluster outskirts. Each data set was individually reprocessed with CIAO version 4.1 using CALDB 4.1.0 and corrected for known time-dependent gain and charge transfer inefficiency problems following techniques similar to those described in the Chandra analysis threads.6 Screening of the event files was also applied to filter out strong background flares. Blank-sky background files, filtered in the same manner as in the Hydra A image and normalized to the count rate of the source image in the 10–12 keV band, were used for background subtraction. We identified and removed the point sources using the CIAO task WAVDETECT, with the detection threshold set to the default value of 10−6. The final, combined exposure time for the two data set is 174.2 ks.
3. THE LARGE-SCALE SHOCK
The observed feature interpreted as a shock front in the X-ray surface brightness of Hydra A surrounds the low-frequency radio lobes, and the correspondence between their shapes supports the interpretation of a cocoon shock of the radio source. Nulsen et al. (2005) estimated an age of the outburst to be ts = 1.4 × 108 yr with a total energy of Es = 9 × 1060 erg. The shock front is clearly aspherical with a complicated three-dimensional shape affected by projection effects. In particular, the Hydra A radio source and cavity system is inclined at about 40° to the plane of the sky, with the northern side lying closer to us (Taylor 1996; Lane et al. 2004; Wise et al. 2007). The cavity and shock geometry as well as the outburst history and energetics have been studied in detail in Wise et al. (2007) and Simionescu et al. (2009a).
3.1. Azimuthal Variations in Shock Strength
We have investigated the azimuthal variations in shock strength and shape by studying surface brightness profiles in different directions. In particular, we have divided the northern cluster semicircle into six sectors of 30° each, and the southern cluster semicircle, where the shock is less evident in the image, into two quadrants of 90° each. Starting from the west with position angle (P.A.) = 0° and counting counterclockwise, the sectors are labeled as 0–30, 30–60, 60–90, 90–120, 120–150, 150–180, 180–270, and 270–360 (see Figure 1). We extracted the background-subtracted, exposure-corrected surface brightness profile along each sector in the energy range 0.5–2.0 keV. The center of the annular regions used to extract each profile was chosen in order to best match the curvature of the radial bins with the shape of the shock front as seen in the image (Figure 1). We then performed fits of a broken power-law density model to each surface brightness profile. This analysis assumes that the radius of curvature in the plane of the sky is the same as that parallel to the line of sight. A hydrodynamic model for the shock was made by initiating an explosion at the center of a hydrostatic, isothermal atmosphere with a power-law density profile (see Nulsen et al. 2005 for more details on the shock model). The power-law index for the density profile of the unshocked gas was determined from the broken power-law fit to the surface brightness profile. The surface brightness profiles and the best fits of broken power-law density models in the various sectors are shown in Figure 2, where the radial axes indicate the distance from the center of curvature of the annular regions used to extract the profiles. However, in the following discussion we indicate the position of the shock in each sector (referred to as the "shock radius") in terms of the distance of the front at mid-aperture of the sector from the cluster center. The best-fit model for each sector is summarized in Table 1 and the corresponding shock front is shown by red arcs in Figure 1. We also investigated the effects of varying the radial binning and center of curvature of the annular regions used to extract the profile. We found that the fit results (i.e., shock radius and Mach number, ) do not depend strongly on the particular choice of the extraction region nor on the radial range of the fit, with systematic variations ≲ 5% in each sector.
Download figure:
Standard image High-resolution imageTable 1. Shock Properties
Sector | Shock Radius | Radius Variation | Density Jump | Mach Number |
---|---|---|---|---|
θ1–θ2 (°) | ('') | – ('') | ||
0–30 | 208 | 201–222 | 1.33 | 1.22 |
30–60 | 257 | 222–303 | 1.45 | 1.30 |
60–90 | 332 | 314–349 | 1.33 | 1.23 |
90–120 | 371 | 364–342 | 1.38 | 1.26 |
120–150 | 301 | 342–271 | 1.46 | 1.32 |
150–180 | 260 | 264–265 | 1.40 | 1.27 |
180–270 | 285 | 260–253 | 1.15 | 1.10 |
270–360 | 205 | 233–202 | 1.35 | 1.23 |
Notes. Results of the fits of a broken power-law density model to the surface brightness profiles extracted along different sectors. The sector aperture (from P.A. θ1 to θ2) is indicated in the first column. The shock radius at the mid-angle of the sector () is indicated in the second column, whereas the shock radii at the starting () and ending () angles are indicated in the third column. Such radial distances are measured from the cluster center. The corresponding density jump and Mach number () are shown in the fourth and fifth columns, respectively. Best-fit statistical errors are on average ∼5%.
Download table as: ASCIITypeset image
The shock is clearly visible as a surface brightness jump in all sectors but the 180–270 sector (i.e., the southeast (SE) quadrant), with a radius varying between ∼205'' in the east–west direction and ∼365'' in the north–south direction. The Mach number varies between 1.20 and 1.32 in the sectors where the density jump is detected, whereas to the southeast (sector 180–270) where the surface brightness discontinuity is less evident we estimate . These results are consistent with Nulsen et al. (2005), and partially in agreement with Simionescu et al. (2009a), who report a detection of the shock front also in the south direction.
3.2. Temperature Variations Across the Shock Front
Our shock model predicts the emission-weighted temperature to rise across the front by ∼8%, 10%, 15%, and 20% for Mach numbers of 1.18, 1.23, 1.33, and 1.40, respectively, reaching its peak at a distance of ∼5%–10% of the shock radius behind the shock and declining below the unshocked temperature values inside radial distances of ∼25% of the shock radius. We used this information to optimize our selection of the pre-shock and post-shock regions in each sector and then extracted the spectra in these regions using the SPECEXTRACT task, which also computes the corresponding event-weighted response matrices. Spectral fitting to a single absorbed apec model was performed in XSPEC version 12.3.1 in the 0.5–8.0 keV energy range. Abundances were measured relative to the abundances of Anders & Grevesse (1989) and a galactic hydrogen column of 4.68 × 1020 cm−2 (Dickey & Lockman 1990) was assumed. As a general result, the post-shock regions are found to be hotter than the corresponding pre-shock regions, although due to the large error bars the temperature of the pre-shock and post-shock gas is consistent. Some examples of these temperature measurements are presented in Section 6.1. A further attempt to bin sectors together in order to increase the statistics did not lead to any clear detections of a temperature jump across the shock front (see Section 6.1). However, previous measurements (David et al. 2001; Simionescu et al. 2009a) indicate that the underlying global temperature profile has a peak close to the inner edge of the shock front, thus complicating the detection of a temperature rise due to the shock itself. To investigate this in more detail, we study the properties of the azimuthally averaged cluster temperature profile in Section 4.
4. GLOBAL CLUSTER TEMPERATURE PROFILE
The azimuthally averaged gas temperature profile was derived from the ACIS-S and ACIS-I data by extracting spectra in the annular regions indicated in Figure 3, and is shown in the left panel of Figure 4. The annular bins and temperature measurements are detailed in Table 2. The two outer temperature points represent the farthest temperature measurements from the cluster center obtained with Chandra data at present. We note that the temperature profile peaks around 180'', just inside the inner edge of the shock front. However, as discussed in Section 6, it is unlikely that this temperature feature is produced by the shock as it is consistent with the general shape of temperature profiles observed for relaxed galaxy clusters.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageTable 2. Temperature Profile Derived by Spectral Fitting
Bin | Radius | kT | χ2/dof | kTmask | χ2/dof |
---|---|---|---|---|---|
No. | ('') | (keV) | (keV) | ||
1 | 0–13 | 2.93+0.04 − 0.04 | 763/560 | 3.11+0.07 − 0.07 | 453/430 |
2 | 13–21 | 3.07+0.04 − 0.04 | 638/590 | 3.26+0.06 − 0.06 | 520/466 |
3 | 21–29 | 3.29+0.05 − 0.04 | 607/598 | 3.34+0.05 − 0.05 | 589/566 |
4 | 29–37 | 3.39+0.05 − 0.05 | 663/600 | 3.55+0.07 − 0.07 | 628/565 |
5 | 37–46 | 3.61+0.06 − 0.06 | 681/622 | 3.76+0.07 − 0.07 | 605/584 |
6 | 46–58 | 3.60+0.06 − 0.06 | 644/607 | 3.73+0.07 − 0.07 | 615/578 |
7 | 58–72 | 3.62+0.07 − 0.06 | 737/623 | 3.81+0.07 − 0.07 | 618/593 |
8 | 72–89 | 3.51+0.07 − 0.07 | 656/627 | 3.73+0.10 − 0.10 | 512/502 |
9 | 89–107 | 3.51+0.07 − 0.07 | 748/629 | 3.85+0.10 − 0.10 | 538/513 |
10 | 107–126 | 3.52+0.07 − 0.07 | 698/634 | 3.82+0.10 − 0.10 | 507/508 |
11 | 126–145 | 3.58+0.07 − 0.07 | 713/648 | 3.78+0.11 − 0.11 | 491/511 |
12 | 145–168 | 3.79+0.07 − 0.07 | 744/674 | 3.91+0.08 − 0.08 | 672/643 |
13 | 168–201 | 4.02+0.08 − 0.08 | 719/709 | 4.02+0.08 − 0.08 | 719/709 |
14 | 201–257 | 3.91+0.08 − 0.08 | 1125/1124 | 3.91+0.08 − 0.08 | 1125/1124 |
15 | 257–339 | 3.70+0.11 − 0.11 | 1252/1177 | 3.70+0.11 − 0.11 | 1252/1177 |
16 | 339–497 | 3.58+0.19 − 0.19 | 1140/1150 | 3.58+0.19 − 0.19 | 1140/1150 |
17 | 497–812 | 2.96+0.30 − 0.29 | 1717/1676 | 2.96+0.30 − 0.29 | 1717/1676 |
Notes. Results of the spectral fitting in concentric 360° annular regions (shown in Figure 3) in the 0.5–8.0 keV energy range using the XSPEC wabs×apec model. The absorbing column density is fixed to the Galactic value (NH = 4.68 × 1020 cm−2), while the temperature (in keV) and metallicity (not reported here) are left as free parameters. Error bars are at the 90% confidence levels on a single parameter of interest. The values of kTmask are measured in the same concentric annular regions after excluding the cool filaments (i.e., by fitting spectra extracted after masking the blue sectors shown in Figure 5, left panel). See Section 4 for details.
Download table as: ASCIITypeset image
On the other hand, a "plateau" is notable in the temperature profile indicating the presence of cool gas in the range of radius ∼70''–150''. This peculiar feature has not been noted previously in the literature. We thus investigated it in more detail thanks to the deep, high-quality Chandra exposure.
4.1. Evidence for Cool Filaments
In order to investigate the origin of the plateau seen in the temperature profile between ∼70'' and 150'', we made several attempts to identify regions of cool gas. Our approach was to measure a new, "undisturbed" temperature profile by masking such regions and then compare it with the global temperature profile measured above. We first masked the cavity regions (as indicated in Wise et al. 2007), finding systematically higher temperatures by ∼0.1–0.2 keV in all bins (except for the very central bin, where the temperature is lower by ∼0.15 keV) up to ∼200''. The temperature profile measured by masking the cavities is therefore shifted up, maintaining its global shape. In particular, the plateau in the temperature profile is still present, indicating that the cavities are not a tracer of the entire amount of cool gas. We then masked the ∼1' bright filament stretching from the inner cavity to the center of a larger outer cavity in the northeast, which was noted by Nulsen et al. (2005). As in the previous attempt, the plateau in the temperature profile remained, indicating that the known filament is not the only repository of the cool gas in Hydra A.
As a diagnostic of the presence of cool gas we finally used the hardness ratio map shown in Figure 5, which was obtained by dividing a smoothed 1.5–7.5 keV image by a smoothed 0.3–1.5 keV image. The dark regions are indicative of low gas temperatures. They agree with the arm-like structures of cooler gas extending toward the north and the south in the temperature map of Simionescu et al. (2009b). Based on a visual inspection of the hardness ratio map, we selected a simple combination of sectors reproducing the shape of the supposedly cool filaments (see the left panel of Figure 5). We then excluded these regions and extracted the spectra in the same radial bins as above (see Figure 3), and generated the temperature profile shown in red in the right panel of Figure 4 (see Table 2). As evident from the comparison of the two profiles, the plateau has largely been removed and the temperature profile is typical of cool core clusters (see also Section 6). This clearly indicates that the masked regions contain cool gas. A detailed spectral analysis of the cool filaments is presented in Section 5.
Download figure:
Standard image High-resolution image5. SPECTRAL PROPERTIES OF THE COOL GAS
We investigate here the spectral properties of the gas which produces the plateau seen in the global cluster temperature profile. In particular, we focus on the gas located in the range of projected distances from the center of 72''–145'' (∼76–152 kpc), which corresponds to the radial range where the plateau in the temperature profile is most evident (i.e., bin numbers 8–11; see Figure 4 and Table 2). We divided the annulus from 72''–145'' into eight sectors, each having an angular width of 45°, obtaining the regions labeled as WNW (west–northwest), NNW (north–northwest), NNE (north–northeast), ENE (east–northeast), ESE (east–southeast), SSE (south–southeast), SSW (south–southwest), and WSW (west–southwest) in the right panel of Figure 5. We extracted the spectra in these sectors and compared two different spectral models. The "1T model" is the absorbed apec model already used above to derive the global temperature profile. The free parameters are the temperature, kT, the metallicity, Z, and the normalization (emission measure, EM). The "2T model" includes a second thermal emission component (apec+apec) and has two additional free parameters: the temperature, kT2, and the normalization, EM2, of the second component (the metallicities of the two components are linked).
The best-fitting parameter values and 90% confidence ranges derived from the fits to the annular spectra in sectors are summarized in Table 3. Although the improvement of adding a second thermal component is formally significant according to the F-test, our results show that the quality of the data is not generally sufficient to demand a model more complex than the 1T model. In fact, in most sectors, the 1T model already produces a very good fit (reduced χ2ν ∼ 1) and therefore a more complicated model appears unnecessary. We also note that the second thermal component is poorly constrained, with temperature errors ≳25% and up to 300%. Only in the SSE sector is the reduced chi squared of the 1T model unacceptable at 90% significance, and the statistical improvement obtained by introducing an additional emission component compared to the single-temperature model is the most significant according to the F-test. The improvement of the 2T model over the 1T model in this sector is also evident from the residuals of the fits in Figure 6. We can therefore conclude that, confirming the hardness ratio map (Figure 5, right panel), we find spectral evidence for multiphase gas in the SSE sector with a hot component at 4.92+0.97 − 0.74 keV and a cool component at 1.59+0.11 − 0.23 keV. Assuming that the two spectral phases are in pressure equilibrium in the same volume, the ratio of the volumes they occupy is estimated as V1/V2 = (EM1/EM2) · (kT1/kT2)2, so the filling factor of the cool gas is ∼0.04.
Download figure:
Standard image High-resolution imageTable 3. Spectral Analysis of the Cool Filaments
1T Model | 2T Model | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Sector | kT (keV) | EM (× 10−3) | Z (solar) | χ2/dof [χ2ν] | kT1 (keV) | EM1(× 10−3) | kT2(keV) | EM2(× 10−3) | Z (solar) | χ2/dof | F stat |
WNW | 3.87+0.10 − 0.10 | 1.60+0.03 − 0.03 | 0.37+0.05 − 0.05 | 505.1/521 [0.97] | 2.17+0.48 − 0.47 | 0.73+0.28 − 0.45 | 6.40+1.89 − 1.74 | 0.94+0.43 − 0.28 | 0.28+0.06 − 0.05 | 484.6/519 | 11.0 |
NNW | 3.24+0.06 − 0.07 | 1.79+0.03 − 0.03 | 0.36+0.04 − 0.04 | 524.9/525 [1.00] | 2.05+0.17 − 0.25 | 1.08+0.24 − 0.34 | 6.74+4.01 − 1.75 | 0.80+0.31 − 0.35 | 0.26+0.04 − 0.04 | 467.6/523 | 32.0 |
NNE | 3.14+0.08 − 0.08 | 1.23+0.03 − 0.03 | 0.41+0.05 − 0.05 | 483.0/460 [1.05] | 2.12+0.23 − 0.21 | 0.81+0.06 − 0.29 | 6.84+1.48 − 2.20 | 0.48+0.26 − 0.11 | 0.30+0.05 − 0.05 | 451.9/458 | 15.8 |
ENE | 3.32+0.08 − 0.08 | 1.15+0.02 − 0.02 | 0.42+0.06 − 0.05 | 455.5/455 [1.00] | 2.24+0.37 − 0.47 | 0.72+0.20 − 0.47 | 6.84+4.62 − 2.72 | 0.48+0.50 − 0.20 | 0.34+0.06 − 0.05 | 433.6/453 | 11.4 |
ESE | 3.45+0.10 − 0.10 | 1.36+0.03 − 0.03 | 0.37+0.05 − 0.05 | 502.2/483 [1.04] | 2.16+0.86 − 0.48 | 0.48+0.89 − 0.35 | 4.49+14.9 − 0.61 | 0.92+0.33 − 0.66 | 0.32+0.05 − 0.05 | 491.2/481 | 5.4 |
SSE | 3.34+0.07 − 0.07 | 1.45+0.03 − 0.03 | 0.38+0.05 − 0.05 | 621.8/491 [1.27] | 1.59+0.11 − 0.23 | 0.44+0.18 − 0.21 | 4.92+0.97 − 0.74 | 1.05+0.17 − 0.15 | 0.28+0.06 − 0.05 | 533.9/489 | 40.3 |
SSW | 4.01+0.12 − 0.12 | 1.18+0.02 − 0.02 | 0.35+0.06 − 0.06 | 458.0/472 [0.97] | 2.39+0.65 − 0.90 | 0.67+0.17 − 0.53 | 8.60+4.98 − 4.12 | 0.58+0.50 − 0.22 | 0.28+0.07 − 0.07 | 438.2/470 | 10.6 |
WSW | 3.92+0.11 − 0.11 | 1.37+0.03 − 0.03 | 0.31+0.06 − 0.05 | 439.4/494 [0.89] | 2.93+0.33 − 1.14 | 1.17+0.07 − 0.85 | 54.4+/ − / | 0.36+0.12 − 0.12 | 0.25+0.06 − 0.05 | 419.5/492 | 11.7 |
Notes. Results of 1T model and 2T model spectral fitting in the 0.5–8.0 keV energy range. The annular sectors of spectra extraction are indicated in the right panel of Figure 5. The normalizations (EMs) are in XSPEC units of 10−14nenpV/4π[DA(1 + z)]2. Errors are at the 90% confidence levels on a single parameter of interest. The last column shows the result of the F-test. See Section 5 for details.
Download table as: ASCIITypeset image
By contrast, we do not find clear spectral signatures of cool gas in the sectors NNW, NNE, and ENE, as expected from a visual inspection of the hardness ratio map. However, the lack of spectral evidence for multiphase gas could be due to the limitations of our data. Indeed, due to the relatively limited spectral resolution of Chandra, the detection of two different thermal components demands a significant temperature separation. The temperature difference required to have a marked effect on a single-phase thermal fit at the 90% confidence level is determined as follows. Using the response matrices, background, and numerical information of a real spectrum extracted in an arbitrary sector, we simulated spectra with two thermal components separated by ΔT around the best-fit value of the single-temperature fit to the real spectrum. We then fitted a single-temperature apec model to the mock two-temperature spectrum. We repeated this exercise, increasing the separation between the two temperature components until χ2ν exceeds the 90% confidence range. For comparison, we performed this procedure by starting from the real spectra in different sectors indicated in Figure 5 and found consistent results. An example of such an analysis is plotted in the left panel of Figure 7, where the dotted and dashed lines are the 90% and 3σ limits. We found that an ∼2.8 keV (3.0 keV) separation is necessary to exclude the presence of single-phase gas at the 90% (3σ) confidence for our data. The fact that the two spectral components detected in sector SSE are separated by ∼3.3 keV is consistent with this result. We note that kT1 found by the 2T model fit in the SSE sector is the lowest among all of the sectors, falling in the energy range which comprises most of the counts and thus is more easily detectable.
Download figure:
Standard image High-resolution imageDespite this, since a temperature higher than 4 keV is not observed at any radius in the cluster (Figure 4), the second thermal component found by the 2T model fit appears unrealistically hot, so we performed a new spectral fit with the 2T model (apec+apec) keeping the temperature of the second thermal component fixed at 4 keV. Such a model has only one additional free parameter than the 1T model: the normalization, EM2, of the second component (the metallicities of the two components are linked). The best-fitting parameter values and 90% confidence levels derived from the fits to the annular spectra in sectors are summarized in Table 4. The F statistics for the improvement over the 1T model, shown in the last column, indicate where the addition of a second thermal component is most significant (sectors SSE, NNW, NNE, and ENE). In agreement with the hardness ratio map (Figure 5), our spectral analysis therefore supports the presence of multiphase gas along the filaments. Interestingly, such cool filaments follow the morphology of the powerful central radio source nicely, although the western part of the southern radio lobe appears devoid of cool gas. The SSW sector is indeed the hottest and is also the region where the presence of a second thermal component is least significant. This sector lies at the location where the southern radio lobe appears to fold back on itself (Lane et al. 2004). The properties of the cavity are also consistent with a sharp bend in the southern jet there (Wise et al. 2007).
Table 4. Spectral Analysis of the Cool Filaments
2T Model with kT2 = 4 keV | |||||
---|---|---|---|---|---|
Sector | kT1 (keV) | EM1(× 10−3) | EM2(× 10−3) | χ2/dof | F Statistics |
WNW | 1.33+0.37 − 0.24 | 0.05+0.05 − 0.04 | 1.55+0.03 − 0.04 | 493.8/520 | 11.9 |
NNW | 1.67+0.16 − 0.11 | 0.39+0.08 − 0.07 | 1.42+0.06 − 0.07 | 478.8/524 | 50.5 |
NNE | 1.73+0.28 − 0.08 | 0.33+0.10 − 0.05 | 0.93+0.04 − 0.10 | 456.9/459 | 26.2 |
ENE | 1.72+0.35 − 0.13 | 0.22+0.09 − 0.05 | 0.95+0.05 − 0.09 | 436.6/454 | 19.7 |
ESE | 2.00+0.52 − 0.35 | 0.27+0.17 − 0.09 | 1.11+0.09 − 0.17 | 492.5/482 | 9.5 |
SSE | 1.33+0.12 − 0.08 | 0.18+0.05 − 0.04 | 1.26+0.03 − 0.03 | 539.1/490 | 75.2 |
SSW | 1.09+0.48 − 0.26 | 0.02+0.03 − 0.01 | 1.17+0.03 − 0.02 | 451.4/471 | 6.9 |
WSW | 1.37+0.63 − 0.35 | 0.04+0.05 − 0.03 | 1.34+0.03 − 0.04 | 433.1/493 | 7.2 |
Notes. Results of 2T model spectral fitting in the 0.5–8.0 keV energy range. The temperature value of the second thermal component is fixed to 4 keV. The annular sectors of spectra extraction are indicated in the right panel of Figure 5. The normalizations (EMs) are in XSPEC units of 10−14nenpV/4π[DA(1 + z)]2. Errors are at the 90% confidence levels on a single parameter of interest. The last column shows the F statistic value testing the improvement of the 2T model over the 1T model in Table 3. See Section 5 for details.
Download table as: ASCIITypeset image
Finally, we also attempted to map the EM distribution by fitting more complicated multiphase spectral models (such as apec+apec+apec+apec, with fixed temperatures and free normalizations). However, we found our data are inadequate due to limited statistics and the spectral resolution of Chandra. These detailed spectral studies will, hopefully, be possible in the future with the spectral capabilities of the International X-ray Observatory (IXO).
6. DISCUSSION
6.1. Scaled Temperature Profile and Shock Front
We show in the left panel of Figure 8 the observed temperature profile of Hydra A (red triangles), scaled by the virial radius (estimated from the relation ; Evrard et al. 1996), overlaid on the scaled temperature profiles of a sample of 12 relaxed clusters observed with Chandra (Vikhlinin et al. 2005). The temperatures are scaled to the emission-weighted cluster temperature 〈TX〉, measured after excluding the central 70 kpc region which is usually affected by radiative cooling. We note that the plateau of cool gas, already discussed in Section 4, stands out from to the typical temperature profile of clusters. In particular, lying just inside the typical temperature peak, such a plateau emphasizes the temperature peak in Hydra A, making it look higher than it actually is. Since this region is also where the shock front is located, it makes it difficult to distinguish between a temperature jump due to the shock and the typical temperature peak of the cluster.
Download figure:
Standard image High-resolution imageIn fact, once we overlay the scaled temperature profile of Hydra A measured after masking the cool filaments (right panel of Figure 4) it is remarkable how well it agrees with the general shape of the temperature profiles observed for relaxed clusters (see the right panel of Figure 8). Therefore, it is unlikely that the temperature peak observed in Hydra A is produced by the shock. In particular, the spectroscopic detection of a temperature rise in the regions immediately inside of the shock front is complicated by the underlying rising temperature profile of the global cluster atmosphere. In Figure 9, we show some examples of such temperature measurements in various sectors (see Section 3.2): the post-shock regions are found to be hotter than the corresponding pre-shock ones, although due to the large error bars, the pre-shock and post-shock temperatures are consistent.
Download figure:
Standard image High-resolution image6.2. Evidence for Gas Dredge-up along the Cool Filaments
The fact that the scaled temperature profile of Hydra A measured after masking the filaments agrees with the general shape of the temperature profiles observed for relaxed clusters (right panel of Figure 8) is a clear indication that these filaments are responsible for the plateau of cool gas. In Section 5, we performed a detailed (projected) spectral analysis and found evidence for multiphase gas, which may have been uplifted with the radio lobes, giving rise to the cooler filaments and plateau. In order to estimate the mass of the cool gas in the filaments, we must correct for the effects of projection along the line of sight. Starting from the annuli used to derive the global temperature profile, we binned them together in order to obtain three shells suitable for a simple deprojection analysis of the cool filaments. In particular, using three annuli: inner (0''–72'' ∼ 0–76 kpc), middle (72''–145'' ∼ 76–152 kpc), and outer (145''–339'' ∼ 152–360 kpc), we performed a deprojection analysis of the absorbed 2T model with the XSPEC projct model. Under the assumption of spherical shells of emission, this model calculates the geometric weighting factors, according to which the emission is redistributed amongst the projected annuli. The outer shell was assumed to include only gas at the ambient temperature, which is linked to the temperature values of the hot components in the inner and middle shells. The detailed results of the deprojection analysis are shown in Table 5.
Table 5. Deprojection Analysis
Shell | kTcool | EMcool | kThot | EMhot | Vtot | fcool | ne, cool | Mcool | fhot | ne, hot | Mhot |
---|---|---|---|---|---|---|---|---|---|---|---|
('') | (keV) | (× 10−3) | (keV) | (× 10−3) | (cm3) | (cm−3) | (M☉) | (cm−3) | (M☉) | ||
Inner | 1.77+0.13 − 0.06 | 3.07 | 4.01+0.04 − 0.04 | 13.3 | 5.35 × 1070 | 0.043 | 0.032 | 6.86 × 1010 | 0.957 | 0.014 | 6.68 × 1011 |
(0–72) | |||||||||||
Middle | 1.58+0.08 − 0.10 | 1.64 | 4.01+0.04 − 0.04 | 10.9 | 3.85 × 1071 | 0.023 | 0.012 | 9.90 × 1010 | 0.977 | 0.0047 | 1.65 × 1012 |
(72–145) | |||||||||||
Outer | ... | 0 | 4.01+0.04 − 0.04 | 12.5 | 5.20 × 1072 | 0 | 0 | 0 | 1 | 0.0136 | 6.59 × 1012 |
(145–339) |
Notes. Results of the deprojection analysis in three concentric 360° annular regions (indicated in the first column) in the 0.5–8.0 keV energy range using the XSPEC projct×wabs×(apec+apec) model. The normalization of the cool component in the outer shell is fixed to zero, and the temperature values of the hot components are linked between all shells. The fit gives χ2/dof = 4993.5/2811. The normalizations, EMcool(hot), are in XSPEC units of 10−14ne, cool(hot)npVcool(hot)/4π[DA(1 + z)]2. Errors are at the 90% confidence levels on a single parameter of interest. The ratio of the relative volumes occupied is estimated as Vcool/Vhot = (EMcool/EMhot) · (kTcool/kThot)2, and the filling factor f is defined as fcool(hot) = Vcool(hot)/Vtot, with Vtot = Vcool + Vhot. The masses are thus estimated from the relation: M = ρfVtot, where we use the conversion from electron number density to gas density ρ = 1.83 μnemH.
Download table as: ASCIITypeset image
We found that the hot phase has a temperature of 4 keV, in agreement with the observed global temperature profile (Figure 4), and that the cool component is, on average, at ∼1.6 keV. By taking into account the relative filling factor, f, calculated under the assumption that the two thermal components are in pressure equilibrium in each shell, in the middle shell we estimated a mass of cool gas of 9.9 × 1010M☉ (and a mass of hot gas of 1.7 × 1012M☉). We note that this estimate is considerably larger than the mass of cool gas reported by Simionescu et al. (2009b). The discrepancy between the two measurements can be ascribed to the fact that we included the mass content of the mushroom-cap (namely, sectors ENE, NNE, NNW in Figure 5). Using the entropy of the cooler phase and comparing to the entropy profile of the Hydra A cluster measured by David et al. (2001), we could estimate where it originated. From the density and temperature values measured in the middle shell (see Table 5), the cool gas has an entropy S = kTcool n−2/3e, cool ∼ 30 keV cm2. As shown by the profile of David et al. (2001), gas at such entropy is located around 30 kpc from the cluster center. Therefore, it is very likely that the cool gas, which is now observed in the filaments at ∼75–150 kpc, was lifted from the central 30 kpc. In this case, the most obvious way to move the gas is by some form of entrainment or dredge up of cool material from the center of the cluster in the rising lobes. Churazov et al. (2000) have indeed shown that hot buoyant bubbles produced by the expanding central radio source can dredge up cool material from the center of the cluster. The mass of cool gas at ∼75–150 kpc is 9.9 × 1010M☉, which is ∼60% of the total mass of gas remaining within 30 kpc (1.7 × 1011M☉; David et al. 2001). Assuming that such a mass was lifted out of the central 30 kpc by a continuous outflow or a series of bursts from the nucleus of the central galaxy over the past 200–500 Myr (which created the cavity system; Wise et al. 2007), it would amount to outflows of a few hundred M☉ yr−1. There would thus be a development of a convectively unstable region that can significantly reduce the net inflow of cooling gas (see also David et al. 2001; Nulsen et al. 2002). Therefore, our results show that the AGN feedback in Hydra A is acting not only by directly re-heating the gas, but also by removing a substantial amount of potential fuel for the supermassive black hole.
The energy required to lift the cool gas gives a lower limit on the amount of AGN outburst energy deposited in the ICM. This value can be estimated by calculating the variation in gravitational potential energy during the lifting process. If we assume that the undisturbed ICM is approximately isothermal with sound speed cs ≈ 1000 km s−1 and is in a hydrostatic configuration with density profile ρ(r), we can calculate this quantity as (Reynolds et al. 2008)
where Mcool is the lifted mass, ρi and ρf are the initial and final densities of the surrounding ICM, and γ = 5/3 is the ratio of specific heat capacities. From the density profile presented by David et al. (2001) we estimated that the energy required to lift the cool gas is ≳ 2.2 × 1060 erg. This value is comparable to the work required to inflate all of the cavities against the surrounding pressure (Wise et al. 2007) and is ∼25% of the total energy of the large-scale shock (Nulsen et al. 2005). Although we find evidence for a more extended gas dredge-up than previously estimated, there is a remarkable correlation between the cool filaments studied here and the metal-rich filaments in the iron-abundance maps measured by Simionescu et al. (2009b) and Kirkpatrick et al. (2009). This is consistent with the idea that Hydra A's powerful radio source is able to lift cool, metal-rich gas from the central region and distribute it throughout the X-ray atmosphere of the cluster. A similar effect is observed in M87 (Simionescu et al. 2008; Werner et al. 2010). We finally note that by summing our estimates of the integrated mass of cool gas, Mcool(< 152 kpc) ∼ 1.7 × 1011M☉, and hot gas, Mhot(< 152 kpc) ∼ 2.3 × 1012M☉, in the inner and middle shells we measure a total mass of gas Mgas(< 152 kpc) ∼ 2.5 × 1012M☉ which is in agreement with the gas mass profile measured by David et al. (2001), so our general picture is self-consistent.
7. SUMMARY
The main results of this work can be summarized as follows.
- 1.We studied the azimuthal properties of the weak (Mach number ∼1.2–1.3), large-scale shock and attempted to measure the temperature jump associated with the shock in different directions. The errors in the temperature measurements are too large to constrain the temperature jump caused by the shock. Furthermore, we note that the detection of a temperature rise in the regions immediately inside of the front is complicated by the underlying rising temperature profile of the cluster atmosphere.
- 2.Our detailed spectral analysis indicates the presence of multiphase gas along soft filaments seen in the hardness ratio map. The cooler gas has a significant impact on the radial temperature profile of the cluster. After correcting for the effect of the cool filaments, Hydra A's temperature profile is consistent with the form observed in relaxed galaxy clusters. Thus, it is unlikely that the observed temperature peak is produced by the shock.
- 3.The cool filaments are direct evidence for substantial dredge-up of low-entropy material by the rising lobes: ∼60% of the gas mass remaining in the central 30 kpc has been lifted along the cool filaments to the observed current position of 75–150 kpc. The outflow amounts to a few hundred M☉ yr−1, which is comparable to the inflow rate formerly estimated for the cooling flow. The energy required to lift the cool gas is ∼25% of the total energy of the large-scale shock and is comparable to the work required to inflate the cavities.
We thank the referee for constructive comments. We also thank A. Vikhlinin for providing the data used to make the plots in Figure 8. M.G. thanks D. Rafferty for precious advice in data reduction, and F. Brighenti for many insightful discussions. M.G. acknowledges support by grants ASI-INAF I/023/05/0 and I/088/06/0 and by Chandra grants GO0-11003X and GO0-11136X. B.R.M acknowledges support from the Natural Sciences and Engineering Research Council of Canada.