Critical Density and Density Parameters

Cosmology's most consequential number is a ratio. Divide the actual energy density of the universe by one particular reference value and the answer tells you the global shape of space, the ultimate fate of the cosmos, and the entire inventory of what the universe is made of. That reference value is the critical density $\rho_c$ : the total energy density that would make space exactly flat at a given moment. It is not a constant of nature but a moving target — because $\rho_c = 3H^2/8\pi G$ tracks the square of the expansion rate, it was colossal in the early fireball and has fallen to about five hydrogen atoms per cubic metre today. Normalizing every energy component by $\rho_c$ converts the first Friedmann equation into a single flat accounting identity: the dimensionless density parameters $\Omega_i = \rho_i/\rho_c$ , plus one curvature term, must add to one. This is the language in which every cosmological result is now quoted — $\Omega_\Lambda \approx 0.69$ , $\Omega_m \approx 0.31$ , $\Omega_b \approx 0.049$ , $\Omega_k \approx 0$ . The sign of $\Omega - 1$ is the sign of the spatial curvature $k$ , so weighing the universe is measuring its geometry. And because each component dilutes at its own rate — matter as $a^{-3}$ , radiation as $a^{-4}$ , the cosmological constant not at all — the $\Omega_i$ shift over cosmic history, carving out radiation, matter, and dark-energy eras. That the total sits so implausibly close to one is the flatness problem, whose resolution is the leading motive for cosmic inflation.

Key takeaways

The critical density $\rho_c = 3H^2/8\pi G$ is the total energy density that makes space exactly flat ( $k=0$ ) at a given epoch; it scales as $H^2$ , so it evolves and is not a universal constant.
Dividing any component's density by $\rho_c$ gives its dimensionless density parameter $\Omega_i = \rho_i/\rho_c$ , the universal currency of cosmology.
The Friedmann equation becomes the exact sum rule $\Omega_m + \Omega_r + \Omega_\Lambda + \Omega_k = 1$ , with curvature entering as a formal density $\Omega_k = -kc^2/(a^2 H^2)$ .
One number fixes the geometry: $\Omega>1$ (closed, $k=+1$ ), $\Omega=1$ (flat, $k=0$ ), $\Omega<1$ (open, $k=-1$ ).
Today $\rho_{c,0} = 1.88\times10^{-26}\,h^2 \approx 8.5\times10^{-27}\ \mathrm{kg\,m^{-3}}$ (h ≈ 0.67) — about five protons per cubic metre — and observations give $\Omega_0 = 1$ to sub-percent precision, $|\Omega_k| \lesssim 0.005$ .
Key idea: Divide the Friedmann equation by $H^2$ : every energy component becomes a dimensionless fraction $\Omega_i = \rho_i/\rho_c$ of the flat-universe critical density, and curvature becomes one more term, so that $\sum_i \Omega_i + \Omega_k = 1$ identically.

Where this lives today

The density parameters are the working output of every precision-cosmology experiment. Planck's map of the cosmic microwave background fixes the whole set at once: the angular scale of the acoustic peaks measures $\Omega_k$ (pinning the first peak at multipole $\ell \approx 220$ , giving $\Omega_k = 0.001 \pm 0.002$ ), the peak heights separate $\Omega_b h^2 \approx 0.0224$ from the total matter, and the late-time decay of gravitational potentials constrains $\Omega_\Lambda$ . Galaxy surveys such as DESI and Rubin/LSST cross-check the same numbers through baryon acoustic oscillations — a geometry-dependent standard ruler — and weak lensing, which weighs $\Omega_m$ directly. The persistent Hubble tension, a 4–5σ disagreement between the local $H_0 \approx 73$ and the CMB-inferred $H_0 \approx 67$ , matters here because $\rho_c \propto H_0^2$ : the two values imply critical densities differing by nearly 20%, so the tension propagates into every derived $\Omega$ . And the nature of the dominant term $\Omega_\Lambda$ — a true cosmological constant ( $w=-1$ ) or a dynamical dark energy ( $w \ne -1$ ) — is the central target of the next decade of surveys, all phrased in the language of density parameters.

