Equation of State in Cosmology

Every smooth ingredient of the cosmos — photons, neutrinos, baryons, cold dark matter, the mysterious dark energy — is modelled in relativistic cosmology as a perfect fluid, an idealisation the cosmological principle makes almost exact. Homogeneity and isotropy forbid any preferred direction, so there is no shear, no viscosity, no bulk flow: a comoving fluid is described completely by two numbers, its energy density ε=ρc2\varepsilon=\rho c^2 and its isotropic pressure pp. Einstein's equations need one more relation to close the system — the thermodynamic link between those two numbers. That relation is the equation of state, and for the components that matter it collapses to a single dimensionless constant ww through p=wεp=w\varepsilon. The reason ww deserves top billing is that it does double duty. Through the fluid equation it fixes how fast a species thins out as space grows, giving the dilution law ε∝a−3(1+w)\varepsilon\propto a^{-3(1+w)}: a universal a−3a^{-3} from volume dilution, times an extra a−3wa^{-3w} measuring the work pressure does during expansion. And through the acceleration equation it fixes whether that species pulls the expansion up or slows it down, because gravity in general relativity responds to ε+3p=(1+3w)ε\varepsilon+3p=(1+3w)\varepsilon. Pressureless dust (w=0w=0) simply spreads over a larger volume; radiation (w=1/3w=1/3) redshifts on top of that, fading one power faster; a cosmological constant (w=−1w=-1) refuses to dilute at all, behaving as an unchanging tension in space itself. Feed three such fluids into the Friedmann equations and the entire arc of cosmic history falls out — a radiation era, a matter era, and the dark-energy era we have just entered. The threshold w=−1/3w=-1/3 separates gravity that attracts from gravity that repels, which is why pinning down the dark-energy ww to a percent is a central goal of modern cosmology.

Key takeaways
  • The cosmological equation of state ties the isotropic pressure pp of a homogeneous cosmic fluid to its energy density ε=ρc2\varepsilon=\rho c^2 through one dimensionless number: p=w ε=wρc2p=w\,\varepsilon=w\rho c^2.
  • Substituted into the fluid (continuity) equation, ww produces the dilution law ε∝a−3(1+w)\varepsilon\propto a^{-3(1+w)}, so different components fade at different rates as the scale factor aa grows.
  • Matter has w=0w=0 (ε∝a−3\varepsilon\propto a^{-3}), radiation has w=1/3w=1/3 (ε∝a−4\varepsilon\propto a^{-4}), and a cosmological constant has w=−1w=-1 (ε\varepsilon constant); the extra power for radiation is the redshift of each photon.
  • Because radiation dilutes faster than matter and matter faster than dark energy, the Universe passes in turn through radiation-, matter-, and dark-energy-dominated eras.
  • The acceleration equation shows expansion speeds up when w<−1/3w<-1/3 — the exact line dividing ordinary attractive gravity from the repulsive gravity of dark energy.
  • Observations favour w≈−1w\approx-1 for dark energy; a phantom value w<−1w<-1 or a time-varying w(a)w(a) remain live possibilities under active test.
  • Key idea: A single number w=p/εw=p/\varepsilon fixes how fast a cosmic component dilutes (ε∝a−3(1+w)\varepsilon\propto a^{-3(1+w)}) and whether it decelerates or accelerates the expansion (threshold w=−1/3w=-1/3).

Where this lives today

The dark-energy equation of state is a headline number of observational cosmology. Type Ia supernovae calibrate the expansion history, baryon acoustic oscillations lay a "standard ruler" across redshift, and the cosmic microwave background anchors the early-time densities; together they currently place ww within a few percent of −1-1, consistent with a pure cosmological constant. The live question is whether ww is exactly −1-1 and constant, or evolves. Surveys fit the Chevallier–Polarski–Linder form w(a)=w0+wa(1−a)w(a)=w_0+w_a(1-a) and hunt for any wa≠0w_a\neq0; in 2024–2025 the DESI spectroscopic survey, combined with supernovae and the CMB, reported a mild preference for evolving dark energy that has not yet reached decisive significance. A confirmed departure from w=−1w=-1 would rule out the cosmological constant and point to a dynamical field — quintessence, whose ww ranges from −1-1 to +1+1 as the field rolls. Euclid, the Vera Rubin Observatory's LSST, and the Nancy Grace Roman telescope are built in large part to measure w(a)w(a) to the sub-percent precision that would settle whether the vacuum is static.

