A&A 484, 317-325 (2008)
DOI: 10.1051/0004-6361:200809493

Cold gas in the Perseus cluster core: excitation of molecular gas in filaments

P. Salomé1,[*] - F. Combes2 - Y. Revaz2,6 - A. C. Edge3 - N. A. Hatch5 - A. C. Fabian4 - R. M. Johnstone4


1 - Institut of Radio Astronomy (IRAM), Domaine Universitaire, 300 rue de la piscine, 38400 St Martin d'Hères, France
2 - Observatoire de Paris, LERMA, 61 Av. de l'Observatoire, 75014 Paris, France
3 - Department of Physics, University of Durham, South Road, Durham DH1 3LEi, UK
4 - Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK
5 - Leiden Observatory, 2300 RA Leiden, The Netherlands
6 - Laboratoire d'Astrophysique, Ecole Polytechnique Fédérale de Lausanne (EPFL), Observatoire, 1290 Sauverny, Switzerland

Received 1 February 2008 / Accepted 13 March 2008

Abstract
We have recently detected CO lines in the well-known filaments around NGC 1275, the galaxy at the centre of the Perseus cluster of galaxies. These previous observations, with the HERA multi-beam array at the IRAM 30 m telescope enabled us to make a large map of the CO(2-1) line and to see hints of molecular gas far away from the cluster centre. To confirm the presence of CO emission lines in the outer filaments and to study the CO(2-1)/CO(1-0) line ratio, we observed seven regions of interest again with the 30 m telescope in both CO(1-0) and CO(2-1). The regions we observed were: the eastern filament, the horseshoe, the northern filament and a southern extension, all selected from H$\alpha $ emission line mapping. Molecular gas is detected in all the observed regions. This result confirms the large extent of the cold molecular gas filaments. We discuss the CO(2-1)/CO(1-0) ratios in the filaments. The eastern filament has optically thick gas, whereas further away, the line ratio increases close to values expected for a warmer optically thin medium. We also show CO(1-0) and CO(2-1) lines in 9 regions closer to the centre. The kinematics of the CO is studied here in more detail and confirms that it follows the motions of the warm H2 gas found in the near-infrared. Finally, we searched for dense gas tracers around 3C 84 and claim here the first detection of HCN(3-2).

Key words: galaxies: elliptical and lenticular, cD - galaxies: cooling flows - galaxies: intergalactic medium - galaxies: individual: NGC 1275

1 Introduction

Over the past decade X-ray observations have shown a lack of cool X-ray emitting gas in cooling flow clusters (e.g. Peterson et al. 2003), and many feedback models have been proposed to explain the required energy injection into the intracluster medium (e.g. Binney & Tabor 1995; Omma & Binney 2004). At the same time, large amounts of molecular gas have been found in many Brightest Cluster Galaxies (BCG) (Edge 2001; Salomé & Combes 2003), providing evidence that some of the ICM may cool to very low temperatures. Detailed studies of this cold molecular gas reservoir are a complementary way to probe feedback processes.

The giant cD galaxy NGC 1275 is the central galaxy of the Perseus cluster (Abell 426) and lies at a redshift of 0.01756. At this distance, 1'' is 350 pc (H0 = 71 km s-1 Mpc-1, $\Omega_M=0.27$, $\Omega_\lambda=0.73$). This object is famous for the huge filamentary structure detected in the optical (Hu et al. 1983; Conselice et al. 2001). These bright H$\alpha $ emitting filaments are known to be directly/indirectly associated with cooling cores (Crawford et al. 1999; Edwards et al. 2007). We have shown previously (Salomé et al. 2006, hereafter S06) that the molecular gas is also detected in CO(2-1) emission in the centre of NGC 1275 with morphology and dynamics identical to that of the H$\alpha $ emission.

We present here the results of follow-up observations aimed at confirming the presence of cold gas associated with the optical filaments further out from the galaxy. We looked for CO(1-0) and CO(2-1) emission lines in 7 regions selected in the optical filaments. We also observed the central region at both 1.3 and 3 mm with high sensitivity receivers. Section 2 describes the observations. Section 3 presents the results and discusses the spectra obtained in the different regions of interest. We then compare the CO results (morphology, kinematics, line ratios) with data at other wavelengths in Sect. 4. Section 5 summarizes our conclusions.