The physical picture

Throw a ball straight up. Whether it falls back or escapes forever is decided by one comparison: is its launch speed above or below escape velocity? The expanding universe poses exactly this question on the largest scale. Expansion is the outward throw; the gravity of all the universe's energy is the pull-back. For a given expansion rate $H$ there is one borderline density at which the throw is precisely at escape speed — the universe decelerates forever but never quite recollapses. That borderline is the critical density $\rho_c = 3H^2/8\pi G$ . Denser than critical ( $\Omega>1$ ) and the universe is bound, recollapsing to a Big Crunch if matter is all there is; less dense ( $\Omega<1$ ) and it coasts to infinity; exactly critical ( $\Omega=1$ ) is the razor's edge, the only flat case.

Figure

The fate fork for a matter-only universe. Ω = ρ/ρ_c decides destiny: Ω over 1 (closed, red) recollapses to a Big Crunch, Ω = 1 (flat, green, ρ = ρ_c = 3H²/8πG) coasts asymptotically to rest, Ω under 1 (open, blue) expands forever. Today Ω_total ≈ 1 to better than one percent.

The catch is that $\rho_c$ itself falls as the universe expands and $H$ drops, so the absolute density is a poor guide — what matters is always the ratio $\Omega = \rho/\rho_c$ . Both $\rho$ and $\rho_c$ can plummet together across billions of years while the geometry (the sign of $k$ ) stays fixed, which is why cosmologists quote $\Omega$ 's rather than kilograms. The bookkeeping becomes perfectly clean once curvature itself is promoted into a fictitious "fluid" $\Omega_k = -kc^2/(a^2 H^2)$ : the Friedmann constraint is then just the statement that all the fractions — real components plus curvature — sum to one. Because the components dilute at different rates, their fractions cross over as space expands. Radiation, redshifting as well as diluting, thins fastest ( $a^{-4}$ ) and ruled the earliest universe; pressureless matter ( $a^{-3}$ ) overtook it at $z \approx 3400$ ; the cosmological constant does not dilute at all ( $a^0$ ) and inevitably came to dominate — the takeover happening only at $z \approx 0.3$ , essentially yesterday in cosmic terms. The second figure is the inventory a cosmologist reaches for: each component's density against the scale factor, with the two crossover epochs that partition cosmic history.

Figure

Cosmic inventory ρᵢ(a) on log axes. Radiation dilutes as a⁻⁴ (red), matter as a⁻³ (blue), dark energy is constant (amber). They cross at radiation–matter equality a ≈ 3×10⁻⁴ (z ≈ 3400) and matter–Λ equality a ≈ 0.77 (z ≈ 0.3), so the universe ran radiation-dominated → matter-dominated → dark-energy-dominated.

The core equations

The starting point is the first Friedmann equation. Setting its curvature term to zero and solving for the density that keeps the equality defines the critical density — the density a flat universe expanding at rate $H$ must have. Everything downstream is a re-expression of this one condition.

\rho_c = \frac{3H^2}{8\pi G} \qquad (1)

where: $\rho_c$ — critical energy density (expressed as a mass density) at a given epoch; $H = \dot a/a$ — Hubble parameter; $G$ — Newton's gravitational constant. Read it as a definition tied to this moment's expansion rate: because $H$ evolves, so does $\rho_c$ , and it was far larger in the past. Inserting today's Hubble constant gives the present-day value, the anchor against which every $\Omega$ is measured. Writing $H_0 = 100\,h\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$ keeps the residual uncertainty in $H_0$ visible through the factor $h^2$ .

\rho_{c,0} = \frac{3H_0^2}{8\pi G} \approx 1.88\times10^{-26}\,h^2\ \mathrm{kg\,m^{-3}} \qquad (2)

where: $H_0$ — Hubble constant today; $h \approx 0.67$ – $0.70$ — dimensionless Hubble parameter. For $h=0.7$ this is $9.2\times10^{-27}\ \mathrm{kg\,m^{-3}}$ , roughly five protons per cubic metre — the emptiness against which "flat" is measured. In astronomer's units it reads $\rho_{c,0} = 2.775\times10^{11}\,h^2\ M_\odot\,\mathrm{Mpc^{-3}}$ . The critical density buys the definition that makes cosmology dimensionless: divide any component's density by $\rho_c$ to get its density parameter. A component with $\Omega_i > 1$ would, by itself, overclose the universe.

