The Cosmic Web and Large-Scale Structure

Point a redshift survey at the sky, plot every galaxy at the distance its recession velocity implies, and a structure leaps out that no one designed: dense knots of clusters strung along luminous filaments, filaments hemming sheet-like walls, and walls wrapped around cavernous, near-empty voids that fill most of the volume. This is the cosmic web, and it spans scales from a few megaparsecs up to more than a hundred. Above roughly $100\,h^{-1}\,\mathrm{Mpc}$ the pattern washes out and the distribution becomes statistically homogeneous, restoring the Cosmological Principle that lets us write a smooth Friedmann background in the first place. The web is the present-day fossil of a nearly featureless beginning. Inflation seeds a Gaussian random field of density contrast $\delta(\mathbf{x})=\rho(\mathbf{x})/\bar\rho-1$ ; gravity then amplifies it. While $|\delta|\ll1$ each Fourier mode grows by the same factor $D(a)$ , so the shape of the matter power spectrum $P(k)$ is frozen and only its amplitude climbs. When modes reach $\delta\sim1$ they collapse — but collapse follows the local tidal tensor, whose three eigenvalues are generically unequal, so matter drains along one axis first (a sheet), then a second (a filament), then the last (a knot). That ordering, made precise by the Zel'dovich approximation, is why the universe looks web-like rather than a scatter of round blobs. Cold dark matter builds the scaffolding, baryons fall into its wells and light up as galaxies, and dark energy freezes further growth once it dominates near $z\sim0.3$ . We read it statistically — through $\xi(r)$ , $P(k)$ , the baryon-acoustic bump near $150\,\mathrm{Mpc}$ , and the galaxy bias $b$ — which turns a spectacular picture into a precision cosmological instrument.

Key takeaways

The cosmic web is the frothy, filamentary arrangement of matter on scales of $\sim1$ – $100\,h^{-1}\,\mathrm{Mpc}$ : clusters at nodes, one-dimensional filaments and two-dimensional walls between them, and vast underdense voids filling most of the volume.
It is the gravitational amplification of tiny Gaussian primordial perturbations, whose statistics are carried entirely by the matter power spectrum $P(k)$ and its Fourier partner, the two-point correlation function $\xi(r)$ .
Anisotropic collapse — the Zel'dovich approximation — makes sheets first, then filaments, then knots, because the tidal deformation tensor generically has three distinct eigenvalues $\lambda_1>\lambda_2>\lambda_3$ .
In the linear era every mode grows by one universal factor $D(a)$ , so $P(k,a)=D^2(a)\,P(k)$ : the amplitude rises but the shape is preserved, which is why the web keeps its statistical pattern while it brightens.
The web is precision cosmology: the shape of $P(k)$ fixes $\Omega_m h$ , its amplitude fixes $\sigma_8\approx0.81$ , the BAO bump at $r\simeq150\,\mathrm{Mpc}$ is a standard ruler, and redshift-space distortions measure the growth rate $f\sigma_8$ .
Galaxies are biased tracers, $\xi_{\mathrm{gal}}\simeq b^2\,\xi_{\mathrm{m}}$ , so raw survey clustering overstates the underlying mass clustering and the bias must be modelled.
Key idea: gravitational instability amplifies a Gaussian primordial density field, and because collapse follows the anisotropic tidal tensor, matter drains sheets → filaments → knots, producing a web whose statistics ( $\xi$ , $P(k)$ , BAO) encode the cosmological parameters.

The core insight

Gravitational instability amplifies a Gaussian random primordial density field; because gravitational collapse follows the anisotropic tidal tensor, matter drains sheets → filaments → knots, producing a web whose statistics ( $\xi$ , $P(k)$ , BAO) encode the cosmological parameters.

Where this lives today

