Galaxy Redshift Surveys

A galaxy redshift survey is a systematic program to measure angular positions and spectroscopic redshifts $z$ for a well-defined sample of galaxies over a chosen sky area and depth. Each redshift, read from the shift of spectral lines, is turned into a line-of-sight distance through a fiducial cosmology, so a two-dimensional catalogue on the sky becomes a three-dimensional map of where matter sits. Those maps delivered one of the twentieth century's great empirical surprises: galaxies are not sprinkled at random but woven into a cosmic web of filaments, sheets, clusters, and cavernous voids — the amplified fossil of the primordial density perturbations, developed by gravity across cosmic time. The power of a survey lies not in any single wall or void but in the statistics of the whole field. The two-point correlation function $\xi(r)$ measures the excess probability, over a random Poisson field, of finding a pair of galaxies separated by comoving distance $r$ ; its Fourier transform is the power spectrum $P(k)$ . On large scales both track the matter power spectrum fixed by early-universe physics and the transfer function, but seen through two filters: galaxies are biased tracers of mass, and the redshift used for distance mixes pure expansion with each galaxy's peculiar velocity. Two features convert these maps into precision cosmology: the baryon acoustic oscillations imprint a fixed comoving scale, the sound horizon at the drag epoch $r_{\mathrm{drag}}\approx147\ \mathrm{Mpc}$ , that acts as a standard ruler for distance and expansion rate; and the anisotropy of clustering in redshift space measures the linear growth rate $f=\mathrm{d}\ln D/\mathrm{d}\ln a$ , a direct test of gravity. From CfA through 2dFGRS, SDSS/BOSS/eBOSS, and now DESI, these datasets anchor the $\Lambda$ CDM concordance model.

Key takeaways

A galaxy redshift survey measures angular positions and spectroscopic redshifts $z$ for large samples, then converts them into a 3D map through a cosmological distance–redshift relation.
The primary products are the two-point correlation function $\xi(r)$ and its Fourier transform, the power spectrum $P(k)$ , quantifying how galaxies cluster relative to random.
Galaxies are biased tracers, so $P_{\mathrm{gal}}(k)=b^{2}P_{\mathrm{m}}(k)$ on large scales, and the redshift mixes expansion with peculiar velocity, adding redshift-space distortions.
The baryon acoustic feature at $r_{\mathrm{drag}}\approx147\ \mathrm{Mpc}$ ( $\approx100\,h^{-1}\mathrm{Mpc}$ ) is a standard ruler constraining $H(z)$ and $D_A(z)$ ; the Kaiser anisotropy yields the growth combination $f\sigma_8(z)$ , testing gravity.
From CfA through 2dFGRS, SDSS/BOSS/eBOSS, to DESI (tens of millions of redshifts), surveys now deliver percent-level cosmological constraints.
Key idea: by tagging every galaxy with a redshift, a survey builds a 3D map whose clustering statistics — $\xi(r)$ and $P(k)$ — carry both the shape of the matter power spectrum (early-universe physics + growth) and standard-ruler/standard-clock features (BAO, RSD) that constrain the cosmic expansion history and the growth of structure.

Where this lives today

Redshift surveys are the front line of observational cosmology. The Dark Energy Spectroscopic Instrument (DESI), which began its main survey in 2021 from Kitt Peak, uses 5000 robotic fibres to record thousands of spectra per exposure and charts the expansion history from $z\approx0.1$ out to $z\approx3.5$ through the Lyman- $\alpha$ forest; its early BAO results already deliver sub-percent distances and have hinted, at modest significance, that dark energy may not be a pure cosmological constant. The complementary Euclid space telescope, launched in 2023, combines BAO and redshift-space distortions with weak lensing. Analysis has moved beyond the bump: "full-shape" fits model the entire anisotropic $P(k,\mu)$ with effective field theory, extracting $H_0$ , $\Omega_m$ , and $f\sigma_8$ at once, while reconstruction runs the density field backwards to sharpen the peak that nonlinear evolution has smeared. The open frontiers are concrete — whether $f\sigma_8(z)$ tracks the general-relativistic prediction or reveals modified gravity, whether the equation of state $w$ deviates from $-1$ , and whether the mild survey–CMB tensions in clustering amplitude survive the next decade of data.

The physical picture