\Omega_i \equiv \frac{\rho_i}{\rho_c} = \frac{8\pi G\,\rho_i}{3H^2} \qquad (3)

where: $\Omega_i$ — dimensionless density parameter of component $i$ ; $\rho_i$ — its energy density as a mass density. The measured budget is $\Omega_\Lambda \approx 0.69$ , $\Omega_m \approx 0.31$ (of which baryons $\Omega_b \approx 0.049$ ), $\Omega_r \approx 9\times10^{-5}$ . Curvature is not a substance, but it enters the Friedmann equation on the same footing as one, so it is convenient to package it as a formal density parameter with a deliberate minus sign that flips the sign relative to $k$ .

\Omega_k \equiv -\frac{kc^2}{a^2 H^2} \qquad (4)

where: $\Omega_k$ — curvature "density" parameter; $k \in \{-1,0,+1\}$ — curvature constant; $a$ — scale factor; $c$ — speed of light. The sign convention means positive curvature ( $k=+1$ , closed) gives $\Omega_k < 0$ , so that $\Omega_{\rm tot} > 1$ — density above critical. Dividing the Friedmann equation through by $H^2$ now collapses it into a single algebraic constraint. This is the master equation of the whole subject: a statement that the components plus curvature partition unity.

\Omega_m + \Omega_r + \Omega_\Lambda + \Omega_k = 1 \qquad (5)

where: $\Omega_m$ — matter (baryonic + dark), $w=0$ ; $\Omega_r$ — radiation (photons + relativistic neutrinos), $w=1/3$ ; $\Omega_\Lambda$ — cosmological constant / dark energy, $w=-1$ ; $\Omega_k$ — curvature. Equivalently $1 - \Omega_{\rm tot} = \Omega_k$ , so measuring the total density measures the geometry. To turn the constraint into an expansion history, feed in how each component scales with $a$ . Every $\Omega_{i,0}$ is a present-day value ( $a_0 = 1$ ), and the whole $a$ -dependence is carried by the powers of $a$ .

H^2(a) = H_0^2\!\left[\,\Omega_{r,0}\,a^{-4} + \Omega_{m,0}\,a^{-3} + \Omega_{k,0}\,a^{-2} + \Omega_{\Lambda,0}\,\right] \qquad (6)

where: $\Omega_{i,0}$ — present-day density parameters. This function, often written $H(a) = H_0\,E(a)$ , is the workhorse for computing distances, ages, and growth in the standard $\Lambda$ CDM model. Each power of $a$ in (6) comes from the equation of state of that component, via the fluid (continuity) equation $\rho_i \propto a^{-3(1+w_i)}$ : radiation ( $w=1/3$ ) gives $a^{-4}$ , matter ( $w=0$ ) gives $a^{-3}$ , a cosmological constant ( $w=-1$ ) gives $a^{0}$ , and the formal curvature "fluid" ( $w=-1/3$ ) gives $a^{-2}$ . Finally, rewriting the curvature term exposes why flatness is dynamically fragile: the deviation of the total density from critical is governed by $\dot a$ , which shrinks during decelerated expansion.

1 - \Omega_{\rm tot}(t) = -\frac{kc^2}{a^2 H^2} = -\frac{kc^2}{\dot a^2} = \Omega_k(t) \qquad (7)

where: $\Omega_{\rm tot}(t)$ — total matter-plus-energy density parameter at time $t$ ; $\dot a = aH$ . Because $kc^2$ is a fixed constant, anything that decreases $\dot a$ drives $\Omega_{\rm tot}$ away from one — the flatness problem, developed below.

Derivation AFrom the Friedmann equation to the density-parameter sum rule

