-
PDF
- Split View
-
Views
-
Cite
Cite
A C Fabian, J S Sanders, G J Ferland, B R McNamara, C Pinto, S A Walker, Hidden cooling flows in clusters of galaxies – III. Accretion on to the central black hole, Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 1, September 2023, Pages 716–730, https://doi.org/10.1093/mnras/stad1870
- Share Icon Share
ABSTRACT
Recently, we have uncovered hidden cooling flows (HCFs) in the X-ray spectra of the central brightest galaxies of 11 clusters, 1 group, and 2 elliptical galaxies. Here, we report such flows in a further 15 objects, consisting of 8 clusters, 3 groups, 3 ellipticals, and 1 Red Nugget. The mass cooling rates are about |$1\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| in the ellipticals, 2 to |$20\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| in the groups, and 20 to |$100\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| in regular clusters. The Red Nugget, MRK 1216, has an HCF of |$10\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|. We review the fate of the cooled gas and investigate how some of it might accrete on to the central black hole. The gas is likely to be very cold and to have fragmented into low-mass stars and smaller objects before being swallowed whole, with little luminous output. If such a scenario is correct and operates at a few |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| then such objects may host the fastest growing black holes in the low-redshift Universe. We briefly discuss the relevance of HCF to the growth of early galaxies and black holes.
1 INTRODUCTION
We have recently found hidden cooling flows (HCFs) in clusters and groups of galaxies, as well as a couple of nearby elliptical galaxies (HCFI and HCFII) (Fabian et al. 2022, 2023), using spectra from the XMM Reflection Grating Spectrometer (RGS). These soft X-ray-emitting flows are hidden within photoelectrically absorbing cold clouds and dust near the centres of the central brightest galaxies. They represent the cooler inner parts of larger, wider-scale cooling flows. Active galactic nucleus (AGN) feedback acts to reduce the main cooling flow in the larger body of these objects but the inner parts drop from direct view behind cold absorbing clouds. The total mass cooling rates can be 20–50 per cent or more of the unabsorbed rates inferred earlier from X-ray imaging studies.
The findings again raise the ‘cooling flow problem’ of what happens to the cooled gas? HCF mass cooling flow rates of tens of solar masses per year in regular clusters and |$1\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| in early-type galaxies (ETGs) lasting |$\sim 8{\rm \, Gyr}$| (since redshift z = 1) means almost |$10^{11}$| and |$10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$|, respectively, of accumulated cooled gas. Where does it go? The issue is not new1 but has largely been ignored for the past two decades, even at the low rates allowed without absorption (see Liu et al. 2019 and Section 2).
We have proposed and discussed several possibilities, namely that (a) the gas cools to invisibility (i.e. so cold that it radiates little), (b) the cooled gas fragments into low-mass stars and substellar objects, and (c) cooled gas is dragged out from the centre by the bubbling action of AGN feedback. (a) and (b) mean that there is increasing unseen mass of gas and/or low-mass stars at the centres of these objects. (c) may be consistent with observed metal abundance profiles. These possibilities are not of course mutually exclusive.
Here, we investigate how much cooled gas can end up in the central black holes. Many of the most massive black holes at low redshift lie in brightest cluster galaxies (BCGs; McConnell & Ma 2013; Bogdán et al. 2018), and we include a couple here including Holm 15A, the central galaxy of A85, which has a black hole of mass 4 × 1010 (Mehrgan et al. 2019). There is some evidence that the black hole to galaxy stellar mass ratio of ETGs has increased significantly from z = 1 to the present day (Farrah et al. 2023). Since massive black holes can swallow stars whole, such accretion need not be luminous.
We now search for hidden cooling flows in eight cool core clusters, three X-ray luminous groups, and four relatively isolated elliptical galaxy including a Red Nugget. They are found in all objects and have the typical mass cooling rates found in HCFI and HCFII. One is the very X-ray luminous cluster, ZW3146, at medium redshift z ∼ 0.3, the results for which compare well with other high-luminosity clusters found at similar and higher redshift. A significant part of its high cooling rate of |$\sim 1000\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| goes into observed normal star formation but it is unlikely to be a long-lived situation. As discussed in HCFII (Fabian et al. 2023), rapid accretion on to the central black hole has the potential to turn it into a luminous quasar, as seen in the Phoenix cluster, perhaps ending in a massive outburst such as has occurred in MS 0735+7421 (McNamara et al. 2005).
We then speculate whether hidden accretion is taking place into the central black holes of HCFs. Hidden in the sense of unobserved because the infalling matter consists of low mass stars, etc., which are swallowed whole without emitting radiation. Finally, we speculate on high-pressure star formation, which occurs in HCF, and discuss its relevance to early galaxy formation and in particular to the origin of ‘Red Nuggets’.
2 SPECTRAL ANALYSIS
The objects and data used are listed in Table 1. SAS 20.0 was used for the data reduction. The spectra were extracted using rgsproc with a 95 per cent extraction in point spread function width (corresponding to 1.7 arcmin) and 95 per cent in pulse-height distribution. To create good time intervals (GTIs), light curves were created for each RGS instrument from CCD 9, with events with flag values of 8 or 16, extracted with a cross-dispersion angle of greater than 1.5 × 10−4, in time bins of 200 s. The GTIs were created when the rate was below 0.3 cts s−1. Background spectra were created with rgsbkgmodel. The spectra and background spectra from RGS1 and RGS2 were combined using rgscombine before spectral fitting. The spectra were then analysed using xspec (Arnaud 1996) over the energy range of 8–22A, where background is minimized.
Observed targets, giving the used source position (deg; J2000), observation identifiers, and average cleaned exposure of the RGS cameras.
Target . | RA . | Dec. . | OBSIDs . | Exposure (ks) . |
---|---|---|---|---|
2A0335 | 54.6691 | 9.9697 | 0 109 870 101 0109870201 0 147 800 201 | 145 |
A85 | 10.4601 | −9.3031 | 0 065 140 101 0723802101 0 723 802 201 | 215 |
A496 | 68.4074 | −13.2619 | 0 135 120 201 0506260301 0 506 260 401 | 162 |
A2597 | 351.3321 | −12.1243 | 0 108 460 201 0147330101 0 723 801 601 0723801701 | 257 |
A2199 | 247.1594 | 39.5512 | 0 008 030 201 0008030301 0 008 030 601 0723801101 0 723 801 201 | 137 |
M87 | 187.7059 | 12.3911 | 0 114 120 101 0200920101 0 803 670 501 0803670601 0 803 671 001 0803671101 | 430 |
NGC1399 | 54.6210 | −35.4505 | 0 012 830 101 0400620101 | 139 |
NGC720 | 28.2519 | −13.7387 | 0 112 300 101 0602010101 | 121 |
NGC1550 | 64.9080 | 2.4101 | 0 152 150 101 0723800401 0 723 800 501 | 200 |
NGC1600 | 67.9156 | −5.0861 | 0 400 490 101 0400490201 | 81 |
NGC3091 | 150.0591 | −19.6364 | 0 041 180 301 0041180701 | 30 |
NGC5813 | 225.2969 | 1.7019 | 0 302 460 101 0554680201 0 554 680 301 0554680401 | 170 |
NGC5846 | 226.6223 | 1.6048 | 0 021 540 101 0021540501 0 045 340 101 0723800101 0 723 800 201 | 228 |
MRK1216 | 127.1964 | −6.9402 | 0 822 960 101 0822960201 | 235 |
ZW3146 | 155.9147 | 4.1866 | 0 108 670 101 0108670401 0 605 540 201 0605540301 | 240 |
Target . | RA . | Dec. . | OBSIDs . | Exposure (ks) . |
---|---|---|---|---|
2A0335 | 54.6691 | 9.9697 | 0 109 870 101 0109870201 0 147 800 201 | 145 |
A85 | 10.4601 | −9.3031 | 0 065 140 101 0723802101 0 723 802 201 | 215 |
A496 | 68.4074 | −13.2619 | 0 135 120 201 0506260301 0 506 260 401 | 162 |
A2597 | 351.3321 | −12.1243 | 0 108 460 201 0147330101 0 723 801 601 0723801701 | 257 |
A2199 | 247.1594 | 39.5512 | 0 008 030 201 0008030301 0 008 030 601 0723801101 0 723 801 201 | 137 |
M87 | 187.7059 | 12.3911 | 0 114 120 101 0200920101 0 803 670 501 0803670601 0 803 671 001 0803671101 | 430 |
NGC1399 | 54.6210 | −35.4505 | 0 012 830 101 0400620101 | 139 |
NGC720 | 28.2519 | −13.7387 | 0 112 300 101 0602010101 | 121 |
NGC1550 | 64.9080 | 2.4101 | 0 152 150 101 0723800401 0 723 800 501 | 200 |
NGC1600 | 67.9156 | −5.0861 | 0 400 490 101 0400490201 | 81 |
NGC3091 | 150.0591 | −19.6364 | 0 041 180 301 0041180701 | 30 |
NGC5813 | 225.2969 | 1.7019 | 0 302 460 101 0554680201 0 554 680 301 0554680401 | 170 |
NGC5846 | 226.6223 | 1.6048 | 0 021 540 101 0021540501 0 045 340 101 0723800101 0 723 800 201 | 228 |
MRK1216 | 127.1964 | −6.9402 | 0 822 960 101 0822960201 | 235 |
ZW3146 | 155.9147 | 4.1866 | 0 108 670 101 0108670401 0 605 540 201 0605540301 | 240 |
Observed targets, giving the used source position (deg; J2000), observation identifiers, and average cleaned exposure of the RGS cameras.
Target . | RA . | Dec. . | OBSIDs . | Exposure (ks) . |
---|---|---|---|---|
2A0335 | 54.6691 | 9.9697 | 0 109 870 101 0109870201 0 147 800 201 | 145 |
A85 | 10.4601 | −9.3031 | 0 065 140 101 0723802101 0 723 802 201 | 215 |
A496 | 68.4074 | −13.2619 | 0 135 120 201 0506260301 0 506 260 401 | 162 |
A2597 | 351.3321 | −12.1243 | 0 108 460 201 0147330101 0 723 801 601 0723801701 | 257 |
A2199 | 247.1594 | 39.5512 | 0 008 030 201 0008030301 0 008 030 601 0723801101 0 723 801 201 | 137 |
M87 | 187.7059 | 12.3911 | 0 114 120 101 0200920101 0 803 670 501 0803670601 0 803 671 001 0803671101 | 430 |
NGC1399 | 54.6210 | −35.4505 | 0 012 830 101 0400620101 | 139 |
NGC720 | 28.2519 | −13.7387 | 0 112 300 101 0602010101 | 121 |
NGC1550 | 64.9080 | 2.4101 | 0 152 150 101 0723800401 0 723 800 501 | 200 |
NGC1600 | 67.9156 | −5.0861 | 0 400 490 101 0400490201 | 81 |
NGC3091 | 150.0591 | −19.6364 | 0 041 180 301 0041180701 | 30 |
NGC5813 | 225.2969 | 1.7019 | 0 302 460 101 0554680201 0 554 680 301 0554680401 | 170 |
NGC5846 | 226.6223 | 1.6048 | 0 021 540 101 0021540501 0 045 340 101 0723800101 0 723 800 201 | 228 |
MRK1216 | 127.1964 | −6.9402 | 0 822 960 101 0822960201 | 235 |
ZW3146 | 155.9147 | 4.1866 | 0 108 670 101 0108670401 0 605 540 201 0605540301 | 240 |
Target . | RA . | Dec. . | OBSIDs . | Exposure (ks) . |
---|---|---|---|---|
2A0335 | 54.6691 | 9.9697 | 0 109 870 101 0109870201 0 147 800 201 | 145 |
A85 | 10.4601 | −9.3031 | 0 065 140 101 0723802101 0 723 802 201 | 215 |
A496 | 68.4074 | −13.2619 | 0 135 120 201 0506260301 0 506 260 401 | 162 |
A2597 | 351.3321 | −12.1243 | 0 108 460 201 0147330101 0 723 801 601 0723801701 | 257 |
A2199 | 247.1594 | 39.5512 | 0 008 030 201 0008030301 0 008 030 601 0723801101 0 723 801 201 | 137 |
M87 | 187.7059 | 12.3911 | 0 114 120 101 0200920101 0 803 670 501 0803670601 0 803 671 001 0803671101 | 430 |
NGC1399 | 54.6210 | −35.4505 | 0 012 830 101 0400620101 | 139 |
NGC720 | 28.2519 | −13.7387 | 0 112 300 101 0602010101 | 121 |
NGC1550 | 64.9080 | 2.4101 | 0 152 150 101 0723800401 0 723 800 501 | 200 |
NGC1600 | 67.9156 | −5.0861 | 0 400 490 101 0400490201 | 81 |
NGC3091 | 150.0591 | −19.6364 | 0 041 180 301 0041180701 | 30 |
NGC5813 | 225.2969 | 1.7019 | 0 302 460 101 0554680201 0 554 680 301 0554680401 | 170 |
NGC5846 | 226.6223 | 1.6048 | 0 021 540 101 0021540501 0 045 340 101 0723800101 0 723 800 201 | 228 |
MRK1216 | 127.1964 | −6.9402 | 0 822 960 101 0822960201 | 235 |
ZW3146 | 155.9147 | 4.1866 | 0 108 670 101 0108670401 0 605 540 201 0605540301 | 240 |
The spectral model used is tbabs(gsmooth*apec + gsmooth(partcov*mlayerz)mkcflow). The intrinsic absorption model mlayerz (see HCFII for details) represents a sequence of interleaved emission and absorption layers with a total column density NH listed in Table 2. tbabs is the Galactic absorption in the direction of the target apec is a constant temperature thermal emission model which represents the outer cluster gas. Its temperature is also used as the hotter temperature in the cooling flow mkcflow model. partcov enables the measurement of the total mass cooling rate of both unabsorbed and absorbed components. A covering fraction of one means that all the cooling flow component is absorbed and if zero then none is absorbed. The model assumes no particular geometry for the absorbed and unabsorbed components. It does assume that all absorbed components are identical. The minimum temperature of the cooling flow model is set at 0.1 keV.
Spectral fitting results. The units of column density NH are |$10^{22}\hbox{${\rm \, cm}^{-2}\, $}$|, the temperature kT of the apec and maximum of mkcflow component kT (which are the same) is in |${\rm \, keV}$|, Z is abundance relative to solar, and |$\dot{M}$| is in |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|. |$\dot{M}_{\rm u}$| is the uncovered rate (i.e. with no absorption). (f) means that a parameter is fixed. All uncertainties correspond to the 90 per cent confidence level.
Cluster . | NH . | kT . | Z . | z . | Norm . | CFrac . | |${N_{\rm H}}^{^{\prime }}$| . | |$\dot{M}$| . | χ2/dof . | |$\dot{M}_{\rm u}$| . |
---|---|---|---|---|---|---|---|---|---|---|
. | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|${\rm \, keV}$|) . | (Z⊙) . | . | . | . | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | . |
2A0335 | 0.26 | |$1.85^{+0.02}_{-0.01}$| | 0.32 ± 0.01 | 0.035 | 4.9e-2 | 1 | 0.32 | |$86^{+9}_{-7}$| | 1359/1224 | 48 ± 4 |
A85 | 0.029 | 3.6 ± 0.4 | |$0.29^{+0.07}_{-0.003}$| | 0.056 | 2e-2 | 0.95 | 1 | |$25^{+60}_{-20}$| | 1406/1327 | |$14.4^{+3}_{-5}$| |
A496 | 0.05 | 2.49 | 0.32 | 0.033 | 2.1e-2 | 0.95 | 2.2 | |$22.7^{+7.5}_{-10}$| | 1406/1299 | <2.1 |
A2597 | 0.023 | 3.3 ± 0.13 | 0.35 ± 0.03 | 0.082 | 1.5e-2 | 0.71 | 5 | |$67^{+139}_{-19}$| | 1457/1364 | |$19.5^{+6}_{-4}$| |
A2199 | 0.008 | |$2.96^{+0.16}_{-0.08}$| | |$0.29^{+0.14}_{-0.02}$| | 0.0296 | 2.42e-2 | 0.95 | 0.46 | |$5.6^{+4.5}_{2.3}$| | 1409/1312 | 2.3 ± 1.5 |
M87 | 0.018 | 1.47 | 0.23 | 4.28e-3 | 7.3e-2 | 1. | 0.2 | |$0.8^{+0.2}_{-0.15}$| | 1875/992 | 0.5 ± 0.04 |
NGC1399 | 0.014 | 1.08 ± 0.01 | 0.27 ± 0.02 | 5.5e-3 | 3.2e-3 | 1. | 2.67 | |$3.3^{+3.4}_{-0.7}$| | 1086/766 | 0.15 ± 0.03 |
NGC720 | 0.14 | 0.61 ± 0.02 | 0.11 ± 0.02 | 0.0065 | 7.9e-4 | 0.95 | 1.13 | |$1.3^{+2.6}_{-0.95}$| | 232/214 | 0.25 ± 0.1 |
NGC1550 | 0.114 | 1.26 ± 0.01 | 0.26 ± 0.02 | 0.0132 | 7.5e-3 | 1. | 1.33 | |$3.6^{+2.4}_{-1.0}$| | 1115/1051 | 0.55 ± 0.2 |
NGC1600 | 0.04 | |$1.33^{+0.02}_{-0.07}$| | |$0.12^{+0.07}_{-0.04}$| | 0.0163 | 1.1e-3 | 0.65 | 1.3 | |$0.81^{+7}_{-0.7}$| | 161/146 | <0.6 |
NGC3091 | 0.013 | 0.01 ± 0.03 | 0.09 ± 0.02 | 0.013 | 1.37e-3 | 0.95 | 3.16 | 8.5 ± 6.6 | 92/78 | <1 |
NGC5813 | 0.043 | 0.72 ± 0.01 | 0.48 ± 0.08 | 0.0063 | 1.8e-3 | 1 | 0.5 | |$2.0^{+5}_{-1}$| | 1227/881 | 2.1 |
NGC5846 | 0.043 | 0.83 ± 0.01 | 0.33 ± 0.02 | 0.0053 | 1.5e-3 | 1 | 3.2 | |$1.29^{+0.3}_{-0.15}$| | 839/617 | 0.8 ± 0.5 |
MRK1216 | 0.04 | 0.7 ± 0.14 | 0.33 ± 0.1 | 0.021 | 4.e-4 | 0.95 | 2 | 9.7 ± 2.7 | 336/336 | 1.3 ± 0.4 |
ZW3146 | 0.024 | 5.7 ± 0.45 | 0.4f | 0.288 | 8.8e-3 | 0.8 | 3 | |$895^{+1000}_{-300}$| | 1196/1169 | |$230^{+50}_{-80}$| |
Cluster . | NH . | kT . | Z . | z . | Norm . | CFrac . | |${N_{\rm H}}^{^{\prime }}$| . | |$\dot{M}$| . | χ2/dof . | |$\dot{M}_{\rm u}$| . |
---|---|---|---|---|---|---|---|---|---|---|
. | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|${\rm \, keV}$|) . | (Z⊙) . | . | . | . | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | . |
2A0335 | 0.26 | |$1.85^{+0.02}_{-0.01}$| | 0.32 ± 0.01 | 0.035 | 4.9e-2 | 1 | 0.32 | |$86^{+9}_{-7}$| | 1359/1224 | 48 ± 4 |
A85 | 0.029 | 3.6 ± 0.4 | |$0.29^{+0.07}_{-0.003}$| | 0.056 | 2e-2 | 0.95 | 1 | |$25^{+60}_{-20}$| | 1406/1327 | |$14.4^{+3}_{-5}$| |
A496 | 0.05 | 2.49 | 0.32 | 0.033 | 2.1e-2 | 0.95 | 2.2 | |$22.7^{+7.5}_{-10}$| | 1406/1299 | <2.1 |
A2597 | 0.023 | 3.3 ± 0.13 | 0.35 ± 0.03 | 0.082 | 1.5e-2 | 0.71 | 5 | |$67^{+139}_{-19}$| | 1457/1364 | |$19.5^{+6}_{-4}$| |
A2199 | 0.008 | |$2.96^{+0.16}_{-0.08}$| | |$0.29^{+0.14}_{-0.02}$| | 0.0296 | 2.42e-2 | 0.95 | 0.46 | |$5.6^{+4.5}_{2.3}$| | 1409/1312 | 2.3 ± 1.5 |
M87 | 0.018 | 1.47 | 0.23 | 4.28e-3 | 7.3e-2 | 1. | 0.2 | |$0.8^{+0.2}_{-0.15}$| | 1875/992 | 0.5 ± 0.04 |
NGC1399 | 0.014 | 1.08 ± 0.01 | 0.27 ± 0.02 | 5.5e-3 | 3.2e-3 | 1. | 2.67 | |$3.3^{+3.4}_{-0.7}$| | 1086/766 | 0.15 ± 0.03 |
NGC720 | 0.14 | 0.61 ± 0.02 | 0.11 ± 0.02 | 0.0065 | 7.9e-4 | 0.95 | 1.13 | |$1.3^{+2.6}_{-0.95}$| | 232/214 | 0.25 ± 0.1 |
NGC1550 | 0.114 | 1.26 ± 0.01 | 0.26 ± 0.02 | 0.0132 | 7.5e-3 | 1. | 1.33 | |$3.6^{+2.4}_{-1.0}$| | 1115/1051 | 0.55 ± 0.2 |
NGC1600 | 0.04 | |$1.33^{+0.02}_{-0.07}$| | |$0.12^{+0.07}_{-0.04}$| | 0.0163 | 1.1e-3 | 0.65 | 1.3 | |$0.81^{+7}_{-0.7}$| | 161/146 | <0.6 |
NGC3091 | 0.013 | 0.01 ± 0.03 | 0.09 ± 0.02 | 0.013 | 1.37e-3 | 0.95 | 3.16 | 8.5 ± 6.6 | 92/78 | <1 |
NGC5813 | 0.043 | 0.72 ± 0.01 | 0.48 ± 0.08 | 0.0063 | 1.8e-3 | 1 | 0.5 | |$2.0^{+5}_{-1}$| | 1227/881 | 2.1 |
NGC5846 | 0.043 | 0.83 ± 0.01 | 0.33 ± 0.02 | 0.0053 | 1.5e-3 | 1 | 3.2 | |$1.29^{+0.3}_{-0.15}$| | 839/617 | 0.8 ± 0.5 |
MRK1216 | 0.04 | 0.7 ± 0.14 | 0.33 ± 0.1 | 0.021 | 4.e-4 | 0.95 | 2 | 9.7 ± 2.7 | 336/336 | 1.3 ± 0.4 |
ZW3146 | 0.024 | 5.7 ± 0.45 | 0.4f | 0.288 | 8.8e-3 | 0.8 | 3 | |$895^{+1000}_{-300}$| | 1196/1169 | |$230^{+50}_{-80}$| |
Spectral fitting results. The units of column density NH are |$10^{22}\hbox{${\rm \, cm}^{-2}\, $}$|, the temperature kT of the apec and maximum of mkcflow component kT (which are the same) is in |${\rm \, keV}$|, Z is abundance relative to solar, and |$\dot{M}$| is in |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|. |$\dot{M}_{\rm u}$| is the uncovered rate (i.e. with no absorption). (f) means that a parameter is fixed. All uncertainties correspond to the 90 per cent confidence level.
Cluster . | NH . | kT . | Z . | z . | Norm . | CFrac . | |${N_{\rm H}}^{^{\prime }}$| . | |$\dot{M}$| . | χ2/dof . | |$\dot{M}_{\rm u}$| . |
---|---|---|---|---|---|---|---|---|---|---|
. | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|${\rm \, keV}$|) . | (Z⊙) . | . | . | . | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | . |
2A0335 | 0.26 | |$1.85^{+0.02}_{-0.01}$| | 0.32 ± 0.01 | 0.035 | 4.9e-2 | 1 | 0.32 | |$86^{+9}_{-7}$| | 1359/1224 | 48 ± 4 |
A85 | 0.029 | 3.6 ± 0.4 | |$0.29^{+0.07}_{-0.003}$| | 0.056 | 2e-2 | 0.95 | 1 | |$25^{+60}_{-20}$| | 1406/1327 | |$14.4^{+3}_{-5}$| |
A496 | 0.05 | 2.49 | 0.32 | 0.033 | 2.1e-2 | 0.95 | 2.2 | |$22.7^{+7.5}_{-10}$| | 1406/1299 | <2.1 |
A2597 | 0.023 | 3.3 ± 0.13 | 0.35 ± 0.03 | 0.082 | 1.5e-2 | 0.71 | 5 | |$67^{+139}_{-19}$| | 1457/1364 | |$19.5^{+6}_{-4}$| |
A2199 | 0.008 | |$2.96^{+0.16}_{-0.08}$| | |$0.29^{+0.14}_{-0.02}$| | 0.0296 | 2.42e-2 | 0.95 | 0.46 | |$5.6^{+4.5}_{2.3}$| | 1409/1312 | 2.3 ± 1.5 |
M87 | 0.018 | 1.47 | 0.23 | 4.28e-3 | 7.3e-2 | 1. | 0.2 | |$0.8^{+0.2}_{-0.15}$| | 1875/992 | 0.5 ± 0.04 |
NGC1399 | 0.014 | 1.08 ± 0.01 | 0.27 ± 0.02 | 5.5e-3 | 3.2e-3 | 1. | 2.67 | |$3.3^{+3.4}_{-0.7}$| | 1086/766 | 0.15 ± 0.03 |
NGC720 | 0.14 | 0.61 ± 0.02 | 0.11 ± 0.02 | 0.0065 | 7.9e-4 | 0.95 | 1.13 | |$1.3^{+2.6}_{-0.95}$| | 232/214 | 0.25 ± 0.1 |
NGC1550 | 0.114 | 1.26 ± 0.01 | 0.26 ± 0.02 | 0.0132 | 7.5e-3 | 1. | 1.33 | |$3.6^{+2.4}_{-1.0}$| | 1115/1051 | 0.55 ± 0.2 |
NGC1600 | 0.04 | |$1.33^{+0.02}_{-0.07}$| | |$0.12^{+0.07}_{-0.04}$| | 0.0163 | 1.1e-3 | 0.65 | 1.3 | |$0.81^{+7}_{-0.7}$| | 161/146 | <0.6 |
NGC3091 | 0.013 | 0.01 ± 0.03 | 0.09 ± 0.02 | 0.013 | 1.37e-3 | 0.95 | 3.16 | 8.5 ± 6.6 | 92/78 | <1 |
NGC5813 | 0.043 | 0.72 ± 0.01 | 0.48 ± 0.08 | 0.0063 | 1.8e-3 | 1 | 0.5 | |$2.0^{+5}_{-1}$| | 1227/881 | 2.1 |
NGC5846 | 0.043 | 0.83 ± 0.01 | 0.33 ± 0.02 | 0.0053 | 1.5e-3 | 1 | 3.2 | |$1.29^{+0.3}_{-0.15}$| | 839/617 | 0.8 ± 0.5 |
MRK1216 | 0.04 | 0.7 ± 0.14 | 0.33 ± 0.1 | 0.021 | 4.e-4 | 0.95 | 2 | 9.7 ± 2.7 | 336/336 | 1.3 ± 0.4 |
ZW3146 | 0.024 | 5.7 ± 0.45 | 0.4f | 0.288 | 8.8e-3 | 0.8 | 3 | |$895^{+1000}_{-300}$| | 1196/1169 | |$230^{+50}_{-80}$| |
Cluster . | NH . | kT . | Z . | z . | Norm . | CFrac . | |${N_{\rm H}}^{^{\prime }}$| . | |$\dot{M}$| . | χ2/dof . | |$\dot{M}_{\rm u}$| . |
---|---|---|---|---|---|---|---|---|---|---|
. | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|${\rm \, keV}$|) . | (Z⊙) . | . | . | . | (|$10^{22}\hbox{${\rm \, cm}^2\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | . |
2A0335 | 0.26 | |$1.85^{+0.02}_{-0.01}$| | 0.32 ± 0.01 | 0.035 | 4.9e-2 | 1 | 0.32 | |$86^{+9}_{-7}$| | 1359/1224 | 48 ± 4 |
A85 | 0.029 | 3.6 ± 0.4 | |$0.29^{+0.07}_{-0.003}$| | 0.056 | 2e-2 | 0.95 | 1 | |$25^{+60}_{-20}$| | 1406/1327 | |$14.4^{+3}_{-5}$| |
A496 | 0.05 | 2.49 | 0.32 | 0.033 | 2.1e-2 | 0.95 | 2.2 | |$22.7^{+7.5}_{-10}$| | 1406/1299 | <2.1 |
A2597 | 0.023 | 3.3 ± 0.13 | 0.35 ± 0.03 | 0.082 | 1.5e-2 | 0.71 | 5 | |$67^{+139}_{-19}$| | 1457/1364 | |$19.5^{+6}_{-4}$| |
A2199 | 0.008 | |$2.96^{+0.16}_{-0.08}$| | |$0.29^{+0.14}_{-0.02}$| | 0.0296 | 2.42e-2 | 0.95 | 0.46 | |$5.6^{+4.5}_{2.3}$| | 1409/1312 | 2.3 ± 1.5 |
M87 | 0.018 | 1.47 | 0.23 | 4.28e-3 | 7.3e-2 | 1. | 0.2 | |$0.8^{+0.2}_{-0.15}$| | 1875/992 | 0.5 ± 0.04 |
NGC1399 | 0.014 | 1.08 ± 0.01 | 0.27 ± 0.02 | 5.5e-3 | 3.2e-3 | 1. | 2.67 | |$3.3^{+3.4}_{-0.7}$| | 1086/766 | 0.15 ± 0.03 |
NGC720 | 0.14 | 0.61 ± 0.02 | 0.11 ± 0.02 | 0.0065 | 7.9e-4 | 0.95 | 1.13 | |$1.3^{+2.6}_{-0.95}$| | 232/214 | 0.25 ± 0.1 |
NGC1550 | 0.114 | 1.26 ± 0.01 | 0.26 ± 0.02 | 0.0132 | 7.5e-3 | 1. | 1.33 | |$3.6^{+2.4}_{-1.0}$| | 1115/1051 | 0.55 ± 0.2 |
NGC1600 | 0.04 | |$1.33^{+0.02}_{-0.07}$| | |$0.12^{+0.07}_{-0.04}$| | 0.0163 | 1.1e-3 | 0.65 | 1.3 | |$0.81^{+7}_{-0.7}$| | 161/146 | <0.6 |
NGC3091 | 0.013 | 0.01 ± 0.03 | 0.09 ± 0.02 | 0.013 | 1.37e-3 | 0.95 | 3.16 | 8.5 ± 6.6 | 92/78 | <1 |
NGC5813 | 0.043 | 0.72 ± 0.01 | 0.48 ± 0.08 | 0.0063 | 1.8e-3 | 1 | 0.5 | |$2.0^{+5}_{-1}$| | 1227/881 | 2.1 |
NGC5846 | 0.043 | 0.83 ± 0.01 | 0.33 ± 0.02 | 0.0053 | 1.5e-3 | 1 | 3.2 | |$1.29^{+0.3}_{-0.15}$| | 839/617 | 0.8 ± 0.5 |
MRK1216 | 0.04 | 0.7 ± 0.14 | 0.33 ± 0.1 | 0.021 | 4.e-4 | 0.95 | 2 | 9.7 ± 2.7 | 336/336 | 1.3 ± 0.4 |
ZW3146 | 0.024 | 5.7 ± 0.45 | 0.4f | 0.288 | 8.8e-3 | 0.8 | 3 | |$895^{+1000}_{-300}$| | 1196/1169 | |$230^{+50}_{-80}$| |
The spectra are shown in the appendix as Figs B1–B15, together with contour plots of absorbed mass cooling rate (|$\dot{M}_{\rm a}$|) versus intrinsic column density NH and covering fraction.
Since the RGS is a slitless spectrometer (den Herder et al. 2001), there is some blurring of the energy scale associated with extended sources. This is included in the spectral model by smoothing the spectral components with separate Gaussian kernels for the outer apec component and inner HCF. When making the contour plots for the less bright objects, we often needed to freeze the smoothing parameters to their best-fitting values in order to have convergence. Detailed spectral results are given in Table 2 and are compared with data from other wavebands in Table 3.
Relevant cluster properties. See subsections of Appendix A for individual object references. A dash indicates lack of data.
Cluster . | L(FIR) . | La . | |$\dot{M}$| . | L(H α) . | MCO . | MBH . |
---|---|---|---|---|---|---|
. | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . |
2A0335 | 4e43 | 2.1e43 | 86 | 8e41 | 1.1e9 | – |
A85 | 2.8e43 | 9.9e42 | 23 | – | – | 4e10 |
A496 | – | 9.6e42 | 23 | 5e40 | – | – |
A2597 | 6.5e43 | 2.1e43 | 67 | 3e42 | 2.3e9 | – |
A2199 | - | 1.5e42 | 5.6 | 3.5e40 | - | 4e9 |
M87 | 5.0e41 | 1.6e41 | 0.8 | 1.9e40 | - | 6.5e9 |
NGC1399 | - | 7.4e41 | 3.3 | 1e39 | 1e9 | |
NGC720 | – | 1.5e41 | 1.0 | – | 1.1e7 | – |
NGC1550 | – | 8.7e41 | 1.5 | – | – | 4.5e9 |
NGC1600 | – | 1.3e41 | 0.8 | 4e39 | – | 1.7e10 |
NGC3091 | – | 1.6e42 | 8.5 | – | – | 3.6e9 |
NGC5813 | 1.1e42 | 5.9e41 | 2.0 | 1.6e40 | – | – |
NGC5846 | 6.2e41 | 2.0e41 | 1.3 | 2.5e40 | 2e6 | – |
MRK1216 | – | 1.3e41 | 9.7 | – | – | 4.9e9 |
ZW3146 | 1.0e45 | 6.3e44 | 1570 | 6e42 | 5e10 | – |
NGC5044 | 3.0e42 | 3.6e42 | 20 | 7.0e40 | 1.5e8 | |
Sersic 159 | 7.3e42 | 2.5e42 | 10 | 2.0e41 | 1.1e9 | |
A262 | 8.0e42 | 2.1e42 | 7 | 9.4e40 | 4.0e8 | |
A2052 | 8.3e42 | 4.4e42 | 15 | 6e40 | 2.8e8 | |
RXJ0821 | 4.5e44 | 7.8e42 | 40 | 3.0e41 | 3.9e10 | |
RXJ1532 | 2.3e45 | 2.0e44 | 1000 | 3e42 | 8.7e10 | |
MACS1931 | 5.6e45 | 4.6e44 | 1000 | 2e42 | 9.0e10 | |
Phoenix Cluster | 3.7e46 | 3.3e44 | 2000 | 8.5e43 | 2e10 | |
M84 | 1.0e42 | 3.3e41 | 2.0 | 4.0e39 | <1.8e7 | |
M49 | 1.2e42 | 2.0e41 | 1.0 | 5.8e39 | <1.4e7 | |
Centaurus | 3.2e42 | 3.6e42 | 15 | 1.7e40 | 1.0e8 | |
Perseus | 5.6e44 | 5.8e42 | 50 | 3.2e42 | 2.0e10 | |
A1835 | 3.2e45 | 5.2e43 | 400 | 4.4e42 | 5.0e10 | |
RXJ1504 | - | 1.9e44 | 520 | 3.2e43 | 1.9e10 |
Cluster . | L(FIR) . | La . | |$\dot{M}$| . | L(H α) . | MCO . | MBH . |
---|---|---|---|---|---|---|
. | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . |
2A0335 | 4e43 | 2.1e43 | 86 | 8e41 | 1.1e9 | – |
A85 | 2.8e43 | 9.9e42 | 23 | – | – | 4e10 |
A496 | – | 9.6e42 | 23 | 5e40 | – | – |
A2597 | 6.5e43 | 2.1e43 | 67 | 3e42 | 2.3e9 | – |
A2199 | - | 1.5e42 | 5.6 | 3.5e40 | - | 4e9 |
M87 | 5.0e41 | 1.6e41 | 0.8 | 1.9e40 | - | 6.5e9 |
NGC1399 | - | 7.4e41 | 3.3 | 1e39 | 1e9 | |
NGC720 | – | 1.5e41 | 1.0 | – | 1.1e7 | – |
NGC1550 | – | 8.7e41 | 1.5 | – | – | 4.5e9 |
NGC1600 | – | 1.3e41 | 0.8 | 4e39 | – | 1.7e10 |
NGC3091 | – | 1.6e42 | 8.5 | – | – | 3.6e9 |
NGC5813 | 1.1e42 | 5.9e41 | 2.0 | 1.6e40 | – | – |
NGC5846 | 6.2e41 | 2.0e41 | 1.3 | 2.5e40 | 2e6 | – |
MRK1216 | – | 1.3e41 | 9.7 | – | – | 4.9e9 |
ZW3146 | 1.0e45 | 6.3e44 | 1570 | 6e42 | 5e10 | – |
NGC5044 | 3.0e42 | 3.6e42 | 20 | 7.0e40 | 1.5e8 | |
Sersic 159 | 7.3e42 | 2.5e42 | 10 | 2.0e41 | 1.1e9 | |
A262 | 8.0e42 | 2.1e42 | 7 | 9.4e40 | 4.0e8 | |
A2052 | 8.3e42 | 4.4e42 | 15 | 6e40 | 2.8e8 | |
RXJ0821 | 4.5e44 | 7.8e42 | 40 | 3.0e41 | 3.9e10 | |
RXJ1532 | 2.3e45 | 2.0e44 | 1000 | 3e42 | 8.7e10 | |
MACS1931 | 5.6e45 | 4.6e44 | 1000 | 2e42 | 9.0e10 | |
Phoenix Cluster | 3.7e46 | 3.3e44 | 2000 | 8.5e43 | 2e10 | |
M84 | 1.0e42 | 3.3e41 | 2.0 | 4.0e39 | <1.8e7 | |
M49 | 1.2e42 | 2.0e41 | 1.0 | 5.8e39 | <1.4e7 | |
Centaurus | 3.2e42 | 3.6e42 | 15 | 1.7e40 | 1.0e8 | |
Perseus | 5.6e44 | 5.8e42 | 50 | 3.2e42 | 2.0e10 | |
A1835 | 3.2e45 | 5.2e43 | 400 | 4.4e42 | 5.0e10 | |
RXJ1504 | - | 1.9e44 | 520 | 3.2e43 | 1.9e10 |
Relevant cluster properties. See subsections of Appendix A for individual object references. A dash indicates lack of data.
Cluster . | L(FIR) . | La . | |$\dot{M}$| . | L(H α) . | MCO . | MBH . |
---|---|---|---|---|---|---|
. | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . |
2A0335 | 4e43 | 2.1e43 | 86 | 8e41 | 1.1e9 | – |
A85 | 2.8e43 | 9.9e42 | 23 | – | – | 4e10 |
A496 | – | 9.6e42 | 23 | 5e40 | – | – |
A2597 | 6.5e43 | 2.1e43 | 67 | 3e42 | 2.3e9 | – |
A2199 | - | 1.5e42 | 5.6 | 3.5e40 | - | 4e9 |
M87 | 5.0e41 | 1.6e41 | 0.8 | 1.9e40 | - | 6.5e9 |
NGC1399 | - | 7.4e41 | 3.3 | 1e39 | 1e9 | |
NGC720 | – | 1.5e41 | 1.0 | – | 1.1e7 | – |
NGC1550 | – | 8.7e41 | 1.5 | – | – | 4.5e9 |
NGC1600 | – | 1.3e41 | 0.8 | 4e39 | – | 1.7e10 |
NGC3091 | – | 1.6e42 | 8.5 | – | – | 3.6e9 |
NGC5813 | 1.1e42 | 5.9e41 | 2.0 | 1.6e40 | – | – |
NGC5846 | 6.2e41 | 2.0e41 | 1.3 | 2.5e40 | 2e6 | – |
MRK1216 | – | 1.3e41 | 9.7 | – | – | 4.9e9 |
ZW3146 | 1.0e45 | 6.3e44 | 1570 | 6e42 | 5e10 | – |
NGC5044 | 3.0e42 | 3.6e42 | 20 | 7.0e40 | 1.5e8 | |
Sersic 159 | 7.3e42 | 2.5e42 | 10 | 2.0e41 | 1.1e9 | |
A262 | 8.0e42 | 2.1e42 | 7 | 9.4e40 | 4.0e8 | |
A2052 | 8.3e42 | 4.4e42 | 15 | 6e40 | 2.8e8 | |
RXJ0821 | 4.5e44 | 7.8e42 | 40 | 3.0e41 | 3.9e10 | |
RXJ1532 | 2.3e45 | 2.0e44 | 1000 | 3e42 | 8.7e10 | |
MACS1931 | 5.6e45 | 4.6e44 | 1000 | 2e42 | 9.0e10 | |
Phoenix Cluster | 3.7e46 | 3.3e44 | 2000 | 8.5e43 | 2e10 | |
M84 | 1.0e42 | 3.3e41 | 2.0 | 4.0e39 | <1.8e7 | |
M49 | 1.2e42 | 2.0e41 | 1.0 | 5.8e39 | <1.4e7 | |
Centaurus | 3.2e42 | 3.6e42 | 15 | 1.7e40 | 1.0e8 | |
Perseus | 5.6e44 | 5.8e42 | 50 | 3.2e42 | 2.0e10 | |
A1835 | 3.2e45 | 5.2e43 | 400 | 4.4e42 | 5.0e10 | |
RXJ1504 | - | 1.9e44 | 520 | 3.2e43 | 1.9e10 |
Cluster . | L(FIR) . | La . | |$\dot{M}$| . | L(H α) . | MCO . | MBH . |
---|---|---|---|---|---|---|
. | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|) . | (|$\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . | (|$\hbox{$\rm \, {\rm M}_{\odot }$}$|) . |
2A0335 | 4e43 | 2.1e43 | 86 | 8e41 | 1.1e9 | – |
A85 | 2.8e43 | 9.9e42 | 23 | – | – | 4e10 |
A496 | – | 9.6e42 | 23 | 5e40 | – | – |
A2597 | 6.5e43 | 2.1e43 | 67 | 3e42 | 2.3e9 | – |
A2199 | - | 1.5e42 | 5.6 | 3.5e40 | - | 4e9 |
M87 | 5.0e41 | 1.6e41 | 0.8 | 1.9e40 | - | 6.5e9 |
NGC1399 | - | 7.4e41 | 3.3 | 1e39 | 1e9 | |
NGC720 | – | 1.5e41 | 1.0 | – | 1.1e7 | – |
NGC1550 | – | 8.7e41 | 1.5 | – | – | 4.5e9 |
NGC1600 | – | 1.3e41 | 0.8 | 4e39 | – | 1.7e10 |
NGC3091 | – | 1.6e42 | 8.5 | – | – | 3.6e9 |
NGC5813 | 1.1e42 | 5.9e41 | 2.0 | 1.6e40 | – | – |
NGC5846 | 6.2e41 | 2.0e41 | 1.3 | 2.5e40 | 2e6 | – |
MRK1216 | – | 1.3e41 | 9.7 | – | – | 4.9e9 |
ZW3146 | 1.0e45 | 6.3e44 | 1570 | 6e42 | 5e10 | – |
NGC5044 | 3.0e42 | 3.6e42 | 20 | 7.0e40 | 1.5e8 | |
Sersic 159 | 7.3e42 | 2.5e42 | 10 | 2.0e41 | 1.1e9 | |
A262 | 8.0e42 | 2.1e42 | 7 | 9.4e40 | 4.0e8 | |
A2052 | 8.3e42 | 4.4e42 | 15 | 6e40 | 2.8e8 | |
RXJ0821 | 4.5e44 | 7.8e42 | 40 | 3.0e41 | 3.9e10 | |
RXJ1532 | 2.3e45 | 2.0e44 | 1000 | 3e42 | 8.7e10 | |
MACS1931 | 5.6e45 | 4.6e44 | 1000 | 2e42 | 9.0e10 | |
Phoenix Cluster | 3.7e46 | 3.3e44 | 2000 | 8.5e43 | 2e10 | |
M84 | 1.0e42 | 3.3e41 | 2.0 | 4.0e39 | <1.8e7 | |
M49 | 1.2e42 | 2.0e41 | 1.0 | 5.8e39 | <1.4e7 | |
Centaurus | 3.2e42 | 3.6e42 | 15 | 1.7e40 | 1.0e8 | |
Perseus | 5.6e44 | 5.8e42 | 50 | 3.2e42 | 2.0e10 | |
A1835 | 3.2e45 | 5.2e43 | 400 | 4.4e42 | 5.0e10 | |
RXJ1504 | - | 1.9e44 | 520 | 3.2e43 | 1.9e10 |
3 THE SPECTRAL RESULTS
As noted in HCFII, χ2-space for the HCF model is often corrugated, which can lead to complex contour plots. We are using a very simplistic model and a real hidden cooling flow is expected to be far more complicated in both space and column density. RGS spectra provide no more than a rough average over the inner arcmin of the target source.
A source like ZW3146, where there is a large continuum fraction, can have a very uncertain abundance Z, with it anticorrelating with the mass cooling rate (see Fig. 1). In this case we fix it at Z = 0.4.

RGS spectrum of ZW3146 with HCF component shown in red and mccall component dotted, mass cooling rate in |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| versus total column density in units of |$10^{22}\hbox{${\rm \, cm}^2\, $}$|, mass cooling rate versus covering fraction of the HCF component, and contours at 68 per cent (red), 90 per cent (green), and 99 per cent (blue).
Of the 15 sources studied here, all but three require a best-fitting covering fraction of 0.95 or more. This emphasizes that they are indeed ‘hidden’. The intrinsic column densities range from 2 × 1021 to |$3\times 10^{22}\hbox{${\rm \, cm}^{-2}\, $}$|.
We also refit the spectra with the covering fraction set to zero, in order to determine the mass cooling rate if there is no absorption, |$\dot{M}_{\rm u}$|. This is listed in the last column of Table 2. As expected, it is generally very low, but quite large for 2A0335. The lowest χ2 value for this no absorption case is however 14 above that for the best-fitting HCF mlayerz model, which is therefore the statistically preferred one. The value in the case of A2597 is about what is expected from the HCF model where the covering fraction is about 70 per cent.
When the temperature of the gas is above about 0.4 keV the fraction of the energy emerging in continuum is about 50 per cent and drops below 10 per cent below 0.2 keV (Fig. 2). Most of the flux below 0.4 keV is absorbed away in our HCF fits, meaning the continuum shape plays a significant role in our spectral fit results.

Flux from cooling flow emerging in lines (blue) and continuum (green) (top panel), as a fraction (bottom panel).
We reduced the energy band of the spectra of several lower temperature objects to 12–20A due to broad excess residuals around 10A. These are likely due to the apec component having a (small) spread in temperature.
The absorbed luminosities (La, Table 3) are all less than the far infrared luminosities, where available. This indicates that the energy lost in the cooling flow to absorption is energetically capable of emerging as radiation from dust in the absorbing gas.
Fig. 3 shows the HCF mass cooling rates (in red) compared (where available) with the ‘classical’ rates from X-ray imaging listed by (Hudson et al. 2010) (in black). The mean ratio of hidden to classical rates is 45 per cent, with a range from 4 to 180 per cent. (Most lie between 17 and 58 per cent.) The HCF rates are about |$1\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| for elliptical galaxies, 2–|$20\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| for brightest group galaxies (BGGs), and about 10–|$100\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| for regular BCGs. There are then a group of more distant, exceptionally X-ray luminous, BCGs with 400 to |$\gt 1000\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|.

Mass cooling rates, classical imaging rate from (Hudson et al. 2010) (black), if available, and spectroscopic HCF rate (red). Objects: (1) 2A0335; (2) A85; (3) A496; (4) A2597; (5) S159, (6) A262, (7) A2052; (8) Cen; (9) Per; (10) A2199; (11) NGC1550; and (12) NGC5044. The average ratio of red (HCF) to black (classical) is 0.45.
We suspect that the last group of rare objects may be highly time variable, with peak luminosity followed by a quasar eruption. The regular clusters and elliptical galaxies generally have low-luminosity nuclei, with radio emission from jets that blow bubbles in the intracluster medium. The bubbles and related activity generally lie outside the inner kpc studied here.
4 THE ACCUMULATION OF COOLED GAS
Over a billion years |$10^9\hbox{$\rm \, {\rm M}_{\odot }$}$| of gas will have cooled in a typical elliptical, up to |$10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| in a BGG and up to |$10^{11}\hbox{$\rm \, {\rm M}_{\odot }$}$| in a BCG. These are large values, the higher end of which exceeds the cold molecular masses observed via CO emission in BCGs (Olivares et al. 2019; Russell et al. 2019). It is possible that the mass of molecular gas has been underestimated due to low abundance and an unseen diffuse component, but this is unlikely to make a very large difference.
In HCFI, we considered the following possibilities: (a) continued cooling to invisibility at 3 K, (b) fragmentation and collapse into substellar objects since the Jeans mass is less than |$0.1\hbox{$\rm \, {\rm M}_{\odot }$}$|, (c) outward dragging of cooled clouds by the bubbling process, or (d) cold front formation. The gas and dark matter peaks may be offset by a kpc or more.
We also flagged the similarity in conditions (e.g. gas pressure) of cooled dusty molecular clouds of a BCG core to those in the Crab Nebula. More detailed observational comparisons are warranted.
It is likely that the cooled material is spread over the innermost few kpc of the core. Clear evidence of the dragging out of dust-enriched material from the centre is provided by the peaks in metal abundance seen |$\sim 10{\rm \, kpc}$| from the centre of low-redshift clusters (Panagoulia, Sanders & Fabian 2015; Lakhchaura, Mernier & Werner 2019; Liu et al. 2019).
Detailed measurements of the mass profile of each separate component in a cool core (black hole, dark matter, stars, gas, etc.) will be invaluable in sorting the possibilities out. We now consider whether some small fraction of the very cold clouds and substellar objects can be swallowed by the central black hole in the next section.
4.1 Accretion of fragmented cold matter by the central black hole
We showed in HCFI (Fabian et al. 2022) that, under the high-pressure conditions of an HCF (|$nT\sim 10^{6.5}\!-\!10^{7.5}\hbox{${\rm \, cm}^{-3}{\rm \, K}\, $}$|) and no heating, the gas cools rapidly (time-scale of tens of years) to ∼3 K. The Jeans mass is below about |$0.1\hbox{$\rm \, {\rm M}_{\odot }$}$| (Jura 1977; Ferland, Fabian & Johnstone 1994) and the gas expected to clump and fragment into low mass stars, brown dwarfs, etc., some of which will fall into the black hole emitting little radiation. Exactly how large a fraction will be swallowed depends on how angular momentum is transported outward. The turbulent viscosity of a luminous accretion disc is absent here and a possible path is that the innermost cooled gas forms a thick disc of low-mass stars and cold gas clouds around the black hole. Dynamical gravitational instabilities such as spiral waves2 and bars within bars transport angular momentum outward in non-spherical systems so that some of the matter falls inward (Shlosman, Frank & Begelman 1989; Hopkins & Quataert 2011; Gualandris et al. 2017) to be swallowed directly by the central black hole without a standard accretion disc forming.3
A very crude estimate of the mass inflow rate may be obtained from an isothermal Bondi flow. This of course assumes the matter is a fluid and ignores rotation but does give some idea of the rate at which matter comes under the gravitational influence of the black hole. This simple rate is
where M, cs, and ρ are the black hole mass |$M_9=10^9\hbox{$\rm \, {\rm M}_{\odot }$}$|, the speed of sound (or random motions), and density of the surrounding gas. Taking |$c_\mathrm{ s}=300\hbox{${\rm \, km}{\rm \, s}^{-1}\, $}$| and ρ equal to the mass density if the medium has a mass in units of |$10^9\hbox{$\rm \, {\rm M}_{\odot }$}$| per sphere of radius |$1{\rm \, kpc}$|, we obtain |$\dot{M}\approx 4 M_9^2\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| and an accretion radius of ∼50 pc.
Hopkins & Quataert (2011) give an analytical estimate of the accretion rate from gravitational torques which agrees with their numerical simulations and find |$1\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| in the middle of the range. The predicted rate has a weak black hole mass dependence of M1/6. Using observations and analytical work, Genzel et al. (2023) show that such torques operating in disc galaxies at z ∼ 2 lead to large-scale inflow on about 10 dynamical times.
In the case of the elliptical and brightest group galaxies studied here, the accretion rate could exceed the HCF mass cooling rate, which would then become the determining rate. We conclude that rates of a few |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| may be possible. We are of course assuming a high efficiency with which the cold matter is swallowed by the black hole.
The possibility thus emerges that the mass of black holes in low redshift elliptical galaxies is increasing due to inflow from HCF at a rate of several |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|. Angular momentum transfer is due to gravitational torques. The black hole mass can thus increase by up to |$\sim 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| since z = 1 and possibly even more for the most massive objects in the BCGs of the most massive clusters. Such objects need not have a luminous AGN, although an ADAF due to a weak gaseous inflow may persist and power jets thus a radio source in these objects. They would be the highest accretion rate black holes in the low-redshift Universe. If the accretion rate continues for several Gyr then this would lead to the most massive black holes appearing now. Examples in our sample include NGC1600 and the central galaxy of A85, Holm 15A, with |$1.7\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| and |$4\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| black holes, respectively. Other BCGs have very high mass black hole including that of the BCG of A1201 for which a gravitationally lensed arc reveals a central mass of |$3\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| (Nightingale et al. 2023).
4.2 Relevance to galaxy formation
The standard model of galaxy formation involves gas falling into dark matter haloes and heated by shocks and compression, gas that can cool quickly (on a dynamical time or less) then leads to star formation enriching the gas with metals and dust and the core of a new galaxy (White & Rees 1978). Further accretion, mergers and feedback later build the outer galaxy. The gas which has a longer cooling time than the local dynamical time, but shorter that the age of the Universe at the time, can form a cooling flow (Nulsen & Fabian 1995). If the conditions such as metal and dust enrichment and especially the pressure are high (|$nT\gt 10^6\hbox{${\rm \, cm}^{-3}{\rm \, K}\, $}$|) then they may resemble the nearby hidden cooling flows discussed here. If at high redshift, the higher temperature of the cosmic microwave background will in turn require a higher Jeans mass. If cloud collapse does lead to large populations of low-mass stars and brown dwarfs then early supermassive black holes can grow by swallowing such fragments whole, independent of the Eddington limit.
4.3 Red Nuggets and MRK 1216
A population of compact ETGs have been identified at redshifts of 2 and above which may be examples of galaxies that did not progress beyond the early core formation galaxy stage (Daddi et al. 2005). These are known as ‘Red Nuggets’ (Damjanov et al. 2009) and have stellar masses of |$1\!-\!2\times 10^{11}\hbox{$\rm \, {\rm M}_{\odot }$}$| and effective radii of only 1–2 kpc. Later some examples have been identified at low redshifts, e.g. NGC1277 (Trujillo et al. 2014) a galaxy unable to grow larger by mergers, or by accretion of cold gas, since it lies in the core of the rich Perseus Cluster. More recently, further examples have been found (Ferré-Mateu et al. 2015) including the isolated rotating ETG MRK 1216 (Ferré-Mateu et al. 2017) that lies at a distance of 94 Mpc and hosts a black hole of mass |$4.5\times 10^9\hbox{$\rm \, {\rm M}_{\odot }$}$| (Walsh et al. 2017).
Werner et al. (2018) noted that MRK 1216 might lie in a halo of mass up to |$10^{13}\hbox{$\rm \, {\rm M}_{\odot }$}$| and so have an X-ray halo. They indeed found extended thermal emission with an X-ray luminosity |$L_{\rm X}=7\times 10^{41}\hbox{${\rm \, erg}{\rm \, s}^{-1}\, $}.$| Buote & Barth (2019) found that its dark matter halo has a high concentration, implying early formation. We have included MRK1216 in our sample and find a significant HCF of |$9.7\pm 2.7 \hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|, larger than the rate of typical ETGs. Ferré-Mateu et al. (2017) show it to have a very bottom-heavy initial mass function (IMF) that is consistent with a significant accumulation of low mass stars and brown dwarfs.
MRK1216 could provide the nearest link between low- and high-redshift HCF and clearly merits deeper study.
4.4 Observational possibilities
Further observations at the time of writing are limited. Hopefully, XRISM will be launched soon and provide new high-resolution, non-dispersive, X-ray spectra of the inner regions of clusters, groups, and ETGs. Its field of view is larger than that of the RGS so can see how any HCF region matches into the rest of the cluster. High spatial resolution X-ray studies await next-generation telescopes such as AXIS. As well as resolving the expected irregular appearance of HCF due to the absorption, it will be particularly helpful for examining the immediate surroundings of the central black hole. The X-IFU of Athena will spectroscopically map HCFs in great detail, as will the Light Element Emission Mapper Probe. JWST may open up the inner regions in the near-IR. Since most of the flow of cooled gas takes place at very low temperatures below |$10{\rm \, K}$|, the bulk of the flow will be inaccessible, except to absorption measurements.
ALMA has opened up molecular absorption studies of cool BCGs using the central radio source as a backlight (David et al. 2014; Tremblay et al. 2016a; Rose et al. 2019, 2023). Four objects in the last study show molecular gas moving towards the central source at |$200\!-\!300\hbox{${\rm \, km}{\rm \, s}^{-1}\, $}$|, each plausibly forming part of an inward cold accretion flow.
5 CONCLUSION
We find that significant cooling flows, closely linked with cold absorbing gas, are common in the brightest galaxies of cool core clusters and groups as well as large elliptical galaxies. The mass cooling rates range from 1 to over |$1000\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|. In most cases they are reduced by AGN feedback to a factor 2–3 times lower than the simple cooling rates derived from X-ray imaging. The gas in the central hidden/absorbed part can cool to below 10 K, collapsing and fragmenting into low mass stars, brown dwarfs, etc., most of which are dragged outward by the bubbling and cold front processes. We speculate that some matter within the inner tens pc may fall into the black hole, with a rate of a few |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| being plausible. Such accretion emits little radiation, although it is likely that some thin plasma is present, possibly in the form of a low-luminosity ADAF, to power the jets usually seen in radio images. If cooled collapsed matter does fall in, then the mass accretion rate can be among the highest in the low redshift Universe.
ACKNOWLEDGEMENTS
BRM acknowledges the Natural Sciences and Engineering Research Council for their support. We thank the referee for a prompt report.
DATA AVAILABILITY
All data used here are available from ESA’s XMM–Newton Science Archive.
Footnotes
We do not repeat here the history of absorption studies in cooling flows, which is discussed in HCFI and HCFII.
A spiral feature is seen at the centre of the Centaurus cluster, see fig. 6 in HCFII.
King (2016) discusses an upper mass limit for a black hole to have a luminous accretion disc of |$\approx 5\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$|. Above that limit, the innermost parts of a luminous turbulent disc would be gravitationally unstable to fragmentation, so preventing the existence of any luminous gaseous disc. In our case we consider that the infalling matter has already fragmented and collapsed at larger radii, irrespective of the black hole mass.
References
APPENDIX A: NOTES ON OBJECTS
A1 2A0335-096
The Chandra X-ray data (Sanders, Fabian & Taylor 2009) of 2A0335 show that the central arcmin consists of six or so blobs of emission. Some are softer than the others. Cavities are seen and in the optical a bright H α nebula is found. ALMA CO data (Vantyghem et al. 2016) of the core reveal |$10^9\hbox{$\rm \, {\rm M}_{\odot }$}$| of cold molecular gas. The paper also discusses the patchy dust extinction reported earlier by Donahue et al. (2007). We use the column density of Sanders et al. (2009) in the spectral fitting of this object.
A2 A85
Holm 15 A, the BCG of A85, hosts one of the most massive black holes known, at |$4\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| (Mehrgan et al. 2019). The cluster is merging with a subcluster seen to the South. Cold fronts and a cavity are seen in Chandra X-ray data presented by (Ichinohe et al. 2015).
A3 A496
The RGS data have been studied by Tamura et al. (2001) in which no evidence for X-ray emitting gas cooler than |$1{\rm \, keV}$| was reported. (No additional absorption was considered.) Chandra X-ray images showing cold fronts have been presented by Dupke & White (2003). Optical emission lines spectra have been studied by Crawford & Fabian (1992).
A4 A2597
A2597 has consistently shown evidence for residual cooling, starting with O vi ultraviolet emission detected with FUSE (Oegerle et al. 2001), together with C iii]. The O vi line originates from thermal gas at around |$7\times 10^5{\rm \, K}$| and could be due to a cooling flow of |$20\pm 15\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$|, or |$75\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| when corrected for dust extinction.
Morris & Fabian (2005) analysed XMM data, including the RGS, and found a possible cooling flow rate of |$\sim 45\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| (corrected to |$H_0=70\hbox{${\rm \, km}\, {\rm \, s}^{-1}{\rm \, Mpc}^{-1}\, $}$|). McNamara et al. (2001) found ghost bubbles in the Chandra X-ray images and Tremblay et al. (2012a, b, 2016a, b) found evidence for some residual cooling as well as CO absorption from ALMA data. They argue for a multiphase outflow or fountain in the cluster core.
A5 A2199
Nulsen et al. (2013) presented the Chandra data for A2199. It shows signs of sloshing and has a complex radio source 3C338.
A6 M87
M87 is the second most luminous galaxy but is situated at the centre of the nearest cluster to us, the Virgo cluster. Its central black hole has a mass of |$6.5\times 10^9\hbox{$\rm \, {\rm M}_{\odot }$}$| as measured by the Event Horizon Telescope (Event Horizon Telescope Collaboration 2019). Werner et al. (2013) studied the Chandra emission in the core of M87, showing that it is highly multiphase and fitting an intrinsic absorption model. Temi, Brighenti & Mathews (2007) report the Spitzer far-infrared (FIR) flux. A power-law component with photon index of 2 was included in the fit to represent the nucleus and jet. Oldham & Auger (2018) found that M87 has a bottom-heavy IMF.
A7 NGC1399
This is the BCG of the Fornax cluster that is the next most distant cluster after Virgo. Su et al. (2019) present analysis of the Chandra X-ray data, revealing a pair of X-ray cavities coincident with radio lobes. Optical spectroscopy studied by Vaughan et al. (2018) shows that IMF in the central region is bottom heavy. They also comment on a filament of ionized gas.
A8 NGC720
NGC720 is an isolated nearby elliptical galaxy. Its Chandra X-ray data have been studied by Buote et al. (2002).
A9 NGC1550
NGC1550 is the brightest galaxy in an X-ray bright group. It has been studied recently with Chandra by Kolokythas et al. (2020). It shows signs of activity near its core due to interactions with its radio source.
A10 NGC1600
NGC1600 is the brightest galaxy in a small group. It hosts a very massive black hole of |$1.6\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$|, as determined by stellar velocity measurements by Thomas et al. (2016). Runge & Walker (2021) present the Chandra data on NGC1600 at the Bondi radius. They find evidence for a multiphase (2-temperature) X-ray emitting gas. H α emission is reported by Trinchieri & di Serego Alighieri (1991)
A11 NGC3091
NGC3091 is a bright elliptical galaxy in the relatively isolated compact group, Hickson 42. It has been selected for not being otherwise special. Our results show that it hosts a modest HCF.
A12 NGC5813
NGC5813 is the brightest galaxy in a subgroup of the larger NGC5846 group. They constitute the third nearest massive cluster/group assembly after the Virgo and Fornax clusters. The Chandra X-ray data of NGC5813 reveal a series of bubbles extending from the BGG (Randall et al. 2015). ALMA CO and other data are reported by Temi et al. (2018).
A13 NGC5846
Chandra X-ray data on NGC5846 have been studied by Machacek et al. (2011). Further X-ray analysis and multiband data on NGC5846 and many other elliptical galaxies are given by Lakhchaura et al. (2018). Radio images of this and some of the above galaxies can be found in Grossová et al. (2022).
A14 MRK1216
The Red Nugget MRK1216 is discussed in Section 4.3.
A15 ZW3146
ZW3146 is a massive X-ray luminous cluster at redshift 0.29. X-ray and ALMA CO and other data are discussed by Vantyghem et al. (2021). The centroid of the |$5\times 10^{10}\hbox{$\rm \, {\rm M}_{\odot }$}$| of molecular gas is offset by 2.6 kpc from the central AGN. The FIR detection is presented by Egami et al. (2006). MUSTANG-2 SZ results are reported by Romero et al. (2020).
APPENDIX B: SPECTRA OF OBJECTS

Top to Bottom: RGS spectrum of 2A0335+096 with HCF component shown in red and mkcflow component dotted, mass cooling rate in |$\hbox{$\hbox{$\rm \, {\rm M}_{\odot }$}{\rm \, yr}^{-1}\, $}$| versus total column density in units of |$10^{22}\hbox{${\rm \, cm}^2\, $}$|, mass cooling rate versus covering fraction of the HCF component. Contours at 68 per cent (red), 90 per cent (green), and 99 per cent (blue).













