advanced18 min studyLesson 851 in Physics

FLRW symmetry and Friedmann dynamics

Evolve a homogeneous cosmic background while keeping geometry, pressure, and thermal assumptions distinct.

On this page

What you’ll understand

  • State the assumptions behind an FLRW background.
  • Derive density scalings from stress-energy conservation.
  • Compute changing fluid fractions and cosmic acceleration.

An expanding universe is specified by a geometry and by the stress-energy that drives it. The Friedmann–Lemaître–Robertson–Walker description assumes spatial homogeneity and isotropy in a preferred family of comoving observers. These are background assumptions, not assertions that galaxies lack structure or that every particle species is in thermal equilibrium. Here and throughout this module, use natural units ℏ=c=kB=1\hbar=c=k_B=1 unless units are explicitly restored. Energy density means energy per physical volume, rather than mass density; this convention fixes the factors of light speed when converting formulas to SI.

Symmetry fixes a form, not an expansion history

Choose cosmic time tt as proper time along comoving worldlines and a dimensionless scale factor a(t)a(t). For comoving radial coordinate χ\chi with dimensions of length, the spatial-curvature parameter KK has dimensions of inverse length squared. A constant rescaling of the comoving coordinate changes the normalization of the scale factor. Observable redshifts depend on ratios of scale factors, so choosing unity at a reference epoch does not select a physical size for the universe.

Spatial symmetry makes the background stress tensor take perfect-fluid form with energy density ρ\rho and isotropic pressure PP. Cold matter has negligible background pressure; freely propagating ultrarelativistic particles have P=ρ/3P=\rho/3. Neither statement requires collisions at every instant. An isotropic collisionless radiation distribution can obey the radiation equation of state after decoupling. Local departures from homogeneity belong to perturbation theory, rather than an additional term inserted arbitrarily into the background metric.

ds2=−dt2+a2(t)[dχ21−Kχ2+χ2dΩ2],H=a˙ads^2=-dt^2+a^2(t)\left[\frac{d\chi^2}{1-K\chi^2}+\chi^2d\Omega^2\right],\qquad H=\frac{\dot a}{a}
The dot denotes cosmic-time differentiation, HH is the expansion rate, and dΩ2d\Omega^2 is the unit-sphere angular metric. The flat examples below set K=0K=0.↔ Scroll wide equations horizontally.

Pressure controls dilution and acceleration

For a component that exchanges no energy with other components, covariant stress-energy conservation reduces to a continuity equation. With constant equation-of-state parameter w=P/ρw=P/\rho, divide by the density and integrate against logarithmic scale factor. Matter then dilutes as inverse volume. Radiation dilutes faster because individual particle energies redshift as well. If components exchange energy, their separate continuity equations require equal-and-opposite source terms; only the total equation retains the source-free form.

Einstein gravity supplies the Friedmann constraint and acceleration equation. Positive density sets the expansion-rate magnitude, while density plus three times pressure controls acceleration. A cosmological constant can be represented as a fluid with negative pressure equal in magnitude to its density. Include it in the total density or as an explicit cosmological-constant term, never both. Expansion does not imply acceleration: a matter universe expands while its scale-factor growth slows.

H2=8πGρ3−Ka2,ρ˙+3H(ρ+P)=0,a¨a=−4πG3(ρ+3P),ρw∝a−3(1+w)H^2=\frac{8\pi G\rho}{3}-\frac{K}{a^2},\qquad \dot\rho+3H(\rho+P)=0,\qquad \frac{\ddot a}{a}=-\frac{4\pi G}{3}(\rho+3P),\qquad \rho_w\propto a^{-3(1+w)}
Here GG is Newton's constant. The density scaling assumes a separately conserved component with constant ww; the gravitational equations use the total density and pressure.↔ Scroll wide equations horizontally.

Fractions evolve even when their labels do not

