Chapter 6

Deceleration Parameter

Now we transition to observations. The universe, it turns out, is large — and a little messy. Measuring a(t)a(t) is no trivial task, especially when dark energy is added to the mix. In simpler times, no one gave dark energy any mind, though, and a(t)a(t) could be determined entirely by local measurements. This is because, without dark energy, curvature is determined by the current matter density (since radiation is so much less important), which leads to only 3 options: open, closed, or flat/coasting universes. These possibilities can be distinguished by measuring the current value of the deceleration of the expansion of space, or a¨\ddot{a}. Without knowing the exact equation for a(t)a(t), the local (t=t0t=t_0) form of the equation can be approximated using a Taylor expansion. That leads to this expression (only valid for tt0t \approx t_0, becoming less correct for other times):

a(t)1+H0(tt0)12q0H02(tt0)2;H0a˙at=t0  &  q0(a¨aa˙2)t=t0.a(t) \approx 1 + H_0 (t-t_0) - \frac{1}{2}q_0 H_0^2 (t-t_0)^2\, ; \\H_0 \equiv \frac{\dot{a}}{a}\bigg\vert_{t=t_0}~~\&~~q_0 \equiv -\left( \frac{\ddot{a}a}{\dot{a}^2} \right)_{t=t_0}\, .

q0q_0 is the deceleration parameter (so named because it is dimensionless and is the negative of acceleration). This definition can also be put in terms of density parameters:

q0=12i=1NΩi,0(1+3wi)=Ωr,0+12Ωm,0ΩΛ,0\begin{array}{lll} q_0 &=& \displaystyle\frac{1}{2}\sum^N_{i=1}\Omega_{i,0} (1+3w_i)\\ &=& \Omega_{\rm r,0} +\frac{1}{2}\Omega_{\rm m,0} - \Omega_{\Lambda,0} \end{array}

Proper Distance

Using this approximation of a(t)a(t) (or comparable expansion for 1/a(t)1/a(t)) and inserting this into the expression for redshift z=1/a(te)1z = 1/a(t_e) - 1, the proper distance can be expressed in terms of zz, H0H_0, and q0q_0:

dp(t0)czH0[11+q02z].d_p(t_0) \approx \frac{cz}{H_0}\left[ 1 - \frac{1+q_0}{2}z \right]\, .

Here we see that Hubble's law (cz=H0dcz=H_0d) gets modified (to second order in zz) by q0q_0 (the classical Hubble law is just the first term in this approximation, and more terms could get added z2,z3,\propto z^2, \propto z^3, etc. For this approximation to be valid, z1z \ll 1 (otherwise the higher order terms aren't negligible).


Measurable Distances

