Combining the precise parallaxes and optical photometry delivered by Gaia's second data release (Gaia DR2) with the photometric catalogues of PanSTARRS-1, 2MASS, and AllWISE, we derive Bayesian stellar parameters, distances, and extinctions for 265 million stars brighter than G=18. Because of the wide wavelength range used, our results substantially improve the accuracy and precision of previous extinction and effective temperature estimates. After cleaning our results for both unreliable input and output data, we retain 137 million stars, for which we achieve a median precision of 5% in distance, 0.20 mag in V-band extinction, and 245 K in effective temperature for G<14, degrading towards fainter magnitudes (12%, 0.20 mag, and 245 K at G=16; 16%, 0.23 mag, and 260 K at G=17, respectively). We find a very good agreement with the asteroseismic surface gravities and distances of 7000 stars in the Kepler, the K2-C3, and the K2-C6 fields, with stellar parameters from the APOGEE survey, as well as with distances to star clusters. Our results are available through the ADQL query interface of the Gaia mirror at the Leibniz-Institut f\"{u}r Astrophysik Potsdam (gaia.aip.de) and as binary tables at data.aip.de. As a first application, in this paper we provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps as a function of distance, and extensive density maps, demonstrating the potential of our value-added dataset for mapping the three-dimensional structure of our Galaxy. In particular, we see a clear manifestation of the Galactic bar in the stellar density distributions, an observation that can almost be considered a direct imaging of the Galactic bar.

There are few warp kinematic models of the Galaxy able to characterise structure and kinematics. These models are necessary to study the lopsidedness of the warp and the twisting of the line-of-nodes of the stellar warp, already seen in gas and dust. We use the \Gaia~Data Release 2 astrometric data up to $G=20$mag to characterise the structure of the Galactic warp, the vertical motions and the dependency on the age. We use two populations up to galactocentric distances of $16$kpc, a young (OB-type) and an old (Red Giant Branch, RGB). We use the nGC3 PCM and LonKin methods based on the Gaia observables, together with 2D projections of the positions and proper motions in the Galactic plane. We confirm the age dependency of the Galactic warp, both in positions and kinematics, being the height of the Galactic warp of about $0.2$kpc for the OB sample and of $1.$kpc for the RGB at a galactocentric distance of $14$kpc. Both methods find that the onset radius is $12\sim 13$kpc for the OB sample and $10\sim 11$kpc for the RGB. From the RGB sample, we find from galactocentric distances larger than $10$kpc the line-of-nodes twists away from the Sun-anticentre line towards galactic azimuths $\sim 180-200^{\circ}$ increasing with radius, though possibly influenced by extinction. The RGB sample reveals a slightly lopsided stellar warp with $\sim 250$pc between the up and down sides. The line of maximum of proper motions in latitude is systematically offset from the line-of-nodes estimated from the spatial data, which our models predict as a kinematic signature of lopsidedness. We also show a prominent wave-like pattern of a bending mode different in the OB and RGB, and substructures that might not be related to the Galactic warp nor to a bending mode. GDR2 triggers the need for complex kinematic models, flexible enough to combine both wave-like patterns and an S-shaped lopsided warp.[abridged]

Context. Open clusters are very good tracers of the evolution of the Galactic disc. Thanks to Gaia, their kinematics can be investigated with an unprecedented precision and accuracy. Aims. The distribution of open clusters in the 6D phase space is revisited with Gaia DR2. Methods. The weighted mean radial velocity of open clusters was determined, using the most probable members available from a previous astrometric investigation that also provided mean parallaxes and proper motions. Those parameters, all derived from Gaia DR2 only, were combined to provide the 6D phase space information of 861 clusters. The velocity distribution of nearby clusters was investigated, as well as the spatial and velocity distributions of the whole sample as a function of age. A high quality subsample was used to investigate some possible pairs and groups of clusters sharing the same Galactic position and velocity. Results. For the high quality sample that has 406 clusters, the median uncertainty of the weighted mean radial velocity is 0.5 km/s. The accuracy, assessed by comparison to ground-based high resolution spectroscopy, is better than 1 km/s. Open clusters nicely follow the velocity distribution of field stars in the close Solar neighbourhood previously revealed by Gaia DR2. As expected, the vertical distribution of young clusters is very flat but the novelty is the high precision to which this can be seen. The dispersion of vertical velocities of young clusters is at the level of 5 km/s. Clusters older than 1 Gyr span distances to the Galactic plane up to 1 kpc with a vertical velocity dispersion of 14 km/s, typical of the thin disc. Five pairs of clusters and one group with five members are possibly physically related. Other binary candidates previously identified turn out to be chance alignment.

