The ESA Gaia mission provides a unique time-domain survey for more than 1.6 billion sources with G ~ 21 mag. We showcase stellar variability across the Galactic colour-absolute magnitude diagram (CaMD), focusing on pulsating, eruptive, and cataclysmic variables, as well as on stars exhibiting variability due to rotation and eclipses. We illustrate the locations of variable star classes, variable object fractions, and typical variability amplitudes throughout the CaMD and illustrate how variability-related changes in colour and brightness induce `motions' using 22 months worth of calibrated photometric, spectro-photometric, and astrometric Gaia data of stars with significant parallax. To ensure a large variety of variable star classes to populate the CaMD, we crossmatch 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. Gaia enables the first investigation of Galactic variable star populations across the CaMD on a similar, if not larger, scale than previously done in the Magellanic Clouds. Despite observed colours not being reddening corrected, we clearly see distinct regions where variable stars occur and determine variable star fractions to within Gaia's current detection thresholds. Finally, we show the most complete description of variability-induced motion within the CaMD to date. 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 previously accessible only for Galactic star clusters or external galaxies.

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.

The goal of this paper is to demonstrate the outstanding quality of the second data release of the Gaia mission and its power for constraining many different aspects of the dynamics of the satellites of the Milky Way. We focus here on determining the proper motions of 75 Galactic globular clusters, nine dwarf spheroidal galaxies, one ultra-faint system, and the Large and Small Magellanic Clouds. Using data extracted from the Gaia archive, we derived the proper motions and parallaxes for these systems, as well as their uncertainties. We demonstrate that the errors, statistical and systematic, are relatively well understood. We integrated the orbits of these objects in three different Galactic potentials, and characterised their properties. We present the derived proper motions, space velocities, and characteristic orbital parameters in various tables to facilitate their use by the astronomical community. Our limited and straightforward analyses have allowed us for example to (i) determine absolute and very precise proper motions for globular clusters; (ii) detect clear rotation signatures in the proper motions of at least five globular clusters; (iii) show that the satellites of the Milky Way are all on high-inclination orbits, but that they do not share a single plane of motion; (iv) derive a lower limit for the mass of the Milky Way of 9.8^{+6.7}_{-2.7} x 10^{11} Msun based on the assumption that the Leo I dwarf spheroidal is bound; (v) derive a rotation curve for the Large Magellanic Cloud based solely on proper motions that is competitive with line-of-sight velocity curves, now using many orders of magnitude more sources; and (vi) unveil the dynamical effect of the bar on the motions of stars in the Large Magellanic Cloud. All these results highlight the incredible power of the Gaia astrometric mission, and in particular of its second data release.

We present the second Gaia data release, Gaia DR2, consisting of astrometry, photometry, radial velocities, and information on astrophysical parameters and variability, for sources brighter than magnitude 21. In addition epoch astrometry and photometry are provided for a modest sample of minor planets in the solar system. A summary of the contents of Gaia DR2 is presented, accompanied by a discussion on the differences with respect to Gaia DR1 and an overview of the main limitations which are still present in the survey. Recommendations are made on the responsible use of Gaia DR2 results. Gaia DR2 contains celestial positions and the apparent brightness in G for approximately 1.7 billion sources. For 1.3 billion of those sources, parallaxes and proper motions are in addition available. The sample of sources for which variability information is provided is expanded to 0.5 million stars. This data release contains four new elements: broad-band colour information in the form of the apparent brightness in the $G_\mathrm{BP}$ (330--680 nm) and $G_\mathrm{RP}$ (630--1050 nm) bands is available for 1.4 billion sources; median radial velocities for some 7 million sources are presented; for between 77 and 161 million sources estimates are provided of the stellar effective temperature, extinction, reddening, and radius and luminosity; and for a pre-selected list of 14000 minor planets in the solar system epoch astrometry and photometry are presented. Finally, Gaia DR2 also represents a new materialisation of the celestial reference frame in the optical, the Gaia-CRF2, which is the first optical reference frame based solely on extragalactic sources. There are notable changes in the photometric system and the catalogue source list with respect to Gaia DR1, and we stress the need to consider the two data releases as independent.

The Gaia spacecraft of the European Space Agency (ESA) has been securing observations of solar system objects (SSOs) since the beginning of its operations. Gaia Data Release 2 (DR2) contains the observations of a selected sample of 14,099 SSOs. These asteroids have been already identified and have been numbered by the Minor Planet Center. Positions are provided for each Gaia observation at CCD level. As additional information, the apparent brightness of SSOs in the unfiltered G band is also provided for selected observations. We explain the processing of SSO data, and describe the criteria we used to select the sample published in Gaia DR2. We then explore the data set to assess its quality. To exploit the epoch astrometry of asteroids in Gaia DR2 it is necessary to take into account the unusual properties of the uncertainty, as the position information is nearly one-dimensional. When this aspect is handled appropriately, an orbit fit can be obtained with post-fit residuals that are overall consistent with the a-priori error model that was used to define individual values of the astrometric uncertainty. The distribution of residuals allowed us to identify possible contaminants in the data set. Photometry in the G band was compared to computed values from reference asteroid shapes and to the flux registered at the corresponding epochs by the red and blue photometers (RP and BP). The overall astrometric performance is close to the expectations, with an optimal range of brightness G~12-17. In this range, the typical transit-level accuracy is well below 1 mas. For fainter asteroids, the growing photon noise deteriorates the performance. Asteroids brighter than G~12 are affected by a lower performance of the processing of their signals. The dramatic improvement brought by Gaia DR2 astrometry of SSOs is demonstrated by preliminary tests on the detection of subtle non-gravitational effects.

To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $\sim$5~kpc to $\sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)

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.

At about 1000 days after the launch of Gaia we present the first Gaia data release, Gaia DR1, consisting of astrometry and photometry for over 1 billion sources brighter than magnitude 20.7. We summarize Gaia DR1 and provide illustrations of the scientific quality of the data, followed by a discussion of the limitations due to the preliminary nature of this release. Gaia DR1 consists of: a primary astrometric data set which contains the positions, parallaxes, and mean proper motions for about 2 million of the brightest stars in common with the Hipparcos and Tycho-2 catalogues and a secondary astrometric data set containing the positions for an additional 1.1 billion sources. The second component is the photometric data set,consisting of mean G-band magnitudes for all sources. The G-band light curves and the characteristics of ~3000 Cepheid and RR Lyrae stars, observed at high cadence around the south ecliptic pole, form the third component. For the primary astrometric data set the typical uncertainty is about 0.3 mas for the positions and parallaxes, and about 1 mas/yr for the proper motions. A systematic component of ~0.3 mas should be added to the parallax uncertainties. For the subset of ~94000 Hipparcos stars in the primary data set, the proper motions are much more precise at about 0.06 mas/yr. For the secondary astrometric data set, the typical uncertainty of the positions is ~10 mas. The median uncertainties on the mean G-band magnitudes range from the mmag level to ~0.03 mag over the magnitude range 5 to 20.7. Gaia DR1 represents a major advance in the mapping of the heavens and the availability of basic stellar data that underpin observational astrophysics. Nevertheless, the very preliminary nature of this first Gaia data release does lead to a number of important limitations to the data quality which should be carefully considered before drawing conclusions from the data.

Gaia is a cornerstone mission in the science programme of the European Space Agency (ESA). The spacecraft construction was approved in 2006, following a study in which the original interferometric concept was changed to a direct-imaging approach. Both the spacecraft and the payload were built by European industry. The involvement of the scientific community focusses on data processing for which the international Gaia Data Processing and Analysis Consortium (DPAC) was selected in 2007. Gaia was launched on 19 December 2013 and arrived at its operating point, the second Lagrange point of the Sun-Earth-Moon system, a few weeks later. The commissioning of the spacecraft and payload was completed on 19 July 2014. The nominal five-year mission started with four weeks of special, ecliptic-pole scanning and subsequently transferred into full-sky scanning mode. We recall the scientific goals of Gaia and give a description of the as-built spacecraft that is currently (mid-2016) being operated to achieve these goals. We pay special attention to the payload module, the performance of which is closely related to the scientific performance of the mission. We provide a summary of the commissioning activities and findings, followed by a description of the routine operational mode. We summarise scientific performance estimates on the basis of in-orbit operations. Several intermediate Gaia data releases are planned and the data can be retrieved from the Gaia Archive, which is available through the Gaia home page at http://www.cosmos.esa.int/gaia.