Begin with the raw idea of clustering. Drop points at random and the chance of a neighbour in a shell of volume $\mathrm{d}V$ at distance $r$ is just proportional to $\mathrm{d}V$ — no memory, no structure. Galaxies are not random: gravity has pulled them together for billions of years, so there is an excess probability of finding a pair at separation $r$ , and that excess is exactly the correlation function $\xi(r)$ . Measuring how $\xi$ falls off with $r$ measures the accumulated memory of gravity acting on the primordial fluctuations. The power spectrum $P(k)$ is the same information sorted by wavelength: large scales ( $k\to0$ ) carry the near-scale-invariant primordial spectrum, the broad turnover near $k\sim0.02\,h\,\mathrm{Mpc}^{-1}$ marks matter–radiation equality, and the small-scale slope encodes the transfer function and the cold-dark-matter content. A redshift is a cheap third dimension: instead of a hard direct distance, one takes a spectrum and reads the shift of known lines. But that shift is contaminated by each galaxy's own motion. Inside a collapsed cluster, virial velocities of $\sim10^{3}\ \mathrm{km\,s^{-1}}$ smear members into radial "fingers" pointing at the observer; on larger scales, coherent infall compresses structure along the line of sight — the Kaiser effect. The distortion is not noise: its amplitude measures the growth rate of structure. The first figure makes these distortions concrete for a single cluster.

Figure

Real space versus redshift space for one cluster: small-scale virial motions stretch it radially into a finger of God pointing at the observer, while large-scale coherent infall squashes structure along the line of sight (the Kaiser effect) — the same peculiar velocities that spoil cz/H₀ as a distance also encode the growth rate f

The BAO feature is a second gift in the same map — a ruler frozen into the galaxy distribution. Before recombination the photon–baryon plasma rang with sound waves; each initial overdensity launched a spherical shell that travelled a fixed comoving distance $r_{\mathrm{drag}}$ before the baryons were released, and galaxies later formed preferentially both at the original peak and on that shell. The result is a faint excess of pairs at that separation — a localized bump in $\xi(r)$ near $100\,h^{-1}\mathrm{Mpc}$ and a train of wiggles in $P(k)$ . Because the scale is set by well-understood linear physics, seeing the bump at redshift $z$ tells you the distance to that redshift. The second figure is the observable itself.

Figure

The galaxy two-point correlation function scaled as r² ξ(r): a smooth broadband decline carries a localized acoustic bump near r ≈ 100 h⁻¹Mpc, the comoving sound-horizon scale r_drag ≈ 147 Mpc frozen into the clustering — the standard ruler first detected by SDSS luminous red galaxies in 2005

The core equations

Everything a survey reports is built on the density contrast — the dimensionless field measuring how far the local galaxy count strays from the mean, and the quantity gravity amplifies.

1+\delta(\mathbf{x})=\frac{n(\mathbf{x})}{\bar n},\qquad \delta(\mathbf{x})=\frac{n(\mathbf{x})-\bar n}{\bar n} \qquad (1)

where: $n(\mathbf{x})$ — local number density of galaxies; $\bar n$ — mean density over the survey; $\delta$ — dimensionless density contrast, with $\langle\delta\rangle=0$ . Reading: $\delta=0$ is the mean, $\delta>0$ an overdensity, $\delta=-1$ a completely empty region — the floor voids press against. The correlation function is defined operationally through pair counts: the fractional excess, over a random field of the same mean density, in the joint probability of finding galaxies in two small volumes a distance $r$ apart.

\mathrm{d}P_{12}=\bar n^{2}\,[\,1+\xi(r)\,]\,\mathrm{d}V_1\,\mathrm{d}V_2 \qquad (2)

where: $\mathrm{d}P_{12}$ — joint probability of a galaxy in each of two volume elements; $\xi(r)$ — two-point correlation function, $\xi=0$ for a random field; $r$ — comoving pair separation. Reading: $\xi>0$ is clustering above random, and the classic $\xi(r)=(r/r_0)^{-\gamma}$ with $r_0\approx5\,h^{-1}\mathrm{Mpc}$ , $\gamma\approx1.8$ summarizes small-scale clustering. The same content in Fourier space defines the power spectrum. Statistical homogeneity forces different modes to be uncorrelated, so the two-point function in $k$ -space is diagonal with amplitude $P(k)$ .

\langle\tilde\delta(\mathbf{k})\,\tilde\delta^{*}(\mathbf{k}')\rangle=(2\pi)^3\,\delta_{\mathrm D}(\mathbf{k}-\mathbf{k}')\,P(k) \qquad (3)

where: $\tilde\delta(\mathbf{k})$ — Fourier transform of $\delta$ ; $\delta_{\mathrm D}$ — Dirac delta enforcing statistical homogeneity; $P(k)$ — power spectrum, depending only on $k=|\mathbf{k}|$ by isotropy. Reading: $P(k)$ is the variance per mode at wavenumber $k$ ; large $P$ means strong clustering on scale $\sim2\pi/k$ . The correlation function and power spectrum are one fact, related by a Fourier (Hankel) transform; isotropy collapses the 3D transform to a single radial integral over the spherical Bessel function $j_0$ .

