The Lambda-CDM Concordance Model

The Lambda-CDM concordance model is the standard six-parameter theory of the whole observable universe, and it is remarkable for how little it asks and how much it delivers. It assumes only that on scales above ∼100\sim100 Mpc the cosmos is statistically homogeneous and isotropic — the cosmological principle — so that its geometry is the FLRW metric and its dynamics the Friedmann equations, sourced by a cosmic fluid of four ingredients: photons and neutrinos (radiation, w=1/3w=1/3), cold dark matter and baryons (pressureless matter, w=0w=0), and a cosmological constant Λ\Lambda (dark energy, w=−1w=-1). Feed those in and the entire expansion history collapses to one dimensionless function E(a)=H(a)/H0E(a)=H(a)/H_0, a sum of power-law terms that dominate in turn. The word concordance is the whole point. Utterly different probes — the acoustic peaks of the cosmic microwave background, the baryon-acoustic ruler in galaxy surveys, the standard-candle distances of Type Ia supernovae, and the light-element abundances of big-bang nucleosynthesis — each sensitive to different physics at different epochs, all converge on the same parameter values: a spatially flat universe that is today about 68%68\% dark energy, 27%27\% cold dark matter, and 5%5\% ordinary baryons, 13.813.8 billion years old. That convergence, not any one measurement, is the empirical triumph. Yet the model is frankly phenomenological. The measured Λ\Lambda is some 1012010^{120} times smaller than the naive quantum-vacuum estimate; dark matter has no confirmed particle; and percent-level data now expose real cracks — a ∼5σ\sim5\sigma tension in H0H_0, and a milder one in the clustering amplitude S8S_8. Lambda-CDM is thus at once the most successful and the most puzzling theory in physics.

Key takeaways
  • Lambda-CDM describes a flat, homogeneous, isotropic universe whose energy budget today is ∼68%\sim68\% cosmological constant Λ\Lambda (dark energy), ∼27%\sim27\% cold (non-relativistic, collisionless) dark matter, and ∼5%\sim5\% baryons, plus a whisker of radiation.
  • It rests on three earlier articles — the FLRW metric, the Friedmann equations, and the density parameters — and is fixed by six numbers: H0H_0, Ωbh2\Omega_b h^2, Ωch2\Omega_c h^2, ΩΛ\Omega_\Lambda (or curvature), the spectral index nsn_s, and the perturbation amplitude AsA_s.
  • The whole expansion history is one equation, H2(a)=H02 (Ωra−4+Ωma−3+Ωka−2+ΩΛ)H^2(a)=H_0^2\,(\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda), whose terms encode the successive radiation-, matter-, and Λ\Lambda-dominated eras.
  • One parameter set simultaneously fits the CMB acoustic peaks, baryon acoustic oscillations, supernova distances, nucleosynthesis abundances, and large-scale structure to percent precision — the concordance.
  • The open frontier is the unexplained nature of Λ\Lambda and dark matter, plus the H0H_0 and S8S_8 tensions that may signal physics beyond the standard model.
  • Key idea: the late-universe expansion history is encoded in the single function E(a)=H(a)/H0=Ωra−4+Ωma−3+Ωka−2+ΩΛE(a)=H(a)/H_0=\sqrt{\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda}; with flatness and the measured densities this one relation reproduces distance, growth, and age data across cosmic time.
The core insight

The entire late-universe expansion history is encoded in one dimensionless function E(a)≡H(a)/H0=Ωra−4+Ωma−3+Ωka−2+ΩΛE(a)\equiv H(a)/H_0=\sqrt{\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda}; with flatness (Ωk=0\Omega_k=0) and the measured density parameters this single equation reproduces redshift-distance, growth, and age data across cosmic history.

Where this lives today

Lambda-CDM is not a settled monument but a working model under intense stress-testing. The Dark Energy Spectroscopic Instrument (DESI), mapping tens of millions of galaxy and quasar redshifts, is measuring the baryon-acoustic ruler to sub-percent precision; its 2024–2025 data hinted that dark energy's equation of state may be evolving rather than exactly w=−1w=-1, which — if it survives — would crack the "constant" in cosmological constant. Euclid and the Vera Rubin Observatory's LSST are charting weak gravitational lensing to test whether structure grew as predicted, sharpening the S8S_8 tension. Most acutely, the H0H_0 tension persists: the SH0ES distance-ladder value (≈73\approx73) and the Planck CMB value (≈67.4\approx67.4) now disagree at roughly 5σ5\sigma, and JWST's independent calibrations have neither dissolved nor confirmed it. Meanwhile the particle identity of cold dark matter stays open, with direct-detection experiments (LZ, XENONnT) squeezing WIMP cross-sections ever tighter. The concordance still fits almost everything — but the frontier is exactly the handful of places where "almost" is starting to show.

The physical picture

