Cooling Flows and Cool-Core Clusters

Galaxy clusters are the largest gravitationally bound structures in the Universe, with masses of $10^{14}$ – $10^{15}\,M_\odot$ set mostly by dark matter. But their baryons are not chiefly stars — they are a diffuse, million-degree atmosphere, the intracluster medium (ICM), radiating X-rays at $T\sim10^7$ – $10^8\,\mathrm{K}$ ( $k_BT\sim1$ – $10\,\mathrm{keV}$ ). Because thermal bremsstrahlung radiates at a rate $\propto n^2\sqrt{T}$ , the densest gas — packed into the core — glows fastest and can radiate its thermal energy away in less than a Hubble time. Hydrostatic reasoning then makes an almost inevitable prediction: as the core loses pressure support, the weight of the overlying atmosphere squeezes it gently inward. This slow, subsonic settling is the cooling flow. The classical model of the 1970s and 80s took the picture to its logical end: gas cooling isobarically to low temperature and condensing at hundreds of solar masses per year onto the central galaxy. Cool-core clusters do show the tell-tale signs — peaked X-ray brightness, short central cooling times — so the story looked secure. Yet it made one sharp, falsifiable prediction: a full cascade of emission lines from gas at every temperature down toward zero. When XMM-Newton's grating spectrometers resolved those lines around 2001, the cascade was missing: gas essentially stopped cooling below about a third of the ambient temperature, and the inferred condensation rates were a tenth of the classical value. This is the cooling-flow problem, and its resolution — finely-tuned mechanical feedback from the central supermassive black hole — has made these systems the definitive laboratory for how black holes regulate the largest galaxies.

Key takeaways

The hot ICM in cluster cores radiates X-rays by thermal bremsstrahlung; where the cooling time drops below the cluster age, gas loses pressure support and settles inward — the cooling flow.
Naive theory predicts a large mass deposition rate, $\dot M\sim100$ – $1000\,M_\odot\,\mathrm{yr}^{-1}$ , feeding runaway condensation onto the central galaxy.
High-resolution X-ray spectroscopy (XMM-Newton RGS, Chandra) finds almost no gas cooling below $\sim T_{\rm vir}/3$ — the cooling-flow problem: the missing sink of cooled gas.
The resolution is feedback: radio jets inflate X-ray cavities that reheat the core, a self-regulated thermostat nearly offsetting radiative losses.
Clusters split into cool-core (peaked density, low central entropy) and non-cool-core systems — a fossil record of their merger histories.
Key idea: Wherever the radiative cooling time of the ICM falls below the cluster age, gas should condense and flow inward; the observed near-absence of such cooling implies a heating source — AGN mechanical feedback — that closely balances radiative losses in a self-regulated thermostat.

Where this lives today

Cooling flows sit at the intersection of X-ray astronomy, black-hole accretion, and galaxy formation. Chandra's arcsecond imaging turned cavities into calorimeters: pressure and volume give the $4PV$ enthalpy, and dividing by a buoyancy time gives a jet power to compare with the cooling luminosity — done now for dozens of clusters. The microcalorimeter on Hitomi measured the turbulent velocity of the Perseus core in 2016, finding it surprisingly quiescent ( $\sim160\,\mathrm{km\,s^{-1}}$ ); its successor XRISM (2023) is now mapping ICM velocities widely. Simulations of chaotic cold accretion test whether the $t_{\rm cool}/t_{\rm ff}\lesssim10$ threshold controls the black hole's fuelling, while ALMA reveals the cold molecular filaments — the "rain" — condensing in cool cores. The physics reaches into cosmology too: the cool-core boost inflates a cluster's X-ray luminosity and must be modelled out of the luminosity-mass relation used to weigh clusters in surveys like eROSITA. How the energy is distributed so gently, and how the dichotomy is set, remain open.

The physical picture