2 Observations

The observations were made with the IRAM 30 m telescope on Pico Veleta, Spain during two different runs in 2006. We used the wobbler switching mode with two 3 mm and two 1.3 mm receivers operating simultaneously. The beam throw was 4 arcmin. Receivers were tuned to the CO(1-0) and CO(2-1) emission lines, redshifted to the velocity of NGC 1275. Frequent pointings were done on 3C 84, the radio source at the centre of NGC 1275. At 3 mm, we used two $512 \times 1$ MHz filter-banks. This gives a total band of $\sim$1300 km s-1 for the CO(1-0) line setup. In addition, we used the two $250 \times 4$ MHz resolution filter-banks for the 1.3 mm receivers, providing a 1 GHz band width, for the 1.3 mm receivers which also corresponds to a 1300 km s-1 bandwidth. The beams of the 30 m telescope, at 3 mm (113.282 GHz) and 1.3 mm (226.559 GHz), are respectively 22 $^{\prime\prime}$ and 11 $^{\prime\prime}$.

Regions close to the centre were very quickly detected (in 10-20 min), while the more distant regions required more time ($\sim$2 h). The signals are expressed here in main beam brightness temperatures. The main-beam efficiency of the 30 m telescope is the ratio of the antenna temperature to the main beam temperature: $T_{\rm A}^*/T_{\rm mb} =
B_{\rm eff}/F_{\rm eff}$ with the ratio of the beam efficiency to the forward efficiency being: $B_{\rm eff}/F_{\rm eff}= 0.75/0.95$ at 3 mm and 0.52/0.91 at 1.3 mm (cf IRAM-30 m site http://www.iram.es/).

The data were calibrated with the MIRA software and reduced with the CLASS90 package. Spiky channels and bad scans were dropped and linear baselines were subtracted for each spectrum. After averaging all the spectra for each line at each position, the data were Hanning smoothed to a 42 km s-1 resolution.

The data are summarized in Tables 1 and 2. The CO(2-1) lines of the filaments use the present observations and the data obtained with HERA (S06): the CO(2-1) spectra were convolved with the beam pattern at 3 mm to compare the CO(1-0) and CO(2-1) temperatures.

Table 1: Results of the observations.

Table 2: Results of the observations.

3 Results

We detected CO in all the regions observed inside the H$\alpha $filaments. Offsets are relative to the 3C 84 position: RA 03:19:48.15, Dec 41:30:42.1 (J2000). In the central region (indicated as centre 1 to 9 in Fig. 1 and Table 2), the 3 mm data suffered from baselines instabilities due to the strong continuum of 3C 84. The 3 mm spectra shown on the left hand side of Fig. 4 are affected by those baseline ripples. The 1.3 mm data, not affected by baseline ripples give a more reliable idea of the line shape. Nevertheless CO lines are detected in all regions except the region outside the H$\alpha $ filaments. This is the region named test, which has no H$\alpha $ emitting counterpart. Our earlier HERA map (S06) showed a hint of a detection at this place, but this was not confirmed by the new observations.

In order to compute the molecular gas content from the integrated CO intensity $I_{{\rm CO}}$ (K km s-1), we used

\begin{displaymath}{I}_{\rm CO}= \int T_{\rm mb}{\rm (CO)} {\rm d}V,
\end{displaymath} (1)

where $T_{\rm mb}{\rm (CO)}$ is the main beam antenna brightness temperature, obtained with the CO emission line and dV is the line width. We then adopted the standard Milky Way conversion factor commonly used for N(H2) estimates:

\begin{displaymath}N({\rm H_2)} = 2.3 \times 10^{20} { I}_{\rm CO} {\rm (molecules/cm^2)}.
\end{displaymath} (2)

From this equation the mass of molecular hydrogen, contained in one beam is:

\begin{displaymath}M_{\rm gas} (M_{\odot}) = 1.36 \times 2.95 \times 10^{-19} { I}_{\rm CO}~
\theta^2 D^2 \frac{N({\rm H}_{2})}{I_{\rm CO}},
\end{displaymath} (3)

$I_{{\rm CO}}$ is the integrated intensity in K  km s-1, $\theta$ is the beamsize of the telescope in arcsec, and D is the distance of the galaxy, taken to be 72.6 Mpc. We also included a factor 1.36 to take into account the He contribution to the gas mass.


  \begin{figure}
\par\includegraphics[origin=rb,width=12cm,angle=-90]{9493fi1.ps}\end{figure} Figure 1: H$\alpha $ image of the filament system in NGC 1275 (Conselice et al. 2001). Overlaid are the regions re-observed with the 30 m telescope. The circles represent the beam size at 3 mm.
Open with DEXTER


  \begin{figure}
\par\begin{tabular}{ccc}
\includegraphics[width=4cm,angle=-90]{9...
...includegraphics[width=4cm,angle=-90]{9493fi19.ps}\\
\end{tabular}
\end{figure} Figure 2: CO(1-0) and CO(2-1) spectra obtained at all the positions observed as indicated at lower right in each diagram. The channel width is 42 km s-1. On the left hand side are the CO(1-0) lines detected with the a100 and b100 receivers. In the middle are the results obtained for the CO(2-1) line with the A230 and B230 receivers. On the right hand side are the CO(2-1) lines computed with both A230 and B230 merged with previous HERA data and smoothed to the 3 mm beam size.
Open with DEXTER

3.1 The filaments

We describe below the regions selected from our HERA map (see S06). These regions were chosen to be distant from the East-West filament seen in CO(2-1) and known for their peculiar optical morphology. The first aim was to confirm the CO detections far away from the central galaxy. The regions called Off 1, Off 2 and Pos 2 trace the northern and southern filaments visible in H$\alpha $ up to a projected distance of 25 kpc north and south of the galaxy. Pos 11 is centred on a looped-back filament seen in optical emission lines that may trace uplifted gas behind a rising cold gas bubble inside the hot intracluster medium (Fabian et al. 2003; Hatch et al. 2006). Finally, closer to the centre, East and Off 3, are two fields in the eastern filament. The East region includes a young star cluster (Shields & Filipenko 1990).

We now compare the velocity shift of our CO detections with the velocities of the warm (2000 K) H2 gas detected by Hatch et al. (2005) and the hot (104 K) optical line emitting gas measured by Hatch et al. (2006).

The eastern filament: refers to positions East and Off 3, at about 8 kpc from the galaxy's centre. It corresponds to regions named A, B, C and D in Hatch et al. (2005). The CO emission lines are centred at -100 km s-1 relative to the systemic velocity. This blueshift is the same as the velocity found for the P$\alpha $, Br$\gamma$ and H2 v=1-0 S(1), S(2)and S(3) lines detected by the above mentioned authors.

The southern filament: refers to position Pos 2 at 25 kpc from the galaxy's centre. It corresponds to regions named SW1 and SW2 in Hatch et al. (2005). Here again, the CO emission is blueshifted relative to the nucleus, like the infrared H2 tracers.

The horseshoe: refers to the position Pos 11, at least 25 kpc from the galaxy's centre. The CO(1-0) line is detected at the limit of the sensitivity reached. There is also a tentative detection of CO(2-1) line which disappears when merged with the HERA data, so this detection is not very strong. Nevertheless, the two emission lines have positive velocities (36 km s-1 and 120 km s-1). This velocity range agrees with that from the infrared H2 tracers by Hatch et al. (2005). It is also the same velocity shift as the hot H$\alpha $ emitting gas found by Hatch et al. (2006) in the loop of the horseshoe filament (see their Fig. 6).

The northern filament: This region refers to the position Off 2, at about 25 kpc from the galaxy centre. The H$\alpha $ emitting gas has a velocity range between +6 and +41 km s-1 relative to the systematic velocity. This is similar to the velocity of the CO line.


  \begin{figure}
\par\begin{tabular}{ccc}
\includegraphics[width=4cm,angle=-90]{9...
...includegraphics[width=4cm,angle=-90]{9493fi34.ps}\\
\end{tabular}
\end{figure} Figure 3: CO(1-0) and CO(2-1) spectra obtained at all the positions near the centre of the galaxy as labeled in the lower right of each diagram. The channel width is 42 km s-1, see Table 2.
Open with DEXTER


  \begin{figure}
\begin{tabular}{ccc}
\includegraphics[width=4cm,angle=-90]{9493f...
...ludegraphics[width=4cm,angle=-90]{9493fi48.ps} & \\
\end{tabular}
\end{figure} Figure 4: CO(1-0) and CO(2-1) spectra obtained at all the positions near the centre of the galaxy as labeled in the lower right of each diagram. The channel width is 42 km s-1, see Table 2.
Open with DEXTER

The tangential filament: This region refers to position Off 1, at least 12.3 kpc from the galaxy centre. It corresponds to a region inside the tangential filament (slit 5) of Hatch et al. (2006). The CO velocity is not well defined here as it is not consistent between the 3 mm and the 1.3 mm data. However, there is a hint of a two-component velocity structure, at $\pm$100 km s-1. This is consistent with the H$\alpha $ velocity found in the same region (distance between $\sim$15 and $\sim$20 kpc in Fig. 7 of Hatch et al. 2006).

The test region: A hint of a CO(2-1) detection was present in the S06 HERA data, in a region that corresponds to no H$\alpha $filament. We re-observed this region and found no CO(1-0) or CO(2-1) emission (i.e. $M_{\rm gas} <10 ^8~M_\odot$). We conclude that there is no, or very little, cold gas outside the H$\alpha $ filaments.

3.2 The central region

We observed 9 regions to map the CO emission close to 3C 84. We present here a fully sampled CO(1-0) map of about $50''\times 50''$ with a spatial resolution of 22''. The CO(2-1) map does not fully sample the same region but has a better resolution of 11''. All the results are described in more detail in Table 2 and Figs. 3 and 4. In the central regions the CO(1-0) and CO(2-1) line profiles are similar. The CO(2-1) line fluxes, when not convolved with the CO(1-0) beam, are 4 times larger than the CO(1-0) fluxes. So within the uncertainties, this is in agreement with an optically thick cold gas. After convolution, all the ratio are close to one. Some regions have a ratio greater than one, but as mentioned before, the CO(1-0) lines are affected by baseline ripples close to 3C 84. The main beam brightness temperatures in CO(1-0) are thus very likely under-estimated in this region.

3.3 Detection of HCN(3-2)

We searched for HCN(1-0), HCN(2-1) and HCN(3-2) emission lines in the centre of NGC 1275. The observations at low frequency were corrupted by strong baseline ripples due to the strong continuum source 3C 84. So we focus on the HCN(3-2) observations at 261 GHz, observed in average weather conditions. The spectrum presented in Fig 5 is the average of data taken during a first run with the two receivers detecting two polarisations (A230, B230) and during a second run with the two other receivers C270 and D270. The signal appears marginally in both data sets and even better in the averaged spectrum. So we claim a detection of HCN(3-2). The S/N ratio on the integrated line is 6.

The HCN molecule has a large dipole moment that requires n(H 2) > 104 cm-3 for significant excitation. This tracer of dense gas regions is thus considered as an indicator of star formation activity (Gao & Solomon, 2004). We computed the apparent CO and HCN luminosities (using Eq. (1) from Downes et al. 1999) from the integrated line flux measured in CO(1-0) and HCN(3-2). We assumed a HCN(3-2)/HCN(1-0) ratio close to one. At densities higher than a few 104 cm-3 for a temperature higher that 40-50 K, radiative transfer models predict a ratio larger than $\sim$0.25. We found $L'_{\rm CO}=5\times 10 ^8$ K km s-1 pc2 and $L'_{\rm
HCN} = 5 \times 10 ^7$ K km s-1 pc2. When compared to the $L_{\rm FIR}$- $L'_{\rm CO}$ and $L_{\rm FIR}$- $L'_{\rm HCN}$ relations as described by Solomon et al. (1992, 1997), our measurements fall in the normal spiral galaxies range and independently predict an $L_{\rm FIR}$of $2\times 10 ^{10}~L_\odot$. Lester et al. (1995) estimated the thermal FIR emission to be $10^{11}~L_\odot$, with the approximation that only 20$\%$ of the total luminosity was non-thermal. This is an order of magnitude higher than the $L_{\rm FIR}$ predicted here. So either the contribution of the synchrotron emission is larger than expected or there is an excess of FIR emission in NGC 1275 compared to the well-known $L_{\rm FIR}$-L $'_{\rm HCN}$ relation for star forming galaxies.

3.4 Molecular gas from the high velocity component

NGC 1275 is surrounded by some star formation regions, with young star clusters, that could be formed in a merger but also in the cooling flow (see S06). Some are associated with the High Velocity System (HVS) at 8200 km s-1, which appears in absorption in front of the optical emission of NGC 1275 (Gillmon et al. 2004). No stellar component is detected, even in the near-infrared, corresponding to the HVS, which could be a gas-rich disrupted dwarf, or tidal debris, as proposed by Hu et al. (1983).

We searched for CO emission at 8280 km s-1 on one position in the HVS (centred on 3C 84's position). The results are shown in Fig. 6. There is no clear detection but a hint of CO(1-0) emission. The line intensity we derive is up to a maximum of the order of 0.4 K km s-1 which corresponds to $\le$10 $^8~M_\odot$ of cold gas. So even if interacting with the ambient medium closely surrounding NGC 1275, the HVS contains a very small amount of cold gas. This mass is comparable or less than what is found in the cold filaments around the cD galaxy. Therefore, if the HVS merges with the cD galaxy, it will not deposit a significant amount of cool gas.

4 Discussion

4.1 CO(2-1)/CO(1-0) line ratios

We have computed the CO(2-1)/CO(1-0) line ratios. The data at 1.3 mm were combined with the HERA data described in S06. We then convolved the CO(2-1) data by the CO(1-0) beam size in order to make a meaningful comparison. We list, in Table 1, the CO(2-1)/CO(1-0) line ratios of the different regions observed. If we compare the main beam temperature in the CO(2-1) spectra with and without the 3 mm beam shape convolution (central and right hand side column of Fig. 2), we can see that we have lost some signal by diluting the 1.3 mm emission with a larger beam size. This means that the molecular gas emission is very likely to come from a region that is smaller than the 1.3 mm beam size (11 $^{\prime\prime}$).

At Pos 11, the tentative detection of CO(2-1) disappears when convolved to the CO(1-0) beam. Regions Off 2, Pos 2, Pos 11 and East show lower main beam temperatures at 1.3 mm when we applied the 3 mm beam convolution. Therefore, the emission from these positions probably comes from regions with an angular size less than 22 $^{\prime\prime}$ (the beam size at 3 mm). Regions Off 1 and Off 3 do not show lower main beam temperatures when convolved with the 3 mm beam and are therefore probably extended by at least 22 $^{\prime\prime}$. We have recently confirmed this last result by mapping the eastern filament with the Plateau de Bure interferometer. This shows that the molecular gas lies in thin ($\le$2'') and elongated structures, exactly coincident with the H$\alpha $ filaments (Salomé et al. 2008).

After convolution to the same beam size, we find the CO(2-1)/CO(1-0) line brightness temperature ratio is close to one in regions East, Off3, Pos2 and close to 3C 84. Such a ratio indicates optically thick CO(1-0) and CO(2-1) emission lines (Eckart et al. 1988). For regions farther out from the galaxy centre like Off2, the ratio is slightly larger than one. However, we sometimes find the line widths smaller in CO(1-0) than in CO(2-1). This means that the apparently high ratio could actually be closer to one if the line widths were constrained to be the same. Note that the Off1 emission-line ratio is about 2.5 while the line shapes are the same at 1.3 mm and 3 mm. So we see hints of higher CO(2-1)/CO(1-0) ratios at larger radius from the galaxy centre. This particular high ratio (larger than one) indicates optically thin emission (Goldsmith et al. 1999).

Closer to the centre, the strong 3C 84 source produced baseline ripples that affected the 3 mm line shapes. However, we estimated CO(2-1)/CO(1-0) line ratios close to one, in agreement with an optically thick cold gas (see Table 2).

4.2 Molecular gas reservoir and star formation

The scenario for the origin of the molecular gas we detect is that the cold gas accumulating towards the centre of the cluster gets dragged out with the bubbles and compressed in curved shocks, or flows back down around the borders of cavities. Feedback from the AGN not only heats the gas (negative feedback), but also compresses it in thin layers, that favors cooling and formation of molecular filaments (positive feedback) at a distance greater than 10 kpc (Revaz et al. 2008).

We expect no significant amounts of dust in the filaments, due to high sputtering rates. The question is then whether, in these conditions, there is still sufficient dust (which acts as a catalyst) to form molecules quickly. It may be that some dust from the ISM of NGC 1275 has been captured inside the uplifted cooler gas and this is sufficient to catalyse the formation of molecules. Dust could also form in interstellar shocks, if present. The presence of dust grains at large radii would then help to form molecules in the outer filaments. A small quantity of dust grains is sufficient. We have proposed further observational investigations in order to constrain the amount of dust inside the filaments of NGC 1275.


  \begin{figure}
\par\includegraphics[width=5cm,angle=-90]{9493fi49.ps}\end{figure} Figure 5: HCN(3-2) emission lines in T$_{\rm mb}$ (mK), observed with the 30 m telescope (290 min). A gaussian fit gives I $_{\rm HCN}$ 2.1 K km s-1 with a peak temperature of 5 mK, a velocity of 120 km s-1 and a linewidth of 400 km s-1. The rms is 1.6 mK with a mean T $_{\rm sys}$ of 445 K. The channel width is 36 km s-1.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=5cm, angle=-90]{9493fi50.ps}\\
\includegraphics[width=5cm, angle=-90]{9493fi51.ps}\\
\end{figure} Figure 6: Observations at 112.175 GHz ( top) and 224.346 GHz ( bottom), the CO(1-0) and CO(2-1) frequencies expected in the HVS restframe at 8280 km s-1. No CO(2-1) is detected. There is a hint of detection in CO(1-0), with $I_{\rm CO} = 0.4 \pm 0.1$ K. km s-1. This would give a mass of about $1 \times 10 ^8~M_\odot$. The channel width is 21 km s-1.
Open with DEXTER

For the assumed standard CO/H2 conversion factor, the amount of molecular gas detected in the different regions observed varies from $7 \times 10 ^7~M_\odot$ to $8 \times 10 ^8~M_\odot$, which when added together gives a total mass of $1.3 \times 10 ^9~M_\odot$. These filaments therefore contain large quantities of cold molecular gas. The regions are also conspicuous in the emission lines of H$\alpha $ and [NII] but the optical line ratios are incompatible with young stars as the source of excitation as might be expected if stars were forming from the molecular gas. The excitation of these lines is not coming from the AGN either, since it does not decrease with distance from the nucleus of NGC 1275 (Hatch et al. 2006). The filaments are emitting extensively in UV and optical lines, much more than in X-rays, and a possible explanation for this is through heating by shocks at the interface between the relativistic jet plasma and the entrained gas at the exterior of the bubbles.

One of the most CO-rich regions is the eastern extension, which contains an unresolved (<1 kpc) young star cluster (Shields et al. 1990). The molecular gas mass of this region is greatly in excess of the mass of the star cluster and confirms a low star formation efficiency within the filaments (see also Salomé et al. 2008, for more details).

5 Conclusions

New CO(1-0) and CO(2-1) observations towards selected optical filaments around NGC 1275 confirm the presence of CO associated with the H$\alpha $ filaments. The gas kinematics are the same in CO and as in H$\alpha $. We detect for the first time, a large amount of cold gas inside the outer filaments (a few 10 $^8~M_\odot$).

The excitation of the gas, derived from the comparison of CO(2-1) and CO(1-0) convolved to the same spatial resolution, indicates a range of densities are present. Optically thick emission is present close to the galaxy centre, whilst further out the emission could be optically thin.

We report a detection of HCN(3-2) emission in the centre, indicating high density gas. The association of the molecular gas with the H$\alpha $filaments is consistent with the scenario in which cooling of the gas is taking place at large radii. This could happen in compressed hot gas at the border of the X-ray cavities. A detailed interpretation awaits specific simulations taking into account the cold, warm and hot gas phases.

Acknowledgements
IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). R.M.J. acknowledges support by STFC and the Royal Society. We thank D. Downes for helpful comments and discussions.

References

 

Copyright ESO 2008