Picture the contents of the universe as a portfolio of assets that decay at different rates as space expands. Radiation loses value fastest, ρ∝a−4\rho\propto a^{-4}, because it dilutes by volume and each photon's wavelength stretches; matter falls off more slowly, ρ∝a−3\rho\propto a^{-3}, from volume dilution alone; and dark energy, ρ∝a0\rho\propto a^{0}, never decays at all. Run the clock forward and whichever component decays slowest must eventually dominate — which guarantees a dark-energy future no matter how the universe began. The first figure plots exactly these three curves and marks where they cross.

Figure
Cosmic energy densities scale differently with the expansion: radiation ρ ∝ a⁻⁴, matter ρ ∝ a⁻³, and dark energy ρ = const. The curves cross so that radiation equals matter near z ≈ 3400 and matter equals Λ near z ≈ 0.30, giving three successive eras: radiation-, matter-, then Λ-dominated
Each component dilutes at its own rate — so the universe passes through three eras 10⁻⁴ 10⁻³ 10⁻² 10⁻¹ 1 scale factor a (log, → cosmic time) energy density ρ / ρ_c,0 (log) today a = 1 ρ ∝ a⁻⁴ radiation ρ ∝ a⁻³ matter ρ = const dark energy Λ radiation = matter, z ≈ 3400 matter = Λ, z ≈ 0.30

The reason this matters is that the ordering of the eras is fixed by the exponents alone: because 4>3>04>3>0, radiation always wins early, matter takes over in the middle, and Λ\Lambda inherits the far future. The concordance model's second face is not a timeline but a statement about agreement — and the cleanest way to see it is a parameter plane. Each observational probe carves out a degeneracy band in the (Ωm,ΩΛ)(\Omega_m,\Omega_\Lambda) plane, and because the bands point in nearly orthogonal directions they pin down a tiny common region. The second figure shows the classic picture.

Figure
The concordance in the density plane: the cosmic microwave background constrains the total density near the flat line Ω_m + Ω_Λ = 1, Type Ia supernovae constrain the accelerating combination Ω_Λ − Ω_m, and baryon acoustic oscillations pin Ω_m ≈ 0.3. The three near-orthogonal bands intersect at the concordance point Ω_m = 0.315, Ω_Λ = 0.685
Three probes, one answer: near-orthogonal bands cross at Ω_m = 0.315, Ω_Λ = 0.685 0 0.5 1.0 0 0.5 1.0 matter density Ω_m dark energy Ω_Λ flat: Ω_m + Ω_Λ = 1 CMB SNe Ia BAO (0.315, 0.685)

Finally, why cold dark matter? Cold means non-relativistic at the epoch of structure formation, so the dark matter has negligible velocity dispersion and does not free-stream out of small overdensities. That lets perturbations on all scales grow gravitationally and build structure hierarchically — small halos first, merging into larger ones — matching the observed galaxy and cluster distribution far better than hot (relativistic) dark matter, whose free-streaming would erase small-scale structure. Cold dark matter turns the smooth early universe into the cosmic web.

The core equations

The geometry is set once and for all by the cosmological principle: the most general homogeneous, isotropic spacetime is the Friedmann-Lemaitre-Robertson-Walker metric, whose only freedom is the scale factor a(t)a(t) and a constant spatial curvature kk. Everything dynamical rides on a(t)a(t).

ds2=−c2 dt2+a2(t)[dr21−kr2+r2 dΩ2](1)\mathrm{d}s^2=-c^2\,\mathrm{d}t^2+a^2(t)\left[\frac{\mathrm{d}r^2}{1-kr^2}+r^2\,\mathrm{d}\Omega^2\right] \qquad (1)

where: a(t)a(t) — scale factor, normalized a(t0)=1a(t_0)=1 today; kk — spatial curvature constant (−1,0,+1-1,0,+1 for open, flat, closed); r,Ωr,\Omega — comoving radial and angular coordinates; cc — speed of light. Reading it: proper distances between comoving observers all scale together as a(t)a(t), so the single number aa carries the entire expansion. The dynamics follow from inserting this metric into Einstein's equations with a perfect-fluid source. Dividing the first Friedmann equation by H02H_0^2 and using how each component dilutes gives the master relation — the workhorse of the whole model.

H2(a)=H02(Ωra−4+Ωma−3+Ωka−2+ΩΛ)(2)H^2(a)=H_0^2\left(\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda\right) \qquad (2)

where: H=a˙/aH=\dot a/a — Hubble parameter; H0H_0 — present value, ≈67\approx67–73 km s−1 Mpc−173\ \mathrm{km\,s^{-1}\,Mpc^{-1}}; Ωi\Omega_i — present density parameter of component ii. Reading it: H2H^2 is a sum of four power laws, each dominating in its era; the whole of observational cosmology is the project of measuring the Ωi\Omega_i. Each term's exponent is not arbitrary — it is fixed by the equation of state wi=pi/(ρic2)w_i=p_i/(\rho_i c^2) through the continuity equation, which tells each decoupled fluid how to dilute as space expands.

