Part I: The Expanding Universe | Chapter 4

Distance Measures in Cosmology

A rigorous derivation of comoving, luminosity, angular diameter, and proper motion distances in an expanding universe with applications to observational cosmology

1. Introduction: Why Multiple Distance Measures?

In Euclidean space, there is only one notion of distance. Given two points, the distance between them is unambiguous. In an expanding universe governed by general relativity, this simplicity evaporates. The fabric of spacetime is dynamic—the metric changes with time—and light signals travel along null geodesics through this evolving geometry. As a result, the “distance” to a distant galaxy depends on how you define the measurement.

The Core Problem

Consider a galaxy at redshift \(z = 2\). The light we observe today was emitted when the universe was roughly 3.3 billion years old. Since that time, the universe has expanded by a factor of \(1 + z = 3\). What do we mean by the “distance” to this galaxy?

The distance the light actually traveled? (Not simply related to any observable.)
The distance inferred from how bright the galaxy appears? (Luminosity distance)
The distance inferred from how large the galaxy appears? (Angular diameter distance)
The “now” distance if we could freeze expansion and measure? (Comoving distance)

Each definition gives a different numerical answer, yet all are physically meaningful and operationally well-defined. The relationships between these distance measures encode the geometry and expansion history of the universe. In a static, flat spacetime they would all agree. Their disagreement is a direct signature of cosmic expansion and spacetime curvature.

We work throughout in the FLRW framework with metric (using conformal notation):

\[ ds^2 = -c^2\,dt^2 + a^2(t)\left[\frac{dr^2}{1 - kr^2} + r^2\,(d\theta^2 + \sin^2\theta\,d\phi^2)\right] \]

where \(a(t)\) is the scale factor (normalized so \(a_0 = 1\) today),\(k = -1, 0, +1\) is the spatial curvature, and \(r\) is the dimensionless comoving radial coordinate.

We also define the Hubble parameter \(H(z)\) via the Friedmann equation:

\[ H(z) = H_0\sqrt{\Omega_r(1+z)^4 + \Omega_m(1+z)^3 + \Omega_k(1+z)^2 + \Omega_\Lambda} \]

where \(\Omega_r, \Omega_m, \Omega_k, \Omega_\Lambda\) are the present-day density parameters for radiation, matter, curvature, and the cosmological constant respectively, satisfying\(\Omega_r + \Omega_m + \Omega_k + \Omega_\Lambda = 1\).

2. Comoving Distance

2.1 Derivation from Null Geodesics

A photon travels along a null geodesic, for which\(ds^2 = 0\). For radial propagation (\(d\theta = d\phi = 0\)), the FLRW metric gives:

\[ 0 = -c^2\,dt^2 + a^2(t)\frac{dr^2}{1 - kr^2} \]

Rearranging and integrating from the emission time \(t_e\) to the observation time\(t_0\) (today):

\[ \int_0^{r}\frac{dr'}{\sqrt{1 - kr'^2}} = c\int_{t_e}^{t_0}\frac{dt'}{a(t')} \]

The right-hand side defines the comoving distance(also called the line-of-sight comoving distance):

\[ \chi \equiv c\int_{t_e}^{t_0}\frac{dt}{a(t)} \]

To convert this to an integral over redshift, use the relation \(1 + z = a_0/a = 1/a\), which gives \(dz = -da/a^2 = -\dot{a}\,dt/a^2\), and therefore\(dt = -dz/[(1+z)H(z)]\). Substituting:

\[ \boxed{\chi(z) = c\int_0^{z}\frac{dz'}{H(z')}} \]

This is the fundamental distance integral of cosmology. All other distances are derived from it.

2.2 Transverse Comoving Distance for Curved Geometries

The left-hand side of our null geodesic equation involves the integral\(\int_0^r dr'/\sqrt{1 - kr'^2}\). This evaluates differently depending on the spatial curvature:

\[ S_k(\chi) = \begin{cases} \frac{1}{\sqrt{|k|}}\sin(\sqrt{|k|}\,\chi) & k > 0 \;\text{(closed)} \\ \chi & k = 0 \;\text{(flat)} \\ \frac{1}{\sqrt{|k|}}\sinh(\sqrt{|k|}\,\chi) & k < 0 \;\text{(open)} \end{cases} \]

In terms of the curvature density parameter, we define \(|k| = |\Omega_k|H_0^2/c^2\), so:

