Cosmic Structure

Thermal Sunyaev-Zel'dovich Cluster Counts as a Cosmology Probe

About one CMB photon in ten thousand that crosses a massive galaxy cluster gets kicked to higher energy by the cluster's 10-million-degree gas, carving a fractional dip of order 10⁻⁴ into the microwave background at frequencies below 217 GHz. That tiny, redshift-independent shadow is the thermal Sunyaev-Zel'dovich (tSZ) effect, and counting the clusters it reveals is one of the sharpest ways to weigh the growth of cosmic structure.

Thermal SZ cluster counts use surveys of these microwave shadows to build a census of galaxy clusters — the most massive gravitationally bound objects in the universe — as a function of mass and redshift. Because cluster abundance depends steeply on the amplitude of matter fluctuations (σ₈) and the matter density (Ω_m), the number of clusters found at each epoch is a direct probe of how gravity has assembled structure over the last 10 billion years.

  • TypeInverse-Compton CMB spectral distortion / cluster survey probe
  • RegimeHot intracluster gas, T_e ~ 10⁷–10⁸ K, kT_e ~ 1–15 keV
  • PredictedSunyaev & Zel'dovich, 1970–1972
  • Null frequency217 GHz (decrement below, increment above)
  • Key relationy = (σ_T k_B / m_e c²) ∫ n_e T_e dl ; Y_SZ ∝ M^5/3 E(z)^2/3
  • Observed inPlanck, ACT, SPT (mm/sub-mm surveys); thousands of clusters

Interactive visualization

Press play, or step through manually. The visualization is yours to drive — try it before reading on.

Open visualization fullscreen ↗

Watch the 60-second explainer

A condensed visual walkthrough — narrated, captioned, under a minute.

What it is: hot gas casting a microwave shadow

A galaxy cluster contains a diffuse atmosphere of ionized gas — the intracluster medium (ICM) — with electron temperatures of 10⁷–10⁸ K (thermal energies kT_e ≈ 1–15 keV) and electron densities around 10⁻³ cm⁻³. As cosmic microwave background (CMB) photons pass through, roughly 1% or less scatter off these fast electrons via inverse Compton scattering. On average the cold 2.725 K photons gain a little energy, so the CMB spectrum seen through the cluster is distorted: photons are shifted from the Rayleigh-Jeans (low-frequency) side to the Wien (high-frequency) side.

  • Below 217 GHz the cluster appears as a decrement — a cold spot against the CMB.
  • Above 217 GHz it becomes an increment — a hot spot.
  • At the crossover null (217 GHz) the effect nearly vanishes for non-relativistic gas.

Crucially, this is a fractional distortion of a fixed-temperature background, so its surface brightness does not dim with distance. That makes tSZ a nearly redshift-independent way to find clusters out to z > 1.5.

The mechanism: the Compton-y parameter

The strength of the distortion at a point on the sky is the dimensionless Compton-y parameter:

y = (σ_T k_B / m_e c²) ∫ n_e T_e dl

where σ_T is the Thomson cross-section, n_e the electron density, T_e the temperature, and the integral runs along the line of sight. Physically, y measures the total thermal energy of electrons integrated through the cluster — essentially the pressure of the gas, P_e = n_e k_B T_e, projected on the sky. Typical cluster central values are y ≈ 10⁻⁵ to 10⁻⁴.

The temperature change in the non-relativistic limit is ΔT/T_CMB = y · f(x), with x = hν/k_B T_CMB and f(x) = x·coth(x/2) − 4, which equals −2 in the Rayleigh-Jeans limit and passes through zero at x ≈ 3.83 (217 GHz). For very hot clusters, relativistic corrections shift the null by a few GHz and modify f(x) at the few-percent level.

Integrating y over the cluster's solid angle gives Y_SZ = ∫ y dΩ, a measure of the cluster's total thermal energy — and hence its mass.

Key quantities: the Y–M scaling relation and a worked example

Cosmology from cluster counts needs to convert the observable Y_SZ into a halo mass. Self-similar collapse predicts a tight power law:

E(z)^(−2/3) · D_A² · Y_500 ∝ M_500^(5/3)

where M_500 is the mass within a radius enclosing 500× the critical density, D_A is the angular-diameter distance, and E(z) = H(z)/H_0 encodes expansion. The 5/3 slope follows because Y ∝ M·T and the virial temperature scales as T ∝ M^(2/3).

  • A M_500 ≈ 5×10¹⁴ M_☉ cluster at z ≈ 0.3 has kT_e ≈ 5 keV and integrated Y_500 ≈ 10⁻⁴ Mpc² — comfortably above Planck's detection floor of a few ×10¹⁴ M_☉.
  • The most massive clusters reach ~2×10¹⁵ M_☉ with kT_e ≈ 12–15 keV.

The abundance of such clusters is set by the halo mass function (Press-Schechter / Tinker), which is exponentially sensitive to σ₈ — a 5% change in σ₈ can change the number of massive clusters by tens of percent.

How it's observed: mm-wave surveys and cluster catalogs

Detecting tSZ requires arcminute-resolution millimeter telescopes with multiple frequency bands to separate the SZ spectrum from the primary CMB, dusty galaxies, and radio sources. Three surveys dominate:

  • Planck (all-sky, 2013–2015): its multi-frequency coverage (100–857 GHz) exploited the 217 GHz null to build a catalog of ~1,650 SZ sources, ~1,200 confirmed clusters, mostly at z < 1.
  • ACT (Atacama Cosmology Telescope, Chile): ~1.4 arcmin resolution, thousands of clusters over ~13,000 deg², reaching z > 1.5.
  • SPT (South Pole Telescope): deep, high-resolution surveys (SPT-SZ, SPTpol, SPT-3G) that pioneered SZ-selected cluster cosmology with samples like the SPT-SZ 2500 deg² catalog.

Clusters are then followed up with optical/IR photometry for redshifts and with X-ray (Chandra, XMM, eROSITA) and weak-lensing data to calibrate mass. The observable-mass calibration — especially from weak gravitational lensing — is now the pacing item for these analyses.

tSZ counts sit alongside several cousins, each with distinct strengths:

  • X-ray counts depend on n_e² and suffer (1+z)⁴ surface-brightness dimming, so they excel at low z but fade at high z where tSZ still shines.
  • Kinetic SZ (kSZ) comes from the bulk peculiar motion of the gas (Doppler, ∝ v/c · τ_e), has a blackbody spectrum with no 217 GHz null, and probes velocities rather than temperature — it is ~10–100× weaker than tSZ.
  • Weak-lensing counts measure mass directly but with low signal-to-noise per object; they are the gold standard for calibrating SZ masses rather than a competing selection.
  • tSZ power spectrum and bispectrum use the diffuse, unresolved SZ background (C_ℓ^yy ∝ σ₈^~8) as a complementary, calibration-light probe.

Because tSZ abundance tracks late-time growth while the primary CMB fixes conditions at z ≈ 1100, comparing them tests whether structure grew as ΛCDM predicts — a powerful lever on dark energy and modified gravity.

Significance, tensions, and open questions

tSZ cluster counts deliver competitive constraints on the combination σ₈(Ω_m/0.27)^0.3 ≈ 0.78. ACT+WMAP found σ₈ ≈ 0.83 and Ω_m ≈ 0.29; SPT and Planck reach similar precision. But a famous puzzle emerged: Planck's SZ cluster counts preferred a lower σ₈ (~0.75) than Planck's own primary-CMB fit (~0.83).

  • The discrepancy hinges on the hydrostatic mass bias (1−b): X-ray masses may underestimate true masses by 20–40% because of non-thermal pressure (turbulence, bulk motions) the assumption of hydrostatic equilibrium ignores.
  • Reconciling SZ counts with the primary CMB requires (1−b) ≈ 0.6, i.e. a ~40% bias — larger than most weak-lensing calibrations find.