Picture the ICM as a vast atmosphere, held up by its own thermal pressure against the cluster's gravity. Radiation is a slow leak in that pressure: every escaping X-ray photon carries off thermal energy, and the densest gas leaks fastest because the emissivity scales as density squared. So the core cools first; its pressure sags, and the weight of the gas above presses it inward. This is not a collapse — the cooling time far exceeds the sound-crossing time, so pressure stays near hydrostatic and the gas merely settles, a subsonic drift of tens of kilometres per second. That slow inward creep is the cooling flow, whose end-point, absent heating, is a torrent of cold gas raining onto the central galaxy. The first figure makes the geometry quantitative: it plots the cooling time as a function of radius and locates the cooling radius where that time equals the cluster's age.

Figure

The cooling time profile of a cool-core cluster: t_cool(r) rises outward and crosses the cluster age t_age ≈ 7.7 Gyr at the cooling radius r_cool ≈ 70 kpc. Inside r_cool the gas radiates within a Hubble time and settles inward at v_in ≈ 25 km/s, deeply subsonic, reaching under 1 Gyr in the core

The naive expectation of that inflow is catastrophe: denser gas cools faster, condenses, and the brightest cluster galaxy should be a blazing star factory drowning in molecular gas. It is not — star-formation rates and cold-gas masses fall an order of magnitude short. Something reheats the core just enough to nearly, but not quite, cancel the cooling. The culprit is the central black hole, which senses the gas through accretion and answers with jets. The jets do not explode the gas; they inflate buoyant bubbles of relativistic plasma that rise like smoke rings, doing $PV$ work and dissipating their energy through weak shocks, sound waves, turbulence and mixing. The second figure shows this thermostat and its energy budget.

Figure

The AGN feedback thermostat: jets from the central black hole inflate a pair of X-ray cavities of relativistic plasma. Each carries enthalpy E = 4PV; for two 15-kpc cavities in a P ≈ 3×10⁻¹⁰ erg cm⁻³ core the total is ≈ 10⁶⁰ erg, and over a buoyant rise time ≈ 10⁸ yr this gives L_cav ≈ 3×10⁴⁴ erg s⁻¹ ≈ L_cool

The energy deposited by this cycle — of order $10^{44}$ – $10^{45}\,\mathrm{erg\,s^{-1}}$ time-averaged — must track the radiated luminosity across 100-kpc scales and gigayear timescales, remarkable given that the black hole is $10^9$ times smaller than the region it heats. The cool-core versus non-cool-core split is the fossil record of this balance: a cluster's central entropy remembers whether a recent merger reheated its core, or whether it has relaxed into a peaked, short-cooling-time state held in check by gentle feedback.

The core equations

The subject rests on one fact — how fast hot, thin plasma radiates. In massive clusters that is thermal bremsstrahlung, whose volume emissivity fixes everything downstream.

\varepsilon_{\rm ff} = 1.4\times10^{-27}\,T^{1/2}\,n_e n_i\,Z^2\,\bar{g}_B\ \ \mathrm{erg\,cm^{-3}\,s^{-1}} \qquad (1)

where: $\varepsilon_{\rm ff}$ — free-free volume emissivity; $T$ — gas temperature in kelvin; $n_e,\,n_i$ — electron and ion number densities; $Z$ — ion charge; $\bar{g}_B$ — frequency-averaged Gaunt factor, of order unity. The lesson is the density-squared scaling: double the density, radiate four times faster — so the dense core dominates the emission. Dividing the gas's thermal energy content by that loss rate gives the cooling time — the field's most important quantity, since it decides where and whether a cooling flow develops.

t_{\rm cool} = \frac{\tfrac{3}{2}(n_e+n_i)k_B T}{n_e n_i \Lambda(T)} \;\approx\; \frac{3}{2}\frac{n\,k_B T}{n_e^2\,\Lambda(T)} \qquad (2)

where: $t_{\rm cool}$ — cooling time; $n=n_e+n_i$ — total particle number density; $\Lambda(T)$ — cooling function (bremsstrahlung plus line emission), units $\mathrm{erg\,cm^3\,s^{-1}}$ ; $k_B$ — Boltzmann constant. Numerator $\propto n$ , denominator $\propto n^2$ , so cooling is short exactly where the gas is dense — the core. Evaluating the bremsstrahlung form at cluster-scale fiducials turns the scaling into a usable number for the densest cores.