\[ d_M(z) = S_k(\chi) = \begin{cases} \dfrac{c}{H_0\sqrt{|\Omega_k|}}\sin\!\left(H_0\sqrt{|\Omega_k|}\;\dfrac{\chi}{c}\right) & \Omega_k < 0 \\[8pt] \chi & \Omega_k = 0 \\[8pt] \dfrac{c}{H_0\sqrt{\Omega_k}}\sinh\!\left(H_0\sqrt{\Omega_k}\;\dfrac{\chi}{c}\right) & \Omega_k > 0 \end{cases} \]

\(d_M\) is the transverse comoving distance, sometimes denoted \(d_C\) or \(r(z)\). For a flat universe, \(d_M = \chi\).

2.3 Relation to Proper Distance

The proper distance at time \(t\) is the physical distance measured on the spatial hypersurface at that instant:

\[ d_{\text{prop}}(t) = a(t)\,\chi \]

Today (\(a_0 = 1\)), \(d_{\text{prop}} = \chi\). At the time of emission,\(d_{\text{prop}}(t_e) = \chi/(1+z)\). The proper distance changes with timeeven though the comoving distance does not—this is the essence of cosmic expansion. Galaxies are not “moving through space”; rather, the space between them stretches.

3. Luminosity Distance

3.1 Definition and Motivation

In static Euclidean space, a source of luminosity \(L\) at distance \(d\)produces a flux \(F = L/(4\pi d^2)\) (the inverse-square law). The luminosity distance is defined to preserve this relation in an expanding universe:

\[ \boxed{d_L \equiv \sqrt{\frac{L}{4\pi F}}} \]

This is the distance an object would need to be at in a static, Euclidean universe to produce the observed flux. For standard candles (objects of known intrinsic luminosity \(L\)), measuring \(F\) directly gives \(d_L\).

3.2 Derivation for a Flat Universe

Consider photons emitted isotropically by a source at comoving distance \(\chi\). Today, these photons are spread over a sphere of proper radius \(a_0\,\chi = \chi\) (since\(a_0 = 1\)). The surface area is \(4\pi\chi^2\). However, two effects reduce the received flux compared to the static case:

Effect 1: Photon Energy Loss

Each photon is redshifted: \(E_{\text{obs}} = E_{\text{emit}}/(1+z)\). This reduces the energy flux by a factor \(1/(1+z)\).

Effect 2: Time Dilation

Photons arrive at a lower rate: \(\Delta t_{\text{obs}} = (1+z)\,\Delta t_{\text{emit}}\). This further reduces the flux by a factor \(1/(1+z)\).

Combining these effects, the observed flux is:

\[ F = \frac{L}{4\pi\chi^2(1+z)^2} \]

Comparing with \(F = L/(4\pi d_L^2)\), we immediately find:

\[ \boxed{d_L = (1+z)\,\chi} \quad \text{(flat universe)} \]

3.3 Curved Spacetimes

For non-flat geometries, the surface area of a sphere at comoving distance \(\chi\)is not \(4\pi\chi^2\) but \(4\pi S_k(\chi)^2\), where\(S_k\) is the curvature function defined in Section 2.2. Therefore:

\[ F = \frac{L}{4\pi\,S_k(\chi)^2\,(1+z)^2} \]

giving the general result:

\[ \boxed{d_L = (1+z)\,S_k(\chi) = (1+z)\,d_M} \]

3.4 Distance Modulus

Astronomers traditionally work with magnitudes rather than fluxes. The distance modulus is defined as:

\[ \boxed{\mu \equiv m - M = 5\log_{10}\!\left(\frac{d_L}{10\;\text{pc}}\right) = 5\log_{10}(d_L) + 25} \]

where \(m\) is the apparent magnitude, \(M\) is the absolute magnitude, and \(d_L\) is measured in Megaparsecs in the final expression.

For Type Ia supernovae (standard candles with \(M \approx -19.3\)), measuring \(m\)gives \(\mu\), hence \(d_L(z)\). Plotting \(\mu\) vs.\(z\) is the famous Hubble diagram, which revealed the accelerating expansion in 1998.

4. Angular Diameter Distance

4.1 Definition and Derivation

The angular diameter distance is defined by the Euclidean relation between an object's physical (proper) size \(\ell\) and its observed angular size \(\delta\theta\):

\[ \boxed{d_A \equiv \frac{\ell}{\delta\theta}} \]