The goal is to show that the first Friedmann equation — apparently a differential equation relating the expansion rate to energy density and curvature — is secretly a simple algebraic identity once every term is measured in units of the critical density. We divide through by $H^2$ , recognize $\rho_c$ in the leading coefficient, and repackage the leftover curvature term as a formal density parameter, arriving at the constraint that all the dimensionless fractions sum to exactly one at every epoch.

Step 1. Start from the first Friedmann equation, $H^2 = \frac{8\pi G}{3}\rho_{\rm tot} - \frac{kc^2}{a^2}$ , where $\rho_{\rm tot} = \rho_m + \rho_r + \rho_\Lambda$ collects every physical component (the cosmological constant written as a density $\rho_\Lambda = \Lambda c^2/8\pi G$ ). Why: this is the exact dynamical equation from general relativity applied to a homogeneous, isotropic universe; we are only going to rearrange it, not approximate.

Step 2. Divide both sides by $H^2$ to make every term dimensionless: $1 = \frac{8\pi G}{3H^2}\rho_{\rm tot} - \frac{kc^2}{a^2 H^2}$ . Why: dividing by $H^2$ is exactly dividing by the thing that sets the critical density, so the natural units for density here are units of $\rho_c$ — the move that turns kilograms into pure numbers.

Step 3. Recognize the leading coefficient as an inverse critical density: since $\rho_c = 3H^2/8\pi G$ , the first term is $\rho_{\rm tot}/\rho_c = \Omega_m + \Omega_r + \Omega_\Lambda$ . Why: this is the payoff of defining $\rho_c$ — the awkward $8\pi G/3H^2$ prefactor is precisely $1/\rho_c$ , so the total density becomes a sum of clean fractions $\Omega_i$ .

Step 4. Define the remaining curvature term as its own parameter, $\Omega_k \equiv -kc^2/(a^2 H^2)$ , choosing the sign so that positive $\Omega_k$ (open, $k=-1$ ) means sub-critical density. Why: with this definition curvature is treated on the same footing as a real component, and the equation reads as a single tidy conservation of unity.

Step 5. Collect the pieces to obtain $\Omega_m + \Omega_r + \Omega_\Lambda + \Omega_k = 1$ , an identity that holds at every epoch provided each $\Omega_i$ is evaluated at that epoch. Why: it is exact, not an approximation — the entire content of the Friedmann equation is compressed into the statement that the fractions always add to one, which is why the deviation of the visible total from one directly measures the curvature.

Go deeper · Derivation B · 1 minThe present-day critical density in numbers

Derivation BThe present-day critical density in numbers

The goal is to turn the compact formula $\rho_{c,0} = 3H_0^2/8\pi G$ into a concrete SI number and then a particle count, so that "the universe is nearly flat" acquires a physical scale. We convert the Hubble constant into inverse seconds, substitute the fundamental constants, and divide by a proton mass to read the critical density as atoms per cubic metre.

Step 1. Convert $H_0$ to SI. Taking $H_0 = 70\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$ and $1\,\mathrm{Mpc} = 3.086\times10^{22}\ \mathrm{m}$ gives $H_0 = (7.0\times10^{4})/(3.086\times10^{22}) = 2.27\times10^{-18}\ \mathrm{s^{-1}}$ . Why: the formula needs $H_0$ in inverse seconds so that $\rho_c$ comes out in kilograms per cubic metre; the megaparsec is an astronomer's convenience, not an SI unit.

Step 2. Square it and assemble the constants: $\rho_{c,0} = \dfrac{3(2.27\times10^{-18})^2}{8\pi(6.674\times10^{-11})}$ , with numerator $3H_0^2 = 1.54\times10^{-35}$ and denominator $8\pi G = 1.68\times10^{-9}$ . Why: keeping the numerator and denominator separate makes the arithmetic auditable and the $h^2$ scaling transparent — every factor of $H_0$ enters squared.

Step 3. Evaluate: $\rho_{c,0} \approx 9.2\times10^{-27}\ \mathrm{kg\,m^{-3}}$ for $h=0.7$ (and $8.5\times10^{-27}$ for $h=0.674$ ), matching $1.88\times10^{-26}\,h^2\ \mathrm{kg\,m^{-3}}$ . Why: this is the anchor value — the same $\rho_{c,0}$ that converts every measured physical density into a density parameter $\Omega_i$ .

