Using 2MASS photometry of Gaia DR2 sources, we present a technique for selecting upper main sequence stars and giants without the need for individual extinction estimates, to a distance of 7 kpc from the Sun. The spatial distribution of the upper main sequence stars clearly shows the nearest spiral arms, while the large-scale kinematics of the two populations perpendicular to the Galactic plane both show for the first time a clear signature of the warp of the Milky Way.

Most of the stars in our Galaxy including our Sun, move in a disk-like component and give the Milky Way its characteristic appearance on the night sky. As in all fields in science, motions can be used to reveal the underlying forces, and in the case of disk stars they provide important diagnostics on the structure and history of the Galaxy. But because of the challenges involved in measuring stellar motions, samples have so far remained limited in their number of stars, precision and spatial extent. This has changed dramatically with the second Data Release of the Gaia mission which has just become available. Here we report that the phase space distribution of stars in the disk of the Milky Way is full of substructure with a variety of morphologies, most of which have never been observed before. This includes shapes such as arches and shells in velocity space, and snail shells and ridges when spatial and velocity coordinates are combined. The nature of these substructures implies that the disk is phase mixing from an out of equilibrium state, and that the Galactic bar and/or spiral structure are strongly affecting the orbits of disk stars. Our analysis of the features leads us to infer that the disk was perturbed between 300 and 900 Myr ago, which matches current estimations of the previous pericentric passage of the Sagittarius dwarf galaxy. The Gaia data challenge the most basic premise of stellar dynamics of dynamical equilibrium, and show that modelling the Galactic disk as a time-independent axisymmetric component is definitively incorrect. These findings mark the start of a new era when, by modelling the richness of phase space substructures, we can determine the gravitational potential of the Galaxy, its time evolution and the characteristics of the perturbers that have most influenced our home in the Universe.

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.

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.

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 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 explore the possibility of detecting and characterizing the warp of the stellar disc of our Galaxy using synthetic Gaia data. The availability of proper motions and, for the brightest stars radial velocities, adds a new dimension to this study. A family of Great Circle Cell Counts (GC3) methods is used. They are ideally suited to find the tilt and twist of a collection of rings, which allow us to detect and measure the warp parameters. To test them, we use random realizations of test particles which evolve in a realistic Galactic potential warped adiabatically to various final configurations. In some cases a twist is introduced additionally. The Gaia selection function, its errors model and a realistic 3D extinction map are applied to mimic three tracer populations: OB, A and Red Clump stars. We show how the use of kinematics improves the accuracy in the recovery of the warp parameters. The OB stars are demonstrated to be the best tracers determining the tilt angle with accuracy better than $\sim0.5$ up to Galactocentric distance of $\sim16$ kpc. Using data with good astrometric quality, the same accuracy is obtained for A type stars up to $\sim13$ kpc and for Red Clump up to the expected stellar cut-off. Using OB stars the twist angle is recovered to within $< 3^\circ$ for all distances.

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 bar-driven dynamics in the inner part of the Milky Way by using invariant manifolds. This theory has been successfully applied to describe the morphology and kinematics of rings and spirals in external galaxies, and now, for the first time, we apply it to the Milky Way. We compute the orbits confined by the invariant manifolds of the unstable periodic orbits located at the ends of the bar. We discuss whether the COBE/DIRBE bar and the Long bar compose a single bar or two independent bars and perform a number of comparisons which, taken together, argue strongly in favour of the former. More specifically, we favour the possibility that the so-called COBE/DIRBE bar is the boxy/peanut bulge of a bar whose outer thin parts are the so-called Long bar. This possibility is in good agreement both with observations of external galaxies, with orbital structure theory and with simulations. We then analyse in detail the morphology and kinematics given by five representative Galactic potentials. Two have a Ferrers bar, two have a quadrupole bar and the last one a composite bar. We first consider only the COBE/DIRBE bar and then extend it to include the effect of the Long bar. We find that the large-scale structure given by the manifolds describes an inner ring, whose size is similar to the near and far 3-kpc arm, and an outer ring, whose properties resemble those of the Galactic Molecular Ring. We also analyse the kinematics of these two structures, under the different galactic potentials, and find they reproduce the main over-densities found in the galactic longitude-velocity CO diagram. Finally, we consider for what model parameters, the global morphology of the manifolds may reproduce the two outer spiral arms. We conclude that this would necessitate either more massive and more rapidly rotating bars, or including in the potential an extra component describing the spiral arms.

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.

In this series of papers, we propose a theory to explain the formation and properties of rings and spirals in barred galaxies. The building blocks of these structures are orbits guided by the manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In this paper we focus on a comparison of the morphology of observed and of theoretical spirals and rings and we also give some predictions for further comparisons. Our theory can account for spirals as well as both inner and outer rings. The model outer rings have the observed $R_1$, $R_1'$, $R_2$, $R_2'$ and $R_1R_2$ morphologies, including the dimples near the direction of the bar major axis. We explain why the vast majority of spirals in barred galaxies are two armed and trailing, and discuss what it would take for higher multiplicity arms to form. We show that the shapes of observed and theoretical spirals agree and we predict that stronger non-axisymmetric forcings at and somewhat beyond corotation will drive more open spirals. We compare the ratio of ring diameters in theory and in observations and predict that more elliptical rings will correspond to stronger forcings. We find that the model potential may influence strongly the numerical values of these ratios.

Two-armed, grand design spirals and inner and outer rings in barred galaxies can be due to orbits guided by the manifolds emanating from the vicinity of the L1 and L2 Lagrangian points, located at the ends of the bar. We first summarise the necessary theoretical background and in particular we describe the dynamics around the unstable equilibrium points in barred galaxy models, and the corresponding homoclinic and heteroclinic orbits. We then discuss two specific morphologies and the circulation of material within the corresponding manifolds. We also discuss the case where mass concentrations at the end of the bar can stabilise the L1 and L2 and the relevance of this work to the gas concentrations in spirals and rings.