ρi(a)=ρi,0 a−3(1+wi),wr=13, wm=0, wΛ=−1(3)\rho_i(a)=\rho_{i,0}\,a^{-3(1+w_i)},\qquad w_r=\tfrac13,\ w_m=0,\ w_\Lambda=-1 \qquad (3)

where: ρi,0\rho_{i,0} — present density of component ii; wiw_i — its (constant) equation-of-state parameter. Reading it: w=1/3w=1/3 gives a−4a^{-4}, w=0w=0 gives a−3a^{-3}, w=−1w=-1 gives a0a^0 — the three slopes of the first figure, read straight off the exponent. The density parameters are ratios to a single benchmark: the critical density that makes the universe spatially flat. Its present value sets the scale for "how much stuff" the cosmos contains.

Ωi≡ρi,0ρc,ρc=3H028πG≈8.5×10−27 kg m−3 (≈5 protons m−3)(4)\Omega_i\equiv\frac{\rho_{i,0}}{\rho_c},\qquad \rho_c=\frac{3H_0^2}{8\pi G}\approx8.5\times10^{-27}\ \mathrm{kg\,m^{-3}}\ (\approx5\ \text{protons m}^{-3}) \qquad (4)

where: GG — Newton's constant; the numerical ρc\rho_c uses H0=67.4 km s−1 Mpc−1H_0=67.4\ \mathrm{km\,s^{-1}\,Mpc^{-1}}. Reading it: with the density parameters so defined, evaluating equation (2) at a=1a=1 forces the closure relation Ωr+Ωm+Ωk+ΩΛ=1\Omega_r+\Omega_m+\Omega_k+\Omega_\Lambda=1, and flat Lambda-CDM sets Ωk=0\Omega_k=0. Whether the universe accelerates is governed by the second Friedmann equation, which depends on the combination ρ+3p/c2\rho+3p/c^2 — pressure gravitates in general relativity, and negative pressure anti-gravitates.

a¨a=−4πG3(ρ+3pc2)+Λc23(5)\frac{\ddot a}{a}=-\frac{4\pi G}{3}\left(\rho+\frac{3p}{c^2}\right)+\frac{\Lambda c^2}{3} \qquad (5)

where: pp — total pressure, pi=wiρic2p_i=w_i\rho_i c^2; a¨\ddot a — cosmic acceleration. Reading it: a¨>0\ddot a>0 requires an effective w<−1/3w<-1/3, so ordinary matter and radiation decelerate the expansion while Λ\Lambda (w=−1w=-1, p=−ρc2p=-\rho c^2) drives it to accelerate once it dominates. To confront data, the expansion history is turned into distances. In a flat universe the comoving distance integrates 1/H1/H, and the observable luminosity and angular-diameter distances differ from it by factors of (1+z)(1+z).

dL(z)=(1+z) c∫0zdz′H(z′),dA(z)=dL(z)(1+z)2,1+z=1a(6)d_L(z)=(1+z)\,c\int_0^z\frac{\mathrm{d}z'}{H(z')},\qquad d_A(z)=\frac{d_L(z)}{(1+z)^2},\qquad 1+z=\frac{1}{a} \qquad (6)

where: zz — redshift; dLd_L — luminosity distance (Type Ia supernovae); dAd_A — angular-diameter distance (CMB, BAO standard rulers). Reading it: supernovae probe dL(z)d_L(z), the acoustic ruler probes dA(z)d_A(z), and both are integrals over the same H(z)H(z) — which is why they constrain the same Ωi\Omega_i. Structure is the last ingredient. On sub-horizon scales in the matter era, the fractional density contrast δ=δρ/ρ\delta=\delta\rho/\rho of cold matter obeys a driven-oscillator equation whose growing solution builds the cosmic web.

δ¨+2Hδ˙−4πGρm δ=0(7)\ddot\delta+2H\dot\delta-4\pi G\rho_m\,\delta=0 \qquad (7)

where: δ\delta — fractional matter overdensity; the 2Hδ˙2H\dot\delta term is Hubble drag; 4πGρm4\pi G\rho_m is self-gravity. Reading it: in the matter era this yields a growing mode δ∝a\delta\propto a on all scales, which only works because cold dark matter is pressureless and cannot resist collapse.

Derivation ABuilding the master Friedmann equation as a sum of power laws

The goal is to show that the single most useful equation in cosmology, H2=H02(Ωra−4+Ωma−3+Ωka−2+ΩΛ)H^2=H_0^2(\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda), is not an assumption but a bookkeeping consequence of two inputs owned by earlier articles: the first Friedmann equation, and the rule that each fluid dilutes as a power of the scale factor set by its equation of state. We assemble it in four moves and check that it closes to unity today.