The physical picture

Picture a comoving box that stretches along with the Universe. When the scale factor doubles, every edge doubles and the volume grows eightfold. Any conserved particles — massive or massless — simply spread over that larger volume, so their number density falls as a−3a^{-3}. That a−3a^{-3} is the common floor shared by every species and is pure geometry. What separates one component from another is what happens to the energy each particle carries as space stretches. A cold, massive particle keeps its rest energy, so matter fades at exactly the volume rate, εm∝a−3\varepsilon_m\propto a^{-3}. A photon does not: its wavelength is stretched by the expansion, λ∝a\lambda\propto a, so each photon loses energy as a−1a^{-1}, and radiation fades one power faster, εr∝a−4\varepsilon_r\propto a^{-4}.

Figure
Doubling the scale factor multiplies the volume by 8: matter density falls by 1/8 as a⁻³, radiation falls by 1/16 as a⁻⁴ because each photon also redshifts by 1/2, and the vacuum energy is unchanged as a⁰ with w = −1
Expand the box a → 2a (volume ×8): each species thins at its own rate scale factor a expand scale factor 2a (volume ×8) matter ρ → ρ/8 (a⁻³) radiation ρ → ρ/16 (a⁻⁴) vacuum ρ → ρ (a⁰, w = −1) · ε(a) = ε₀ a⁻³⁽¹⁺ʷ⁾

The parameter ww is really a bookkeeping of how expensive pressure is for the fluid. A positive pressure (w>0w>0) pushes outward on the receding boundary and does work, draining energy faster than volume alone — which is exactly why radiation fades faster than matter. A negative pressure (w<0w<0) is a tension: the fluid resists being stretched, and the expansion does work on it. For w=−1w=-1 that inflow of work precisely balances the volume growth, so the density never changes — the vacuum gains energy at just the rate needed to fill the new space, which is why a cosmological constant behaves like a fixed property of spacetime rather than a diluting substance. The same ww decides the fate of the expansion. In general relativity pressure gravitates, entering the source of the acceleration equation with a factor of three: the active gravitational density is ε+3p=(1+3w)ε\varepsilon+3p=(1+3w)\varepsilon. Ordinary matter and radiation have 1+3w>01+3w>0 and decelerate the expansion, while a fluid with w<−1/3w<-1/3 makes that combination negative — its gravity repels — and the expansion accelerates. Track all three species at once and their crossings lay out cosmic history: the steep radiation line dominates first, then matter, then the flat dark-energy line takes over near the present day.

Figure
Energy density of each species versus scale factor on log axes: radiation ρ ∝ a⁻⁴ (slope −4), matter ρ ∝ a⁻³ (slope −3), and dark energy ρ ∝ a⁰ (flat). Radiation and matter are equal at a ≈ 2.9×10⁻⁴ (z ≈ 3400); matter and dark energy are equal at a ≈ 0.77 (z ≈ 0.3)
The punchline: three slopes, three eras — who dominates when −5 −4 −3 −2 −1 0 log₁₀ a (a = 1 today, at right edge) log₁₀ (ρ / ρ_crit,0) today radiation ρ ∝ a⁻⁴ matter ρ ∝ a⁻³ dark energy ρ ∝ a⁰ radiation = matter a ≈ 2.9×10⁻⁴, z ≈ 3400 matter = Λ a ≈ 0.77, z ≈ 0.3

The core equations

The whole subject rests on one modelling choice: that a comoving cosmic fluid is barotropic, its pressure a function of its energy density alone. For the dominant components that function is linear, and the definition below is the entire content of the equation of state.

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

where: pp — isotropic pressure of the fluid; ε\varepsilon — energy density, equal to ρc2\rho c^2; ρ\rho — mass-equivalent density ε/c2\varepsilon/c^2; ww — dimensionless equation-of-state parameter; cc — speed of light. Reading: ww is just the pressure measured in units of the energy density — dimensionless, positive or negative, and everything else in the article follows from its value. Local energy conservation in an expanding spacetime, ∇μTμν=0\nabla_\mu T^{\mu\nu}=0, gives the fluid or continuity equation. It says the energy density falls both because the volume grows and because pressure does work on the expanding boundary.

ε˙+3 a˙a (ε+p)=0(2)\dot{\varepsilon} + 3\,\frac{\dot a}{a}\,(\varepsilon + p) = 0 \qquad (2)

