Observations of Field and Cluster RR Lyrae with Spitzer. Towards High Precision Distances with Population II Stellar Tracers

TL;DR

This paper discusses calibrating RR Lyrae stars' period-luminosity-metallicity relation using Spitzer observations to improve distance measurements with Population II stellar tracers.

Contribution

It introduces a new calibration of RR Lyrae stars' relation using high-precision infrared data from Spitzer, enhancing distance measurement accuracy.

Findings

Calibrated RR Lyrae period-luminosity-metallicity relation.

Improved distance estimates to Galactic halo and globular clusters.

Enhanced understanding of Population II stellar tracers.

Abstract

We present our project to calibrate the RR Lyrae period-luminosity-metallicity relation using a sample of Galactic calibrators in the halo and globular clusters.

Figures

Equations

$$M_{\lambda }=a_{\lambda }+b_{\lambda }\cdot \log P+c_{\lambda }\cdot \text{[Fe/H]}$$

$$m_{\lambda ,i}-{\mu }_i=a_{\lambda }+b_{\lambda }\cdot \log P_i+c_{\lambda }\cdot {\text{[Fe/H]}}_i+A_{V,i}\cdot [A_{\lambda }/A_V]$$

Full text

1 Introduction

The availability of reliable, high-precision, RR Lyræ distances could be a game-changer not just for tracing old stellar populations in the Milky Way and other Local Group galaxies, but also for renewed efforts to establish an independent cosmological distance scale entirely based on Population II indicators 2016ApJ…832..210B . This will be crucial to resolve the current $\geq 3\sigma$ tension existing between different measurements of the Hubble constant and other cosmological parameters. Measurements obtained with methods based on Population I stellar indicators, anchored on Cepheids, are at odd with determinations based on the Cosmic Background Radiation and Baryonic Acoustic Oscillations 2016ApJ…826…56R . If confirmed, this tension could be an indication of unknown physics, missing in the current formulation of $\Lambda$CDM cosmological models.

Historically, the adoption of RR Lyræ as precision standard candles has been limited by the lack of a proper period-luminosity (PL) relation. Their visible magnitude only shows a shallow dependence from their pulsation period, and the magnitude vs. metallicity relation commonly used to derive their photometric parallaxes is characterized by a large intrinsic scatter ($>5$%). The situation however changes at near- and mid-infrared wavelengths (NIR and MIR thereafter), for three main reasons: (i) in the NIR the slope in the RR Lyræ PL relation steadily increases (see Figure 1); (ii) the intrinsic scatter due to evolutionary effects and temperature dependence is greatly reduced in the infrared, because of the narrowing of the instability strip width, and because in the MIR the brightness variations are mainly driven by the radius changes; (iii) in the NIR the effect of reddening is one order of magnitude less than in the visible, with an extra factor of 3 reduction in the MIR 2016MmSAIt.XXX..XXXB . These predictions were recently confirmed, observationally, both in NIR 2004ApJ…610..269D and MIR 2011ApJ…738..185K ; 2013ApJ…776..135M bands.

With the upcoming availability of large aperture infrared space telescopes (JWST in the NIR and MIR, and WFIRST that will operate at wavelengths up to 2 $\mu$m), it is time to provide a robust calibration of RR Lyræ PL relations, and characterize their dependence from metallicity. This paper presents a summary of our calibration efforts based on the MIR observation of Galactic RR Lyræ with the Spitzer Space Telescope InfraRed Array Camera (IRAC) 2004ApJS..154…10F , part of the Carnegie RR Lyræ Program (CRRP), as well as ground-based observations at visible and NIR wavelengths.

2 Calibrating the Period-Luminosity-Metallicity (PLZ) Relations

The absolute brightness of RR Lyræ can be approximated by a linearized expression of their pulsation period $P$ and iron content [Fe/H] (as proxy for their detailed abundances):

[TABLE]

where the dependence on metallicity has been factorized in the $c_{\lambda}$ term. Detailed time-dependent hydrodynamic models (such as 2015ApJ…808…50M , see also M. Marconi contribution to this volume) show that the metallicity slope $c_{\lambda}$ is only weakly dependent from wavelength, since the metallicity dependence of these stars’ bolometric corrections effectively cancels-out the corresponding changes in effective temperature, resulting in a minimal change in color as [Fe/H] is increased. Furthermore, in the MIR the stellar photosphere is well described by the Rayleigh-Jeans approximation, and lacks significant spectral features. For this reasons the PLZ metallicity slope approaches a constant value $c$. In this formulation, the period slope $b_{\lambda}$ is assumed independent from metallicity, a result supported by models; in the MIR this term converges to the value expected for brightness variations driven by changes in radius ($b_{\lambda}\simeq-2.3$ for fundamental mode pulsators, see Figure 1).

To be truly useful, RR Lyræ PLZs need however to be calibrated observationally. This can be achieved by observing a number of stars in the field and globular clusters (GCs) with known true distance modulus $\mu_{i}$, having a broad range of periods $P_{i}$ and metallicity [Fe/H]i. The average magnitudes $m_{\lambda,i}$ of each of these calibrators need to be measured at multiple wavelengths, in order to characterize their individual extinction $A_{V,i}$ (after assuming an appropriate reddening law $[A_{\lambda}/A_{V}]$), so that:

[TABLE]

