Lecture 5 Observational Techniques: The Distance Ladder

Updates - 27/4/20:

  • Clarified in Sec. 5.7 that you are not expected to memorise the numbers in the Key Project error budget.
  • Clarified in Sec. 5.7 that the details of the Bayesian analysis of Gaia parallex uncertainties is included only for those who want to look into it further.

In Chapters 2, 3 and 4, we mainly considered the theoretical aspects of Cosmology. In this section we will investigate the link between the models that describe the expanding Universe and the observations that allow us to measure their parameters.

5.1 Equations to Observations

As we have seen in the previous sections, the Friedman equation is the foundation of Cosmology: \[\begin{equation} \left(\dfrac{\dot{a}}{a}\right)^2 + \dfrac{kc^2}{a^2} = \dfrac{8\pi G}{3}\rho \tag{5.1} \end{equation}\] In Section 3.1, we showed that the scale factor, \(a\), evolves as \[\begin{equation} a(t) = \left(\dfrac{t}{t_0}\right)^{2/3} \tag{5.2} \end{equation}\] in a matter-dominated Universe, leading us to the Hubble parameter, \(H(t)\): \[\begin{equation} H(t) \equiv \dfrac{\dot{a}}{a} = \dfrac{2}{3t} \tag{5.3} \end{equation}\]

An accurate and precise measurement of \(H(t)\) will tell us how \(a\) changes over time. Additionally, measuring \(H(t_0)\), more commonly written as \(H_0\) (the Hubble constant), gives us an estimate of the age of the Universe.

The longest established route for measuring \(H_0\) is via the distance ladder. The distance ladder ties together various standard candles in order to estimate the distances of increasingly distant objects.

5.2 Standard candle measurements of \(H_0\)

Standard candles are astronomical objects of known intrinsic brightness. Examples of standard candles are objects such as type Ia supernovae, or Cepheid variable stars.

Cepheids are standard candles as they have a well-defined relationship between their intrinsic brightness (or absolute magnitude) and their period of variation. This relation is known as the Leavitt law (Leavitt and Pickering 1912), but you may also see it referred to as the period-luminosity relation. Figure 5.1 shows the Leavitt laws for a selection of Cepheids in the Small Magellanic Cloud.

Leavitt laws for Cepheids in the Small Magellanic Cloud observed at optical (purple, bottom) to mid-infrared (pink, top) wavelengths. There is a clear relation between period and luminosity, with the relation tightening and steepening as we move to longer wavelengths. The relations are offset so these changes can be seen. Adapted from @2016Scowcroft.

Figure 5.1: Leavitt laws for Cepheids in the Small Magellanic Cloud observed at optical (purple, bottom) to mid-infrared (pink, top) wavelengths. There is a clear relation between period and luminosity, with the relation tightening and steepening as we move to longer wavelengths. The relations are offset so these changes can be seen. Adapted from Scowcroft et al. (2016).

By comparing the apparent magnitude of a standard candle with it’s absolute magnitude we can obtain its distance: \[\begin{equation} \mu = m - M = 5 \log(d) - 5 \tag{5.4} \end{equation}\] where \(\mu\) is the distance modulus, \(m\) is the apparent magnitude, \(M\) is the absolute magnitude, and \(d\) is the distance in parsecs. In the case of Cepheids, we observe the star’s apparent magnitude and calculate its absolute magnitude using the Leavitt law. The Leavitt law takes the form \[\begin{equation} M_{\lambda} = a_{\lambda} \log(P) + b_{\lambda} \tag{5.5} \end{equation}\]

where \(M_{\lambda}\) is the absolute magnitude at wavelength \(\lambda\), and \(a_{\lambda}\) and \(b_{\lambda}\) are the slope and intercept of the relation. Figure 5.1 shows how the slope (\(a\)) and intercept (\(b\)) change with wavelength, with \(a\) becoming increasingly negative as we move to longer wavelengths.