Consider an object of proper size \(\ell\) at comoving distance \(\chi\), oriented perpendicular to the line of sight. The proper size at the time of emission was related to its comoving size by \(\ell = a(t_e)\,S_k(\chi)\,\delta\theta\). Since\(a(t_e) = 1/(1+z)\):

\[ \delta\theta = \frac{\ell}{a(t_e)\,S_k(\chi)} = \frac{\ell\,(1+z)}{S_k(\chi)} \]

Therefore:

\[ \boxed{d_A = \frac{S_k(\chi)}{1+z} = \frac{d_M}{1+z} = \frac{d_L}{(1+z)^2}} \]

For a flat universe, \(d_A = \chi/(1+z)\).

4.2 The Etherington Reciprocity Theorem

The relation \(d_L = (1+z)^2\,d_A\) is known as the Etherington reciprocity theorem (also called the distance-duality relation). It was proven by I. M. H. Etherington in 1933 and holds forany metric theory of gravity in which photons travel on null geodesics and photon number is conserved. It requires:

Photons propagate along null geodesics (guaranteed by the equivalence principle).
Photon number is conserved (no absorption or emission along the line of sight).
The spacetime is described by a metric theory (Riemannian geometry).

Violations of this relation would signal exotic physics: photon-axion oscillations, varying fundamental constants, or a breakdown of Riemannian geometry. Current observations confirm it to \(\sim 2\%\) accuracy.

4.3 The Angular Diameter Distance Maximum

One of the most counterintuitive results in cosmology: \(d_A\) has a maximum at some redshift \(z_{\max}\)and decreases for higher redshifts. This means distant objects can appear largeron the sky than closer ones!

To see why, note that \(d_A = \chi(z)/(1+z)\) for a flat universe. At low \(z\),\(\chi\) grows roughly linearly and dominates, so \(d_A\) increases. At high \(z\), the \(1/(1+z)\) factor wins because \(\chi\)approaches the particle horizon (a finite limit), while \(1/(1+z) \to 0\). The competition produces a maximum.

For the standard \(\Lambda\)CDM cosmology (\(\Omega_m = 0.315\),\(\Omega_\Lambda = 0.685\)), the maximum occurs at:

\[ z_{\max} \approx 1.6, \qquad d_A^{\max} \approx 1760\;\text{Mpc} \]

This has a profound observational consequence: the cosmic microwave background (CMB), emitted at\(z \approx 1100\), has \(d_A \approx 12.8\;\text{Mpc}\)—much smaller than the maximum. The sound horizon at recombination (\(\approx 147\;\text{Mpc}\)comoving) subtends an angle of about \(1°\) on the sky, which sets the angular scale of the first acoustic peak in the CMB power spectrum.

5. Proper Motion Distance and Comoving Volume

5.1 Proper Motion Distance

The proper motion distance (also called the transverse comoving distance) is defined through the transverse proper velocity. If an object at redshift \(z\) has a transverse proper velocity \(v_\perp\) and an observed angular velocity \(d\theta/dt_0\), then:

\[ d_M \equiv \frac{v_\perp}{d\theta/dt_0} = (1+z)\,d_A = \frac{d_L}{1+z} = S_k(\chi) \]

For a flat universe, \(d_M = \chi\). The proper motion distance is numerically equal to the transverse comoving distance and is the fundamental building block from which all other distances are constructed.

5.2 Parallax Distance

One can also define a parallax distance \(d_P\)based on the trigonometric parallax. For a baseline \(b\) and measured parallax angle\(\pi\), the parallax distance is \(d_P = b/\pi\). In a curved, expanding spacetime, this differs from the other distance measures. For small angles in a flat universe:

\[ d_P = \frac{d_M}{1 + z} = d_A \]

In general, the parallax distance involves the geometry of the full past light cone and can differ from \(d_A\) in curved spacetimes.

5.3 Comoving Volume Element

The comoving volume element is essential for counting objects (galaxies, quasars, supernovae) as a function of redshift. In a solid angle \(d\Omega\) and redshift interval \(dz\):

\[ \boxed{\frac{dV}{dz\,d\Omega} = \frac{d_M^2\,c}{H(z)}} \]

Derivation: The comoving volume of a shell between\(\chi\) and \(\chi + d\chi\) subtending solid angle \(d\Omega\) is:

\[ dV = S_k(\chi)^2\,d\chi\,d\Omega = d_M^2\,d\chi\,d\Omega \]

