Astronomy & Astrophysics manuscript no. GaiaCons_EyerEtAl c ESO 2020 April 17, 2020

Gaia Data Release 2 Variable stars in the colour-absolute magnitude diagram

Gaia Collaboration, L. Eyer1, L. Rimoldini2, M. Audard1, R.I. Anderson3, 1, K. Nienartowicz2, F. Glass1, O. Marchal4, M. Grenon1, N. Mowlavi1, B. Holl1, G. Clementini5, C. Aerts6, 7, T. Mazeh8, D.W. Evans9, L. Szabados10, A.G.A. Brown11, A. Vallenari12, T. Prusti13, J.H.J. de Bruijne13, C. Babusiaux4, 14, C.A.L. Bailer-Jones15, M. Biermann16, F. Jansen17, C. Jordi18, S.A. Klioner19, U. Lammers20, L. Lindegren21, X. Luri18, F. Mignard22, C. Panem23, D. Pourbaix24, 25, S. Randich26, P. Sartoretti4, H.I. Siddiqui27, C. Soubiran28, F. van Leeuwen9, N.A. Walton9, F. Arenou4, U. Bastian16, M. Cropper29, R. Drimmel30, D. Katz4, M.G. Lattanzi30, J. Bakker20, C. Cacciari5, J. Castañeda18, L. Chaoul23, N. Cheek31, F. De Angeli9, C. Fabricius18, R. Guerra20, E. Masana18, R. Messineo32, P. Panuzzo4, J. Portell18, M. Riello9, G.M. Seabroke29, P. Tanga22, F. Thévenin22, G. Gracia-Abril33, 16, G. Comoretto27, M. Garcia-Reinaldos20, D. Teyssier27, M. Altmann16, 34, R. Andrae15, I. Bellas-Velidis35, K. Benson29, J. Berthier36, R. Blomme37, P. Burgess9, G. Busso9, B. Carry22, 36, A. Cellino30, M. Clotet18, O. Creevey22, M. Davidson38, J. De Ridder6, L. Delchambre39, A. Dell’Oro26, C. Ducourant28, J. Fernández-Hernández40, M. Fouesneau15, Y. Frémat37, L. Galluccio22, M. García-Torres41, J. González-Núñez31, 42, J.J. González-Vidal18, E. Gosset39, 25, L.P. Guy2, 43, J.-L. Halbwachs44, N.C. Hambly38, D.L. Harrison9, 45, J. Hernández20, D. Hestroffer36, S.T. Hodgkin9, A. Hutton46, G. Jasniewicz47, A. Jean-Antoine-Piccolo23, S. Jordan16, A.J. Korn48, A. Krone-Martins49, A.C. Lanzafame50, 51, T. Lebzelter52, W. Löffler16, M. Manteiga53, 54, P.M. Marrese55, 56, J.M. Martín-Fleitas46, A. Moitinho49, A. Mora46, K. Muinonen57, 58, J. Osinde59, E. Pancino26, 56, T. Pauwels37, J.-M. Petit60, A. Recio-Blanco22, P.J. Richards61, A.C. Robin60, L.M. Sarro62, C. Siopis24, M. Smith29, A. Sozzetti30, M. Süveges15, J. Torra18, W. van Reeven46, U. Abbas30, A. Abreu Aramburu63, S. Accart64, G. Altavilla55, 56, 5, M.A. Álvarez53, R. Alvarez20, J. Alves52, A.H. Andrei65, 66, 34, E. Anglada Varela40, E. Antiche18, T. Antoja13, 18, B. Arcay53, T.L. Astraatmadja15, 67, N. Bach46, S.G. Baker29, L. Balaguer-Núñez18, P. Balm27, C. Barache34, C. Barata49, D. Barbato68, 30, F. Barblan1, P.S. Barklem48, D. Barrado69, M. Barros49, M.A. Barstow70, S. Bartholomé Muñoz18, J.-L. Bassilana64, U. Becciani51, M. Bellazzini5, A. Berihuete71, S. Bertone30, 34, 72, L. Bianchi73, O. Bienaymé44, S. Blanco-Cuaresma1, 28, 74, T. Boch44, C. Boeche12, A. Bombrun75, R. Borrachero18, D. Bossini12, S. Bouquillon34, G. Bourda28, A. Bragaglia5, L. Bramante32, M.A. Breddels76, A. Bressan77, N. Brouillet28, T. Brüsemeister16, E. Brugaletta51, B. Bucciarelli30, A. Burlacu23, D. Busonero30, A.G. Butkevich19, R. Buzzi30, E. Caffau4, R. Cancelliere78, G. Cannizzaro79, 7, T. Cantat-Gaudin12, 18, R. Carballo80, T. Carlucci34, J.M. Carrasco18, L. Casamiquela18, M. Castellani55, A. Castro-Ginard18, P. Charlot28, L. Chemin81, A. Chiavassa22, G. Cocozza5, G. Costigan11, S. Cowell9, F. Crifo4, M. Crosta30, C. Crowley75, J. Cuypers†37, C. Dafonte53, Y. Damerdji39, 82, A. Dapergolas35, P. David36, M. David83, P. de Laverny22, F. De Luise84, R. De March32, D. de Martino85, R. de Souza86, A. de Torres75, J. Debosscher6, E. del Pozo46, M. Delbo22, A. Delgado9, H.E. Delgado62, S. Diakite60, C. Diener9, E. Distefano51, C. Dolding29, P. Drazinos87, J. Durán59, B. Edvardsson48, H. Enke88, K. Eriksson48, P. Esquej89, G. Eynard Bontemps23, C. Fabre90, M. Fabrizio55, 56, S. Faigler8, A.J. Falcão91, M. Farràs Casas18, L. Federici5, G. Fedorets57, P. Fernique44, F. Figueras18, F. Filippi32, K. Findeisen4, A. Fonti32, E. Fraile89, M. Fraser9, 92, B. Frézouls23, M. Gai30, S. Galleti5, D. Garabato53, F. García-Sedano62, A. Garofalo93, 5, N. Garralda18, A. Gavel48, P. Gavras4, 35, 87, J. Gerssen88, R. Geyer19, P. Giacobbe30, G. Gilmore9, S. Girona94, G. Giuffrida56, 55, M. Gomes49, M. Granvik57, 95, A. Gueguen4, 96, A. Guerrier64, J. Guiraud23, R. Gutiérrez-Sánchez27, 4 87, 35 16, 15 4 48 76 4 19 arXiv:1804.09382v2 [astro-ph.SR] 16 Apr 2020 R. Haigron , D. Hatzidimitriou , M. Hauser , M. Haywood , U. Heiter , A. Helmi , J. Heu , T. Hilger , D. Hobbs21, W. Hofmann16, G. Holland9, H.E. Huckle29, A. Hypki11, 97, V. Icardi32, K. Janßen88, G. Jevardat de Fombelle2, P.G. Jonker79, 7, Á.L. Juhász10, 98, F. Julbe18, A. Karampelas87, 99, A. Kewley9, J. Klar88, A. Kochoska100, 101, R. Kohley20, K. Kolenberg102, 6, 74, M. Kontizas87, E. Kontizas35, S.E. Koposov9, 103, G. Kordopatis22, Z. Kostrzewa-Rutkowska79, 7, P. Koubsky104, S. Lambert34, A.F. Lanza51, Y. Lasne64, J.-B. Lavigne64, Y. Le Fustec105, C. Le Poncin-Lafitte34, Y. Lebreton4, 106, S. Leccia85, N. Leclerc4, I. Lecoeur-Taibi2, H. Lenhardt16, F. Leroux64, S. Liao30, 107, 108, E. Licata73, H.E.P. Lindstrøm109, 110, T.A. Lister111, E. Livanou87, A. Lobel37, M. López69, D. Lorenz52, S. Managau64, R.G. Mann38, G. Mantelet16, J.M. Marchant112, M. Marconi85, S. Marinoni55, 56, G. Marschalkó10, 113, D.J. Marshall114, M. Martino32, G. Marton10, N. Mary64, D. Massari76, G. Matijevicˇ88, P.J. McMillan21, S. Messina51, D. Michalik21, N.R. Millar9, D. Molina18, R. Molinaro85, L. Molnár10,

Article number, page 1 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

P. Montegriffo5, R. Mor18, R. Morbidelli30, T. Morel39, S. Morgenthaler115, D. Morris38, A.F. Mulone32, T. Muraveva5, I. Musella85, G. Nelemans7, 6, L. Nicastro5, L. Noval64, W. O’Mullane20, 43, C. Ordénovic22, D. Ordóñez-Blanco2, P. Osborne9, C. Pagani70, I. Pagano51, F. Pailler23, H. Palacin64, L. Palaversa9, 1, A. Panahi8, M. Pawlak116, 117, A.M. Piersimoni84, F.-X. Pineau44, E. Plachy10, G. Plum4, E. Poggio68, 30, E. Poujoulet118, A. Prša101, L. Pulone55, E. Racero31, S. Ragaini5, N. Rambaux36, M. Ramos-Lerate119, S. Regibo6, C. Reylé60, F. Riclet23, V. Ripepi85, A. Riva30, A. Rivard64, G. Rixon9, T. Roegiers120, M. Roelens1, M. Romero-Gómez18, N. Rowell38, F. Royer4, L. Ruiz-Dern4, G. Sadowski24, T. Sagristà Sellés16, J. Sahlmann20, 121, J. Salgado122, E. Salguero40, N. Sanna26, T. Santana-Ros97, M. Sarasso30, H. Savietto123, M. Schultheis22, E. Sciacca51, M. Segol124, J.C. Segovia31, D. Ségransan1, I-C. Shih4, L. Siltala57, 125, A.F. Silva49, R.L. Smart30, K.W. Smith15, E. Solano69, 126, F. Solitro32, R. Sordo12, S. Soria Nieto18, J. Souchay34, A. Spagna30, F. Spoto22, 36, U. Stampa16, I.A. Steele112, H. Steidelmüller19, C.A. Stephenson27, H. Stoev127, F.F. Suess9, J. Surdej39, E. Szegedi-Elek10, D. Tapiador128, 129, F. Taris34, G. Tauran64, M.B. Taylor130, R. Teixeira86, D. Terrett61, P. Teyssandier34, W. Thuillot36, A. Titarenko22, F. Torra Clotet131, C. Turon4, A. Ulla132, E. Utrilla46, S. Uzzi32, M. Vaillant64, G. Valentini84, V. Valette23, A. van Elteren11, E. Van Hemelryck37, M. van Leeuwen9, M. Vaschetto32, A. Vecchiato30, J. Veljanoski76, Y. Viala4, D. Vicente94, S. Vogt120, C. von Essen133, H. Voss18, V. Votruba104, S. Voutsinas38, G. Walmsley23, M. Weiler18, O. Wertz134, T. Wevers9, 7, Ł. Wyrzykowski9, 116, A. Yoldas9, M. Žerjal100, 135, H. Ziaeepour60, J. Zorec136, S. Zschocke19, S. Zucker137, C. Zurbach47, and T. Zwitter100 (Affiliations can be found after the references)

Received ; accepted

ABSTRACT

