Free coursewareRewritten from MIT OpenCourseWareCC BY-NC-SA 4.0
Open CosmologyPhysics & cosmology: four years, plus a graduate year

How bright is a distant supernova?

The formula the 1998 supernova teams fitted to their data — the measured intensity of a standard candle as a function of its redshift, for a universe with matter, radiation, vacuum energy and curvature.

You need this first

Full reading path (32 pages)
  1. Vectors (8.01)
  2. Position and velocity in one dimension (8.01)
  3. Acceleration in one dimension (8.01)
  4. Motion in two dimensions and projectiles (8.01)
  5. The homogeneous expanding universe (8.286)
  6. The Doppler effect and special relativity (8.286)
  7. Cosmological redshift (8.286)
  8. Newton's laws of motion (8.01)
  9. The universal law of gravitation (8.01)
  10. Why a static universe cannot work (8.286)
  11. The expanding sphere and the equation for a(t) (8.286)
  12. The Friedmann equation (8.286)
  13. Open, closed and flat universes (8.286)
  14. The closed universe as a cycloid (8.286)
  15. The age of the universe and the Big Bang singularity (8.286)
  16. The age of a closed universe from H and Ω (8.286)
  17. The open universe and the age of any matter universe (8.286)
  18. Non-Euclidean geometry (8.286)
  19. The metric of a sphere and of a closed universe (8.286)
  20. Curvature and expansion — the Robertson–Walker metric (8.286)
  21. The geometry of open universes (8.286)
  22. The spacetime interval (8.286)
  23. The spacetime metric of the universe (8.286)
  24. The geodesic equation (8.286)
  25. Simplifying the geodesic equation (8.286)
  26. Falling into a black hole (8.286)
  27. Mass, energy and the dilution of radiation (8.286)
  28. Matter–radiation equality and the pressure of light (8.286)
  29. The horizon distance (8.286)
  30. The radiation-dominated universe (8.286)
  31. The cosmological constant and vacuum energy (8.286)
  32. The age of the universe with vacuum energy (8.286)
On this page · 5 sections
  1. The metric in terms of an angle
  2. Spreading the light over a sphere
  3. Rewriting it in terms astronomers measure
  4. Flat and open universes
  5. Summary

The second calculation of Lecture 19 is the one behind the discovery of acceleration itself. In 1998 two teams used type Ia supernovae as standard candles to trace the expansion rate back in time. Looking back about 5–6 billion years, they found the expansion was slower then than now.

What they actually measured was simply each supernova's apparent brightness against its redshift: for a standard candle, brightness stands in for distance. So the question for us is: as a function of the cosmological parameters, what intensity should we expect for a source of known power at redshift zz? Guth stresses that the formula we are about to derive is exactly the one the teams fitted to their data.

Nothing here is really new — the class has calculated intensities before — but now the equations appear in full, with vacuum energy, and allowing for curvature. In the end, Guth says, the curvature of our universe turns out to be negligible as far as anyone can tell: it is not there at the level of one part in 100, though it could still be there at about one part in 1,000, and people keep looking.

The metric in terms of an angle#

Start from the Robertson–Walker metric:

Robertson–Walker metric
ds2=−c2dt2+a2(t){dr21−kr2+r2(dθ2+sin⁡2θ dϕ2)}ds^2 = -c^2dt^2 + a^2(t)\left\{\frac{dr^2}{1 - kr^2} + r^2\left(d\theta^2 + \sin^2\theta\,d\phi^2\right)\right\}
In words

ds squared equals minus c squared dt squared plus a squared (t) dr squared over 1 minus kr squared plus r squared (d theta squared plus sine squared theta d phi squared )

Source: 8.286 Lecture 19 transcript, pp. 13–14 (board)

Light from the source to us travels radially, so it pays to simplify the radial part. For a flat universe nothing needs doing; for open and closed universes the best variable differs. Guth takes the closed case (k>0k > 0) and introduces an angle ψ\psi — which, traced back, is the angle from the ww-axis used when the closed universe was first built as a 3-sphere in four dimensions:

Radial angle
sin⁡ψ=k r⟹dψ=k drcos⁡ψ=k dr1−kr2\sin\psi = \sqrt k\,r \quad\Longrightarrow\quad d\psi = \frac{\sqrt k\,dr}{\cos\psi} = \frac{\sqrt k\,dr}{\sqrt{1 - kr^2}}
In words

sine psi equals the square root of k r which gives d psi equals the square root of k dr over cos psi equals the square root of k dr over the square root of 1 minus kr squared

Source: 8.286 Lecture 19 transcript, p. 14 (board)

The radial term of the metric is exactly the square of this, so