Step 1. Start from the first Friedmann equation with curvature and cosmological constant, H2=8πG3ρ−kc2a2+Λc23H^2=\tfrac{8\pi G}{3}\rho-\tfrac{kc^2}{a^2}+\tfrac{\Lambda c^2}{3}, where ρ\rho is the total density of all matter and radiation. Why: this is the tttt Einstein equation for the FLRW metric — it is where geometry (HH, kk) is tied to content (ρ\rho, Λ\Lambda), and everything else is algebra on top of it.

Step 2. Split ρ\rho into decoupled components and use the continuity equation result ρi=ρi,0 a−3(1+wi)\rho_i=\rho_{i,0}\,a^{-3(1+w_i)}, so radiation (w=1/3w=1/3) scales as a−4a^{-4}, matter (w=0w=0) as a−3a^{-3}, and a constant-Λ\Lambda vacuum (w=−1w=-1) as a0a^{0}. Why: each component redshifts independently once it stops exchanging energy with the others, and the exponent is dictated entirely by ww — this is what makes the eras separate cleanly.

Step 3. Absorb the curvature and Λ\Lambda terms into the same language by defining Ωi=8πGρi,0/3H02\Omega_i=8\pi G\rho_{i,0}/3H_0^2, together with ΩΛ=Λc2/3H02\Omega_\Lambda=\Lambda c^2/3H_0^2 and Ωk=−kc2/H02\Omega_k=-kc^2/H_0^2, all normalized to the critical density ρc=3H02/8πG\rho_c=3H_0^2/8\pi G. Why: writing every term as a dimensionless fraction of the critical density lets curvature and dark energy sit in the same sum as real matter, and turns the equation into a statement about pure numbers.

Step 4. Substitute and factor out H02H_0^2 to obtain H2=H02(Ωra−4+Ωma−3+Ωka−2+ΩΛ)H^2=H_0^2(\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda); evaluating at a=1a=1 where H=H0H=H_0 forces Ωr+Ωm+Ωk+ΩΛ=1\Omega_r+\Omega_m+\Omega_k+\Omega_\Lambda=1. Why: the closure condition is not an extra assumption but an automatic identity of the construction, and setting Ωk=0\Omega_k=0 is exactly the flatness that the CMB measures to better than a percent.

Go deeper · Derivation B · 1 minThe onset of acceleration and the redshift of Lambda-domination
Derivation BThe onset of acceleration and the redshift of Lambda-domination

The goal is to pin down two distinct late-time milestones using nothing but the master equation and its acceleration counterpart: the redshift zmΛz_{m\Lambda} at which dark energy first out-densities matter, and the earlier redshift zaccz_\mathrm{acc} at which the expansion actually starts to accelerate. The two are not the same, and seeing why is the sharpest way to understand what Λ\Lambda does.

Step 1. For the equality of densities, set the matter and dark-energy terms of E2(a)E^2(a) equal, Ωma−3=ΩΛ\Omega_m a^{-3}=\Omega_\Lambda, neglecting the tiny radiation and curvature terms. Solving gives amΛ=(Ωm/ΩΛ)1/3a_{m\Lambda}=(\Omega_m/\Omega_\Lambda)^{1/3}, hence 1+zmΛ=(ΩΛ/Ωm)1/31+z_{m\Lambda}=(\Omega_\Lambda/\Omega_m)^{1/3}. Why: this is the moment the dark-energy floor rises above the diluting matter curve — the crossing marked in the first figure — after which Λ\Lambda dominates the budget forever.

Step 2. Acceleration is a different question: it is governed by a¨\ddot a, not by which density is larger. From the second Friedmann equation, in scaled units the matter term contributes −12Ωma−3-\tfrac12\Omega_m a^{-3} (because ρ+3p/c2=ρ\rho+3p/c^2=\rho for pressureless matter) while Λ\Lambda contributes +ΩΛ+\Omega_\Lambda. Why: pressure gravitates, so it is the combination ρ+3p/c2\rho+3p/c^2 that curves the expansion — and only Λ\Lambda, with its large negative pressure, can flip the sign of a¨\ddot a.

Step 3. Acceleration begins where a¨=0\ddot a=0, i.e. 12Ωma−3=ΩΛ\tfrac12\Omega_m a^{-3}=\Omega_\Lambda, giving aacc=(Ωm/2ΩΛ)1/3a_\mathrm{acc}=(\Omega_m/2\Omega_\Lambda)^{1/3} and 1+zacc=(2ΩΛ/Ωm)1/31+z_\mathrm{acc}=(2\Omega_\Lambda/\Omega_m)^{1/3}. The factor of 22 pushes this to a higher redshift than equality. Why: acceleration only needs Λ\Lambda to beat half the matter density (the effective gravitating matter is ρm\rho_m, but the threshold for a¨>0\ddot a>0 is set by ΩΛ>12Ωma−3\Omega_\Lambda>\tfrac12\Omega_m a^{-3}), so the cosmos starts speeding up before dark energy formally dominates.