SN 1572 (Tycho Brahe's supernova) clearly belongs to the Ia (thermonuclear) type. It was produced by the explosion of a white dwarf in a binary system. Its remnant has been the first of this type to be explored in search of a possible surviving companion, the mass donor that brought the white dwarf to the point of explosion. A high peculiar motion with respect to the stars at the same location in the Galaxy, mainly due to the orbital velocity at the time of the explosion, is a basic criterion for the detection of such companions. Radial velocities from the spectra of the stars close to the geometrical center of Tycho's supernova remnant, plus proper motions of the same stars, obtained by astrometry with the {\it Hubble Space Telescope}, have been used so far. In addition, a detailed chemical analysis of the atmospheres of a sample of candidate stars had been made. However, the distances to the stars, remained uncertain. Now, the Second {\it Gaia} Data Release (DR2) provides unprecedent accurate distances and new proper motions for the stars can be compared with those made from the {\it HST}. We consider the Galactic orbits that the candidate stars to SN companion would have in the future. We do this to explore kinematic peculiarity. We also locate a representative sample of candidate stars in the Toomre diagram. Using the new data, we reevaluate here the status of the candidates suggested thus far, as well as the larger sample of the stars seen in the central region of the remnant.

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.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

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)

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

The Young Local Associations constitute an excellent sample for the study of a variety of astrophysical topics, especially the star formation process in low-density environments. Data from the Gaia mission allows us to undertake studies of the YLAs with unprecedented accuracy. We determine the dynamical age and place of birth of a set of associations in a uniform and dynamically consistent manner. There are nine YLAs in our sample $\epsilon$ Chamaeleontis, TW Hydrae, $\beta$ Pictoris, Octans, Tucana-Horologium, Columba, Carina, Argus and AB Doradus. We designed a method for deriving the dynamical age of the YLAs based on the orbital integration. It involves a strategy to account for the effect of observational errors and we tested it using mock YLAs. Finally, we applied it to our set of nine YLAs with astrometry from the first Gaia data release and complementary on-ground radial velocities from the literature. Our orbital analysis yields a first estimate of the dynamical age of 3$^{+9}_{-0}$ Myr, 13$^{+7}_{-0}$ Myr and 5$^{+23}_{-0}$ Myr for $\epsilon$ Cha, $\beta$ Pict and Tuc-Hor, respectively. For four other associations (Oct, Col, Car and Arg), we provide a lower limit for the dynamical age. Our rigorous error treatment indicates that TW Hya and AB Dor deserve further study. Our dynamical ages are compatible spectroscopic and isochrone fitting ages obtained elsewhere. From the orbital analysis, we suggest a scenario with two episodes of star formation: one ~40 Myr ago in the first quadrant that gave birth to $\epsilon$ Cha, TW Hya and $\beta$ Pic, and another 5-15 Myr ago close to the Sun that formed Tuc-Hor, Col, and Car. Future Gaia data will provide the necessary accuracy to improve the present results, especially for the controversial age determinations, and additional evidence for the proposed scenario, once a complete census of YLAs and better membership can be obtained.

The stellar [$\alpha$/Fe] abundance is sometimes used as a proxy for stellar age, following standard chemical evolution models for the Galaxy, as seen by different observational results. In this work we show that the Open Cluster NGC6705/M11 has a significant $\alpha$-enhancement [$\alpha$/Fe]$>0.1$ dex, despite its young age ($\sim$300 Myr), challenging the current paradigm. We use high resolution (R$>65,000$) high signal-to-noise ($\sim$70) spectra of 8 Red Clump stars, acquired within the OCCASO survey. We determine very accurate chemical abundances of several $\alpha$ elements, using an equivalent width methodology (Si, Ca and Ti), and spectral synthesis fits (Mg and O). We obtain [Si/Fe]=$0.13\pm0.05$, [Mg/Fe]=$0.14\pm0.07$, [O/Fe]=$0.17\pm0.07$, [Ca/Fe]=$0.06\pm0.05$ and [Ti/Fe]=$0.03\pm0.03$. Our results place these cluster within the group of young [$\alpha$/Fe]-enhanced field stars recently found by several authors in the literature. The ages of our stars have an uncertainty of around 50 Myr, much more precise than for field stars. By integrating the cluster's orbit in several non-axisymmetric Galactic potentials, we establish the M11's most likely birth radius to lie between 6.8-7.5 kpc from the Galactic center, not far from its current position. With the robust Open Cluster age scale, our results prove that a moderate [$\alpha$/Fe]-enhancement is no guarantee for a star to be old, and that not all $\alpha$-enhanced stars can be explained with an evolved blue straggler scenario. Based on our orbit calculations, we further argue against a Galactic bar origin of M11.

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.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