Context. The ESA Gaia mission provides a unique time-domain survey for more than 1.6 billion sources with G . 21 mag. Aims. We showcase stellar variability in the Galactic colour-absolute magnitude diagram (CaMD). We focus on pulsating, eruptive, and cataclysmic variables, as well as on stars that exhibit variability that is due to rotation and eclipses. Methods. We describe the locations of variable star classes, variable object fractions, and typical variability amplitudes throughout the CaMD and show how variability-related changes in colour and brightness induce ‘motions’. To do this, we use 22 months of calibrated photometric, spectro-photometric, and astrometric Gaia data of stars with a significant parallax. To ensure that a large variety of variable star classes populate the CaMD, we crossmatched Gaia sources with known variable stars. We also used the statistics and variability detection modules of the Gaia variability pipeline. Corrections for interstellar extinction are not implemented in this article. Results. Gaia enables the first investigation of Galactic variable star populations in the CaMD on a similar, if not larger, scale as was previously done in the Magellanic Clouds. Although the observed colours are not corrected for reddening, distinct regions are visible in which variable stars occur. We determine variable star fractions to within the current detection thresholds of Gaia. Finally, we report the most complete description of variability-induced motion within the CaMD to date. Conclusions. Gaia enables novel insights into variability phenomena for an unprecedented number of stars, which will benefit the understanding of stellar astrophysics. The CaMD of Galactic variable stars provides crucial information on physical origins of variability in a way that has previously only been accessible for Galactic star clusters or external galaxies. Future Gaia data releases will enable significant improvements over this preview by providing longer time series, more accurate astrometry, and additional data types (time series BP and RP spectra, RVS spectra, and radial velocities), all for much larger samples of stars. Key words. stars: general – Stars: variables: general – Stars: oscillations – binaries: eclipsing – Surveys – Methods: data analysis

1. Introduction Proper and orbital motions of stars further enable mass measure- ments in multiple stellar systems, as well as the investigation of The ESA space mission Gaia (Gaia Collaboration et al. 2016b) cluster membership. has been conducting a unique survey since the beginning of its operations (end of July 2014). Its uniqueness derives from sev- Thirdly, Gaia data are homogeneous throughout the entire eral aspects that we list in the following paragraphs. sky, since they are observed with a single set of instruments and Firstly, Gaia delivers nearly simultaneous measurements in are not subject to the Earth’s atmosphere or seasons. All-sky sur- the three observational domains on which most stellar astronom- veys cannot be achieved using a single ground-based telescope; ical studies are based: astrometry, photometry, and spectroscopy surveys using multiple sites and telescopes and instruments re- (Gaia Collaboration et al. 2016a; van Leeuwen et al. 2017). As quire cross-calibration, which unavoidably introduces systemat- consequence of the spin of the spacecraft, it takes about 80 sec- ics and reduces precision because of the increased scatter. Thus, onds for sources to be measured from the first to the last CCD Gaia will play an important role as a standard source in cross- during a single field-of-view transit. calibrating heterogeneous surveys and instruments, much like Secondly, the Gaia data releases provide accurate astromet- the Hipparcos mission (Perryman et al. 1997; ESA 1997) did in ric measurements for an unprecedented number of objects. In the past. Of course, Gaia represents a quantum leap from Hip- particular, trigonometric parallaxes carry invaluable information, parcos in many regards, including an increase of four orders of since they are required to infer stellar luminosities, which form magnitude in the number of objects observed, additional types of the basis of the understanding of much of stellar astrophysics. observations (spectrophotometry and spectroscopy), and signif-

Article number, page 2 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars icantly improved sensitivity and precision for all types of mea- (e.g. Hermes et al. 2017). The results we provide in Sects. 3 and 4 surements. on the variability fractions and levels are representative of milli- Fourthly, there are unprecedented synergies for calibrating mag level variability and not of micro-mag levels, as are found distance scales using the dual astrometric and time-domain ca- in space asteroseismic data. pabilities of Gaia (e.g. Eyer et al. 2012). Specifically, Gaia will Any of these asteroseismic surveys can benefit from Gaia enable the discovery of unrivalled numbers of standard candles astrometry, however, so that distances and luminosities can be residing in the Milky Way, and anchor Leavitt laws (period- derived, as De Ridder et al. (2016) and Huber et al. (2017) re- luminosity relations) to trigonometric parallaxes (see Gaia Col- ported with Gaia DR1 data. Gaia will also contribute to these laboration et al. 2017; Casertano et al. 2017, for two examples surveys with its photometry, and some surveys will also bene- based on the first Gaia data release). fit from the Gaia radial velocities (depending on their operating Variable stars have for a long time been recognised to offer magnitude range). crucial insights into stellar structure and evolution. Similarly, the Stellar variability comprises a great variety of observable Hertzsprung-Russell diagram (HRD) provides an overview of all features that are due to different physical origins. Figure 1 shows stages of stellar evolution, and together with its empirical cousin, the updated variability tree (Eyer & Mowlavi 2008), which pro- the colour-magnitude diagram (CMD), it has shaped stellar as- vides a useful overview of the various types of variability and trophysics like no other diagram. Henrietta Leavitt (1908) was their known causes. The variability tree has four levels: the dis- one of the first to note the immense potential of studying vari- tinction of intrinsic versus extrinsic variability, the separation able stars in populations, where distance uncertainties did not into major types of objects (asteroid, stars, and AGN), the phys- introduce significant scatter. Soon thereafter, Leavitt & Picker- ical origin of the variability, and the class name. In this article, ing (1912) discovered the period-luminosity relation of Cepheid we follow the classical distinction of the different causes of the variables, which has become a cornerstone of stellar physics and variability phenomena: variability induced by pulsation, rota- cosmology. It appears that Eggen (1951, his Fig. 42) was the first tion, eruption, eclipses, and cataclysmic events. A large number to use (photoelectric) observations of variable stars (in this case, of variability types can be identified in the Gaia data even now, classical Cepheids) to constrain regions where Cepheids occur in as described in the subsequent sections. the HRD; these regions are today referred to as instability strips. We here provide an overview of stellar variability in the Eggen further illustrated how Cepheids change their locus in the CaMD, building on the astrometric and photometric data of colour-absolute magnitude diagram (CaMD) during the course the second Gaia data release (DR2). Future Gaia DRs will en- of their variability, thus developing a time-dependent CMD for able much more detailed investigations of this kind using longer variable stars. Kholopov (1956) and Sandage (1958) later illus- temporal baselines, greater number of observations, and added trated the varying locations of variable stars in the HRD using classes of variable stars (such as eclipsing binaries, which will classical Cepheids located within star clusters. By combining the be published in DR3). different types of Gaia time-series data with Gaia parallaxes, we This paper is structured as follows. Section 2 shows the lo- are now in a position to construct time-dependent CaMD towards cation of different variability types in the CaMD, making use any direction in the Milky Way, building on previous work based of known objects from the literature that are published in Gaia on Hipparcos (Eyer et al. 1994; Eyer & Grenon 1997), but on a DR2, but without any further analysis of the Gaia data. Section 3 much larger scale. presents the fraction of variables as a function of colour and ab- Many variability (ground- and space-based) surveys have ex- solute magnitude, obtained by processing the Gaia time series ploited the power of identifying variable stars in stellar pop- for the detection of variability (Eyer et al. 2018). Section 4 inves- ulations at similar distances, for example, in star clusters or tigates the variability level in the CaMD by employing statistics nearby galaxies such as the Magellanic Clouds. Ground-based and classification results (some of which are related to unpub- microlensing surveys such as the Optical Gravitational Lensing lished Gaia time series). Section 5 shows the motion of known Experiment (OGLE; e.g. Udalski et al. 2015), the Expérience variables stars in the CaMD, that is, a time-dependent CaMD, pour la Recherche d’Objets Sombres (EROS Collaboration et al. which also includes sources that are not available in the DR2 1999), and the Massive Compact Halo Objects project (MA- archive but are online material. Section 6 summarises our results CHO; Alcock et al. 1993) deserve a special mention in this re- and presents an outlook to future Gaia DRs. Further informa- gard. The data will continue to grow with the next large multi- tion on the literature cross-match and on the selection criteria epoch surveys such as the Zwicky Transient Facility (Bellm applied to our data samples can be found in Appendices A and 2014) and the Large Synoptic Survey Telescope (LSST Science B, respectively. Collaboration et al. 2009) from the ground, and the Transit- ing Exoplanet Survey Satellite (TESS; Ricker et al. 2015) and 2. Location of variability types in the CaMD PLATO (Rauer et al. 2014) from space. Another ground-breaking observational trend has been the The precision of the location in the CaMD depends on the pre- long-term high-precision high-cadence uninterrupted space pho- cision of the colour on the one hand and on the determination tometry with CoRoT/BRITE (Auvergne et al. 2009; Pablo et al. of the absolute magnitude on the other. The precision of the 2016, with time bases of up to five months) and Kepler/K2 absolute magnitude of variable stars depends on the photomet- (Gilliland et al. 2010; Howell et al. 2014, with time bases of up to ric precision, the number of measurements, the amplitude of four years and three months, respectively) provided entirely new variability, and the relative parallax precision σ$/$. The up- insights into micro-magnitude level variability of stars, with pe- per limits of σ$/$ employed in this article vary between 5 and riodicities ranging from minutes to years. These missions opened 20%, which means that the uncertainty of the absolute magni- up stellar interiors from the detection of solar-like oscillations of tude that is solely due to the parallax uncertainty can be as large −1 thousands of sun-like stars and red giants (e.g. Bedding et al. as 5 (ln 10) σ$/$ ≈ 0.43 mag. 2011; Chaplin & Miglio 2013; Hekker & Christensen-Dalsgaard As we determined the colour as a function of integrated BP 2017, for reviews), as well as hundreds of intermediate-mass and integrated RP spectro-photometric measurements with tight stars (e.g. Aerts 2015; Bowman 2017) and compact pulsators constraints on the precision of these quantities (see Appendix B),

Article number, page 3 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

Variability Tree

Extrinsic Intrinsic Asteroids AGN

Rotation Stars Stars Eclipse

Microlensing Eclipse Rotation Secular Eclipse Eruptive Cataclysmic Pulsation PV Tel He star Photom. Period Asteroid RCB DY Per BLAP ZAND occultation V361 Hya Symbiotic UG (EC14026) Eclipsing Planetary short period sdB Dwarf novae L M binary transits FKCOM SN V1093 Her UV N (PG1716+426 / Betsy) Irregulars Miras Single red giants BY Dra Supernovae long period sdB Ceti Novae LPV PG 1159 SR Red dwarfs WR ELL (K-M stars) (DO,V GW Vir) ZZ Ceti SARV OSARG Semi- EA EW He/C/O-WDs regulars (DAV) H-WDs Small ampl. β Per / α Vir LBV ACYG V777 Her red var. EB RS CVn S Dor α Cygni (DBV) He-WDs SXA Binary red giants FU Solar-like GCAS RV PMS SX Arietis ACV Be stars SPBe RV Tau MS (B0-A7) with RR CW Credit: Eyer et al. (2018) 2 λ Eri strong B fields α Canum Venaticorum RR Lyrae CEP (W Vir / BL Her) MS (B8-A7) with Type II Ceph. Adapted from: Eyer & Mowlavi (2008) strong B fields Hot OB Supergiants PMS δ Cepheids δ Scuti BCEP GDOR roAp γ Doradus β Cephei SPB SXPHE DST Slowly δ Scuti pulsating B stars SX Phoenicis