t_{\rm cool} \approx 8.5\times10^{10}\,\mathrm{yr}\left(\frac{n_e}{10^{-3}\,\mathrm{cm^{-3}}}\right)^{-1}\left(\frac{T}{10^{8}\,\mathrm{K}}\right)^{1/2} \qquad (3)

where: $n_e$ — electron density, reaching $\sim10^{-2}$ – $10^{-1}\,\mathrm{cm^{-3}}$ in cool-core centres; $T$ — gas temperature. At the outskirts the cooling time exceeds the age of the Universe; a hundredfold-denser core cools in under a gigayear — the origin of the cool-core phenomenon. The classical model turns a measured cooling luminosity into a predicted condensation rate: every radiated erg comes from gas losing its enthalpy as it settles isobarically.

\dot{M} = \frac{2\,\mu m_H\,L_{\rm cool}}{5\,k_B T} \qquad (4)

where: $\dot{M}$ — mass deposition rate; $L_{\rm cool}$ — luminosity radiated within the cooling region; $\mu\approx0.6$ — mean molecular weight of the ionized ICM; $m_H$ — hydrogen mass. For luminous cool cores this returns $\sim100$ – $1000\,M_\odot\,\mathrm{yr}^{-1}$ — the celebrated, badly over-stated prediction. The entropy index compares cores across clusters of different mass; it is conserved adiabatically, changing only when gas gains or loses heat.

K \equiv \frac{k_B T}{n_e^{2/3}} \qquad (5)

where: $K$ — entropy index (units $\mathrm{keV\,cm^2}$ ), with $K\propto e^{2s/3k_B}$ ; $n_e$ — electron density. A low central $K_0\lesssim30\,\mathrm{keV\,cm^2}$ flags a cool core; a high, flat central entropy flags a merger-heated cluster. On the heating side, the currency is the enthalpy of a jet-inflated cavity — the energy a buoyant relativistic-plasma bubble releases as it rises, giving up its internal energy and its inflation work.

E_{\rm cav} = \frac{\gamma}{\gamma-1}\,P\,V \;=\; 4\,P V \quad(\gamma=4/3) \qquad (6)

where: $P$ — ambient ICM pressure at the cavity; $V$ — cavity volume; $\gamma=4/3$ — adiabatic index of a relativistic plasma. The $4PV$ is $3PV$ of internal energy plus $1PV$ inflation work; divided by a buoyant rise time it gives the power $L_{\rm cav}$ that offsets $L_{\rm cool}$ . Finally, whether cold gas condenses is set not by the cooling time alone but by its ratio to the local free-fall time — the precipitation criterion tying cluster thermodynamics to the black hole's fuel supply.

\frac{t_{\rm cool}}{t_{\rm ff}} \lesssim 10 \;\Rightarrow\; \text{thermal instability, cold-gas condensation} \qquad (7)

where: $t_{\rm ff}=\sqrt{2r/g}$ — free-fall time at radius $r$ ; $t_{\rm cool}$ — local cooling time. Where the ratio dips below $\sim10$ , unstable blobs cool and rain onto the central galaxy before buoyancy mixes them away, feeding chaotic cold accretion. For context, the black hole is bounded by its Eddington luminosity, though here it operates far below it in a jet-dominated, radiatively inefficient mode.

L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T}\approx1.3\times10^{38}\left(\frac{M_{\rm BH}}{M_\odot}\right)\,\mathrm{erg\,s^{-1}} \qquad (8)

where: $M_{\rm BH}\sim10^9$ – $10^{10}\,M_\odot$ — central black hole mass; $\sigma_T$ — Thomson cross section; $m_p$ — proton mass. A $10^{10}\,M_\odot$ hole has $L_{\rm Edd}\sim10^{48}\,\mathrm{erg\,s^{-1}}$ , so jet powers of $\sim10^{44}$ – $10^{45}$ mean accretion at $\sim10^{-3}$ – $10^{-4}$ of Eddington — the low, hot-mode regime.

