Stellar variability studies are now reaching a completely new level thanks to ESA's Gaia mission, which enables us to locate many variable stars in the Hertzsprung-Russell diagram and determine the various instability strips/bands. Furthermore, this mission also allows us to detect, characterise and classify many millions of new variable stars thanks to its very unique nearly simultaneous multi-epoch survey with different instruments (photometer, spectro-photometer, radial velocity spectrometer). An overview of what can be found in literature in terms of mostly data products by the Gaia consortium is given. This concerns the various catalogues of variable stars derived from the Gaia time series and also the location and motion of variable stars in the observational Hertzsprung-Russell diagram. In addition, we provide a list of a few thousands of variable white dwarf candidates derived from the DR2 published data, among them probably many hundreds of new pulsating white dwarfs. On a very different topic, we also show how Gaia allows us to reveal the 3D structures of and around the Milky Way thanks to the RR Lyrae stars.

More than half a million of the 1.69 billion sources in Gaia Data Release 2 (DR2) are published with photometric time series that exhibit light variations during the 22 months of observation. An all-sky classification of common high-amplitude pulsators (Cepheids, long-period variables, Delta Scuti / SX Phoenicis, and RR Lyrae stars) is provided for stars with brightness variations greater than 0.1 mag in G band. A semi-supervised classification approach was employed, firstly training multi-stage random forest classifiers with sources of known types in the literature, followed by a preliminary classification of the Gaia data and a second training phase that included a selection of the first classification results to improve the representation of some classes, before the improved classifiers were applied to the Gaia data. Dedicated validation classifiers were used to reduce the level of contamination in the published results. A relevant fraction of objects were not yet sufficiently sampled for reliable Fourier series decomposition, consequently classifiers were based on features derived from statistics of photometric time series in the G, BP, and RP bands, as well as from some astrometric parameters. The published classification results include 195,780 RR Lyrae stars, 150,757 long-period variables, 8550 Cepheids, and 8882 Delta Scuti / SX Phoenicis stars. All of these results represent candidates whose completeness and contamination are described as a function of variability type and classification reliability. Results are expressed in terms of class labels and classification scores, which are available in the vari_classifier_result table of the Gaia archive.

The Gaia second Data Release (DR2) presents a first mapping of full-sky RR Lyrae stars and Cepheids observed by the spacecraft during the initial 22 months of science operations. The Specific Object Study (SOS) pipeline, developed to validate and fully characterise Cepheids and RR Lyrae stars (SOS Cep&RRL) observed by Gaia, has been presented in the documentation and papers accompanying the Gaia first Data Release. Here we describe how the SOS pipeline was modified to allow for processing the Gaia multiband (G, G_BP and G_RP) time series photometry of all-sky candidate variables and produce specific results for confirmed RR Lyrae stars and Cepheids that are published in the DR2 catalogue. The SOS Cep&RRL processing uses tools such as the period-amplitude and the period-luminosity relations in the G band. For the analysis of the Gaia DR2 candidates we also used tools based on the G_BP and G_RP photometry, such as the period-Wesenheit relation in (G,G_RP). Multiband time series photometry and characterisation by the SOS Cep&RRL pipeline are published in Gaia DR2 for 150,359 such variables (9,575 classified as Cepheids and 140,784 as RR Lyrae stars) distributed all over the sky. The sample includes variables in 87 globular clusters and 14 dwarf galaxies. To the best of our knowledge, as of 25 April 2018, variability of 50,570 of these sources (350 Cepheids and 50,220 RR Lyrae stars) is not known in the literature, hence likely they are new discoveries by Gaia. An estimate of the interstellar absorption is published for 54,272 fundamental-mode RR Lyrae stars from a relation based on the G-band amplitude and the pulsation period. Metallicities derived from the Fourier parameters of the light curves are also released for 64,932 RR Lyrae stars and 3,738 fundamental-mode classical Cepheids with period shorter than 6.3 days.