Closed metric with angle ψ
ds2=−c2dt2+a~2(t){dψ2+sin⁡2ψ(dθ2+sin⁡2θ dϕ2)},a~(t)=a(t)kds^2 = -c^2dt^2 + \tilde a^2(t)\left\{d\psi^2 + \sin^2\psi\left(d\theta^2 + \sin^2\theta\,d\phi^2\right)\right\}, \qquad \tilde a(t) = \frac{a(t)}{\sqrt k}
In words

ds squared equals minus c squared dt squared plus a tilde squared (t) d psi squared plus sine squared psi (d theta squared plus sine squared theta d phi squared ) , a tilde (t) equals a(t) over the square root of k

Source: 8.286 Lecture 19 transcript, p. 14 (board)

The tilde only distinguishes the new scale factor from the old one; starting from this metric, you could drop it. A student asks whether aa and a~\tilde a have different units. They do. In Guth's convention rr is a coordinate distance in notches, so aa is in metres per notch; ψ\psi is an angle and naturally dimensionless, so a~\tilde a is simply a length.

Spreading the light over a sphere#

Put the source at the centre (ψ=0\psi = 0), and imagine a sphere around it, drawn in comoving coordinates so that our detector, of area AA, lies on its surface at ψ=ψd\psi = \psi_d. The fraction of the light that hits the detector is the detector's area divided by the sphere's area — and the sphere's area must be computed from the metric.

On the sphere only θ\theta and ϕ\phi vary. The angular part of the metric is the usual one for a sphere, multiplied by the square of the sphere's radius, which is therefore a~sin⁡ψ\tilde a\sin\psi. We evaluate it today, t0t_0, because the detector has area AA today and the comparison only makes sense if all distances are measured at the same time. So the fraction is

Fraction of the light reaching the detector
A4π a~2(t0) sin⁡2ψd\frac{A}{4\pi\,\tilde a^2(t_0)\,\sin^2\psi_d}
In words

A over 4 pi a tilde squared (t sub 0 ) sine squared psi sub d

Source: 8.286 Lecture 19 transcript, pp. 15–16 (board)

The redshift costs two more factors. Any clock at the source appears slowed by 1+zs1 + z_s, where zsz_s is the redshift of the source. One clock is the photon's frequency, so each photon arrives with its energy reduced by 1+zs1 + z_s; the other is the rate at which photons arrive, reduced by the same factor:

Redshift
1+zs=a(t0)a(te)1 + z_s = \frac{a(t_0)}{a(t_e)}
In words

1 plus z sub s equals a(t sub 0 ) over a(t sub e )

Source: 8.286 Lecture 19 transcript, p. 16 (board)

With PP the power emitted by the source, the power received and the measured intensity JJ (power per unit area of detector) are

Received power and intensity
Prec=P(1+zs)2 A4π a~2(t0)sin⁡2ψd,J=PrecA=P4π(1+zs)2 a~2(t0) sin⁡2ψdP_\text{rec} = \frac{P}{(1+z_s)^2}\,\frac{A}{4\pi\,\tilde a^2(t_0)\sin^2\psi_d}, \qquad J = \frac{P_\text{rec}}{A} = \frac{P}{4\pi(1+z_s)^2\,\tilde a^2(t_0)\,\sin^2\psi_d}
In words

P sub rec equals P over (1 plus z sub s ) squared A over 4 pi a tilde squared (t sub 0 ) sine squared psi sub d , J equals P sub rec over A equals P over 4 pi (1 plus z sub s ) squared a tilde squared (t sub 0 ) sine squared psi sub d

Source: 8.286 Lecture 19 transcript, pp. 16–17 (board)

Rewriting it in terms astronomers measure#

That is effectively the answer, but a~(t0)\tilde a(t_0) and ψd\psi_d mean nothing to an astronomer. Both can be traded for measured quantities.

The scale factor. From the definition of the curvature fraction, Ωk,0=−kc2/(a2(t0)H02)\Omega_{k,0} = -kc^2/(a^2(t_0)H_0^2):

ã(t₀) from the curvature fraction
a~(t0)=c H0−1−Ωk,0\tilde a(t_0) = \frac{c\,H_0^{-1}}{\sqrt{-\Omega_{k,0}}}
In words

a tilde (t sub 0 ) equals c H sub 0 to the power minus 1 over the square root of minus omega sub k,0

Source: 8.286 Lecture 19 transcript, p. 17 (board)

For a closed universe k>0k > 0, so Ωk,0\Omega_{k,0} is negative and the square root is real. [The transcript reads as if the square root multiplies; from the definition of Ωk,0\Omega_{k,0} it divides.]