Define each density fraction by dividing its density by the instantaneous critical density 3H2/(8πG)3H^2/(8\pi G). In a flat universe these fractions sum to unity. They are not permanent mixture weights: radiation loses importance relative to matter as expansion proceeds. A convenient calculation first evolves each density relative to one reference critical density, adds those contributions to obtain the new squared expansion rate, and only then renormalizes the fractions at the new epoch.

For a single constant-pressure-ratio fluid in a flat expanding universe, integration gives a power-law scale factor when w>−1w>-1. This solution chooses an origin of cosmic time and assumes that fluid dominates throughout the interval. Applying it across matter-radiation equality without matching the changing mixture produces a false time-temperature history. Similarly, a successful homogeneous solution does not prove the initial state was thermal, identify dark matter, or explain why the large-scale background was smooth.

Follow a changing matter-radiation mixture

Consider a stipulated flat universe containing only separately conserved cold matter and radiation. At a=1a=1, their fractions are 3/43/4 and 1/41/4, with expansion rate HrefH_{\rm ref}. At a=2a=2, find the expansion rate, both new fractions, and the deceleration parameter q=−a¨/(aH2)q=-\ddot a/(aH^2). These inputs describe a toy universe, not a fitted cosmic epoch.

  1. Evolve the reference-normalized density contributions separately. Matter contributes (3/4)2−3=3/32(3/4)2^{-3}=3/32 and radiation contributes (1/4)2−4=1/64(1/4)2^{-4}=1/64 to H2/Href2H^2/H_{\rm ref}^2. Their sum is 7/647/64. Taking the positive square root selects the expanding branch.

  2. Normalize each contribution by the new total, rather than preserving the reference fractions. The matter fraction is (6/64)/(7/64)=6/7(6/64)/(7/64)=6/7, and the radiation fraction is 1/71/7. The faster dilution of radiation explains the shift.

  3. Use the acceleration equation to obtain q=Ωm/2+Ωrq=\Omega_m/2+\Omega_r. Substitution gives 3/7+1/7=4/73/7+1/7=4/7. Its positive sign means decelerating expansion even though the expansion rate remains positive.

The result

The new rate is H/Href=7/8≃0.3307H/H_{\rm ref}=\sqrt7/8\simeq0.3307, the fractions are 6/76/7 and 1/71/7, and q=4/7q=4/7. Radiation contributes less to the density yet more acceleration per unit density than cold matter.

Explore the model

Interactive lab

Interactive lab · matrix

Separate a scale parameter from a volume claim

Reset and set both controls to diagonal scale a=1a=1 and shear b=0b=0; the selected final determinant is 11. Reset and set both controls to a=2a=2, b=0b=0; it is 22. Geometry begins at identity after Reset and interpolates during playback, while readouts describe the selected final matrix. This two-dimensional one-axis area scaling is a bookkeeping analogy, not isotropic three-dimensional FLRW expansion.

Separate a scale parameter from a volume claimThe selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA geometric scale has meaning only after the transformation and dimension are specified; this lab's determinant is not the cosmic volume factor a3a^3.

Check your understanding

For a separately conserved ultrarelativistic component, why does density scale as a−4a^{-4}?

Try it yourself

A different flat expanding toy universe contains one separately conserved fluid with w=−1/3w=-1/3 and positive density. If its scale factor doubles, by what factor does its density change? Determine whether its expansion accelerates and find the power-law exponent after choosing the zero-size time as the origin.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Use continuity first, then solve the flat Friedmann equation with the resulting density scaling.

Compare with the solution

The density scales as a−2a^{-2}, so doubling the scale factor reduces it to one quarter. Since ρ+3P=0\rho+3P=0, the acceleration vanishes. The flat Friedmann equation gives a˙2\dot a^2 constant and positive, hence a∝ta\propto t on the expanding branch. This fluid has negative pressure but sits exactly at the boundary between acceleration and deceleration; negative pressure alone is insufficient for acceleration.

Take it with you

  • FLRW symmetry does not establish thermal equilibrium.
  • Continuity fixes dilution only after energy exchange and the equation of state are specified.
  • Recompute density fractions using the critical density at the epoch of interest.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.