The cosmic web is now a working instrument of precision cosmology, not just a spectacular map. The Dark Energy Spectroscopic Instrument (DESI) is measuring tens of millions of galaxy and quasar redshifts to pin the baryon-acoustic ruler across $0<z<3.5$ , and its first results already sharpen — and mildly stir — the constraints on whether dark energy is a pure cosmological constant. Euclid and the Rubin Observatory's LSST add weak-lensing tomography, mapping the web's mass directly through the tiny shear it imprints on background galaxy shapes and combining clustering with lensing in " $3\times2$ pt" analyses that break the $\sigma_8$ – $\Omega_m$ degeneracy. Redshift-space distortions — Fingers-of-God stretching in clusters and Kaiser flattening of infalling sheets — become a measurement of the growth rate $f\sigma_8$ , testing whether cosmic-scale gravity really is General Relativity. On the theory side, billion-particle N-body and hydrodynamic simulations (Millennium, IllustrisTNG) grow synthetic webs matched against the surveys, while eigenvalue classifiers (the tidal T-web, the velocity-shear V-web) label every point as void, sheet, filament, or knot to study how environment shapes a galaxy's color, spin, and star formation. The sharpest open problem is the mild but persistent $S_8=\sigma_8\sqrt{\Omega_m/0.3}$ tension: lensing surveys tend to find the low-redshift web slightly smoother than the value extrapolated from the Planck CMB, and whether that gap is systematics, new physics, or a fluke is unresolved.

The physical picture

Begin with an almost-smooth fluid carrying tiny ripples in density. Gravity is a runaway amplifier: a slightly overdense patch pulls in more matter, deepens its well, and pulls in still more. But a real patch is never a perfect sphere — it is a squashed ellipsoid, and it collapses first along its shortest axis, flattening into a pancake or wall; then along the next, pinching into a filament; and only last along the third, drawing up into a knot where clusters form. Run that everywhere at once and the sheets, filaments, and knots stitch themselves into one connected network. The voids are not passive gaps but active drivers: underdense regions expand faster than the mean and sweep matter outward onto their shared boundaries. Where two voids' shells meet you get a wall; where three walls meet, a filament; where filaments cross, a cluster. The web is the drainage pattern of matter flowing out of expanding voids.

Figure

A schematic slice of the cosmic web: expanding voids drain matter onto the walls and filaments between them, clusters sit at the nodes where filaments cross, and above about 100 h⁻¹ Mpc the pattern averages out to homogeneity

Everything about the early web is Gaussian, and a Gaussian random field is fully specified by its two-point statistics — the power spectrum $P(k)$ , or equivalently the correlation function $\xi(r)$ . The abundance of peaks, the coherence length of clustering, and the frozen acoustic wiggle all follow from $P(k)$ and the growth factor; nonlinear evolution then adds the non-Gaussianity that higher-order statistics capture. The most useful single summary is the correlation function — the excess probability, over random, of finding a second galaxy a distance $r$ from a first. The next figure is that curve, where the abstract statistics turn into real numbers a survey can measure.

Figure

The galaxy two-point correlation function ξ(r): a power law ξ = (r/r₀)^−γ with correlation length r₀ ≈ 5 h⁻¹ Mpc and slope γ ≈ 1.8, the baryon-acoustic bump near r ≈ 100 h⁻¹ Mpc (≈ 150 Mpc), and ξ → 0 above the homogeneity scale

Scale sets the regime. Below the nonlinear scale — about $8\,h^{-1}\,\mathrm{Mpc}$ today, where the rms fluctuation $\sigma_8\approx0.81$ reaches order unity — structure has collapsed and virialized into halos; above it, perturbations are still linear and the web is a low-contrast lattice preserving its initial phases. The dividing scale grows with time as more of the universe crosses into the nonlinear regime.

The core equations

The starting variable is the density contrast, the dimensionless deviation of the local matter density from the cosmic mean. Everything downstream — growth, collapse, the correlation function — is written in terms of this one field.

\delta(\mathbf{x}) \equiv \frac{\rho(\mathbf{x}) - \bar\rho}{\bar\rho} \qquad (1)

where: $\rho(\mathbf{x})$ — local matter density; $\bar\rho$ — mean cosmic matter density; $\delta$ — dimensionless density contrast, bounded below by $-1$ (total void) and unbounded above. A region with $\delta>0$ is overdense and will deepen; a region with $\delta<0$ is a void that will empty further. Because the field is random, we characterize it statistically. The two-point correlation function is the excess probability of a pair at separation $r$ ; for a homogeneous field it is the Fourier transform of the power spectrum, so $\xi$ and $P$ carry identical information in real and reciprocal space.

\xi(r) \equiv \langle \delta(\mathbf{x})\,\delta(\mathbf{x}+\mathbf{r}) \rangle = \frac{1}{(2\pi)^3}\int P(k)\,e^{i\mathbf{k}\cdot\mathbf{r}}\,\mathrm{d}^3k \qquad (2)