Fig. 1. Updated version of the variability tree presented in Eyer & Mowlavi (2008), separated according to the cause of variability phenomena: variability induced by pulsations, rotation, eruptions, eclipses, and cataclysmic events. there are parts of the CaMD that are not explored here. For ex- multiplicity, rotation, inclination of the rotation axis, and chem- ample, the faint end of the main sequence presented in Fig. 9 ical composition. In this work, we do not correct for such phe- of Gaia Collaboration et al. (2018) does not fulfil the condition nomena and instead rely on the apparent magnitudes and colours on the precision in BP, so our diagrams do not include L and measured by Gaia, computing ‘absolute’ magnitudes using Gaia T brown dwarfs (which are fainter than MG ∼ 14 mag). If we parallaxes. We note that interstellar extinction and reddening can cross-match the Gaia data (conditional on the selection of Ap- be significant at the considered distances (up to 1 kpc), in partic- pendix B) with the catalogue of M dwarfs of Lépine & Gaidos ular for objects in the Galactic plane. This leads to distortions of (2011), only a few M6, M7, and M8 dwarf stars are found. certain observed features, such as the long tail in the red clump, In Fig. 2 we introduce the Gaia CaMD, which is displayed which extends to redder and fainter magnitudes. as a background in subsequent figures. For clarity, we note ba- The stellar variability aspects covered in the second Data Re- sic astronomical features such as the main sequence, the red lease of Gaia include a limited number of variability classes clump (and its long tail due to interstellar extinction), the hor- (Holl et al. 2018): long-period variables, Cepheids, RR Lyrae izontal branch, the extreme horizontal branch (see D’Cruz et al. stars, SX Phoenicis/δ Scuti stars, and rotation-modulated solar- 1996, for its physical origin), the red giant branch, the asymp- like variability (i.e. all late-type BY Draconis stars). Short- totic giant branch, the white dwarf sequence, the subdwarfs, the timescale variability (within one day) was explored regardless supergiants, and the binary sequence. There are additional sub- of the physical origin of the variability (Roelens et al. 2018), tle features above and below the red clump that are described in although stars that are classified as eclipsing binaries were re- Fig. 10 of Gaia Collaboration et al. (2018) and are known as the moved as planned to first appear in the third Data Release of asymptotic giant branch bump and the red giant branch bump, Gaia. The stars presented in this section are solely based on the respectively. On the right-hand side of Fig. 2, we also note the cross-match with known objects in the literature. The list of vari- typical limiting distance that can be reached because of the se- ability types presented here is not meant to be comprehensive. lection of σ$/$, up to 1 kpc (which was the largest distance we Figures 3–7 illustrate the locations of known variable stars considered for background stars). from catalogues in the literature that are cross-matched with the Several effects can influence the average location of a star in Gaia data. We indicate these targets according to their known the CaMD (in both axes), including interstellar extinction, stellar variability type published in the literature (the references are

Article number, page 4 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

Fig. 2. CaMD with its most striking known features (see text). The points in grey denote objects with parallax greater than 1 mas, with relative parallax precision better than 20% and other criteria described in Appendix B. listed in Table A.1), and only the stars that satisfy the selec- Future Gaia data releases will provide a more homogeneous tion criteria described in Appendix B are kept. Each of these variability classification that will rely primarily on the results figures includes as reference the location and density (in grey of the variability processing (Holl et al. 2018). (b) The CaMDs scale) of all stars, regardless of stellar variability, that satisfy the are not corrected for extinction, which leads to increased scat- astrometric and photometric criteria of Appendix B with the ad- ter, in particular for objects that primarily reside in heavily at- ditional constraint of a minimum parallax of 1 mas (i.e. within 1 tenuated areas such as the Galactic disc and the Galactic bulge. kpc to the Sun). This radius seems a good compromise between (c) The cross-match of sources can be erroneous when stars are a large number of stars and a limited effect of interstellar matter. located in crowded regions or have high proper motion, espe- Variable stars whose variability type was previously known are cially if the positions of stars in the published catalogues are not represented by combinations of symbols and colours. Following sufficiently precise or if proper motion information is not avail- the structure of the variability tree in Fig. 1, we show in sepa- able. (d) Some variability types like magnetically active stars rate figures the CaMDs of stars whose variability is induced by (e.g. RS CVn stars) exhibit different observational phenomena, different causes, such as pulsations, rotation, eruptions, eclipses, such as rotational modulation variations as well as flares. To and cataclysmic events. avoid overcrowding the CaMD diagrams, these types are repre- Several caveats apply to Figs. 3–7 and should be kept in mind sented in only one of the relevant diagrams (e.g. with rotational for their interpretation. (a) The quality of catalogues published or eruptive variables). Furthermore, we note that the time sam- in the literature can be rather different, in part because variability pling and the waveband coverage of a given survey might favour is often classified without knowledge of a parallax. To reduce the the detection of only some of these aspects. (e) Gaia represents impact of misclassified objects on these figures, we selected sub- a milestone for space astrometry and photometry. Nevertheless, sets of all available catalogues as reference for specific variable some sources can be affected by problems such as corrupt mea- star classes, depending on their agreement with the expected lo- surements so that their location in the CaMD may be incorrect cations in the CaMD. In certain cases, we have excluded sources (Arenou et al. 2018). However, we stress that these problems are from the literature by choice of specific catalogues (Table A.1) limited to a small fraction of sources so that most known vari- and by using the Gaia astrometry and the multi-band photomet- able classes are recovered as expected. The cyclic approach of ric time-series data for occasional cuts in magnitude or colour.

Article number, page 5 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl the Gaia data processing and analysis will allow us to correct 12. Rapidly oscillating Am and Ap stars, which are chemically for these unexpected features in the future data releases. peculiar A stars that exhibit multiperiodic non-radial pres- sure modes in the period range of about 5 − 20 min. 13. V361 Hydrae (or EC 14026) stars, which are subdwarf B 2.1. Pulsating variable stars stars on the extreme horizontal branch that pulsate in pres- Figure 3 shows the positions of different classes of pulsating sure modes with very short periods of ∼ 1 − 10 min. variable stars based on the Gaia data and can be compared to its 14. V1093 Her (or PG 1716) stars, which are subdwarf B stars on theoretical counterpart in recent textbooks on asteroseismology the extreme horizontal branch that pulsate in gravity modes (Fig. 1.12 in Aerts et al. 2010) and on pulsating stars (Catelan with periods of about 1 − 3 h. & Smith 2015). We refer to these books for further details of 15. ZZ Ceti stars, which are white dwarfs featuring fast non- specific variability classes. Here, we only consider the following radial gravity-mode pulsations with periods of 0.5 − 25 min. types of pulsating variable stars: The CaMD of pulsating stars carries a great deal of infor- 1. Long-period variables, which are red giant stars that populate mation, much of which has shaped the understanding of stellar the reddest and brightest regions of the CaMD. They include structure and evolution and can be found in textbooks. Briefly Miras, semi-regular variables, slow irregular variables, and summarised, we note the following particularly interesting fea- small-amplitude red giants. tures of Fig. 3. 2. α Cygni stars, which are luminous supergiant stars that pul- sate in non-radial modes. They are particularly affected by – Extinction affects variability classes that belong to differ- interstellar extinction as they are young massive stars that ent populations unequally, as expected. Stars located away reside in the Galactic disc, so that their position in Fig. 3 from the Galactic disc are much less reddened and thus must be treated with caution. clump more clearly. This effect is particularly obvious when 3. δ Scuti stars, which are Population-I stars of spectral types A RR Lyrae stars and classical Cepheids are compared, which and F with short periods (< 0.3 d) that dominantly pulsate in both occupy the same instability strip, and it cannot be ex- pressure modes, but may also reveal low-order gravity modes plained by the known fact that the classical instability strip of low amplitude. becomes wider in colour at higher luminosity (e.g. see An- 4. SX Phoenicis stars, which are Population-II high-amplitude derson et al. 2016; Marconi et al. 2005; Bono et al. 2000, and δ Scuti stars. references therein). 5. γ Doradus stars, which are stars with spectral types A and – Interstellar reddening blurs the boundaries between variabil- F with periods from 0.3 to 3 d that dominantly pulsate in ity classes. Correcting for interstellar extinction will be cru- high-order gravity modes, but may also reveal low-amplitude cial to delineate the borders of the instability strips in the pressure modes. CaMD, as well as to deduce their purity in terms of the frac- 6. RR Lyrae stars (Bailey’s type ab and c), which are tion of stars that exhibit pulsations while residing in such Population-II horizontal branch stars with periods from 0.2 regions. to 1 d that pulsate in pressure mode. C-type RR Lyrae stars – Practical difficulties involved in separating variable star are bluer than ab-type stars. 7. Slowly pulsating B (SPB) stars, which are non-radial multi- classes in the way required to construct Fig. 3 include a) that periodic gravity-mode pulsators of spectral type B with peri- variable stars are often subject to multiple types of variability ods typically in the range from 0.4 to 5 d. (e.g. γ Doradus/δ Scuti, β Cephei/SPB hybrid pulsators, pul- 8. β Cephei stars, which are late-O to early-B spectral type sating stars in eclipsing binary systems, or pulsating white non-supergiant stars with dominant low-order pressure and dwarfs that exhibit eruptions), and b) that naming conven- gravity modes that feature periods in the range from 0.1 to tions are often historical or purely based on light-curve mor- 0.6 d. Several of them have been found to also exhibit low- phology, so that they do not account for different evolution- amplitude high-order gravity modes as in the SPB stars (e.g. ary scenarios (e.g. type-II Cepheids). With additional data Stankov & Handler 2005). The β Cephei stars are located in and a fully homogeneous variable star classification based the Galactic disc so that their CaMD position is easily af- on Gaia alone, such ambiguities will be resolved in the fu- fected by interstellar extinction. ture unless they are intrinsically connected to the nature of 9. Classical Cepheids (prototype δ Cephei), which are evolved the variability. Population-I (young intermediate-mass) stars featuring ra- – We note multiple groups of ZZ Ceti stars along the white dial p-mode pulsations with periods of approximately 1 − dwarf sequence. The most prominent of these is located at 100 d. Cepheids can be strongly affected by interstellar ex- GBP − GRP ' 0 and MG ' 12, as reported in Fontaine & tinction as they reside in the Galactic disc and can be ob- Brassard (2008). served at great distances. 10. Type-II Cepheids, which are Population-II stars pulsating in p-mode that were historically thought to be identical to 2.2. Variability due to rotation and eclipses classical Cepheids. Type-II Cepheids consist of three differ- Figure 4 shows stars whose variability is induced by rotation. ent sub-classes (separated by period) that are commonly re- There are three primary categories: spotted stars, stars deformed ferred to as BL Herculis, W Virginis, and RV Tauri stars. by tidal interactions, and objects whose variability is due to light Their evolutionary scenarios differ significantly, although the reflected by a companion. Following the nomenclature in the lit- three sub-classes together define a tight period-luminosity re- erature (Table A.1), we list the following variability classes sep- lation. arately, although we note occasional overlaps among the defi- 11. PV Telescopii stars, which include the sub-classes V652 Her, nitions of these variability classes. The following types are in- V2076 Oph, and FQ Aqr (Jeffery 2008). They are rare cluded in Fig. 4: hydrogen-deficient supergiant stars that cover a wide range of spectral types and exhibit complex light and radial veloc- 1. α2 Canum Venaticorum stars, which are highly magnetic ity variations. variable Bp and Ap MS stars.