This mass-calibration uncertainty is the central open problem. Ongoing weak-lensing programs (DES, HSC, KiDS, Euclid, Rubin/LSST) and next-generation SZ surveys (SPT-3G, Simons Observatory, CMB-S4) aim to pin down (1−b), sharpening whether the σ₈ tension is a systematic or a genuine crack in ΛCDM.

Thermal SZ compared with related cluster-detection and CMB probes
ProbeSignal / observableRedshift dependencePrimary cosmology sensitivity
Thermal SZ countsCMB spectral distortion, Compton-y ∝ ∫ n_e T_e dlSurface brightness independent of zσ₈, Ω_m via cluster abundance growth
X-ray cluster countsBremsstrahlung, L_X ∝ n_e² Λ(T)Dims as (1+z)⁴ (cosmological dimming)σ₈, Ω_m; strong at low z
Kinetic SZDoppler shift ∝ (v/c) τ_eIndependent of zPeculiar velocities, reionization
Weak-lensing countsShear of background galaxiesDepends on source-lens geometryσ₈, Ω_m; direct mass, low S/N per cluster
Primary CMB (Planck)Acoustic peaks at z ≈ 1100Single early-universe snapshotΩ_m, σ₈ extrapolated from z ≈ 1100

Frequently asked questions

Why doesn't the thermal SZ signal get fainter for more distant clusters?

The tSZ effect is a fractional distortion of the CMB, whose temperature is essentially the same (scaling as (1+z) along with the photons) everywhere. Because it is a change in the ratio of scattered to background photons rather than a flux from the cluster itself, its surface brightness is independent of distance. This is why SZ surveys can detect massive clusters at z > 1.5 as easily as at z = 0.3, unlike X-ray surveys that dim as (1+z)⁴.

What is the significance of the 217 GHz null frequency?

217 GHz is where the thermal SZ distortion crosses zero for non-relativistic gas: below it clusters appear as CMB decrements (cold spots), above it as increments (hot spots). Observing at and around the null lets experiments like Planck cleanly separate the SZ signal from the primary CMB and from foregrounds, since the SZ signal has a unique, well-defined spectral shape that other components lack.

How do cluster counts constrain σ₈ and Ω_m?

The number of clusters above a given mass is set by the halo mass function, which is exponentially sensitive to σ₈ (the amplitude of matter fluctuations on 8 Mpc/h scales) and depends on Ω_m through the growth of structure and the volume element. Counting clusters as a function of mass and redshift therefore measures how fast gravity assembled structure, tightly constraining the degenerate combination σ₈(Ω_m/0.27)^0.3 ≈ 0.78.

What is the Y–M scaling relation and why does it matter?

Y_SZ, the integrated Compton-y, measures a cluster's total thermal energy and scales with mass as roughly E(z)^(−2/3) D_A² Y_500 ∝ M_500^(5/3). This tight, low-scatter relation lets surveys convert an easily-measured SZ signal into a halo mass. Its calibration — the normalization and any mass bias — is the dominant systematic in tSZ cosmology, because a biased mass scale shifts the inferred σ₈.

What is the hydrostatic mass bias and why does it cause tension?

Cluster masses derived from X-ray or SZ data often assume the gas is in hydrostatic equilibrium, but real clusters have turbulence and bulk motions (non-thermal pressure) that make the gas 'lighter' than it appears, so masses are underestimated by a factor (1−b) ≈ 0.8. Planck's SZ counts preferred a lower σ₈ than its primary-CMB fit; reconciling them needs (1−b) ≈ 0.6 (a ~40% bias), larger than most lensing measurements support — hence the ongoing tension.

How does thermal SZ differ from kinetic SZ?

Thermal SZ arises from the random thermal motions of hot electrons (∝ n_e T_e) and produces the characteristic spectral distortion with a 217 GHz null. Kinetic SZ arises from the cluster's bulk peculiar velocity (∝ v/c · τ_e), keeps a blackbody spectrum with no null, and is 10–100× weaker. Thermal SZ measures gas pressure and mass; kinetic SZ measures line-of-sight velocity, probing cosmic flows and reionization.