Step 4. Convert to a number density by dividing by the proton mass: $n = \rho_{c,0}/m_p = (9.2\times10^{-27})/(1.67\times10^{-27}) \approx 5.5\ \mathrm{m^{-3}}$ . Why: about five or six hydrogen atoms per cubic metre is a strikingly small density — and most of it is not atoms at all but dark energy and dark matter, which is what makes the near-perfect flatness so remarkable.

Go deeper · Derivation C · 1 minWhy cosmic flatness is unstable — the flatness problem

Derivation CWhy cosmic flatness is unstable — the flatness problem

The goal is to show that in an ordinary matter- or radiation-dominated universe $\Omega_{\rm tot} = 1$ is an unstable equilibrium: any tiny early deviation from flatness is amplified with time, so today's near-flatness implies an absurd fine-tuning at early epochs. We track the deviation $1 - \Omega_{\rm tot}$ using the constancy of $kc^2$ , show it grows whenever $\dot a$ falls, quantify the required tuning, and then show how accelerated expansion reverses the instability.

Step 1. Rewrite the deviation from flatness using $aH = \dot a$ : $1 - \Omega_{\rm tot} = \Omega_k = -kc^2/(a^2 H^2) = -kc^2/\dot a^2$ . Why: expressing everything through $\dot a$ isolates the only time-dependent quantity, because $kc^2$ is a fixed constant set once and for all by the geometry.

Step 2. In a decelerating universe $\dot a$ decreases with time. For radiation $a \propto t^{1/2}$ so $\dot a \propto t^{-1/2}$ , and for matter $a \propto t^{2/3}$ so $\dot a \propto t^{-1/3}$ — both fall as $t$ grows. Why: gravity decelerates ordinary expansion, and a shrinking $\dot a$ sits in the denominator of the deviation, so it drives $|1-\Omega_{\rm tot}|$ upward.

Step 3. Therefore $|1 - \Omega_{\rm tot}| = |k|c^2/\dot a^2$ grows monotonically: any nonzero deviation is amplified, so $\Omega_{\rm tot}=1$ is an unstable fixed point. Why: this is the crux — flatness is not self-maintaining but self-destroying under deceleration, so finding the universe flat today demands an explanation.

Step 4. Run the number backward. Today $|\Omega_0 - 1| \lesssim 0.005$ , which extrapolated to the Planck era requires $|\Omega - 1| \lesssim 10^{-60}$ — a fine-tuning of sixty decimal places. Why: no natural initial condition explains such a coincidence, so the near-flatness is a genuine problem, not a triviality.

Step 5. Inflation reverses the sign of $\ddot a$ : during accelerated expansion $\dot a$ grows, so $|1-\Omega_{\rm tot}| = |k|c^2/\dot a^2$ shrinks, driving $\Omega_{\rm tot} \to 1$ . Why: a burst of accelerated expansion makes flatness a dynamical attractor rather than a fine-tuning, which is why the flatness problem is one of the founding motivations for inflationary cosmology.

Limiting cases

The total density parameter $\Omega_{\rm tot} = \rho/\rho_c$ compared to unity fixes both the sign of the spatial curvature and, for a matter/radiation-only cosmos, the ultimate fate:

Case	Result	Meaning
Closed, positive curvature	$\Omega_{\rm tot} > 1 \iff k = +1$	Density above critical; space is a finite 3-sphere. Without $\Lambda$ the expansion halts and reverses into a Big Crunch.
Flat, zero curvature	$\Omega_{\rm tot} = 1 \iff k = 0$	Density exactly critical; Euclidean geometry, infinite volume. The observed universe sits here to sub-percent precision.
Open, negative curvature	$\Omega_{\rm tot} < 1 \iff k = -1$	Density below critical; hyperbolic geometry, infinite volume. Matter-only case expands forever, $a \propto t$ asymptotically.
Dark-energy limit	$\Omega_\Lambda \to 1$ , $a \to \infty$	$H \to H_0\sqrt{\Omega_\Lambda}$ ; de Sitter exponential expansion, every other $\Omega_i \to 0$ .

