Gaia Data Release 2 - Variable stars in the colour-absolute magnitude diagram (original) (raw)
A&A 623, A110 (2019)
Variable stars in the colour-absolute magnitude diagram★,★★
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, R. Haigron4, D. Hatzidimitriou87,35, M. Hauser16,15, M. Haywood4, U. Heiter48, A. Helmi76, J. Heu4, T. Hilger19, 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. Matijevič88, P. J. McMillan21, S. Messina51, D. Michalik21, N. R. Millar9, D. Molina18, R. Molinaro85, L. Molnár10, 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
1Department of Astronomy, University of Geneva, Chemin des Maillettes 51, 1290 Versoix, Switzerland
e-mail: laurent.eyer@unige.ch
2Department of Astronomy, University of Geneva, Chemin d’Ecogia 16, 1290 Versoix, Switzerland
3 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany
4GEPI, Observatoire de Paris, Université PSL, CNRS, 5 Place Jules Janssen, 92190 Meudon, France
5 INAF – Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Piero Gobetti 93/3, 40129 Bologna, Italy
6Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
7Department of Astrophysics/IMAPP, Radboud University, PO Box 9010, 6500 GL Nijmegen, The Netherlands
8School of Physics and Astronomy, Tel Aviv University, Tel Aviv 6997801, Israel
9Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
10Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, Konkoly Thege Miklós út 15-17, 1121 Budapest, Hungary
11Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands
12 INAF – Osservatorio astronomico di Padova, Vicolo Osservatorio 5, 35122 Padova, Italy
13Science Support Office, Directorate of Science, European Space Research and Technology Centre (ESA/ESTEC), Keplerlaan 1, 2201AZ, Noordwijk, The Netherlands
14Université Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
15 Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
16Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstr. 12-14, 69120 Heidelberg, Germany
17Mission Operations Division, Operations Department, Directorate of Science, European Space Research and Technology Centre (ESA/ESTEC), Keplerlaan 1, 2201 AZ, Noordwijk, The Netherlands
18Institut de Ciències del Cosmos, Universitat de Barcelona (IEEC-UB), Martí i Franquès 1, 08028 Barcelona, Spain
19Lohrmann Observatory, Technische Universität Dresden, Mommsenstraße 13, 01062 Dresden, Germany
20European Space Astronomy Centre (ESA/ESAC), Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
21Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 43, 22100 Lund, Sweden
22Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange, Bd de l’Observatoire, CS 34229, 06304 Nice Cedex 4, France
23 CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France
24Institut d’Astronomie et d’Astrophysique, Université Libre de Bruxelles CP 226, Boulevard du Triomphe, 1050 Brussels, Belgium
25F.R.S.-FNRS, Rue d’Egmont 5, 1000 Brussels, Belgium
26 INAF – Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy
27Telespazio Vega UK Ltd for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
28Laboratoire d’astrophysique de Bordeaux, Université Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France
29Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK
30 INAF – Osservatorio Astrofisico di Torino, Via Osservatorio 20, 10025 Pino Torinese (TO), Italy
31Serco Gestión de Negocios for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
32ALTEC S.p.a, Corso Marche, 79, 10146 Torino, Italy
33Gaia DPAC Project Office, ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
34SYRTE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, LNE, 61 avenue de l’Observatoire 75014 Paris, France
35National Observatory of Athens, I. Metaxa and Vas. Pavlou, Palaia Penteli, 15236 Athens, Greece
36IMCCE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Lille, 77 av. Denfert-Rochereau, 75014 Paris, France
37 Royal Observatory of Belgium, Ringlaan 3, 1180 Brussels, Belgium
38Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK
39Institut d’Astrophysique et de Géophysique, Université de Liège, 19c, Allée du 6 Août, 4000 Liège, Belgium
40ATG Europe for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
41Área de Lenguajes y Sistemas Informáticos, Universidad Pablo de Olavide, Ctra. de Utrera, km 1. 41013, Sevilla, Spain
42ETSE Telecomunicación, Universidade de Vigo, Campus Lagoas-Marcosende, 36310 Vigo, Galicia, Spain
43 Large Synoptic Survey Telescope, 950 N. Cherry Avenue, Tucson, AZ 85719, USA
44Observatoire Astronomique de Strasbourg, Université de Strasbourg, CNRS, UMR 7550, 11 rue de l’Université, 67000 Strasbourg, France
45Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambride CB3 0HA, UK
46Aurora Technology for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
47Laboratoire Univers et Particules de Montpellier, Université Montpellier, Place Eugène Bataillon, CC72, 34095 Montpellier Cedex 05, France
48Department of Physics and Astronomy, Division of Astronomy and Space Physics, Uppsala University, Box 516, 75120 Uppsala, Sweden
49CENTRA, Universidade de Lisboa, FCUL, Campo Grande, Edif. C8, 1749-016 Lisboa, Portugal
50Università 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
52University of Vienna, Department of Astrophysics, Türkenschanzstraße 17, A1180 Vienna, Austria
53CITIC – Department of Computer Science, University of A Coruña, Campus de Elviña S/N, 15071, A Coruña, Spain
54CITIC – 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
57Department of Physics, University of Helsinki, PO Box 64, 00014 Helsinki, Finland
58 Finnish Geospatial Research Institute FGI, Geodeetinrinne 2, 02430 Masala, Finland
59Isdefe for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
60Institut UTINAM UMR6213, CNRS, OSU THETA Franche-Comté Bourgogne, Université Bourgogne Franche-Comté, 25000 Besançon, France
61STFC, Rutherford Appleton Laboratory, Harwell, Didcot, OX11 0QX, UK
62Departamento de Inteligencia Artificial, UNED, c/ Juan del Rosal 16, 28040 Madrid, Spain
63Elecnor 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
67Department of Terrestrial Magnetism, Carnegie Institution for Science, 5241 Broad Branch Road, NW, Washington, DC 20015-1305, USA
68Dipartimento di Fisica, Università di Torino, Via Pietro Giuria 1, 10125 Torino, Italy
69Departamento de Astrofísica, Centro de Astrobiología (CSIC-INTA), ESA-ESAC. Camino Bajo del Castillo s/n. 28692 Villanueva de la Cañada, Madrid, Spain
70Leicester Institute of Space and Earth Observation and Department of Physics and Astronomy, University of Leicester, University Road, Leicester LE1 7RH, UK
71Departamento 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
73EURIX S.r.l., Corso Vittorio Emanuele II 61, 10128, Torino, Italy
74 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA
75HE Space Operations BV for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
76Kapteyn 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
78University of Turin, Department of Computer Sciences, Corso Svizzera 185, 10149 Torino, Italy
79SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584CA, Utrecht, The Netherlands
80Departamento de Matemática Aplicada y Ciencias de la Computación, Universidad de Cantabria, ETS Ingenieros de Caminos, Canales y Puertos, Avda. de los Castros s/n, 39005 Santander, Spain
81Unidad de Astronomía, Universidad de Antofagasta, Avenida Angamos 601, Antofagasta 1270300, Chile
82 CRAAG – Centre de Recherche en Astronomie, Astrophysique et Géophysique, Route de l’Observatoire Bp 63 Bouzareah 16340 Algiers, Algeria
83University of Antwerp, Onderzoeksgroep Toegepaste Wiskunde, Middelheimlaan 1, 2020 Antwerp, Belgium
84 INAF – Osservatorio Astronomico d’Abruzzo, Via Mentore Maggini, 64100 Teramo, Italy
85 INAF – Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131, Napoli, Italy
86Instituto de Astronomia, Geofìsica e Ciências Atmosféricas, Universidade de São Paulo, Rua do Matão, 1226, Cidade Universitaria, 05508-900 São Paulo, SP, Brazil
87Department of Astrophysics, Astronomy and Mechanics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos, 15783 Athens, Greece
88 Leibniz Institute for Astrophysics Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
89RHEA 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
91UNINOVA – CTS, Campus FCT-UNL, Monte da Caparica, 2829-516 Caparica, Portugal
92School of Physics, O’Brien Centre for Science North, University College Dublin, Belfield, Dublin 4, Ireland
93Dipartimento di Fisica e Astronomia, Università di Bologna, Via Piero Gobetti 93/2, 40129 Bologna, Italy
94 Barcelona Supercomputing Center – Centro Nacional de Supercomputación, c/ Jordi Girona 29, Ed. Nexus II, 08034 Barcelona, Spain
95Department of Computer Science, Electrical and Space Engineering, Luleå University of Technology, Box 848, 981 28 Kiruna, Sweden
96Max Planck Institute for Extraterrestrial Physics, High Energy Group, Gießenbachstraße, 85741 Garching, Germany
97Astronomical Observatory Institute, Faculty of Physics, Adam Mickiewicz University, Słoneczna 36, 60-286 Poznań, Poland
98 Eötvös Loránd University, Egyetem tér 1-3, 1053 Budapest, Hungary
99 American Community Schools of Athens, 129 Aghias Paraskevis Ave. & Kazantzaki Street, Halandri, 15234 Athens, Greece
100Faculty of Mathematics and Physics, University of Ljubljana, Jadranska ulica 19, 1000 Ljubljana, Slovenia
101Villanova University, Department of Astrophysics and Planetary Science, 800 E Lancaster Avenue, Villanova PA 19085, USA
102Physics Department, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerp, Belgium
103McWilliams Center for Cosmology, Department of Physics, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA
104Astronomical Institute, Academy of Sciences of the Czech Republic, Fričova 298, 25165 Ondřejov, 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
107Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Rd, 200030 Shanghai, PR China
108School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, PR China
109Niels 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
112Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK
113 Baja Observatory of University of Szeged, Szegedi út III/70, 6500 Baja, Hungary
114Laboratoire 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
117Institute 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
119Vitrociset 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, Baltimore, MD 21218, USA
122QUASAR Science Resources for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
123Fork Research, Rua do Cruzado Osberno, Lt. 1, 9 esq., Lisboa, Portugal
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, Spain
127 Fundación Galileo Galilei – INAF, Rambla José Ana Fernández Pérez 7, 38712 Breña Baja, Santa Cruz de Tenerife, Spain
128INSA for ESA/ESAC, Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
129Departamento Arquitectura de Computadores y Automática, Facultad de Informática, Universidad Complutense de Madrid, C/ Prof. José García Santesmases s/n, 28040 Madrid, Spain
130H H Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, UK
131 Institut d’Estudis Espacials de Catalunya (IEEC), Gran Capita 2-4, 08034 Barcelona, Spain
132Applied Physics Department, Universidade de Vigo, 36310 Vigo, Spain
133Stellar Astrophysics Centre, Aarhus University, Department of Physics and Astronomy, 120 Ny Munkegade, Building 1520, 8000 Aarhus C, Denmark
134Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
135Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611 Australia
136Sorbonne Universités, UPMC Université Paris 6 et CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis bd. Arago, 75014 Paris, France
137Department of Geosciences, Tel Aviv University, Tel Aviv 6997801, Israel
Received: 25 April 2018
Accepted: 3 November 2018
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
© ESO 2019
1. Introduction
The ESA space mission Gaia (Gaia Collaboration 2016a) has been conducting a unique survey since the beginning of its operations (end of July 2014). Its uniqueness derives from several aspects that we list in the following paragraphs.
Firstly, Gaia delivers nearly simultaneous measurements in the three observational domains on which most stellar astronomical studies are based: astrometry, photometry, and spectroscopy (Gaia Collaboration 2016b; van Leeuwen et al. 2017). As consequence of the spin of the spacecraft, it takes about 80 s for sources to be measured from the first to the last CCD during a single field-of-view transit.
Secondly, the Gaia data releases provide accurate astrometric measurements for an unprecedented number of objects. In particular, trigonometric parallaxes carry invaluable information, since they are required to infer stellar luminosities, which form the basis of the understanding of much of stellar astrophysics. Proper and orbital motions of stars further enable mass measurements in multiple stellar systems, as well as the investigation of cluster membership.
Thirdly, Gaia data are homogeneous throughout the entire sky, since they are observed with a single set of instruments and are not subject to the Earth’s atmosphere or seasons. All-sky surveys cannot be achieved using a single ground-based telescope; surveys using multiple sites and telescopes and instruments require cross-calibration, which unavoidably introduces systematics and reduces precision because of the increased scatter. Thus, Gaia will play an important role as a standard source in cross-calibrating heterogeneous surveys and instruments, much like the HIPPARCOS mission (Perryman et al. 1997; ESA 1997) did in the past. Of course, Gaia represents a quantum leap from HIPPARCOS in many regards, including an increase of four orders of magnitude in the number of objects observed, additional types of observations (spectrophotometry and spectroscopy), and significantly improved sensitivity and precision for all types of measurements.
Fourthly, there are unprecedented synergies for calibrating distance scales using the dual astrometric and time-domain capabilities of Gaia (e.g. Eyer et al. 2012). Specifically, Gaia will enable the discovery of unrivalled numbers of standard candles residing in the Milky Way, and anchor Leavitt laws (period-luminosity relations) to trigonometric parallaxes (see Gaia Collaboration 2017; Casertano et al. 2017, for two examples based on the first Gaia data release).
Variable stars have for a long time been recognised to offer crucial insights into stellar structure and evolution. Similarly, the Hertzsprung–Russell diagram (HRD) provides an overview of all stages of stellar evolution, and together with its empirical cousin, the colour-magnitude diagram (CMD), it has shaped stellar astrophysics like no other diagram. Henrietta Leavitt (1908) was one of the first to note the immense potential of studying variable stars in populations, where distance uncertainties did not introduce significant scatter. Soon thereafter, Leavitt & Pickering (1912) discovered the period-luminosity relation of Cepheid variables, which has become a cornerstone of stellar physics and cosmology. It appears that Eggen (1951, his Fig. 42) was the first to use (photoelectric) observations of variable stars (in this case, classical Cepheids) to constrain regions where Cepheids occur in the HRD; these regions are today referred to as instability strips. Eggen further illustrated how Cepheids change their locus in the colour-absolute magnitude diagram (CaMD) during the course of their variability, thus developing a time-dependent CMD for variable stars. Kholopov (1956) and Sandage (1958) later illustrated the varying locations of variable stars in the HRD using classical Cepheids located within star clusters. By combining the different types of Gaia time-series data with Gaia parallaxes, we are now in a position to construct time-dependent CaMD towards any direction in the Milky Way, building on previous work based on HIPPARCOS (Eyer et al. 1994; Eyer & Grenon 1997), but on a much larger scale.
Many variability (ground- and space-based) surveys have exploited the power of identifying variable stars in stellar populations at similar distances, for example, in star clusters or nearby galaxies such as the Magellanic Clouds. Ground-based microlensing surveys such as the Optical Gravitational Lensing Experiment (OGLE; e.g. Udalski et al. 2015), the Expérience pour la Recherche d’Objets Sombres (EROS Collaboration 1999), and the Massive Compact Halo Objects project (MACHO; Alcock et al. 1993) deserve a special mention in this regard. The data will continue to grow with the next large multi-epoch surveys such as the Zwicky Transient Facility (Bellm 2014) and the Large Synoptic Survey Telescope (LSST Science Collaboration 2009) from the ground, and the Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015) and PLATO (Rauer et al. 2014) from space.
Another ground-breaking observational trend has been the long-term high-precision high-cadence uninterrupted space photometry with CoRoT/BRITE (Auvergne et al. 2009; Pablo et al. 2016, with time bases of up to five months) and Kepler/K2 (Gilliland et al. 2010; Howell et al. 2014, with time bases of up to four years and three months, respectively) provided entirely new insights into micro-magnitude level variability of stars, with periodicities ranging from minutes to years. These missions opened up stellar interiors from the detection of solar-like oscillations of thousands of sun-like stars and red giants (e.g. Bedding et al. 2011; Chaplin & Miglio 2013; Hekker & Christensen-Dalsgaard 2017, for reviews), as well as hundreds of intermediate-mass stars (e.g. Aerts 2015; Bowman 2017) and compact pulsators (e.g. Hermes et al. 2017). The results we provide in Sects. 3 and 4 on the variability fractions and levels are representative of milli-mag level variability and not of micro-mag levels, as are found in space asteroseismic data.
Any of these asteroseismic surveys can benefit from Gaia astrometry, however, so that distances and luminosities can be derived, as De Ridder et al. (2016) and Huber et al. (2017) reported with Gaia DR1 data. Gaia will also contribute to these surveys with its photometry, and some surveys will also benefit from the Gaia radial velocities (depending on their operating magnitude range).
Stellar variability comprises a great variety of observable features that are due to different physical origins. Figure 1 shows the updated variability tree (Eyer & Mowlavi 2008), which provides a useful overview of the various types of variability and their known causes. The variability tree has four levels: the distinction of intrinsic versus extrinsic variability, the separation into major types of objects (asteroid, stars, and AGN), the physical origin of the variability, and the class name. In this article, we follow the classical distinction of the different causes of the variability phenomena: variability induced by pulsation, rotation, eruption, eclipses, and cataclysmic events. A large number of variability types can be identified in the Gaia data even now, as described in the subsequent sections.
We here provide an overview of stellar variability in the CaMD, building on the astrometric and photometric data of the second Gaia data release (DR2). Future Gaia DRs will enable much more detailed investigations of this kind using longer temporal baselines, greater number of observations, and added classes of variable stars (such as eclipsing binaries, which will be published in DR3).
This paper is structured as follows. Section 2 shows the location of different variability types in the CaMD, making use of known objects from the literature that are published in Gaia DR2, but without any further analysis of the Gaia data. Section 3 presents the fraction of variables as a function of colour and absolute magnitude, obtained by processing the Gaia time series for the detection of variability (Eyer et al. 2018). Section 4 investigates the variability level in the CaMD by employing statistics and classification results (some of which are related to unpublished Gaia time series). Section 5 shows the motion of known variables stars in the CaMD, that is, a time-dependent CaMD, which also includes sources that are not available in the DR2 archive but are online material. Section 6 summarises our results and presents an outlook to future Gaia DRs. Further information on the literature cross-match and on the selection criteria applied to our data samples can be found in Appendices A and B, respectively.
2. Location of variability types in the CaMD
The precision of the location in the CaMD depends on the precision of the colour on the one hand and on the determination of the absolute magnitude on the other. The precision of the absolute magnitude of variable stars depends on the photometric precision, the number of measurements, the amplitude of variability, and the relative parallax precision σ ϖ ∕ϖ. The upper limits of σ ϖ ∕ϖ employed in this article vary between 5 and 20%, which means that the uncertainty of the absolute magnitude that is solely due to the parallax uncertainty can be as large as 5 (ln10)−1σϖ/ϖ≈0.43 mag.
As we determined the colour as a function of integrated BP and integrated RP spectro-photometric measurements with tight constraints on the precision of these quantities (see Appendix B), there are parts of the CaMD that are not explored here. For example, the faint end of the main sequence presented in Fig. 9 of Gaia Collaboration (2018) does not fulfil the condition on the precision in BP, so our diagrams do not include L and T brown dwarfs (which are fainter than M G ~ 14 mag). If we cross-match the Gaia data (conditional on the selection of Appendix B) with the catalogue of M dwarfs of Lépine & Gaidos (2011), only a few M6, M7, and M8 dwarf stars are found.
In Fig. 2 we introduce the Gaia CaMD, which is displayed as a background in subsequent figures. For clarity, we note basic astronomical features such as the main sequence, the red clump (and its long tail due to interstellar extinction), the horizontal branch, the extreme horizontal branch (see D’Cruz et al. 1996, for its physical origin), the red giant branch, the asymptotic giant branch, the white dwarf sequence, the subdwarfs, the supergiants, and the binary sequence. There are additional subtle features above and below the red clump that are described in Fig. 10 of Gaia Collaboration (2018) and are known as the asymptotic giant branch bump and the red giant branch bump, respectively. On the right-hand side of Fig. 2, we also note the typical limiting distance that can be reached because of the selection of σ ϖ ∕ϖ, up to 1 kpc (which was the largest distance we considered for background stars).
Several effects can influence the average location of a star in the CaMD (in both axes), including interstellar extinction, stellar multiplicity, rotation, inclination of the rotation axis, and chemical composition. In this work, we do not correct for such phenomena and instead rely on the apparent magnitudes and colours measured by Gaia, computing “absolute” magnitudes using Gaia parallaxes. We note that interstellar extinction and reddening can be significant at the considered distances (up to 1 kpc), in particular for objects in the Galactic plane. This leads to distortions of certain observed features, such as the long tail in the red clump, which extends to redder and fainter magnitudes.
The stellar variability aspects covered in the second Data Release of Gaia include a limited number of variability classes (Holl et al. 2018): long-period variables, Cepheids, RR Lyrae stars, SX Phoenicis/δ Scuti stars, and rotation-modulated solar-like variability (i.e. all late-type BY Draconis stars). Short-timescale variability (within one day) was explored regardless of the physical origin of the variability (Roelens et al. 2018), although stars that are classified as eclipsing binaries were removed as planned to first appear in the third Data Release of Gaia. The stars presented in this section are solely based on the cross-match with known objects in the literature. The list of variability types presented here is not meant to be comprehensive.
Figures 3–7 illustrate the locations of known variable stars from catalogues in the literature that are cross-matched with the Gaia data. We indicate these targets according to their known variability type published in the literature (the references are listed in Table A.1), and only the stars that satisfy the selection criteria described in Appendix B are kept. Each of these figures includes as reference the location and density (in grey scale) of all stars, regardless of stellar variability, that satisfy the astrometric and photometric criteria of Appendix B with the additional constraint of a minimum parallax of 1 mas (i.e. within 1 kpc to the Sun). This radius seems a good compromise between a large number of stars and a limited effect of interstellar matter. Variable stars whose variability type was previously known are represented by combinations of symbols and colours. Following the structure of the variability tree in Fig. 1, we show in separate figures the CaMDs of stars whose variability is induced by different causes, such as pulsations, rotation, eruptions, eclipses, and cataclysmic events.
Several caveats apply to Figs. 3–7 and should be kept in mind for their interpretation. (a) The quality of catalogues published in the literature can be rather different, in part because variability is often classified without knowledge of a parallax. To reduce the impact of misclassified objects on these figures, we selected subsets of all available catalogues as reference for specific variable star classes, depending on their agreement with the expected locations in the CaMD. In certain cases, we have excluded sources from the literature by choice of specific catalogues (Table A.1) and by using the Gaia astrometry and the multi-band photometric time-series data for occasional cuts in magnitude or colour. Future Gaia data releases will provide a more homogeneous variability classification that will rely primarily on the results of the variability processing (Holl et al. 2018). (b) The CaMDs are not corrected for extinction, which leads to increased scatter, in particular for objects that primarily reside in heavily attenuated areas such as the Galactic disc and the Galactic bulge. (c) The cross-match of sources can be erroneous when stars are located in crowded regions or have high proper motion, especially if the positions of stars in the published catalogues are not sufficiently precise or if proper motion information is not available. (d) Some variability types like magnetically active stars (e.g. RS CVn stars) exhibit different observational phenomena, such as rotational modulation variations as well as flares. To avoid overcrowding the CaMD diagrams, these types are represented in only one of the relevant diagrams (e.g. with rotational or eruptive variables). Furthermore, we note that the time sampling and the waveband coverage of a given survey might favour the detection of only some of these aspects. (e) Gaia represents a milestone for space astrometry and photometry.Nevertheless, some sources can be affected by problems such as corrupt measurements so that their location in the CaMD may be incorrect (Arenou et al. 2018). However, we stress that these problems are limited to a small fraction of sources so that most known variable classes are recovered as expected. The cyclic approach of the Gaia data processing and analysis will allow us to correct for these unexpected features in the future data releases.
Fig. 1Updated 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. |
---|
Fig. 2CaMD 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. |
---|
2.1. Pulsating variable stars
Figure 3 shows the positions of different classes of pulsating variable stars based on the Gaia data and can be compared to its theoretical counterpart in recent textbooks on asteroseismology (Fig. 1.12 in Aerts et al. 2010) and on pulsating stars (Catelan & Smith 2015). We refer to these books for further details of specific variability classes. Here, we only consider the following types of pulsating variable stars:
- 1.
Long-period variables, which are red giant stars that populate the reddest and brightest regions of the CaMD. They include Miras, semi-regular variables, slow irregular variables, and small-amplitude red giants. - 2.
_α_Cygni stars, which are luminous supergiant stars that pulsate in non-radial modes. They are particularly affected by interstellar extinction as they are young massive stars that reside in the Galactic disc, so that their position inFig. 3 must be treated with caution. - 3.
δ Scuti stars, which are Population-I stars of spectral types A and F with short periods (<0.3 d) that dominantly pulsate in pressure modes, but may also reveal low-order gravity modes of low amplitude. - 4.
SX Phoenicis stars, which are Population-II high-amplitude δ Scuti stars. - 5.
γ Doradus stars, which are stars with spectral types A and F with periods from 0.3 to 3 d that dominantly pulsate in high-order gravity modes, but may also reveal low-amplitude pressure modes. - 6.
RR Lyrae stars (Bailey’s type ab and c), which are Population-II horizontal branch stars with periods from 0.2 to 1 d that pulsate in pressure mode. C-type RR Lyrae stars are bluer than ab-type stars. - 7.
Slowly pulsating B (SPB) stars, which are non-radial multi-periodic gravity-mode pulsators of spectral type B with periods typically in the range from 0.4 to 5 d. - 8.
β_Cephei stars, which are late-O to early-B spectral type non-supergiant stars with dominant low-order pressure and gravity modes that feature periods in the range from 0.1 to 0.6 d. Several of them have been found to also exhibit low-amplitude high-order gravity modes as in the SPB stars (e.g. Stankov & Handler 2005). The_β Cephei stars are located in the Galactic disc so that their CaMD position is easily affected by interstellar extinction. - 9.
Classical Cepheids (prototype δ Cephei), which are evolved Population-I (young intermediate-mass) stars featuring radial _p_-mode pulsations with periods of approximately 1_–_100 d. Cepheids can be strongly affected by interstellar extinction as they reside in the Galactic disc and can be observed at great distances. - 10.
Type-II Cepheids, which are Population-II stars pulsating in _p_-mode that were historically thought to be identical to classical Cepheids. Type-II Cepheids consist of three different sub-classes (separated by period) that are commonly referred to as BL Herculis, W Virginis, and RV Tauri stars. Their evolutionary scenarios differ significantly, although the three sub-classes together define a tight period-luminosity relation. - 11.
PV Telescopii stars, which include the sub-classes V652 Her, V2076 Oph, and FQ Aqr (Jeffery 2008). They are rare hydrogen-deficient supergiant stars that cover a wide range of spectral types and exhibit complex light and radial velocity variations. - 12.
Rapidly oscillating Am and Ap stars, which are chemically peculiar A stars that exhibit multiperiodic non-radial pressure modes in the period range of about 5_–_20 min. - 13.
V361 Hydrae (or EC 14026) stars, which are subdwarf B stars on the extreme horizontal branch that pulsate in pressure modes with very short periods of ~1_–_10 min. - 14.
V1093 Her (or PG 1716) stars, which are subdwarf B stars on the extreme horizontal branch that pulsate in gravity modes with periods of about 1_–_3 h. - 15.
ZZ Ceti stars, which are white dwarfs featuring fast non-radial gravity-mode pulsations with periods of 0.5_–_25 min.
The CaMD of pulsating stars carries a great deal of information, much of which has shaped the understanding of stellar structure and evolution and can be found in textbooks. Briefly summarised, we note the following particularly interesting features of Fig. 3.
– Extinction affects variability classes that belong to different populations unequally, as expected. Stars located away from the Galactic disc are much less reddened and thus clump more clearly. This effect is particularly obvious when RR Lyrae stars and classical Cepheids are compared, which both occupy the same instability strip, and it cannot be explained by the known fact that the classical instability strip becomes wider in colour at higher luminosity (e.g. see Anderson et al. 2016; Marconi et al. 2005; Bono et al. 2000, and references therein).
– Interstellar reddening blurs the boundaries between variability classes. Correcting for interstellar extinction will be crucial to delineate the borders of the instability strips in the CaMD, as well as to deduce their purity in terms of the fraction of stars that exhibit pulsations while residing in such regions.
– Practical difficulties involved in separating variable star classes in the way required to construct Fig. 3 include (a) that variable stars are often subject to multiple types of variability (e.g. γ Doradus/δ Scuti, β Cephei/SPB hybrid pulsators, pulsating stars in eclipsing binary systems, or pulsating white dwarfs that exhibit eruptions), and (b) that naming conventions are often historical or purely based on light-curve morphology, so that they do not account for different evolutionary scenarios (e.g. type-II Cepheids). With additional data and a fully homogeneous variable star classification based on Gaia alone, such ambiguities will be resolved in the future unless they are intrinsically connected to the nature of the variability.
– We note multiple groups of ZZ Ceti stars along the white dwarf sequence. The most prominent of these is located at _G_BP −_G_RP ≃ 0 and M G ≃ 12, as reported in Fontaine & Brassard (2008).
Fig. 3Known 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 Table 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 _G_BP measurements introduces artificial cuts in the distributions of low-mass main-sequence stars and red (super)giants. |
---|
Fig. 4Same as Fig. 3, but for rotational-induced variability types. |
---|
Fig. 5Same 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. |
---|
Fig. 7Same as Fig. 3, but for cataclysmic variables and some sub-types. |
---|
2.2. Variability due to rotation and eclipses
Figure 4 shows stars whose variability is induced by rotation. There are three primary categories: spotted stars, stars deformed by tidal interactions, and objects whose variability is due to light reflected by a companion. Following the nomenclature in the literature (Table A.1), we list the following variability classes separately, although we note occasional overlaps among the definitions of these variability classes. The following types are included in Fig. 4:
- 1.
_α_2 Canum Venaticorum stars, which are highly magnetic variable Bp and Ap MS stars. - 2.
Spotted stars, which show rotational modulation variability from spots. - 3.
BY Draconis stars, which are main-sequence stars with late spectral types (K and M) that exhibit quasi-periodic light curves due to spots and chromospheric activity. - 4.
RS Canum Venaticorum stars, which are spotted stars whose rotation-induced variability is frequently accompanied by other phenomena, such as eclipses and flares. - 5.
Ellipsoidal variables, which show variability (without eclipses) due to orbital motion of a star that is distorted by a stellar companion; - 6.
solar-like stars with magnetic activity. Stars of this type in Fig. 4 are limited to a catalogue focused on the Pleiades, which explains the thin distribution of the main sequence. We can see a hint of the binary sequence. - 7.
SX Arietis stars, which are similar to _α_2 Canum Venaticorum stars but have a higher temperature. We note that some overlap of the two distributions occurs for these two variability types. - 8.
Binary systems with a strong reflection component in the light curve with re-radiation of the hotter star’s light from the cooler companion’s surface. - 9.
FK Comae Berenices stars, which are spotted giant stars.
Figure 4 shows the following properties, among other things.
– RS Canum Venaticorum stars are significantly brighter than BY Draconis stars near the bottom of the main sequence (at cool temperatures).
– The reflection binary class is primarily present among very compact (subdwarf) stars; there is a cluster near absolute mag 4, _G_BP –_G_RP approximately − 0.4 mag.
– There seems to be a dearth of rotational spotted variables around _G_BP –_G_RP ~ 0.4, which corresponds with the transition region of stars with a radiative versus convective outer envelope.
– SX Arietis stars form a fairly well-defined hot-temperature envelope of the most luminous _α_2 Canum Venaticorum variables.
Figure 5 shows eclipsing binary systems as well as stars that have been identified to host exoplanets through the transit method. Symbols distinguish the following sub-classes:
- 1.
Eclipsing binaries of type EA; the prototype is Algol. Binaries with spherical or slightly ellipsoidal components with well-separated, nearly constant light curves in between minima. Secondary minima can be absent. - 2.
Eclipsing binaries of type EB; the prototype is β Lyrae. Binaries with continuously changing light curves and not clearly defined onsets or ends of eclipses. Secondary minima are always present, but can be significantly less deep than primary minima. - 3.
Eclipsing binaries of type EW; the prototype is W Ursae Majoris. The components are nearly or actually in contact and minima are virtually equally strong. Onsets and ends of minima are not well defined. - 4.
Stars known to exhibit exoplanetary transits from the literature.
Based on Fig. 5, we observe that EA stars are present almost throughout the CaMD. There are groups of EB stars that are overluminous compared to the white dwarf sequence, which are likely white dwarf stars with main-sequence companions. Moreover, the majority of the stars that host exoplanets are identified by the Kepler spacecraft, and only very few of them have detectable transits in the Gaia data because of the different photometric precision and time sampling.
2.3. Eruptive and cataclysmic variables
Figure 6 focuses on eruptive variable stars. As for the rotationally induced variables, we adopt the nomenclature from the literature (see Table A.1), which includes partially overlapping definitions. The following types are considered.
- 1.
S Doradus stars, also known as luminous blue variables, are massive evolved stars that feature major and irregular photometric variations that are due to heavy mass loss by a radiation-driven wind. - 2.
R Coronae Borealis stars, which are carbon-rich supergiants that emit obscuring material and as a consequence have drastic rapid dimming phases. - 3.
Wolf-Rayet (WR) stars, which are the almost naked helium core that is left over from originally very high mass evolved stars. They feature strong emission lines of hydrogen, nitrogen, carbon, or oxygen. WR stars undergo very fast mass loss and can be significantly dust-attenuated. - 4.
γ Cassiopeiae stars and stars with B spectral types that exhibit hydrogen emission lines, that is, Be stars. These are emitting shell stars. During their “eruptive” phenomena, they become brighter. - 5.
Flare stars, which are magnetically active stars that display flares. This category incudes many subtypes of magnetically active stars, such as UV Ceti-type, RS CVn-type, and T Tauri stars. - 6.
UV Cetistars, which usually are K-M dwarfs that show flares. - 7.
T Tauri stars (classical and weak lined), which are young pre-main sequence stars that either accrete strongly (classical) or show little sign of accretion (weak lined). These stars show variability that is due either to magnetic activity (e.g. rotational modulation and flares) or accretion (quasi-periodic, episodic, or stochastic variations), in addition to pulsations that may also occur in some of them.
In Fig. 6, we notice the absence of eruptive variables among hot main-sequence stars (non-supergiants). This region is populated by pulsating stars, such as γ Doradus and δ Scuti stars, cf. Fig. 3. Moreover, WR stars, R Coronae Borealis stars, and S Doradus stars are among the most luminous stars in this diagram.
Figure 7 illustrates three types of cataclysmic variables.
- 1.
Cataclysmic variables (generic class), typically novae and dwarf novae involving a white dwarf. Many of these stars are situated between the main and white dwarf sequences. - 2.
U Geminorum stars, which are dwarf novae that in principle consist of a white dwarf with a red dwarf companion that experiences mass transfer. - 3.
Z Andromedae stars, which are symbiotic binary stars composed of a giant and a white dwarf.
Further information on cataclysmic variables can be found for example in Warner (2003) and Hellier (2001).
In Fig. 7, we note a clump of cataclysmic variables located in the ZZ Ceti variability strip near G ~ 12 and _G_BP –_G_RP ~ 0.1. The most significant clump of cataclysmic variables is near G ~ 4 and _G_BP − _G_RP ~ 0.1 mag. These are probably binary systems with stars from the extreme horizontal branch and the main sequence.
3. Variable object fractions in the CaMD
The different types of brightness variations as presented in the CaMD may strongly depend on the colour and absolute magnitude, as seen in Sect. 2, because they are driven by different physical mechanisms. Similarly, the variable object fraction, which is defined as the number of variable objects per colour-magnitude bin divided by the total number of objects in the same bin, is expected to depend on the location in the CaMD. The variable object fraction waspreviously determined based on variable objects detected in the HIPPARCOS time series (ESA 1997), for example. Here we significantly expand this investigation using 13.5 million stars with heliocentric distances of up to 1 kpc that satisfy the astrometric and photometric selection criteria listed in Appendix B as well as (a) at least 20 observations in the G, _G_BP, and _G_RP bands, and (b) a relative parallax uncertainty of <5%. In order to reduce the number of objects that are affected by significant extinction, stars at low Galactic latitudes (from −5 to 5°) are excluded. This effectively reduces the number of disc variables such as classical Cepheids and β Cephei stars.
Figure 8 illustrates this _Gaia_-based high-resolution map of the variable object fraction in the CaMD at the precision level of approximately 5–10 mmag. Variability is identified in about 9% of the stars based on a supervised classification of Gaia sources. This method depends heavily on the selection of the training set of constant and variable objects. Minor colour-coded features can be due to training-set related biases. The detection of variability further depends upon the amplitude of the variables, their apparent magnitude distribution, and the instrumental precision. The accuracy of the fraction of variables is also affected by the number of sources per bin of absolute magnitude and colour, which can be as low as one in the tails of the two-dimensional source number density distribution.
Figure 8 contains many informative features, despite possible biases. Future data releases will significantly improve upon Fig. 8 by correcting for reddening and extinction and using a larger number of objects with more accurate source classifications. For the time being, we remark the following.
– The classical instability strip is clearly visible with a variability in about 50–60% of the stars (although extinction limits the precision of this estimate).
– For evolved stars, red giants, and asymptotic giant branch stars, we find that higher luminosity and redder colour implies a higher probability of variability.
– The red clump has a very low fraction of variable stars in the Gaia data. Kepler photometry of red clump stars has revealed complex variability at the micro-mag level that has been used extensively for asteroseismology, cf. Sect. 1 and references therein.
– The classical ZZ Ceti location is extremely concentrated in colour and magnitude, with variability in about half of the stars. The concentration is due to the partial ionisation of hydrogen in the outer envelope of white dwarfs, which is developed only in extremely narrow ranges of effective temperatures (see Fontaine & Brassard 2008).
– Extreme horizontal branch stars show a high probability of variability.
– The hottest and most luminous main-sequence stars are very frequently variable.
– There is a clear gradient towards larger fractions of variables above the low-mass main-sequence stars.
4. Variability amplitudes in the CaMD
Figure 9 shows variability amplitudes as a function of position in the CaMD. Here, we quantify variability amplitudes using the _G_-band inter-quartile range (IQR). Objects are selected according to the general criteria described in Appendix B, with stricter conditions on the parallax (greater than 1 mas) and its relative precision (better than 5%). To prevent the false impression that faint (and very bright) sources have intrinsically higher amplitudes, we corrected for the instrumental spread of the IQR as a function of the median G magnitude. This correction was determined using sources that are classified as constant in the all-sky classification (Rimoldini et al. 2018) and subtracted in quadrature from the measured IQR. Instead of plotting individual data points in Fig. 9, we show the (colour-coded) mean of the corrected _G_-band IQR of sources within each square bin measuring 0.02 mag in both colour and magnitude after trimming the top and bottom 5%. This binning was applied to each variability type individually, and cuts were applied to select minimum classification probabilities per type to minimise incorrect classifications. We emphasise the location of variable object classes that feature large amplitudes by plotting classes with higher IQR on top of variability classes with lower IQR.
Figure 9 contains the following stellar variability types based on the all-sky classification (Rimoldini et al. 2018): _α_2 Canum Venaticorum, α Cygni, β Cephei, cataclysmic, classicalCepheids, δ Scuti, γ Cassiopeiae, γ Doradus, Mira, ellipsoidal, RR Lyrae of Bailey’s type ab and c, semiregular, slowly pulsating B stars, solar-like variability due to magnetic activity (flares, spots, and rotational modulation), SX Arietis, and SX Phoenicis. We did not include other classes (listed in Eyer et al. 2018)for clarity or because there were too few objects. We note that any specific selection criteria applied to the objects shown in Fig. 9 introduce biases that can highlight or diminish the prominence of certain phenomena. Nevertheless, Fig. 9 provides a first detailed illustration of some of the most important amplitude-related variability features in the CaMD. A number of clumps and instability regions are visible in Fig. 9, which are related to the variability classes described in Sects. 2 and 3. We notice the following features.
– The classical instability strip that contains classical Cepheids and RR Lyrae stars is not very prominent, although some clumps (in red or cyan) are apparent.
– The instability regions linked to SPB stars and β Cephei stars are broad and uniform.
– Higher amplitude variations are clearly correlated with redder colours for long-period variables.
– Variables with the highest amplitude (IQR > 0.1 mag) occur in several regions in the CaMD, including the classical instability strip, long-period variables, below the red clump, above the main sequence of low-mass stars (in correspondence of the observed gradient in the fraction of variables), and between the white dwarf sequence and the main sequence.
– Significant amplitudes of >0.04 mag are found very frequently among the coolest white dwarfs.
– The stars between the main sequence and the white dwarfs sequence feature large variability amplitudes and extend into the clump of ZZ Ceti stars in the white dwarf sequence. This intermediate region is populated in particular by the high-amplitude cataclysmic variables, cf. Fig. 7. A closeup view of the white dwarf sequence is shown in Fig. 10, which represents all classified variables within 200 pc. Each object is plotted without binning to emphasise the variability of the ZZ Ceti stars.
Fig. 8Variable 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, _G_BP, and _G_RP 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°) 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 _G_BP measurements introduces artificial cuts in the distributions of low-mass main-sequence stars and red (super)giants. |
---|
Fig. 9Amplitude 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. |
---|
Fig. 10Same as Fig. 9, but focusing on the white dwarf sequence and plotting all classified variables within 200 pc with a relative uncertainty for the parallax better than 5%. A close inspection of this sequence reveals amplitudes at the level of 40 mmag in various regions. |
---|
5. Variability-induced motion in the CaMD
In this section, we visualise the variability-induced motion of stars in the time-dependent CaMD using all-sky measurements made in the G, _G_BP, and _G_RP passbands. Gaia data are uniquely suited to create this time-dependent CaMD, since the different data types (astrometric and photometricin three bands) are acquired in a quasi-simultaneous fashion at many epochs that are distributed over a multi-year time span. The first of such representations, although much less detailed, was presented for individual classical Cepheids in the Milky Way (Eggen 1951) and in Galactic star clusters (Kholopov 1956; Sandage 1958). Similarly minded representations in the literature were based on data from the SDSS (Ivezić et al. 2003, mostly Galactic objects), EROS (Spano et al. 2009, LMC objects), and, very recently, the HST observations of M51 (Conroy et al. 2018).
Figure 11 illustrates the variability-induced motion of stars in the CaMD. As elsewhere in this paper, no correction for interstellar extinction is applied. Individual stars are shown by differently (arbitrarily) coloured lines that connect successive absolute G magnitudes and _G_BP − _G_RP measurements, that is, the observations are ordered in time as opposed to variability phase. This choice was made to avoid uncertainties related to phase-folding the relatively sparsely sampled light curves based on 22 or fewer months of observations and to include both periodic and non-periodic variable objects.
Figure 11 is limited to a subset of all available variable stars in order to avoid overcrowding the diagram. As a preview for future data releases, we include here the variability-induced motions of some stars whose time series and variability types are not published in DR2 (but which are available as online material). Figure 11 includes the following variability types as defined in Sect. 2: _α_2 Canum Venaticorum variables, B-type emission line/γ Cassiopeiae stars, cataclysmic variables, classical and type-II Cepheids, δ Scuti stars, eclipsing binaries, RR Lyrae stars, long-period variables, and SX Phoenicis stars. All sources shown satisfy the general criteria described in Appendix B and typically have at least ten available observations1. We further prioritised the selection of objects featuring wider ranges of variations in the G band (with a minimum of about 0.1 mag)2. The number of sources shown for each variability type ranges from a few to several tens and was selected to ensure clarity in case of high source density or overlapping variability classes in certain regions of the CaMD. In order to limit the effect of outlying values, time-series data are filtered by operators as described in Holl et al. (2018), and the 10 % of the brightest and faintest observations in the _G_RP band are excluded for sources with _G_BP − _G_RP lower than 1.5 mag. Non-variable objects are shown as a grey background to provide a visual reference for the variable object locations in the CaMD. These stars satisfy the criteria described in Appendix B as well as the stricter condition of ϖ > 1 mas. Stars whose variability is caused by different physical effects exhibit different motions within the time-dependent CaMD. We briefly summarise the different motions seen in Fig.11 as follows.
- 1.
Pulsating stars, including long-period variables, Cepheids, RR Lyrae, and δ Scuti/SX Phoenicis stars, exhibit a similar behaviour. These stars are bluer when brighter in G, which illustrates that brightness variations of pulsating stars are 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 generally ordered in phase, so that an overall “envelope” of many cycles is revealed. - 2.
The motions of eclipsing binary systems in the CaMD depend on the colour difference between the two stars. The components of eclipsing systems of the EW type have similar mass (and colour), leading to a rather vertically aligned motion (parallel to the absolute magnitude axis). For eclipsing 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 primary star. The slope of the motion of eclipsing binaries in the CaMD is very different (much steeper) than the one of pulsating stars. - 3.
Rotationally induced variables (shown here: _α_2 Canum Venaticorum stars) feature small amplitudes in absolute _G_and are rather horizontal in the CaMD. This is as expected from starspots, which have a lower temperature than the surroundings and hence absorb the light at bluer frequencies to re-emit it by back-warming effect at redder frequencies. Therefore, the magnitude change in a broad band like_G_is smaller than it would be if measured in narrow bands. - 4.
Eruptive stars (shown here: γ Cassiopeiae and Be-type stars) become redder when brighter because of additional extinction during their eruptive phase. The slopes oftheir motions in the CaMD therefore have the opposite sign with respect to the sign of pulsating stars. - 5.
The variability of cataclysmic variables (shown here: novae) primarily features strong outbursts in the ultraviolet and blue part of the spectrum that are understood to be caused by mass transfer from donor stars in binary systems. These outbursts very significantly change the colour of the system towards bluer values.
The current version of Fig. 11 represents a first step towards a more global description of stellar variability. The motions described by the variable stars in the time-dependent CaMD provide new perspectives on the data that can be exploited as variable star classification attributes to appreciably improve the classification results. Gaia data will definitively help identify misclassifications and problems in published catalogues, thanks to its astrometry and the quasi-simultaneous measurements.
In future Gaia data releases, there will be more data points per source, which will enable us to refine Fig. 11. In particular, periods can be determined with an accuracy inversely proportional to the total Gaia time base for periodic objects. In this way, the motion in the CaMD can be represented more precisely by connecting points that are sorted in phase (rather than in time). This leads to Lissajous-type configurations for pulsators. For sufficiently bright stars,radial velocity time series will add a third and unprecedented dimension to Fig. 11.
An animated version of Fig. 11 is provided online and at https://www.cosmos.esa.int/web/gaia/gaiadr2_cu7. We provide material at the CDS that includes the time series in the G, BP, and RP bands of the selected field-of-view transits for 224 sources that are not published in Gaia DR2, but are plotted in Fig. 11.
Fig. 11Motions of selected variable stars in the CaMD, highlighted by segments connecting successive absolute G magnitudes and _G_BP − _G_RP 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. An animated version of this figure is available online. |
---|
6. Conclusions
The Gaia mission enables a comprehensive description of phenomena 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 investigate the purity of instability strips and how sharply such regions are truly defined or how they depend on chemical composition. Of course, additional work is required to this end, and accurately correcting for reddening and extinction will be crucial. The (time-dependent) CaMD will play an important role for improving the variable star classification by providing additional 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 interrelations 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 observed properties of variable stars in the CaMD, such as instability strip boundaries or period-luminosity relations, provide crucial input and constraints for models describing pulsational instability, convection, and stellar structure in general.
Future Gaia data releases will further surpass the variability 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.
Movie
Movie of Fig. 11.Access here
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. We gratefully acknowledge Mark Taylor for creating the astronomy-oriented data handling and visualization software TOPCAT (Taylor 2005).
Appendix A: Literature for each variability type
See Table A.1 for details on the references from the literature regarding the objects included in Figs. 3–7 and 11.
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).
Appendix B: Selection criteria
Astrometric and photometric conditions are applied to all CaMDs for improved accuracy of the star locations in such diagrams. Astrometric constraints include limits on the number of visibility periods (observation groups separated from other groups by at least four days) per source used in the secondary astrometric solution (Gaia Collaboration 2018), the excess astrometric noise of the source postulated to explain the scatter of residuals in the astrometric solution for that source (Gaia Collaboration 2018), and the relative parallax precision (herein set to 20% although in other applications it was set to 5 or 10%):
- 1.
visibility_periods_used > 5; - 2.
astrometric_excess_noise < 0.5 mas; - 3.
parallax > 0 mas; - 4.
parallax_over_error > 5.
Photometric conditions set limits for each source on the relative precisions of the mean fluxes in the _G_BP, _G_RP, and G bands, as well as on the mean flux excess in the _G_BP and _G_RP bands with respectto the G band as a function of colour (Evans et al. 2018):
- 5.
phot_bp_mean_flux_error/phot_bp_mean_flux < 0.05; - 6.
phot_rp_mean_flux_error/phot_rp_mean_flux < 0.05; - 7.
phot_g_mean_flux_error/phot_g_mean_flux < 0.02; - 8.
(phot_bp_mean_flux + phot_rp_mean_flux) / {phot_g_mean_flux * [1.2 + 0.03 * (phot_bp_mean_mag - phot_rp_mean_mag)2]} < 1.2.
The ADQL query to select a sample of sources that satisfy all of the above listed criteria follows.
SELECT TOP 10 source_id
FROM gaiadr2.gaia_source
WHERE visibility_periods_used > 5
AND astrometric_excess_noise < 0.5
AND parallax > 0
AND parallax_over_error > 5
AND phot_bp_mean_flux_over_error > 20
AND phot_rp_mean_flux_over_error > 20
AND phot_g_mean_flux_over_error > 50
AND phot_bp_rp_excess_factor < 1.2*(1.2+0.03* power(phot_bp_mean_mag-phot_rp_mean_mag,2))
References
- Abbas, M., Grebel, E. K., Martin, N. F., et al. 2014, AJ, 148, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Aerts, C. 2015, in New Windows on Massive Stars, eds. G. Meynet, C. Georgy, J. Groh, & P. Stee, IAU Symp., 307, 154 [NASA ADS] [Google Scholar]
- Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, Astronomy and Astrophysics Library (Dordrecht: Springer Science+Business Media B.V.)[Google Scholar]
- Alcock, C., Allsman, R. A., Axelrod, T. S., et al. 1993, in Sky Surveys. Protostars to Protogalaxies, ed. B. T. Soifer, ASP Conf. Ser., 43, 291 [NASA ADS] [Google Scholar]
- Alfonso-Garzón, J., Domingo, A., Mas-Hesse, J. M., & Giménez, A. 2012, A&A, 548, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Anderson, R. I., Saio, H., Ekström, S., Georgy, C., & Meynet, G. 2016, A&A, 591, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Arenou, F., Luri, X., Babusiaux, C., et al. 2018, A&A, 616, A17 (Gaia 2 SI)[Google Scholar]
- Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Bellm, E. 2014, in The Third Hot-wiring the Transient Universe Workshop, eds. P. R. Wozniak, M. J. Graham, A. A. Mahabal, & R. Seaman, 27[Google Scholar]
- Bono, G., Castellani, V., & Marconi, M. 2000, ApJ, 529, 293 [NASA ADS] [CrossRef] [Google Scholar]
- Bowman, D. M. 2017, Amplitude Modulation of Pulsation Modes in Delta Scuti Stars, Springer Thesis series, Jeremiah Horrocks Institute, University of Central Lancashire, Preston (Ph.D. Thesis) [NASA ADS] [Google Scholar]
- Braga, V. F., Stetson, P. B., Bono, G., et al. 2016, AJ, 152, 170 [NASA ADS] [CrossRef] [Google Scholar]
- Casertano, S., Riess, A. G., Bucciarelli, B., & Lattanzi, M. G. 2017, A&A, 599, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Catelan, M., & Smith, H. A. 2015, Pulsating Stars (Weinheim: Wiley-VCH Verlag GmbH) [CrossRef] [Google Scholar]
- Chang, S.-W., Byun, Y.-I., & Hartman, J. D. 2015, ApJ, 814, 35 [NASA ADS] [CrossRef] [Google Scholar]
- Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353 [NASA ADS] [CrossRef] [Google Scholar]
- Conroy, C., Strader, J., van Dokkum, P., et al. 2018, ApJ, 864, 111 [NASA ADS] [CrossRef] [Google Scholar]
- D’Cruz, N. L., Dorman, B., Rood, R. T., & O’Connell, R. W. 1996, ApJ, 466, 359 [NASA ADS] [CrossRef] [Google Scholar]
- De Ridder, J., Molenberghs, G., Eyer, L., & Aerts, C. 2016, A&A, 595, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Debosscher, J., Sarro, L. M., Aerts, C., et al. 2007, A&A, 475, 1159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Debosscher, J., Blomme, J., Aerts, C., & De Ridder J. 2011, A&A, 529, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013a, ApJ, 763, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Drake, A. J., Catelan, M., Djorgovski, S. G., et al. 2013b, ApJ, 765, 154 [NASA ADS] [CrossRef] [Google Scholar]
- Drake, A. J., Gänsicke, B. T., Djorgovski, S. G., et al. 2014a, MNRAS, 441, 1186 [NASA ADS] [CrossRef] [Google Scholar]
- Drake, A. J., Graham, M. J., Djorgovski, S. G., et al. 2014b, ApJS, 213, 9 [NASA ADS] [CrossRef] [Google Scholar]
- Eggen, O. J. 1951, ApJ, 113, 367 [NASA ADS] [CrossRef] [Google Scholar]
- EROS Collaboration (Derue, F., et al.) 1999, A&A, 351, 87 [NASA ADS] [Google Scholar]
- ESA, ed. 1997, The Hipparcos and Tycho catalogues. Astrometric and Photometric Star Catalogues derived from the ESA Hipparcos Space Astrometry Mission, ESA SP, 1200[Google Scholar]
- Evans, D., Riello, M., De Angeli, F., et al. 2018, A&A, 616, A4 (Gaia 2 SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eyer, L., & Grenon, M. 1997, in Hipparcos – Venice ’97, eds. R. M. Bonnet, E. Høg, P. L. Bernacca, L. Emiliani, A. Blaauw, C. Turon, J. Kovalevsky, L. Lindegren, H. Hassan, M. Bouffard, B. Strim, D. Heger, M. A. C. Perryman, & L. Woltjer, ESA SP, 402, 467–472 [NASA ADS] [Google Scholar]
- Eyer, L., & Mowlavi, N. 2008, J. Phys. Conf. Ser., 118, 012010 [NASA ADS] [CrossRef] [Google Scholar]
- Eyer, L., Grenon, M., Falin, J.-L., Froeschle, M., & Mignard, F. 1994, Sol. Phys., 152, 91 [NASA ADS] [CrossRef] [Google Scholar]
- Eyer, L., Palaversa, L., Mowlavi, N., et al. 2012, Ap&SS, 341, 207 [NASA ADS] [CrossRef] [Google Scholar]
- Eyer, L., Guy,L., Distefano, E., et al. 2018, Gaia DR2 documentation Chapter 7: Variability, Tech. rep[Google Scholar]
- Fontaine, G., & Brassard, P. 2008, PASP, 120, 1043 [NASA ADS] [CrossRef] [Google Scholar]
- Gaia Collaboration (Prusti, T., et al.) 2016a, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2016b, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Clementini, G., et al.) 2017, A&A, 605, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Babusiaux, C., et al) 2018, A&A, 616, A10 (Gaia 2 SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASP, 122, 131 [NASA ADS] [CrossRef] [Google Scholar]
- Hartman, J. D., Bakos, G. Á., Kovács, G., & Noyes, R. W. 2010, MNRAS, 408, 475 [NASA ADS] [CrossRef] [Google Scholar]
- Hekker, S., & Christensen-Dalsgaard, J. 2017, A&ARv, 25, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hellier, C. 2001, Cataclysmic Variable Stars (Dordrecht: Springer)[Google Scholar]
- Hermes, J. J., Kawaler, S. D., Romero, A. D., et al. 2017, ApJ, 841, L2 [NASA ADS] [CrossRef] [Google Scholar]
- Holl, B., Audard, M., Nienartowicz, K., et al. 2018, A&A, 618, A30 (Gaia 2 SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398 [NASA ADS] [CrossRef] [Google Scholar]
- Huber, D., Zinn, J., Bojsen-Hansen, M., et al. 2017, ApJ, 844, 102 [NASA ADS] [CrossRef] [Google Scholar]
- Ivezić, Ž., Lupton, R. H., Anderson, S., et al. 2003, Mem. Soc. Astron. It., 74, 978 [NASA ADS] [Google Scholar]
- Jeffery, C. S. 2008, IBVS Stars, 5817[Google Scholar]
- Kahraman Aliçavuş, F., Niemczura, E., De Cat, P., et al. 2016, MNRAS, 458, 2307 [NASA ADS] [CrossRef] [Google Scholar]
- Kholopov, P. N. 1956, Peremennye Zvezdy, 11, 325 [NASA ADS] [Google Scholar]
- Leavitt, H. S., & Pickering, E. C. 1912, Harvard College Observatory Circ., 173, 1[Google Scholar]
- Lépine, S., & Gaidos, E. 2011, AJ, 142, 138 [NASA ADS] [CrossRef] [Google Scholar]
- LSST Science Collaboration (Abell, P. A., et al.) 2009, ArXiv e-prints [arXiv:0912.0201][Google Scholar]
- Marconi, M., Musella, I., & Fiorentino, G. 2005, ApJ, 632, 590 [NASA ADS] [CrossRef] [Google Scholar]
- Mróz, P., Udalski, A., Poleski, R., et al. 2015, Acta Astron., 65, 313 [NASA ADS] [Google Scholar]
- Niemczura, E. 2003, A&A, 404, 689 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pablo, H., Whittaker, G. N., Popowicz, A., et al. 2016, PASP, 128, 125001 [NASA ADS] [CrossRef] [Google Scholar]
- Palaversa, L., Ivezić, Ž., Eyer, L., et al. 2013, AJ, 146, 101 [CrossRef] [Google Scholar]
- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 [NASA ADS] [Google Scholar]
- Pigulski, A., Pojmański, G., Pilecki, B., & Szczygieł, D. M. 2009, Acta Astron., 59, 33 [NASA ADS] [Google Scholar]
- Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [NASA ADS] [CrossRef] [Google Scholar]
- Reinhold, T., & Gizon, L. 2015, A&A, 583, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Richards, J. W., Starr, D. L., Miller, A. A., et al. 2012, ApJS, 203, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003 [NASA ADS] [CrossRef] [Google Scholar]
- Rimoldini, L., Holl, B., Audard, M., et al. 2018, A&A, submitted (Gaia 2 SI) [arXiv:1811.03919][Google Scholar]
- Roelens, M., Eyer, L., Mowlavi, N., et al. 2018, A&A, 620, A197 (Gaia 2 SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- ROTSE Collaboration (Kinemuchi, K., et al.) 2006, AJ, 132, 1202 [NASA ADS] [CrossRef] [Google Scholar]
- Sandage, A. 1958, ApJ, 128, 150 [NASA ADS] [CrossRef] [Google Scholar]
- Shibayama, T., Maehara, H., Notsu, S., et al. 2013, ApJS, 209, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Spano, M., Mowlavi, N., Eyer, L., & Burki, G. 2009, AIP Conf. Ser., 1170, 324 [NASA ADS] [CrossRef] [Google Scholar]
- Stankov, A., & Handler, G. 2005, ApJS, 158, 193 [NASA ADS] [CrossRef] [Google Scholar]
- Süveges, M., Sesar, B., Váradi, M., et al. 2012, MNRAS, 424, 2528 [NASA ADS] [CrossRef] [Google Scholar]
- Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, ASP Conf. Ser., 347, 29[Google Scholar]
- Torrealba, G., Catelan, M., Drake, A. J., et al. 2015, MNRAS, 446, 2251 [NASA ADS] [CrossRef] [Google Scholar]
- Udalski, A., Szymański, M. K., & Szymański, G. 2015, Acta Astron., 65, 1 [NASA ADS] [Google Scholar]
- Uytterhoeven, K., Moya, A., Grigahcène, A., et al. 2011, A&A, 534, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Leeuwen, F., Evans, D. W., De Angeli, F., et al. 2017, A&A, 599, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Walkowicz, L. M., Basri, G., Batalha, N., et al. 2011, AJ, 141, 50 [NASA ADS] [CrossRef] [Google Scholar]
- Warner, B. 2003, Cataclysmic Variable Stars (Cambridge: Cambridge University Press)[Google Scholar]
- Watson, C. L., Henden, A. A., & Price, A. 2006, Soc. Astron. Sci. Ann. Symp., 25, 47 [NASA ADS] [Google Scholar]
- Watson, C. L., Henden, A. A., & Price, A. 2015, CDS/ADC Collection of Electronic Catalogues, 1, 2027[Google Scholar]
- Watson, C. L., Wils, P., & Ochsenbein, O. 2016, VizieR Online Data Catalogue, VSX Version 2016-07-04[Google Scholar]
- Woźniak, P. R., Williams, S. J., Vestrand, W. T., & Gupta, V. 2004, AJ, 128, 2965 [NASA ADS] [CrossRef] [Google Scholar]
- Wu, C.-J., Ip, W.-H., & Huang, L.-C. 2015, ApJ, 798, 92 [NASA ADS] [CrossRef] [Google Scholar]
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 observations is removed for cataclysmic variables.
2
A minimum range in the _G_band is not required for _α_2 Canum Venaticorum stars and cataclysmic variables as their variability may be small in the “white”_G_band.
All Tables
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).
All Figures
Fig. 1Updated 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. | |
---|---|
In the text |
Fig. 2CaMD 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. | |
---|---|
In the text |
Fig. 3Known 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 Table 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 _G_BP measurements introduces artificial cuts in the distributions of low-mass main-sequence stars and red (super)giants. | |
---|---|
In the text |
Fig. 5Same 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. | |
---|---|
In the text |
Fig. 8Variable 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, _G_BP, and _G_RP 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°) 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 _G_BP measurements introduces artificial cuts in the distributions of low-mass main-sequence stars and red (super)giants. | |
---|---|
In the text |
Fig. 9Amplitude 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. | |
---|---|
In the text |
Fig. 10Same as Fig. 9, but focusing on the white dwarf sequence and plotting all classified variables within 200 pc with a relative uncertainty for the parallax better than 5%. A close inspection of this sequence reveals amplitudes at the level of 40 mmag in various regions. | |
---|---|
In the text |
Fig. 11Motions of selected variable stars in the CaMD, highlighted by segments connecting successive absolute G magnitudes and _G_BP − _G_RP 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. An animated version of this figure is available online. | |
---|---|
In the text |