In this paper we present building blocks which can explain the formation and properties both of spirals and of inner and outer rings in barred galaxies. We first briefly summarise the main results of the full theoretical description we have given elsewhere, presenting them in a more physical way, aimed to an understanding without the requirement of extended knowledge of dynamical systems or of orbital structure. We introduce in this manner the notion of manifolds, which can be thought of as tubes guiding the orbits. The dynamics of these manifolds can govern the properties of spirals and of inner and outer rings in barred galaxies. We find that the bar strength affects how unstable the L1 and L2 Lagrangian points are, the motion within the 5A5A5Amanifold tubes and the time necessary for particles in a manifold to make a complete turn around the galactic centre. We also show that the strength of the bar, or, to be more precise, of the non-axisymmetric forcing at and somewhat beyond the corotation region, determines the resulting morphology. Thus, less strong bars give rise to R1 rings or pseudorings, while stronger bars drive R2, R1R2 and spiral morphologies. We examine the morphology as a function of the main parameters of the bar and present descriptive two dimensional plots to that avail. We also derive how the manifold morphologies and properties are modified if the L1 and L2 Lagrangian points become stable. Finally, we discuss how dissipation affects the manifold properties and compare the manifolds in gas-like and in stellar cases. Comparison with observations, as well as clear predictions to be tested by observations will be given in an accompanying paper.

We study the role of the unstable equilibrium points in the transfer of matter in a galaxy using the potential of a rotating triaxial system. In particular, we study the neighbourhood of these points for energy levels and for main model parameters where the zero velocity curves just open and form a bottleneck in the region. For these energies, the transfer of matter from the inner to the outer parts and vice versa starts being possible. We study how the dynamics around the unstable equilibrium points is driven, by performing a partial normal form scheme and by computing the invariant manifolds of periodic orbits and quasi-periodic orbits using the reduced Hamiltonian. In particular, we compute some homoclinic and heteroclinic orbits playing a crucial role. Our results also show that in slow rotating and/or axisymmetric systems the hyperbolic character of the equilibrium points is cancelled, so that no transfer of matter is possible through the bottleneck.

We propose a theory to explain the formation of spiral arms and of all types of outer rings in barred galaxies, extending and applying the technique used in celestial mechanics to compute transfer orbits. Thus, our theory is based on the chaotic orbital motion driven by the invariant manifolds associated to the periodic orbits around the hyperbolic equilibrium points. In particular, spiral arms and outer rings are related to the presence of heteroclinic or homoclinic orbits. Thus, R1 rings are associated to the presence of heteroclinic orbits, while R1R2 rings are associated to the presence of homoclinic orbits. Spiral arms and R2 rings, however, appear when there exist neither heteroclinic nor homoclinic orbits. We examine the parameter space of three realistic, yet simple, barred galaxy models and discuss the formation of the different morphologies according to the properties of the galaxy model. The different morphologies arise from differences in the dynamical parameters of the galaxy.

We propose a theory to explain the formation of both spirals and rings in barred galaxies using a common dynamical framework. It is based on the orbital motion driven by the unstable equilibrium points of the rotating bar potential. Thus, spirals, rings and pseudo-rings are related to the invariant manifolds associated to the periodic orbits around these equilibrium points. We examine the parameter space of three barred galaxy models and discuss the formation of the different morphological structures according to the properties of the bar model. We also study the influence of the shape of the rotation curve in the outer parts, by making families of models with rising, flat or falling rotation curves in the outer parts. The differences between spiral and ringed structures arise from differences in the dynamical parameters of the host galaxies.

We propose a new theory to explain the formation of spiral arms and of all types of outer rings in barred galaxies. We have extended and applied the technique used in celestial mechanics to compute transfer orbits. Thus, our theory is based on the chaotic orbital motion driven by the invariant manifolds associated to the periodic orbits around the hyperbolic equilibrium points. In particular, spiral arms and outer rings are related to the presence of heteroclinic or homoclinic orbits. Thus, R1 rings are associated to the presence of heteroclinic orbits, while R1R2 rings are associated to the presence of homoclinic orbits. Spiral arms and R2 rings, however, appear when there exist neither heteroclinic nor homoclinic orbits. We examine the parameter space of three realistic, yet simple, barred galaxy models and discuss the formation of the different morphologies according to the properties of the galaxy model. The different morphologies arise from differences in the dynamical parameters of the galaxy.

In this and in a previous paper (Romero-Gomez et al. 2006) we propose a theory to explain the formation of both spirals and rings in barred galaxies using a common dynamical framework. It is based on the orbital motion driven by the unstable equilibrium points of the rotating bar potential. Thus, spirals, rings and pseudo-rings are related to the invariant manifolds associated to the periodic orbits around these equilibrium points. We examine the parameter space of three barred galaxy models and discuss the formation of the different morphological structures according to the properties of the bar model. We also study the influence of the shape of the rotation curve in the outer parts, by making families of models with rising, flat, or falling rotation curves in the outer parts. The differences between spiral and ringed structures arise from differences in the dynamical parameters of the host galaxies. The results presented here will be discussed and compared with observations in a forthcoming paper.

We propose a new theory for the formation of rR_1 ring structures, i.e. for ring structures with both an inner and an outer ring, the latter having the form of ``8''. We propose that these rings are formed by material from the stable and unstable invariant manifolds associated with the Lyapunov orbits around the equilibrium points of a barred galaxy. We discuss the shape and velocity structure of the rings thus formed and argue that they are in agreement with the observed properties of rR_1 structures.