where: ε˙\dot\varepsilon — time rate of change of the energy density; a˙/a=H\dot a/a=H — the Hubble parameter; ε+p\varepsilon+p — the enthalpy density, which is the combination that sources dilution beyond the pure a−3a^{-3} of volume. Reading: if pressure were zero the equation would give simple volume dilution; the pressure term is precisely the "extra" physics that distinguishes species. Combining (1) and (2) with ww held constant integrates to a single power law — the most-used result in the field. It is the master formula from which every scaling is read.

ε(a)=ε0 a−3(1+w)(3)\varepsilon(a) = \varepsilon_0\,a^{-3(1+w)} \qquad (3)

where: ε0\varepsilon_0 — present-day energy density, evaluated at a0=1a_0=1; aa — scale factor normalised to unity today; ww — assumed constant. Reading: the exponent −3(1+w)-3(1+w) carries all the information. Set w=0w=0 and get a−3a^{-3} (matter); set w=1/3w=1/3 and get a−4a^{-4} (radiation); set w=−1w=-1 and the exponent vanishes, so ε\varepsilon is constant. The second Friedmann or acceleration equation shows how the same ww controls dynamics. Because pressure gravitates in general relativity, the source term is ε+3p\varepsilon+3p, not ε\varepsilon alone.

a¨a=−4πG3c2 (ε+3p)=−4πG3c2 ε (1+3w)(4)\frac{\ddot a}{a} = -\frac{4\pi G}{3c^2}\,(\varepsilon + 3p) = -\frac{4\pi G}{3c^2}\,\varepsilon\,(1+3w) \qquad (4)

where: a¨\ddot a — second time derivative of the scale factor; GG — Newton's constant; 1+3w1+3w — the sign-determining factor. Reading: since ε>0\varepsilon>0 and G>0G>0, the sign of a¨\ddot a is set entirely by −(1+3w)-(1+3w). Ordinary fluids decelerate; a fluid with 1+3w<01+3w<0 makes gravity push outward. Setting a¨=0\ddot a=0 in (4) isolates the single value of ww that marks the boundary between a slowing and a speeding-up universe.

w<−13 ⟹ a¨>0(5)w < -\tfrac{1}{3} \ \Longrightarrow\ \ddot a > 0 \qquad (5)

where: a¨>0\ddot a>0 — accelerating expansion; −1/3-1/3 — the critical equation of state at which gravity neither attracts nor repels on the largest scales. Reading: this is the operational definition of dark energy — any component negative enough in pressure to cross this line. For a barotropic fluid the adiabatic sound speed comes straight from the slope of p(ε)p(\varepsilon), which is another face of ww.

cs2=c2 ∂p∂ε∣S=w c2(w=const)(6)c_s^2 = c^2\,\frac{\partial p}{\partial \varepsilon}\bigg|_{S} = w\,c^2 \quad (w=\text{const}) \qquad (6)

where: csc_s — sound speed in the fluid; ∂p/∂ε∣S\partial p/\partial\varepsilon|_S — the derivative at fixed entropy. Reading: radiation (w=1/3w=1/3) has cs=c/3c_s=c/\sqrt{3}, the speed of sound in a photon gas; matter has cs≈0c_s\approx0. A negative ww gives a formally imaginary sound speed — a signal that dark energy is not a simple barotropic fluid but usually a field with its own effective sound speed. Because Λ\LambdaCDM contains several species at once, the total is a sum of terms, each scaling with its own ww, and this is what actually enters the Friedmann equation.

ρ(a)=ρr,0 a−4+ρm,0 a−3+ρΛ,0(7)\rho(a) = \rho_{r,0}\,a^{-4} + \rho_{m,0}\,a^{-3} + \rho_{\Lambda,0} \qquad (7)

where: ρr,0\rho_{r,0} — present radiation density (w=1/3w=1/3); ρm,0\rho_{m,0} — present matter density (w=0w=0); ρΛ,0\rho_{\Lambda,0} — present dark-energy density (w=−1w=-1). Reading: the dominant term changes with aa because the exponents differ — the ordering of exponents is the ordering of the cosmic eras.

Derivation AThe dilution law from the fluid equation