Since \(d\chi = c\,dz/H(z)\) from the definition of comoving distance:

\[ dV = d_M^2\,\frac{c}{H(z)}\,dz\,d\Omega \]

The total comoving volume out to redshift \(z\) over the full sky (\(4\pi\) sr) is:

\[ V(z) = 4\pi\int_0^z \frac{d_M^2(z')\,c}{H(z')}\,dz' \]

6. Summary of Distance Measures

Distance	Symbol	Formula (flat)	Observable
Comoving	\(\chi\)	\(c\int_0^z dz'/H(z')\)	BAO ruler
Luminosity	\(d_L\)	\((1+z)\,\chi\)	Standard candle flux
Angular diameter	\(d_A\)	\(\chi/(1+z)\)	Standard ruler angle
Proper motion	\(d_M\)	\(\chi = d_L/(1+z)\)	Transverse velocity
Proper (today)	\(d_{\text{prop}}\)	\(\chi\)	—

Key Relationships

\[ d_L = (1+z)\,d_M = (1+z)^2\,d_A \]

All distances agree to first order for \(z \ll 1\): \(d \approx cz/H_0\) (Hubble's law). They diverge at cosmological redshifts, revealing the expansion history.

7. Applications

7.1 The Cosmic Distance Ladder

Measuring distances in cosmology requires a chain of calibrated methods, each valid over a different range:

Nearby (\(z \lesssim 0.01\))

Parallax: Gaia measures stellar parallaxes to \(\sim 10\;\mu\)as accuracy.
Cepheid variables: Period-luminosity relation (Leavitt law) calibrates \(d_L\) to \(\sim 40\) Mpc.
TRGB: Tip of the Red Giant Branch provides an independent calibration.

Cosmological (\(z \gtrsim 0.01\))

Type Ia SNe: Standardizable candles measure \(d_L(z)\) to \(z \sim 2.3\).
BAO: Standard ruler (\(\approx 147\) Mpc) measures \(d_A(z)\) and \(H(z)\).
CMB: Sound horizon angular size gives \(d_A(z \approx 1100)\).

7.2 The Hubble Tension

One of the most pressing problems in modern cosmology is the disagreement between two methods of measuring \(H_0\):

Early Universe (Planck CMB)

\(H_0 = 67.4 \pm 0.5\;\text{km/s/Mpc}\)

Derived from fitting the \(\Lambda\)CDM model to CMB anisotropies at\(z \approx 1100\) and extrapolating to today.

Late Universe (SH0ES Cepheids + SNe)

\(H_0 = 73.0 \pm 1.0\;\text{km/s/Mpc}\)

From the local distance ladder: Cepheids calibrate Type Ia supernovae, giving\(d_L(z)\) directly at low \(z\).

The tension stands at \(\sim 5\sigma\). Possible explanations include early dark energy, modified gravity, new neutrino physics, or systematic errors in the distance ladder. Resolving this tension is a central goal of DESI, Euclid, and the Vera Rubin Observatory.

8. Historical Context

Key Figures

Edwin Hubble (1929): Discovered the linear velocity-distance relation \(v = H_0 d\) using Cepheid-calibrated distances to “nebulae” (galaxies), establishing that the universe is expanding. His original value of \(H_0 \approx 500\;\text{km/s/Mpc}\) was far too large due to calibration errors in the distance ladder.

Allan Sandage (1958–2000s): Hubble's scientific heir at Mt. Wilson/Palomar. Spent decades refining\(H_0\), championing values near \(50\;\text{km/s/Mpc}\) against Gérard de Vaucouleurs' higher estimates (\(\sim 100\)). This “Hubble wars” lasted from the 1970s through the 1990s.

Adam Riess & the High-z Supernova Search Team (1998): Discovered that Type Ia supernovae at \(z \sim 0.5\) are fainter than expected in a decelerating universe, implying \(d_L\) is larger than predicted—the universe's expansion is accelerating.

Saul Perlmutter & the Supernova Cosmology Project (1999): Independently confirmed the accelerating expansion using a different sample of Type Ia SNe. The \(d_L(z)\) data from both teams required\(\Omega_\Lambda \sim 0.7\), resurrecting Einstein's cosmological constant.

Brian Schmidt (1998): Led the High-z team alongside Riess. Together with Perlmutter, the three shared the 2011 Nobel Prize in Physics for the discovery of the accelerating expansion of the universe through observations of distant supernovae.

The measurement of cosmological distances has been one of the great intellectual achievements of the 20th century. From Hubble's first estimates using Cepheids to the precision measurements by Planck and DESI, the distance measures derived in this chapter underpin our entire quantitative understanding of the cosmos.

9. Python Simulation: Distance Measures in \(\Lambda\)CDM

The following simulation computes and plots all four distance measures as a function of redshift for the standard \(\Lambda\)CDM cosmology, along with the Hubble diagram (distance modulus vs. redshift). We use the trapezoidal rule for numerical integration—no external dependencies beyond NumPy.

Python

script.py129 lines

import numpy as np

# ============================================================
# Cosmological Parameters (Planck 2018)
# ============================================================
H0_km = 67.4          # Hubble constant [km/s/Mpc]
H0 = H0_km * 1e3 / 3.0857e22  # Convert to [1/s]
c = 2.998e5           # Speed of light [km/s]
Omega_r = 9.1e-5      # Radiation density parameter
Omega_m = 0.315       # Matter density parameter
Omega_L = 0.685       # Dark energy (cosmological constant)
Omega_k = 1.0 - Omega_r - Omega_m - Omega_L  # Curvature (~ 0)

# ============================================================
# Hubble parameter H(z) / H0
# ============================================================
def E(z):
    """Dimensionless Hubble parameter E(z) = H(z)/H0"""
    return np.sqrt(Omega_r*(1+z)**4 + Omega_m*(1+z)**3
                   + Omega_k*(1+z)**2 + Omega_L)

# ============================================================
# Numerical integration via trapezoidal rule
# ============================================================
def trapz_cumulative(y, x):
    """Cumulative trapezoidal integration."""
    result = np.zeros(len(x))
    for i in range(1, len(x)):
        result[i] = result[i-1] + 0.5*(y[i-1] + y[i])*(x[i] - x[i-1])
    return result

# ============================================================
# Compute distance measures
# ============================================================
N = 2000
z_max = 5.0
z = np.linspace(0, z_max, N)
dH = c / H0_km   # Hubble distance [Mpc]

# Comoving distance: chi = (c/H0) * integral_0^z dz'/E(z')
integrand = 1.0 / E(z)
chi = dH * trapz_cumulative(integrand, z)  # [Mpc]

# For flat universe (Omega_k ~ 0): d_M = chi
d_M = chi.copy()

# Luminosity distance
d_L = (1 + z) * d_M

# Angular diameter distance
d_A = d_M / (1 + z)

# Proper distance (today) = chi (same as comoving for a_0 = 1)
d_prop = chi.copy()

# Distance modulus (only for z > 0)
mu = np.full_like(z, np.nan)
mask = z > 0.001
mu[mask] = 5.0 * np.log10(d_L[mask]) + 25.0  # d_L in Mpc

# ============================================================
# Print key values
# ============================================================
print("=" * 65)
print("DISTANCE MEASURES IN LCDM COSMOLOGY")
print("Planck 2018: H0=67.4 km/s/Mpc, Om=0.315, OL=0.685")
print("=" * 65)
print(f"{'z':>6}  {'chi [Mpc]':>12}  {'d_L [Mpc]':>12}  {'d_A [Mpc]':>12}  {'mu [mag]':>10}")
print("-" * 65)

z_values = [0.1, 0.5, 1.0, 1.5, 2.0, 3.0, 5.0]
for zv in z_values:
    idx = np.argmin(np.abs(z - zv))
    chi_v = chi[idx]
    dL_v = d_L[idx]
    dA_v = d_A[idx]
    mu_v = mu[idx] if not np.isnan(mu[idx]) else 0.0
    print(f"{zv:6.1f}  {chi_v:12.1f}  {dL_v:12.1f}  {dA_v:12.1f}  {mu_v:10.2f}")

# Find d_A maximum
idx_max = np.argmax(d_A[1:]) + 1
print(f"\nAngular diameter distance maximum:")
print(f"  z_max = {z[idx_max]:.2f},  d_A_max = {d_A[idx_max]:.1f} Mpc")

print("\n" + "=" * 65)
print("TABULAR OUTPUT FOR PLOTTING")
print("=" * 65)
print(f"{'z':>8}  {'chi':>12}  {'d_L':>12}  {'d_A':>12}  {'d_prop':>12}")

# Output every 50th point for the chart
for i in range(0, N, 50):
    print(f"{z[i]:8.3f}  {chi[i]:12.2f}  {d_L[i]:12.2f}  {d_A[i]:12.2f}  {d_prop[i]:12.2f}")
# Always include last point
print(f"{z[-1]:8.3f}  {chi[-1]:12.2f}  {d_L[-1]:12.2f}  {d_A[-1]:12.2f}  {d_prop[-1]:12.2f}")

print("\n" + "=" * 65)
print("HUBBLE DIAGRAM DATA (distance modulus vs z)")
print("=" * 65)
print(f"{'z':>8}  {'mu':>10}")
for i in range(20, N, 50):
    if not np.isnan(mu[i]):
        print(f"{z[i]:8.3f}  {mu[i]:10.2f}")
print(f"{z[-1]:8.3f}  {mu[-1]:10.2f}")

# ============================================================
# Comparison: Empty universe (Milne model) vs LCDM
# ============================================================
print("\n" + "=" * 65)
print("COMPARISON: LCDM vs EMPTY (MILNE) UNIVERSE")
print("=" * 65)
def E_milne(z):
    return (1+z)  # Omega_k = 1, all others = 0

integrand_milne = 1.0 / E_milne(z)
chi_milne = dH * trapz_cumulative(integrand_milne, z)
d_L_milne = (1+z) * chi_milne
mu_milne = np.full_like(z, np.nan)
mask2 = z > 0.001
mu_milne[mask2] = 5.0*np.log10(d_L_milne[mask2]) + 25.0

print(f"{'z':>8}  {'mu_LCDM':>10}  {'mu_Milne':>10}  {'Delta_mu':>10}")
for zv in [0.1, 0.3, 0.5, 0.7, 1.0, 1.5, 2.0]:
    idx = np.argmin(np.abs(z - zv))
    dm = mu[idx] - mu_milne[idx] if not np.isnan(mu[idx]) else 0
    print(f"{zv:8.1f}  {mu[idx]:10.2f}  {mu_milne[idx]:10.2f}  {dm:10.3f}")

print("\nPositive Delta_mu means SNe are FAINTER than expected in an")
print("empty universe => the expansion is ACCELERATING (dark energy!)")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Reading the Output

All distances agree at low z: For \(z \ll 1\), all distance measures converge to \(d \approx cz/H_0\) (the Hubble law regime).
\(d_L\) is always the largest: The luminosity distance grows as \((1+z)\chi\), making distant objects appear fainter than their comoving distance would suggest.
\(d_A\) turns over: The angular diameter distance reaches a maximum near \(z \approx 1.6\) and decreases at higher redshifts.
Distance modulus residuals: The difference \(\Delta\mu\) between\(\Lambda\)CDM and an empty (Milne) universe is positive at \(z \gtrsim 0.3\), meaning supernovae are fainter than expected—this is the signature of accelerating expansion discovered in 1998.

References

Hogg, D. W. (2000). “Distance Measures in Cosmology,” arXiv:astro-ph/9905116. — The classic pedagogical reference for cosmological distance measures, which this chapter closely follows.
Weinberg, S. (1972). Gravitation and Cosmology, Wiley. — The standard graduate text deriving all distance measures from first principles.
Etherington, I. M. H. (1933). “On the Definition of Distance in General Relativity,” Phil. Mag. 15, 761–773. — The original proof of the reciprocity theorem \(d_L = (1+z)^2 d_A\).
Riess, A. G. et al. (1998). “Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant,” Astron. J. 116, 1009. — The High-z team's discovery of cosmic acceleration via \(d_L(z)\) measurements.
Perlmutter, S. et al. (1999). “Measurements of \(\Omega\) and \(\Lambda\) from 42 High-Redshift Supernovae,” Astrophys. J. 517, 565. — The Supernova Cosmology Project's independent confirmation.
Planck Collaboration (2020). “Planck 2018 results. VI. Cosmological parameters,” Astron. Astrophys. 641, A6. — Definitive measurement of cosmological parameters used throughout this chapter.
Riess, A. G. et al. (2022). “A Comprehensive Measurement of the Local Value of the Hubble Constant,” Astrophys. J. Lett. 934, L7. — The SH0ES team's latest \(H_0 = 73.0 \pm 1.0\) km/s/Mpc measurement.

← Cosmological Redshift Cosmological Horizons →