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.
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 z? 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.
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 )
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>0) and introduces an angle ψ — which, traced back, is the angle from the w-axis used when the closed universe was first built as a 3-sphere in four dimensions:
Radial angle
sinψ=kr⟹dψ=cosψkdr=1−kr2kdr
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
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
The tilde only distinguishes the new scale factor from the old one; starting from this metric, you could drop it. A student asks whether a and a~ have different units. They do. In Guth's convention r is a coordinate distance in notches, so a is in metres per notch; ψ is an angle and naturally dimensionless, so a~ is simply a length.
Put the source at the centre (ψ=0), and imagine a sphere around it, drawn in comoving coordinates so that our detector, of area A, lies on its surface at ψ=ψ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 θ and ϕ 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ψ. We evaluate it today, t0, because the detector has area A 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
4πa~2(t0)sin2ψdA
In words
A over 4 pi a tilde squared (t sub 0 ) sine squared psi sub d
The redshift costs two more factors. Any clock at the source appears slowed by 1+zs, where zs is the redshift of the source. One clock is the photon's frequency, so each photon arrives with its energy reduced by 1+zs; the other is the rate at which photons arrive, reduced by the same factor:
Redshift
1+zs=a(te)a(t0)
In words
1 plus z sub s equals a(t sub 0 ) over a(t sub e )
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
For a closed universe k>0, so Ω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 it divides.]
The angle. Light travels on null geodesics. For a radial ray, ds2=0:
Radial light ray
−c2dt2+a~2(t)dψ2=0⟹dtdψ=a~(t)c
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)
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=∫tet0cdt/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) and differentiate:
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)
The Friedmann equation from the first half of the lecture, H2=H02F(x)/x4 with x=1/(1+z), supplies H(z). Substituting it and the expression for a~(t0):
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
Put this ψ and a~(t0) into the formula for J, and the expected intensity depends only on the present Ω's, H0, 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.
Guth works out only the closed case. The other two follow from it (this section is the page's own addition):
Open (k<0, Ωk,0>0): the same steps with sinhψ=∣k∣r give the same formula with sinψd replaced by sinhψd, and −Ωk,0 by Ωk,0.
Flat (Ωk,0→0): ψd becomes tiny and sinψd≈ψd, so a~(t0)sinψd→(c/H0)∫0zsdz/E(z), where E(z)=H(z)/H0 is the square root in the integral above.
Astronomers package all three cases into a luminosity distancedL, defined so that J=P/(4πdL2) looks like the inverse-square law of a static, flat space. Comparing with the formula above, dL=(1+zs)a~(t0)sinψd.
Loading graph…
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.5. Take H0=70 km/s/Mpc and compare two flat universes: (A) matter only, (B) Ωm=0.3, Ω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, Ωvac,0=0.70 (so Ωk,0=−0.01), changes model (B).
(A) Matter only
E(z)=(1+z)3/2, so ∫00.5dz/E=2(1−1/1.5)=0.3670. Then H0dL/c=1.5×0.3670=0.5505. With c/H0=4283 Mpc, dL=2358 Mpc.
(B) Matter + vacuum energy
E(z)=0.3(1+z)3+0.7. Numerically, ∫00.5dz/E=0.4410, so H0dL/c=0.6615 and dL=2833 Mpc.
How much fainter
J∝1/dL2, so JB/JA=(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 radians — a small angle, as the closed universe is huge compared with the distance to the source. a~(t0)=(c/H0)/0.1, so H0dL/c=1.5×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.
The 1998 teams measured the intensity of type Ia supernovae against redshift; for standard candles, intensity measures distance.
In a closed universe, with sinψ=kr and a~=a/k, the light spreads over a sphere of area 4πa~2(t0)sin2ψd, and redshift costs two factors of 1+zs: J=P/[4π(1+zs)2a~2(t0)sin2ψd].
a~(t0)=cH0−1/−Ωk,0, and the null-geodesic condition dψ/dt=c/a~ gives ψs=a~(t0)1∫0zscdz/H(z), with 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.
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”.
What did the 1998 supernova teams actually measure?
Why does the received power carry a factor 1/(1+z)2?
In the closed-universe metric written with the angle ψ, what is the area of the sphere around the source on which the detector sits, measured today?
For a radial light ray in ds2=−c2dt2+a~2(t)dψ2+…, what is dψ/dt?
At z=0.5, a flat universe with Ωm=0.3 and Ωvac=0.7 gives a luminosity distance about 20% larger than a flat matter-only universe with the same H0. How bright does a standard candle look compared with the matter-only prediction?