The single most consequential fact about any cosmic species is how quickly its energy density fades as the Universe grows, and that fact is not inserted by hand — it is forced by energy conservation the moment we adopt a constant equation of state. The goal here is to take the perfect-fluid continuity equation, feed in p=wεp=w\varepsilon with ww fixed, integrate cleanly, and arrive at the closed-form power law ε(a)=ε0 a−3(1+w)\varepsilon(a)=\varepsilon_0\,a^{-3(1+w)}, so that the three headline scalings of matter, radiation, and the cosmological constant can be read directly off one exponent rather than memorised separately.

Step 1. Start from the fluid equation ε˙+3(a˙/a)(ε+p)=0\dot\varepsilon+3(\dot a/a)(\varepsilon+p)=0 and substitute the equation of state, using ε+p=ε(1+w)\varepsilon+p=\varepsilon(1+w), to get ε˙+3(a˙/a) ε(1+w)=0\dot\varepsilon+3(\dot a/a)\,\varepsilon(1+w)=0. Why: this is the only place the physics of the fluid enters; once ww is chosen, the equation is a first-order linear relation between ε˙\dot\varepsilon and a˙\dot a.

Step 2. Separate the variables by dividing through by ε\varepsilon and rewriting the time derivative as a scale-factor derivative, (a˙/a) dt=da/a(\dot a/a)\,\mathrm{d}t=\mathrm{d}a/a, which gives

dεε=−3(1+w) daa.\frac{\mathrm{d}\varepsilon}{\varepsilon} = -3(1+w)\,\frac{\mathrm{d}a}{a}.

Why: separating variables is legitimate only because ww is constant, so the factor −3(1+w)-3(1+w) can be pulled outside the integral — a time-varying ww would leave an integral ∫(1+w) dln⁡a\int(1+w)\,\mathrm{d}\ln a in the exponent instead.

Step 3. Integrate both sides, giving ln⁡ε=−3(1+w)ln⁡a+const\ln\varepsilon=-3(1+w)\ln a+\text{const}. Why: both integrands are exact logarithmic derivatives, so the integration is immediate and the constant of integration is all that remains to be fixed.

Step 4. Exponentiate and pin the constant with present-day values, ε(a0)=ε0\varepsilon(a_0)=\varepsilon_0 at a0=1a_0=1, to obtain

ε(a)=ε0 a−3(1+w).\varepsilon(a) = \varepsilon_0\,a^{-3(1+w)}.

Why: normalising at today is the natural choice because all densities are quoted relative to their present values, and it makes aa dimensionless with a=1a=1 now.

Step 5. Read off the special cases from the exponent: w=0w=0 gives a−3a^{-3} (matter, pure volume dilution); w=1/3w=1/3 gives a−4a^{-4} (radiation); w=−1w=-1 gives a0a^{0}, a constant density (cosmological constant); and w=1w=1 gives the steepest causal fall-off a−6a^{-6}. Why: this is the payoff — a single formula reproduces every canonical scaling, so the whole zoo of components is encoded in one number.

Go deeper · Derivation B · 1 minWhy radiation carries the extra factor of one over a
Derivation BWhy radiation carries the extra factor of one over a

The dilution law tells us radiation scales as a−4a^{-4}, one power steeper than matter, but the formula alone does not explain where that extra power comes from. The goal of this derivation is to rebuild εr∝a−4\varepsilon_r\propto a^{-4} from microphysics — counting photons and tracking their individual energies separately — so the abstract exponent becomes a physical story of volume dilution times cosmological redshift, and then to confirm w=1/3w=1/3 independently from the kinetic pressure of a photon gas.

Step 1. Count the photons. After decoupling their number is conserved, so their density falls purely with the growing volume, nγ∝V−1∝a−3n_\gamma\propto V^{-1}\propto a^{-3}. Why: number conservation is the same geometric dilution matter feels, so this factor is the common a−3a^{-3} floor and cannot be what distinguishes radiation.

Step 2. Track each photon's energy. Its wavelength is stretched by the expansion, λ∝a\lambda\propto a, so its energy Eγ=hc/λ∝a−1E_\gamma=hc/\lambda\propto a^{-1}. Why: this is the cosmological redshift acting on a single quantum — the piece that has no analogue for a cold massive particle, which keeps its rest energy.

Step 3. Multiply the two effects: εr=nγEγ∝a−3⋅a−1=a−4\varepsilon_r=n_\gamma E_\gamma\propto a^{-3}\cdot a^{-1}=a^{-4}, matching the dilution law at w=1/3w=1/3 since −3(1+1/3)=−4-3(1+1/3)=-4. Why: energy density is number density times energy per particle, so the two independent powers simply add in the exponent.