\xi(r)=\int\frac{\mathrm{d}^3k}{(2\pi)^3}\,P(k)\,e^{i\mathbf{k}\cdot\mathbf{r}}=\int_0^\infty\frac{k^2\,\mathrm{d}k}{2\pi^2}\,P(k)\,\frac{\sin kr}{kr} \qquad (4)

where: $\dfrac{\sin kr}{kr}=j_0(kr)$ — spherical Bessel function surviving the angular integration. Reading: $P(k)$ makes the BAO wiggles and equality turnover explicit; $\xi(r)$ localizes the BAO as a single bump near $100\,h^{-1}\mathrm{Mpc}$ and keeps window effects additive. Because a real catalogue lives inside a ragged mask, $\xi(r)$ is estimated by comparing observed pairs against a synthetic random catalogue of the same geometry — the Landy–Szalay estimator, a minimum-variance, edge-corrected combination.

\hat\xi(r)=\frac{DD(r)-2\,DR(r)+RR(r)}{RR(r)} \qquad (5)

where: $DD$ — data–data pairs at separation $r$ ; $DR$ — data–random pairs; $RR$ — random–random pairs from a catalogue matching the survey geometry. Reading: the numerator vanishes when data look random ( $\hat\xi\to0$ ), and the $-2DR+RR$ combination cancels leading edge and integral-constraint biases. The most important complication is that the map is built in redshift space. Linear peculiar-velocity infall boosts modes along the line of sight, giving the anisotropic Kaiser power spectrum — the equation that turns a distortion into a growth measurement.

P^{s}(k,\mu)=\big(b+f\mu^{2}\big)^{2}\,P_{\mathrm{m}}(k) \qquad (6)

where: $b$ — linear galaxy bias, $\delta_{\mathrm{gal}}=b\,\delta_{\mathrm{m}}$ ; $f=\mathrm{d}\ln D/\mathrm{d}\ln a\approx\Omega_m(a)^{0.55}$ — growth rate; $\mu$ — cosine of the angle between $\mathbf{k}$ and the line of sight; $P_{\mathrm{m}}$ — real-space matter power spectrum. Reading: line-of-sight modes ( $\mu=1$ ) are boosted by $(b+f)^2$ , transverse modes ( $\mu=0$ ) only by $b^2$ ; the anisotropy measures $f$ . What surveys report from that anisotropy is the bias-robust growth combination $f\sigma_8$ , together with geometric BAO distances; folding the growth rate into the clustering amplitude makes the poorly known bias drop out.

f\sigma_8(z)=\frac{\mathrm{d}\ln D}{\mathrm{d}\ln a}\,\sigma_8\,D(z),\qquad D_V(z)=\Big[\,c\,z\,\frac{(1+z)^2 D_A^2(z)}{H(z)}\,\Big]^{1/3} \qquad (7)

where: $\sigma_8$ — rms linear density fluctuation in $8\,h^{-1}\mathrm{Mpc}$ spheres today; $D(z)$ — growth factor normalized to $D(0)=1$ ; $D_V$ — spherically-averaged BAO distance; $D_A(z)$ — angular diameter distance; $H(z)$ — Hubble parameter. Reading: $f\sigma_8(z)$ is the survey's growth report card, $D_V(z)/r_{\mathrm{drag}}$ its geometric one when the signal cannot split radial from transverse. Finally, every reported number carries an error budget set by the finite number of modes — cosmic variance on large scales, Poisson shot noise from discrete tracers — the formula that drives survey design toward larger volumes and denser sampling.

\sigma^2_{P}(k)=\frac{2}{N_k}\left[P(k)+\frac{1}{\bar n}\right]^2 \qquad (8)

where: $N_k$ — number of independent $\mathbf{k}$ modes in the bin; $1/\bar n$ — Poisson shot-noise floor. Reading: where $\bar n P\gg1$ the error is set by mode count (survey volume), so shrinking error bars means surveying more volume, not observing longer.

Derivation AFrom catalogue to power spectrum — the FKP estimator

The goal here is to bridge the clean definition of $P(k)$ , valid for an infinite continuous field, and the messy reality of a survey: a finite set of discrete galaxies seen only inside a complicated mask and radial selection, contaminated by shot noise. We want one estimator that is unbiased, reduces to the definition in the ideal limit, and weights each region optimally so the error bar is as small as the data allow. The Feldman–Kaiser–Peacock (FKP) construction achieves this by building a weighted overdensity field, Fourier transforming, squaring, and subtracting the Poisson self-pair term.