Common trapThe critical density is a fixed constant and Ω = 1 means the universe is static.

Both halves are wrong. The critical density is not a constant of nature: $\rho_c = 3H^2/8\pi G$ tracks the square of the expansion rate, so it was vastly larger in the early universe and shrinks as $H$ falls. For the same reason the $\Omega_i$ are functions of time — the familiar figures $\Omega_m \approx 0.31$ , $\Omega_\Lambda \approx 0.69$ are present-day values $\Omega_{i,0}$ , evaluated at $a=1$ , and were completely different in the past (radiation once dominated). And $\Omega_{\rm tot}=1$ says nothing about the universe being static: it means space is flat, Euclidean three-geometry, while that flat space expands perfectly happily.

Worked examples

1 · Weighing the universe in kilograms and atoms

Take $H_0 = 70\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$ ( $h=0.7$ ). Find the present critical density in SI units, express it as protons per cubic metre, and then give the physical mass density of matter and of baryons using $\Omega_m = 0.31$ and $\Omega_b = 0.049$ .

Step 1. Convert the Hubble constant: $H_0 = (7.0\times10^{4}\ \mathrm{m\,s^{-1}})/(3.086\times10^{22}\ \mathrm{m}) = 2.27\times10^{-18}\ \mathrm{s^{-1}}$ . This inverse-second form is what the formula demands.

Step 2. Insert into $\rho_{c,0} = 3H_0^2/8\pi G$ : numerator $3H_0^2 = 1.54\times10^{-35}$ , denominator $8\pi G = 1.68\times10^{-9}$ , giving $\rho_{c,0} = \mathbf{9.2\times10^{-27}\ kg\,m^{-3}}$ , an energy density $\rho_{c,0}c^2 \approx 5\ \mathrm{GeV\,m^{-3}}$ .

Step 3. As a number density, $n = \rho_{c,0}/m_p = (9.2\times10^{-27})/(1.67\times10^{-27}) = \mathbf{5.5\ m^{-3}}$ — between five and six hydrogen atoms per cubic metre.

Step 4. Multiply by the density parameters. Matter: $\rho_m = \Omega_m\,\rho_{c,0} = 0.31\times9.2\times10^{-27} = \mathbf{2.9\times10^{-27}\ kg\,m^{-3}}$ . Baryons: $\rho_b = \Omega_b\,\rho_{c,0} = \mathbf{4.5\times10^{-28}\ kg\,m^{-3}}$ , about 0.27 protons per cubic metre — ordinary atoms are only one part in twenty of the critical density, the rest being dark matter and dark energy.

2 · Geometry and the epoch of dark-energy takeover

A universe is measured to have $\Omega_m = 0.31$ , $\Omega_\Lambda = 0.69$ , $\Omega_r = 9.0\times10^{-5}$ . (a) What is its spatial geometry? (b) At what redshift did matter and radiation have equal density? (c) At what redshift did dark energy overtake matter?

Step 1. Apply the sum rule (5): $\Omega_k = 1 - (\Omega_m + \Omega_\Lambda + \Omega_r) = 1 - (0.31 + 0.69 + 0.00009) = -0.00009$ . To measurement precision $\Omega_k \approx 0$ , so the universe is spatially flat, $k=0$ — this is essentially our own universe.

Step 2. Radiation-matter equality occurs where $\Omega_{r,0}\,a^{-4} = \Omega_{m,0}\,a^{-3}$ , i.e. $a_{\rm eq} = \Omega_{r,0}/\Omega_{m,0} = (9.0\times10^{-5})/0.31 = 2.9\times10^{-4}$ . The redshift is $z_{\rm eq} = 1/a_{\rm eq} - 1 = \mathbf{3400}$ . Before this, radiation ruled; after, matter did.

Step 3. Matter- $\Lambda$ equality occurs where $\Omega_{m,0}\,a^{-3} = \Omega_{\Lambda,0}$ , i.e. $a = (\Omega_{m,0}/\Omega_{\Lambda,0})^{1/3} = (0.31/0.69)^{1/3} = 0.77$ . The redshift is $z = 1/0.77 - 1 = \mathbf{0.3}$ . Dark energy took over only recently — at a lookback time of roughly four billion years, well after the Sun formed.

