The ΛCDM Paradigm: Foundations and Tensions

The ΛCDM paradigm stands as the cornerstone of modern cosmology, a model of remarkable simplicity and predictive power that encapsulates the thermal history of the universe in just six free parameters. Yet, beneath this veneer of success lie growing tensions—most prominently in the Hubble constant ($H_0$) and the structure growth parameter ($S_8$)—that challenge the paradigm’s completeness. This chapter constructs the ΛCDM model from first principles, derives its observable predictions, and subjects it to rigorous statistical scrutiny against the latest cosmological datasets. We will see how the model’s elegance coexists with profound puzzles, setting the stage for the explorations of dark matter and dark energy that follow.

The Friedmann-Lemaître-Robertson-Walker Framework

The cosmological principle—the assumption of homogeneity and isotropy on large scales—dictates that the spacetime metric must be of the Friedmann-Lemaître-Robertson-Walker (FLRW) form. In comoving coordinates $(t, \mathbf{x})$, the line element is given by

$$ds^2 = -c^2 dt^2 + a^2(t) \left[ \frac{dr^2}{1 - k r^2} + r^2 (d\theta^2 + \sin^2\theta \, d\phi^2) \right],$$

where $a(t)$ is the scale factor (normalized to unity at the present epoch, $a(t_0)=1$), and $k \in \{-1, 0, +1\}$ encodes the spatial curvature. The Einstein field equations,

$$G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu},$$

when applied to the FLRW metric and a perfect fluid energy-momentum tensor $T^{\mu}_{\ \nu} = \text{diag}(-\rho c^2, p, p, p)$, yield the Friedmann equations:

$$H^2(t) \equiv \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho(t) - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}, \tag{1}$$

$$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3} \left( \rho + \frac{3p}{c^2} \right) + \frac{\Lambda c^2}{3}. \tag{2}$$

Here $H(t)$ is the Hubble parameter, $\rho$ the total energy density, and $p$ the pressure. The cosmological constant $\Lambda$ is interpreted as a fluid with equation of state $p_\Lambda = -\rho_\Lambda c^2$, where $\rho_\Lambda = \Lambda c^2 / (8\pi G)$. Equation (1) is a direct consequence of the $00$-component of the field equations, while (2) follows from the spatial components or, equivalently, from energy-momentum conservation $\nabla_\mu T^{\mu\nu}=0$, which gives the fluid equation:

$$\dot{\rho} + 3H \left( \rho + \frac{p}{c^2} \right) = 0. \tag{3}$$

For a component with constant equation of state $w \equiv p/(\rho c^2)$, integration of (3) yields $\rho(a) = \rho_0 a^{-3(1+w)}$. This scaling law underlies the cosmic inventory: non-relativistic matter ($w=0$, $\rho_m \propto a^{-3}$), radiation ($w=1/3$, $\rho_r \propto a^{-4}$), and dark energy with $w \approx -1$ ($\rho_\Lambda = \text{const}$). The curvature term scales as $a^{-2}$, making it dynamically insignificant at late times if not exactly zero.

The critical density $\rho_c(t) = 3H^2(t)/(8\pi G)$ defines the density parameter $\Omega_i(t) = \rho_i(t)/\rho_c(t)$. At present, the total density parameter $\Omega_{\text{tot}} = \sum_i \Omega_i$ determines curvature via $\Omega_k = 1 - \Omega_{\text{tot}} = -kc^2/(a_0^2 H_0^2)$. Observations of the cosmic microwave background (CMB) anisotropies, particularly the angular scale of the first acoustic peak, constrain $\Omega_k$ to be extremely close to zero, implying a spatially flat universe to within $0.2\%$ (Planck Collaboration 2018). Hence, we set $k=0$ henceforth, simplifying (1) to

$$H^2(a) = H_0^2 \left[ \Omega_r a^{-4} + \Omega_m a^{-3} + \Omega_\Lambda \right], \tag{4}$$

with $\Omega_m = \Omega_b + \Omega_c$ comprising baryons ($\Omega_b$) and cold dark matter ($\Omega_c$), and $\Omega_r$ including photons and three massless neutrino species (effectively contributing $\Delta N_{\text{eff}} \approx 0.027$ to the effective number of relativistic degrees of freedom). The sum of density parameters equals unity: $\Omega_r + \Omega_m + \Omega_\Lambda = 1$.

Thermodynamic Composition of the Universe