Derivation AThe cooling radius and why the inflow is a flow, not a collapse

The goal is to locate where cooling matters and to show the resulting inflow is slow and subsonic — a settling, not a collapse, which is what makes "cooling flow" the right name. We fix the cooling radius where the local cooling time equals the cluster age, use mass continuity to get the inflow speed, and compare it to the sound speed.

Step 1. The cooling-time profile follows from equation (3): $t_{\rm cool}(r)\propto \sqrt{T(r)}/n_e(r)$ . With the centrally peaked density of a relaxed cool core it plunges inward, crossing the cluster age $t_{\rm age}\approx7.7\,\mathrm{Gyr}$ at $r_{\rm cool}\sim50$ – $100\,\mathrm{kpc}$ — the cooling radius. Why: outside $r_{\rm cool}$ gas cannot radiate within a Hubble time, so all the action is confined to $r<r_{\rm cool}$ .

Step 2. Mass continuity for steady spherical inflow, $\dot M = 4\pi r^2 \rho\, v_{\rm in}$ , gives $v_{\rm in}\sim r_{\rm cool}/t_{\rm cool}\approx25\,\mathrm{km\,s^{-1}}$ for $r_{\rm cool}\approx70\,\mathrm{kpc}$ , $t_{\rm cool}\approx2.7\,\mathrm{Gyr}$ . Why: this turns the cooling time into a physical drift speed to weigh against the atmosphere's dynamical speed.

Step 3. Compare to $c_s=\sqrt{\gamma k_B T/\mu m_H}\approx950\,\mathrm{km\,s^{-1}}$ at $k_BT\approx3.4\,\mathrm{keV}$ : the Mach number is $\approx0.03$ , deeply subsonic, so $t_{\rm cool}\gg r/c_s$ and pressure stays hydrostatic. Why: sound re-levels pressure far faster than radiation drains it, so the gas settles gently — a flow, not a collapse.

Step 4. The classical dynamics are thus self-consistent. The error the observations exposed is thermodynamic, not dynamical: the assumption that cooling proceeds unopposed to low temperature. Why: this points the finger at a missing heat source rather than at the hydrostatic reasoning.

Go deeper · Derivation B · 1 minAGN feedback energy balance from cavity enthalpy

Derivation BAGN feedback energy balance from cavity enthalpy

The goal is to show, with numbers, that jet-inflated cavities carry enough energy to offset a cooling core's radiative losses — the quantitative case for mechanical AGN feedback. We compute a relativistic bubble's enthalpy, estimate its buoyant deposition time, and compare the resulting power to the cooling luminosity; the two match to order of magnitude.

Step 1. A cavity is relativistic plasma ( $\gamma=4/3$ ) in pressure balance with the ICM; its energy is the enthalpy $E_{\rm cav}=\gamma/(\gamma-1)\,PV=4PV$ — internal energy $3PV$ plus the $PV$ inflation work against the ambient pressure. Why: enthalpy, not just internal energy, is the right measure because the bubble both stores heat and displaces gas, and the ICM recovers both as the wake refills.

Step 2. The bubble releases this on its buoyant rise time, $t_{\rm buoy}\sim R\sqrt{SC/2gV}\sim10^7$ – $10^8\,\mathrm{yr}$ (comparable to the sound-crossing time), giving power $L_{\rm cav}=4PV/t_{\rm buoy}$ . Why: it is the power, not the total energy, that must match the continuous radiative drain.

Step 3. Numbers: two cavities of radius $R=15\,\mathrm{kpc}$ in $P\approx3\times10^{-10}\,\mathrm{erg\,cm^{-3}}$ give $V=\tfrac{4}{3}\pi R^3\approx4\times10^{68}\,\mathrm{cm^3}$ , total enthalpy $\approx10^{60}\,\mathrm{erg}$ , and over $t_{\rm buoy}\approx10^{8}\,\mathrm{yr}$ , $L_{\rm cav}\approx3\times10^{44}\,\mathrm{erg\,s^{-1}}$ . Why: this is the same order as a luminous cool core's $L_{\rm cool}\sim10^{44}\,\mathrm{erg\,s^{-1}}$ .

