The Cosmological Fluid and Equation of State

The universe is a preposterously complicated object — a hundred billion galaxies of stars, gas, dust, black holes, streaming neutrinos, and a dark sector nobody has held in their hand. Yet cosmology extracts a working model of the whole thing from a single act of courage: on scales larger than about 100 Mpc, throw away every detail and treat each ingredient as a perfect fluid, a smooth continuum with nothing but an energy density and a pressure. The Cosmological Principle licenses the move — homogeneity and isotropy forbid shear, viscosity, and heat flow, so the stress-energy tensor has no choice but the perfect-fluid form, its four-velocity locked to the comoving observers of the FLRW metric. A perfect fluid is two numbers, ρc2\rho c^2 and pp, and to make them dynamical you need one relation between them: the equation of state. For cosmology it is astonishingly simple — linear, p=w ρc2p = w\,\rho c^2, a single dimensionless constant ww per species. That one number does two enormous jobs. Fed through the conservation law ∇μTμν=0\nabla_\mu T^{\mu\nu}=0, it dictates how a fluid's density dilutes as space stretches, ρ∝a−3(1+w)\rho \propto a^{-3(1+w)}. Fed through the acceleration equation, it decides whether the fluid's gravity halts the expansion or drives it faster. Matter (w=0w=0), radiation (w=1/3w=1/3), and vacuum (w=−1w=-1) fall out as three values, and because they dilute at different rates the cosmos automatically sorts itself into a radiation era, a matter era, and the dark-energy era we live in now. This article is where thermodynamics, conservation laws, and the dynamics of spacetime collapse into one number.

Key takeaways
  • The contents of the universe are modeled on the largest scales as non-interacting perfect fluids, each fixed by an energy density ρc2\rho c^2 and an isotropic pressure pp tied together by an equation of state p=w ρc2p = w\,\rho c^2.
  • Local energy conservation ∇μTμν=0\nabla_\mu T^{\mu\nu}=0 in an FLRW spacetime gives the cosmological continuity equation ρ˙+3H(ρ+p/c2)=0\dot\rho + 3H(\rho + p/c^2)=0, which for constant ww integrates to the scaling law ρ∝a−3(1+w)\rho \propto a^{-3(1+w)}.
  • The three canonical values are w=0w=0 (pressureless matter, ρ∝a−3\rho\propto a^{-3}), w=1/3w=1/3 (radiation, ρ∝a−4\rho\propto a^{-4}), and w=−1w=-1 (a cosmological constant, ρ=\rho= const).
  • Pressure gravitates: the source of the acceleration equation is ρc2+3p\rho c^2 + 3p, so a¨>0\ddot a>0 requires w<−1/3w<-1/3 — only sufficiently negative pressure produces cosmic repulsion.
  • Because each species dilutes at its own rate, the energy budget is dominated successively by radiation, then matter, then dark energy — the origin of the cosmic eras, with matter–radiation equality near z≈3400z\approx3400.
  • Key idea: one dimensionless number ww in p=wρc2p = w\rho c^2 controls everything — it fixes how a fluid dilutes via ρ∝a−3(1+w)\rho \propto a^{-3(1+w)}, and whether that fluid decelerates (w>−1/3w>-1/3) or accelerates (w<−1/3w<-1/3) the cosmic expansion.
The core insight

One dimensionless number ww in p=wρc2p = w\rho c^2 controls everything: it fixes how a fluid's density dilutes with expansion via ρ∝a−3(1+w)\rho \propto a^{-3(1+w)}, and whether that fluid decelerates (w>−1/3w>-1/3) or accelerates (w<−1/3w<-1/3) the cosmic expansion.

Where this lives today

The equation of state is not a settled number but a central measurement of twenty-first-century cosmology. The dark-energy surveys — the Dark Energy Survey, ESA's Euclid (launched 2023), and the Vera Rubin Observatory's LSST — exist to answer one question: is ww exactly −1-1, a true cosmological constant, or slightly different, and does it drift with time? The tools are geometric: Type Ia supernovae trace H(z)H(z), baryon acoustic oscillations imprint a standard ruler in the galaxy distribution, and weak lensing weighs the growth of structure — each folding ww into an observable through ρ∝a−3(1+w)\rho\propto a^{-3(1+w)}. Current constraints hover around w=−1.0±0.03w=-1.0\pm0.03, consistent with a bare cosmological constant, though 2024 DESI results hinted at a mild preference for an evolving ww — which, if it survives, would rule out the constant and demand a dynamical field. The abstract barotropic parameter of this article is, quite literally, the quantity billion-dollar observatories are built to pin down.