Step 1. Build the weighted overdensity field $F(\mathbf{x})=w(\mathbf{x})\,[\,n(\mathbf{x})-\alpha\,n_{\mathrm{r}}(\mathbf{x})\,]$ , where $n_{\mathrm{r}}$ is a dense random catalogue tracing the survey selection and $\alpha=N_{\mathrm{data}}/N_{\mathrm{rand}}$ normalizes it. Why: subtracting the mean field encoded by the randoms turns raw counts into a fluctuation field with zero mean inside the mask, so the estimator responds to clustering, not to the survey shape.

Step 2. Choose the optimal weight $w(\mathbf{x})=1/[\,1+\bar n(\mathbf{x})\,P_0\,]$ , with $P_0$ the power amplitude at the scale of interest. Why: it balances the two noise sources — dense regions ( $\bar n P_0\gg1$ ) count equally, sparse shot-noise-dominated regions are down-weighted — so no part of the survey injects excess noise.

Step 3. Fourier transform and average the squared amplitude over a shell in $|\mathbf{k}|$ : $\hat P(k)\propto\langle|\tilde F(\mathbf{k})|^2\rangle$ . Why: by equation (3) the ensemble-averaged squared amplitude is the power spectrum, so shell-averaging over independent modes trades the ensemble average we cannot do for a mode average we can.

Step 4. Subtract the shot-noise bias $P_{\mathrm{shot}}=\int w^2\bar n\,\mathrm{d}^3x\big/\int w^2\bar n^2\,\mathrm{d}^3x$ , leaving $\hat P(k)$ convolved with the window $|\tilde W(\mathbf{k})|^2$ . Why: discrete self-pairs add a constant positive offset to $|\tilde F|^2$ unrelated to clustering; it must be removed, and the residual window convolution forward-modelled rather than deconvolved. *In the full-sky, high-density limit $w\to\mathrm{const}$ , $\alpha\to0$ , the estimator collapses to $\langle|\tilde\delta|^2\rangle=P(k)$ from equation (3), confirming it reduces to the clean definition in the ideal case.*

Go deeper · Derivation B · 1 minThe Kaiser factor and the growth rate

Derivation BThe Kaiser factor and the growth rate

The aim is to show precisely how a distortion becomes a measurement — how plotting galaxies at the wrong radial position, because their redshift includes a peculiar velocity, imprints a specific, calculable anisotropy on the clustering that is proportional to the growth rate $f$ . We start from the coordinate map between real and redshift space, linearize it, invoke linear-theory continuity to relate velocity to density, and read off the $(b+f\mu^2)$ boost; then we expand its square in Legendre multipoles to isolate a ratio that depends only on $\beta=f/b$ , delivering the growth rate independently of the clustering amplitude.

Step 1. Map real to redshift space along the line of sight, $s=r+v_\parallel/(aH)$ , and linearize the Jacobian: $\delta^{s}(\mathbf{k})=\delta(\mathbf{k})-(ik_\parallel/aH)\,v_\parallel(\mathbf{k})$ . Why: the redshift is displaced from the true position only along the line of sight, so the distortion appears with an explicit $k_\parallel=\mu k$ — the geometric origin of the $\mu$ -dependence.

Step 2. Use linear continuity, $\theta\equiv\nabla\cdot\mathbf{v}=-aHf\,\delta_{\mathrm{m}}$ , to trade velocity for density, and combine with $\delta_{\mathrm{gal}}=b\,\delta_{\mathrm{m}}$ to get $\delta^{s}_{\mathrm{gal}}=(b+f\mu^2)\,\delta_{\mathrm{m}}$ . Why: in linear theory the velocity is sourced by the same growing mode as the density, with proportionality exactly $f$ — so the distortion carries direct information about how fast structure grows.

Step 3. Square and factor out $b^2$ : $P^s(k,\mu)=b^2(1+\beta\mu^2)^2 P_{\mathrm{m}}(k)$ with $\beta\equiv f/b$ , then expand in even Legendre polynomials, $P^s=\sum_{\ell=0,2,4}P_\ell^s\mathcal{L}_\ell(\mu)$ . Why: only even multipoles survive by the $\mu\to-\mu$ symmetry, and the multipole moments are the compressed data vector surveys fit.