This is great, but proper distances are not what we measure as they are instantaneous distances — light travels at a finite speed, during which aa changes, affecting the ways in which we can estimate distances with observations. There are various ways to estimate distances, the best of which are geometric (e.g., parallax), assuming we know (or can precisely calculate) the size of the short side of a triangle (Earth's semimajor axis in the case of parallax). These cases are rare in practice; most are based are knowing the intrinsic brightness, or luminosity, of an object and measuring its flux.

Luminosity Distance

Distances can be estimated if we know how intrinsically luminous an object is, i.e., if it is a standard candle. This distance is just that assumed by the inverse square law:

dL(L4πf)1/2.d_L \equiv \left( \frac{L}{4\pi f} \right)^{1/2}\, .

The observed flux from a source, however, depends on the curvature of the universe, or Sκ(r)S_\kappa(r). This is because the surface area of a sphere is different depending on the curvature (Asph=4πSκ(r)2A_{\rm sph} = 4\pi S_\kappa(r)^2) AND because light is redshifted, losing energy (1+z)1\propto (1+z)^{-1}, AND because the distance between photons increase with the expansion of space, also (1+z)1\propto (1+z)^{-1}. So the actual measured flux (energy/time/area) is then

f=L4πSκ(r)2(1+z)2f = \frac{L}{4\pi S_\kappa(r)^2(1+z)^2}

so dL=Sκ(r)(1+z)\boxed{d_L = S_\kappa(r)(1+z)} by definition. In a flat universe (like ours), Sκ(r)=r=dp(t0)S_\kappa(r) = r = d_p(t_0), so dL=dp(t0)(1+z)d_L = d_p(t_0) (1+z). Using the approximate expression for dp(t0)d_p(t_0) above,

dLczH0(1+1q02z)\boxed{d_L \approx \frac{cz}{H_0} \left(1 + \frac{1-q_0}{2}z \right)}

Angular-diameter Distance

Distances can also be estimated if we know the intrinsic size of an object, i.e., it is a standard ruler. This distance is just that assumed by the small-angle formula:

dAδθd_A \equiv \frac{\ell}{\delta \theta}

where \ell is the length of the object and δθ\delta \theta is its observed angular size. The angular size depends on when the light was emitted, and =a(te)Sκ(r)δθ\ell = a(t_e) S_\kappa(r) \delta \theta, so dA=Sκ(r)/(1+z)\boxed{d_A = S_\kappa (r) / (1+z)}. In a flat universe, Sκ(r)=r=dp(t0)S_\kappa(r) = r = d_p(t_0), so dA=dp(t0)/(1+z)d_A = d_p(t_0) / (1+z). Using the approximate expression for dp(t0)d_p(t_0) above,

dAczH0(13+q02z)\boxed{d_A \approx \frac{cz}{H_0} \left(1 - \frac{3+q_0}{2}z \right)}

As can be seen comparing dLd_L and dAd_A in a flat universe, they are related to each other by

dA(1+z)=dp(t0)=dL/(1+z)d_A (1+z) = d_p(t_0) = d_L / (1+z)

Luminosity distances, but using magnitudes (ugh!)

Astronomers still use magnitudes, first defined by ancient Greek astronomers to rank star brightnesses. Due to the semi-logarithmic sensitivity of human eyes, this means you have to deal with a garbage relation between fluxes (a sane unit of measure) and observed magnitudes, which are properly called apparent magnitudes:

m2.5log10(f/fx)m \equiv -2.5 \log_{10}(f/f_x)

where fx=2.53×108 W m2f_x = 2.53 \times 10^{-8}~{\rm W~m}^{-2} is just a constant defined by convention (and can be different in different magnitude systems). Similarly, a corresponding equation for luminosities can be written, which are called absolute magnitudes:

M2.5log10(L/Lx)M \equiv-2.5 \log_{10}(L/L_x)

The constant Lx=78.7 LL_x = 78.7~{\rm L}_\odot is defined as the luminosity a star (or whatever) with a flux of fxf_x would have if its luminosity distance dL=10 pcd_L = 10~{\rm pc}. Using our definition of dLd_L above,

M=2.5log10(4πdL2f/(4π(10 pc2fx)=2.5log10([dL/(10 pc)]2[f/fx])=2.5log10(f/fx)5log10(dL/10 pc)=m5log10(dL/10 pc)M = -2.5\log_{10}(4\pi d_L^2 f / (4\pi (10~{\rm pc}^2 f_x) = -2.5\log_{10}([d_L/(10~{\rm pc})]^2[f/f_x]) \\ = -2.5\log_{10}(f/f_x) - 5\log_{10}(d_L/10~{\rm pc}) = m - 5\log_{10}(d_L/10~{\rm pc})

Since galaxies are separated by Mpc, not pc, it's convenient to change the units of dLd_L, equivalent to multiplying the second term by a factor of log10(105)=5\log_{10}(10^5) = 5, so to keep the equation balanced, we get a version of this equation called the distance modulus:

mM=5log10(dL1 Mpc)+25m-M = 5\log_{10}\left(\frac{d_L}{1~{\rm Mpc}} \right) + 25

This is nothing more than the inverse square law, but in a semi-logarithmic form for no good reason. Think of it like a secret handshake in a secret society — an overly-complicated way to say hi to a fellow astronomer. Plugging this into our approximate expression for dLd_L above yields the truly awful equation (but the one used IRL):

mM43.235log10(H068 km s1 Mpc1)+5log10z+1.086(1q0)zm-M \approx 43.23 - 5\log_{10}\left(\frac{H_0}{68~{\rm km~s}^{-1}~{\rm Mpc}^{-1}} \right) + 5\log_{10}z+1.086(1-q_0)z