Draft version August 28, 2020 Typeset using LATEX twocolumn style in AASTeX63 Near-infrared Census of RR Lyrae variables in the Messier 3 globular cluster and the Period–Luminosity Relations Anupam Bhardwaj,1, ∗ Marina Rejkuba,2 Richard de Grijs,3,4,5 Gregory J. Herczeg,1 Harinder P. Singh,6 Shashi Kanbur,7 and Chow-Choong Ngeow8 1Kavli Institute for Astronomy and Astrophysics, Peking University, Yi He Yuan Lu 5, Hai Dian District, Beijing 100871, China 2European Southern Observatory, Karl-Schwarzschild-Straße 2, 85748, Garching, Germany 3Department of Physics and Astronomy, Macquarie University, Balaclava Road, Sydney, NSW 2109, Australia 4Research Centre for Astronomy, Astrophysics and Astrophotonics, Macquarie University, Balaclava Road, Sydney, NSW 2109, Australia 5International Space Science Institute – Beijing, 1 Nanertiao, Zhongguancun, Hai Dian District, Beijing 100190, China 6Department of Physics and Astrophysics, University of Delhi, Delhi-110007, India 7State University of New York, Oswego, NY 13126, USA 8Graduate Institute of Astronomy, National Central University, 300 Jhongda Road, 32001 Jhongli, Taiwan (Received 31 July 2020; Revised August 28, 2020; Accepted August 28, 2020) Submitted to AJ ABSTRACT We present new near-infrared (JHKs) time-series observations of RR Lyrae variables in the Messier 3 (NGC 5272) globular cluster using the WIRCam instrument at the 3.6-m Canada France Hawaii Telescope. Our observations cover a sky area of ∼ 21′ × 21′ around the cluster center and provide an average of twenty epochs of homogeneous JHKs-band photometry. New homogeneous photometry is used to estimate robust mean magnitudes for 175 fundamental-mode (RRab), 47 overtone-mode (RRc), and 11 mixed-mode (RRd) variables. Our sample of 233 RR Lyrae variables is the largest thus far obtained in a single cluster with time-resolved, multi-band near-infrared photometry. Near-infrared to optical amplitude ratios for RR Lyrae in Messier 3 exhibit a systematic increase moving from RRc to short-period (P < 0.6 days) and long-period (P & 0.6 days) RRab variables. We derive JHKs-band Period–Luminosity relations for RRab, RRc, and the combined sample of variables. Absolute calibrations based on the theoretically predicted Period–Luminosity–Metallicity relations for RR Lyrae stars yield a distance modulus, µ = 15.041 ± 0.017 (statistical) ± 0.036 (systematic) mag, to Messier 3. When anchored to trigonometric parallaxes for nearby RR Lyrae stars from the Hubble Space Telescope and the Gaia mission, our distance estimates are consistent with those resulting from the theoretical calibrations, albeit with relatively larger systematic uncertainties. 1. INTRODUCTION metallicity relation with negligible dependence on pulsa- RR Lyrae (RRL) variables are low-mass (0.5 . tion periods unlike classical Cepheids (Bono et al.2003). arXiv:2008.11745v1 [astro-ph.SR] 26 Aug 2020 The reason for this different behavior is that the bolo- M/M ⊙ . 0.8), old (> 10 Gyr) stars that are located in a region between the cross-section of the horizontal branch metric correction’s sensitivity to effective temperature and the classical “instability strip” in the Hertzsprung– becomes significant only at longer wavelengths (R-band Russell diagram. These horizontal branch stars pulsate onwards, Catelan et al. 2004). Indeed, RRL exhibit well during their central helium burning evolutionary phase, defined Period–Luminosity relations (PLRs) at infrared wavelengths, first demonstrated in pioneering work by similar to intermediate-mass (3 . M/M ⊙ . 10) classi- cal Cepheids. RRL follow a visual (V -band) magnitude– Longmore et al. (1986), which makes them excellent dis- tance indicators (see recent reviews, Beaton et al. 2018; Bhardwaj 2020). RRL play a key role in our understand- Corresponding author: Anupam Bhardwaj ing of stellar evolution and pulsation (Catelan 2009), [email protected]; [email protected] and as stellar population tracers for Galactic archaeol- ∗ IAU Gruber Foundation Fellow 2020 2 Bhardwaj A. et al. ogy and the cosmic distance scale (Kunder et al. 2018; models, Denissenkov et al. (2017) found a good agree- Beaton et al. 2018). ment with observed properties of RRL and non-variable Globular clusters (GCs) typically host a rich and ho- horizontal branch stars for a distance modulus and red- mogeneous population of RRL stars. Messier 3 (M3 dening of µ = 15.02 mag and E(B − V )=0.013 mag, or NGC 5272), located at a distance of ∼ 10 kpc, respectively. hosts one of the largest samples of RRL with a domi- RRL as distance indicators have gained signif- nant population of fundamental mode RRL (RRab) vari- icance with increasing NIR observations especially ables (Clement et al. 2001). M3 has a mean metallicity in GCs (Sollima et al. 2006; Coppola et al. 2011; of [Fe/H]∼ −1.5 dex (Harris 2010) and the observed Stetson et al. 2014; Braga et al. 2015; Navarrete et al. period distribution of its RRL population exhibits a 2015; Braga et al. 2018). These horizontal branch vari- sharp peak at fundamental pulsator period of 0.55 day ables can complement the tip of the red giant branch (Jurcsik et al. 2017), indicating that it is a typical stars to provide an absolute primary calibration for the Oosterhoff I type (OoI, Oosterhoff 1939; Fabrizio et al. population II distance ladder (Beaton et al. 2016). The 2019) cluster. While multiple stellar populations have homogeneous population of RRL in different GCs offers been detected along the red giant branch of M3 (e.g., the possibility to derive their PLRs and estimate the Massari et al. 2016; Lee & Sneden 2020), no signifi- dependence on metal abundance (Sollima et al. 2006). cant variation has been detected in the iron abundance While theoretical models predict a significant metallicity (Sneden et al. 2004). Furthermore, helium enhancement coefficient of the RRL Period–Luminosity–Metallicity (∆Y . 0.02) has been suggested to explain observed (PLZ) relation (e.g., Catelan et al. 2004; Marconi et al. properties of horizontal branch stars (Dalessandro et al. 2015), it is still a topic of active debate considering 2013; Valcarce et al. 2016; Denissenkov et al. 2017). the paucity of RRL with both high-resolution spectro- Insignificant interstellar reddening (VandenBerg et al. scopic metallicities and precise parallaxes suitable to es- 2016) and the close proximity of M3 motivated sev- tablish an empirical calibration (Muraveva et al. 2018; eral detailed long-term photometric investigations at Neeley et al. 2019; Bhardwaj 2020). Therefore, NIR ob- optical wavelengths (Bakos et al. 2000; Cacciari et al. servations of abundant RRL in M3 will not only be use- 2005; Benko 2006; Jurcsik et al. 2012, 2017, and ref- ful for studies of the distance scale but also complement erences therein). Optical photometry has been used optical and ultraviolet data for a rigorous comparison to investigate the Blazhko effect, multimode pulsations with evolutionary and pulsation models. and period doubling in M3 RRL variables (Jurcsik et al. In this work, we present NIR time-series observations 2015; Jurcsik 2019). The RRL population in M3 clus- of RRL in M3 for the largest sample of variables in an ter has also been explored at ultraviolet wavelengths individual GC. The paper is organized as follows. In Sec- (Siegel et al. 2015). At near-infrared (NIR) wave- tion 2, we describe the observations, the data reduction lengths, Longmore et al. (1990) derived Ks-band PLRs and the photometric calibrations. The NIR light curves for RRL in GCs including 49 variables in the outer re- and pulsation properties of the RRL are discussed in gion of M3. Apart from that, NIR photometry of RRL Section 3. We discuss the JHKs-band PLRs for M3 in M3 has been limited to a sample of 7 RRL in the RRL in Section 4 and estimate a robust distance to the inner region of the cluster (Butler 2003). cluster. The results are summarized in Section 5. M3 has been the subject of several theoretical stud- ies aimed at reproducing the observed pulsation prop- 2. OBSERVATIONS, DATA REDUCTION, AND erties and, in particular, the period distribution of PHOTOMETRIC CALIBRATION its RRL population (Catelan 2004; Castellani et al. 2.1. Observations and data reduction 2005; Fadeyev 2019). Catelan (2004) showed that Our NIR observations were obtained using the the predicted period distribution based on canon- WIRCam instrument (Puget et al. 2004) mounted on ical horizontal branch models is inconsistent with the 3.6-m Canada France Hawaii Telescope (CFHT) on observations, while Castellani et al. (2005) suggested the summit of Mauna Kea in Hawaii during 4 nights be- that a bimodal mass distribution would be required tween 26 and 29 May 2019. WIRcam is an array of four to reproduce the period distribution with canoni- 2048 × 2048 HgCdTe HAWAII-RG2 detectors arranged cal models. Marconi & Degl’Innocenti (2007) accu- in 2 × 2 grid with gaps of 45′′ between adjacent detec- rately modeled optical light curves of M3 RRL using tors. The pixel scale of each detector is 0.3′′ pixel−1 nonlinear pulsation models with [Fe/H]∼ −1.34 dex resulting in a field of view of ∼ 21′ × 21′. We requested (Carretta & Gratton 1997), and estimated a distance JHKs time-series observations in queue mode centered modulus of 15.10 ± 0.10 mag. Using horizontal branch on the M3 cluster center, and obtained 22 epochs in J JHKs observations of RR Lyrae stars in M3 3 Table 1. Log of NIR observations. J-band H-band Ks-band Date MJD Airmass IQ Nf MJD Airmass IQ Nf MJD Airmass IQ Nf ET days arcsec days arcsec days arcsec sec 2019-05-26 58629.2882