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.
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.
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.
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\).
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.
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.
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 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.
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.
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).
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.
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.