The angle. Light travels on null geodesics. For a radial ray, ds2=0ds^2 = 0:

Radial light ray
−c2dt2+a~2(t) dψ2=0⟹dψdt=ca~(t)-c^2dt^2 + \tilde a^2(t)\,d\psi^2 = 0 \quad\Longrightarrow\quad \frac{d\psi}{dt} = \frac{c}{\tilde a(t)}
In words

minus c squared dt squared plus a tilde squared (t) d psi squared equals 0 which gives d psi over dt equals c over a tilde (t)

Source: 8.286 Lecture 19 transcript, p. 18 (board)

Earlier in the course this rule was motivated by intuition; here it follows directly from the metric, even in a curved universe. Integrating from emission to today gives ψd=∫tet0c dt/a~(t)\psi_d = \int_{t_e}^{t_0} c\,dt/\tilde a(t).

To integrate over redshift instead of time — redshift is how astronomers label the past — use 1+z=1/x=a~(t0)/a~(t)1 + z = 1/x = \tilde a(t_0)/\tilde a(t) and differentiate:

Changing variable to redshift
dz=−a~(t0)a~2(t) a~˙ dt=−a~(t0) H(t) dta~(t)⟹ψs=1a~(t0)∫0zsc dzH(z)dz = -\frac{\tilde a(t_0)}{\tilde a^2(t)}\,\dot{\tilde a}\,dt = -\tilde a(t_0)\,H(t)\,\frac{dt}{\tilde a(t)} \quad\Longrightarrow\quad \psi_s = \frac{1}{\tilde a(t_0)}\int_0^{z_s}\frac{c\,dz}{H(z)}
In words

dz equals minus a tilde (t sub 0 ) over a tilde squared (t) a tilde dot dt equals minus a tilde (t sub 0 ) H(t) dt over a tilde (t) which gives psi sub s equals 1 over a tilde (t sub 0 ) the integral from 0 to z sub s of c dz over H(z)

Source: 8.286 Lecture 19 transcript, p. 18 (board)

The Friedmann equation from the first half of the lecture, H2=H02F(x)/x4H^2 = H_0^2F(x)/x^4 with x=1/(1+z)x = 1/(1+z), supplies H(z)H(z). Substituting it and the expression for a~(t0)\tilde a(t_0):

Angle to a source at redshift z_s
ψ(zs)=∣Ωk,0∣∫0zsdzΩm,0(1+z)3+Ωrad,0(1+z)4+Ωvac,0+Ωk,0(1+z)2\psi(z_s) = \sqrt{|\Omega_{k,0}|}\int_0^{z_s}\frac{dz}{\sqrt{\Omega_{m,0}(1+z)^3 + \Omega_{\text{rad},0}(1+z)^4 + \Omega_{\text{vac},0} + \Omega_{k,0}(1+z)^2}}
In words

psi (z sub s ) equals the square root of | omega sub k,0 | the integral from 0 to z sub s of dz over the square root of omega sub m,0 (1 plus z) cubed plus omega sub rad ,0 (1 plus z) to the power 4 plus omega sub vac ,0 plus omega sub k,0 (1 plus z) squared

Source: 8.286 Lecture 19 transcript, p. 19 (board)

Put this ψ\psi and a~(t0)\tilde a(t_0) into the formula for JJ, and the expected intensity depends only on the present Ω's, H0H_0, the redshift of the source and its intrinsic power. That, Guth says, is exactly the formula the supernova teams used in 1998 — and to fit their data they needed a very significant vacuum energy, of order 60 or 70 percent.

Flat and open universes#