Moving to longer wavelengths to measure Cepheid distances is advantageous for several reasons. First, effects of reddening and extinction on apparent magnitudes are dramatically reduced compared to optical wavelengths. Second, the Leavitt law has smaller dispersion – i.e. the standard deviation, \(\sigma\) of the points around the Leavitt law is smaller, and the amplitude of their light curves are smaller (see Figure 5.2). This means that distances obtained from redder wavelengths are more accurate than those from shorter wavelength observations.

Multiwavelength light curves for Milky Way Cepheid CD Cyg. The light curves become more sinusoidal and smaller in amplitude as you move from blue to red wavelengths.

Figure 5.2: Multiwavelength light curves for Milky Way Cepheid CD Cyg. The light curves become more sinusoidal and smaller in amplitude as you move from blue to red wavelengths.

5.3 Calibration of standard candles

Standard candles are useful as they are objects of known intrinsic luminosity. But how do we know their intrinsic luminosity? To determine the intrinsic luminosity of a standard candle, we must have a calibration sample for which we know the objects distance in addition to their apparent magnitudes. In the case of variable stars, this is done using parallax measurements of objects within the Milky way.

5.4 Parallax calibration of standard candles

Parallax is a geometric distance determination. Distances are estimated by measuring how the positions of objects change over time with respect to a reference frame of distant, background objects.