The thermal history of the universe is governed by the interplay between expansion and particle interactions. The photon-baryon plasma, tightly coupled via Thomson scattering prior to recombination, supports acoustic oscillations whose imprint is seen in the CMB power spectrum. To model this, we must characterize each component’s equation of state and evolution.

Radiation Domination and Neutrino Decoupling

At early times ($a \ll 10^{-4}$), radiation dominates. The photon temperature $T_\gamma \propto a^{-1}$; the CMB today has $T_{\gamma,0} = 2.72548 \pm 0.00057$ K (Fixsen 2009). Neutrinos decouple at $T \sim 2-3$ MeV, after which their momenta redshift as $a^{-1}$. Electron-positron annihilation heats the photons relative to neutrinos, yielding $T_{\nu,0} = (4/11)^{1/3} T_{\gamma,0}$. The effective number of relativistic species is $N_{\text{eff}} = 3.045$ in the standard model (de Salas et al. 2016), contributing to $\Omega_r$ as

$$\Omega_r h^2 = \Omega_\gamma h^2 \left(1 + 0.2271 N_{\text{eff}}\right),$$

where $\Omega_\gamma h^2 = 2.469 \times 10^{-5} (T_{\gamma,0}/2.7255\ \text{K})^4$ for $h = H_0/(100\ \text{km/s/Mpc})$. This radiation density affects the expansion rate during the radiation era, thereby shifting the sound horizon at recombination—a key standard ruler for BAO and CMB analyses.

Matter Domination and Recombination

As the universe expands, matter ($\Omega_m$) eventually dominates at $z_{\text{eq}} \approx 3400$, where $\rho_m = \rho_r$. The matter density includes baryons $\Omega_b h^2$ and cold dark matter $\Omega_c h^2$, both scaling as $a^{-3}$. Baryons remain ionized until recombination, when the rate of electron-proton recombination exceeds the expansion rate. The residual ionization fraction $x_e$ determines the Thomson scattering optical depth

$$\tau = \int_{t_{\text{rec}}}^{t_0} n_e(t) \sigma_T c \, dt,$$

where $n_e$ is the free electron density and $\sigma_T$ the Thomson cross-section. The peak of the visibility function at $z_{\text{rec}} \approx 1089$ (Planck 2018) defines the surface of last scattering. Precise modeling of recombination, including radiative transfer effects and two-photon processes, is essential for predicting the CMB damping tail (Seager et al. 1999). Codes like *RECFAST* and its successors *HyRec* and *CosmoRec* compute $x_e(z)$ to better than $0.1\%$ accuracy, enabling sub-percent cosmological parameter inference.

Dark Energy and Late-Time Acceleration

The discovery of cosmic acceleration from Type Ia supernovae (Riess et al. 1998; Perlmutter et al. 1999) introduced the need for a component with $w < -1/3$. In ΛCDM, this is the cosmological constant, with $w = -1$ exactly. Its energy density remains constant while matter dilutes, leading to a recent transition to acceleration at $z_{\text{acc}} \approx 0.7$. Alternative models parametrize dark energy with a time-varying equation of state $w(a) = w_0 + w_a (1-a)$, but current data are consistent with $w = -1$ to within $5\%$ (Planck 2018). The flatness condition then forces $\Omega_\Lambda = 1 - \Omega_m - \Omega_r \approx 0.69$.

Linear Perturbation Theory and CMB Anisotropies

The observed CMB temperature fluctuations $\Delta T/T \sim 10^{-5}$ are the linear perturbations $\Theta(\mathbf{k}, \hat{n})$ on the FLRW background. In the tight-coupling regime before recombination, photons and baryons behave as a single fluid with sound speed $c_s = c/\sqrt{3(1+R)}$, where $R \equiv 3\rho_b/(4\rho_\gamma)$ is the baryon-to-photon momentum density ratio. The photon temperature perturbation in Fourier space obeys, in the absence of gravity,

$$\ddot{\Theta} + \frac{\dot{a}}{a} \dot{\Theta} + k^2 c_s^2 \Theta = 0, \tag{5}$$

which describes acoustic oscillations. Gravity from dark matter perturbations introduces a driving term, enhancing compressional peaks. The solution for the photon-baryon fluid yields a characteristic scale: the sound horizon at recombination,

$$r_s(z_*) = \int_{z_*}^{\infty} \frac{c_s(z)}{H(z)} \, dz, \tag{6}$$

wh