Article number, page 6 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

Fig. 3. Known pulsating variable stars retrieved from published catalogues are placed in the observational CaMD, with symbols and colours representing types as shown in the legend (see Appendix A.1 for the references from the literature per type). All stars satisfy the selection criteria described in Appendix B. The background points in grey denote a reference subset of objects with a stricter constraint on parallax ($ > 1 mas), which limits the sample size, extinction, and reddening. The effects of interstellar matter and other phenomena (see text) are not corrected for. The condition on the relative precision of GBP measurements introduces artificial cuts in the distributions of low-mass main-sequence stars and red (super)giants.

2. Spotted stars, which show rotational modulation variability – RS Canum Venaticorum stars are significantly brighter than from spots. BY Draconis stars near the bottom of the main sequence (at 3. BY Draconis stars, which are main-sequence stars with late cool temperatures). spectral types (K and M) that exhibit quasi-periodic light – The reflection binary class is primarily present among very curves due to spots and chromospheric activity. compact (subdwarf) stars; there is a cluster near absolute 4. RS Canum Venaticorum stars, which are spotted stars whose mag 4, GBP − GRP ∼ −0.4 mag. rotation-induced variability is frequently accompanied by – There seems to be a dearth of rotational spotted variables − ∼ other phenomena, such as eclipses and flares. around GBP GRP 0.4, which corresponds with the tran- sition region of stars with a radiative versus convective outer 5. Ellipsoidal variables, which show variability (without envelope. eclipses) due to orbital motion of a star that is distorted by a – SX Arietis stars form a fairly well-defined hot-temperature stellar companion. envelope of the most luminous alpha2 Canum Venaticorum 6. Solar-like stars with magnetic activity. Stars of this type in variables. Fig. 4 are limited to a catalogue focused on the Pleiades, which explains the thin distribution of the main sequence. Figure 5 shows eclipsing binary systems as well as stars We can see a hint of the binary sequence. that have been identified to host exoplanets through the transit 7. SX Arietis stars, which are similar to α2 Canum Venaticorum method. Symbols distinguish the following sub-classes: stars but have a higher temperature. We note that some over- 1. Eclipsing binaries of type EA; the prototype is Algol. Bina- lap of the two distributions occurs for these two variability ries with spherical or slightly ellipsoidal components with types. well-separated, nearly constant light curves in between min- 8. Binary systems with a strong reflection component in the ima. Secondary minima can be absent. light curve with re-radiation of the hotter star’s light from 2. Eclipsing binaries of type EB; the prototype is β Lyrae. Bina- the cooler companion’s surface. ries with continuously changing light curves and not clearly 9. FK Comae Berenices stars, which are spotted giant stars. defined onsets or ends of eclipses. Secondary minima are al- ways present, but can be significantly less deep than primary Figure 4 shows the following properties, among other things. minima.

Article number, page 7 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

3. Eclipsing binaries of type EW; the prototype is W Ursae Ma- 3. Z Andromedae stars, which are symbiotic binary stars com- joris. The components are nearly or actually in contact and posed of a giant and a white dwarf. minima are virtually equally strong. Onsets and ends of min- ima are not well defined. Further information on cataclysmic variables can be found for 4. Stars known to exhibit exoplanetary transits from the litera- example in Warner (2003) and Hellier (2001). ture. In Fig. 7, we note a clump of cataclysmic variables located in the ZZ Ceti variability strip near G ∼ 12 and GBP − GRP ∼ Based on Fig. 5, we observe that EA stars are present al- 0.1. The most significant clump of cataclysmic variables is near most throughout the CaMD. There are groups of EB stars that are G ∼ 4 and GBP − GRP ∼ 0.1 mag. These are probably binary overluminous compared to the white dwarf sequence, which are systems with stars from the extreme horizontal branch and the likely white dwarf stars with main-sequence companions. More- main sequence. over, the majority of the stars that host exoplanets are identified by the Kepler spacecraft, and only very few of them have de- tectable transits in the Gaia data because of the different photo- 3. Variable object fractions in the CaMD metric precision and time sampling. The different types of brightness variations as presented in the CaMD may strongly depend on the colour and absolute mag- 2.3. Eruptive and cataclysmic variables nitude, as seen in Sect. 2, because they are driven by differ- ent physical mechanisms. Similarly, the variable object fraction, Figure 6 focuses on eruptive variable stars. As for the rotation- which is defined as the number of variable objects per colour- ally induced variables, we adopt the nomenclature from the liter- magnitude bin divided by the total number of objects in the same ature (see Table A.1), which includes partially overlapping defi- bin, is expected to depend on the location in the CaMD. The vari- nitions. The following types are considered. able object fraction was previously determined based on variable 1. S Doradus stars, also known as luminous blue variables, are objects detected in the Hipparcos time series (ESA 1997), for massive evolved stars that feature major and irregular pho- example. Here we significantly expand this investigation using tometric variations that are due to heavy mass loss by a 13.5 million stars with heliocentric distances of up to 1 kpc that radiation-driven wind. satisfy the astrometric and photometric selection criteria listed 2. R Coronae Borealis stars, which are carbon-rich supergiants in Appendix B as well as (a) at least 20 observations in the G, that emit obscuring material and as a consequence have dras- GBP, and GRP bands, and (b) a relative parallax uncertainty of tic rapid dimming phases. < 5% . In order to reduce the number of objects that are affected 3. Wolf-Rayet (WR) stars, which are the almost naked helium by significant extinction, stars at low Galactic latitudes (from −5 core that is left over from originally very high mass evolved to 5 deg) are excluded. This effectively reduces the number of stars. They feature strong emission lines of hydrogen, nitro- disc variables such as classical Cepheids and β Cephei stars. gen, carbon, or oxygen. WR stars undergo very fast mass loss Fig. 8 illustrates this Gaia-based high-resolution map of the and can be significantly dust-attenuated. variable object fraction in the CaMD at the precision level of ap- 4. γ Cassiopeiae stars and stars with B spectral types that ex- proximately 5–10 mmag. Variability is identified in about 9% of hibit hydrogen emission lines, that is, Be stars. These are the stars based on a supervised classification of Gaia sources. emitting shell stars. During their ‘eruptive’ phenomena, they This method depends heavily on the selection of the training set become brighter. of constant and variable objects. Minor colour-coded features 5. Flare stars, which are magnetically active stars that display can be due to training-set related biases. The detection of vari- flares. This category incudes many subtypes of magnetically ability further depends upon the amplitude of the variables, their active stars, such as UV Ceti-type, RS CVn-type, and T Tauri apparent magnitude distribution, and the instrumental precision. stars. The accuracy of the fraction of variables is also affected by the 6. UV Ceti stars, which usually are K-M dwarfs that show number of sources per bin of absolute magnitude and colour, flares. which can be as low as one in the tails of the two-dimensional 7. T Tauri stars (classical and weak lined), which are young pre- source number density distribution. main sequence stars that either accrete strongly (classical) or Figure 8 contains many informative features, despite possi- show little sign of accretion (weak lined). These stars show ble biases. Future data releases will significantly improve upon variability that is due either to magnetic activity (e.g. rota- Fig. 8 by correcting for reddening and extinction and using a tional modulation and flares) or accretion (quasi-periodic, larger number of objects with more accurate source classifica- episodic, or stochastic variations), in addition to pulsations tions. For the time being, we remark the following. that may also occur in some of them. – The classical instability strip is clearly visible with a variabil- In Fig. 6, we notice the absence of eruptive variables among ity in about 50-60% of the stars (although extinction limits hot main-sequence stars (non-supergiants). This region is popu- the precision of this estimate). lated by pulsating stars, such as γ Doradus and δ Scuti stars, cf. – For evolved stars, red giants, and asymptotic giant branch Fig. 3. Moreover, WR stars, R Coronae Borealis stars, and S Do- stars, we find that higher luminosity and redder colour im- radus stars are among the most luminous stars in this diagram. plies a higher probability of variability. Figure 7 illustrates three types of cataclysmic variables. – The red clump has a very low fraction of variable stars in the 1. Cataclysmic variables (generic class), typically novae and Gaia data. Kepler photometry of red clump stars has revealed dwarf novae involving a white dwarf. Many of these stars complex variability at the micro-mag level that has been used are situated between the main and white dwarf sequences. extensively for asteroseismology, cf. Sect. 1 and references 2. U Geminorum stars, which are dwarf novae that in principle therein. consist of a white dwarf with a red dwarf companion that – The classical ZZ Ceti location is extremely concentrated in experiences mass transfer. colour and magnitude, with variability in about half of the

Article number, page 8 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

Fig. 4. Same as Fig. 3, but for rotational-induced variability types.

Fig. 5. Same as Fig. 3, but for eclipsing binaries (of types EA, EB, and EW) and known host-stars that show exoplanet transits. As expected, eclipsing binaries can be anywhere in the CaMD, which explains why they are the main source of contamination of pulsating stars, for instance.

Article number, page 9 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

Fig. 6. Same as Fig. 3, but for eruptive variability types.

Fig. 7. Same as Fig. 3, but for cataclysmic variables and some sub-types.

Article number, page 10 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

Fig. 8. Variable object fraction in the CaMD shown as a colour scale as labelled. This figure is not based on variable objects from the literature. Instead, variability is detected directly using Gaia data and employing supervised classification for sources with at least 20 observations in the G, GBP, and GRP bands. All objects satisfy the selection criteria described in Appendix B, but with more restrictive constraints on the parallax precision (parallax_over_error > 20) and on the parallax value ($ > 1 mas), which limits the sample (size, extinction, and reddening). In order to reduce the extinction effect, objects at low Galactic latitudes (from −5 to 5 deg) are excluded. About 9% of the 13.5 million stars that satisfy these criteria are variable. Some of the bins (especially the outlying ones) can contain only a few or even single sources. The condition on the relative precision of GBP measurements introduces artificial cuts in the distributions of low-mass main-sequence stars and red (super)giants.