Schematic of the parallax method (not to scale). Distances are determined by measuring the apparent angular shift of nearby objects compared to an seemingly static background of distant objects. Observations are separated by six months (i.e. opposite sides of the Earth's orbit) in order to maximise the baseline of the observations.

Figure 5.3: Schematic of the parallax method (not to scale). Distances are determined by measuring the apparent angular shift of nearby objects compared to an seemingly static background of distant objects. Observations are separated by six months (i.e. opposite sides of the Earth’s orbit) in order to maximise the baseline of the observations.

Figure 5.3 illustrates how parallax distance determinations are made. The positions of target stars relative to distant, background objects are measured. Observations are separated by six months, in order to maximise the baseline of the observations. The parallax angle, \(\varpi = \left(\theta / 2\right)\) is used to calculate the distance via simple trigonometry:

\[\begin{equation} d = \dfrac{r}{\tan(\varpi)} \tag{5.6} \end{equation}\]

As the parallax angle, \(\varpi\), is extremely small, and the baseline length is much smaller than the distance of the target object \(\left(\text{i.e. }r \ll d \right)\), Eq. (5.6) can be simplified to:

\[\begin{equation} d = \dfrac{1}{\varpi} \tag{5.7} \end{equation}\]

where \(d\) is the distance of the target in parsecs.

In the case of the distance ladder, parallax distances of Cepheids are used to fix the intercept, or zero-point, of the Leavitt law. Until recently only 10 Cepheids with high-precision parallaxes were available for this calibration. However, the Gaia mission promises to dramatically improve this situation. Gaia is a European Space Agency mission which is measuring the parallaxes of over one billion stars in the Milky Way and its nearest neighbours to extremely high precision. The precision of Gaia is such that if we put it on top of Buckingham Palace it could resolve a human hair on top of the Empire State Building, over 5,500 km away. Gaia’s parallax catalogue will contain many thousands of Cepheids that could potentially be included in the Leavitt law calibration.

5.5 The Distance Ladder

With the Leavitt law calibrated, we can now measure distances to objects outside of our own Galaxy. However, although Cepheids are bright stars, we cannot observe them in very distant galaxies. In order to measure distances to the furthest objects we must use different techniques, each calibrated to our parallax zero-point. Tying together different methods to estimate distances to progressively further objects is known as the distance ladder. Figure 5.4 shows how the techniques tie together to measure \(H_0\).

The distance ladder. The bottom rung of the distance ladder is established by measuring parallax distances to Cepheids within our own Galaxy. Cepheids are measured in increasingly distant galaxies in order to calibrate the absolute scale of techniques such as the Tully Fisher relation and Type Ia supernovae (SN Ia). The Hubble constant is measured using the recession velocities (redshift) and distances of SN Ia in the furthest galaxies, where contribution of peculiar velocities is minimized.

Figure 5.4: The distance ladder. The bottom rung of the distance ladder is established by measuring parallax distances to Cepheids within our own Galaxy. Cepheids are measured in increasingly distant galaxies in order to calibrate the absolute scale of techniques such as the Tully Fisher relation and Type Ia supernovae (SN Ia). The Hubble constant is measured using the recession velocities (redshift) and distances of SN Ia in the furthest galaxies, where contribution of peculiar velocities is minimized.

To get an accurate measurement of \(H_0\) we must use objects that are as far away as possible. If we were to use only galaxies that were nearby in our estimation then the velocities we measure would be dominated by peculiar motions. Peculiar motions are caused by gravitational interaction with other nearby objects. For very distant galaxies – those deemed to be in the Hubble flow – the galaxy’s velocity is dominated by its recession velocity, and the contribution from the peculiar velocity is negligible.

5.6 Measurements of \(H_0\) from the distance ladder

The distance ladder is one of the oldest and most frequently used techniques to measure \(H_0\). Figure 5.5 shows how our measurements of \(H_0\) have converged over time, from the early measurements in the 1930’s, to the high-precision measurements of the 2010’s and 2020’s. From Figure 5.5 it looks like we may have reached a consensus on the value of \(H_0\). To confirm that, we need to consider the uncertainties on our measurements.

Measurements of $H_0$ as a function of time. Hubble's original measurement [@1929Hubble] was around 500 km s$^{-1}$ Mpc$^{-1}$. Over the last one-hundred years, we appear to have reached a consensus that $H_{0}$ is approximately 70 km s$^{-1}$ Mpc$^{-1}$, although recent high precision measurements have brought that consensus into question.

Figure 5.5: Measurements of \(H_0\) as a function of time. Hubble’s original measurement (Hubble 1929) was around 500 km s\(^{-1}\) Mpc\(^{-1}\). Over the last one-hundred years, we appear to have reached a consensus that \(H_{0}\) is approximately 70 km s\(^{-1}\) Mpc\(^{-1}\), although recent high precision measurements have brought that consensus into question.

5.7 Uncertainties in the distance ladder

There are two types of uncertainties that can be present in a result – systematic and random. A comparison of the two and their effects on the precision and accuracy of a measurement are shown in Table 5.1 and Figure 5.6.

Table 5.1: Comparison of systematic and random uncertainties
Systematic Random
Affects all data points the same way Affects each point differently
Often assumed to have Gaussian distribution
Causes offsets Causes dispersion
Affects accuracy Affects precision
Doesn’t decrease with \(N_{\text{obs}}\) Decreases with \(\sqrt{N_{\text{obs}}}\)
The effect of systematic and random uncertainties on accuracy and precision.

Figure 5.6: The effect of systematic and random uncertainties on accuracy and precision.

The distance ladder ties together several techniques to measure distances to the furthest galaxies, hence uncertainties on each “rung” of the ladder propagate through to our final measurements. Large uncertainties at the bottom will lead to unreliable measurements at the top, in a similar way to having an unstable foundation for the bottom of a ladder leads to an unstable, wobbly ladder that you really wouldn’t want to stand at the top of.

A major step forward in measuring a high-precision value of \(H_0\) was made in 2001 with the publication of the results from the HST Key Project to Measure the Hubble Constant Freedman et al. (2001), who found \(H_0\) = 71 \(\pm\) 2 (rand.) \(\pm\) 6 (sys.) km s\(^{-1}\) Mpc\(^{-1}\). Table 5.2 shows how different components of systematic uncertainty combine in their final result.

You are not expected to memorise these numbers. Table 5.2 is provided so you can see the variety of effects that contribute to \(H_0\) measurements. Some are obvious; the calibration of the camera is something you would expect to see in the error budget for an observational experiment. Others are less obvious,For example, the bias in the Leavitt Law comes from the fact that the Cepheids observed may not fully sample the true distribution of the Leavitt law. This means that when we derive a distance to a galaxy using this biased sample, that distance measurement will also be biased.

Table 5.2: Overall systematic uncertainties in the Hubble constant. From Freedman et al. (2001).
Source of uncertainty Method of estimation Error (%)
LMC zero point Error on mean from Cepheids, TRGB, SN 1987A, red clump, eclipsing binaries \(\pm\) 5
HST camera calibration Tie-in to Galactic star clusters \(\pm\) 3.5
Reddening Limits from photometry \(\pm\) 1
Metallicity Observational and theoretical constraints \(\pm\) 4
Bias in Leavitt law Short-end period cut-off \(\pm\) 1
Crowding Artificial star experiments +5, -0
Large scale flows Limits from SN Ia, CMB \(\pm\) 5

In the years since the HST Key Project, these uncertainties have been further reduced. One of the highest precision measurements of \(H_0\) to date comes from the SH\(_{0}\)ES experiment (Riess et al. 2016), who found \(H_0= 73.24 \pm 1.74\) km s\(^{-1}\)~Mpc\(^{-1}\). Figure 5.7 shows the breakdown of their uncertainty budget.

Uncertainty budget from the *SH$_{0}$ES* experiment [@2016Riess].

Figure 5.7: Uncertainty budget from the SH\(_{0}\)ES experiment (Riess et al. 2016).

One of the largest uncertainties remaining in distance ladder measurements of \(H_0\) comes from the bottom rung of the ladder. In the Key Project, the ladder was tied to our nearest galaxy, the Large Magellanic Cloud (LMC). SH\(_{0}\)ES reduced this uncertainty by tying their calibration to the Cepheids in the LMC (which can all be assumed to be at the same distance), Cepheids with parallaxes in the Milky Way (MW), and Cepheids in the mega-maser host galaxy NGC4258.

:::fyi Astrophysical masers are sources of stimulated microwave emission, similar to lasers. Mega-masers are astrophysical masers with very high isotropic luminosities. Mega-masers are excellent geometric distance indicators. By assuming that the maser is travelling round a galaxy with a Keplerian orbit we can calculate the intrinsic size of the system. Comparing this with its apparent size gives us an estimate of the distance. :::

Until recently, only a handful of MW Cepheids had high-precision parallax measurements. The ESA Gaia mission will change this in the coming years by measuring high-precision parallaxes for over 1~billion stars in the MW. The Gaia sample will include several thousand Cepheids which can then be used to anchor the distance ladder. Figure 5.8 shows the sky as seen by Gaia. Each star in this image has a parallax measurement, including those in the Magellanic Clouds (seen to the bottom right of the Galactic Plane).

The *Gaia* sky in colour. Credit: ESA/Gaia/DPAC, CC BY-SA 3.0 IGO.

Figure 5.8: The Gaia sky in colour. Credit: ESA/Gaia/DPAC, CC BY-SA 3.0 IGO.

In order to properly calibrate the Leavitt law and the uncertainties it contributes to the distance ladder, the Gaia Cepheid parallaxes must be treated extremely carefully. When using small samples of Cepheids with reasonably small uncertainties, it is appropriate to estimate the parallax distance (\(r\)) as \[\begin{equation} r = \dfrac{1}{\varpi} \tag{5.8} \end{equation}\]

where \(\varpi\) is the parallax angle. However, if we wish to ensure an unbiased calibration of the Leavitt law, we must not truncate the Cepheid sample with limits on (a) parallax uncertainty (\(\sigma_{\varpi}\)) or (b) apparent magnitude. As there is an inverse relation between parallax and distance, the probability distribution of distances for a given parallax angle is skewed, becoming more so as the fractional uncertainty \({f}_{\mathrm{obs}} = \sigma_{\varpi}/\varpi\) increases. This is illustrated in Figure 5.9.

Left panel is the probability $P(\varpi | r,{\sigma }_{\varpi })$ as a function of $\varpi$. The middle panel is $P(\varpi | r,{\sigma }_{\varpi })$ as a function of $r$ for an object with ${f}_{obs}=0.2$, and the right panel is as the middle panel but for ${f}_{obs}=0.5$, where. ${f}_{obs} = \sigma_{\varpi}/\varpi$. The shaded areas indicate the 2$\sigma$ credible interval about $1/r$ and the corresponding transformed credible interval. As $f_{obs}$ increases, the probability distribution for $r$ becomes more skewed, such that stars with larger parallax uncertainties become more likely to have overestimated $r$ than to have underestimated $r$. From @2016BJ

Figure 5.9: Left panel is the probability \(P(\varpi | r,{\sigma }_{\varpi })\) as a function of \(\varpi\). The middle panel is \(P(\varpi | r,{\sigma }_{\varpi })\) as a function of \(r\) for an object with \({f}_{obs}=0.2\), and the right panel is as the middle panel but for \({f}_{obs}=0.5\), where. \({f}_{obs} = \sigma_{\varpi}/\varpi\). The shaded areas indicate the 2\(\sigma\) credible interval about \(1/r\) and the corresponding transformed credible interval. As \(f_{obs}\) increases, the probability distribution for \(r\) becomes more skewed, such that stars with larger parallax uncertainties become more likely to have overestimated \(r\) than to have underestimated \(r\). From Astraatmadja and Bailer-Jones (2016a)

This skewness must be taken into account in order to prevent systematic uncertainties being introduced into the Leavitt law calibration. Therefore, rather than use the simple parallax equation given in Eq. (5.8), we must use a full Bayesian analysis to estimate distances. A full discussion of the application of Bayesian statistics to parallax measurements is beyond the scope of this course. However, if you are interested in the methodology, the series of papers by Astraatmadja & Bailer-Jones gives an excellent introduction to the topic (Bailer-Jones 2015; Astraatmadja and Bailer-Jones 2016a; Astraatmadja and Bailer-Jones 2016b).

References

Astraatmadja, Tri L., and Coryn A. L. Bailer-Jones. 2016a. “Estimating Distances from Parallaxes. II. Performance of Bayesian Distance Estimators on a Gaia-like Catalogue.” ApJ 832 (2): 137.

Astraatmadja, Tri L., and Coryn A. L. Bailer-Jones. 2016a. “Estimating Distances from Parallaxes. II. Performance of Bayesian Distance Estimators on a Gaia-like Catalogue.” ApJ 832 (2): 137.

2016b. “Estimating Distances from Parallaxes. III. Distances of Two Million Stars in the Gaia DR1 Catalogue.” ApJ 833 (1): 119.

Bailer-Jones, Coryn A. L. 2015. “Estimating Distances from Parallaxes.” PASP 127 (956): 994.

Freedman, Wendy L, Barry F Madore, Brad K Gibson, Laura Ferrarese, Daniel D Kelson, Shoko Sakai, Jeremy R Mould, et al. 2001. “Final Results from the Hubble Space Telescope Key Project to Measure the Hubble Constant.” ApJ 553 (1): 47.

Hubble, Edwin. 1929. “A Relation between Distance and Radial Velocity among Extra-Galactic Nebulae.” Proceedings of the National Academy of Science 15: 168.

Leavitt, Henrietta S., and Edward C. Pickering. 1912. “Periods of 25 Variable Stars in the Small Magellanic Cloud.” Harvard College Observatory Circular 173: 1.

Riess, Adam G., Lucas M. Macri, Samantha L. Hoffmann, Dan Scolnic, Stefano Casertano, Alexei V. Filippenko, Brad E. Tucker, et al. 2016. “A 2.4% Determination of the Local Value of the Hubble Constant.” ApJ 826 (1): 56.

Scowcroft, Victoria, Wendy L. Freedman, Barry F. Madore, Andy Monson, S. E. Persson, Jeff Rich, Mark Seibert, and Jane R. Rigby. 2016. “The Carnegie Hubble Program: The Distance and Structure of the SMC as Revealed by Mid-infrared Observations of Cepheids.” ApJ 816: 49.