where: $\xi(r)$ — excess pair probability over random at separation $r$ ; $P(k)$ — matter power spectrum, the variance of Fourier mode $\mathbf{k}$ ; $\langle\cdot\rangle$ — ensemble average over realizations of the field. A positive $\xi$ means clustering (pairs closer than random), and $\xi(r)\to0$ marks the scale of statistical homogeneity. Isotropy lets the angular integral in equation (2) be done once and for all, replacing the plane wave by a spherical Bessel function $j_0(kr)=\sin(kr)/(kr)$ . This is the form a survey pipeline actually evaluates.

\xi(r) = \frac{1}{2\pi^2}\int_0^\infty P(k)\,\frac{\sin(kr)}{kr}\,k^2\,\mathrm{d}k \qquad (3)

where: $\sin(kr)/(kr)$ — $j_0(kr)$ , from averaging $e^{i\mathbf{k}\cdot\mathbf{r}}$ over directions; $k^2\,\mathrm{d}k$ — the radial part of $\mathrm{d}^3k$ . Setting $r=0$ recovers the total variance $\sigma^2=\langle\delta^2\rangle$ . To ask how much structure lives at scale $k$ , use the variance per logarithmic interval — the dimensionless power spectrum — and smooth the field on a physical radius $R$ with a top-hat window. Smoothing on $R=8\,h^{-1}\,\mathrm{Mpc}$ defines the standard amplitude $\sigma_8$ .

\Delta^2(k)\equiv\frac{k^3P(k)}{2\pi^2},\qquad \sigma^2(R)=\frac{1}{2\pi^2}\int_0^\infty P(k)\,W^2(kR)\,k^2\,\mathrm{d}k,\quad W(x)=\frac{3(\sin x-x\cos x)}{x^3} \qquad (4)

where: $\Delta^2(k)$ — variance of $\delta$ per unit $\ln k$ ; $R$ — top-hat smoothing radius; $W(kR)$ — Fourier transform of a spherical top-hat; $\sigma_8\approx0.81$ — the rms mass fluctuation in $8\,h^{-1}\,\mathrm{Mpc}$ spheres. The condition $\Delta^2(k)\sim1$ marks the scale at which clustering becomes nonlinear. The observed shape of $P(k)$ is assembled from three physically distinct pieces: a nearly scale-invariant primordial spectrum, a transfer function that bends it at sub-horizon scales, and the growth factor that sets its present amplitude.

P(k) = A\,k^{n_s}\,T^2(k)\,D^2(a) \qquad (5)

where: $A$ — normalization, fixed by $\sigma_8$ or the CMB amplitude $A_s$ ; $n_s\approx0.965$ — primordial spectral index; $T(k)$ — transfer function encoding horizon-crossing and radiation-era physics; $D(a)$ — linear growth factor. The turnover of $T(k)$ near the horizon scale at matter–radiation equality fixes the combination $\Omega_m h$ , one of the web's headline outputs. Galaxies do not trace mass perfectly: they form in the peaks of the field, so their clustering is amplified by a bias factor. On large, linear scales the amplification is a simple squared multiplier.

\xi_{\mathrm{gal}}(r) \simeq b^2\,\xi_{\mathrm{m}}(r) \qquad (6)

where: $b$ — linear bias factor, $b>1$ for massive galaxies and clusters; $\xi_{\mathrm{gal}},\xi_{\mathrm{m}}$ — galaxy and matter correlation functions. Since the observable is $\xi_{\mathrm{gal}}$ , the mass clustering $\xi_{\mathrm{m}}$ is recoverable only once $b$ is known. To describe the collapse itself, the Zel'dovich approximation trades the Eulerian density for a Lagrangian displacement: every fluid element is pushed off its initial grid position $\mathbf{q}$ by a displacement proportional to the growth factor.

\mathbf{x}(\mathbf{q},t) = \mathbf{q} - D(t)\,\nabla_{\mathbf{q}}\Psi(\mathbf{q}) \qquad (7)

where: $\mathbf{q}$ — initial (Lagrangian) comoving position; $\mathbf{x}$ — present (Eulerian) comoving position; $\Psi$ — displacement potential with $\nabla^2\Psi=\delta_0$ ; $D(t)$ — growth factor setting the displacement amplitude. This single kinematic rule, extrapolated from linear theory, already produces sheets, filaments, and knots. Mass conservation turns the displacement map into a density: the density is the inverse Jacobian of the map, which factorizes into the three eigenvalues of the deformation tensor.