The physical picture

Hold two pictures side by side. In the first, a comoving box of space is stretched uniformly in all directions and you watch whatever is inside. Cold matter is a fixed count of slow particles: the box grows, the number stays the same, so density drops purely as the volume, a−3a^{-3}. Radiation is a gas of photons: the same volume dilution applies, but every photon is also redshifted, its energy bled away as a−1a^{-1} — one extra power, giving a−4a^{-4}. Vacuum energy is the strange one: its density per unit volume never changes, so stretching empty space just makes more of the same energy, a0a^{0}. Three ingredients, three fates, one expansion. The second picture is the master relation itself, on a single axis. The number w=p/ρc2w=p/\rho c^2 is the pressure a fluid carries per unit of its energy density, and where a fluid sits on that axis fixes both its dilution law and its gravitational sign. This is the counter-intuitive heart of relativistic cosmology: pressure gravitates. Positive pressure — radiation pushing outward — adds to the attractive source and slows the expansion. Only tension, pressure below the threshold w=−1/3w=-1/3, gravitates repulsively and drives acceleration.

Figure
The master relation places every cosmic fluid on one axis. The equation-of-state number w fixes both the dilution ρ ∝ a⁻³⁽¹⁺ʷ⁾ and, through ä > 0 ⟺ w < −1/3, whether the fluid halts the expansion or drives it faster. Matter sits at w = 0, radiation at w = 1/3, a cosmological constant at w = −1, and the curvature-like threshold at w = −1/3
One number w fixes both the dilution law and the sign of gravity accelerates · repulsive · w below −1/3 decelerates · attractive · w above −1/3 w = −1 vacuum / Λ ρ = const w = −1/3 curvature-like ρ ∝ a⁻² w = 0 matter · ρ ∝ a⁻³ w = 1/3 radiation · ρ ∝ a⁻⁴ w = 1 stiff · ρ ∝ a⁻⁶ dilution: ρ ∝ a⁻³⁽¹⁺ʷ⁾ the universe accelerates when w falls below −1/3

Because the three fluids dilute at different rates, their relative importance marches monotonically with the scale factor. Run the clock backward and radiation, falling as a−4a^{-4}, wins; run it forward and vacuum energy, refusing to dilute, wins in the end. The crossover points are the boundaries of the cosmic eras, set entirely by today's densities and the exponents −3(1+w)-3(1+w). The second figure is what a cosmologist actually plots: the energy density of each component against the scale factor, on log-log axes where the power laws are straight lines and the era boundaries are where those lines cross.