Step 4. Because heating must track cooling over gigayears, feedback is high-duty-cycle: cooling feeds the hole, jets fire, cooling is suppressed, accretion drops, jets fade, cooling resumes. The hole accretes at $\sim10^{-3}$ – $10^{-4}\,L_{\rm Edd}$ , channelling output into jets rather than light. Why: this closes the thermostat and explains why cool-core AGN are radio-loud but optically feeble.

Go deeper · Derivation C · 1 minThermal instability and the precipitation threshold

Derivation CThermal instability and the precipitation threshold

The goal is to derive when cold gas condenses from the hot atmosphere, linking Field's linear thermal-instability criterion to the empirical precipitation threshold that modern simulations converge on. We write the net loss function, apply the isobaric criterion, then let gravity compete through the free-fall time.

Step 1. Define the net volumetric loss $\mathcal{L}=n_e^2\Lambda(T)-\mathcal{H}$ (cooling minus heating $\mathcal{H}$ ); global balance means $\mathcal{L}=0$ on average. Why: feedback offsets cooling on average, so the question is the stability of that balance, not whether cooling occurs.

Step 2. Perturb at constant pressure: Field's criterion $(\partial\mathcal{L}/\partial T)_P<0$ means a slightly cooler blob radiates faster and runs away; in the line-cooling regime, where $\Lambda(T)$ rises as gas cools, it is readily met. Why: constant pressure is right because a condensing blob stays in pressure balance, contracting and densifying as it cools.

Step 3. Gravity opposes the runaway: a cool, dense blob sinks on the free-fall time $t_{\rm ff}=\sqrt{2r/g}$ ; if it falls and mixes before cooling, it is reabsorbed. Why: instability alone is insufficient — the blob must cool faster than buoyancy carries it into hotter gas.

Step 4. The contest is set by the ratio $t_{\rm cool}/t_{\rm ff}$ . Cold molecular gas, H $\alpha$ filaments and star formation appear where $\min(t_{\rm cool}/t_{\rm ff})\lesssim10$ — the precipitation regime that feeds chaotic cold accretion. Why: this links 100-kpc thermodynamics to sub-parsec fuelling, so the atmosphere regulates its own black hole.

Limiting cases

The behaviour of the ICM core depends on how the cooling time compares to the cluster age and on whether bremsstrahlung or line emission dominates the losses.

Case	Result	Meaning
Bremsstrahlung regime ( $k_BT\gtrsim2\,\mathrm{keV}$ )	$t_{\rm cool}\propto \sqrt{T}/n$	Hot clusters cool by free-free emission; cooling time grows with $T$ , so only the dense core cools appreciably.
Line-cooling regime ( $k_BT\lesssim2\,\mathrm{keV}$ )	$\Lambda(T)$ rises as metal lines dominate	In groups and cooler gas, line emission steepens $\Lambda(T)$ ; cooling accelerates below $\sim1\,\mathrm{keV}$ , sharpening the truncation.
Cool-core limit ( $t_{\rm cool}\ll t_{\rm Hubble}$ )	$t_{\rm cool}(r{\to}0)\lesssim1\,\mathrm{Gyr},\ K_0\lesssim30\,\mathrm{keV\,cm^2}$	Peaked density, low central entropy; a cooling flow is expected but held in check by feedback.
Non-cool-core limit	$t_{\rm cool}(r{\to}0)\gtrsim t_{\rm Hubble},\ K_0\gtrsim30\,\mathrm{keV\,cm^2}$	Flat, high core entropy — merger-disturbed; no short-cooling-time region, so no cooling flow.

Worked examples

1 · Cooling time and classical deposition rate of a cool core (fully worked)

