Data Analysis problem set
AdvancedIOAA 2017 | Phuket, Thailand | Data Analysis D1 | 75 pt

2017 D1 | Calibrating the distance ladder to the LMC

An accurate trigonometric parallax calibration for Galactic Cepheids has long been sought but is hard to achieve. All known classical Cepheids are more than 250 pc away, so a 10% distance needs parallax uncertainties of 0.2 milliarcsec or better, which requires space-based observations. The Fine Guidance Sensor on HST measured parallaxes better than 10% for 10 Cepheids, with periods from 3.7 to 35.6 days and distances of about 300 to 560 pc.

Table 1 gives periods PP, mean magnitudes in V, K and I, extinctions AVA_V and AKA_K, and parallaxes with their uncertainties. The apparent magnitudes have negligible uncertainty.

(D1.1) The period–luminosity (PL) relation is L∝PβL\propto P^\beta; in practice it is written in terms of absolute magnitude (the "PL relation"). Use Table 1 to plot a suitable linear graph to derive the PL relation in V and in K, on separate graphs. Determine the slope of the best-fit line. (You may find Δ(log⁡10x)≈Δx/(xln⁡10)\Delta(\log_{10}x)\approx\Delta x/(x\ln10) useful.) (36.5)

The reddening is a function of wavelength and varies between Cepheids, so the PL relation is really a period–luminosity–colour relation. The reddening-free "Wesenheit" magnitude removes it using the star's own colour:

WVI=V−AVE(V−I)(V−I)=V−RV(V−I),RV=2.45 here.W_{VI} = V - \frac{A_V}{E(V-I)}(V-I) = V - R_V(V-I), \qquad R_V = 2.45 \text{ here}.

(D1.2) Plot and derive the reddening-free PL relation using WVIW_{VI}. Estimate its slope and the slope's uncertainty. (14.5)

(D1.3) Use the PL relations from (D1.1) and (D1.2) with the LMC Cepheids in Table 2 (with extinction-corrected ⟨Vcorr⟩\langle V_{\rm corr}\rangle and ⟨WVI⟩\langle W_{VI}\rangle). Estimate the distance modulus of each star, then the distance to the LMC (in pc) and its standard deviation for each band. Are the two distances statistically different (YES/NO)? Are the standard deviations different (YES/NO)? Which band is more accurate for the LMC distance? (24)

Table 1 | Galactic Cepheids with HST parallaxes
StarP (d)⟨V⟩⟨K⟩AVA_VAKA_K⟨I⟩π (mas)σπ (mas)
RT Aur3.7285.4643.9250.200.024.7782.400.19
FF Aql4.4715.3723.4650.640.084.5102.810.18
X Sgr7.0134.5562.5570.580.073.6613.000.18
ζ Gem10.1513.9112.0970.060.013.0852.780.18
l Car35.5513.7321.0710.520.062.5572.010.20
Table 2 | Cepheids in the Large Magellanic Cloud
StarP (d)⟨V_corr⟩ (mag)⟨W_VI⟩ (mag)
HV121992.6316.0814.56
HV122032.9515.9314.40
HV128169.1014.3012.80
HV89930.9013.0710.97
HV225739.3612.8610.54

Saved in this browser only, no account needed.

Reading resources

Material that covers the physics behind this problem.