Step 4. Cross-check the value of ww from pressure. Kinetic theory of an isotropic photon gas gives p=ε/3p=\varepsilon/3 directly, because a massless gas exerts a pressure equal to one third of its energy density. Why: an entirely separate argument — radiation pressure — lands on w=1/3w=1/3, confirming that the fluid picture and the microscopic picture are the same physics seen two ways.

Go deeper · Derivation C · 1 minThe acceleration threshold at w equals minus one third
Derivation CThe acceleration threshold at w equals minus one third

Cosmic acceleration is not a small correction but a change of sign, and the equation of state controls that sign exactly. The goal here is to obtain the acceleration equation from the two Friedmann equations, expose the source term ε+3p\varepsilon+3p in which pressure gravitates with a factor of three, and substitute the equation of state to show that the expansion switches from decelerating to accelerating precisely at w=−1/3w=-1/3 — the dividing line any candidate for dark energy must cross.

Step 1. Write the first Friedmann equation as the energy constraint,

(a˙a)2=8πG3c2 ε−kc2a2,\left(\frac{\dot a}{a}\right)^2 = \frac{8\pi G}{3c^2}\,\varepsilon - \frac{kc^2}{a^2},

with kk the spatial curvature. Why: this is the 0000 Einstein equation and fixes the expansion rate from the energy density; differentiating it is the route to the acceleration.

Step 2. Differentiate with respect to time and eliminate ε˙\dot\varepsilon using the fluid equation ε˙=−3H(ε+p)\dot\varepsilon=-3H(\varepsilon+p), which makes the curvature term drop out and leaves

a¨a=−4πG3c2 (ε+3p).\frac{\ddot a}{a} = -\frac{4\pi G}{3c^2}\,(\varepsilon+3p).

Why: the curvature term cancels because the same kk appears in both the constraint and its derivative, isolating a clean statement about a¨\ddot a that depends only on the fluid content.

Step 3. Substitute the equation of state p=wεp=w\varepsilon to get

a¨a=−4πG3c2 ε (1+3w).\frac{\ddot a}{a} = -\frac{4\pi G}{3c^2}\,\varepsilon\,(1+3w).

Why: this reduces the dynamics to a single parameter — the entire question of the Universe's fate collapses onto the sign of 1+3w1+3w.

Step 4. Analyse the sign. With ε>0\varepsilon>0 and G>0G>0, the expansion decelerates for w>−1/3w>-1/3, coasts at w=−1/3w=-1/3 (a¨=0\ddot a=0), and accelerates for w<−1/3w<-1/3. Why: at the threshold the attractive gravity of the energy density is exactly cancelled by the repulsive contribution of its negative pressure, which is why w=−1/3w=-1/3 is the defining line for dark energy.

Limiting cases

The power of the linear equation of state is that a short table of ww values captures every fluid the standard model needs, plus the effective fluids that mimic curvature and stiff matter.

CaseResultMeaning
Non-relativistic matter (dust)w=0,  ε∝a−3w=0,\ \ \varepsilon\propto a^{-3}Cold dark matter and baryons; thermal pressure negligible next to rest energy, so pure volume dilution.
Radiation / relativistic speciesw=13,  ε∝a−4w=\tfrac{1}{3},\ \ \varepsilon\propto a^{-4}Photons and (early) neutrinos; the extra a−1a^{-1} is the cosmological redshift of each quantum.
Cosmological constant / vacuumw=−1,  ε=constw=-1,\ \ \varepsilon=\text{const}Density independent of aa; drives constant HH and de Sitter expansion — the endpoint of Λ\LambdaCDM.
Spatial curvature (effective)w=−13,  εk∝a−2w=-\tfrac{1}{3},\ \ \varepsilon_k\propto a^{-2}Not a real fluid, but the curvature term scales like one; sits exactly on the acceleration boundary.
Stiff / kination fluidw=1,  ε∝a−6w=1,\ \ \varepsilon\propto a^{-6}The steepest causal equation of state (cs=cc_s=c); the kinetic energy of a fast-rolling scalar field behaves this way.
Phantom energyw<−1,  ε growsw<-1,\ \ \varepsilon\ \text{grows}Density increases with expansion; if real, ends in a finite-time "Big Rip". Allowed by some data fits, physically pathological.

Worked examples