A cool-core cluster has central electron density $n_e=0.02\,\mathrm{cm^{-3}}$ , temperature $k_BT=3.4\,\mathrm{keV}$ (so $T\approx4\times10^{7}\,\mathrm{K}$ ), and a cooling luminosity $L_{\rm cool}=3\times10^{44}\,\mathrm{erg\,s^{-1}}$ within its cooling radius. Find the central cooling time and the classical mass deposition rate, then say what the spectroscopy actually shows.

Step 1. Central cooling time from equation (3), scaling the fiducial by the density and temperature ratios: $t_{\rm cool}=8.5\times10^{10}\,\mathrm{yr}\,(0.02/10^{-3})^{-1}(4\times10^{7}/10^{8})^{1/2}=8.5\times10^{10}\times(1/20)\times0.632\approx\mathbf{2.7\times10^{9}\,yr}$ — well under a Hubble time, so this is a genuine cool core.

Step 2. Classical mass deposition rate from equation (4). Write $k_BT=3.4\,\mathrm{keV}=5.45\times10^{-9}\,\mathrm{erg}$ , $\mu=0.6$ , $m_H=1.67\times10^{-24}\,\mathrm{g}$ : $\dot M = 2(0.6)(1.67\times10^{-24})(3\times10^{44})/[5(5.45\times10^{-9})]\approx2.2\times10^{28}\,\mathrm{g\,s^{-1}}$ .

Step 3. Convert to solar masses per year: $\dot M\approx2.2\times10^{28}\times(3.156\times10^{7})/(1.989\times10^{33})\approx\mathbf{350\,M_\odot\,yr^{-1}}$ — the classical prediction, enough to build a giant elliptical's stellar mass in a few gigayears.

Step 4. The reality check: XMM-Newton RGS spectra of such cores show essentially no gas cooling below $\sim T_{\rm vir}/3$ , and the inferred rate is a tenth of this, $\sim\mathbf{35\,M_\odot\,yr^{-1}}$ . The missing factor of ten is the cooling-flow problem; suppressing the other $\sim90\%$ is what AGN feedback must supply. Result: t_cool ≈ 2.7 Gyr, classical Ṁ ≈ 350 M⊙/yr, observed ≈ 35 M⊙/yr.

2 · Can cavities offset the cooling? (energy balance)

The same core hosts a pair of X-ray cavities, each of radius $R=15\,\mathrm{kpc}$ , in ambient pressure $P\approx3\times10^{-10}\,\mathrm{erg\,cm^{-3}}$ , rising on $t_{\rm buoy}\approx10^{8}\,\mathrm{yr}$ . Does their power match the cooling luminosity of Example 1?

Step 1. Cavity volume: $R=15\,\mathrm{kpc}=4.63\times10^{22}\,\mathrm{cm}$ , so $V=\tfrac{4}{3}\pi R^3\approx4.2\times10^{68}\,\mathrm{cm^3}$ per cavity.

Step 2. Enthalpy per cavity from equation (6): $E_{\rm cav}=4PV=4(3\times10^{-10})(4.2\times10^{68})\approx5.0\times10^{59}\,\mathrm{erg}$ ; for the pair, $E_{\rm tot}\approx\mathbf{1.0\times10^{60}\,erg}$ .

Step 3. Power: $L_{\rm cav}=E_{\rm tot}/t_{\rm buoy}=1.0\times10^{60}/(3.16\times10^{15}\,\mathrm{s})\approx\mathbf{3\times10^{44}\,erg\,s^{-1}}$ .

Step 4. Compare: $L_{\rm cav}\approx3\times10^{44}\,\mathrm{erg\,s^{-1}}\approx L_{\rm cool}$ from Example 1. The cavities carry just enough enthalpy, released on the buoyancy time, to offset the radiative losses — the heart of the feedback solution. For the extreme system MS 0735.6+7421 the same calculation with $\sim100$ -kpc cavities gives $\sim10^{62}\,\mathrm{erg}$ , among the most energetic AGN outbursts known. Result: E_tot ≈ 10⁶⁰ erg, L_cav ≈ 3×10⁴⁴ erg/s ≈ L_cool — heating balances cooling.