Step 4. Sanity check with (6): at $a=1$ , $H^2/H_0^2 = 9\times10^{-5} + 0.31 + 0.69 = 1.00$ , as it must, and as $a\to\infty$ the bracket tends to $\Omega_\Lambda$ (the de Sitter limit). Flat geometry, $z_{\rm eq} \approx 3400$ , dark-energy takeover at $z \approx 0.3$ .

Physics of the portrait

The critical density is where Newtonian intuition and relativistic cosmology meet in a single formula. The number $3H^2/8\pi G$ can be guessed from an energy argument about an expanding ball of dust, yet its true meaning — the density that flattens space — belongs to general relativity, and its measured value took a century of instruments to pin down.

Origins · 1 mincritical density and how this was discovered

critical density

the flat-universe density ρ_c = 3H²/8πG and the Ω budget · relativistic cosmology of Friedmann, Einstein, and de Sitter

The idea that a single density separates an eternally expanding universe from one destined to recollapse emerged from Alexander Friedmann's 1922 solutions of Einstein's field equations, which first exhibited expanding and recollapsing cosmologies side by side. Einstein had introduced the cosmological constant $\Lambda$ in 1917 to force a static universe, a move he reportedly regretted once Hubble's 1929 discovery of expansion made a dynamic cosmos unavoidable. In 1932 Einstein and Willem de Sitter distilled the picture into the flat model $\Omega_m = 1$ , $k=0$ — the textbook standard for fifty years — precisely because it sits on the critical borderline, and $\Omega$ became the organizing variable of the effort to "weigh the universe." Two developments overturned the simple $\Omega_m = 1$ expectation: Guth's 1981 inflation, which predicted $\Omega_{\rm tot} = 1$ to extraordinary precision and solved the flatness problem, and the 1998 supernova discovery of accelerating expansion, which revealed a budget filled not by matter but predominantly by dark energy, $\Omega_\Lambda \approx 0.7$ . Planck closed the loop a century after Friedmann: space is flat to a fraction of a percent, and the once-abstract $\rho_c$ is now the precisely known anchor of the entire cosmic inventory.

History

1922

— Friedmann derives the expanding-universe equations and the density that separates their fates

1932

— Einstein and de Sitter propose the flat, critical-density Ω_m = 1 model

1981

— Guth's inflation predicts Ω_total = 1 and resolves the flatness problem

1998

— Type Ia supernovae reveal accelerating expansion and Ω_Λ ≈ 0.7

2018

— Planck measures Ω_k = 0.001 ± 0.002, confirming space is flat

Going further

Texts. Ryden, Introduction to Cosmology, 2nd ed., ch. 5–6 (critical density, density parameters, the flatness problem) · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., ch. 29 · Dodelson & Schmidt, Modern Cosmology, 2nd ed., ch. 1–2 · Mo, van den Bosch & White, Galaxy Formation and Evolution, ch. 3 · Binney & Tremaine, Galactic Dynamics, 2nd ed., cosmology appendix. The numerical anchors used throughout — $\rho_{c,0} = 1.88\times10^{-26}\,h^2\ \mathrm{kg\,m^{-3}} = 2.775\times10^{11}\,h^2\ M_\odot\,\mathrm{Mpc^{-3}}$ , $\Omega_\Lambda \approx 0.69$ , $\Omega_m \approx 0.31$ , $\Omega_b \approx 0.049$ , $\Omega_r \approx 9\times10^{-5}$ , $z_{\rm eq} \approx 3400$ , $\Omega_k = 0.001\pm0.002$ — are Planck 2018 (Results VI: Cosmological Parameters) values.

Prerequisites. The Friedmann Equations · The Hubble Constant · Spatial Curvature of the Universe · The Expansion of the Universe and the Scale Factor · Equation of State in Cosmology. Leads to. The FLRW Metric · The Fluid and Acceleration Equations · Cosmological Redshift.