1 · When did radiation give way to matter? (fully worked)

Take the present density parameters Ωm,0=0.31\Omega_{m,0}=0.31 (matter, w=0w=0) and Ωr,0=9.1×10−5\Omega_{r,0}=9.1\times10^{-5} (radiation, w=1/3w=1/3, photons plus relativistic neutrinos). Find the scale factor and redshift of matter–radiation equality, and the photon temperature there, given T0=2.725T_0=2.725 K.

Step 1. Equality is where the two densities cross. Using the dilution law for each, Ωm,0 a−3=Ωr,0 a−4\Omega_{m,0}\,a^{-3}=\Omega_{r,0}\,a^{-4}. The a−3a^{-3} cancels on both sides, leaving Ωm,0=Ωr,0 a−1\Omega_{m,0}=\Omega_{r,0}\,a^{-1}, so

aeq=Ωr,0Ωm,0=9.1×10−50.31=2.9×10−4.a_{\rm eq}=\frac{\Omega_{r,0}}{\Omega_{m,0}}=\frac{9.1\times10^{-5}}{0.31}=\mathbf{2.9\times10^{-4}}.

Step 2. Convert to redshift with 1+z=1/a1+z=1/a: 1+zeq=1/(2.9×10−4)=34071+z_{\rm eq}=1/(2.9\times10^{-4})=3407, so zeq≈3.4×103\mathbf{z_{\rm eq}\approx3.4\times10^{3}}. Before this redshift radiation dominated the expansion; after it, matter did — the transition that lets structure begin to grow.

Step 3. The photon temperature scales as T∝a−1=1+zT\propto a^{-1}=1+z, so Teq=T0(1+zeq)=2.725×3407=9.3×103 KT_{\rm eq}=T_0(1+z_{\rm eq})=2.725\times3407=\mathbf{9.3\times10^{3}\ K}, corresponding to a thermal energy kBTeq≈0.80k_BT_{\rm eq}\approx0.80 eV. Equality falls at aeq≈2.9×10−4a_{\rm eq}\approx2.9\times10^{-4}, zeq≈3400z_{\rm eq}\approx3400, Teq≈9300T_{\rm eq}\approx9300 K — comfortably before recombination at z≈1100z\approx1100, so the Universe was already matter-dominated when the CMB was released.

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

Take a flat Λ\LambdaCDM model with Ωm,0=0.31\Omega_{m,0}=0.31 and ΩΛ,0=0.69\Omega_{\Lambda,0}=0.69, neglecting radiation. Find the present deceleration parameter q0q_0, and the redshift at which acceleration began.

Step 1. The deceleration parameter is q≡−a¨ a/a˙2q\equiv-\ddot a\,a/\dot a^{2}, and for a mix of fluids the acceleration equation gives q=12∑iΩi(1+3wi)q=\tfrac12\sum_i\Omega_i(1+3w_i). Today, with matter (w=0w=0) and Λ\Lambda (w=−1w=-1),

q0=12[Ωm,0(1)+ΩΛ,0(−2)]=12[0.31−1.38]=−0.54.q_0=\tfrac12\big[\Omega_{m,0}(1)+\Omega_{\Lambda,0}(-2)\big]=\tfrac12\big[0.31-1.38\big]=\mathbf{-0.54}.

The negative value is the direct statement that the Universe is accelerating now.

Step 2. Acceleration begins when a¨=0\ddot a=0, i.e. when the total ε+3p\varepsilon+3p changes sign: ρm−2ρΛ=0\rho_m-2\rho_\Lambda=0. Using the scalings, Ωm,0 a−3=2 ΩΛ,0\Omega_{m,0}\,a^{-3}=2\,\Omega_{\Lambda,0}, so

aacc=(Ωm,02 ΩΛ,0)1/3=(0.311.38)1/3=0.61.a_{\rm acc}=\left(\frac{\Omega_{m,0}}{2\,\Omega_{\Lambda,0}}\right)^{1/3}=\left(\frac{0.31}{1.38}\right)^{1/3}=0.61.

Step 3. In redshift, 1+zacc=1/0.61=1.651+z_{\rm acc}=1/0.61=1.65, so zacc≈0.6\mathbf{z_{\rm acc}\approx0.6}. The cosmos was decelerating until z≈0.6z\approx0.6 and has been accelerating since, with q0≈−0.54q_0\approx-0.54 today. Note that matter and dark energy are only equal somewhat later, at z≈0.3z\approx0.3 (from Ωm,0a−3=ΩΛ,0\Omega_{m,0}a^{-3}=\Omega_{\Lambda,0}) — acceleration starts before Λ\Lambda formally dominates, because pressure enters with a factor of three.