Step 4. With Planck values Ωm=0.315\Omega_m=0.315, ΩΛ=0.685\Omega_\Lambda=0.685: equality at zmΛ≈0.30z_{m\Lambda}\approx0.30 and acceleration onset at zacc≈0.63z_\mathrm{acc}\approx0.63 — both firmly in the reach of supernova surveys, and consistent with the observed turnover from deceleration to acceleration. Why: putting real numbers on abstract thresholds is what let the 1998 supernova data confirm that we live after zaccz_\mathrm{acc}, in an accelerating universe. *Numerically: (ΩΛ/Ωm)1/3−1=0.296(\Omega_\Lambda/\Omega_m)^{1/3}-1=0.296 and (2ΩΛ/Ωm)1/3−1=0.632(2\Omega_\Lambda/\Omega_m)^{1/3}-1=0.632, confirming that acceleration precedes Lambda-domination, as it must.*

Go deeper · Derivation C · 1 minThe age of the flat Lambda-CDM universe in closed form
Derivation CThe age of the flat Lambda-CDM universe in closed form

The goal is to integrate the master equation from the big bang to today and obtain the present age t0t_0 as a clean closed form, not a black-box number. Ignoring radiation (negligible for the age integral) and curvature, the matter-plus-Λ\Lambda Friedmann equation integrates to an inverse hyperbolic sine — which is why the age of the universe has an exact analytic expression at all.

Step 1. Write cosmic time as t0=∫01da/a˙=1H0∫01da/[aΩma−3+ΩΛ]t_0=\int_0^1 \mathrm{d}a/\dot a=\tfrac{1}{H_0}\int_0^1 \mathrm{d}a/[a\sqrt{\Omega_m a^{-3}+\Omega_\Lambda}], using a˙=aH\dot a=aH and the two-component master equation. Why: the age is just the time to grow the scale factor from 00 to 11, and H(a)H(a) tells us how fast aa grows at each stage.

Step 2. Clear the negative power by multiplying numerator and denominator by a3/2a^{3/2}, giving t0=1H0∫01a da/Ωm+ΩΛa3t_0=\tfrac{1}{H_0}\int_0^1 \sqrt{a}\,\mathrm{d}a/\sqrt{\Omega_m+\Omega_\Lambda a^3}. Why: pulling the a−3a^{-3} out of the square root turns a messy integrand into one that a single substitution can flatten.

Step 3. Substitute u=a3/2u=a^{3/2}, so du=32a1/2da\mathrm{d}u=\tfrac32 a^{1/2}\mathrm{d}a and the integral becomes 23H0∫01du/Ωm+ΩΛu2\tfrac{2}{3H_0}\int_0^1 \mathrm{d}u/\sqrt{\Omega_m+\Omega_\Lambda u^2} — a standard arcsinh form. Why: the integrand is now 1/const+const⋅u21/\sqrt{\text{const}+\text{const}\cdot u^2}, whose antiderivative is an inverse hyperbolic sine, so the answer is exact.

Step 4. Evaluate to get t0=23H0ΩΛ arcsinhΩΛ/Ωmt_0=\tfrac{2}{3H_0\sqrt{\Omega_\Lambda}}\,\mathrm{arcsinh}\sqrt{\Omega_\Lambda/\Omega_m}; with H0=67.4 km s−1 Mpc−1H_0=67.4\ \mathrm{km\,s^{-1}\,Mpc^{-1}} (Hubble time 1/H0=14.51/H_0=14.5 Gyr) and Ωm=0.315\Omega_m=0.315 this gives t0≈13.8t_0\approx13.8 Gyr. Why: the closed form makes the age's dependence transparent — it scales as 1/H01/H_0 (so the H0H_0 tension is also a mild age tension) and rises with ΩΛ\Omega_\Lambda because dark energy stretches the expansion timeline. *Numerically: 230.685×14.5 Gyr×arcsinh(1.475)=13.80\tfrac{2}{3\sqrt{0.685}}\times14.5\ \mathrm{Gyr}\times\mathrm{arcsinh}(1.475)=13.80 Gyr, matching Planck's t0=13.797±0.023t_0=13.797\pm0.023 Gyr to better than 0.1%0.1\%.*

Limiting cases

The expansion factor E2(a)=Ωra−4+Ωma−3+Ωka−2+ΩΛE^2(a)=\Omega_r a^{-4}+\Omega_m a^{-3}+\Omega_k a^{-2}+\Omega_\Lambda collapses to single-component behaviour in each era, and the resulting a(t)a(t) laws are the backbone of the standard timeline.