1+\delta(\mathbf{x},t) = \prod_{i=1}^{3}\frac{1}{1 - D(t)\,\lambda_i(\mathbf{q})} \qquad (8)

where: $\lambda_i$ — eigenvalues of the deformation tensor $\partial_i\partial_j\Psi$ , ordered $\lambda_1\ge\lambda_2\ge\lambda_3$ . The density diverges when $D(t)\,\lambda_1=1$ : the largest eigenvalue collapses first, and because the three are generically unequal the collapse is a sheet, not a sphere.

Derivation ALinear growth preserves the shape of the power spectrum

The goal is to show why the web keeps the same statistical pattern while it brightens: in the linear regime the equation of motion for each Fourier mode carries no dependence on the wavenumber, so every mode is amplified by one and the same time-dependent factor. That single fact means the power spectrum grows purely in amplitude, as $D^2(t)$ , while its dependence on scale — the shape encoding the primordial spectrum and transfer function — is frozen. We start from the linearized fluid equations for pressureless self-gravitating matter in an expanding background and follow the mode $\delta_{\mathbf{k}}$ .

Step 1. For pressureless (cold) matter the linearized equations combine into $\ddot\delta_{\mathbf{k}}+2H\dot\delta_{\mathbf{k}}-4\pi G\bar\rho\,\delta_{\mathbf{k}}=0$ , where the Hubble term $2H\dot\delta$ is the drag from cosmic expansion. Why: the would-be Jeans pressure term $c_s^2k^2/a^2$ vanishes for cold matter, and with it goes the only place $k$ could enter the equation.

Step 2. Because no $k$ appears, the time dependence separates from the spatial pattern: $\delta_{\mathbf{k}}(t)=[D(t)/D(t_i)]\,\delta_{\mathbf{k}}(t_i)$ , with the growing mode $D(t)$ obeying $\ddot D+2H\dot D-4\pi G\bar\rho\,D=0$ . Why: a $k$ -independent linear operator cannot mix scales, so it acts as a universal multiplier common to every mode.

Step 3. The power spectrum is the variance of the modes, $P(k,t)=\langle|\delta_{\mathbf{k}}(t)|^2\rangle=[D(t)/D(t_i)]^2\,P(k,t_i)$ , so the amplitude scales as $D^2$ while the shape $P(k,t_i)$ is untouched. Why: squaring the universal multiplier gives a single overall factor that pulls out of the ensemble average, leaving the scale dependence exactly as it was.

Step 4. In a matter-dominated Einstein–de Sitter universe the growing mode is $D(t)\propto a(t)\propto t^{2/3}$ , so $\delta\propto a$ ; once dark energy dominates, $D$ grows more slowly than $a$ and structure formation stalls. Why: this is the standard normalization used to extrapolate linear theory forward, and it fixes the "clock" $D$ that appears in the Zel'dovich map.

Go deeper · Derivation B · 1 minA power-law spectrum makes a power-law correlation function

Derivation BA power-law spectrum makes a power-law correlation function

The aim is to explain why the observed galaxy correlation function is so nearly a straight line on a log–log plot, $\xi(r)\propto r^{-1.8}$ , by showing that a scale-free power spectrum $P(k)\propto k^n$ transforms into a scale-free $\xi(r)\propto r^{-(n+3)}$ . The demonstration is a change of variables inside the isotropic transform of equation (3): once the power law is inserted, all $r$ -dependence factors out, leaving a pure number times a fixed power of $r$ — and the exponent connects a $k$ -space spectral index to a real-space clustering slope.

Step 1. Insert $P(k)=A\,k^n$ into the isotropic transform, giving $\xi(r)=\dfrac{A}{2\pi^2}\displaystyle\int_0^\infty k^{n}\,\dfrac{\sin(kr)}{kr}\,k^2\,\mathrm{d}k$ . Why: a scale-free spectrum is the one input for which the transform can be done without knowing any physical scale, isolating the pure geometry.

