Draft version January 23, 2019 Preprint typeset using LATEX style emulateapj v. 01/23/15 PHYSICAL PROPERTIES OF GALACTIC RV TAURI STARS FROM GAIA DR2 DATA A. Bodi´ 1,2 and L.L. Kiss1,3 Draft version January 23, 2019 ABSTRACT We present the first period-luminosity and period-radius relation of Galactic RV Tauri variable stars. We have surveyed the literature for all variable stars belonging to this class and compiled the full set of their photometric and spectroscopic measurements. We cross-matched the final list of stars with the Gaia DR2 database and took the parallaxes, G-band magnitudes and effective temperatures to calculate the distances, luminosities and radii using a probabilistic approach. As it turned out, the sample was very contaminated and thus we restricted our study to those objects for which the RV Tau-nature was securely confirmed. We found that several stars are located outside the red edge of the classical instability strip, which implies a wider pulsational region for RV Tau stars. The period-luminosity relation of galactic RV Tauri stars is steeper than that of the shorter-period Type II Cepheids, in agreement with previous result obtained for the Magellanic Clouds and globular clusters. The median masses of RVa and RVb stars were calculated to be 0.45-0.52 M and 0.83 M , respectively. Keywords: stars: binaries { stars: oscillations { stars: variables: Cepheids { stars: AGB and post-AGB { stars: fundamental parameters 1. INTRODUCTION photometric sub-classes, respectively. RV Tauri-type variables form the long-period exten- Previous studies on the period-luminosity (PL) rela- sion of Population II Cepheids, which are metal-poor tions of RV Tauri stars were almost exclusively based on low- and intermediate mass F-, G- and K-type supergiant various samples of Population II Cepheids in the Magel- stars, older than the classical Cepheids (see Wallerstein lanic Clouds or globular clusters. There were hints of a (2002) for a review). RV Tauris are among the most lu- different slope for longer-period Type II Cepheids (Mc- Namara 1995), which was found to depend on the wave- minous stars (103 - 104 L ), having already left the red length of the observations, with negligible effects in the giant (RGB) or the asymptotic giants branch (AGB) and JHK bands (Matsunaga et al. 2006). Recently, Groe- rapidly evolving through the post-RGB/AGB instability S newegen & Jurkovi´c (2017a) and Manick et al. (2017) strip to become planetary nebulae (Manick et al. 2018; presented supporting indications of a steeper RV Tau PL Kamath, Wood & Van Winckel 2015, 2014; Jura 1986). relation for the dusty objects in the Magellanic Clouds Hence, they provide important information about the not and the Milky Way, respectively. so well-known late-phases of stellar evolution, in which Until now no PL relation has been published for pulsations and mass-loss processes can interact and in- nearby, bright and in all other aspects well-observed fluence stellar evolution. Galactic RV Tauri stars. Gaia DR2 has opened, for the During their evolution, post-AGB stars cross the in- first time, the possibility of a geometric distance mea- stability strip where they become unstable against ra- surement of Galactic RV Tau-stars. The main inspira- dial pulsation. The observed periods of RV Tau stars tion of our study is to compare the RV Tau populations usually fall between 20 and 90 days (Soszy´nskiet al. of the Milky Way and the Magellanic Clouds and the 2010, 2008). The main characteristics of the light curve universality of their PL relations. is the presence of alternating minima (i.e. every second minima are shallower). The periodicity is not strict, as 2. DATA AND METHODS the cycle-to-cycle variations can be quite significant, and To identify all the known stars in the Galaxy that are in some cases, it has been shown to be caused by low- arXiv:1901.01409v2 [astro-ph.SR] 22 Jan 2019 thought to be of RV Tau-type, we searched the catalog dimensional chaos (e.g. Buchler et al. 1996, Plachy, B´odi of variable star index (VSX database4), the General Cat- & Koll´ath 2018). In addition to the Cepheid-like pulsa- alog of Variable Stars (GCVS; Samus et al. 2017) and tions, some RV Tau stars show long-term modulation of the SIMBAD (Wenger et al. 2000) database. Then, we the mean brightness with periods of 700-2500 days, asso- cross-correlated this sample with the Gaia archive5 (Gaia ciated with time-variable dust obscuration (Kiss & B´odi Collaboration et al. 2016, 2018) and downloaded all 2017). The absence or presence of the slow modulation is available measurements for stars that have relative par- the basis for classifying the stars into the RVa and RVb allax error (σπ/π) smaller than 0.2 (well-determined par- [email protected] allaxes; Astraatmadja & Bailer-Jones 2016). We have 1 Konkoly Observatory, Research Centre for Astronomy and found 56 stars. We note that in most cases extinction Earth Sciences, Hungarian Academy of Sciences, H-1121 Bu- dapest, Konkoly Thege M. ´ut15-17, Hungary values are missing in the Gaia table. This turned out to 2 MTA CSFK Lend¨ulet Near-Field Cosmology Research be a consequence of the data filtering done by the Gaia Group 3 Sydney Institute for Astronomy, School of Physics A29, Uni- 4 http://www.aavso.org/vsx/ versity of Sydney, NSW 2006, Australia 5 https://gea.esac.esa.int/archive/ 2 Bodi´ & Kiss team (see Sect. 6.5, Eqs. (8-11) and Figs. 31 and 19d in Andrae et al. 2018), which practically removed all stars 9000 20 between the upper main sequence and the RGB/AGB in 10 the Hertzsprung{Russell-diagram. Because of this and 8000 the strong degeneracy between extinction and tempera- tures, the Gaia T values were also found to be system- 0 eff 7000 0 5000 lit Gaia atically biased (see details below). Teff −Teff Bailer-Jones (2015) showed that distance estimation Gaia eff 6000 from parallaxes becomes an inference problem when mea- T surement errors are present. Traditional inverse ap- proach gives an incorrect (symmetric) error estimate, 5000 which can be avoided by using a properly normalized prior. Astraatmadja & Bailer-Jones (2016) investigated 4000 the performance of various priors for estimating distances and found that the exponentially decreasing space den- 3000 sity (EDSD) prior performs well with a length scale of 3000 4000 5000 6000 7000 8000 9000 10000 lit 1.35 kpc. To determine Gaia distances we followed the Teff prescription of the EDSD method. To calculate the absolute magnitudes and luminosities, Figure 1. A comparison of effective temperatures obtained from we used photometric measurements, extinction and bolo- Gaia and those derived using spectroscopy. The black dashed line metric correction values. 2MASS J,H,Ks and Johnson shows the 1:1 relation. The increasing deviation is clearly visible as the literature Teff increases. The insert shows the distribution V,I band photometric values were taken from the SIM- of the differences. BAD catalog. However, during the calculation of lumi- nosities, the Gaia magnitudes were preferred, if avail- When we plotted the Hertzsprung{Russell-diagram able. Extinctions AV were taken from the combination and the period-luminosity relation, we found that the of 3D reddening maps by Marshall et al. (2006), Green observational scatter was huge. The initial sample con- et al. (2015), and Drimmel, Cabrera-Lavers & L´opez- tained several low-luminosity stars, and also objects lo- Corredoira (2003) as implemented in the python package cated far from the theoretical instability strip. We in- of mwdust (Bovy et al. 2016). Bolometric corrections terpreted the large scatter as due to contamination of (BC) were linearly interpolated from Teff , logg, [Fe/H] misclassified objects, hence we performed a strict revi- and AV in the MIST/C3K grid (C. Conroy et al., in sion of the sample as follows. 6 preparation ). Teff , logg and [Fe/H] values were taken First, we thoroughly surveyed the literature for all from the literature or the Gaia archive. Missing logg and stars in the initial sample to remove all objects for which [Fe/H] values were filled by 0.0, assuming the errors to be there was the slightest doubt about the RV Tau-nature. 0.4 and 0.5, respectively. These parameters have subtle Second, we took into account the systematic investiga- influence on the BCs. tion of misclassified RV Tauris by Zsoldos (1991). Third, Absolute magnitude, luminosity and radius values were we also excluded stars with poorly determined pulsation determined in a probabilistic approach using the direct periods (meaning that we may have left out genuine RV mode of the slightly modified isoclassify code of Hu- Tauris, too, which need further photometric observations ber et al. (2017), which uses a Monte-Carlo sampling to measure the periods). scheme and derives posterior distributions of all param- In the next step, we reviewed the AAVSO and ASAS eters. The parameter estimation was performed star by light curves for each star remaining and checked the vari- star as follows. First, the distance is determined, then ability characteristics by visual inspection. Stars with the AV is estimated from the reddening map, which is few observations and those that do not show the alter- used to calculate the absolute magnitude in an available nating behavior clearly (i.e. every minima have almost photometric band. Using BC, luminosity is determined, the same depth) were filtered out. During the inspec- which is then converted into radius using the effective tion, we found incorrectly determined magnitudes in the temperature and the Stefan-Boltzmann law. The final commonly used catalogues.