CaseResultMeaning
Radiation era (a≪aeqa\ll a_{eq})a(t)∝t1/2a(t)\propto t^{1/2}the Ωra−4\Omega_r a^{-4} term dominates; sets nucleosynthesis and the early sound horizon
Matter era (aeq≪a≪aΛa_{eq}\ll a\ll a_\Lambda)a(t)∝t2/3a(t)\propto t^{2/3}the Ωma−3\Omega_m a^{-3} term dominates; recombination and structure growth, δ∝a\delta\propto a
Dark-energy era (a≫aΛa\gg a_\Lambda)a(t)∝eH0ΩΛ ta(t)\propto e^{H_0\sqrt{\Omega_\Lambda}\,t}the constant ΩΛ\Omega_\Lambda term dominates; de Sitter acceleration, the future attractor
Flat geometryΩk=1−Ωr−Ωm−ΩΛ≈0\Omega_k=1-\Omega_r-\Omega_m-\Omega_\Lambda\approx0CMB peak positions force ∣Ωk∣≲0.005\lvert\Omega_k\rvert\lesssim0.005; Lambda-CDM assumes exact flatness
Empty de Sitter (Ωm→0\Omega_m\to0)t0→∞t_0\to\infty, pure exponentialthe Λ\Lambda-only limit of the age formula; sets the asymptotic Hubble rate H∞=H0ΩΛH_\infty=H_0\sqrt{\Omega_\Lambda}

Worked examples

1 · Radiation-matter equality and the critical density (fully worked)

Using Ωm=0.315\Omega_m=0.315, Ωr=9.2×10−5\Omega_r=9.2\times10^{-5}, and H0=67.4 km s−1 Mpc−1H_0=67.4\ \mathrm{km\,s^{-1}\,Mpc^{-1}}, find the redshift of matter-radiation equality and the present critical density, and translate the latter into particles per cubic metre.

Step 1. Equality is where the two power-law terms balance, Ωra−4=Ωma−3\Omega_r a^{-4}=\Omega_m a^{-3}, so aeq=Ωr/Ωma_{eq}=\Omega_r/\Omega_m and 1+zeq=Ωm/Ωr1+z_{eq}=\Omega_m/\Omega_r. Plugging in, 1+zeq=0.315/(9.2×10−5)≈34201+z_{eq}=0.315/(9.2\times10^{-5})\approx3420, so zeq≈3400z_{eq}\approx3400 — comfortably before recombination at z≈1090z\approx1090, which is why the scales relevant for galaxies entered the horizon during matter domination and could grow.

Step 2. The critical density is ρc=3H02/8πG\rho_c=3H_0^2/8\pi G. Convert H0=67.4 km s−1 Mpc−1=2.18×10−18 s−1H_0=67.4\ \mathrm{km\,s^{-1}\,Mpc^{-1}}=2.18\times10^{-18}\ \mathrm{s^{-1}}, then ρc=3(2.18×10−18)2/[8π(6.674×10−11)]=8.5×10−27 kg m−3\rho_c=3(2.18\times10^{-18})^2/[8\pi(6.674\times10^{-11})]=\mathbf{8.5\times10^{-27}\ kg\,m^{-3}}.

Step 3. Divide by a proton mass 1.67×10−27 kg1.67\times10^{-27}\ \mathrm{kg}: that is ≈5.1\approx5.1 protons per cubic metre of critical density. Since baryons are only Ωb≈0.049\Omega_b\approx0.049 of it, ordinary matter amounts to about 0.250.25 protons per cubic metre — roughly one hydrogen atom in four cubic metres. The universe is staggeringly empty, yet that trickle of baryons plus the dark matter is enough to grow every galaxy. z_eq ≈ 3400, ρ_c ≈ 8.5×10⁻²⁷ kg m⁻³ ≈ 5 protons m⁻³.

2 · When did the cosmos start to accelerate (fully worked)

With Ωm=0.315\Omega_m=0.315 and ΩΛ=0.685\Omega_\Lambda=0.685, find the redshift at which dark energy overtook matter and the redshift at which the expansion began to accelerate. Then confirm the present age.

Step 1. Matter-Λ\Lambda equality: 1+zmΛ=(ΩΛ/Ωm)1/3=(0.685/0.315)1/3=(2.175)1/3=1.2961+z_{m\Lambda}=(\Omega_\Lambda/\Omega_m)^{1/3}=(0.685/0.315)^{1/3}=(2.175)^{1/3}=1.296, so zmΛ≈0.30z_{m\Lambda}\approx0.30. In cosmic time that is only about 3.53.5 billion years ago — dark energy is a recent tenant.