Figure
Energy density of each component versus scale factor on log-log axes, in units of today's critical density. The straight lines are the power laws ρ ∝ a⁻⁴ (radiation), a⁻³ (matter), a⁰ (Λ). They cross at matter–radiation equality near z ≈ 3400 and matter–Λ equality near z ≈ 0.3, sorting cosmic history into three eras
The three fluids dilute at different rates — cosmic history sorts into eras radiation era matter era Λ era 10⁻⁵ 10⁻⁴ 10⁻³ 10⁻² 10⁻¹ 1 scale factor a (log scale; a = 1 today) energy density (log, units of today's ρ_crit) ρ ∝ a⁻⁴ ρ ∝ a⁻³ ρ_Λ = const matter–radiation equality z ≈ 3400 matter–Λ equality z ≈ 0.3 today

The core equations

Everything begins with the shape of the stress-energy tensor. Homogeneity and isotropy leave one option — no momentum flux, no shear, no anisotropic stress — so the source of gravity is built from only the fluid four-velocity and the metric, the most general stress-energy compatible with the Cosmological Principle.

Tμν=(ρ+pc2)uμuν+p gμν(1)T^{\mu\nu} = \left(\rho + \frac{p}{c^2}\right) u^\mu u^\nu + p\, g^{\mu\nu} \qquad (1)

where: ρ\rho — rest-frame mass-energy density (so ρc2\rho c^2 is the energy density); pp — isotropic pressure in the fluid rest frame; uμu^\mu — fluid four-velocity, normalized uμuμ=−c2u^\mu u_\mu = -c^2; gμνg^{\mu\nu} — inverse FLRW metric. Reading: in the comoving rest frame uμ=(c,0,0,0)u^\mu=(c,0,0,0) this collapses to T00=−ρc2T^0{}_0=-\rho c^2 and Tij=p δjiT^i{}_j = p\,\delta^i_j — energy density on the time-time slot, isotropic pressure on the spatial diagonal, nothing off-diagonal. To close the system you need one thermodynamic relation between pp and ρ\rho. For the barotropic fluids of cosmology it is linear, and the single slope ww is its entire content.

p=w ρc2(2)p = w\,\rho c^2 \qquad (2)

where: ww — dimensionless equation-of-state parameter, w=p/(ρc2)w = p/(\rho c^2), constant (or slowly varying) per species. Reading: ww is a pure ratio of pressure to energy density; standard cosmology lives at w∈{0, 1/3, −1}w \in \{0,\,1/3,\,-1\}. Demanding local energy-momentum conservation — automatic once Einstein's equation holds — and taking the time component in FLRW gives the cosmological continuity equation, the fluid's equation of motion in an expanding space.

ρ˙+3a˙a(ρ+pc2)=0(3)\dot\rho + 3\frac{\dot a}{a}\left(\rho + \frac{p}{c^2}\right) = 0 \qquad (3)

where: a(t)a(t) — FLRW scale factor; a˙/a≡H\dot a/a \equiv H — Hubble parameter; ρ˙\dot\rho — time derivative along comoving worldlines. Reading: density changes for two reasons — dilution (3Hρ3H\rho) and the work done by pressure (3Hp/c23Hp/c^2). It is nothing but the first law of thermodynamics d(ρc2V)=−p dV\mathrm{d}(\rho c^2 V) = -p\,\mathrm{d}V for an adiabatically expanding comoving volume V∝a3V\propto a^3. Inserting the equation of state and integrating (Derivation A) turns the differential law into the algebraic scaling law every cosmology calculation begins from.

ρ(a)=ρ0(aa0)−3(1+w)(4)\rho(a) = \rho_0\left(\frac{a}{a_0}\right)^{-3(1+w)} \qquad (4)

where: ρ0\rho_0 — present density of the component; a0a_0 — present scale factor (usually set to 1). Reading: one exponent, −3(1+w)-3(1+w), does all the work — w=0⇒a−3w=0\Rightarrow a^{-3}, w=1/3⇒a−4w=1/3\Rightarrow a^{-4}, w=−1⇒a0w=-1\Rightarrow a^0 — exactly the exponent of each term of the Friedmann equation. The dynamical payoff is the acceleration equation, the time-time projection of Einstein's equation. Its source is not the energy density alone but ρc2+3p\rho c^2 + 3p — the reason pressure gravitates.

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

where: a¨\ddot a — second time derivative of the scale factor; GG — Newton's constant. Reading: positive pressure (w>0w>0) deepens the deceleration; a¨\ddot a flips sign only when the pressure is negative enough to make 1+3w<01+3w<0. That sign flip is the most important inequality in cosmology; it defines dark energy.

a¨>0  ⟺  1+3w<0  ⟺  w<−13(6)\ddot a > 0 \iff 1+3w < 0 \iff w < -\frac{1}{3} \qquad (6)

where: a¨>0\ddot a>0 — the expansion is speeding up. Reading: a fluid accelerates the universe only if w<−1/3w<-1/3; a cosmological constant at w=−1w=-1 clears the bar comfortably, curvature at w=−1/3w=-1/3 sits exactly on it, and all ordinary matter and radiation fail it.

Derivation AIntegrating the fluid equation to the density scaling law

The goal of this derivation is to turn the cosmological continuity equation — a first-order differential equation relating a density to the expansion rate — into the closed-form power law ρ∝a−3(1+w)\rho \propto a^{-3(1+w)} that underlies every term of the Friedmann equation. The only physical input beyond energy conservation is that ww is constant for the species in question; everything else is a clean separation of variables and one integration, with the constant fixed by matching to today's density.

Step 1. Start from the continuity equation ρ˙+3a˙a(ρ+p/c2)=0\dot\rho + 3\frac{\dot a}{a}\left(\rho + p/c^2\right)=0, the ν=0\nu=0 component of ∇μTμν=0\nabla_\mu T^{\mu\nu}=0 in FLRW. Why this is the right starting point: it is the exact statement of local energy conservation, the 3a˙/a3\dot a/a arising from the spatial trace of the connection Γi0j=(a˙/a)δji.\Gamma^i{}_{0j}=(\dot a/a)\delta^i_j.

Step 2. Insert the equation of state p/c2=wρp/c^2 = w\rho so the bracket becomes ρ(1+w)\rho(1+w), giving ρ˙=−3a˙a ρ(1+w)\dot\rho = -3\frac{\dot a}{a}\,\rho(1+w). Why this simplifies things: it replaces two unknown functions ρ\rho and pp by a single one, because the equation of state ties them rigidly together.

Step 3. Divide by ρ\rho to separate variables, ρ˙/ρ=−3(1+w) a˙/a\dot\rho/\rho = -3(1+w)\,\dot a/a, and recognize both sides as time derivatives of logarithms. Why the constancy of ww matters here: only if ww does not vary can the right side be integrated in closed form — a time-varying ww leaves an integral ∫w dln⁡a\int w\,\mathrm d\ln a that does not collapse to a power law.

Step 4. Integrate to ln⁡ρ=−3(1+w)ln⁡a+const\ln\rho = -3(1+w)\ln a + \text{const} and exponentiate, then match the constant to the present epoch by demanding ρ=ρ0\rho=\rho_0 at a=a0a=a_0. Why the boundary condition is chosen this way: densities are measured today, so anchoring the power law at a0a_0 makes ρ0\rho_0 the observable normalization that surveys report.

Step 5. The result is ρ(a)=ρ0 (a/a0)−3(1+w)\rho(a)=\rho_0\,(a/a_0)^{-3(1+w)}. Setting w=0,13,−1w=0,\tfrac13,-1 returns a−3,a−4,a0a^{-3},a^{-4},a^0 respectively. Why this is the whole story: substituting these three exponents into the Friedmann equation reproduces the entire radiation-to-matter-to-vacuum chronology with no further physics. *Verified: reduces to MTW §27 and Carroll Eq. (8.36); w=0⇒a−3w=0\Rightarrow a^{-3}, w=1/3⇒a−4w=1/3\Rightarrow a^{-4}, w=−1⇒w=-1\Rightarrow constant.*

Go deeper · Derivation B · 1 minThe acceleration equation and the threshold for cosmic repulsion
Derivation BThe acceleration equation and the threshold for cosmic repulsion

The goal here is to derive the acceleration equation a¨/a=−4πG3ρ(1+3w)\ddot a/a = -\tfrac{4\pi G}{3}\rho(1+3w) directly from the first Friedmann equation and the continuity equation, and from it read off the exact equation of state at which the universe stops decelerating and begins to accelerate. This is the calculation that converts the abstract slogan "pressure gravitates" into the sharp, testable threshold w<−1/3w<-1/3 separating ordinary gravitating matter from dark energy, using nothing but differentiation of a relation we already have.

Step 1. Take the first Friedmann equation (a˙/a)2=8πG3ρ−kc2/a2\left(\dot a/a\right)^2 = \tfrac{8\pi G}{3}\rho - kc^2/a^2 and differentiate once with respect to cosmic time. Why start from the Friedmann constraint: it already encodes Einstein's equation for the FLRW metric, so its time derivative must contain the dynamical content of the acceleration.

Step 2. The differentiation yields 2a˙a(a¨a−a˙2a2)=8πG3ρ˙+2kc2a˙/a32\frac{\dot a}{a}\left(\frac{\ddot a}{a}-\frac{\dot a^2}{a^2}\right) = \tfrac{8\pi G}{3}\dot\rho + 2kc^2\dot a/a^3. Why keep the curvature term: it looks like it should matter, but it cancels in the next step against the a˙2/a2\dot a^2/a^2 piece from the Friedmann relation.

Step 3. Eliminate ρ˙\dot\rho using the continuity equation ρ˙=−3a˙aρ(1+w)\dot\rho = -3\frac{\dot a}{a}\rho(1+w), and substitute the Friedmann expression back for a˙2/a2\dot a^2/a^2. Why this is the crux: it injects the equation of state — the factor (1+w)(1+w) — into what was a purely kinematic identity, so the pressure enters the dynamics precisely here.

Step 4. After cancellation the curvature term and the a˙2/a2\dot a^2/a^2 term disappear, leaving a¨/a=−4πG3ρ(1+3w)\ddot a/a = -\tfrac{4\pi G}{3}\rho(1+3w). Why the curvature drops out: acceleration is a local dynamical quantity sourced by the fluid, independent of the global spatial geometry, which only sets the expansion rate, not its rate of change.

Step 5. The sign of a¨\ddot a is the sign of −(1+3w)-(1+3w), so a¨>0\ddot a>0 demands 1+3w<01+3w<0, i.e. w<−1/3w<-1/3; at w=−1w=-1 this gives the de Sitter result a¨/a=+8πG3ρΛ>0\ddot a/a = +\tfrac{8\pi G}{3}\rho_\Lambda>0. Why −1/3-1/3 and not 00: the "33" is the three spatial dimensions over which pressure acts as a source, so it takes pressure equal to −13-\tfrac13 of the energy density just to make the fluid gravitationally weightless. *Verified: matches Carroll Eq. (8.23) and Hartle Ch. 18; w=−1w=-1 gives exponential de Sitter expansion.*

Go deeper · Derivation C · 1 minWhy radiation carries one extra power of the scale factor
Derivation CWhy radiation carries one extra power of the scale factor

This derivation explains why radiation dilutes as a−4a^{-4} while matter dilutes as a−3a^{-3} — a single extra power of the scale factor whose origin is the cosmological redshift of each individual quantum. We reach the same w=1/3w=1/3 two independent ways: by counting photons in a comoving box and tracking their energies, and by taking the trace of the electromagnetic stress-energy tensor, which must vanish for a massless field. The agreement ties thermodynamics to field theory.

Step 1. In a comoving region the photon number NN is conserved and the proper volume grows as V∝a3V\propto a^3, so the number density falls as n=N/V∝a−3n = N/V \propto a^{-3} — the same volume dilution matter feels. Why number is conserved: away from interactions, photons are neither created nor destroyed as space expands, so only the geometry dilutes them.

Step 2. Each photon's wavelength stretches with the expansion, λ∝a\lambda\propto a, so its energy E=hc/λ∝a−1E=hc/\lambda\propto a^{-1} — this is the cosmological redshift acting on every quantum. Why matter does not share this factor: a massive particle's energy is dominated by its rest mass mc2mc^2, which does not redshift, so cold matter loses only the volume factor.

Step 3. Multiply density by mean energy: ρrc2=n⟨E⟩∝a−3⋅a−1=a−4\rho_r c^2 = n\langle E\rangle \propto a^{-3}\cdot a^{-1} = a^{-4}, and matching to −3(1+w)=−4-3(1+w)=-4 gives w=1/3w=1/3. Why this is the physical content of w=1/3w=1/3: it is literally the statement that radiation suffers one redshift factor on top of the volume dilution shared with matter.

Step 4. Independently, the electromagnetic stress-energy is trace-free, Tμμ=0T^\mu{}_\mu=0; for a perfect fluid the trace is −ρc2+3p-\rho c^2 + 3p, so 3p=ρc23p=\rho c^2 and p=13ρc2p=\tfrac13\rho c^2, i.e. w=1/3w=1/3 exactly. Why the two routes must agree: the trace-free condition is masslessness expressed in field theory, and the redshift argument is masslessness expressed in kinematics — one physics, two languages. *Verified: photon-gas pressure p=u/3p=u/3 with u=ρc2u=\rho c^2 gives w=1/3w=1/3, consistent with Stefan–Boltzmann u∝T4u\propto T^4 and T∝a−1T\propto a^{-1}.*

Limiting cases

The equation-of-state catalogue is short, and every entry is a value of ww read straight into the exponent −3(1+w)-3(1+w).

CaseResultMeaning
Pressureless matter (dust)w=0,ρm∝a−3w=0,\quad \rho_m \propto a^{-3}Cold dark matter and baryons; non-relativistic, negligible pressure. Density falls purely from volume dilution.
Radiation / relativistic speciesw=13,ρr∝a−4w=\tfrac{1}{3},\quad \rho_r \propto a^{-4}Photons and relativistic neutrinos. The extra a−1a^{-1} over matter is the cosmological redshift of each quantum.
Cosmological constant / vacuumw=−1,ρΛ=constw=-1,\quad \rho_\Lambda = \text{const}Dark energy in ΛCDM. Density independent of aa; drives exponential de Sitter expansion.
Spatial curvature (effective)w=−13,ρk∝a−2w=-\tfrac{1}{3},\quad \rho_k \propto a^{-2}The curvature term in the Friedmann equation behaves like a fluid with w=−1/3w=-1/3 — the exact boundary between deceleration and acceleration.
Stiff fluidw=1,ρ∝a−6w=1,\quad \rho \propto a^{-6}Maximal causal ww (sound speed equals cc); the kinetic energy of a free scalar field. Dilutes fastest of all.
Phantom energyw<−1,ρ grows with aw<-1,\quad \rho \text{ grows with } aDensity increases as space expands; violates the dominant energy condition and ends in a finite-time "Big Rip".

Worked examples

1 · Matter–radiation equality, in redshift and temperature (fully worked)

Take today's density parameters Ωm,0=0.315\Omega_{m,0}=0.315 for matter and Ωr,0=9.2×10−5\Omega_{r,0}=9.2\times10^{-5} for radiation (photons plus three relativistic neutrino species), and the CMB temperature T0=2.725T_0=2.725 K. At what redshift were matter and radiation equally important, and how hot was the universe then?

Step 1. Write each density with its scaling law relative to today: ρm(a)=ρm,0 a−3\rho_m(a)=\rho_{m,0}\,a^{-3} and ρr(a)=ρr,0 a−4\rho_r(a)=\rho_{r,0}\,a^{-4}, with a0=1a_0=1. Because the two exponents differ by one, their ratio is a single power of aa: ρr/ρm=(Ωr,0/Ωm,0) a−1\rho_r/\rho_m = (\Omega_{r,0}/\Omega_{m,0})\,a^{-1}.

Step 2. Equality is where the ratio equals one, so aeq=Ωr,0/Ωm,0=9.2×10−5/0.315=2.92×10−4a_{\rm eq}=\Omega_{r,0}/\Omega_{m,0}=9.2\times10^{-5}/0.315 = 2.92\times10^{-4}. In redshift, 1+zeq=1/aeq=Ωm,0/Ωr,0=3.42×1031+z_{\rm eq}=1/a_{\rm eq}=\Omega_{m,0}/\Omega_{r,0}=3.42\times10^{3}, so zeq≈3400z_{\rm eq}\approx 3400.

Step 3. The photon temperature scales as T∝a−1=(1+z)T\propto a^{-1}=(1+z), so Teq=T0(1+zeq)=2.725×3423=9.33×103T_{\rm eq}=T_0(1+z_{\rm eq})=2.725\times3423 = 9.33\times10^{3} K, equivalently kBTeq=0.80 eVk_BT_{\rm eq}=\mathbf{0.80\ eV}. That is the crossover the density plot marks: before z≈3400z\approx3400 the universe was radiation-dominated; after it, matter-dominated. Radiation's one extra power of aa is what puts equality at a redshift of thousands.

2 · Is the universe accelerating today? The deceleration parameter (fully worked)

Using the flat ΛCDM budget Ωm,0=0.315\Omega_{m,0}=0.315, ΩΛ,0=0.685\Omega_{\Lambda,0}=0.685, Ωr,0=9.2×10−5\Omega_{r,0}=9.2\times10^{-5} with w=0,−1,1/3w=0,-1,1/3 respectively, compute the present deceleration parameter q0≡−a¨ a/a˙2q_0\equiv-\ddot a\,a/\dot a^2 and decide the sign of a¨\ddot a.

Step 1. Divide the acceleration equation a¨/a=−4πG3∑iρi(1+3wi)\ddot a/a = -\tfrac{4\pi G}{3}\sum_i\rho_i(1+3w_i) by H2=8πG3ρtotH^2=\tfrac{8\pi G}{3}\rho_{\rm tot} (flat universe). The prefactors combine to 12\tfrac12 and each ρi/ρtot\rho_i/\rho_{\rm tot} becomes Ωi\Omega_i, giving q0=12∑iΩi(1+3wi)q_0=\tfrac12\sum_i\Omega_i(1+3w_i).

Step 2. Evaluate term by term: matter contributes Ωm(1+0)=0.315\Omega_m(1+0)=0.315; radiation contributes Ωr(1+1)=1.84×10−4\Omega_r(1+1)=1.84\times10^{-4} (negligible); dark energy contributes ΩΛ(1−3)=−2×0.685=−1.370\Omega_\Lambda(1-3)=-2\times0.685=-1.370.

Step 3. Sum and halve: q0=12(0.315+0.000184−1.370)=12(−1.055)=−0.53q_0=\tfrac12(0.315+0.000184-1.370)=\tfrac12(-1.055)=\mathbf{-0.53}. A negative q0q_0 means a¨>0\ddot a>0: the expansion is accelerating today. The dark-energy term, with its factor (1+3w)=−2(1+3w)=-2, overwhelms the positive matter contribution — a numerical demonstration that w<−1/3w<-1/3 tips the universe into acceleration. Set ΩΛ=0\Omega_\Lambda=0 and you get q0=+0.16>0q_0=+0.16>0, a decelerating matter universe; the accelerating cosmos exists only because ~70% of the budget sits below the w=−1/3w=-1/3 threshold.

Physics of the portrait

The cosmological fluid was assembled by a handful of relativists between the wars — but one figure realized that the most mysterious term in Einstein's equation is itself a fluid, a vacuum with negative pressure.

Origins · 1 minGeorges Lemaître and how this was discovered
GL
Georges Lemaître
1894–1966 · Belgian priest, physicist, and father of the expanding-universe fluid model

Georges Lemaître took holy orders and a doctorate in physics in the same years, and used both to overturn the static cosmos of his elders. In 1927, working from Einstein's equations for a homogeneous isotropic fluid, he derived an expanding-universe solution and tied it to the observed recession of galaxies — two years before Hubble's law, in a paper so obscure it went unnoticed until Eddington championed it. But his deepest contribution to the equation of state came in 1934. Where Einstein had bolted the cosmological constant Λ\Lambda onto his equations in 1917 as a device to keep the universe static, Lemaître recognized what it physically was: a fluid filling all of space with constant energy density and a pressure equal to the negative of that density, p=−ρc2p=-\rho c^2 — the value w=−1w=-1. "Everything happens as though the energy in vacuo would be different from zero," he wrote, the first clear statement that empty space itself gravitates and, because its pressure is negative, gravitates repulsively. This is precisely the dark-energy component supernova surveys confirmed sixty-four years later. The modern picture — every cosmic ingredient a perfect fluid indexed by a single number ww, the strangest of all the vacuum at w=−1w=-1 — is the direct inheritance of his insight.

History
1917
— Einstein adds the cosmological constant Λ to his equations to permit a static universe
1922
— Friedmann derives the expanding-universe equations for a homogeneous perfect fluid
1927
— Lemaître links the expansion to galaxy recession and rediscovers the dynamical solution
1934
— Lemaître and Tolman recast Λ as a vacuum energy with pressure p = −ρc², the case w = −1
1998
— Supernova surveys find the expansion accelerating, measuring dark energy at w ≈ −1

Going further

Texts. Misner, Thorne & Wheeler, Gravitation, §27 (Idealized Cosmologies) · Wald, General Relativity, Ch. 5 · Carroll, Spacetime and Geometry, §8.1–8.3 · Hartle, Gravity, Ch. 18 · Schutz, A First Course in General Relativity, Ch. 12 · Dodelson & Schmidt, Modern Cosmology, Ch. 2 · Weinberg, Cosmology, Ch. 1.

Provenance. The equation-of-state values (w=0,1/3,−1w=0,1/3,-1) and scaling laws are verified against MTW §27, Wald Ch. 5, Carroll §8.1–8.3, and Hartle Ch. 18. The continuity equation is the ν=0\nu=0 component of ∇μTμν=0\nabla_\mu T^{\mu\nu}=0 in FLRW; the acceleration equation is cross-checked against Carroll Eq. (8.23). A verifier should recompute the integration of the fluid equation (ρ∝a−3(1+w)\rho\propto a^{-3(1+w)}), the curvature cancellation giving a¨/a=−4πG3ρ(1+3w)\ddot a/a=-\tfrac{4\pi G}{3}\rho(1+3w), and the photon-gas w=1/3w=1/3 from the trace-free electromagnetic stress-energy. Numerical anchors are standard Planck-era figures: matter–radiation equality zeq≈3400z_{\rm eq}\approx3400, kBTeq≈0.8k_BT_{\rm eq}\approx0.8 eV (from Ωm,0=0.315\Omega_{m,0}=0.315, Ωr,0=9.2×10−5\Omega_{r,0}=9.2\times10^{-5}, T0=2.725T_0=2.725 K); matter–Λ equality z≈0.3z\approx0.3; q0≈−0.53q_0\approx-0.53.

Prerequisites. The Cosmological Principle · The FLRW Metric · The Scale Factor and Hubble Parameter · Comoving and Proper Coordinates. Leads to. The Friedmann Equations · The Radiation-Dominated Era · The Matter-Dominated Era · The Cosmological Constant and Dark Energy · Critical Density and the Density Parameters.

Ask the tutor: Ask me to show how the same ww that fixes a fluid's dilution law also decides, through a¨/a=−4πG3ρ(1+3w)\ddot a/a=-\tfrac{4\pi G}{3}\rho(1+3w), whether it halts the cosmos or drives it apart.

Why does pressure gravitate in general relativity but not in Newtonian gravity?

In Newton's theory the source of gravity is mass density alone; pressure never appears. In general relativity the source is the full stress-energy tensor, and its active gravitational part — the combination that governs how a cloud of matter accelerates inward — is ρc2+3p\rho c^2 + 3p, sometimes called the Tolman or Komar mass density. The factor of three counts the three spatial directions in which pressure acts as a momentum flux, each contributing equally to the curvature that bends worldlines together. Pressure is energy of motion or of fields, and in relativity all energy gravitates. This is why a hot gas gravitates more strongly than cold dust of the same rest-mass energy, and why only tension can reverse the sign and push the expansion faster.

Is the cosmological continuity equation independent of the Friedmann equations?

Not entirely — it is a consequence of them, but a useful one. The two Friedmann equations plus the continuity equation form a set of three, but only two are independent: the contracted Bianchi identity guarantees that energy-momentum conservation (∇μTμν=0\nabla_\mu T^{\mu\nu}=0, which gives continuity) follows automatically once Einstein's equation holds. In practice cosmologists use the first Friedmann equation together with continuity as the working pair, deriving the acceleration equation when needed, because continuity most transparently encodes the fluid's thermodynamics. You can regard it either as an independent physical law (local energy conservation) or as a derived identity — both views are correct.

Can the equation-of-state parameter w be less than −1?

Mathematically yes, and it has a name: phantom energy. For w<−1w<-1 the scaling law ρ∝a−3(1+w)\rho\propto a^{-3(1+w)} has a positive exponent, so the energy density grows as the universe expands — the vacuum becomes more potent the emptier space gets. This violates the dominant energy condition and leads to a "Big Rip", a finite-time singularity in which the accelerating expansion tears apart galaxies, then stars, then atoms. Current data sit at w=−1w=-1 within a few percent and cannot exclude a slightly phantom value. Most theorists regard w<−1w<-1 with suspicion because the simplest models producing it have instabilities, but it remains an active possibility the surveys are built to test.

Why is it legitimate to treat galaxies and photons as a smooth fluid at all?

Because cosmology only ever asks about scales far larger than the lumps. Averaged over a box 100 Mpc on a side — big enough to hold thousands of galaxies — the matter distribution is smooth to about one part in a hundred thousand, the level of the CMB anisotropies. On those scales the Cosmological Principle holds, and a homogeneous isotropic distribution has no preferred direction, no bulk flow relative to comoving observers, and no shear, forcing the stress-energy into the perfect-fluid form. The fluid description is not an approximation to apologize for; it is the leading term of a systematic expansion, with the lumpiness treated as small perturbations on top — the subject of structure formation.

What is the difference between ρ and ρc², and why does it matter here?

Throughout this article ρ\rho is a mass density (kg/m³) and ρc2\rho c^2 the corresponding energy density (J/m³); the equation of state p=wρc2p=w\rho c^2 compares pressure to energy density because pressure itself has units of energy density. The distinction is easy to lose because many texts set c=1c=1 and write ρ\rho for the energy density directly. Keeping c2c^2 explicit matters most in the acceleration equation, whose source ρc2+3p\rho c^2+3p adds an energy density to three times a pressure — both must share units for the sum to make sense.

How do astronomers actually measure w from the sky?

They never measure ww directly; they measure the expansion history and infer ww through the scaling law. Type Ia supernovae give luminosity distances versus redshift, hence H(z)H(z); baryon acoustic oscillations provide a fixed comoving ruler (∼\sim150 Mpc) whose apparent size calibrates the same H(z)H(z); and weak lensing measures how structure grew, which is sensitive to how dark energy suppressed that growth. Each observable folds in ρΛ(a)∝a−3(1+w)\rho_\Lambda(a)\propto a^{-3(1+w)}, so fitting them jointly pins down ww and its possible time variation, usually parametrized w(a)=w0+wa(1−a)w(a)=w_0+w_a(1-a). The precision now distinguishes w=−1w=-1 from w=−0.95w=-0.95, which is why any hint of w≠−1w\neq-1 makes headlines.

Before this

Where this leads