3 · Reading ww off a scaling law (quick)

A hypothetical component dilutes as ε∝a−5\varepsilon\propto a^{-5}. What is its equation of state? Match exponents in the dilution law: −3(1+w)=−5-3(1+w)=-5, so 1+w=5/31+w=5/3 and w=2/3\mathbf{w=2/3}. Its sound speed is cs=2/3 c≈0.82 c<cc_s=\sqrt{2/3}\,c\approx0.82\,c<c, so it is causal, lying between radiation (w=1/3w=1/3) and the stiff limit (w=1w=1); it decelerates the expansion since w>−1/3w>-1/3.

Physics of the portrait

The equation of state is where thermodynamics meets gravitation, and no one saw that junction more clearly, or earlier, than the Belgian priest-physicist who first read Einstein's cosmological constant as a pressure.

Origins · 1 minGeorges Lemaître and how this was discovered
GL
Georges Lemaître
1894–1966 · relativistic cosmology and the vacuum equation of state

Georges Lemaître trained as a civil engineer, served in the First World War, was ordained a Catholic priest, and then reshaped cosmology. In 1927, two years before Hubble's paper, he derived an expanding-universe solution of Einstein's equations and predicted a linear velocity–distance relation from the redshifts of galaxies — the relation now called the Hubble–Lemaître law. In 1931 he proposed that the expansion traced back to a dense initial state, the "primeval atom," the conceptual seed of the Big Bang. But his sharpest contribution to this subject came in 1934. Where Einstein had introduced the cosmological constant Λ\Lambda in 1917 as a reluctant geometric device to keep the Universe static, Lemaître reinterpreted it physically: Λ\Lambda is the energy of the vacuum, with density ρvac=Λc2/8πG\rho_{\rm vac}=\Lambda c^2/8\pi G and — crucially — a pressure equal to minus that energy density. In his words, one must "associate a pressure p=−ρc2p=-\rho c^2 to the density of energy of vacuum." That is exactly the equation of state w=−1w=-1, written down six decades before observations would demand it. Lemaître thus supplied the modern reading of dark energy: not empty geometry but a fluid with strongly negative pressure whose gravity repels. He lived to hear, days before his death in 1966, that the cosmic microwave background had been found — evidence for the hot, dense beginning he had championed.

History
1917
— Einstein inserts the cosmological constant Λ to permit a static universe
1922
— Friedmann derives the expanding-universe equations from general relativity
1927
— Lemaître links the redshifts of galaxies to cosmic expansion
1934
— Lemaître identifies Λ as vacuum energy with pressure p = −ρc², the case w = −1
1998
— supernova surveys find cosmic acceleration, requiring a component with w < −1/3

Going further

Texts. Ryden, Introduction to Cosmology, 2nd ed., ch. 5 (fluid and acceleration equations) · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., ch. 29 · Dodelson, Modern Cosmology, 2nd ed., ch. 2 · Mo, van den Bosch & White, Galaxy Formation and Evolution, ch. 3 · Peebles, Principles of Physical Cosmology, ch. 5. The key results of this article are: the dilution law ε∝a−3(1+w)\varepsilon\propto a^{-3(1+w)} from integrating the fluid equation with constant ww; the three canonical scalings a−3a^{-3}, a−4a^{-4}, and constant for w=0,1/3,−1w=0,1/3,-1; the acceleration threshold w=−1/3w=-1/3 from equation (4); and the sound speed cs2=wc2c_s^2=wc^2, giving c/3c/\sqrt{3} for radiation. The worked numbers use standard Λ\LambdaCDM parameters Ωm,0=0.31\Omega_{m,0}=0.31, ΩΛ,0=0.69\Omega_{\Lambda,0}=0.69, Ωr,0=9.1×10−5\Omega_{r,0}=9.1\times10^{-5}, yielding matter–radiation equality zeq≈3400z_{\rm eq}\approx3400, acceleration onset zacc≈0.6z_{\rm acc}\approx0.6, and present deceleration parameter q0≈−0.54q_0\approx-0.54. The identification of vacuum pressure as p=−ρc2p=-\rho c^2 traces to Lemaître (1934).