Step 4. Integrate $(1+\beta\mu^2)^2$ against each $\mathcal{L}_\ell$ to get $P_0^s=(1+\tfrac23\beta+\tfrac15\beta^2)b^2P_{\mathrm{m}}$ and $P_2^s=(\tfrac43\beta+\tfrac47\beta^2)b^2P_{\mathrm{m}}$ , so $P_2^s/P_0^s$ depends only on $\beta$ . Why: the common $b^2P_{\mathrm{m}}$ prefactor cancels in the ratio, isolating $\beta=f/b$ ; combined with the amplitude it yields $f\sigma_8$ without needing $b$ and $\sigma_8$ separately. *Direct integration gives the coefficients $\{1,\tfrac23,\tfrac15\}$ for the monopole and $\{\tfrac43,\tfrac47\}$ for the quadrupole; as a check, the ratio $\to0$ as $\beta\to0$ , correctly recovering isotropy when there are no peculiar velocities.*

Go deeper · Derivation C · 1 minBAO as a standard ruler — reading distance and expansion rate

Derivation CBAO as a standard ruler — reading distance and expansion rate

The purpose is to explain how one fixed length — the comoving sound horizon $r_{\mathrm{drag}}$ imprinted on the galaxy distribution — is converted into two independent cosmological distances by measuring the acoustic feature separately across the sky and along the line of sight. Its apparent transverse angle measures the angular diameter distance while its apparent radial extent in redshift measures the Hubble rate; when the data are too noisy to separate the two, a single volume-averaged distance is measured instead. This is the geometric heart of why redshift surveys constrain dark energy.

Step 1. For the transverse mode, a comoving separation $r_\perp$ subtends angle $\Delta\theta$ through the angular diameter distance: $r_\perp=(1+z)\,D_A(z)\,\Delta\theta=D_M(z)\,\Delta\theta$ . Why: $D_A$ is defined so that a known transverse length maps to an observed angle — exactly what an on-sky pair separation measures.

Step 2. For the radial mode, a comoving separation $r_\parallel$ corresponds to a redshift interval, $r_\parallel=c\,\Delta z/H(z)$ . Why: line-of-sight separations are observed as redshift differences, and the conversion factor is exactly the Hubble rate — so the radial ruler is a direct clock on $H(z)$ .

Step 3. Set both equal to the known ruler, $r_\perp=r_\parallel=r_{\mathrm{drag}}$ , giving $\Delta\theta=r_{\mathrm{drag}}/D_M(z)$ and $\Delta z=r_{\mathrm{drag}}\,H(z)/c$ . Why: the same acoustic scale seen in two directions yields two independent observables, $D_M/r_{\mathrm{drag}}$ and $H\,r_{\mathrm{drag}}$ — why anisotropic BAO breaks degeneracies isotropic BAO cannot.

Step 4. When signal-to-noise permits only one number, combine into $D_V(z)=[\,c\,z\,D_M^2(z)/H(z)\,]^{1/3}$ . Why: the volume-average weights the two-to-one ratio of transverse to radial modes, and with $r_{\mathrm{drag}}\approx147\ \mathrm{Mpc}$ calibrated by the CMB it fixes absolute distances, not mere ratios. *For fiducial flat $\Lambda$ CDM ( $\Omega_m=0.31$ , $H_0=67.4$ ) at $z=0.35$ the formulas give $D_V=1379\ \mathrm{Mpc}$ , in excellent agreement with the Eisenstein et al. (2005) SDSS detection $D_V\approx1370\ \mathrm{Mpc}$ , with $r_{\mathrm{drag}}\approx147\ \mathrm{Mpc}$ from Planck 2018.*

Limiting cases

The clustering statistics behave very differently across scale and across the transition into the nonlinear regime; the table collects the four regimes a survey analyst must keep separate.

Case	Result	Meaning
Large scales (linear)	$P_{\mathrm{gal}}(k)=b^2\,P_{\mathrm{m}}(k),\quad k\lesssim0.1\,h\,\mathrm{Mpc}^{-1}$	Constant linear bias; clustering tracks the linear matter power spectrum directly, and the BAO wiggles live here.
Redshift space, linear (Kaiser)	$P^{s}(k,\mu)=\big(b+f\mu^{2}\big)^{2}P_{\mathrm{m}}(k)$	Coherent infall boosts line-of-sight ( $\mu=1$ ) modes; the anisotropy yields $f=\mathrm{d}\ln D/\mathrm{d}\ln a$ .
Small scales (nonlinear)	$\xi(r)\propto r^{-\gamma},\ \gamma\approx1.8,\ r_0\approx5\,h^{-1}\mathrm{Mpc}$	Classic power-law clustering; nonlinear collapse and virial fingers-of-God dominate, needing halo-model or simulation treatment.
BAO scale	$r_{\mathrm{BAO}}\simeq r_{\mathrm{drag}}\approx147\ \mathrm{Mpc}\ (\approx100\,h^{-1}\mathrm{Mpc})$	A localized bump in $\xi(r)$ , robust to broadband systematics, making it the cleanest standard ruler in a survey.