Physics of the portrait

Cooling flows are inseparable from the career of one physicist: Andrew Fabian, who built the classical model, championed it for two decades, then led the community in dismantling it when the spectra disagreed.

Origins · 1 minAndrew Fabian and how this was discovered

Andrew Fabian

b. 1948 · X-ray astronomy and cluster cooling flows

Andrew Fabian, at Cambridge's Institute of Astronomy, has shaped the study of hot cosmic atmospheres for half a century. In 1977, with Paul Nulsen, he formulated the steady-state cooling-flow model, and his authoritative 1994 Annual Reviews article, "Cooling Flows in Clusters of Galaxies," codified the field and its prediction of large mass deposition rates. Fabian was also central to the model's undoing: when XMM-Newton's Reflection Grating Spectrometer resolved the soft-X-ray lines of cool cores around 2001, the expected cascade of gas below $\sim T_{\rm vir}/3$ was absent, and Peterson, Fabian and collaborators quantified the deficit. Rather than defend the old picture, Fabian pivoted to its resolution, showing with deep Chandra imaging of Perseus that the central AGN inflates buoyant cavities and drives sound waves and weak shocks that reheat the gas. His work knit together the thermodynamics of gas on 100-kpc scales with accretion physics on sub-parsec scales — a model of a scientist following the data past his own most famous idea.

History

1965

— Field derives the isobaric thermal-instability criterion for a radiatively cooling gas

1977

— Cowie and Binney, and Fabian and Nulsen, formulate the steady cooling-flow model

1994

— Fabian's Annual Reviews article codifies the classical cooling-flow picture

2001

— XMM-Newton grating spectra reveal the missing cool gas: the cooling-flow problem

2005

— Chandra imaging of MS 0735.6+7421 reveals a roughly 10⁶² erg AGN outburst

2007

— McNamara and Nulsen review establishes AGN mechanical feedback as the resolution

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §27–28 (clusters, ICM, X-ray emission) · Binney and Tremaine, Galactic Dynamics, 2nd ed. · Mo, van den Bosch and White, Galaxy Formation and Evolution (ICM thermodynamics and feedback) · Dodelson, Modern Cosmology · Fabian, ARA&A 32, 277 (1994), "Cooling Flows in Clusters of Galaxies" · McNamara and Nulsen, ARA&A 45, 117 (2007), "Heating Hot Atmospheres with AGN" · Peterson and Fabian, Phys. Rep. 427, 1 (2006), "X-ray spectroscopy of cooling clusters." The $\sqrt{T}/n$ cooling-time scaling gives $v_{\rm in}\approx25\,\mathrm{km\,s^{-1}}$ and Mach $\approx0.03$ for the subsonic inflow. The cooling-time coefficient $8.5\times10^{10}\,\mathrm{yr}$ (at $n_e=10^{-3}\,\mathrm{cm^{-3}}$ , $T=10^{8}\,\mathrm{K}$ ) uses $\bar g_B\approx1.1$ ; Example 1 yields $t_{\rm cool}\approx2.7\,\mathrm{Gyr}$ and $\dot M\approx350\,M_\odot\,\mathrm{yr^{-1}}$ . The $4PV$ cavity enthalpy gives $E_{\rm tot}\approx10^{60}\,\mathrm{erg}$ and $L_{\rm cav}\approx3\times10^{44}\,\mathrm{erg\,s^{-1}}$ in Example 2.

Prerequisites. Thermal Bremsstrahlung · Hydrostatic Equilibrium · Supermassive Black Holes · Active Galactic Nuclei: Unified Model. Leads to. Relativistic Jets and Launching Mechanisms · Advection-Dominated Accretion Flows · M–sigma Relation and Black Hole–Galaxy Coevolution · Eddington Luminosity and Accretion Limits.

Ask the tutor: Ask me to walk through why XMM-Newton's missing cool gas forced the field to invoke black-hole feedback — or to work any cavity-power or cooling-time estimate with you step by step.

