Abstract
In the centres of many galaxy clusters, the hot (approximately 107 kelvin) intracluster medium can become dense enough that it should cool on short timescales1,2. However, the low measured star formation rates in massive central galaxies3,4,5,6 and the absence of soft X-ray lines from the cooling gas7,8,9 suggest that most of this gas never cools. This is known as the cooling flow problem. The latest observations suggest that black hole jets are maintaining the vast majority of gas at high temperatures10,11,12,13,14,15,16. A cooling flow has yet to be fully mapped through all the gas phases in any galaxy cluster. Here we present observations of the Phoenix cluster17 using the James Webb Space Telescope to map the [Ne vi] λ 7.652-μm emission line, enabling us to probe the gas at 105.5 kelvin on large scales. These data show extended [Ne vi] emission that is cospatial with the cooling peak in the intracluster medium, the coolest gas phases and the sites of active star formation. Taken together, these imply a recent episode of rapid cooling, causing a short-lived spike in the cooling rate, which we estimate to be 5,000–23,000 solar masses per year. These data provide a large-scale map of gas at temperatures between 105 kelvin and 106 kelvin in a cluster core, and highlight the critical role that black hole feedback has in not only regulating cooling but also promoting it18.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
27,99 € / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
199,00 € per year
only 3,90 € per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout




Similar content being viewed by others
Data availability
Data for JWST programme ID 2439 are publicly available through the Mikulski Archive for Space Telescopes of the Space Telescope Science Institute. Supplementary data used in our analysis from HST are also available through the Mikulski Archive for Space Telescopes. Chandra data are available from the Chandra Data Archive (https://cxc.harvard.edu/cda/).
Code availability
JWST, HST and Chandra data were reduced using the publicly available reduction pipeline codes provided by the Space Telescope Science Institute and Chandra X-ray Center. The Likelihood Optimization of gas Kinematics code used in our analysis is available through GitHub at https://github.com/Michael-Reefe/Loki.jl, and Cloudy is available at https://nublado.org. We also provide our customized driver scripts for data reduction and analysis through GitHub at https://github.com/Michael-Reefe/Reefe2024_code_supplements.
References
Fabian, A. C., Nulsen, P. E. J. & Canizares, C. R. Cooling flows in clusters of galaxies. Nature 310, 733–740 (1984).
Fabian, A. C. Cooling flows in clusters of galaxies. Annu. Rev. Astron. Astrophys. 32, 277–318 (1994).
McNamara, B. R. & O’Connell, R. W. Star formation in cooling flows in clusters of galaxies. Astron. J. 98, 2018–2043 (1989).
Allen, S. W. Starbursts in cooling flows: blue continua and emission-line nebulae in central cluster galaxies. Mon. Not. R. Astron. Soc. 276, 947–960 (1995).
Hicks, A. K. & Mushotzky, R. Star formation rates in cooling flow clusters: a UV pilot study with archival XMM-Newton optical monitor data. Astrophys. J. Lett. 635, L9–L12 (2005).
McDonald, M., Gaspari, M., McNamara, B. R. & Tremblay, G. R. Revisiting the cooling flow problem in galaxies, groups, and clusters of galaxies. Astrophys. J. 858, 45 (2018).
Canizares, C. R., Markert, T. H. & Donahue, M. E. in Cooling Flows in Clusters and Galaxies (ed. Fabian, A. C.) 63 (1988).
David, L. P. et al. A high-resolution study of the Hydra A cluster with Chandra: comparison of the core mass distribution with theoretical predictions and evidence for feedback in the cooling flow. Astrophys. J. 557, 546–559 (2001).
Peterson, J. R. et al. High-resolution X-ray spectroscopic constraints on cooling-flow models for clusters of galaxies. Astrophys. J. 590, 207–224 (2003).
McNamara, B. R. & Nulsen, P. E. J. Heating hot atmospheres with active galactic nuclei. Annu. Rev. Astron. Astrophys. 45, 117–175 (2007).
Fabian, A. C. Observational evidence of active galactic nuclei feedback. Annu. Rev. Astron. Astrophys. 50, 455–489 (2012).
Hlavacek-Larrondo, J. et al. X-ray cavities in a sample of 83 SPT-selected clusters of galaxies: tracing the evolution of AGN feedback in clusters of galaxies out to z = 1.2. Astrophys. J. 805, 35 (2015).
Gaspari, M., Melioli, C., Brighenti, F. & D’Ercole, A. The dance of heating and cooling in galaxy clusters: three-dimensional simulations of self-regulated active galactic nuclei outflows. Mon. Not. R. Astron. Soc. 411, 349–372 (2011).
Prasad, D., Sharma, P. & Babul, A. Cool core cycles: cold gas and AGN jet feedback in cluster cores. Astrophys. J. 811, 108 (2015).
Li, Y., Ruszkowski, M. & Bryan, G. L. AGN heating in simulated cool-core clusters. Astrophys. J. 847, 106 (2017).
Yang, H. Y. K., Gaspari, M. & Marlow, C. The impact of radio AGN bubble composition on the dynamics and thermal balance of the intracluster medium. Astrophys. J. 871, 6 (2019).
McDonald, M. et al. A massive, cooling-flow-induced starburst in the core of a luminous cluster of galaxies. Nature 488, 349–352 (2012).
McNamara, B. R. et al. A mechanism for stimulating AGN feedback by lifting gas in massive galaxies. Astrophys. J. 830, 79 (2016).
Revaz, Y., Combes, F. & Salomé, P. Formation of cold filaments in cooling flow clusters. Astron. Astrophys. 477, L33–L36 (2008).
McDonald, M. et al. Anatomy of a cooling flow: the feedback response to pure cooling in the core of the Phoenix cluster. Astrophys. J. 885, 63 (2019).
Russell, H. R. et al. Alma observations of massive molecular gas filaments encasing radio bubbles in the Phoenix cluster. Astrophys. J. 836, 130 (2017).
Gaspari, M. et al. Shaken snow globes: kinematic tracers of the multiphase condensation cascade in massive galaxies, groups, and clusters. Astrophys. J. 854, 167 (2018).
McDonald, M. et al. Deep Chandra, HST-COS, and Megacam observations of the Phoenix cluster: extreme star formation and AGN feedback on hundred kiloparsec scales. Astrophys. J. 811, 111 (2015).
Pinto, C. et al. AGN feedback in the Phoenix cluster. Mon. Not. R. Astron. Soc. 480, 4113–4123 (2018).
Oegerle, W. R. et al. FUSE observations of cooling-flow gas in the galaxy clusters A1795 and A2597. Astrophys. J. 560, 187–193 (2001).
Bregman, J. N., Fabian, A. C., Miller, E. D. & Irwin, J. A. On VI observations of galaxy clusters: evidence for modest cooling flows. Astrophys. J. 642, 746–751 (2006).
Fabian, A. C. et al. Hidden cooling flows in clusters of galaxies. Mon. Not. R. Astron. Soc. 515, 3336–3345 (2022).
Sparks, W. B. et al. Hundred thousand degree gas in the Virgo cluster of galaxies. Astrophys. J. Lett. 750, L5 (2012).
Anderson, M. E. & Sunyaev, R. Searching for FUV line emission from 107 K gas in massive elliptical galaxies and galaxy clusters as a tracer of turbulent velocities. Mon. Not. R. Astron. Soc. 459, 2806–2821 (2016).
Canning, R. E. A. et al. Detection of optical coronal emission from 106-K gas in the core of the Centaurus cluster. Mon. Not. R. Astron. Soc. 411, 411–421 (2011).
Chatzikos, M. et al. Implications of coronal line emission in NGC 4696*. Mon. Not. R. Astron. Soc. 446, 1234–1244 (2015).
McDonald, M. et al. The state of the warm and cold gas in the extreme starburst at the core of the Phoenix galaxy cluster (SPT-CLJ2344-4243). Astrophys. J. 784, 18 (2014).
Pope, E. C. D., Babul, A., Pavlovski, G., Bower, R. G. & Dotter, A. Mass transport by buoyant bubbles in galaxy clusters. Mon. Not. R. Astron. Soc. 406, 2023–2037 (2010).
Gaspari, M., Ruszkowski, M. & Sharma, P. Cause and effect of feedback: multiphase gas in cluster cores heated by AGN jets. Astrophys. J. 746, 94 (2012).
Wittor, D. & Gaspari, M. Dissecting the turbulent weather driven by mechanical AGN feedback. Mon. Not. R. Astron. Soc. 498, 4983–5002 (2020).
Ferland, G. J. et al. CLOUDY 90: numerical simulation of plasmas and their spectra. Publ. Astron. Soc. Pac. 110, 761–778 (1998).
Chatzikos, M. et al. The 2023 release of Cloudy. Rev. Mex. Astron. Astrofis. 59, 327–343 (2023).
Saccheo, I. et al. The WISSH quasars project. XI. The mean spectral energy distribution and bolometric corrections of the most luminous quasars. Astron. Astrophys. 671, A34 (2023).
Indriolo, N., Geballe, T. R., Oka, T. & McCall, B. J. H+3 in diffuse interstellar clouds: a tracer for the cosmic-ray ionization rate. Astrophys. J. 671, 1736–1747 (2007).
Gaspari, M. Shaping the X-ray spectrum of galaxy clusters with AGN feedback and turbulence. Mon. Not. R. Astron. Soc. 451, L60–L64 (2015).
Jenkins, E. B. A unified representation of gas-phase element depletions in the interstellar medium. Astrophys. J. 700, 1299–1348 (2009).
Gaspari, M., Temi, P. & Brighenti, F. Raining on black holes and massive galaxies: the top-down multiphase condensation model. Mon. Not. R. Astron. Soc. 466, 677–704 (2017).
Prasad, D., Sharma, P., Babul, A., Voit, G. M. & O’Shea, B. W. Cool-core cycles and Phoenix. Mon. Not. R. Astron. Soc. 495, 594–599 (2020).
Fabian, A. C. et al. The relationship between the optical Hα filaments and the X-ray emission in the core of the Perseus cluster. Mon. Not. R. Astron. Soc. 344, L48–L52 (2003).
Gendron-Marsolais, M. et al. Revealing the velocity structure of the filamentary nebula in NGC 1275 in its entirety. Mon. Not. R. Astron. Soc. 479, L28–L33 (2018).
McNamara, B. R. & Nulsen, P. E. J. Mechanical feedback from active galactic nuclei in galaxies, groups and clusters. New J. Phys. 14, 055023 (2012).
Fabian, A. C. et al. Hidden cooling flows in clusters of galaxies II: a wider sample. Mon. Not. R. Astron. Soc. 521, 1794–1807 (2023).
Spilker, J. S. et al. Spatial variations in aromatic hydrocarbon emission in a dust-rich galaxy. Nature 618, 708–711 (2023).
Bezanson, J., Edelman, A., Karpinski, S. & Shah, V. B. Julia: a fresh approach to numerical computing. SIAM Rev. 59, 65–98 (2017).
Smith, J. D. T. et al. The mid-infrared spectrum of star-forming galaxies: global properties of polycyclic aromatic hydrocarbon emission. Astrophys. J. 656, 770–791 (2007).
Fitzpatrick, E. L. & Massa, D. An analysis of the shapes of ultraviolet extinction curves. I. The 2,175 angstrom bump. Astrophys. J. 307, 286 (1986).
Donnan, F. R. et al. The obscured nucleus and shocked environment of VV 114E revealed by JWST/MIRI spectroscopy. Mon. Not. R. Astron. Soc. 519, 3691–3705 (2023).
Labiano, A. et al. Wavelength calibration and resolving power of the JWST MIRI medium resolution spectrometer. Astron. Astrophys. 656, A57 (2021).
Tsuchikawa, T. et al. A systematic study of silicate absorption features in heavily obscured AGNs observed by Spitzer/IRS. Astron. Astrophys. 651, A117 (2021).
Veilleux, S. et al. Spitzer Quasar and Ulirg Evolution Study (QUEST). IV. Comparison of 1 Jy ultraluminous infrared galaxies with Palomar-Green quasars. Astrophys. J. Supp. 182, 628–666 (2009).
Markwardt, C. B. Non-linear least-squares fitting in IDL with MPFIT. in Astronomical Data Analysis Software and Systems XVIII (eds Bohlender, D. A. et al.) 251 (2009).
Mogensen, P. K. & Riseth, A. N. Optim: a mathematical optimization package for Julia. J. Open Source Softw. 3, 615 (2018).
Del Zanna, G., Dere, K. P., Young, P. R. & Landi, E. CHIANTI—an atomic database for emission lines. XVI. version 10, further extensions. Astrophys. J. 909, 38 (2021).
Fielding, D. B. & Bryan, G. L. The structure of multiphase galactic winds. Astrophys. J. 924, 82 (2022).
Meena, B. et al. Investigating the narrow-line region dynamics in nearby active galaxies. Astrophys. J. 943, 98 (2023).
Acknowledgements
This work is based on observations with the NASA/ESA/CSA JWST obtained from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract number NAS 5-03127. Support for programme number JWST-GO-02439.001-A was provided through a grant from the Space Telescope Science Institute under NASA contract number NAS 5-03127. This work is also based in part on observations made with the NASA/ESA HST obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy under NASA contract number NAS 5-26555. These observations are associated with programme number GO15315. Support for this work was also provided by NASA through Chandra award number GO7-18124 issued by Chandra, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of NASA under contract number NAS8-03060. M.R. acknowledges support from the National Science Foundation Graduate Research Fellowship under grant number 2141064. M.G. acknowledges support from the ERC Consolidator Grant BlackHoleWeather (101086804). M. Chatzikos acknowledges support from NSF (1910687) and NASA (19-ATP19-0188 and 22-ADAP22-0139). H.R. acknowledges an Anne McLaren Fellowship from the University of Nottingham. We thank the members of the MIRI/MRS instrument team, particularly D.L., for providing advice and guidance in reducing and cleaning the MIRI/MRS data.
Author information
Authors and Affiliations
Contributions
M.R. reduced the MIRI/MRS data; wrote and ran the spectral fitting and PSF decomposition code; performed the main analysis of the ionization, kinematics and cooling; and wrote the paper. M.M. provided substantial guidance in the main analyses, interpretation of the results and writing the paper; he also wrote the original proposal for the MIRI/MRS data. M. Chatzikos performed the Cloudy simulations and wrote the description of these processes in Methods. J.S. performed a sanity check on our data reduction and spectral modelling by running an independent pipeline (Q3Dfit). K.S. formulated strong lensing models to estimate the gravitational potential, which we used in our kinematic modelling. R.M., S.V., S.W.A., M.B., M. Calzadilla, R.C., B.F., M.G., J.H.-L., B.M., H.R. and T.S. provided feedback on the paper drafts.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 The field of view of the MIRI/MRS instrument relative to the size of the central galaxy in the Phoenix Cluster.
The field of view of each of the four MIRI channels shown over a 3-color image of the core of the Phoenix cluster using data from HST/ACS in the F475W, F775W, and F850LP filters. The dust-obscured AGN lies at the center of the brightest cluster galaxy. The angular and physical scales are annotated in the bottom-left.
Extended Data Fig. 2 The full mid-infrared spectrum of the central galaxy in the Phoenix Cluster.
A spectrum covering MIRI channels 2–4 (excluding the long wavelength end of 4C), and integrated over the full channel 2 FOV. The top panel shows the data in black and the model in orange, while the bottom panel shows the residuals. The model is also decomposed into components: the gray solid lines are the thermal dust continua, the green line is the QSO PSF model, the blue line is the PAH emission, the purple lines are the emission lines, and the dotted gray line shows the extinction profile. The sum of the thermal dust continua and QSO PSF is also shown by the thick gray line. The emission lines are labeled at the top of the plot, with vertical dashed lines showing their locations in the spectrum. The boundaries between channels and bands are also shown at the top of the plot with blue arrows. The translucent orange region shows the range of models produced in the different bootstrapping iterations.
Extended Data Fig. 3 A zoomed-in emission line spectrum of the central galaxy in the Phoenix Cluster.
A subset of the MIRI channel 3 spectrum zoomed in on the [Ne vi] line. This spectrum includes only a single spaxel north of the nucleus. The formatting of everything is identical to Extended Data Fig. 2, except that there is no translucent orange region since this fit has not been bootstrapped. The two distinct kinematic components of the emission lines can be seen in purple.
Extended Data Fig. 4 A series of flux maps for [Ne vi] over the MIRI channel 3 FOV showing the decomposition into various spatial and kinematic components.
a-c show the QSO and host galaxy components of the flux. a, the total observed flux, b, the flux from the QSO that has been dispersed according to the PSF model, and c, the QSO-subtracted flux that is attributed to the host galaxy. d-f show the further decomposition of the host galaxy flux into 2 distinct kinematic components, now with an S/N cut such that only spaxels with a detection of S/N ≥ 3 are shown. d, the combined flux from both kinematic components, making it identical to c except for the SNR cut. e and f, the fluxes from each individual kinematic component, sorted in order of decreasing flux. The color scales are shown on the right of each row and are the same for each panel in that row. The physical scale in kpc and angular scale in arcsec are annotated in the bottom left of each panel, and the FWHM of the PSF at the wavelength of [Ne vi] is shown in the bottom right of each panel.
Extended Data Fig. 5 A series of flux maps for [Ne ii] over the MIRI channel 4 FOV showing the decomposition into various spatial and kinematic components.
Same as Extended Data Fig. 4, but for [Ne ii] λ12.813 μm.
Extended Data Fig. 6 A series of maps for [Ne v] over the MIRI channel 4 FOV showing the decomposition into various spatial and kinematic components.
Same as Extended Data Fig. 4, but for [Ne v] λ14.322 μm.
Extended Data Fig. 7 A series of kinematic maps for [Ne vi] over the MIRI channel 3 FOV showing the decomposition into various kinematic components.
a-c show the line velocity widths. a, the W80 of the full line profile, b and c, the FWHMs of the individual kinematic components. d-f show the LOS velocities. d, the median LOS velocity of the full line profile, e and f, the peak velocities of the individual kinematic components. All maps have an S/N ≥ 3 cut. The color scales are shown on the right of each row and are the same for each panel in that row. The physical scale in kpc and angular scale in arcsec are annotated in the bottom left of each panel, and the FWHM of the PSF at the wavelength of [Ne vi] is shown in the bottom right of each panel.
Extended Data Fig. 8 A series of kinematic maps for [Ne ii] over the MIRI channel 4 FOV showing the decomposition into various kinematic components.
Same as Extended Data Fig. 7, but for [Ne ii] λ12.813 μm.
Extended Data Fig. 9 A series of kinematic maps for [Ne v] over the MIRI channel 4 FOV showing the decomposition into various kinematic components.
Same as Extended Data Fig. 7, but for [Ne v] λ14.322 μm.
Supplementary information
Supplementary Information
Supplementary Methods, Figs. 1–3, Tables 1 and 2, and References.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Reefe, M., McDonald, M., Chatzikos, M. et al. Directly imaging the cooling flow in the Phoenix cluster. Nature 638, 360–364 (2025). https://doi.org/10.1038/s41586-024-08369-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-024-08369-x