Common trapRedshift is not distance, and galaxies are not mass

Two seductive shortcuts wreck a clustering analysis. The first is treating $cz/H_0$ as a clean comoving distance: the measured redshift sums the cosmological (Hubble) shift and the Doppler shift from peculiar velocity, so mapping with $cz/H_0$ places every galaxy in redshift space, systematically deformed — Kaiser squashing on large scales, radial fingers-of-God on small scales. These distortions carry the growth rate but must be modelled, never assumed away. The second shortcut is reading the clustering amplitude as the matter amplitude. Galaxies are biased, $\delta_{\mathrm{gal}}=b\,\delta_{\mathrm{m}}$ only at leading order, so $P_{\mathrm{gal}}=b^2P_{\mathrm{m}}$ hides an unknown $b$ degenerate with $\sigma_8$ . This is why surveys report the bias-cancelling $f\sigma_8$ and the geometric BAO scale, not the raw amplitude — and why cross-correlation with lensing is needed to pin $b$ down.

Worked examples

1 · The BAO standard ruler at z = 0.35 (fully worked)

Adopt flat $\Lambda$ CDM with $\Omega_m=0.31$ , $\Omega_\Lambda=0.69$ , $H_0=67.4\ \mathrm{km\,s^{-1}Mpc^{-1}}$ , $c=2.998\times10^{5}\ \mathrm{km\,s^{-1}}$ , $r_{\mathrm{drag}}=147\ \mathrm{Mpc}$ . At the SDSS luminous-red-galaxy effective redshift $z=0.35$ , compute the expansion rate, the comoving and angular diameter distances, the BAO distance $D_V$ , and the apparent size of the acoustic feature.

Step 1. Evaluate the dimensionless expansion rate $E(z)=\sqrt{\Omega_m(1+z)^3+\Omega_\Lambda}$ . With $(1.35)^3=2.460$ , $E(0.35)=\sqrt{0.31\times2.460+0.69}=\sqrt{1.4527}=1.205$ , so the Hubble rate is $H(0.35)=H_0E=67.4\times1.205=\mathbf{81.2\ km\,s^{-1}Mpc^{-1}}$ .

Step 2. The Hubble distance is $D_H=c/H_0=4448\ \mathrm{Mpc}$ . Integrating $D_C=D_H\int_0^{0.35}\mathrm{d}z'/E(z')$ numerically gives $\int=0.3204$ , so $D_C=D_M=\mathbf{1425\ Mpc}$ (flat, hence $D_M=D_C$ ), and $D_A=D_M/(1+z)=\mathbf{1056\ Mpc}$ .

Step 3. Assemble $D_V=[\,c\,z\,D_M^2/H(z)\,]^{1/3}$ . The radial factor is $cz/H=1292\ \mathrm{Mpc}$ and $D_M^2=2.031\times10^6\ \mathrm{Mpc^2}$ , so $D_V=(2.624\times10^9)^{1/3}=\mathbf{1379\ Mpc}$ — in excellent agreement with the Eisenstein et al. (2005) measurement $D_V\approx1370\ \mathrm{Mpc}$ .

Step 4. Predict the feature size. Transverse: $\Delta\theta=r_{\mathrm{drag}}/D_M=147/1425=0.1032\ \mathrm{rad}=\mathbf{5.9^\circ}$ . Radial: $\Delta z=r_{\mathrm{drag}}H/c=147\times81.2/(2.998\times10^5)=\mathbf{0.040}$ . The dimensionless ruler measurement is $D_V/r_{\mathrm{drag}}=\mathbf{9.38}$ . So the same $147\ \mathrm{Mpc}$ shell appears as a $\sim6^\circ$ ring on the sky and a $\Delta z\approx0.04$ shell in depth — two independent handles on cosmic geometry from one feature.

2 · Redshift-space anisotropy and the growth rate

A BOSS-like sample at $z=0.5$ has linear bias $b=2$ . Using $f\approx\Omega_m(a)^{0.55}$ with $\Omega_m=0.31$ , $\Omega_\Lambda=0.69$ , find the growth rate, the Kaiser parameter $\beta$ , and the quadrupole-to-monopole ratio of the redshift-space power spectrum, then estimate $f\sigma_8$ .

