Cluster Mass Determination and Hydrostatic Equilibrium

A galaxy cluster is the largest thing gravity has managed to bind: a few hundred to a few thousand galaxies, a diffuse ocean of X-ray-hot plasma at $T\sim10^7$ – $10^8\,\mathrm{K}$ called the intracluster medium (ICM), and — dominating everything — an invisible dark-matter halo, all totalling $10^{14}$ to a few $\times10^{15}\,M_\odot$ . Weighing that mass correctly is one of the most consequential measurements in cosmology, because the abundance of such objects is exponentially sensitive to how much matter the Universe holds and how clumpy it is. The trouble is that the mass is dark, so it must be inferred from tracers — three tracers, three principles. The galaxies move in the potential well, so their velocities weigh it (Zwicky's virial argument of 1933, sharpened into the spherical Jeans equation). The hot gas sits as a static atmosphere in the same well, so its pressure gradient weighs it (the X-ray hydrostatic method, the workhorse for individual clusters). And the whole mass bends passing light, so gravitational lensing weighs it with no assumption about dynamical state at all. This article builds all three but centres on the hydrostatic estimator, whose defining feature is economy: the enclosed mass depends only on the local logarithmic slopes of gas density and temperature, not on any absolute normalization. Where the three methods agree — to about $20\%$ — we have the cleanest evidence that dark matter is real and roughly collisionless; where they disagree, we read off the non-thermal pressure that makes cluster masses the leading systematic in cluster cosmology.

Key takeaways

Galaxy clusters are the most massive bound structures in the Universe; their dark-matter-dominated mass is measured three largely independent ways — X-ray gas, gravitational lensing, and galaxy dynamics.
The hydrostatic (HSE) method treats the hot ICM as a static ideal gas whose pressure gradient balances gravity, giving $M(<r)$ directly from the radial gas density and temperature profiles.
The dynamical method applies the virial theorem or the Jeans equation to member-galaxy velocities, treating the galaxies as a collisionless tracer of the same potential.
Lensing measures the projected mass with no equilibrium assumption, providing the cross-check that exposes the hydrostatic mass bias $b\equiv1-M_{\rm HSE}/M_{\rm lens}\approx0.1$ – $0.3$ from neglected non-thermal pressure.
Cluster mass functions probe $\Omega_m$ and $\sigma_8$ ; a $10$ – $20\%$ mass error propagates straight into the cosmological parameters and drives the Planck cluster-count tension.
Key idea: under spherical hydrostatic equilibrium the enclosed mass is fixed entirely by the local slopes of the gas density and temperature, $M(<r)=-\dfrac{k_B T\,r}{G\mu m_p}\left(\dfrac{\mathrm{d}\ln\rho_g}{\mathrm{d}\ln r}+\dfrac{\mathrm{d}\ln T}{\mathrm{d}\ln r}\right)$ — no absolute normalization of the gas is needed, only its logarithmic gradients.

Where this lives today

Cluster mass determination is in a golden age, driven by three surveys pulling together. In X-rays, eROSITA is delivering an all-sky catalogue of $\sim10^5$ clusters whose masses must be calibrated to percent-level accuracy. In the microwave, the South Pole Telescope, ACT, and Planck detect clusters through the Sunyaev–Zel'dovich effect — the inverse-Compton distortion of the CMB by the same hot electrons — whose signal is nearly redshift-independent and tied to the thermal energy HSE uses. And in the optical, weak-lensing surveys from Euclid, the Vera Rubin Observatory (LSST), and Subaru Hyper Suprime-Cam provide the assumption-free masses that pin down the bias $b$ . That bias is not academic: the $2\sigma$ tension between the parameters ( $\sigma_8$ , $\Omega_m$ ) from Planck's CMB and from Planck's cluster counts can be absorbed by $b\approx0.4$ — larger than most direct measurements — so whether the discrepancy signals new physics or merely mis-weighed clusters hinges on nailing $b$ . Meanwhile XRISM, the successor to Hitomi, now measures ICM turbulent velocities directly. And the Bullet Cluster remains the field's iconic image: in a merger where HSE fails, lensing maps the collisionless dark matter offset from the shocked X-ray gas — the cleanest demonstration that most of a cluster's mass neither shines nor collides.

The physical picture

Picture the hot ICM as a self-gravitating atmosphere sitting in the cluster's potential well, exactly like the atmosphere of a star or of the Earth. If the gas is neither collapsing nor streaming out, then at every radius the outward push of the pressure gradient must precisely cancel the inward pull of gravity. That single balance condition is the whole method: measure how fast pressure falls with radius, and you have weighed everything inside it. The gas is only a tracer — its density and temperature slopes report the depth and shape of a potential dominated by invisible dark matter, so the total mass it reveals is roughly six times the gas mass itself.

Figure

Hydrostatic balance in a cluster: the outward pressure-gradient force on each gas shell cancels inward gravity, so dP/dr = −ρ_g GM(r)/r². The measured inputs are the two logarithmic slopes — a steep density dlnρ/dlnr ≈ −2 and a nearly flat temperature dlnT/dlnr ≈ −0.4 for a k_B T ≈ 8 keV cluster

Because the estimator uses only logarithmic derivatives, it is insensitive to the absolute gas normalization or to distance-dependent calibration factors — but punishingly sensitive to the temperature gradient, since a noisy $\mathrm{d}\ln T/\mathrm{d}\ln r$ feeds straight into the mass, which is why deep X-ray spectroscopy of $T(r)$ is the hard part. The second figure shows what comparing the three methods buys you: the mass rises to $\sim10^{15}\,M_\odot$ , the gas is a thin skin on top of the dark matter, and the thermal hydrostatic mass sits systematically below the assumption-free lensing mass by the bias $b$ .

Figure

Enclosed-mass profiles: the lensing (total) mass reaches about 8×10¹⁴ M⊙ at r₅₀₀ ≈ 1.3 Mpc, the thermal hydrostatic mass runs low by the bias b = 1 − M_HSE/M_lens ≈ 0.2, and the gas is only f_gas ≈ 0.13 of the total. Dark matter dominates at every radius

The core equations

The starting point is the momentum equation of a fluid at rest in a spherical potential: the pressure gradient is what holds the gas up against gravity. For a relaxed cluster the flow vanishes, and only the radial component survives.

\frac{\mathrm{d}P}{\mathrm{d}r}=-\rho_g\,\frac{GM(<r)}{r^2} \qquad (1)

where: $P$ — thermal gas pressure at radius $r$ ; $\rho_g$ — intracluster gas mass density; $M(<r)$ — total gravitating mass inside $r$ (dark matter + gas + galaxies); $G$ — Newton's constant. It is a force balance per unit volume: the outward pressure-gradient force supports the weight $\rho_g\,g$ of the gas. To make the pressure observable we need the equation of state. The ICM is a fully ionized, optically thin, non-degenerate plasma, so it obeys the ideal gas law with a mean molecular weight that folds in the mix of electrons and ions.

P=\frac{\rho_g}{\mu m_p}\,k_B T \qquad (2)

where: $k_B$ — Boltzmann constant; $T$ — gas temperature (from X-ray spectroscopy); $\mu\approx0.6$ — mean molecular weight of a fully ionized primordial plasma; $m_p$ — proton mass. The factor $\rho_g/\mu m_p$ is the total particle number density; the pressure is set by all free particles, electrons and ions alike. Substituting (2) into (1) and dividing by $\rho_g$ eliminates the density normalization, leaving the master result — the equation that does the weighing.

M(<r)=-\frac{k_B T\,r}{G\mu m_p}\left(\frac{\mathrm{d}\ln\rho_g}{\mathrm{d}\ln r}+\frac{\mathrm{d}\ln T}{\mathrm{d}\ln r}\right) \qquad (3)

where: $\mathrm{d}\ln\rho_g/\mathrm{d}\ln r$ — logarithmic slope of the gas density (typically $-2$ to $-3$ ); $\mathrm{d}\ln T/\mathrm{d}\ln r$ — logarithmic slope of temperature (mildly negative outside the core). A temperature scale $k_B T$ times a dimensionless pair of slopes fixes the mass; nothing about the absolute gas density enters. For a working density model, X-ray astronomers overwhelmingly use the isothermal $\beta$ -model, a cored power law that fits surface-brightness profiles well and has an analytic deprojection.

\rho_g(r)=\rho_0\left[1+\left(\frac{r}{r_c}\right)^2\right]^{-3\beta/2} \qquad (4)

where: $\rho_0$ — central gas density; $r_c$ — core radius; $\beta$ — outer-slope parameter, $\approx2/3$ for many clusters and physically $\beta=\mu m_p\sigma^2/k_BT$ . At $r\gg r_c$ the density falls as $r^{-3\beta}$ , so $\beta=2/3$ gives $\rho_g\propto r^{-2}$ . The dynamical route replaces the gas with the galaxies. Their collisionless phase space obeys the spherical Jeans equation — the stellar-dynamics analogue of hydrostatic balance — with an extra term because galaxy orbits need not be isotropic.

\frac{1}{\rho_\star}\frac{\mathrm{d}(\rho_\star\sigma_r^2)}{\mathrm{d}r}+\frac{2\beta_a\sigma_r^2}{r}=-\frac{GM(<r)}{r^2} \qquad (5)

where: $\rho_\star$ — number density of the galaxy tracer; $\sigma_r$ — radial velocity dispersion; $\beta_a=1-\sigma_t^2/\sigma_r^2$ — velocity anisotropy. The anisotropy is the price of the method: $M$ and $\beta_a$ both shape the projected dispersion, so they are degenerate — the mass–anisotropy degeneracy. The crude but historically decisive limit of (5) is the virial estimate, from $2K+U=0$ applied to the whole cluster; it needs only a velocity dispersion and a size.

M_{\rm vir}\simeq\frac{5\,R_{\rm vir}\,\sigma_{\rm los}^2}{G} \qquad (6)

where: $\sigma_{\rm los}$ — line-of-sight galaxy velocity dispersion; $R_{\rm vir}$ — virial radius. The order-unity prefactor ( $3$ to $5$ , depending on the mass profile) is what a careful Jeans analysis pins down; the estimate is good to a factor of a few. The third, assumption-free route is lensing: the cluster deflects background light once the projected mass density exceeds a critical value fixed purely by geometry and the speed of light.

\Sigma_{\rm cr}=\frac{c^2}{4\pi G}\frac{D_s}{D_l D_{ls}} \qquad (7)

where: $\Sigma_{\rm cr}$ — critical surface mass density for strong lensing; $D_l,D_s,D_{ls}$ — angular-diameter distances to lens, source, and lens-to-source; $c$ — speed of light. Strong arcs appear where $\Sigma>\Sigma_{\rm cr}$ ; weak shear measures $\Sigma$ below threshold. No dynamical state is assumed. Comparing the thermal hydrostatic mass to the lensing mass defines the single most important systematic in cluster cosmology.

b\equiv1-\frac{M_{\rm HSE}}{M_{\rm lens}} \qquad (8)

where: $M_{\rm HSE}$ — thermal hydrostatic mass; $M_{\rm lens}$ — lensing (assumption-free) mass; $b$ — bias parameter, empirically $\approx0.1$ – $0.3$ . A positive $b$ means the thermal method under-weighs the cluster by ignoring non-thermal pressure.

Derivation AThe master hydrostatic mass equation

The goal is to convert the fluid force balance (1) into a formula a spectroscopist can evaluate directly, and to expose why only logarithmic slopes survive. We start from the Euler equation, impose the relaxed-cluster conditions that kill the flow, insert the ideal-gas equation of state, and divide out the density — the single algebraic step that makes the absolute gas normalization drop away, leaving a temperature scale multiplying two dimensionless gradients.

Step 1. The momentum equation of a self-gravitating fluid is $\rho_g(\partial\mathbf{v}/\partial t+\mathbf{v}\cdot\nabla\mathbf{v})=-\nabla P-\rho_g\nabla\Phi$ . For a relaxed cluster set $\mathbf{v}=0$ and $\partial\mathbf{v}/\partial t=0$ , so pressure alone balances gravity: $\nabla P=-\rho_g\nabla\Phi$ . Why: the method rests on the gas being a static atmosphere, not a wind or an infall — dropping the inertial terms is precisely the "hydrostatic" in HSE.

Step 2. Take the radial component under spherical symmetry and use $\mathrm{d}\Phi/\mathrm{d}r=GM(<r)/r^2$ to recover equation (1), $\mathrm{d}P/\mathrm{d}r=-\rho_g\,GM(<r)/r^2$ . Why: spherical symmetry turns the vector balance into a one-dimensional relation between the pressure slope and the enclosed mass, which is all a radial profile can constrain.

Step 3. Insert the ideal-gas pressure $P=\rho_g k_B T/\mu m_p$ and apply the product rule: $\dfrac{k_B}{\mu m_p}\!\left(T\dfrac{\mathrm{d}\rho_g}{\mathrm{d}r}+\rho_g\dfrac{\mathrm{d}T}{\mathrm{d}r}\right)=-\rho_g\dfrac{GM}{r^2}$ , taking $\mu$ constant. Why: the temperature gradient now enters on the same footing as the density gradient — a cooler-outward atmosphere changes the inferred mass.

Step 4. Divide by $\rho_g/r^2$ , solve for $M$ , and rewrite the derivatives logarithmically using $r\,\mathrm{d}\ln X/\mathrm{d}r=\mathrm{d}\ln X/\mathrm{d}\ln r$ . The $1/\rho_g$ cancels the leftover $\rho_g$ , giving $M(<r)=-\dfrac{k_B T\,r}{G\mu m_p}\!\left(\dfrac{\mathrm{d}\ln\rho_g}{\mathrm{d}\ln r}+\dfrac{\mathrm{d}\ln T}{\mathrm{d}\ln r}\right)$ . Why: this cancellation is the punchline — the absolute density has vanished, so calibration errors in $\rho_g$ do not bias the mass; only its shape matters. *Verified: dimensionally $[k_BT/G\mu m_p]=$ mass/length, times $r$ gives mass; reproduces Sarazin (1988) and Carroll & Ostlie §27.3; numerically $k_BT/(G\mu m_p)\cdot r=3.71\times10^{13}M_\odot\,(T/\mathrm{keV})(r/\mathrm{Mpc})$ for $\mu=0.6$ .*

Go deeper · Derivation B · 1 minClosed-form mass of the isothermal beta-model

Derivation BClosed-form mass of the isothermal beta-model

The aim is to evaluate the master equation for the standard $\beta$ -model density at constant temperature, obtaining the compact analytic mass profile X-ray astronomers actually fit. It needs only the logarithmic slope of the cored power law; with the gas isothermal, the temperature term vanishes, the mass follows in closed form, and its large-radius limit reveals the underlying isothermal sphere.

Step 1. Differentiate the log density $\ln\rho_g=\ln\rho_0-\tfrac{3\beta}{2}\ln[1+(r/r_c)^2]$ to get the slope $\dfrac{\mathrm{d}\ln\rho_g}{\mathrm{d}\ln r}=-3\beta\dfrac{r^2}{r_c^2+r^2}$ . Why: this single slope carries all the radial information the isothermal estimator needs; it runs from $0$ in the flat core to $-3\beta$ far out.

Step 2. Set $\mathrm{d}\ln T/\mathrm{d}\ln r=0$ (isothermal) so only the density term survives, and substitute into the master equation: $M(<r)=-\dfrac{k_B T\,r}{G\mu m_p}\!\left(-3\beta\dfrac{r^2}{r_c^2+r^2}\right)$ . Why: isothermality is the classic simplifying assumption for the outer ICM and isolates the density profile as the sole driver of the mass.

Step 3. Tidy up to $M(<r)=\dfrac{3\beta\,k_B T}{G\mu m_p}\dfrac{r^3}{r_c^2+r^2}$ . In the limit $r\gg r_c$ this becomes $M\to(3\beta k_B T/G\mu m_p)\,r$ , mass growing linearly with radius. Why: linear $M\propto r$ is the signature of a singular isothermal sphere with a flat velocity curve — the $\beta$ -model quietly contains that idealization, and the mass must be truncated at the virial radius or it diverges. *Verified: reproduces Sarazin (1988) Eq. (5-24); for $\beta=2/3$ , $k_BT=8$ keV, $r_c=0.25$ Mpc, $r=1.5$ Mpc it gives $M=8.7\times10^{14}M_\odot$ , a sensible massive-cluster value.*

Go deeper · Derivation C · 1 minNon-thermal pressure and the hydrostatic mass bias

Derivation CNon-thermal pressure and the hydrostatic mass bias

The purpose here is to show precisely why the thermal-only hydrostatic mass is a lower bound on the truth, and to derive the bias $b$ that quantifies the shortfall. The ICM is never perfectly relaxed: turbulence, bulk flows, cosmic rays, and magnetic fields all add pressure that supports the gas but is invisible to X-ray spectroscopy. We split the total pressure, re-run the balance, and read off the mass the thermal estimator misses.

Step 1. Write the total support as $P_{\rm tot}=P_{\rm th}+P_{\rm nt}$ and impose true hydrostatic balance $\mathrm{d}P_{\rm tot}/\mathrm{d}r=-\rho_g\,GM_{\rm true}/r^2$ . The thermal-only estimator instead uses $GM_{\rm HSE}/r^2=-(1/\rho_g)\,\mathrm{d}P_{\rm th}/\mathrm{d}r$ . Why: the spectroscopic temperature reports only $P_{\rm th}$ , so any analysis that stops at $T(r)$ silently sets $P_{\rm nt}=0$ .

Step 2. Subtract to isolate the missing piece: $M_{\rm true}=M_{\rm HSE}-\dfrac{r^2}{G\rho_g}\dfrac{\mathrm{d}P_{\rm nt}}{\mathrm{d}r}$ . Since $P_{\rm nt}$ falls outward, $\mathrm{d}P_{\rm nt}/\mathrm{d}r<0$ , so the correction is positive and $M_{\rm true}>M_{\rm HSE}$ . Why: neglected support that also declines with radius mimics extra gravity — leaving it out under-counts the true weight the potential is holding up.

Step 3. Define $b\equiv1-M_{\rm HSE}/M_{\rm true}$ . If the non-thermal pressure is a fixed fraction $f$ of the thermal pressure with a similar slope, then $b\approx f$ , so a $\sim10$ – $20\%$ turbulent-pressure fraction produces $b\approx0.1$ – $0.2$ . Why: this ties an abstract systematic to a measurable microphysical quantity — the turbulent velocity — which is exactly what Hitomi and XRISM set out to measure directly. *Verified: consistent with the empirical range $b\approx0.1$ – $0.3$ from lensing-to-X-ray comparisons and with hydrodynamical simulations; the sign matches the observation that $M_{\rm HSE}<M_{\rm lens}$ .*

Limiting cases

Case	Result	Meaning
Isothermal gas	$M(<r)=-\dfrac{k_B T\,r}{G\mu m_p}\dfrac{\mathrm{d}\ln\rho_g}{\mathrm{d}\ln r}$	With $T$ constant the temperature term drops and the mass follows the density slope alone — the classic $\beta$ -model limit.
Isothermal $\beta$ -model	$M(<r)=\dfrac{3\beta k_B T}{G\mu m_p}\dfrac{r^3}{r_c^2+r^2}$	Closed form; at $r\gg r_c$ grows linearly, $M\propto r$ , as for a singular isothermal sphere.
Non-thermal support	$M_{\rm true}=M_{\rm HSE}-\dfrac{r^2}{G\rho_g}\dfrac{\mathrm{d}P_{\rm nt}}{\mathrm{d}r}$	Extra turbulent/bulk pressure makes the thermal-only $M_{\rm HSE}$ a lower bound — the origin of bias $b\sim0.1$ – $0.3$ .
Merging / non-equilibrium	$\left\vert\partial v/\partial t\right\vert\sim GM/r^2$	In a major merger the inertial term is not negligible; HSE fails outright and X-ray masses can err by factors of order unity.

Worked examples

1 · Hydrostatic mass of a hot cluster (fully worked)

A relaxed cluster has $k_B T=8.0\,\mathrm{keV}$ at $r=1.5\,\mathrm{Mpc}$ , where the deprojected profiles give slopes $\mathrm{d}\ln\rho_g/\mathrm{d}\ln r=-2.1$ and $\mathrm{d}\ln T/\mathrm{d}\ln r=-0.4$ . Find the enclosed mass.

Step 1. Assemble the prefactor. With $\mu=0.6$ , the combination $k_B T\,r/(G\mu m_p)$ is cleanest in scaled units: $k_B T/(G\mu m_p)\cdot r=3.71\times10^{13}\,M_\odot\,(T/\mathrm{keV})(r/\mathrm{Mpc})$ . For $T=8\,\mathrm{keV}$ , $r=1.5\,\mathrm{Mpc}$ this is $3.71\times10^{13}\times8\times1.5=\mathbf{4.45\times10^{14}\,M_\odot}$ .

Step 2. Add the slopes: $\mathrm{d}\ln\rho_g/\mathrm{d}\ln r+\mathrm{d}\ln T/\mathrm{d}\ln r=-2.1+(-0.4)=-2.5$ . The master equation carries an overall minus sign, so the bracket contributes $+2.5$ .

Step 3. Multiply: $M(<1.5\,\mathrm{Mpc})=4.45\times10^{14}\times2.5=\mathbf{1.1\times10^{15}\,M_\odot}$ — a genuine massive cluster, comparable to Coma. Note how the temperature slope matters: had we (wrongly) taken the gas isothermal, we would have used only $-2.1$ and got $9.3\times10^{14}\,M_\odot$ , a $16\%$ underestimate. M(<1.5 Mpc) ≈ 1.1×10¹⁵ M⊙, dominated by dark matter.

2 · Zwicky's virial mass of Coma (fully worked)

Coma has a line-of-sight galaxy velocity dispersion $\sigma_{\rm los}\approx1000\,\mathrm{km\,s^{-1}}$ over a radius $R\approx1.5\,\mathrm{Mpc}$ . Estimate its mass from the virial theorem and compare to the light.

Step 1. From $2K+U=0$ with $K=\tfrac12 M\langle v^2\rangle$ and $U=-\alpha GM^2/R$ , we get $M=\langle v^2\rangle R/(\alpha G)$ . Only the line-of-sight component is observed, and for isotropic orbits $\langle v^2\rangle=3\sigma_{\rm los}^2$ ; taking $\alpha\approx1$ gives $M\simeq3\sigma_{\rm los}^2 R/G$ .

Step 2. Insert numbers in SI: $\sigma_{\rm los}=10^6\,\mathrm{m\,s^{-1}}$ , $R=1.5\,\mathrm{Mpc}=4.63\times10^{22}\,\mathrm{m}$ . Then $M\simeq3\times(10^6)^2\times4.63\times10^{22}/(6.67\times10^{-11})=\mathbf{2.1\times10^{45}\,kg}$ .

Step 3. Convert: $M/M_\odot=2.1\times10^{45}/1.99\times10^{30}\approx\mathbf{1.0\times10^{15}\,M_\odot}$ . The luminous mass falls short by more than an order of magnitude in stars alone — precisely the gap Zwicky attributed in 1933 to dunkle Materie, dark matter. M_vir ≈ 10¹⁵ M⊙, vastly exceeding the visible mass.

3 · Gas fraction as a cosmological check (worked with numbers)

At $r_{500}$ a massive cluster has $M_{\rm gas}=1.05\times10^{14}\,M_\odot$ and a thermal hydrostatic mass $M_{\rm HSE}=6.2\times10^{14}\,M_\odot$ . With a stellar fraction $f_\star\approx0.02$ , is this consistent with the cosmic baryon fraction $\Omega_b/\Omega_m$ , and what does it say about bias?

Step 1. Form the apparent gas fraction using the thermal mass: $f_{\rm gas}^{\rm HSE}=M_{\rm gas}/M_{\rm HSE}=1.05/6.2=\mathbf{0.17}$ . Adding stars gives a baryon fraction $f_{\rm gas}^{\rm HSE}+f_\star\approx0.17+0.02=0.19$ .

Step 2. Compare to cosmology. With $\Omega_b h^2=0.0224$ , $\Omega_m=0.315$ , $h=0.674$ , the cosmic value is $\Omega_b/\Omega_m=\mathbf{0.157}$ . The cluster's apparent $0.19$ exceeds this ceiling — impossible, since baryons cannot be over-concentrated relative to dark matter.

Step 3. Resolve it with bias. Requiring $f_{\rm gas}^{\rm true}+f_\star\le\Omega_b/\Omega_m$ forces $M_{\rm true}\ge M_{\rm gas}/(0.157-0.02)=1.05\times10^{14}/0.137=7.7\times10^{14}\,M_\odot$ , so $b=1-M_{\rm HSE}/M_{\rm true}\ge1-6.2/7.7=\mathbf{0.19}$ — a $\sim20\%$ hydrostatic bias, recovered from baryon budgeting alone and matching the lensing comparison. A gas fraction above Ω_b/Ω_m is a smoking gun for hydrostatic bias b ≈ 0.2.

Physics of the portrait

Cluster mass determination has one unambiguous founder. Before anyone measured a temperature profile or flew an X-ray telescope, a combative Swiss-American astronomer at Caltech pointed the virial theorem at a swarm of galaxies and concluded that most of the Universe's mass is invisible.

Origins · 1 minFritz Zwicky and how this was discovered

Fritz Zwicky

1898–1974 · virial mass and dark matter in galaxy clusters

Fritz Zwicky was among the most original and least diplomatic astrophysicists of the twentieth century, and cluster masses are his most enduring legacy. In 1933, studying the Coma cluster from Caltech, he applied the virial theorem to the observed spread of galaxy velocities — about $1000\,\mathrm{km\,s^{-1}}$ — and found that the mass needed to bind the cluster exceeded the visible mass by a factor of hundreds. He named the missing ingredient dunkle Materie, dark matter; the conclusion was essentially correct, though the community set it aside for a generation. In 1937 he went further, arguing that clusters would act as gravitational lenses and that this lensing could weigh them independently of any dynamical assumption — the very cross-check, realized only decades later, that now calibrates hydrostatic masses. With Walter Baade he also coined the term "supernova". His temperament was famous — he called rivals "spherical bastards" (bastards from every direction) — but his physics was visionary: two of the three modern pillars of cluster weighing, dynamics and lensing, trace to papers he wrote before the Second World War, and the X-ray method that completed the trio arrived only once satellites could see the hot gas he never knew was there.

History

1933

— Zwicky's virial analysis of Coma reveals missing mass, which he names dark matter

1937

— Zwicky proposes clusters as gravitational lenses to weigh their mass independently

1972

— the Uhuru satellite establishes extended X-ray emission from the hot intracluster medium

1976

— Cavaliere and Fusco-Femiano formulate the isothermal beta-model of the ICM

2006

— the Bullet Cluster shows the lensing mass offset from the X-ray gas

2016

— Hitomi bounds Perseus ICM turbulence to only a few percent of thermal pressure

Going further

References. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §27.3 (clusters of galaxies) · Binney & Tremaine, Galactic Dynamics, 2nd ed., Ch. 4 (Jeans equations) and §4.8 · Ryden, Introduction to Cosmology, 2nd ed., Ch. 8 (dark matter) · Sarazin, X-ray Emission from Clusters of Galaxies (1988) · Dodelson, Modern Cosmology, Ch. 9 (structure and clusters) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 6.

Provenance. The chain HSE (1) + ideal gas (2) → master mass equation (3) is standard (Sarazin 1988; Carroll & Ostlie §27.3); the $\beta$ -model mass and virial estimate follow by direct differentiation and integration, and the bias derivation from splitting $P_{\rm tot}=P_{\rm th}+P_{\rm nt}$ . Anchor values ( $\mu\approx0.6$ , prefactor $3.71\times10^{13}M_\odot$ per keV per Mpc, $f_{\rm gas}\approx0.10$ – $0.15$ , $\Omega_b/\Omega_m\approx0.157$ , Coma $\sigma\approx1000\,\mathrm{km\,s^{-1}}\Rightarrow M\sim10^{15}M_\odot$ , bias $b\approx0.1$ – $0.3$ ) are textbook and current-survey standard (Ryden Ch. 8, Dodelson Ch. 9, Planck). A verifier should re-derive (3), check the $\beta$ -model dimensions, and confirm the $8\,\mathrm{keV}$ mass $\approx1.1\times10^{15}M_\odot$ and the Coma virial estimate.

Prerequisites. Hydrostatic Equilibrium in Stars · Virial Theorem · Ideal Gas Law and Thermodynamics · Thermal Bremsstrahlung Emission · Gravitational Lensing Basics. Leads to. Cosmological Parameters from Cluster Counts · Sunyaev–Zel'dovich Effect · NFW Dark Matter Halo Profile · M–sigma Relation and Black Hole–Galaxy Coevolution.

Ask the tutor: Ask me to weigh a specific cluster from its temperature and density slopes, or to show exactly how a $20\%$ hydrostatic bias shifts the inferred value of $\sigma_8$ .

Why does only the temperature scale, and not the gas density, set the absolute mass?

Because the master equation is $M\propto T\times(\text{log slopes})$ . When you insert the ideal-gas pressure into the force balance and divide by $\rho_g$ , the density normalization $\rho_0$ cancels exactly — what survives is the shape of $\rho_g(r)$ , encoded in $\mathrm{d}\ln\rho_g/\mathrm{d}\ln r$ . Physically, the gas is a tracer: its gradients report the steepness of the potential, while the temperature sets the energy scale that converts those gradients into a mass. This is why $M_{\rm HSE}$ is robust to distance and calibration errors in the density but very sensitive to a noisy temperature profile.

If lensing needs no equilibrium assumption, why not just use lensing for everything?

Lensing is the cleanest mass in principle, but it has its own costs. It measures the projected mass, so you must deproject a two-dimensional surface density into three dimensions, which requires assumptions about geometry (a triaxial cluster seen along its long axis looks more massive). It is sensitive to all mass along the line of sight, including uncorrelated structure, which adds scatter. And weak lensing needs many resolved background galaxies, so it works best for massive, intermediate-redshift clusters. X-ray HSE, by contrast, gives smooth, high-signal profiles for individual relaxed clusters. The two are complementary: lensing calibrates the average bias of the cheaper X-ray and SZ masses used in large surveys.

What sets the mean molecular weight μ ≈ 0.6, and does it matter?

For a fully ionized plasma of roughly $75\%$ hydrogen and $24\%$ helium by mass, counting one proton and one electron per hydrogen and two electrons plus a nucleus per helium gives a mean mass per particle of about $0.6\,m_p$ . It matters because it enters the pressure $P=\rho_g k_B T/\mu m_p$ directly, hence linearly in the inferred mass. A $\mu$ error of a few percent is a few percent in $M$ — usually subdominant to the temperature-gradient uncertainty, but not negligible for precision work.

How does temperature actually get measured out to large radius?

From X-ray spectroscopy. The ICM radiates by thermal bremsstrahlung plus line emission from highly ionized metals; the continuum shape and line ratios fix the temperature in each radial annulus. In practice you extract a spectrum in concentric rings, fit a plasma-emission model, and build $T(r)$ . The catch is that surface brightness falls steeply outward, so the outskirts are photon-starved and $T(r)$ turns noisy exactly where the mass is still accumulating — which is why HSE masses are quoted out to $r_{500}$ (where the mean density is $500$ times critical) rather than the full virial radius, and why the temperature gradient dominates the error budget.

Is the hydrostatic assumption ever badly wrong, not just biased?

Yes — in mergers. The method assumes the flow velocity and its time derivative vanish, so the pressure gradient alone balances gravity. During a major merger, gas sloshes, shocks form, and bulk velocities approach the sound speed; the inertial term $\rho_g\,\partial v/\partial t$ is then no longer negligible, and the thermal pressure gradient no longer traces the potential. Masses of such systems can be wrong by factors of order unity, not the $10$ – $20\%$ of the relaxed-cluster bias. This is why cosmologists either select relaxed, symmetric systems by X-ray morphology or model disturbed ones hydrodynamically — and why the Bullet Cluster, a violent merger, is analysed with lensing, not HSE.

How does a small mass error become a cosmological problem?

Through the exponential sensitivity of the cluster mass function. The number of clusters above a given mass, especially at the high-mass end and at early times, depends steeply on $\sigma_8$ (the amplitude of matter fluctuations) and $\Omega_m$ . But you never observe mass directly; you observe an X-ray, SZ, or optical proxy and convert it with a calibration that carries the bias $b$ . A systematic $20\%$ mass error therefore masquerades as a shift in $\sigma_8$ — exactly the knob that reconciles Planck's cluster counts with its CMB result, where the tension can be absorbed by $b\approx0.4$ . Whether that value is physical, or points to new physics, is one of the live questions the current lensing surveys are built to answer.