stars. The concentration is due to the partial ionisation of points in Fig. 9, we show the (colour-coded) mean of the cor- hydrogen in the outer envelope of white dwarfs, which is rected G-band IQR of sources within each square bin measuring developed only in extremely narrow ranges of effective tem- 0.02 mag in both colour and magnitude after trimming the top peratures (see Fontaine & Brassard 2008). and bottom 5%. This binning was applied to each variability type – Extreme horizontal branch stars show a high probability of individually, and cuts were applied to select minimum classifica- variability. tion probabilities per type to minimise incorrect classifications. – The hottest and most luminous main-sequence stars are very We emphasise the location of variable object classes that feature frequently variable. large amplitudes by plotting classes with higher IQR on top of – There is a clear gradient towards larger fractions of variables variability classes with lower IQR. above the low-mass main-sequence stars. Figure 9 contains the following stellar variability types based on the all-sky classification (Rimoldini et al. 2018): α2 Canum Venaticorum, α Cygni, β Cephei, cataclysmic, clas- 4. Variability amplitudes in the CaMD sical Cepheids, δ Scuti, γ Cassiopeiae, γ Doradus, Mira, ellip- soidal, RR Lyrae of Bailey’s type ab and c, semiregular, slowly Figure 9 shows variability amplitudes as a function of position pulsating B stars, solar-like variability due to magnetic activ- in the CaMD. Here, we quantify variability amplitudes using the ity (flares, spots, and rotational modulation), SX Arietis, and G-band inter-quartile range (IQR). Objects are selected accord- SX Phoenicis. We did not include other classes (listed in Eyer ing to the general criteria described in Appendix B, with stricter et al. 2018) for clarity or because there were too few objects. conditions on the parallax (greater than 1 mas) and its relative We note that any specific selection criteria applied to the objects precision (better than 5%). To prevent the false impression that shown in Fig. 9 introduce biases that can highlight or diminish faint (and very bright) sources have intrinsically higher ampli- the prominence of certain phenomena. Nevertheless, Fig. 9 pro- tudes, we corrected for the instrumental spread of the IQR as a vides a first detailed illustration of some of the most important function of the median G magnitude. This correction was deter- amplitude-related variability features in the CaMD. A number mined using sources that are classified as constant in the all-sky of clumps and instability regions are visible in Fig. 9, which are classification (Rimoldini et al. 2018) and subtracted in quadra- related to the variability classes described in Sects. 2 and 3. We ture from the measured IQR. Instead of plotting individual data notice the following features.

Article number, page 11 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

Fig. 9. Amplitude of variability in the CaMD based on a selection of classified variables within 1 kpc and with a relative uncertainty for the parallax of 5%. The colour scale shows the corrected G-band IQR (see text) with a cut-off at 0.1 mag to emphasise the low- and mid-level variability. The background points in grey represent classified constant stars. All objects satisfy the selection criteria described in Appendix B, in addition to the stricter conditions on parallax and its precision, as mentioned above. The effects of interstellar extinction are not corrected for.

– The classical instability strip that contains classical Cepheids 5. Variability-induced motion in the CaMD and RR Lyrae stars is not very prominent, although some clumps (in red or cyan) are apparent. In this section, we visualise the variability-induced motion of – The instability regions linked to SPB stars and β Cephei stars stars in the time-dependent CaMD using all-sky measurements are broad and uniform. made in the G, GBP, and GRP passbands. Gaia data are uniquely – Higher amplitude variations are clearly correlated with red- suited to create this time-dependent CaMD, since the different der colours for long-period variables. data types (astrometric and photometric in three bands) are ac- – Variables with the highest amplitude (IQR > 0.1 mag) oc- quired in a quasi-simultaneous fashion at many epochs that are cur in several regions in the CaMD, including the classical distributed over a multi-year time span. The first of such repre- instability strip, long-period variables, below the red clump, sentations, although much less detailed, was presented for indi- above the main sequence of low-mass stars (in correspon- vidual classical Cepheids in the Milky Way (Eggen 1951) and in dence of the observed gradient in the fraction of variables), Galactic star clusters (Kholopov 1956; Sandage 1958). Similarly and between the white dwarf sequence and the main se- minded representations in the literature were based on data from quence. the SDSS (Ivezic´ et al. 2003, mostly Galactic objects), EROS – Significant amplitudes of > 0.04 mag are found very fre- (Spano et al. 2009, LMC objects), and, very recently, the HST quently among the coolest white dwarfs. – The stars between the main sequence and the white dwarfs observations of M51 (Conroy et al. 2018). sequence feature large variability amplitudes and extend into Figure 11 illustrates the variability-induced motion of stars the clump of ZZ Ceti stars in the white dwarf sequence. This in the CaMD. As elsewhere in this paper, no correction for inter- intermediate region is populated in particular by the high- stellar extinction is applied. Individual stars are shown by differ- amplitude cataclysmic variables, cf. Fig. 7. A closeup view ently (arbitrarily) coloured lines that connect successive absolute of the white dwarf sequence is shown in Fig. 10, which rep- G magnitudes and GBP − GRP measurements, that is, the obser- resents all classified variables within 200 pc. Each object is vations are ordered in time as opposed to variability phase. This plotted without binning to emphasise the variability of the choice was made to avoid uncertainties related to phase-folding ZZ Ceti stars. the relatively sparsely sampled light curves based on 22 or fewer

Article number, page 12 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

dominated by the effect of change in temperature rather than radius. For the longest-period variables, the 22-month time span of the Gaia data is similar to the pulsation cycle, so that in some cases, loop-like shapes are apparent. For variable stars with shorter periods (e.g. RR Lyrae stars or classical Cepheids), successive measurements in time are not gener- ally ordered in phase, so that an overall ‘envelope’ of many cycles is revealed. 2. The motions of eclipsing binary systems in the CaMD de- pends on the colour difference between the two stars. The components of eclipsing systems of the EW type have sim- ilar mass (and colour), leading to a rather vertically aligned motion (parallel to the absolute magnitude axis). For eclips- ing binary systems with stars of different mass (and thus colour) close to the main sequence, the deepest eclipse is usually slightly redder, since the secondary (less massive and redder) component eclipses part of the light of the pri- mary star. The slope of the motion of eclipsing binaries in Fig. 10. Same as Fig. 9, but focusing on the white dwarf sequence and the CaMD is very different (much steeper) than the one of plotting all classified variables within 200 pc with a relative uncertainty for the parallax better than 5%. A close inspection of this sequence re- pulsating stars. veals amplitudes at the level of 40 mmag in various regions. 3. Rotationally induced variables (shown here: α2 Canum Ve- naticorum stars) feature small amplitudes in absolute G and are rather horizontal in the CaMD. This is as expected from months of observations and to include both periodic and non- starspots, which have a lower temperature than the surround- periodic variable objects. ings and hence absorb the light at bluer frequencies to re- Figure 11 is limited to a subset of all available variable stars emit it by back-warming effect at redder frequencies. There- in order to avoid overcrowding the diagram. As a preview for fore, the magnitude change in a broad band like G is smaller future data releases, we include here the variability-induced mo- than it would be if measured in narrow bands. tions of some stars whose time series and variability types are not published in DR2 (but which are available as online material). 4. Eruptive stars (shown here: γ Cassiopeiae and Be-type stars) Figure 11 includes the following variability types as defined in become redder when brighter because of additional extinc- Sec. 2: α2 Canum Venaticorum variables, B-type emission line tion during their eruptive phase. The slopes of their motions /γ Cassiopeiae stars, cataclysmic variables, classical and type-II in the CaMD therefore have the opposite sign with respect to Cepheids, δ Scuti stars, eclipsing binaries, RR Lyrae stars, long- the sign of pulsating stars. period variables, and SX Phoenicis stars. All sources shown sat- 5. The variability of cataclysmic variables (shown here: novae) isfy the general criteria described in Appendix B and typically primarily features strong outbursts in the ultraviolet and blue 1 have at least ten available observations . We further prioritised part of the spectrum that are understood to be caused by mass the selection of objects featuring wider ranges of variations in transfer from donor stars in binary systems. These outbursts 2 the G band (with a minimum of about 0.1 mag) . The number of very significantly change the colour of the system towards sources shown for each variability type ranges from a few to sev- bluer values. eral tens and was selected to ensure clarity in case of high source density or overlapping variability classes in certain regions of the The current version of Fig. 11 represents a first step towards CaMD. In order to limit the effect of outlying values, time-series a more global description of stellar variability. The motions de- data are filtered by operators as described in Holl et al. (2018), scribed by the variable stars in the time-dependent CaMD pro- and the 10 % of the brightest and faintest observations in the vide new perspectives on the data that can be exploited as vari- G band are excluded for sources with G − G lower than RP BP RP able star classification attributes to appreciably improve the clas- 1.5 mag. Non-variable objects are shown as a grey background sification results. Gaia data will definitively help identify mis- to provide a visual reference for the variable object locations classifications and problems in published catalogues, thanks to in the CaMD. These stars satisfy the criteria described in Ap- its astrometry and the quasi-simultaneous measurements. pendix B as well as the stricter condition of $ > 1 mas. Stars whose variability is caused by different physical effects exhibit In future Gaia data releases, there will be more data points different motions within the time-dependent CaMD. We briefly per source, which will enable us to refine Fig. 11. In particular, summarise the different motions seen in Fig.11 as follows. periods can be determined with an accuracy inversely propor- tional to the total Gaia time base for periodic objects. In this way, 1. Pulsating stars, including long-period variables, Cepheids, the motion in the CaMD can be represented more precisely by RR Lyrae, and δ Scuti/SX Phoenicis stars, exhibit a similar connecting points that are sorted in phase (rather than in time). behaviour. These stars are bluer when brighter in G, which this leads to Lissajous-type configurations for pulsators. For suf- illustrates that brightness variations of pulsating stars are ficiently bright stars, radial velocity time series will add a third and unprecedented dimension to Fig. 11. 1 The minimum number of observations per source is increased to 20 in the case of long-period variables, but the condition on the number of An animated version of Fig. 11 is provided at https:// observations is removed for cataclysmic variables. www.cosmos.esa.int/web/gaia/gaiadr2_cu7. We provide 2 A minimum range in the G band is not required for α2 Canum Venati- online material that includes the time series in the G, BP, and RP corum stars and cataclysmic variables as their variability may be small bands of the selected field-of-view transits for 224 sources that in the ‘white’ G band. are not published in Gaia DR2, but are plotted in Fig. 11.

Article number, page 13 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

Fig. 11. Motions of selected variable stars in the CaMD, highlighted by segments connecting successive absolute G magnitudes and GBP − GRP measurements in time with the same colour for the same source. Preferential directions and amplitudes of magnitude and colour variations can be inferred as a function of variability type (α2 Canum Venaticorum, Be-type and γ Cassiopeiae, cataclysmic, classical and type-II Cepheid, δ Scuti and SX Phoenicis, eclipsing binary, long period, and RR Lyrae), as labelled in the figure. For clarity of visualisation, the selection of eclipsing binaries (and partially other types) was adjusted to minimise the overlap with other types. Selection criteria of all sources represented in colour or grey are the same as in Fig. 3. Additional conditions are described in the text.

Article number, page 14 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