The Gaia DR2 sample of short-timescale variable candidates results from the investigation of the first 22 months of Gaia photometry for a subsample of sources at the Gaia faint end. For this exercise, we limited ourselves to the case of suspected rapid periodic variability. Our study combines fast-variability detection through variogram analysis, high-frequency search by means of least-squares periodograms, and empirical selection based on the investigation of specific sources seen through the Gaia eyes (e.g. known variables or visually identified objects with peculiar features in their light curves). The progressive definition and validation of this selection criterion also benefited from supplementary ground-based photometric monitoring of a few preliminary candidates, performed at the Flemish Mercator telescope (Canary Islands, Spain) between August and November 2017. We publish a list of 3,018 short-timescale variable candidates, spread throughout the sky, with a false-positive rate up to 10-20% in the Magellanic Clouds, and a more significant but justifiable contamination from longer-period variables between 19% and 50%, depending on the area of the sky. Although its completeness is limited to about 0.05%, this first sample of Gaia short-timescale variables recovers some very interesting known short-period variables, such as post-common envelope binaries or cataclysmic variables, and brings to light some fascinating, newly discovered variable sources. In the perspective of future Gaia data releases, several improvements of the short-timescale variability processing are considered, by enhancing the existing variogram and period-search algorithms or by classifying the identified candidates. Nonetheless, the encouraging outcome of our Gaia DR2 analysis demonstrates the power of this mission for such fast-variability studies, and opens great perspectives for this domain of astrophysics.

Gaia DR2 provides a unique all-sky catalogue of 550'737 variable stars, of which 151'761 are long-period variable (LPV) candidates with G variability amplitudes larger than 0.2 mag (5-95% quantile range). About one-fifth of the LPV candidates are Mira candidates, the majority of the rest are semi-regular variable candidates. For each source, G, BP , and RP photometric time-series are published, together with some LPV-specific attributes for the subset of 89'617 candidates with periods in G longer than 60 days. We describe this first Gaia catalogue of LPV candidates, and present various validation checks. Various samples of LPVs were used to validate the catalogue: a sample of well-studied very bright LPVs with light curves from the AAVSO that are partly contemporaneous with Gaia light curves, a sample of Gaia LPV candidates with good parallaxes, the ASAS_SN catalogue of LPVs, and the OGLE catalogues of LPVs towards the Magellanic Clouds and the Galactic bulge. The analyses of these samples show a good agreement between Gaia DR2 and literature periods. The same is globally true for bolometric corrections of M-type stars. The main contaminant of our DR2 catalogue comes from young stellar objects (YSOs) in the solar vicinity (within ~1 kpc), although their number in the whole catalogue is only at the percent level. A cautionary note is provided about parallax-dependent LPV attributes published in the catalogue. This first Gaia catalogue of LPVs approximately doubles the number of known LPVs with amplitudes larger than 0.2 mag, despite the conservative candidate selection criteria that prioritise low contamination over high completeness, and despite the limited DR2 time coverage compared to the long periods characteristic of LPVs. It also contains a small set of YSO candidates, which offers the serendipitous opportunity to study these objects at an early stage of the Gaia data releases.

The Gaia Data Release 2 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability studies.

The second Gaia data release is expected to contain data products from about 22 months of observation. Based on these data, we aim to provide an advance publication of a full-sky Gaia map of RR Lyrae stars. Although comprehensive, these data still contain a significant fraction of sources which are insufficiently sampled for Fourier series decomposition of the periodic light variations. The challenges in the identification of RR Lyrae candidates with (much) fewer than 20 field-of-view transits are described. General considerations of the results, their limitations, and interpretation are presented together with prospects for improvement in subsequent Gaia data releases.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