Guth works out only the closed case. The other two follow from it (this section is the page's own addition):

  • Open (k<0k < 0, Ωk,0>0\Omega_{k,0} > 0): the same steps with sinh⁡ψ=∣k∣ r\sinh\psi = \sqrt{|k|}\,r give the same formula with sin⁡ψd\sin\psi_d replaced by sinh⁡ψd\sinh\psi_d, and −Ωk,0\sqrt{-\Omega_{k,0}} by Ωk,0\sqrt{\Omega_{k,0}}.
  • Flat (Ωk,0→0\Omega_{k,0} \to 0): ψd\psi_d becomes tiny and sin⁡ψd≈ψd\sin\psi_d \approx \psi_d, so a~(t0)sin⁡ψd→(c/H0)∫0zsdz/E(z)\tilde a(t_0)\sin\psi_d \to (c/H_0)\int_0^{z_s}dz/E(z), where E(z)=H(z)/H0E(z) = H(z)/H_0 is the square root in the integral above.

Astronomers package all three cases into a luminosity distance dLd_L, defined so that J=P/(4πdL2)J = P/(4\pi d_L^2) looks like the inverse-square law of a static, flat space. Comparing with the formula above, dL=(1+zs) a~(t0)sin⁡ψdd_L = (1 + z_s)\,\tilde a(t_0)\sin\psi_d.

Computed from the formulas: the intensity of a standard candle in a flat universe with matter and vacuum energy (Ω_vac = 1 − Ω_m, radiation neglected), divided by the intensity the same candle would have in a flat matter-only universe with the same H₀. The matter-only case has the closed form H₀d_L/c = 2(1+z)(1 − 1/√(1+z)); the other is the integral above, evaluated numerically. For Ω_m = 0.3 the candle looks 69% as bright at z = 0.5 and 58% at z = 1: the supernovae were fainter, so farther away, than a decelerating universe allows.
Worked example

A supernova at z = 0.5

A type Ia supernova is seen at z=0.5z = 0.5. Take H0=70H_0 = 70 km/s/Mpc and compare two flat universes: (A) matter only, (B) Ωm=0.3\Omega_m = 0.3, Ωvac=0.7\Omega_\text{vac} = 0.7. Find the luminosity distance in each, and how much fainter the supernova looks in (B). Then check how much a slight closed curvature, Ωm,0=0.31\Omega_{m,0} = 0.31, Ωvac,0=0.70\Omega_{\text{vac},0} = 0.70 (so Ωk,0=−0.01\Omega_{k,0} = -0.01), changes model (B).

(A) Matter only

E(z)=(1+z)3/2E(z) = (1+z)^{3/2}, so ∫00.5dz/E=2(1−1/1.5)=0.3670\int_0^{0.5}dz/E = 2(1 - 1/\sqrt{1.5}) = 0.3670. Then H0dL/c=1.5×0.3670=0.5505H_0d_L/c = 1.5 \times 0.3670 = 0.5505. With c/H0=4283c/H_0 = 4283 Mpc, dL=2358d_L = 2358 Mpc.

(B) Matter + vacuum energy

E(z)=0.3(1+z)3+0.7E(z) = \sqrt{0.3(1+z)^3 + 0.7}. Numerically, ∫00.5dz/E=0.4410\int_0^{0.5}dz/E = 0.4410, so H0dL/c=0.6615H_0d_L/c = 0.6615 and dL=2833d_L = 2833 Mpc.

How much fainter

J∝1/dL2J \propto 1/d_L^2, so JB/JA=(0.5505/0.6615)2=0.69J_B/J_A = (0.5505/0.6615)^2 = 0.69. The supernova delivers only 69% of the light expected in a decelerating matter-only universe.

A little curvature

Now ψd=0.01∫00.5dz/E=0.0440\psi_d = \sqrt{0.01}\int_0^{0.5}dz/E = 0.0440 radians — a small angle, as the closed universe is huge compared with the distance to the source. a~(t0)=(c/H0)/0.1\tilde a(t_0) = (c/H_0)/0.1, so H0dL/c=1.5×sin⁡(0.0440)/0.1=0.6603H_0d_L/c = 1.5 \times \sin(0.0440)/0.1 = 0.6603, only 0.2% below the flat value. Curvature at the 1% level hardly changes the brightness at this redshift; the vacuum energy is what matters.

Summary#

  • The 1998 teams measured the intensity of type Ia supernovae against redshift; for standard candles, intensity measures distance.
  • In a closed universe, with sin⁡ψ=k r\sin\psi = \sqrt k\,r and a~=a/k\tilde a = a/\sqrt k, the light spreads over a sphere of area 4πa~2(t0)sin⁡2ψd4\pi\tilde a^2(t_0)\sin^2\psi_d, and redshift costs two factors of 1+zs1+z_s: J=P/[4π(1+zs)2a~2(t0)sin⁡2ψd]J = P/[4\pi(1+z_s)^2\tilde a^2(t_0)\sin^2\psi_d].
  • a~(t0)=cH0−1/−Ωk,0\tilde a(t_0) = cH_0^{-1}/\sqrt{-\Omega_{k,0}}, and the null-geodesic condition dψ/dt=c/a~d\psi/dt = c/\tilde a gives ψs=1a~(t0)∫0zsc dz/H(z)\psi_s = \frac{1}{\tilde a(t_0)}\int_0^{z_s}c\,dz/H(z), with H(z)H(z) from the Ω's.
  • This is the formula fitted to the supernova data; it required a vacuum energy of order 60–70% of the critical density.

Symbols on this page

Show the 28 symbols used here
θ\theta theta
Polar angle: measured from the north pole (0) to the south pole (π). Say it: “theta” (sounds like “thay-ta or thee-ta”).
π\pi pi
The ratio of a circle's circumference to its diameter, π ≈ 3.14159. Say it: “pi”.
ϕ\phi phi
Azimuthal angle (longitude), running from 0 to 2π around the axis. Say it: “phi” (sounds like “fye”).
ψ\psi psi
The extra angle of the 3-sphere: how far round the closed universe a point is. Say it: “psi” (sounds like “sigh”).
Ω\Omega capital omega
The density parameter Ω = ρ/ρ_c. Subscripts pick one ingredient: Ω_m matter, Ω_vac (or Ω_Λ) vacuum energy, Ω_r radiation, Ω_k curvature. Say it: “omega”.
≈\approx approximately equal
Is approximately equal to. Say it: “is approximately”.
∝\propto proportional to
Is proportional to: doubling one doubles the other. Say it: “is proportional to”.
→\to arrow
Goes to (tends to), or becomes. Say it: “goes to”.
⇒\Rightarrow implies
Implies: what follows is a consequence of what came before. Say it: “which gives”.
×\times times
Times (multiplication). Between two vectors: the cross product. Say it: “times”.
∫\int integral
Integral: the sum of infinitely many small pieces, such as the area under a curve. Say it: “the integral of”.
sin⁡\sin sine
Sine of an angle. Say it: “sine”.
cos⁡\cos cosine
Cosine of an angle. Say it: “cos”.
sinh⁡\sinh hyperbolic sine
Hyperbolic sine, (eˣ − e⁻ˣ)/2. Say it: “sinh” (sounds like “shine or sinch”).
a˙\dot a dot
Rate of change with time: ȧ = da/dt. Say it: “dot”.
a~\tilde a tilde
A modified version of a quantity (ã, k̃); its meaning is given where it is introduced. Say it: “tilde”.
aa a
The scale factor a(t): how much the universe has expanded. Physical distance = a × comoving distance (Guth's units: metres per notch). Say it: “a”.
HH H
The Hubble parameter H = ȧ/a: the expansion rate of the universe at a given time. Say it: “H”.
zz z
Redshift: 1 + z = a(observed)/a(emitted), how much the universe has stretched since the light left. Say it: “z”.
cc c
The speed of light, c ≈ 2.998 × 10⁸ m/s. Say it: “c”.
kk k
The curvature constant: k > 0 closed, k = 0 flat, k < 0 open (Guth's units: 1/notch²). Say it: “k”.
H0H_{0} Hubble constant
The Hubble constant: today's value of H, about 67–73 km/s/Mpc. Say it: “H nought”.
t0t_{0} today
Today: the present time, about 13.8 billion years after the Big Bang in the real universe (simple models on some pages give other ages). Say it: “t nought”.
Ωm\Omega_{m} matter density parameter
The matter density as a fraction of the critical density, ρ_m/ρ_c; ',0' means today (about 0.3). Say it: “omega m”.
Ωvac\Omega_{\text{vac}} vacuum density parameter
The vacuum-energy density as a fraction of the critical density; ',0' means today (about 0.7). Say it: “omega vac”.
Ωr\Omega_{r} radiation density parameter
The radiation density as a fraction of the critical density (about 9 × 10⁻⁵ today). Say it: “omega r”.
Ωk,0\Omega_{k,0} curvature density parameter
The curvature term written like a density parameter, −kc²/(a²H²); ',0' means today. Zero for a flat universe. Say it: “omega k”.
dLd_{L} luminosity distance
The luminosity distance: the distance that would make the observed brightness obey the inverse-square law. Say it: “d L”.

Tap or click a symbol in a formula to see the same explanation. All symbols: symbols and notation.

Where this leads

Pages that build on this

Check your understanding

  1. What did the 1998 supernova teams actually measure?
  2. Why does the received power carry a factor 1/(1+z)21/(1+z)^2?
  3. In the closed-universe metric written with the angle ψ\psi, what is the area of the sphere around the source on which the detector sits, measured today?
  4. For a radial light ray in ds2=−c2dt2+a~2(t) dψ2+…ds^2 = -c^2dt^2 + \tilde a^2(t)\,d\psi^2 + \ldots, what is dψ/dtd\psi/dt?
  5. At z=0.5z = 0.5, a flat universe with Ωm=0.3\Omega_m = 0.3 and Ωvac=0.7\Omega_\text{vac} = 0.7 gives a luminosity distance about 20% larger than a flat matter-only universe with the same H0H_0. How bright does a standard candle look compared with the matter-only prediction?
0 of 5 answered

Flashcards

1 / 4