6. Conclusions We gratefully acknowledge Mark Taylor for creating the astronomy-oriented data handling and visualization software TOPCAT (Taylor 2005). The Gaia mission enables a comprehensive description of phe- nomena related to stellar variability. We here focused on stellar variability in the CaMD and showed locations that are occupied by different variability types as well as variable object fractions, variability amplitudes, and variability-induced motions that are described by different variability classes in the CaMD. The wealth of information related to variable stars that is contained in Gaia DR2 is unprecedented for the Milky Way. The CaMD can provide guidance for further detailed studies, which can focus on individual regions or clumps, for instance, to in- vestigate the purity of instability strips and how sharply such regions are truly defined or how they depend on chemical com- position. Of course, additional work is required to this end, and accurately correcting for reddening and extinction will be cru- cial. The (time-dependent) CaMD will play an important role for improving the variable star classification by providing addi- tional attributes, such as the expected direction of variability for specific variable classes, and for illustrating stellar variability to non-expert audiences. The CaMD of variable stars can further point out interre- lations between variability phenomena that are otherwise not easily recognised, and it might be able to identify new types of variability. Detailed follow-up observations from the ground will help correct previous misclassifications and enable in-depth studies of peculiar and particularly interesting objects. Based on the variable stars that reside in the Milky Way, as presented here, it will be possible to obtain data with particularly high S/N, for example through high-resolution spectroscopy. Finally, the ob- served properties of variable stars in the CaMD, such as instabil- ity strip boundaries or period-luminosity relations, provide cru- cial input and constraints for models describing pulsational in- stability, convection, and stellar structure in general. Future Gaia data releases will further surpass the variabil- ity content of this second data release3. By the end of mission, Gaia data are expected to comprise many tens of millions of variable celestial objects, including many additional variability types, as well as time series BP and RP spectra. Eventually, time series of radial velocities and spectra from the radial velocity spectrometer will be published for subsets of variables. Finally, the variability classification of future Gaia data will also make use of unsupervised clustering techniques aimed at discovering entirely new (sub-)clusters and classes of variable phenomena. Acknowledgements. We would like to thank Laurent Rohrbasser for tests done on the representation of time series. This work has made use of data from the ESA space mission Gaia, processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC has been provided by national institutions, some of which participate in the Gaia Multilateral Agreement, which include, for Switzerland, the Swiss State Secretariat for Education, Research and Innovation through the ESA Prodex program, the “Mesures d’accompagnement”, the “Activités Nationales Complémentaires”, the Swiss National Science Foundation, and the Early Postdoc.Mobility fellowship; for Belgium, the BELgian federal Science Policy Office (BELSPO) through PRODEX grants; for Italy, Istituto Nazionale di Astrofisica (INAF) and the Agenzia Spaziale Italiana (ASI) through grants I/037/08/0, I/058/10/0, 2014- 025-R.0, and 2014-025-R.1.2015 to INAF (PI M.G. Lattanzi); for France, the Centre National d’Etudes Spatiales (CNES). Part of this research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Advanced Grant agreements N◦670519: MAMSIE “Mixing and Angular Momentum tranSport in MassIvE stars”). This research has made use of NASA’s Astrophysics Data System, the VizieR catalogue access tool, CDS, Strasbourg, France, and the International Variable Star Index (VSX) database, operated at AAVSO, Cambridge, Massachusetts, USA.

3 cf. https://www.cosmos.esa.int/web/gaia/release

Article number, page 15 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

Appendix A: Literature for each variability type Bellm, E. 2014, in The Third Hot-wiring the Transient Universe Workshop, ed. P. R. Wozniak, M. J. Graham, A. A. Mahabal, & R. Seaman, 27–33 See Table A.1 for details on the references from the literature Bono, G., Castellani, V., & Marconi, M. 2000, ApJ, 529, 293 regarding the objects included in Figs. 3–7 and 11. Bowman, D. M. 2017, Amplitude Modulation of Pulsation Modes in Delta Scuti Stars, Jeremiah Horrocks Institute, University of Central Lancashire, Preston, UK. (PhD Thesis) (Springer Thesis series) Braga, V. F., Stetson, P. B., Bono, G., et al. 2016, AJ, 152, 170 Appendix B: Selection criteria Casertano, S., Riess, A. G., Bucciarelli, B., & Lattanzi, M. G. 2017, A&A, 599, A67 Astrometric and photometric conditions are applied to all Catelan, M. & Smith, H. A. 2015, Pulsating Stars (Wiley-VCH) CaMDs for improved accuracy of the star locations in such di- Chang, S.-W., Byun, Y.-I., & Hartman, J. D. 2015, ApJ, 814, 35 agrams. Astrometric constraints include limits on the number Chaplin, W. J. & Miglio, A. 2013, ARA&A, 51, 353 of visibility periods (observation groups separated from other Conroy, C., Strader, J., van Dokkum, P., et al. 2018, ArXiv e-prints D’Cruz, N. L., Dorman, B., Rood, R. T., & O’Connell, R. W. 1996, ApJ, 466, groups by at least four days) per source used in the secondary 359 astrometric solution (Gaia Collaboration et al. 2018), the excess De Ridder, J., Molenberghs, G., Eyer, L., & Aerts, C. 2016, A&A, 595, L3 astrometric noise of the source postulated to explain the scatter Debosscher, J., Blomme, J., Aerts, C., & De Ridder, J. 2011, A&A, 529, A89 of residuals in the astrometric solution for that source (Gaia Col- Debosscher, J., Sarro, L. M., Aerts, C., et al. 2007, A&A, 475, 1159 Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013a, ApJ, 763, 32 laboration et al. 2018), and the relative parallax precision (herein Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013b, ApJ, 765, 154 set to 20% although in other applications it was set to 5 or 10%): Drake, A. J., Gänsicke, B. T., Djorgovski, S. G., et al. 2014a, MNRAS, 441, 1186 Drake, A. J., Graham, M. J., Djorgovski, S. G., et al. 2014b, ApJS, 213, 9 1. visibility_periods_used > 5; Eggen, O. J. 1951, ApJ, 113, 367 2. astrometric_excess_noise < 0.5 mas; EROS Collaboration, Derue, F., Afonso, C., et al. 1999, A&A, 351, 87 ESA, ed. 1997, ESA Special Publication, Vol. 1200, The Hipparcos and Ty- 3. parallax > 0 mas; cho catalogues. Astrometric and photometric star catalogues derived from the 4. parallax_over_error > 5. ESA Hipparcos Space Astrometry Mission Evans, D., Riello, M., De Angeli, F., & et al. 2018, e-print Photometric conditions set limits for each source on the relative Eyer, L. & Grenon, M. 1997, in ESA Special Publication, Vol. 402, Hipparcos - precisions of the mean fluxes in the GBP, GRP, and G bands, as Venice ’97, ed. R. M. Bonnet, E. Høg, P. L. Bernacca, L. Emiliani, A. Blaauw, well as on the mean flux excess in the G and G bands with C. Turon, J. Kovalevsky, L. Lindegren, H. Hassan, M. Bouffard, B. Strim, BP RP D. Heger, M. A. C. Perryman, & L. Woltjer, 467–472 respect to the G band as a function of colour (Evans et al. 2018): Eyer, L., Grenon, M., Falin, J.-L., Froeschle, M., & Mignard, F. 1994, Sol. Phys., 152, 91 5. phot_bp_mean_flux_error / phot_bp_mean_flux < 0.05; Eyer, L., Guy, L., Distefano, E., et al. 2018, Gaia DR2 documentation Chapter 6. phot_rp_mean_flux_error / phot_rp_mean_flux < 0.05; 7: Variability, Tech. rep. 7. phot_g_mean_flux_error / phot_g_mean_flux < 0.02; Eyer, L. & Mowlavi, N. 2008, Journal of Physics Conference Series, 118, 012010 Eyer, L., Palaversa, L., Mowlavi, N., et al. 2012, Ap&SS, 341, 207 8.( phot_bp_mean_flux + phot_rp_mean_flux) / Fontaine, G. & Brassard, P. 2008, PASP, 120, 1043 {phot_g_mean_flux * [1.2 + 0.03 * (phot_bp_mean_mag Gaia Collaboration, Babusiaux, C., van Leeuwen, F., Barstow, M., & et al. 2018, - phot_rp_mean_mag)2]} < 1.2. e-print Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016a, A&A, 595, A2 The ADQL query to select a sample of sources that satisfy Gaia Collaboration, Clementini, G., Eyer, L., & et al. 2017, A&A, 605, A79 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016b, A&A, 595, A1 all of the above listed criteria follows. Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASP, 122, 131 SELECT TOP 10 source_id Hartman, J. D., Bakos, G. Á., Kovács, G., & Noyes, R. W. 2010, MNRAS, 408, FROM gaiadr2.gaia_source 475 WHERE visibility_periods_used > 5 Hekker, S. & Christensen-Dalsgaard, J. 2017, A&A Rev., 25, 1 Hellier, C. 2001, Cataclysmic Variable Stars (Springer) AND astrometric_excess_noise < 0.5 Hermes, J. J., Kawaler, S. D., Romero, A. D., et al. 2017, ApJ, 841, L2 AND parallax > 0 Holl, B., Audard, M., Nienartowicz, K., et al. 2018, A&A, 618, A30 AND parallax_over_error > 5 Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398 AND phot_bp_mean_flux_over_error > 20 Huber, D., Zinn, J., Bojsen-Hansen, M., et al. 2017, ApJ, 844, 102 Ivezic,´ Ž., Lupton, R. H., Anderson, S., et al. 2003, Mem. Soc. Astron. Italiana, AND phot_rp_mean_flux_over_error > 20 74, 978 AND phot_g_mean_flux_over_error > 50 Jeffery, C. S. 2008, Information Bulletin on Variable Stars, 5817 AND phot_bp_rp_excess_factor < 1.2*(1.2+0.03* Kahraman Aliçavu¸s,F., Niemczura, E., De Cat, P., et al. 2016, MNRAS, 458, power(phot_bp_mean_mag-phot_rp_mean_mag,2)) 2307 Kholopov, P. N. 1956, Peremennye Zvezdy, 11, 325 Kinemuchi, K., Smith, H. A., Wo´zniak,P. R., McKay, T. A., & ROTSE Collab- oration. 2006, AJ, 132, 1202 References Leavitt, H. S. & Pickering, E. C. 1912, Harvard College Observatory Circular, 173, 1 Abbas, M., Grebel, E. K., Martin, N. F., et al. 2014, AJ, 148, 8 Lépine, S. & Gaidos, E. 2011, AJ, 142, 138 Aerts, C. 2015, in IAU Symposium, Vol. 307, New Windows on Massive Stars, LSST Science Collaboration, Abell, P. A., Allison, J., et al. 2009, ArXiv e-prints ed. G. Meynet, C. Georgy, J. Groh, & P. Stee, 154–164 Marconi, M., Musella, I., & Fiorentino, G. 2005, ApJ, 632, 590 Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, Mróz, P., Udalski, A., Poleski, R., et al. 2015, Acta Astron., 65, 313 Astronomy and Astrophysics Library (Springer Science+Business Media Niemczura, E. 2003, A&A, 404, 689 B.V.) Pablo, H., Whittaker, G. N., Popowicz, A., et al. 2016, PASP, 128, 125001 Alcock, C., Allsman, R. A., Axelrod, T. S., et al. 1993, in Astronomical Society Palaversa, L., Ivezic,´ Ž., Eyer, L., et al. 2013, AJ, 146, 101 of the Pacific Conference Series, Vol. 43, Sky Surveys. Protostars to Proto- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323 galaxies, ed. B. T. Soifer, 291 Pigulski, A., Pojmanski,´ G., Pilecki, B., & Szczygieł, D. M. 2009, Acta Astron., Alfonso-Garzón, J., Domingo, A., Mas-Hesse, J. M., & Giménez, A. 2012, 59, 33 A&A, 548, A79 Rauer, H., Catala, C., Aerts, C., et al. 2014, Experimental Astronomy, 38, 249 Anderson, R. I., Saio, H., Ekström, S., Georgy, C., & Meynet, G. 2016, A&A, Reinhold, T. & Gizon, L. 2015, A&A, 583, A65 591 Richards, J. W., Starr, D. L., Miller, A. A., et al. 2012, Astrophys. J. Suppl. Arenou, F., Luri, X., Babusiaux, C., & et al. 2018, e-print Series, 203, 32 Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, Journal of Astronomical Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608 Telescopes, Instruments, and Systems, 1, 014003