Prerequisites. The Fluid and Acceleration Equations · The Scale Factor and Cosmic Expansion · The FLRW Metric · Homogeneity and Isotropy in Cosmology. Leads to. The Friedmann Equations · Critical Density and the Density Parameters · Cosmological Redshift.

Ask the tutor: Ask me to slide ww from +1+1 down to −2-2 and watch the dilution law, the sound speed, and the fate of the expansion change at each step.

Why does a single number w capture so much physics?

Because the cosmological principle strips a cosmic fluid down to two quantities — energy density and isotropic pressure — and a barotropic relation ties them together. With no shear, heat flow, or bulk velocity, p(ε)p(\varepsilon) is the entire thermodynamic content. For the dominant components that relation is linear, p=wεp=w\varepsilon, and ww enters both governing equations: the fluid equation (how fast the density dilutes) and the acceleration equation (whether gravity attracts or repels). Two of the most important questions about a species — how it fades, and what it does to the expansion — are answered by the same number.

How can the vacuum keep a constant energy density while space expands?

Follow the fluid equation. It reads ε˙=−3H(ε+p)\dot\varepsilon=-3H(\varepsilon+p), and for w=−1w=-1 the pressure is p=−εp=-\varepsilon, so ε+p=0\varepsilon+p=0 and therefore ε˙=0\dot\varepsilon=0 exactly. Physically, the expansion does work against the vacuum's tension, and that work is precisely what is needed to fill each new parcel of volume at the same density. Energy does not appear from nowhere in any suspicious sense — in general relativity total energy is not globally conserved in an expanding spacetime, and the fluid equation is the correct local statement of energy–momentum conservation. Constant density is what it demands when pressure equals minus the energy density.

Does negative pressure mean the dark energy has negative mass?

No. The energy density of dark energy is positive; only its pressure is negative. The confusion comes from the acceleration equation, whose source is ε+3p\varepsilon+3p, not ε\varepsilon. For a cosmological constant this is ε+3(−ε)=−2ε\varepsilon+3(-\varepsilon)=-2\varepsilon, negative even though ε>0\varepsilon>0. So the "active gravitational density" driving the expansion is negative, producing repulsion — but it is built from a positive energy and a more-negative pressure, not from any negative mass.

Why is there a factor of three in front of the pressure?

It comes from the structure of the stress-energy tensor in Einstein's equations. The acceleration equation is sourced by ρ+3p/c2\rho+3p/c^2 because the mass density ρ\rho contributes once while the isotropic pressure contributes through all three spatial components — one factor per spatial direction. This is a purely relativistic effect: in Newtonian gravity pressure does not gravitate at all. It is why radiation, despite a modest energy density, decelerated the early Universe strongly, and why negative pressure can flip the sign of gravity.

What is the difference between w = −1/3 and w = −1?

They mark two different thresholds. The value w=−1/3w=-1/3 is the acceleration boundary: any fluid with w<−1/3w<-1/3 makes ε+3p<0\varepsilon+3p<0 and drives a¨>0\ddot a>0 — the minimum requirement to be called dark energy. The value w=−1w=-1 is the cosmological constant, the special case where the density is exactly constant and the equation of state matches vacuum energy. A component could accelerate the Universe with, say, w=−0.8w=-0.8 (still evolving, still diluting slowly) without being a true cosmological constant. Distinguishing "w<−1/3w<-1/3" from "w=−1w=-1 exactly" is the whole observational program of dark-energy surveys.

Can w change over cosmic time, and how would we know?

Yes — nothing forbids a time-varying equation of state, and a rolling scalar field (quintessence) generically produces one. When ww is not constant the dilution law generalises to ε∝exp⁡ ⁣[3 ⁣∫(1+w(a)) dln⁡a]\varepsilon\propto\exp\!\big[3\!\int(1+w(a))\,\mathrm{d}\ln a\big], reducing to the power law only when ww is fixed. Surveys parametrise the evolution as w(a)=w0+wa(1−a)w(a)=w_0+w_a(1-a) and measure (w0,wa)(w_0,w_a) from the expansion history in supernovae, baryon acoustic oscillations, and the CMB. A significant wa≠0w_a\neq0 would prove dark energy evolves and rule out a pure cosmological constant. As of the mid-2020s the evidence is suggestive but not decisive, which is why the next generation of surveys is being built.

Before this

Where this leads