Step 2. Substitute the dimensionless variable $u=kr$ , so $k=u/r$ and $\mathrm{d}k=\mathrm{d}u/r$ , which pulls the whole $r$ -dependence outside: $\xi(r)=\dfrac{A}{2\pi^2\,r^{n+3}}\displaystyle\int_0^\infty u^{n+1}\sin u\,\mathrm{d}u$ . Why: every power of $k$ and the volume element together contribute $r^{-(n+3)}$ , and the remaining $u$ -integral is a pure number, convergent for $-3<n<0$ in the distributional sense.

Step 3. Reading off the exponent gives $\xi(r)\propto r^{-(n+3)}\equiv(r/r_0)^{-\gamma}$ with $\gamma=n+3$ . Observations give $\gamma\approx1.8$ , hence an effective $n\approx-1.2$ near the nonlinear scale. Why: this is the quantitative bridge between the theorist's spectral index and the observer's clustering slope, and it explains the famous $\xi\propto r^{-1.8}$ .

Step 4. The normalization is fixed by the correlation length, defined by $\xi(r_0)=1$ , which for galaxies is $r_0\approx5\,h^{-1}\,\mathrm{Mpc}$ . Why: $r_0$ marks where clustering falls to order unity — the coherence scale of the web's dense parts — and separates the strongly clustered interior from the near-random large scales.

Go deeper · Derivation C · 1 minZel'dovich collapse: why pancakes form first

Derivation CZel'dovich collapse: why pancakes form first

The purpose is to prove that the first structure to appear at any point is a two-dimensional sheet — a Zel'dovich pancake — rather than a spherical clump, and that this follows inevitably from the anisotropy of the tidal field. We take the density as the inverse Jacobian of the Lagrangian-to-Eulerian map, diagonalize the deformation tensor so the density becomes a product over its three eigenvalues, and ask when that product first diverges. Because a Gaussian field almost never gives three equal eigenvalues, the generic outcome is sequential collapse: sheet, then filament, then knot.

Step 1. Mass conservation writes the density as the inverse Jacobian of the map, $1+\delta=|\det(\partial\mathbf{x}/\partial\mathbf{q})|^{-1}$ , and from the Zel'dovich map $\partial x_i/\partial q_j=\delta_{ij}-D(t)\,\partial_i\partial_j\Psi$ . Why: the amount of mass is fixed in Lagrangian coordinates, so compressing the coordinate volume by the Jacobian is exactly what raises the density.

Step 2. Diagonalizing $\partial_i\partial_j\Psi$ in its principal axes gives $1+\delta=\prod_{i=1}^3(1-D\lambda_i)^{-1}$ with eigenvalues $\lambda_1\ge\lambda_2\ge\lambda_3$ . Why: in the principal frame the deformation is a product of three independent one-dimensional stretches, so the density factorizes cleanly.

Step 3. As $D(t)$ grows, the factor $1-D\lambda_1$ reaches zero first because $\lambda_1$ is the largest (and positive in an overdensity), so the density diverges when $D(t)\,\lambda_1=1$ . Only one factor blows up, meaning matter has collapsed along a single axis into a two-dimensional sheet; later $\lambda_2$ and then $\lambda_3$ collapse, giving a filament and then a knot. Why: a caustic in one direction is a pancake by definition, and the strict ordering of eigenvalues fixes the sheet → filament → knot sequence.

Step 4. For a Gaussian field the joint distribution of $(\lambda_1,\lambda_2,\lambda_3)$ — computed by Doroshkevich in 1970 — assigns essentially zero probability to three equal eigenvalues, so isotropic spherical collapse is a measure-zero special case and the web, not a field of round blobs, is the expected morphology. Why: this is what makes filamentarity generic rather than accidental, turning a statement about one point into a statement about the whole universe.

Limiting cases

The web wears a different face in each regime of density contrast and epoch; the table collects the four that matter.