Step 2. Acceleration onset: 1+zacc=(2ΩΛ/Ωm)1/3=(4.35)1/3=1.6321+z_\mathrm{acc}=(2\Omega_\Lambda/\Omega_m)^{1/3}=(4.35)^{1/3}=1.632, so zacc≈0.63z_\mathrm{acc}\approx0.63. Note zacc>zmΛz_\mathrm{acc}>z_{m\Lambda}: the universe began speeding up before dark energy formally dominated the budget, because acceleration needs Λ\Lambda only to beat half the matter density.

Step 3. Age check: t0=23H0ΩΛ arcsinhΩΛ/Ωmt_0=\tfrac{2}{3H_0\sqrt{\Omega_\Lambda}}\,\mathrm{arcsinh}\sqrt{\Omega_\Lambda/\Omega_m}. With 1/H0=14.51/H_0=14.5 Gyr, ΩΛ/Ωm=1.475\sqrt{\Omega_\Lambda/\Omega_m}=1.475, and arcsinh(1.475)=1.181\mathrm{arcsinh}(1.475)=1.181: t0=23(0.828)(14.5)(1.181)=13.8 Gyrt_0=\tfrac{2}{3(0.828)}(14.5)(1.181)=\mathbf{13.8\ Gyr}, matching Planck's 13.79713.797 Gyr. The three numbers hang together: a universe that began accelerating at z≈0.63z\approx0.63, became Λ\Lambda-dominated at z≈0.30z\approx0.30, and is 13.813.8 Gyr old. z_mΛ ≈ 0.30, z_acc ≈ 0.63, t₀ ≈ 13.8 Gyr.

Physics of the portrait

The concordance model has no single author — it is the residue of a century of argument between theory and telescopes — but if one physicist embodies its theoretical spine it is Jim Peebles, who spent five decades turning cosmology from speculation into a precision science and shared the 2019 Nobel Prize for exactly that.

Origins · 1 minP. J. E. Peebles and how this was discovered
PJ
P. J. E. Peebles
b. 1935, Winnipeg, Canada · physical cosmology and the cold-dark-matter paradigm

Phillip James Edwin Peebles is, more than anyone, the architect of the theoretical framework that became Lambda-CDM. In 1965, working with Robert Dicke at Princeton, he predicted the temperature of a relic radiation background from the hot big bang — days before Penzias and Wilson's accidental detection of the cosmic microwave background was connected to it, a discovery Peebles's group immediately interpreted. Over the following decades he did the quantitative heavy lifting a working model demands: he computed primordial nucleosynthesis abundances, developed the statistical theory of galaxy clustering (the two-point correlation function is largely his language), and in the early 1980s helped establish the cold-dark-matter paradigm, showing that non-relativistic collisionless matter grows perturbations bottom-up in a way that matches the galaxy distribution. His textbooks trained generations. Notably, Peebles long resisted the cosmological constant, only to see the 1998 supernova results force Λ\Lambda back into the model he had built; he has since been the field's most eloquent skeptic, insisting that Lambda-CDM's empirical success does not excuse its unexplained ingredients. The 2019 Nobel Prize in Physics honoured his "theoretical discoveries in physical cosmology" — the quiet mathematics behind the concordance.

History
1917
— Einstein introduces the cosmological constant Λ into general relativity
1933
— Zwicky infers unseen dark matter in the Coma cluster
1965
— Penzias and Wilson detect the microwave background, interpreted by Dicke and Peebles
1998
— the supernova teams discover cosmic acceleration, reviving Λ
2013
— the Planck satellite fixes the six Lambda-CDM parameters to percent precision

Going further

Texts. Misner, Thorne & Wheeler, Gravitation (1973), Ch. 27 · Wald, General Relativity (1984), Ch. 5 · Carroll, Spacetime and Geometry (2004), Ch. 8 · Hartle, Gravity (2003), Ch. 18–19 · Schutz, A First Course in General Relativity (2009), Ch. 12 · Dodelson & Schmidt, Modern Cosmology (2nd ed., 2020) · Planck Collaboration, Planck 2018 Results VI: Cosmological Parameters, A&A 641, A6 (2020). All numerical values — density parameters, H0H_0, t0t_0, zeqz_{eq}, and zrecz_{rec} — are from Planck 2018 (A&A 641, A6). Key cross-checks: the master equation follows from substituting ρi∝a−3(1+wi)\rho_i\propto a^{-3(1+w_i)} into the first Friedmann equation; the quoted Ωi\Omega_i yield zmΛ=0.30z_{m\Lambda}=0.30, zacc=0.63z_\mathrm{acc}=0.63, and zeq≈3400z_{eq}\approx3400; the arcsinh age integral gives t0≈13.8t_0\approx13.8 Gyr for H0=67.4 km s−1 Mpc−1H_0=67.4\ \mathrm{km\,s^{-1}\,Mpc^{-1}}; and ρc≈8.5×10−27 kg m−3≈5\rho_c\approx8.5\times10^{-27}\ \mathrm{kg\,m^{-3}}\approx5 protons m⁻³.