NGC 7067 is a young open cluster located in the direction between the first and the second Galactic quadrants and close to the Perseus spiral arm. This makes it useful for studies of the nature of the Milky Way spiral arms. Stromgren photometry taken with the Wide Field Camera at the Isaac Newton Telescope allowed us to compute individual physical parameters for the observed stars and hence to derive cluster's physical parameters. Spectra from the 1.93-m telescope at the Observatoire de Haute-Provence helped to check and improve the results. We obtained photometry for 1233 stars, individual physical parameters for 515 and spectra for 9 of them. The 139 selected cluster members lead to a cluster distance of 4.4+/-0.4 kpc, with an age below log10(t(yr))=7.3 and a present Mass of 1260+/-160Msun. The morphology of the data reveals that the centre of the cluster is at (ra,dec)=(21:24:13.69,+48:00:39.2) J2000, with a radius of 6.1arcsec. Stromgren and spectroscopic data allowed us to improve the previous parameters available for the cluster in the literature.

We map the kinematics of stars in simulated galaxy disks with spiral arms using the velocity ellipsoid vertex deviation (l$_v$). We use test particle simulations, and for the first time, fully self-consistent high resolution N-body models. We compare our maps with the Tight Winding Approximation model analytical predictions. We see that for all barred models spiral arms rotate closely to a rigid body manner and the vertex deviation values correlate with the density peaks position bounded by overdense and underdense regions. In such cases, vertex deviation sign changes from negative to positive when crossing the spiral arms in the direction of disk rotation, in regions where the spiral arms are in between corotation (CR) and the Outer Lindblad Resonance (OLR). By contrast, when the arm sections are inside the CR and outside the OLR, l$_v$ changes from negative to positive.We propose that measurements of the vertex deviations pattern can be used to trace the position of the main resonances of the spiral arms. We propose that this technique might exploit future data from Gaia and APOGEE surveys. For unbarred N-body simulations with spiral arms corotating with disk material at all radii, our analysis suggests that no clear correlation exists between l$_v$ and density structures.

High resolution N-body simulations using different codes and initial condition techniques reveal two different behaviours for the rotation frequency of transient spiral arms like structures. Whereas unbarred disks present spiral arms nearly corotatingwith disk particles, strong barred models (bulged or bulge-less) quickly develop a bar-spiral structure dominant in density, with a pattern speed almost constant in radius. As the bar strength decreases the arm departs from bar rigid rotation and behaves similar to the unbarred case. In strong barred models we detect in the frequency space other subdominant and slower modes at large radii, in agreement with previous studies, however we also detect them in the configuration space. We propose that the distinctive behaviour of the dominant spiral modes can be exploited in order to constraint the nature of Galactic spiral arms by the astrometric survey GAIA and by 2-D spectroscopic surveys like CALIFA and MANGA in external galaxies.

With the advent of the Gaia data, the unprecedented kinematic census of great part of the Milky Way disc will allow us to characterise the local kinematic groups and new groups in different disc neighbourhoods. First, we show here that the models predict a stellar kinematic response to the spiral arms and bar strongly dependent on disc position. For example, we find that the kinematic groups induced by the spiral arm models change significantly if one moves only ~ 0.6 kpc in galactocentric radius, but ~ 2 kpc in azimuth. There are more and stronger groups as one approaches the spiral arms. Depending on the spiral pattern speed, the kinematic imprints are more intense in nearby vicinities or far from the Sun. Secondly, we present a preliminary study of the kinematic groups observed by RAVE. This sample will allow us, for the first time, to study the dependence on Galactic position of the (thin and thick) disc moving groups. In the solar neighbourhood, we find the same kinematics groups as detected in previous surveys, but now with better statistics and over a larger spatial volume around the Sun. This indicates that these structures are indeed large scale kinematic features.

We study the spiral arm influence on the solar neighbourhood stellar kinematics. As the nature of the Milky Way (MW) spiral arms is not completely determined, we study two models: the Tight-Winding Approximation (TWA) model, which represents a local approximation, and a model with self-consistent material arms named PERLAS. This is a mass distribution with more abrupt gravitational forces. We perform test particle simulations after tuning the two models to the observational range for the MW spiral arm properties. We explore the effects of the arm properties and find that a significant region of the allowed parameter space favours the appearance of kinematic groups. The velocity distribution is mostly sensitive to the relative spiral arm phase and pattern speed. In all cases the arms induce strong kinematic imprints for pattern speeds around 17 km/s/kpc (close to the 4:1 inner resonance) but no substructure is induced close to corotation. The groups change significantly if one moves only ~0.6 kpc in galactocentric radius, but ~2 kpc in azimuth. The appearance time of each group is different, ranging from 0 to more than 1 Gyr. Recent spiral arms can produce strong kinematic structures. The stellar response to the two potential models is significantly different near the Sun, both in density and kinematics. The PERLAS model triggers more substructure for a larger range of pattern speed values. The kinematic groups can be used to reduce the current uncertainty about the MW spiral structure and to test whether this follows the TWA. However, groups such as the observed ones in the solar vicinity can be reproduced by different parameter combinations. Data from velocity distributions at larger distances are needed for a definitive constraint.