Spacetime Diagram

The Davis & Lineweaver (2004) view of an expanding universe. Past and future light cones, particle horizon, event horizon, the Hubble sphere, and galaxy worldlines — all on one plot, for any wCDM cosmology you like.

The Spacetime Diagram — Comoving distance vs Cosmic time After Davis & Lineweaver (2004)
01 — Cosmological Model Live updates — diagram redraws as you drag
Preset
Hubble Constant H₀ km/s/Mpc
Matter Density ΩM
Vacuum Energy ΩΛ
Curvature ΩK flat
auto: 1 − Ωᵐ − ΩΛ
Dark Energy EoS w
02 — Display Controls Toggle layers on / off
Emission redshift z₂ amber dot on diagram
Reading the Diagram What every line is telling you

Particle horizon

The comoving distance light has had time to travel since the Big Bang. At each cosmic time on the y-axis, anything outside this curve is causally disconnected from us — light from it hasn’t reached us yet.

Event horizon

The comoving distance from which light emitted now could ever reach us in the future. In an accelerating universe this is finite (about 16–17 Glyr for Planck 2018) — galaxies beyond it will redshift out of contact forever.

Hubble sphere

The surface where the recession velocity vrec = H(t) D equals c. In comoving coordinates the Hubble sphere shrinks in the de Sitter era — one of the more counter-intuitive features Davis & Lineweaver use this diagram to clarify.

Light cones

The past light cone (solid) traces the worldline of every photon arriving at us right now. The future light cone (dotted) traces every photon we emit today — it asymptotes to the event horizon, never crossing it.

Worldlines

Galaxies sit at constant comoving distance, so their worldlines are vertical lines in this view. We mark the worldlines for objects we observe today at z = 1, 3, 10, and the cosmic microwave background last-scattering surface at z ≈ 1100.

“Now” line and us

The horizontal dashed line at t = t₀ is the present epoch. The white star at (0, t₀) is our worldline crossing “now”. Following the past light cone down from that point traces back every direction in the sky.

Method & References

Everything is computed from the wCDM expansion function E(z) = √[ΩM(1+z)3 + ΩK(1+z)2 + ΩΛ(1+z)3(1+w)]. The diagram requires four integrals, each evaluated with Simpson’s rule:

Particle horizon today D_p(0) = (c/H₀) · ∫₀ dz/E(z) Event horizon today D_e(0) = (c/H₀) · ∫-10 dz/E(z) Comoving D_C(z) from us = (c/H₀) · ∫₀z dz/E(z) Cosmic time t(a) = (1/H₀) · ∫₀a da′/[a′ E(a′)]

For the high-z particle-horizon integral the calculator substitutes u = ln(1+z) so Simpson’s rule samples both the low-z bulk and the high-z tail evenly — without that change, the linear-z rule undersamples z < 5 and the particle horizon comes out roughly 3× too large. Cross-check (Planck 2018): Dp(0) ≈ 46.8 Glyr, De(0) ≈ 16.7 Glyr, DH(0) = c/H₀ ≈ 14.5 Glyr.

  • T. M. Davis & C. H. Lineweaver, “Expanding Confusion: Common Misconceptions of Cosmological Horizons and the Superluminal Expansion of the Universe,” PASA 21, 97 (2004). arXiv:astro-ph/0310808
  • D. W. Hogg, “Distance measures in cosmology,” (1999). Used for the transverse-distance / sinh-sin curvature handling. arXiv:astro-ph/9905116
  • E. F. Bunn & D. W. Hogg, “The kinematic origin of the cosmological redshift,” AJP 77, 688 (2009). For follow-up on the kinematic / GR interpretation of the recession-velocity surface. arXiv:0808.1081
  • See also the Cosmology Calculator for the full set of distance / time / observational outputs (DC, DA, DL, μ, batch processing, reverse lookups, etc.).