In this article we present an overview of the ESA Gaia mission and of the unprecedented impact that Gaia will have on the field of variable star research. We summarise the contents and impact of the first Gaia data release on the description of variability phenomena, with particular emphasis on pulsating star research. The Tycho-Gaia astrometric solution, although limited to 2.1 million stars, has been used in many studies related to pulsating stars. Furthermore a set of 3,194 Cepheids and RR Lyrae stars with their times series have been released. Finally we present the plans for the ongoing study of variable phenomena with Gaia and highlight some of the possible impacts of the second data release on variable, and specifically, pulsating stars.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Tens of millions of new variable objects are expected to be identified in over a billion time series from the Gaia mission. Crossmatching known variable sources with those from Gaia is crucial to incorporate current knowledge, understand how these objects appear in the Gaia data, train supervised classifiers to recognise known classes, and validate the results of the Variability Processing and Analysis Coordination Unit (CU7) within the Gaia Data Analysis and Processing Consortium (DPAC). The method employed by CU7 to crossmatch variables for the first Gaia data release includes a binary classifier to take into account positional uncertainties, proper motion, targeted variability signals, and artefacts present in the early calibration of the Gaia data. Crossmatching with a classifier makes it possible to automate all those decisions which are typically made during visual inspection. The classifier can be trained with objects characterized by a variety of attributes to ensure similarity in multiple dimensions (astrometry, photometry, time-series features), with no need for a-priori transformations to compare different photometric bands, or of predictive models of the motion of objects to compare positions. Other advantages as well as some disadvantages of the method are discussed. Implementation steps from the training to the assessment of the crossmatch classifier and selection of results are described.

The ESA Gaia mission provides a unique time-domain survey for more than one billion sources brighter than G=20.7 mag. Gaia offers the unprecedented opportunity to study variability phenomena in the Universe thanks to multi-epoch G-magnitude photometry in addition to astrometry, blue and red spectro-photometry, and spectroscopy. Within the Gaia Consortium, Coordination Unit 7 has the responsibility to detect variable objects, classify them, derive characteristic parameters for specific variability classes, and provide global descriptions of variable phenomena. We describe the variability processing and analysis that we plan to apply to the successive data releases, and we present its application to the G-band photometry results of the first 14 months of Gaia operations that comprises 28 days of Ecliptic Pole Scanning Law and 13 months of Nominal Scanning Law. Out of the 694 million, all-sky, sources that have calibrated G-band photometry in this first stage of the mission, about 2.3 million sources that have at least 20 observations are located within 38 degrees from the South Ecliptic Pole. We detect about 14% of them as variable candidates, among which the automated classification identified 9347 Cepheid and RR Lyrae candidates. Additional visual inspections and selection criteria led to the publication of 3194 Cepheid and RR Lyrae stars, described in Clementini et al. (2016). Under the restrictive conditions for DR1, the completenesses of Cepheids and RR Lyrae stars are estimated at 67% and 58%, respectively, numbers that will significantly increase with subsequent Gaia data releases. Data processing within the Gaia Consortium is iterative, the quality of the data and the results being improved at each iteration. The results presented in this article show a glimpse of the exceptional harvest that is to be expected from the Gaia mission for variability phenomena. [abridged]

We present an overview of the Specific Objects Study (SOS) pipeline developed within the Coordination Unit 7 (CU7) of the Gaia Data Processing and Analysis Consortium (DPAC), the coordination unit charged with the processing and analysis of variable sources observed by Gaia, to validate and fully characterise Cepheids and RR Lyrae stars observed by the spacecraft. We describe how the SOS for Cepheids and RR Lyrae stars (SOS Cep&RRL) was specifically tailored to analyse Gaia's G-band photometric time-series with a South Ecliptic Pole (SEP) footprint, which covers an external region of the Large Magellanic Cloud (LMC). G-band time-series photometry and characterization by the SOS Cep&RRL pipeline (mean magnitude and pulsation characteristics) are published in Gaia Data Release 1 (Gaia DR1) for a total sample of 3,194 variable stars, 599 Cepheids and 2,595 RR Lyrae stars, of which 386 (43 Cepheids and 343 RR Lyrae stars) are new discoveries by Gaia. All 3,194 stars are distributed over an area extending 38 degrees on either side from a point offset from the centre of the LMC by about 3 degrees to the north and 4 degrees to the east. The vast majority, but not all, are located within the LMC. The published sample also includes a few bright RR Lyrae stars that trace the outer halo of the Milky Way in front of the LMC.