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 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.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

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 ESA Gaia satellite scans the whole sky with a temporal sampling ranging from seconds and hours to months. Each time a source passes within the Gaia field of view, it moves over 10 CCDs in 45 s and a lightcurve with 4.5 s sampling (the crossing time per CCD) is registered. Given that the 4.5 s sampling represents a virtually unexplored parameter space in optical time domain astronomy, this data set potentially provides a unique opportunity to open up the fast transient sky. We present a method to start mining the wealth of information in the per CCD Gaia data. We perform extensive data filtering to eliminate known on-board and data processing artefacts, and present a statistical method to identify sources that show transient brightness variations on ~2 hours timescales. We illustrate that by using the Gaia photometric CCD measurements, we can detect transient brightness variations down to an amplitude of 0.3 mag on timescales ranging from 15 seconds to several hours. We search an area of ~23.5 square degrees on the sky, and find four strong candidate fast transients. Two candidates are tentatively classified as flares on M-dwarf stars, while one is probably a flare on a giant star and one potentially a flare on a solar type star. These classifications are based on archival data and the timescales involved. We argue that the method presented here can be added to the existing Gaia Science Alerts infrastructure for the near real-time public dissemination of fast transient events.

We investigate the capabilities of the ESA Gaia mission for detecting and character- izing short timescale variability, from tens of seconds to a dozen hours. We assess the efficiency of the variogram analysis, for both detecting short timescale variability and estimating the underlying characteristic timescales from Gaia photometry, through extensive light-curve simulations for various periodic and transient short timescale variable types. We show that, with this approach, we can detect fast periodic variabil- ity, with amplitudes down to a few millimagnitudes, as well as some M dwarf flares and supernovae explosions, with limited contamination from longer timescale variables or constant sources. Timescale estimates from the variogram give valuable informa- tion on the rapidity of the underlying variation, which could complement timescale estimates from other methods, like Fourier-based periodograms, and be reinvested in preparation of ground-based photometric follow-up of short timescale candidates evi- denced by Gaia. The next step will be to find new short timescale variable candidates from real Gaia data, and to further characterize them using all the Gaia information, including color and spectrum.

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.

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.

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]

Combined studies of variable stars and stellar clusters open great horizons, and they allow us to improve our understanding of stellar cluster formation and stellar evolution. In that prospect, the Gaia mission will provide astrometric, photometric, and spectroscopic data for about one billion stars of the Milky Way. This will represent a major census of stellar clusters, and it will drastically increase the number of known variable stars. In particular, the peculiar Gaia scanning law offers the opportunity to investigate the rather unexplored domain of short timescale variability (from tens of seconds to a dozen of hours), bringing invaluable clues to the fields of stellar physics and stellar aggregates. We assess the Gaia capabilities in terms of short timescale variability detection, using extensive light-curve simulations for various variable object types. We show that Gaia can detect periodic variability phenomena with amplitude variations larger than a few millimagnitudes. Additionally, we plan to perform subsequent follow-up of variables stars detected in clusters by Gaia to better characterize them. Hence, we develop a pipeline for the analysis of high cadence photometry from ground-based telescopes such as the 1.2m Euler telescope (La Silla, Chile) and the 1.2m Mercator telescope (La Palma, Canary Islands).

We investigate the radial velocity (RV) variability and spectroscopic binarity of 19 Galactic long-period ($P_{\rm{puls}} \gtrsim 10$ d) classical Cepheid variable stars whose trigonometric parallaxes are being measured using the Hubble Space Telescope and Gaia. Our primary objective is to constrain possible parallax error due to undetected orbital motion. Using $>1600$ high-precision RVs measured between 2011 and 2016, we find no indication of orbital motion on $\lesssim 5$ yr timescales for 18 Cepheids and determine upper limits on allowed configurations for a range of input orbital periods. The results constrain the unsigned parallax error due to orbital motion to $< 2 \%$ for 16 stars, and $< 4 \%$ for 18. We improve the orbital solution of the known binary YZ Carinae and show that the astrometric model must take into account orbital motion to avoid significant error ($\sim \pm 100 \mu$arcsec). We further investigate long-timescale ($P_{\rm{orb}} > 10$ yr) variations in pulsation-averaged velocity $v_\gamma$ via a template fitting approach using both new and literature RVs. We discover the spectroscopic binarity of XZ Car and CD Cyg, find first tentative evidence for AQ Car, and reveal KN Cen's orbital signature. Further (mostly tentative) evidence of time-variable $v_\gamma$ is found for SS CMa, VY Car, SZ Cyg, and X Pup. We briefly discuss considerations regarding a vetting process of Galactic Leavitt law calibrators and show that light contributions by companions are insignificant for most distance scale applications.

Classical Cepheid variable stars are crucial calibrators of the cosmic distance scale thanks to a relation between their pulsation periods and luminosities. Their archetype, {\delta} Cephei, is an important calibrator for this relation. In this paper, we show that {\delta} Cephei is a spectroscopic binary based on newly-obtained high-precision radial velocities. We combine these new data with literature data to determine the orbit, which has period 2201 days, semi-amplitude 1.5 km/s, and high eccentricity (e = 0.647). We re-analyze Hipparcos intermediate astrometric data to measure {\delta} Cephei's parallax ($\varpi = 4.09 \pm 0.16$ mas) and find tentative evidence for an orbital signature, although we cannot claim detection. We estimate that Gaia will fully determine the astrometric orbit. Using the available information from spectroscopy, velocimetry, astrometry, and Geneva stellar evolution models ($M_{\delta Cep} ~ 5.0 - 5.25 M_\odot$), we constrain the companion mass to within $0.2 < M_2 < 1.2 M_\odot$. We discuss the potential of ongoing and previous interactions between the companion and {\delta} Cephei near pericenter passage, informing reported observations of circumstellar material and bow-shock. The orbit may have undergone significant changes due to a Kozai-Lidov mechanism driven by the outer (visual and astrometric) companion HD 213307. Our discovery of {\delta} Cephei's nature as a spectroscopic binary exposes a hidden companion and reveals a rich and dynamical history of the archetype of classical Cepheid variables.