Article number, page 16 of 18 Gaia Collaboration et al.: Gaia DR2: CaMD of Variable Stars

Table A.1. Literature references of stars as a function of variability type and the corresponding number of sources depicted in Figs. 3–7, after selections based on reliability, photometric accuracy, and astrometric parameters (Appendix B). Figure 11 includes only subsets of variability types and of sources per type.

Variability Type Reference # Sources Pulsating α Cygni Hip97, VSX16 17 β Cephei PDC05 20 Cepheid ASA09, Hip97, INT12 155 δ Scuti ASA09, Hip97, JD07, Kep11b, Kep11c, SDS12 724 γ Doradus FKA16, Kep11b, Kep11c, VSX16 561 Long Period Variable ASA12, Hip97, INT12, Kep11b, NSV04 5221 PV Telescopii VSX16 3 Rapidly Oscillating Am star VSX16 8 Rapidly Oscillating Ap star VSX16 25 RR Lyrae, fundamental mode (RRab) ASA09, ASA12, Cat13a, Cat13b, Cat14b, Cat15, Hip97, 1676 INT12, LIN13, NSV06, VFB16, VSX16 RR Lyrae, first overtone (RRc) ASA09, ASA12, Cat13b, Cat14b, Hip97, INT12, Kep11b, 611 LIN13, MA14, VFB16, VSX16 RV Tauri ASA12, Hip97, VSX16 48 Slowly Pulsating B star IUE03, Hip97, PDC05 78 SX Phoenicis ASA12, Hip97, VSX16 41 Type-II Cepheid ASA12, Cat14b, Hip97, VSX16 21 V361 Hya (also EC 14026) VSX16 41 V1093 Her (also PG 1716) VSX16 1 ZZ Ceti VSX16 61 Rotational α2 Canum Venaticorum Hip97, VSX16 598 Binary with Reflection VSX16 27 BY Draconis VSX16 713 Ellipsoidal ASA12, Cat14b, Hip97, Kep11b, VSX16 398 FK Comae Berenices Hip97 3 Rotating Spotted Kep15b 16 593 RS Canum Venaticorum ASA12, Cat14b, Hip97, VSX16 1381 Solar-Like Variations HAT10 176 SX Arietis Hip97, VSX16 14 Eclipsing EA, β Persei (Algol) ASA09, Cat14b, Hip97, LIN13, VSX16 8123 EB, β Lyrae ASA09, Cat14b, Hip97, LIN13, VSX16 3096 EW, W Ursae Majoris ASA09, Hip97, VSX16 3248 Exoplanet JS15 278 Eruptive B-type emission-line star ASA12, VSX16 86 Classical T Tauri Star VSX16 75 Flares (UV, BY, TTS) Kep11a, Kep13, Kep15a, MMT15 478 γ Cassiopeiae Hip97, VSX16 84 R Coronae Borealis VSX16 4 S Doradus ASA12, INT12 2 T Tauri Star (TTS) VSX16 173 UV Ceti INT12, VSX16 425 Weak-lined T Tauri Star VSX16 119 Wolf-Rayet INT12, VSX16 15 Cataclysmic Cataclysmic Variable (generic) Cat14a, OGL15, VSX16 132 U Geminorum INT12, VSX16 4 Z Andromedae INT12, VSX16 5

Notes. ASA09: Pigulski et al. (2009); ASA12: Richards et al. (2012); Cat13a: Drake et al. (2013a); Cat13b: Drake et al. (2013b); Cat14a: Drake et al. (2014a); Cat14b: Drake et al. (2014b); Cat15: Torrealba et al. (2015); FKA16: Kahraman Aliçavu¸set al. (2016); HAT10: Hartman et al. (2010); Hip97: ESA (1997); INT12: Alfonso-Garzón et al. (2012); IUE03: Niemczura (2003); JD07: Debosscher et al. (2007); JS15: J. South- worth, http://www.astro.keele.ac.uk/jkt/tepcat/observables.html (as of Aug. 2015); Kep11a: Walkowicz et al. (2011); Kep11b: Debosscher et al. (2011); Kep11c: Uytterhoeven et al. (2011); Kep13: Shibayama et al. (2013); Kep15a: Wu et al. (2015); Kep15b: Reinhold & Gizon (2015); LIN13: Palaversa et al. (2013); MA14: Abbas et al. (2014); MMT15: Chang et al. (2015); NSV04: Wo´zniaket al. (2004); NSV06: Kinemuchi et al. (2006); OGL15: Mróz et al. (2015); PDC05: P. De Cat, http://www.ster.kuleuven.ac.be/~peter/Bstars/ (as of Jan. 2005); SDS12: Süveges et al. (2012); VFB16: Braga et al. (2016); VSX16: Watson et al. (2016, 2015, 2006).

Article number, page 17 of 18 A&A proofs: manuscript no. GaiaCons_EyerEtAl