Equation 2 can then be fit to derive the PLZ coefficients $a_{\lambda}$, $b_{\lambda}$ and $c_{\lambda}$ (and the nuisance parameter $A_{V,i}$ for each star) for each band in which the calibrators have been observed. To minimize the uncertainty in the best fit PLZ coefficients, the sample of RR Lyræ calibrators in our program have been selected so that they will have Gaia parallax determined better than $\sim 1$%. For all stars we are also acquiring a homogeneous set of metallicity (with precision better than $\sim 0.1$ dex), using optical and NIR high resolution spectroscopic observations.

3 Calibrating the Period Slope with Galactic Globular Clusters

From Equation 2 it is clear that even small uncertainties in the observed quantities (especially [Fe/H], rarely known better than $\sim 0.2$ dex), and variations in the reddening law will introduce a large scatter, and possibly a period-dependent bias, in the observed PLZ. This makes it difficult to simultaneously fit all PLZ coefficients, and the extinction $A_{V}$, using an heterogeneous sample of calibrator RR Lyræ located at different distance, and having different metallicity. The period slope $b_{\lambda}$, however, can be determined independently from all other parameters by observing a sample of single-metallicity GCs. The advantage of using GCs is of course that the distance of the cluster from Earth is much larger than the relative spatial distribution of the individual stars within the GC, and variable extinction in front of the cluster only contributes to a small unbiased scatter around the overall RR Lyræ PL: for RR Lyræ in a GC the period slope $b_{\lambda}$ can be accurately fit as the slope of the PL relation constructed with the reddened apparent magnitudes.

Figure 2 (adapted from 2017ApJ…XXX…..N ) presents an example based on our observations of M4 (NGC 6121). The right panel shows the PLZ for first overtone and fundamental mode RR Lyræ, derived for IRAC 3.6 and 4.5 $\mu$m bands. The period slope we measured in the NIR 2015ApJ…799..165B and MIR 2017ApJ…XXX…..N are shown in Figure 1, and are in good agreement with their theoretical values. Repeating this measurement for the other single-metallicity clusters observed as part of the CRRP program will allow us to refine the fitted values of $b_{\lambda}$, and check that this term is indeed independent from [Fe/H]. Once the $b_{\lambda}$ are determined, we will use our sample of field calibrator RR Lyræ to fit the remaining parameters $a_{\lambda}$ and $c_{\lambda}$.

The left panel of Figure 2 assesses how well, if the PLZ is fully calibrated, it can be used to accurately derive the distance and average reddening of GCs. For this purpose we have used the photometric data of all RR Lyræ in M4, at $RIJHK$ and IRAC MIR bands, fitting at each wavelength the cluster reddened distance modulus by assuming the validity of the theoretical PLZ from 2015ApJ…808…50M . We have then derived the average distance modulus and extinction, adopting the peculiar reddening law for M4 derived by 2012AJ….144…25H . The result is consistent with independent measurements of the cluster’s distance and average $A_{V}$, and shows that this method is capable of providing accurate distances for individual Galactic GCs with a total uncertainty of just a few percent.

4 High Precision Individual RR Lyræ Distances

The theoretical PLZ relations derived by 2015ApJ…808…50M can be used to estimate the effect of metallicity on individual RR Lyræ distances. This is shown in Figure 3, where the distance modulus and extinction of each calibrator in our sample was fit using the theoretical PLZ relations (one example in the left panel, with the best fit confidence level derived using the Feldman and Cousin method 1998PhRvD..57.3873F , appropriate when a fitting parameter is bound by a physical boundary, e.g. $A_{V}\geq 0$). The central panels show the distribution of the calibrators’ magnitudes (derived using their best fit distance and extinction) in the $I$ and IRAC 3.6 $\mu$m bands, binned based on their [Fe/H]. In both wavelengths the dispersion is of the order of $\pm 0.13$ magnitudes. Once the metallicity term $c_{\lambda}\cdot\textrm{[Fe/H]}$ is subtracted from the stars’ absolute magnitude, the spread is significantly reduced (right). This is especially true in the MIR bands, where temperature and evolutionary effects are vanishing, leading to drastically smaller intrinsic scatter in the resulting PLZ (once metallicity effects are removed, we found $\sigma\simeq 0.02$ mag in the IRAC bands).

To test the validity of the theoretical PLZ derived by 2015ApJ…808…50M and 2017ApJ…XXX…..N we have compared our best fit distance moduli with the Gaia first data release (DR1, 2016A&A…595A…2G ) distances. Figure 4 shows that our results are within $1\sigma$ of the Gaia measurements, and that the residuals are not biased in either period, metallicity or the value of extinction derived with our best fitting procedure.

A detailed analysis of these results is provided in 2017ApJ…XXX…..N . Here we want to highlight how precisely calibrating the metallicity term in the RR Lyræ PLZ is essential to have reliable distances for individual stars, a necessary condition for using them as accurate tracers of old stellar populations. The typical $\sim 0.2$ dex uncertainty in most spectrophotometric determination of metallicity, just by itself, would result in a distance uncertainty of almost $\sim 4$ %, stressing the need of accurately measuring the individual star’s metallicity to obtain precise distances. Including the metallicity term is also important for adopting RR Lyræ as a new anchor for the distance scale: the reduced scatter in the MIR PLZ enables the determination of distances with accuracy comparable to Cepheids, fulfilling the promise of an independent route to the cosmological distance scale based on these Population II pulsators.

{acknowledgement}

Acknowledgments: This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the JPL-Caltech under a contract with NASA. This work also made use of data from the ESA mission Gaia, processed by the Gaia Data Processing and Analysis Consortium (DPAC).