Why is the inflow called a "flow" rather than a collapse?

Because it is deeply subsonic and quasi-static. The cooling time is a few gigayears, but sound crosses the same region in only $\sim10^8$ years, so pressure disturbances re-level the atmosphere far faster than radiation drains it — the gas never builds a deficit large enough to fall freely. The settling speed at the cooling radius is only $\sim25\,\mathrm{km\,s^{-1}}$ against a sound speed near $950\,\mathrm{km\,s^{-1}}$ , a Mach number around $0.03$ . The gas creeps inward while staying almost perfectly hydrostatic — exactly what "flow" conveys. A collapse would need cooling to outrun the sound-crossing time, which does not happen on cluster scales.

What is the difference between a cool-core and a non-cool-core cluster?

It is a difference in central thermodynamic state, best read off the entropy index $K=k_BT/n_e^{2/3}$ . A cool core has sharply peaked density, an inward-dipping temperature, a short central cooling time ( $\lesssim1\,\mathrm{Gyr}$ ), and low central entropy ( $K_0\lesssim30\,\mathrm{keV\,cm^2}$ ) — everything a cooling flow needs, kept in check by feedback. A non-cool-core cluster has a flat, high-entropy core and no short-cooling-time region, so no cooling flow. The leading explanation is history: a major merger shock-heats and mixes the core, erasing the entropy peak, while relaxed clusters settle into the cool-core state. Roughly half of nearby clusters are cool-core.

If gas is not actually cooling and piling up, where does the radiated energy come from?

From the AGN, in a steady state. The core does radiate the X-rays we see, so on short timescales some gas is losing its enthalpy. But that energy is continuously replenished by mechanical heating from the black hole — jet-inflated cavities, weak shocks, sound waves and turbulence — rather than by a bottomless supply of gas cooling to zero. Time-averaged, heating nearly equals cooling, so only a small residual condenses (the observed $\sim10\%$ ). The classical flow drew its energy from gas cooling all the way down; the real system draws most of it from the black hole and returns the gas to the hot phase before it fully cools.

Why is the cavity energy 4PV and not just its internal energy?

Because a buoyant bubble gives back to the ICM both the heat it contains and the work it did to make room for itself. For a relativistic plasma with $\gamma=4/3$ , the internal energy is $U=PV/(\gamma-1)=3PV$ . Inflating the cavity against the ambient pressure required a further $PV$ of work, stored in the displaced atmosphere and recovered as the wake refills. The sum $3PV+PV=4PV$ is the enthalpy — the right budget for a bubble transporting energy through a pressurized medium. Using only the internal energy would undercount the heating and the inferred jet power by a quarter.

What does the precipitation criterion have to do with feeding the black hole?

Everything — it bridges cluster-scale gas and sub-parsec accretion. Where the ratio of cooling time to free-fall time drops below about ten, thermally unstable blobs cool faster than buoyancy can carry them into hotter gas and dissolve them. They become cold clouds and filaments (seen in H $\alpha$ , CO and dust) that rain toward the centre. Because cold, low-angular-momentum clouds shed angular momentum efficiently, they are captured stochastically at rates far above the smooth Bondi estimate — chaotic cold accretion. That fuelling powers the next burst of jets, which reheat the core, raise the cooling time, and shut off the rain until it resumes. The threshold thus sets the duty cycle of the thermostat.

How do we actually know the classical cooling rates are wrong?

From high-resolution X-ray spectroscopy. The classical isobaric cooling flow predicts that gas passing through every temperature emits a full ladder of lines, including strong low-ionization lines such as Fe XVII from gas well below the ambient temperature. When XMM-Newton's grating spectrometers resolved these lines in cool cores around 2001, the lines from gas below roughly a third of the virial temperature were largely missing. The spectra fit gas cooling only partway and then being reheated, not cascading to zero. Fitting the residual lines gives condensation rates about ten times below the classical value — a direct, model-independent sign that cooling is truncated.