Step 1. The matter fraction at scale factor $a$ is $\Omega_m(a)=\Omega_m(1+z)^3/E^2(z)$ . At $z=0.5$ , $(1.5)^3=3.375$ , so $\Omega_m(1+z)^3=1.046$ , $E^2=1.046+0.69=1.736$ , and $\Omega_m(a)=1.046/1.736=0.603$ . Hence $f=0.603^{0.55}=\mathbf{0.76}$ (compared with $f=0.31^{0.55}=0.53$ today).

Step 2. The Kaiser distortion parameter is $\beta=f/b=0.76/2=\mathbf{0.38}$ . This sets the strength of the redshift-space anisotropy relative to the isotropic clustering.

Step 3. Use the multipole coefficients from the Kaiser derivation. Monopole: $1+\tfrac23\beta+\tfrac15\beta^2=1+0.253+0.029=1.281$ . Quadrupole: $\tfrac43\beta+\tfrac47\beta^2=0.507+0.082=0.587$ . Their ratio is $P_2^s/P_0^s=\mathbf{0.46}$ — a large, readily measurable anisotropy that pins $\beta$ , and through $f=\beta b$ the growth rate.

Step 4. Fold in the amplitude. With $\sigma_8=0.81$ today and growth factor $D(0.5)\approx0.74$ , the amplitude at that epoch is $\sigma_8(z)=\sigma_8 D(z)=0.60$ , so $f\sigma_8(0.5)=0.76\times0.60=\mathbf{0.46}$ — consistent with the BOSS measurement $f\sigma_8\approx0.46$ near $z=0.5$ , a clean test that general relativity's growth prediction holds there.

Physics of the portrait

Galaxy redshift surveys are, at bottom, statistical mechanics applied to the universe — and the person who built that language, more than any other, is Jim Peebles. Every symbol above traces to his program of turning the galaxy distribution from a gallery of curiosities into a rigorously computable random field.

Origins · 1 minJim Peebles and how this was discovered

Jim Peebles

b. 1935 · statistical cosmology of large-scale structure

Phillip James Edwin "Jim" Peebles, born in Winnipeg in 1935 and a Princeton physicist for his entire career, is the principal architect of modern physical cosmology, and the theory of galaxy clustering is one of his deepest imprints. Arriving in Robert Dicke's group, he helped predict the cosmic microwave background just before its 1965 discovery, then spent the 1970s asking how a nearly smooth early universe came to be full of galaxies. His answer was to treat the galaxy distribution as a realization of a random field, characterized by its statistical moments — above all the two-point correlation function $\xi(r)$ and its Fourier partner $P(k)$ . The 1980 monograph The Large-Scale Structure of the Universe codified this framework, deriving the correlation-function hierarchy, the growth of perturbations, and the power-law clustering $\xi\propto r^{-1.8}$ that early surveys were already measuring. An early advocate of cold dark matter and of a cosmological constant when both were unfashionable, he supplied the perturbation theory that lets a measured $P(k)$ be read as a statement about $\Omega_m$ , the primordial spectrum, and the transfer function. When the 2005 SDSS detection of the baryon acoustic peak turned his abstract $\xi(r)$ into a precision ruler, it vindicated a program four decades old. Peebles received the 2019 Nobel Prize in Physics "for theoretical discoveries in physical cosmology."

History

1980

— Peebles publishes The Large-Scale Structure of the Universe, codifying the correlation-function formalism

1986

— de Lapparent, Geller and Huchra publish the first CfA slice, revealing walls and voids in the galaxy distribution

1987

— Kaiser derives the linear redshift-space distortion boost, turning a distortion into a growth measurement

2005

— Eisenstein and collaborators detect the baryon acoustic peak in SDSS luminous red galaxies

2021

— DESI begins its five-year campaign to map tens of millions of galaxy and quasar redshifts

Going further

Texts. Dodelson, Modern Cosmology, 2nd ed., Ch. 9 (Power Spectra) and Ch. 11 (Probes of Structure) · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 27 (Large-Scale Structure) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 6 (clustering statistics) · Peebles, The Large-Scale Structure of the Universe (1980), the founding correlation-function text · Binney & Tremaine, Galactic Dynamics, 2nd ed., §1.1 · Ryden, Introduction to Cosmology, 2nd ed., Ch. 12. The angular integral over the 3D Fourier transform yields $j_0(kr)=\sin(kr)/(kr)$ , giving the Hankel relation between $\xi(r)$ and $P(k)$ ; the Kaiser multipole coefficients are $(1+\tfrac23\beta+\tfrac15\beta^2)$ for the monopole and $(\tfrac43\beta+\tfrac47\beta^2)$ for the quadrupole. The BAO numbers ( $D_V=1379\ \mathrm{Mpc}$ , $\Delta z=0.040$ , $\Delta\theta=5.9^\circ$ at $z=0.35$ ) use flat $\Lambda$ CDM ( $\Omega_m=0.31$ , $H_0=67.4$ ); the growth values ( $f=0.76$ , $\beta=0.38$ , $f\sigma_8\approx0.46$ at $z=0.5$ ) use $f\approx\Omega_m(a)^{0.55}$ (Linder 2005). The standard anchors are $r_{\mathrm{drag}}\approx147\ \mathrm{Mpc}$ (Planck 2018), $r_0\approx5\,h^{-1}\mathrm{Mpc}$ , and $\gamma\approx1.8$ .