Prerequisites. The FLRW Metric · The Friedmann Equations · Critical Density and the Density Parameters · The Cosmological Equation of State · The Cosmological Principle. Leads to. The Cosmological Constant and Dark Energy · The Dark-Energy-Dominated Era · The Matter-Dominated Era · Cosmological Redshift.

Ask the tutor: Ask me to walk you through why acceleration starts at a higher redshift than dark-energy domination, or to derive any of the era-transition redshifts from the Ωi\Omega_i you choose.

Is Lambda-CDM the same thing as the big bang theory?

Not quite — it is the big bang theory made quantitative and specific. "Big bang" names the broad idea that the universe expanded and cooled from a hot dense state. Lambda-CDM is the particular six-parameter model that fills in the details: a spatially flat FLRW universe with specified amounts of radiation, baryons, cold dark matter, and a cosmological constant, evolving by the Friedmann equations. You can accept the big bang framework and still argue about the ingredients; Lambda-CDM is the current best-fit choice of those ingredients. It is falsifiable in a way "the big bang" is not — which is exactly why the H0H_0 and S8S_8 tensions are taken so seriously.

What are the six parameters, and why exactly six?

The standard set is Ωbh2\Omega_b h^2 (baryons), Ωch2\Omega_c h^2 (cold dark matter), θs\theta_s or H0H_0 (the sound-horizon angle, i.e. the expansion rate), τ\tau (reionization optical depth), AsA_s (primordial fluctuation amplitude), and nsn_s (spectral tilt). "Exactly six" is an empirical claim, not a theorem: these six suffice to fit essentially all current data. The model assumes flatness (Ωk=0\Omega_k=0), a pure cosmological constant (w=−1w=-1), and standard neutrinos, each of which could be promoted to a free parameter. Testing "extensions to Lambda-CDM" means letting one of those fixed assumptions float and asking whether the data demand it.

If dark energy is a constant, where did all of it come from at late times?

None of it was "created" — the dark-energy density has been the same at every epoch. What changes is everything else. Matter and radiation densities fall as the universe expands (ρ∝a−3\rho\propto a^{-3} and a−4a^{-4}), while ρΛ\rho_\Lambda stays fixed. Early on, the enormous matter and radiation densities dwarfed the constant Λ\Lambda term, so it was dynamically irrelevant. As space expanded, the diluting components dropped below the fixed Λ\Lambda floor, and once they did, Λ\Lambda took over. So dark energy's "emergence" is really the fading of its competitors, not the arrival of new stuff — a subtle but important distinction.

Why is the cosmological constant such an embarrassment if the model works?

Because Lambda-CDM measures the value of Λ\Lambda superbly but offers no explanation for it. If Λ\Lambda is the energy of the quantum vacuum, dimensional estimates from quantum field theory overshoot the measured value by up to 1012010^{120} — often called the worst prediction in physics. Worse, there is a coincidence problem: the matter and dark-energy densities are within a factor of two of each other right now, even though their ratio spans many orders of magnitude over cosmic history, so we seem to live at a special moment. The model works as a description but fails as an explanation — which is why theorists keep probing modified gravity, quintessence, and anthropic arguments.

How can four completely different measurements agree — isn't that circular?

It would be circular if one measurement calibrated another, but they do not. The CMB acoustic peaks are physics at z≈1090z\approx1090; big-bang nucleosynthesis is physics in the first few minutes; Type Ia supernovae are late-time standard candles at z≲2z\lesssim2; baryon acoustic oscillations are a geometric ruler in the galaxy distribution. Each rests on independent physics and independent systematics, and each carves out a differently-oriented band in parameter space. That they intersect in one small region — for instance, Ωbh2\Omega_b h^2 from nucleosynthesis matching Ωbh2\Omega_b h^2 from CMB peaks measured 380,000380{,}000 years later — is a genuine, non-circular consistency test. The tensions that do exist (like H0H_0) are notable precisely because concordance is otherwise so tight.

Does the H0 tension mean Lambda-CDM is wrong?

Not necessarily, but it is the model's most serious challenge. The tension is that the local distance ladder (Cepheids calibrating Type Ia supernovae) gives H0≈73H_0\approx73, while extrapolating the CMB forward within Lambda-CDM gives H0≈67.4H_0\approx67.4, disagreeing at about 5σ5\sigma. Three possibilities remain live: an unidentified systematic, a statistical fluke (unlikely at 5σ5\sigma), or new physics — for example early dark energy that shrinks the sound horizon, or a non-standard neutrino sector. Because the CMB value is model-dependent (you must assume Lambda-CDM to extrapolate) while the local value is nearly model-independent, a real tension would point to a flaw in the model between recombination and today. It is the single most-watched number in cosmology.

Before this

Where this leads