Ask the tutor: Ask me to walk through why the same $\Omega$ that weighs the universe also measures the shape of space, or to work out any equality redshift or cosmic fate from a given density budget.

Is the critical density the density the universe actually has?

Only in our universe, and only because observations turned out that way. The critical density is defined abstractly as $\rho_c = 3H^2/8\pi G$ — the density that would make space flat at the current expansion rate — for any universe, whatever its real density. The measurements tell us our universe's actual total density equals its critical density to better than a percent, so $\Omega_{\rm tot} \approx 1$ and space is flat. A different universe could have twice or half the critical density and be closed or open. So $\rho_c$ is a yardstick, and the remarkable empirical fact is that our universe measures exactly one yardstick.

Why does curvature get written as a density parameter Ω_k if it is not a substance?

Purely as a bookkeeping convenience, and a very good one. Curvature is not a form of energy — it is a property of geometry — but in the Friedmann equation the term $-kc^2/a^2$ sits algebraically alongside the real densities. Dividing through by $H^2$ turns it into $\Omega_k = -kc^2/(a^2 H^2)$ , which scales as $a^{-2}$ , exactly as a hypothetical fluid with $w = -1/3$ would. Packaging it this way lets you write the clean identity $\sum_i \Omega_i + \Omega_k = 1$ and treat all four terms uniformly. The minus sign is chosen so a closed universe ( $k=+1$ ) has $\Omega_k < 0$ and total density above critical.

How can Ω_m be 0.31 today if radiation once dominated the universe?

Because the density parameters are ratios of two quantities that both change with time, at different rates. The matter density falls as $a^{-3}$ while radiation falls faster, as $a^{-4}$ , because radiation both dilutes and redshifts. Running the clock backward, radiation grows relative to matter, and before $a \approx 3\times10^{-4}$ (redshift about 3400) it dominated — $\Omega_r$ near one, $\Omega_m$ small. The values $\Omega_m \approx 0.31$ , $\Omega_r \approx 9\times10^{-5}$ are strictly today's numbers, the present-day $\Omega_{i,0}$ ; quoting an $\Omega$ without an epoch always means the present.

If the universe is flat, why do we say dark energy and matter each have an Ω less than one?

Because flatness is a statement about the total, not about any single component. Flat means $\Omega_{\rm tot} = \Omega_m + \Omega_r + \Omega_\Lambda = 1$ (with $\Omega_k = 0$ ), and that budget is shared: matter about 0.31, dark energy about 0.69, radiation a whisker. No individual component reaches one, but together they exactly fill the critical density. If any single component had $\Omega > 1$ on its own, the others would have to be negative to keep the total at one, which is unphysical. So "flat" and "each $\Omega_i < 1$ " are perfectly consistent — one talks about the sum, the other about the parts.

Does a closed universe (Ω > 1) have to end in a Big Crunch?

Not anymore, and this is a genuinely modern twist. In the old matter-only picture geometry was destiny: $\Omega > 1$ meant closed and doomed to recollapse. But a cosmological constant breaks that link. Dark energy does not dilute as space expands, so once it dominates it drives ever-accelerating expansion regardless of curvature. A closed universe with enough $\Omega_\Lambda$ expands forever, its matter thinning out while the constant $\Lambda$ term takes over. So the sign of the curvature no longer fixes the fate — the density parameters as a set do. Fate and geometry have been decoupled.

Why does the Hubble tension change the critical density?

Because $\rho_c$ depends on the Hubble constant squared: $\rho_{c,0} = 3H_0^2/8\pi G$ . The competing values — the local distance-ladder $H_0 \approx 73\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$ and the CMB-inferred $H_0 \approx 67$ — differ by about 9%, so their implied critical densities differ by $(73/67)^2 - 1 \approx 19\%$ . Since every density parameter is a physical density divided by $\rho_c$ , a shift in $\rho_{c,0}$ reshuffles the whole inventory: the same measured matter density gives a different $\Omega_m$ depending on which $H_0$ you trust. That is why the tension ripples through the entire cosmic budget rather than staying a peripheral disagreement.