Prerequisites. Density Perturbations and the Density Contrast · Transfer Function and the Matter Power Spectrum · Linear Growth of Structure · Baryon Acoustic Oscillations. Leads to. The Sound Horizon and the Standard Ruler · The ΛCDM Concordance Model · Dark Energy and the Cosmological Constant.

Ask the tutor: Want to walk through how DESI turns a fibre spectrum into a point on the $f\sigma_8(z)$ curve, or work a BAO distance at a redshift of your choice?

Why measure a correlation function at all, instead of just looking at the map?

Because individual walls and voids are not predictable, but their statistics are. The initial conditions are a random field, so theory does not forecast that a particular filament sits at a particular place — only statistical measures like $\xi(r)$ and $P(k)$ can be compared between theory and data. Averaging also beats down noise: any single structure is one draw from the distribution, whereas $\xi(r)$ combines every galaxy pair at a given separation into a clean signal. The map is useful for intuition, but the science lives in the moments of the field.

What is galaxy bias, and why does it complicate everything?

Galaxies form in the peaks of the dark-matter field, so they over-represent overdense regions: $\delta_{\mathrm{gal}}=b\,\delta_{\mathrm{m}}$ with $b$ typically $1$ – $2.5$ for survey samples. The trouble is that $b$ is degenerate with the clustering amplitude — a survey measures $b^2\sigma_8^2$ , so a highly biased population with low $\sigma_8$ looks identical to a weakly biased one with high $\sigma_8$ . That is why surveys report the BAO scale (immune to a constant $b$ ) and $f\sigma_8$ (where the anisotropy ratio cancels $b$ ), and cross-correlate with gravitational lensing, which responds to mass directly, to pin $b$ down.

How can peculiar velocities both distort the map and help us?

They do both at once. Because a galaxy's redshift includes its line-of-sight peculiar velocity, plotting by redshift misplaces galaxies radially and deforms structures — that is the distortion. But the deformation is not arbitrary: in linear theory the infall velocity is sourced by the same growing density mode, with a proportionality equal to the growth rate $f$ . So the anisotropy that spoils a naive distance map is a direct readout of how fast structure is growing — one of the few clean, geometry-independent tests of gravity on cosmological scales.

Why is the BAO feature such a trustworthy standard ruler?

Because its length is fixed by simple early-universe physics and it sits at a scale where astrophysics is quiet. The acoustic scale $r_{\mathrm{drag}}$ is the comoving distance sound travelled in the photon–baryon plasma before baryons decoupled, computed from linear theory and calibrated to sub-percent precision by the CMB. In the galaxy distribution it appears as a sharp bump near $100\,h^{-1}\mathrm{Mpc}$ , far larger than the scales where galaxy formation, nonlinear collapse, and bias operate — so those effects shift the broadband clustering but leave the position of the bump almost untouched. A known scale whose location is robust to systematics is exactly what a standard ruler requires.

Correlation function or power spectrum — which should I use?

They hold the same information as a Fourier pair, so the choice is practical. The power spectrum $P(k)$ separates scales cleanly and makes the BAO wiggles and equality turnover explicit, but the survey window convolves it, mixing modes. The correlation function $\xi(r)$ localizes the BAO as a single bump and treats the geometry additively rather than as a convolution, but its bins are correlated and it spreads each scale across a range of $r$ . Most modern analyses fit both as a consistency check; any discrepancy flags an unmodelled systematic.

How do surveys deal with the ragged shape of the observed sky region?

Through random catalogues. A survey never covers a clean box: it has holes around bright stars, a boundary set by the telescope footprint, and a radial selection that thins with distance. Analysts generate a dense synthetic catalogue of "random" points sharing the survey's angular mask and radial selection but carrying no clustering, then measure clustering as the excess of real pairs over random pairs — the Landy–Szalay estimator for $\xi(r)$ , the FKP-weighted field for $P(k)$ . The randoms encode the geometry, so what remains is the genuine cosmological signal, not the shape of the observation.