Case	Result	Meaning
Linear regime	$\vert\delta\vert\ll1,\ \ \delta(\mathbf{x},t)=D(t)\,\delta_0(\mathbf{x})$	Every mode grows self-similarly and $P(k,t)=D^2(t)P(k)$ ; the web is a faint, symmetric lattice preserving its initial phases.
Mildly nonlinear (Zel'dovich)	$\mathbf{x}=\mathbf{q}-D(t)\,\nabla_q\Psi(\mathbf{q})$	Anisotropic collapse; sheets form first at caustics where the map turns multivalued. Excellent until the first shell crossing.
Strongly nonlinear	$\delta\gtrsim1,\ \ r\lesssim8\,h^{-1}\mathrm{Mpc}$	Virialized halos and clusters; needs N-body simulation. $\xi(r)\simeq(r/r_0)^{-\gamma}$ with $r_0\approx5\,h^{-1}\mathrm{Mpc}$ , $\gamma\approx1.8$ .
Homogeneity scale	$r\gtrsim100\,h^{-1}\mathrm{Mpc}\ \Rightarrow\ \xi(r)\to0$	The distribution is statistically homogeneous and isotropic, recovering the Cosmological Principle that the smooth Friedmann background assumes.

Common trapThe correlation function describes statistics, not shells around galaxies

Two traps hide in the same curve. First, $\xi(r)$ and $P(k)$ are statistical descriptors of the whole density field, not maps of any one object: the baryon-acoustic bump near $150\,\mathrm{Mpc}$ is a slight excess in the ensemble pair count at that separation, not a physical ring drawn around each galaxy. Reading a peak in $\xi(r)$ as a structure around individual galaxies is a category error. Second, galaxies are biased tracers, so the measured $\xi_{\mathrm{gal}}$ is roughly $b^2$ times the matter correlation $\xi_{\mathrm{m}}$ ; with a typical $b\sim1.5$ the raw survey clustering overstates the underlying mass fluctuations by more than a factor of two. Any parameter — $\sigma_8$ especially — read straight off galaxy clustering without dividing out $b^2$ will be biased high. The bias must be measured (from redshift-space distortions, lensing, or higher-order statistics), not assumed away.

Worked examples

1 · From galaxy clustering to matter clustering (fully worked)

A survey measures $\xi_{\mathrm{gal}}(r)=(r/5\,h^{-1}\mathrm{Mpc})^{-1.8}$ for a population with linear bias $b=1.5$ . What is the matter correlation function at $r=10\,h^{-1}\mathrm{Mpc}$ , and what does the answer say about how much galaxies exaggerate clustering?

Step 1. Evaluate the galaxy correlation at the requested separation: $\xi_{\mathrm{gal}}(10)=(10/5)^{-1.8}=2^{-1.8}$ . Since $2^{1.8}=3.48$ , this is $\xi_{\mathrm{gal}}(10)=\mathbf{0.287}$ — the pair count there is about 29% above random.

Step 2. Divide out the bias using equation (6), $\xi_{\mathrm{m}}=\xi_{\mathrm{gal}}/b^2=0.287/1.5^2=0.287/2.25$ . The matter correlation is $\xi_{\mathrm{m}}(10)=\mathbf{0.128}$ .

Step 3. Interpret it: the underlying mass is clustered less than the galaxies by exactly the factor $b^2=2.25$ at every separation on these linear scales. So a naive reader who took galaxy clustering for mass clustering would overstate $\xi$ by a factor of $2.25$ — and because $\sigma\propto\sqrt{\xi}$ in the same crude sense, would overstate the amplitude $\sigma_8$ by about $b=1.5$ . $\xi_{\mathrm{m}}(10)=0.128$ , i.e. $b^2=2.25$ below the galaxy value.

2 · When does a Zel'dovich pancake collapse? (fully worked)

A protostructure has a largest tidal eigenvalue whose linear-theory value today (with the growth factor normalized to $D_0=1$ at $a=1$ ) is $D_0\lambda_1=2$ . In an Einstein–de Sitter universe, where $D\propto a$ , at what scale factor and redshift does its first caustic form?

Step 1. The first caustic condition from equation (8) is $D(a)\,\lambda_1=1$ . With $D\propto a$ normalized so $D_0=1$ , write $D(a)=a$ , and $\lambda_1=D_0\lambda_1/D_0=2$ in the same units.

Step 2. Solve $a\,\lambda_1=1\Rightarrow a=1/\lambda_1=1/2$ . The sheet forms at scale factor $\mathbf{a=0.5}$ .

Step 3. Convert to redshift, $1+z=1/a=2$ , so $\mathbf{z=1}$ . A region whose linearly-extrapolated eigenvalue reaches $2$ by today actually crossed its first caustic long ago, at $z=1$ — the sheet is already in place while linear theory, blind to the collapse, keeps reporting a finite $\delta$ . This is exactly why the web is fully assembled by $z\sim1$ even though linear extrapolation of $\delta$ never diverges. First pancake at $a=0.5$ , $z=1$ .

3 · Reading a spectral index off a clustering slope

A galaxy sample shows $\xi(r)\propto r^{-1.8}$ near $r_0$ . What effective power-spectrum index $n$ does it have there, and how does that compare with the primordial slope?

Step 1. Use the scale-free result $\gamma=n+3$ from the second derivation, with $\gamma=1.8$ : $n=\gamma-3=\mathbf{-1.2}$ .

Step 2. Compare with the primordial spectrum $P(k)\propto k^{n_s}$ , $n_s\approx0.965$ . On the largest scales the spectrum rises with $k$ ; near the nonlinear scale the transfer function and nonlinear processing have driven the effective index down to $-1.2$ , far redder. Effective $n\approx-1.2$ versus primordial $n_s\approx0.965$ — the spectrum reddens toward small scales.

Physics of the portrait

The cosmic web has many parents, but one figure turned it from a picture into a prediction. Before 1970 the growth of structure was understood only in the linear, spherical idealization; the insight that collapse is generically flat — a pancake — and that this single idea organizes the entire large-scale morphology belongs to Yakov Zel'dovich.

Origins · 1 minYakov Zel'dovich and how this was discovered

Yakov Zel'dovich

1914–1987 · gravitational instability and the pancake theory of cosmic structure

Yakov Borisovich Zel'dovich was a self-taught polymath who moved from combustion and detonation theory through the Soviet nuclear program to become, in his final decades, one of the founders of physical cosmology. His 1970 paper introduced what is now called the Zel'dovich approximation: instead of tracking the Eulerian density, follow each fluid element as it is displaced from its initial position by an amount proportional to the growth factor, $\mathbf{x}=\mathbf{q}-D(t)\nabla_q\Psi$ . Extrapolated straight out of linear theory, this deceptively simple rule captures the essential nonlinearity of collapse — the formation of caustics — and predicts that matter first flattens into two-dimensional sheets, the "Zel'dovich pancakes." When redshift surveys in the 1980s revealed sheets, filaments, and voids rather than isolated clumps, they vindicated the picture. With collaborators he built much of the modern theory of structure formation: the statistics of peaks in Gaussian fields, the eigenvalue statistics of the tidal tensor (with Doroshkevich), the Sunyaev–Zel'dovich effect of clusters on the microwave background, and early work on baryon acoustic oscillations — a range that led Hawking to quip he had once assumed "Zel'dovich" was a collective pen-name like Bourbaki.

History

1922

— Friedmann derives the expanding-universe solutions that fix the smooth background the web perturbs

1970

— Zel'dovich publishes the approximation predicting sheet-first (pancake) collapse

1970

— Doroshkevich computes the eigenvalue statistics of the Gaussian tidal tensor

1980

— Peebles' "The Large-Scale Structure of the Universe" codifies the ξ(r) description

1986

— the CfA redshift survey reveals the Great Wall and enormous voids

2005

— SDSS and 2dFGRS detect the baryon-acoustic peak in the galaxy correlation function

Going further

Texts. Ryden, Introduction to Cosmology, 2nd ed., Ch. 11–12 · Dodelson & Schmidt, Modern Cosmology, 2nd ed., Ch. 8, 9, 11 · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 30 · Binney & Tremaine, Galactic Dynamics, 2nd ed., Ch. 9 · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 4, 6, 7 · Peebles, The Large-Scale Structure of the Universe (1980) · Zel'dovich, A&A 5, 84 (1970). The key relations used throughout are: the Fourier-pair relation between $\xi(r)$ and $P(k)$ ; the scale-free result $\xi\propto r^{-(n+3)}$ ; the Zel'dovich density $1+\delta=\prod(1-D\lambda_i)^{-1}$ with first-caustic condition $D\lambda_1=1$ ; and the $k$ -independence of the linear growth equation giving $P\propto D^2$ . Standard numbers quoted: $\sigma_8\approx0.81$ , $n_s\approx0.965$ , $r_0\approx5\,h^{-1}\mathrm{Mpc}$ , $\gamma\approx1.8$ , BAO scale $\approx150\,\mathrm{Mpc}\ (\approx100\,h^{-1}\mathrm{Mpc})$ , homogeneity scale $\approx100\,h^{-1}\mathrm{Mpc}$ .

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

Ask the tutor: Ask me to walk through why a squashed ellipsoid of dark matter always collapses into a pancake first, or to turn a measured $\xi(r)$ into $\sigma_8$ , $\Omega_m$ , and the bias $b$ .

Is the cosmic web made of dark matter or of galaxies?

Both, in layers. The load-bearing structure is cold dark matter: because it is collisionless and pressureless, its perturbations grow on all scales and collapse first, forming the filaments and halos whose gravity defines the web's shape. Baryons then fall into those dark-matter potential wells, cool, and condense into the galaxies we actually see, so galaxies light up the densest parts of a scaffold they did not build. Diffuse gas — the warm-hot intergalactic medium — traces the filaments between galaxies and is only now being detected directly. So when a redshift survey shows the web, it is showing biased galaxy tracers sitting on an underlying dark-matter network that carries most of the mass.

Why do sheets form before filaments and knots, not all at once?

Because gravitational collapse is set by the tidal deformation tensor, and that tensor almost never has three equal eigenvalues. In the principal axes the density is $1+\delta=\prod_i(1-D\lambda_i)^{-1}$ , and as the growth factor $D$ rises the term with the largest eigenvalue $\lambda_1$ hits its caustic first. When it does, matter has collapsed along one axis only — a two-dimensional sheet. The next eigenvalue collapses later to make a filament, and the smallest last to make a knot. Spherical collapse, with all three eigenvalues equal, is a measure-zero coincidence for a Gaussian field, so the generic, expected outcome is the ordered sheet → filament → knot sequence that weaves the web.

What is the difference between the correlation function and the power spectrum?

They are the same information in two languages. The correlation function $\xi(r)$ lives in real space and answers "how much more likely than random is a galaxy pair separated by $r$ ?" The power spectrum $P(k)$ lives in Fourier space and answers "how much variance sits in modes of wavenumber $k$ ?" For a statistically homogeneous field the two are an exact Fourier-transform pair, equation (2). Cosmologists use whichever is convenient: $P(k)$ separates scales cleanly and has nearly diagonal errors, favored for extracting parameters, while $\xi(r)$ isolates the localized baryon-acoustic bump and resists some systematics. For a Gaussian field, neither contains information the other lacks.

Why is 100 h⁻¹ Mpc special — what happens above it?

Above roughly $100\,h^{-1}\,\mathrm{Mpc}$ the correlation function $\xi(r)$ drops to essentially zero, meaning galaxy positions on those scales are statistically independent of one another. This is the homogeneity scale, and it matters because the whole framework of Friedmann cosmology rests on the Cosmological Principle — that the universe is homogeneous and isotropic when averaged on large enough scales. The web is the perturbation; the smooth expanding background is what survives the averaging. If clustering persisted to arbitrarily large scales the smooth-background description would fail and there would be no well-defined mean density to perturb around. Surveys confirm the transition to homogeneity near this scale, which is one of the deepest observational checks of the standard model.

Does dark energy change the shape of the cosmic web?

Not the shape so much as the timing. Dark energy became dynamically important around $z\sim0.3$ , and its accelerated expansion suppresses the growth of new structure by stretching space faster than gravity can pull matter together. The growth factor $D(a)$ then grows more slowly than the scale factor $a$ and eventually nearly freezes. So the web is essentially assembled by $z\sim1$ and changes little afterward — voids keep emptying and the largest filaments sharpen slightly, but no fundamentally new large-scale pattern forms. This freezing is itself a probe: comparing $f\sigma_8$ from redshift-space distortions across epochs constrains the dark-energy equation of state and tests General Relativity.

How do we know the initial field was Gaussian and nearly scale-invariant?

From the cosmic microwave background. The anisotropies Planck maps are a snapshot of the same primordial field at $z\approx1100$ , before nonlinear collapse scrambled it. Their statistics are Gaussian to exquisite precision — primordial non-Gaussianity is consistent with zero — and their spectrum is a nearly scale-invariant power law with $n_s\approx0.965$ , slightly red as simple inflation predicts. Because linear growth preserves both the shape of the spectrum and the Gaussianity of the field, the low-redshift web inherits both on large scales, and the match between the CMB-predicted $P(k)$ and the galaxy-measured $P(k)$ across a factor of a thousand in redshift is a triumph of the standard model. The non-Gaussianity of today's web is generated by nonlinear gravity, not present at the start.