Rimoldini, L., Holl, B., Audard, M., et al. 2018, arXiv e-prints, (ESA/ESTEC), Keplerlaan 1, 2201 AZ, Noordwijk, The Nether- arXiv:1811.03919 lands Roelens, M., Eyer, L., Mowlavi, N., et al. 2018, arXiv e-prints, 18 Institut de Ciències del Cosmos, Universitat de Barcelona (IEEC- arXiv:1805.00747 UB), Martí i Franquès 1, 08028 Barcelona, Spain Sandage, A. 1958, ApJ, 128, 150 19 Lohrmann Observatory, Technische Universität Dresden, Momm- Shibayama, T., Maehara, H., Notsu, S., et al. 2013, ApJS, 209, 5 senstraße 13, 01062 Dresden, Germany Spano, M., Mowlavi, N., Eyer, L., & Burki, G. 2009, in American Institute of 20 Physics Conference Series, Vol. 1170, American Institute of Physics Confer- European Space Astronomy Centre (ESA/ESAC), Camino bajo del ence Series, ed. J. A. Guzik & P. A. Bradley, 324–326 Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de Stankov, A. & Handler, G. 2005, ApJS, 158, 193 la Cañada, 28692 Madrid, Spain Süveges, M., Sesar, B., Váradi, M., et al. 2012, MNRAS, 424, 2528 21 Lund Observatory, Department of Astronomy and Theoretical Taylor, M. B. 2005, in Astronomical Society of the Pacific Conference Se- Physics, Lund University, Box 43, 22100 Lund, Sweden ries, Vol. 347, Astronomical Data Analysis Software and Systems XIV, ed. 22 Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, P. Shopbell, M. Britton, & R. Ebert, 29 Laboratoire Lagrange, Bd de l’Observatoire, CS 34229, 06304 Torrealba, G., Catelan, M., Drake, A. J., et al. 2015, MNRAS, 446, 2251 Nice Cedex 4, France Udalski, A., Szymanski,´ M. K., & Szymanski,´ G. 2015, Acta Astron., 65, 1 23 Uytterhoeven, K., Moya, A., Grigahcène, A., et al. 2011, A&A, 534, A125 CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 van Leeuwen, F., Evans, D. W., De Angeli, F., et al. 2017, A&A, 599, A32 Toulouse Cedex 9, France 24 Walkowicz, L. M., Basri, G., Batalha, N., et al. 2011, AJ, 141, 50 Institut d’Astronomie et d’Astrophysique, Université Libre de Warner, B. 2003, Cataclysmic Variable Stars (Cambridge University Press) Bruxelles CP 226, Boulevard du Triomphe, 1050 Brussels, Belgium Watson, C. L., Henden, A. A., & Price, A. 2006, Society for Astronomical Sci- 25 F.R.S.-FNRS, Rue d’Egmont 5, 1000 Brussels, Belgium ences Annual Symposium, 25, 47 26 INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, Watson, C. L., Henden, A. A., & Price, A. 2015, CDS/ADC Collection of Elec- 50125 Firenze, Italy tronic Catalogues, 1, 2027 27 Telespazio Vega UK Ltd for ESA/ESAC, Camino bajo del Castillo, Watson, C. L., Wils, P., & Ochsenbein, O. 2016, VizieR Online Data Catalogue, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, VSX Version 2016-07-04 28692 Madrid, Spain Wo´zniak,P. R., Williams, S. J., Vestrand, W. T., & Gupta, V. 2004, AJ, 128, 2965 28 Wu, C.-J., Ip, W.-H., & Huang, L.-C. 2015, ApJ, 798, 92 Laboratoire d’astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France 29 Mullard Space Science Laboratory, University College London, 1 Holmbury St Mary, Dorking, Surrey RH5 6NT, United Kingdom Department of Astronomy, University of Geneva, Chemin des 30 Maillettes 51, 1290 Versoix, Switzerland INAF - Osservatorio Astrofisico di Torino, via Osservatorio 20, 2 10025 Pino Torinese (TO), Italy Department of Astronomy, University of Geneva, Chemin d’Ecogia 31 16, 1290 Versoix, Switzerland Serco Gestión de Negocios for ESA/ESAC, Camino bajo del 3 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain Garching, Germany 32 4 ALTEC S.p.a, Corso Marche, 79,10146 Torino, Italy GEPI, Observatoire de Paris, Université PSL, CNRS, 5 Place Jules 33 Janssen, 92190 Meudon, France Gaia DPAC Project Office, ESAC, Camino bajo del Castillo, s/n, 5 INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain Bologna, via Piero Gobetti 93/3, 40129 Bologna, Italy 34 6 Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, SYRTE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, LNE, 61 avenue de l’Observatoire 75014 Paris, France 3001 Leuven, Belgium 35 7 Department of Astrophysics/IMAPP, Radboud University, P.O.Box National Observatory of Athens, I. Metaxa and Vas. Pavlou, Palaia Penteli, 15236 Athens, Greece 9010, 6500 GL Nijmegen, The Netherlands 36 8 School of Physics and Astronomy, Tel Aviv University, Tel Aviv IMCCE, Observatoire de Paris, Université PSL, CNRS, Sorbonne 6997801, Israel Université, Univ. Lille, 77 av. Denfert-Rochereau, 75014 Paris, 9 France Institute of Astronomy, University of Cambridge, Madingley Road, 37 Cambridge CB3 0HA, United Kingdom Royal Observatory of Belgium, Ringlaan 3, 1180 Brussels, Bel- 10 gium Konkoly Observatory, Research Centre for Astronomy and Earth 38 Sciences, Hungarian Academy of Sciences, Konkoly Thege Miklós Institute for Astronomy, University of Edinburgh, Royal Observa- tory, Blackford Hill, Edinburgh EH9 3HJ, United Kingdom út 15-17, 1121 Budapest, Hungary 39 11 Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Institut d’Astrophysique et de Géophysique, Université de Liège, 19c, Allée du 6 Août, B-4000 Liège, Belgium Leiden, The Netherlands 40 12 INAF - Osservatorio astronomico di Padova, Vicolo Osservatorio ATG Europe for ESA/ESAC, Camino bajo del Castillo, s/n, Urban- 5, 35122 Padova, Italy izacion Villafranca del Castillo, Villanueva de la Cañada, 28692 13 Madrid, Spain Science Support Office, Directorate of Science, European Space 41 Research and Technology Centre (ESA/ESTEC), Keplerlaan 1, Área de Lenguajes y Sistemas Informáticos, Universidad Pablo de Olavide, Ctra. de Utrera, km 1. 41013, Sevilla, Spain 2201AZ, Noordwijk, The Netherlands 42 14 Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France ETSE Telecomunicación, Universidade de Vigo, Campus Lagoas- 15 Marcosende, 36310 Vigo, Galicia, Spain Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidel- 43 berg, Germany Large Synoptic Survey Telescope, 950 N. Cherry Avenue, Tucson, 16 AZ 85719, USA Astronomisches Rechen-Institut, Zentrum für Astronomie der Uni- 44 versität Heidelberg, Mönchhofstr. 12-14, 69120 Heidelberg, Ger- Observatoire Astronomique de Strasbourg, Université de Stras- many bourg, CNRS, UMR 7550, 11 rue de l’Université, 67000 Stras- 17 bourg, France Mission Operations Division, Operations Department, Directorate 45 of Science, European Space Research and Technology Centre Kavli Institute for Cosmology, University of Cambridge, Mading- ley Road, Cambride CB3 0HA, United Kingdom 46 Aurora Technology for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 47 Laboratoire Univers et Particules de Montpellier, Université Mont- pellier, Place Eugène Bataillon, CC72, 34095 Montpellier Cedex 05, France 48 Department of Physics and Astronomy, Division of Astronomy and Article number, page 18 of 18 Space Physics, Uppsala University, Box 516, 75120 Uppsala, Swe- den 49 CENTRA, Universidade de Lisboa, FCUL, Campo Grande, Edif. C8, 1749-016 Lisboa, Portugal 50 Università di Catania, Dipartimento di Fisica e Astronomia, Sezione Astrofisica, Via S. Sofia 78, 95123 Catania, Italy 51 INAF - Osservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy 52 University of Vienna, Department of Astrophysics, Türkenschanzs- traße 17, A1180 Vienna, Austria 53 CITIC – Department of Computer Science, University of A Coruña, Campus de Elviña S/N, 15071, A Coruña, Spain 54 CITIC – Astronomy and Astrophysics, University of A Coruña, Campus de Elviña S/N, 15071, A Coruña, Spain 55 INAF - Osservatorio Astronomico di Roma, Via di Frascati 33, 00078 Monte Porzio Catone (Roma), Italy 56 Space Science Data Center - ASI, Via del Politecnico SNC, 00133 Roma, Italy 57 University of Helsinki, Department of Physics, P.O. Box 64, 00014 Helsinki, Finland 58 Finnish Geospatial Research Institute FGI, Geodeetinrinne 2, 02430 Masala, Finland 59 Isdefe for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 60 Institut UTINAM UMR6213, CNRS, OSU THETA Franche- Comté Bourgogne, Université Bourgogne Franche-Comté, 25000 Besançon, France 61 STFC, Rutherford Appleton Laboratory, Harwell, Didcot, OX11 0QX, United Kingdom 62 Dpto. de Inteligencia Artificial, UNED, c/ Juan del Rosal 16, 28040 Madrid, Spain 63 Elecnor Deimos Space for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 64 Thales Services for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France 65 ON/MCTI-BR, Rua Gal. José Cristino 77, Rio de Janeiro, CEP 20921-400, RJ, Brazil 66 OV/UFRJ-BR, Ladeira Pedro Antônio 43, Rio de Janeiro, CEP 20080-090, RJ, Brazil 67 Department of Terrestrial Magnetism, Carnegie Institution for Sci- ence, 5241 Broad Branch Road, NW, Washington, DC 20015-1305, USA 68 Università di Torino, Dipartimento di Fisica, via Pietro Giuria 1, 10125 Torino, Italy 69 Departamento de Astrofísica, Centro de Astrobiología (CSIC- INTA), ESA-ESAC. Camino Bajo del Castillo s/n. 28692 Vil- lanueva de la Cañada, Madrid, Spain 70 Leicester Institute of Space and Earth Observation and Depart- ment of Physics and Astronomy, University of Leicester, University Road, Leicester LE1 7RH, United Kingdom 71 Departamento de Estadística, Universidad de Cádiz, Calle República Árabe Saharawi s/n. 11510, Puerto Real, Cádiz, Spain 72 Astronomical Institute Bern University, Sidlerstrasse 5, 3012 Bern, Switzerland (present address) 73 EURIX S.r.l., Corso Vittorio Emanuele II 61, 10128, Torino, Italy 74 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA 75 HE Space Operations BV for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 76 Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands 77 SISSA - Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy 78 University of Turin, Department of Computer Sciences, Corso Svizzera 185, 10149 Torino, Italy 79 SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584CA, Utrecht, The Netherlands 80 Dpto. de Matemática Aplicada y Ciencias de la Computación, Univ. de Cantabria, ETS Ingenieros de Caminos, Canales y Puertos, Avda. de los Castros s/n, 39005 Santander, Spain 81 Unidad de Astronomía, Universidad de Antofagasta, Avenida Ang- amos 601, Antofagasta 1270300, Chile 82 CRAAG - Centre de Recherche en Astronomie, Astrophysique et Géophysique, Route de l’Observatoire Bp 63 Bouzareah 16340 Al- giers, Algeria 83 University of Antwerp, Onderzoeksgroep Toegepaste Wiskunde, Middelheimlaan 1, 2020 Antwerp, Belgium 84 INAF - Osservatorio Astronomico d’Abruzzo, Via Mentore Mag- gini, 64100 Teramo, Italy 85 INAF - Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131, Napoli, Italy 86 Instituto de Astronomia, Geofìsica e Ciências Atmosféricas, Uni- versidade de São Paulo, Rua do Matão, 1226, Cidade Universitaria, 05508-900 São Paulo, SP, Brazil 87 Department of Astrophysics, Astronomy and Mechanics, National and Kapodistrian University of Athens, Panepistimiopolis, Zo- grafos, 15783 Athens, Greece 88 Leibniz Institute for Astrophysics Potsdam (AIP), An der Stern- warte 16, 14482 Potsdam, Germany 89 RHEA for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 90 ATOS for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France 91 UNINOVA - CTS, Campus FCT-UNL, Monte da Caparica, 2829- 516 Caparica, Portugal 92 School of Physics, O’Brien Centre for Science North, University College Dublin, Belfield, Dublin 4, Ireland 93 Dipartimento di Fisica e Astronomia, Università di Bologna, Via Piero Gobetti 93/2, 40129 Bologna, Italy 94 Barcelona Supercomputing Center - Centro Nacional de Supercom- putación, c/ Jordi Girona 29, Ed. Nexus II, 08034 Barcelona, Spain 95 Department of Computer Science, Electrical and Space Engineer- ing, Luleå University of Technology, Box 848, S-981 28 Kiruna, Sweden 96 Max Planck Institute for Extraterrestrial Physics, High Energy Group, Gießenbachstraße, 85741 Garching, Germany 97 Astronomical Observatory Institute, Faculty of Physics, Adam Mickiewicz University, Słoneczna 36, 60-286 Poznan,´ Poland 98 Eötvös Loránd University, Egyetem tér 1-3, H-1053 Budapest, Hungary 99 American Community Schools of Athens, 129 Aghias Paraskevis Ave. & Kazantzaki Street, Halandri, 15234 Athens, Greece 100 Faculty of Mathematics and Physics, University of Ljubljana, Jad- ranska ulica 19, 1000 Ljubljana, Slovenia 101 Villanova University, Department of Astrophysics and Planetary Science, 800 E Lancaster Avenue, Villanova PA 19085, USA 102 Physics Department, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerp, Belgium 103 McWilliams Center for Cosmology, Department of Physics, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA 104 Astronomical Institute, Academy of Sciences of the Czech Repub- lic, Fricovaˇ 298, 25165 Ondrejov,ˇ Czech Republic 105 Telespazio for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France 106 Institut de Physique de Rennes, Université de Rennes 1, 35042 Rennes, France 107 Shanghai Astronomical Observatory, Chinese Academy of Sci- ences, 80 Nandan Rd, 200030 Shanghai, China 108 School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China 109 Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen Ø, Denmark 110 DXC Technology, Retortvej 8, 2500 Valby, Denmark 111 Las Cumbres Observatory, 6740 Cortona Drive Suite 102, Goleta, CA 93117, USA 112 Astrophysics Research Institute, Liverpool John Moores Univer- sity, 146 Brownlow Hill, Liverpool L3 5RF, United Kingdom 113 Baja Observatory of University of Szeged, Szegedi út III/70, 6500 Baja, Hungary 114 Laboratoire AIM, IRFU/Service d’Astrophysique - CEA/DSM - CNRS - Université Paris Diderot, Bât 709, CEA-Saclay, 91191 Gif- sur-Yvette Cedex, France 115 École polytechnique fédérale de Lausanne, SB MATHAA STAP, MA B1 473 (Bâtiment MA), Station 8, CH-1015 Lausanne, Switzerland 116 Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warszawa, Poland 117 Institute of Theoretical Physics, Faculty of Mathematics and Physics, Charles University in Prague, Czech Republic 118 AKKA for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France 119 Vitrociset Belgium for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 120 HE Space Operations BV for ESA/ESTEC, Keplerlaan 1, 2201AZ, Noordwijk, The Netherlands 121 Space Telescope Science Institute, 3700 San Martin Drive, Balti- more, MD 21218, USA 122 QUASAR Science Resources for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 123 Fork Research, Rua do Cruzado Osberno, Lt. 1, 9 esq., Lisboa, Por- tugal 124 APAVE SUDEUROPE SAS for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France 125 Nordic Optical Telescope, Rambla José Ana Fernández Pérez 7, 38711 Breña Baja, Spain 126 Spanish Virtual Observatory 127 Fundación Galileo Galilei - INAF, Rambla José Ana Fernández Pérez 7, 38712 Breña Baja, Santa Cruz de Tenerife, Spain 128 INSA for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain 129 Dpto. Arquitectura de Computadores y Automática, Facultad de In- formática, Universidad Complutense de Madrid, C/ Prof. José Gar- cía Santesmases s/n, 28040 Madrid, Spain 130 H H Wills Physics Laboratory, University of Bristol, Tyndall Av- enue, Bristol BS8 1TL, United Kingdom 131 Institut d’Estudis Espacials de Catalunya (IEEC), Gran Capita 2-4, 08034 Barcelona, Spain 132 Applied Physics Department, Universidade de Vigo, 36310 Vigo, Spain 133 Stellar Astrophysics Centre, Aarhus University, Department of Physics and Astronomy, 120 Ny Munkegade, Building 1520, DK- 8000 Aarhus C, Denmark 134 Argelander-Institut fü Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany 135 Research School of Astronomy and Astrophysics, Australian Na- tional University, Canberra, ACT 2611 Australia 136 Sorbonne Universités, UPMC Univ. Paris 6 et CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis bd. Arago, 75014 Paris, France 137 Department of Geosciences, Tel Aviv University, Tel Aviv 6997801, Israel