We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

The Sagittarius stream is one of the best tools that we currently have to estimate the mass and shape of our Galaxy. However, assigning membership and obtaining the phase-space distribution of the stars that form the tails is quite challenging. Our goal is to produce a catalogue of RR Lyrae stars of Sagittarius and obtain an empiric measurement of the trends along the stream in sky position, distance and tangential velocities. We generate two initial samples from the Gaia DR2 RR Lyrae catalogue: one, selecting only the stars within \pm20deg of the orbital plane of Sagittarius (Strip) and the other, the result of applying the Pole Count Map (nGC3) algorithm. We then use the model-independent, deterministic method developed in this work to remove most of the contamination by detecting and isolating the stream in distance and proper motions. The output is two empiric catalogues: the Strip sample (higher-completeness, lower-purity) which contains 11 677 stars, and the nGC3 sample (higher-purity, lower-completeness) with 6 608 stars. We characterise the changes along the stream in all the available dimensions, the 5 astrometric ones plus the metallicity, covering more than 2pi rad in the sky and obtain new estimates for the apocentres and the mean [Fe/H] of the RR Lyrae population. Also, we show the first map of the two components of the tangential velocity, thanks to the combination of distances and proper motions. Finally, we detect the bifurcation in the leading arm and report no significant difference between the two branches, either in metallicity, kinematics or distance. We provide the largest sample of RR Lyrae candidates of Sagittarius, which can be used as an input for a spectroscopic follow-up or as a reference for the new generation of models of the stream through the interpolators in distance and velocity that we have constructed.

We aim to measure the proper motion along the Sagittarius stream that is the missing piece to determine its full 6D phase space coordinates. We conduct a blind search of over-densities in proper motion from Gaia DR2 in a broad region around the Sagittarius stream by applying wavelet transform techniques. We find that for most of the sky patches, the highest intensity peaks delineate the path of the Sagittarius stream. The 1500 peaks identified depict a continuous sequence spanning almost $2\pi$ in the sky, only obscured when the stream crosses the Galactic disk. Altogether, around $100\,000$ stars potentially belong to the stream as indicated by a coarse inspection of the colour-magnitude diagrams. From these stars, we determine the proper motion along the Sagittarius stream, making it the proper motion sequence with the largest span and continuity ever measured for a stream. A first comparison with existing N-body models of the stream reveals some discrepancies, especially near the pericentre of the trailing arm and an overestimation of the total proper motion for the leading arm. Our study can be the starting point for determining the variation of the population of stars along the stream, the distance to the stream with red clump stars, and the solar motion. It will also allow a much better measurement of the Milky Way potential.

We analyzed the velocity space of the thin and thick-disk Gaia white dwarf population within 100 pc looking for signatures of the Hercules stellar stream. We aimed to identify those objects belonging to the Hercules stream and, by taking advantage of white dwarf stars as reliable cosmochronometers, to derive a first age distribution. We applied a kernel density estimation to the $UV$ velocity space of white dwarfs. For the region where a clear overdensity of stars was found, we created a 5-D space of dynamic variables. We applied a hierarchichal clustering method, HDBSCAN, to this 5-D space, identifying those white dwarfs that share similar kinematic characteristics. Finally, under general assumptions and from their photometric properties, we derived an age estimate for each object. The Hercules stream was firstly revealed as an overdensity in the $UV$ velocity space of the thick-disk white dwarf population. Three substreams were then found: Hercules $a$ and Hercules $b$, formed by thick-disk stars with an age distribution peaked $4\,$Gyr in the past and extended to very old ages; and Hercules $c$, with a ratio of 65:35 thin:thick stars and a more uniform age distribution younger than 10 Gyr.

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.

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.

We develop a new theoretical framework to generate Besan\c{c}on Galaxy Model fast approximate simulations (BGM FASt) to address fundamental questions of the Galactic structure and evolution performing multi-parameter inference. As a first application of our strategy we simultaneously infer the IMF, the star formation history and the stellar mass density in the Solar Neighbourhood. The BGM FASt strategy is based on a reweighing scheme, that uses a specific pre-sampled simulation, and on the assumption that the distribution function of the generated stars in the Galaxy can be described by an analytical expression. To validate BGM FASt we execute a set of tests. Finally, we use BGM FASt with an approximate Bayesian computation algorithm to obtain the posterior PDF of the inferred parameters, by comparing synthetic versus Tycho-2 colour-magnitude diagrams. Results: The validation shows a very good agreement between BGM FASt and the standard BGM, with BGM FASt being $\approx 10^4$ times faster. By analysing Tycho-2 data we obtain a thin disc star formation history decreasing in time and a present rate of $1.2 \pm 0.2 M_\odot/yr$. The resulting total stellar mass density in the Solar Neighbourhood is $0.051_{-0.005}^{+0.002} M_\odot/pc^3$ and the local dark matter density is $0.012 \pm 0.001 M_\odot/pc^3$. For the composite IMF we obtain a slope of $\alpha_2={2.1}_{-0.3}^{+0.1}$ in the mass range between $0.5 M_\odot$ and $1.53M_\odot$. The results of the slope at the high mass range are trustable up to $4M_\odot$ and highly depend on the choice of the extinction map (obtaining $\alpha_3={2.9}_{-0.2}^{+0.2}$ and $\alpha_3={3.7}_{-0.2}^{+0.2}$ respectively, for two different extinction maps). Systematic uncertainties are not included. Conclusions: The good performance of BGM FASt demonstrates that it is a very valuable tool to perform multi-parameter inference using Gaia data releases.

The evolution of the Milky Way disk, which contains most of the stars in the Galaxy, is affected by several phenomena. For example, the bar and the spiral arms of the Milky Way induce radial migration of stars and can trap or scatter stars close to orbital resonances. External perturbations from satellite galaxies can also have a role, causing dynamical heating of the Galaxy, ring-like structures in the disk and correlations between different components of the stellar velocity. These perturbations can also cause 'phase wrapping' signatures in the disk, such as arched velocity structures in the motions of stars in the Galactic plane. Some manifestations of these dynamical processes have already been detected, including kinematic substructure in samples of nearby stars, density asymmetries and velocities across the Galactic disk that differ from the axisymmetric and equilibrium expectations, especially in the vertical direction, and signatures of incomplete phase mixing in the disk. Here we report an analysis of the motions of six million stars in the Milky Way disk. We show that the phase-space distribution contains different substructures with various morphologies, such as snail shells and ridges, when spatial and velocity coordinates are combined. We infer that the disk must have been perturbed between 300 million and 900 million years ago, consistent with estimates of the previous pericentric passage of the Sagittarius dwarf galaxy. Our findings show that the Galactic disk is dynamically young and that modelling it as time-independent and axisymmetric is incorrect.

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.

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.

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of 150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions.

Gaia DR2 has delivered full-sky 6-D measurements for millions of stars, and the quest to understand the dynamics of our Galaxy has entered a new phase. Our aim is to reveal and characterize the kinematic sub-structure of the different Galactic neighbourhoods, to form a picture of their spatial evolution that can be used to infer the Galactic potential, its evolution and its components. We take ~5 million stars in the Galactic disk from the Gaia DR2 catalogue and build the velocity distribution of many different Galactic Neighbourhoods distributed along 5 kpc in Galactic radius and azimuth. We decompose their distribution of stars in the V_R-V_phi plane with the wavelet transformation and asses the statistical significance of the structures found. We detect many kinematic sub-structures (arches and more rounded groups) that diminish their azimuthal velocity as a function of Galactic radius in a continuous way, connecting volumes up to 3 kpc apart in some cases. The decrease rate is, on average, of ~23 km/s/kpc. In azimuth, the kinematic sub-structures present much smaller variations. We also observe a duality in their behaviour: some conserve their vertical angular momentum with radius (e.g., Hercules), while some seem to have nearly constant kinetic energy (e.g., Sirius). These two trends are consistent with the approximate predictions of resonances and of phase mixing, respectively. Besides, the overall spatial evolution of Hercules is consistent with being related to the Outer Lindblad Resonance of the Bar. We also detect structures without apparent counterpart in the vicinity of the Sun. The various trends observed and their continuity with radius and azimuth allows for future work to deeply explore the parameter space of the models. Also, the characterization of extrasolar moving groups opens the opportunity to expand our understanding of the Galaxy beyond the Solar Neighbourhood.

Before the publication of the Gaia DR2 we confirmed with RAVE and TGAS an observation recently made with the GALAH survey by Quillen ey al. concerning the Coma Berenices moving group in the Solar neighbourhood, namely that it is only present at negative Galactic latitudes. This allowed us to show that it is coherent in vertical velocity, providing a first evidence for incomplete vertical phase-mixing. We estimated for the first time from dynamical arguments that the moving group must have formed at most ~ 1.5 Gyr ago, and related this to a pericentric passage of the Sagittarius dwarf satellite galaxy. The present note is a rewritten version of the original arXiv post on this result now also including a confirmation of our finding with Gaia DR2.

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.

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)

We explore the connections between stellar age, chemistry, and kinematics across a Galactocentric distance of $7.5 < R\,(\mathrm{kpc}) < 9.0$, using a sample of $\sim 12\,000$ intermediate-mass (FGK) turnoff stars observed with the RAdial Velocity Experiment (RAVE) survey. The kinematics of this sample are determined using radial velocity measurements from RAVE, and parallax and proper motion measurements from the Tycho-Gaia Astrometric Solution (TGAS). In addition, ages for RAVE stars are determined using a Bayesian method, taking TGAS parallaxes as a prior. We divide our sample into young ($0 < \tau < 3$ Gyr) and old ($8 < \tau < 13$ Gyr) populations, and then consider different metallicity bins for each of these age groups. We find significant differences in kinematic trends of young and old, metal-poor and metal-rich, stellar populations. In particular, we find a strong metallicity dependence in the mean Galactocentric radial velocity as a function of radius ($\partial {V_{\rm R}}/\partial R$) for young stars, with metal-rich stars having a much steeper gradient than metal-poor stars. For $\partial {V_{\phi}}/\partial R$, young, metal-rich stars significantly lag the LSR with a slightly positive gradient, while metal-poor stars show a negative gradient above the LSR. We interpret these findings as correlations between metallicity and the relative contributions of the non-axisymmetries in the Galactic gravitational potential (the spiral arms and the bar) to perturb stellar orbits.

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.

We present DECam imaging for the peculiar Galactic globular clusters NGC1851, NGC1904 (M79), NGC2298 and NGC2808. Our deep photometry reveals that all the clusters have an important contribution of stars beyond their King tidal radii and present tails with different morphologies. We have also explored the surroundings of the clusters where the presence of the Canis Major overdensity and/or the low Galactic latitude Monoceros ring at d~8kpc is evident. A second stellar system is found at d~17kpc and spans at least 18deg x 15deg in the sky. As one of the possible scenarios to explain that feature, we propose that the unveiled system is part of Monoceros explained as a density wave moving towards the outer Milky Way. Alternatively, the unveiled system might be connected with other known halo substructures or associated with the progenitor dwarf galaxy of NGC1851 and NGC1904, which are widely considered accreted globular clusters.

We use the Gaia-TGAS data to compare the transverse velocities in Galactic longitude (coming from proper motions and parallaxes) in the Milky Way disk for negative and positive longitudes as a function of distance. The transverse velocities are strongly asymmetric and deviate significantly from the expectations for an axisymmetric Galaxy. The value and sign of the asymmetry changes at spatial scales of several tens of degrees in Galactic longitude and about 0.5 kpc in distance. The asymmetry is statistically significant at 95% confidence level for 57% of the region probed, which extends up to ~1.2 kpc. A percentage of 24% of the region studied shows absolute differences at this confidence level larger than 5 km/s and 7% larger than 10 km/s. The asymmetry pattern shows mild variations in the vertical direction and with stellar type. A first qualitative comparison with spiral arm models indicates that the arms are unlikely to be the main source of the asymmetry. We briefly discuss alternative origins. This is the first time that global all-sky asymmetries are detected in the Milky Way kinematics, beyond the local neighbourhood, and with a purely astrometric sample.

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.

We explore the correlations between velocity and metallicity and the possible distinct chemical signatures of the velocity over-densities of the local Galactic neighbourhood. We use the large spectroscopic survey RAVE and the Geneva Copenhagen Survey. We compare the metallicity distribution of regions in the velocity plane ($v_R,v_\phi$) with that of their symmetric counterparts ($-v_R,v_\phi$). We expect similar metallicity distributions if there are no tracers of a sub-population (e.g., a dispersed cluster, accreted stars), if the disk of the Galaxy is axisymmetric, and if the orbital effects of the spiral arms and the bar are weak. We find that the metallicity-velocity space of the solar neighbourhood is highly patterned. A large fraction of the velocity plane shows differences in the metallicity distribution when comparing symmetric $v_R$ regions. The typical differences in the median metallicity are of $0.05$ dex with a statistical significance of at least $95\%$, and with values up to $0.6$ dex. For low azimuthal velocity $v_\phi$, stars moving outwards in the Galaxy have on average higher metallicity than those moving inwards. These include stars in the Hercules and Hyades moving groups and other velocity branch-like structures. For higher $v_\phi$, the stars moving inwards have higher metallicity than those moving outwards. The most likely interpretation of the metallicity asymmetry is that it is due to the orbital effects of the bar and the radial metallicity gradient of the disk. We present a simulation that supports this idea. We have also discovered a positive gradient in $v_\phi$ with respect to metallicity at high metallicities, apart from the two known positive and negative gradients for the thick and thin disks, respectively.

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.

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.

Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795\AA) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities from 520,781 spectra of 457,588 unique stars are presented, of which 255,922 stellar observations have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. For our main DR5 catalog, stellar parameters (effective temperature, surface gravity, and overall metallicity) are computed using the RAVE DR4 stellar pipeline, but calibrated using recent K2 Campaign 1 seismic gravities and Gaia benchmark stars, as well as results obtained from high-resolution studies. Also included are temperatures from the Infrared Flux Method, and we provide a catalogue of red giant stars in the dereddened color $(J-Ks)_0$ interval (0.50,0.85) for which the gravities were calibrated based only on seismology. Further data products for sub-samples of the RAVE stars include individual abundances for Mg, Al, Si, Ca, Ti, Fe, and Ni, and distances found using isochrones. Each RAVE spectrum is complemented by an error spectrum, which has been used to determine uncertainties on the parameters. The data can be accessed via the RAVE Web site or the Vizier database.

We model the effects of the spiral arms of the Milky Way on the disk stellar kinematics in the Gaia observable space. We also estimate the Gaia capabilities of detecting the predicted signatures. We use both controlled orbital integrations in analytic potentials and self-consistent simulations. We introduce a new strategy, which consists of comparing the stellar kinematics of symmetric Galactic longitudes (+l and -l), in particular the median transverse velocity (from parallaxes and proper motions). This approach does not require the assumption of an axisymmetric model. The typical differences between the transverse velocity in symmetric longitudes in the models are of ~2 km/s, but can be larger than 10 km/s for certain longitudes and distances. The kinematic differences for +l and -l show trends that depend on the properties of spiral arms. Thus, this method can be used to quantify the importance of the effects of spiral arms on the orbits of stars in the different disk regions, and to constrain the location of the arms, main resonances and, thus, pattern speed. Moreover, the method allows us to test the dynamical nature of the spiral structure (e.g. grand design versus transient multiple arms). We estimate the number of stars of each spectral type that Gaia will observe in representative Galactic longitudes, their errors in distance and transverse velocity, and the error in computing the median velocity as a function of distance. We will be able to measure the median transverse velocity with Gaia data, with precision smaller than ~1 km/s up to distances of ~4-6 kpc for certain giant stars, and up to ~2-4 kpc and better kinematic precision (<= 0.5 km\s) for certain sub-giants and dwarfs. These are enough to measure the typical signatures seen in the models. We also show that a similar strategy can be used with line-of-sight velocities from other upcoming spectroscopic surveys.(Abridged)

We use self-consistent numerical simulations of the evolution and disruption of the Sun's birth cluster in the Milky Way potential to investigate the present-day phase space distribution of the Sun's siblings. The simulations include the gravitational N-body forces within the cluster and the effects of stellar evolution on the cluster population. In addition the gravitational forces due to the Milky Way potential are accounted for in a self-consistent manner. Our aim is to understand how the astrometric and radial velocity data from the Gaia mission can be used to pre-select solar sibling candidates. We vary the initial conditions of the Sun's birth cluster, as well as the parameters of the Galactic potential. We show that the disruption time-scales of the cluster are insensitive to the details of the non-axisymmetric components of the Milky Way model and we make predictions, averaged over the different simulated possibilities, about the number of solar siblings that should appear in surveys such as Gaia or GALAH. We find a large variety of present-day phase space distributions of solar siblings, which depend on the cluster initial conditions and the Milky Way model parameters. We show that nevertheless robust predictions can be made about the location of the solar siblings in the space of parallaxes ($\varpi$), proper motions ($\mu$) and radial velocities ($V_\mathrm{r}$). By calculating the ratio of the number of simulated solar siblings to that of the number of stars in a model Galactic disk, we find that this ratio is above 0.5 in the region given by: $\varpi \geq 5$mas, $4 \leq \mu \leq 6$masyr$^{-1}$, and $-2\leq V_\mathrm{r} \leq 0$kms$^{-1}$. Selecting stars from this region should increase the probability of success in identifying solar siblings through follow up observations [Abridged].

We present a method to identify Ultra Faint Dwarf Galaxy (UFDG) candidates in the halo of the Milky Way using the future Gaia catalogue and we explore its detection limits and completeness. The method is based on the Wavelet Transform and searches for over-densities in the combined space of sky coordinates and proper motions, using kinematics in the search for the first time. We test the method with a Gaia mock catalogue that has the Gaia Universe Model Snapshot (GUMS) as a background, and use a library of around 30 000 UFDGs simulated as Plummer spheres with a single stellar population. For the UFDGs we use a wide range of structural and orbital parameters that go beyond the range spanned by real systems, where some UFDGs may remain undetected. We characterize the detection limits as function of the number of observable stars by Gaia in the UFDGs with respect to that of the background and their apparent sizes in the sky and proper motion planes. We find that the addition of proper motions in the search improves considerably the detections compared to a photometric survey at the same magnitude limit. Our experiments suggest that Gaia will be able to detect UFDGs that are similar to some of the known UFDGs even if the limit of Gaia is around 2 magnitudes brighter than that of SDSS, with the advantage of having a full-sky catalogue. We also see that Gaia could even find some UFDGs that have lower surface brightness than the SDSS limit.

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE dataset. We find clear asymmetries in the $v_R$-$v_\phi$ velocity distributions of thin and thick disk stars: here are more stars moving radially outwards for low azimuthal velocities and more radially inwards for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick disk stars should no longer be described by Schwarzschild's, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully.

The standard errors of the end-of-mission Gaia astrometry have been re-assessed after conclusion of the in-orbit commissioning phase of the mission. An analytical relation is provided for the parallax standard error as function of Gaia G magnitude (and V-I colour) which supersedes the pre-launch relation provided in de Bruijne (2012).

In this first paper we simulate the population of disc Red Clump stars to be observed by Gaia. We generate a set of test particles and we evolve it in a 3D barred Milky Way like galactic potential. We assign physical properties of the Red Clump trace population and a realistic 3D interstellar extinction model. We add Gaia observational constraints and an error model according to the pre-commissioning scientific performance assessments. We present and analyse two mock catalogues, offered to the community, that are an excellent test bed for testing tools being developed for the future scientific exploitation of Gaia data. The first catalogue contains stars up to Gaia G 20, while the second is the subset containing Gaia radial velocity data with a maximum error of sigmaVr=10 kms. Here we present first attempts to characterise the density structure of the Galactic bar in the Gaia space of observables. The Gaia large errors in parallax and the high interstellar extinction in the inner parts of the Galactic disc prevent us to model the bar overdensity. This result suggests the need to combine Gaia and IR data to undertake such studies. We find that IR photometric distances for this Gaia sample allow us to recover the Galactic bar orientation angle with an accuracy of approximately 5 degrees.

Aims: We investigate whether the cylindrical (galactocentric) radial velocity gradient of ~ -3 km/s/kpc, directed radially from the Galactic center and recently observed in the stars of the solar neighborhood with the RAVE survey, can be explained by the resonant effects of the bar near the solar neighborhood. Methods: We compared the results of test particle simulations of the Milky Way with a potential that includes a rotating bar with observations from the RAVE survey. To this end we applied the RAVE selection function to the simulations and convolved these with the characteristic RAVE errors. We explored different "solar neighborhoods" in the simulations, as well as different bar models Results: We find that the bar induces a negative radial velocity gradient at every height from the Galactic plane, outside the outer Lindblad resonance and for angles from the long axis of the bar compatible with the current estimates. The selection function and errors do not wash away the gradient, but often make it steeper, especially near the Galactic plane, because this is where the RAVE survey is less radially extended. No gradient in the vertical velocity ispresent in our simulations, from which we may conclude that this cannot be induced by the bar.

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.

Non-axisymmetries in the Galactic potential (spiral arms and bar) induce kinematic groups such as the Hercules stream. Assuming that Hercules is caused by the effects of the Outer Lindblad Resonance of the Galactic bar, we model analytically its properties as a function of position in the Galaxy and its dependence on the bar's pattern speed and orientation. Using data from the RAVE survey we find that the azimuthal velocity of the Hercules structure decreases as a function of Galactocentric radius, in a manner consistent with our analytical model. This allows us to obtain new estimates of the parameters of the Milky Way's bar. The combined likelihood function of the bar's pattern speed and angle has its maximum for a pattern speed of $\Omega_b=(1.89 \pm 0.08) \times \Omega_0$ where $\Omega_0$ is the local circular frequency. Assuming a Solar radius of 8.05 kpc and a local circular velocity of 238 km/s, this corresponds to $\Omega_b=56\pm 2$ km/s/kpc. On the other hand, the bar's orientation $\phi_b$ cannot be constrained with the available data. In fact, the likelihood function shows that a tight correlation exists between the pattern speed and the orientation, implying that a better description of our best fit results is given by the linear relation $\Omega_b/\Omega_0=1.91 + 0.0044[\phi_b(\deg)-48]$, with standard deviation of 0.02. For example, for an angle of $\phi_b=30\deg$ the pattern speed is $54.0\pm0.5$ km/s/kpc. These results are not very sensitive to the other Galactic parameters such as the circular velocity curve or the peculiar motion of the Sun, and are robust to biases in distance.

Aims: To study the imprints of a rotating bar on the kinematics of stars in the thin and thick disks throughout the Galaxy. Methods: We perform test particle numerical simulations of the thin and thick disks in a 3D Galactic potential that includes a halo, a bulge, thin and thick disks, and a Ferrers bar. We analyze the resulting velocity distributions of populations corresponding to both disks, for different positions in the Galaxy and for different structural parameters of the bar. Results: We find that the velocity distributions of the disks are affected by the bar, and that strong transient effects are present for approximately 10 bar rotations after this is introduced adiabatically. On long (more realistic) timescales, the effects of the bar are strong on the kinematics of thin disk stars, and weaker on those in the thick disk, but in any case significant. Furthermore, we find that it is possible to trace the imprints of the bar also vertically and at least up to z ~ 1 kpc for the thin disk and z ~ 2 kpc for the thick disk.

We analyse the kinematics of disc stars observed by the RAVE survey in and beyond the Solar neighbourhood.We detect significant overdensities in the velocity distributions using a technique based on the wavelet transform.We find that the main local kinematic groups are large scale features, surviving at least up to ~1 kpc from the Sun in the direction of anti-rotation, and also at ~700 pc below the Galactic plane.We also find that for regions located at different radii than the Sun, the known groups appear shifted in the velocity plane. For example, the Hercules group has a larger azimuthal velocity for regions inside the Solar circle and a lower value outside. We have also discovered a new group at (U, V) = (92,-22) km/s in the Solar neighbourhood and confirmed the significance of other previously found groups. Some of these trends detected for the first time are consistent with dynamical models of the effects of the bar and the spiral arms. More modelling is required to definitively characterise the non-axisymmetric components of our Galaxy using these groups.

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 present a method to constrain the potential of the non-axisymmetric components of the Galaxy using the kinematics of stars in the solar neighborhood. The basic premise is that dynamical substructures in phase-space (i.e. due to the bar and/or spiral arms) are associated with families of periodic or irregular orbits, which may be easily identified in orbital frequency space. We use the "observed" positions and velocities of stars as initial conditions for orbital integrations in a variety of gravitational potentials. We then compute their characteristic frequencies, and study the structure present in the frequency maps. We find that the distribution of dynamical substructures in velocity- and frequency-space is best preserved when the integrations are performed in the "true" gravitational potential.

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.

We study the phase space available to the local stellar distribution using a Galactic potential consistent with several recent observational constraints. We find that the induced phase space structure has several observable consequences. The spiral arm contribution to the kinematic structure in the solar neighborhood may be as important as the one produced by the Galactic bar. We suggest that some of the stellar kinematic groups in the solar neighborhood, like the Hercules structure and the kinematic branches, can be created by the dynamical resonances of self-gravitating spiral arms and not exclusively by the Galactic bar. A structure coincident with the Arcturus kinematic group is developed when a hot stellar disk population is considered, which introduces a new perspective on the interpretation of its extragalactic origin. A bar-related resonant mechanism can modify this kinematic structure.We show that particles in the dark matter disk-like structure predicted by recent LCDM galaxy formation experiments, with similar kinematics to the thick disk, are affected by the same resonances, developing phase space structures or dark kinematic groups that are independent of the Galaxy assembly history and substructure abundance. We discuss the possibility of using the stellar phase space groups as constraints to non-axisymmetric models of the Milky Way structure.

Context. Recent studies have suggested that moving groups have a dynamic or "resonant" origin. Under this hypothesis, these kinematic structures become a powerful tool for studying the large-scale structure and dynamics of the Milky Way. Aims. We aim to characterize these structures in the U-V-age-[Fe/H] space and establish observational constraints that will allow us to study their origin and evolution. Methods. We apply multiscale techniques -wavelet denoising (WD)- to an extensive compendium of more than 24000 stars in the solar neighbourhood with the best available astrometric, photometric and spectroscopic data. Results. We confirm that the dominant structures in the U-V plane are the branches of Sirius, Coma Berenices, Hyades-Pleiades and Hercules, which are nearly equidistant in this kinematic plane and show a negative slope. The abrupt drops in the velocity distribution are characterized. We find a certain dependence of these kinematic structures on Galactic position with a significant change of contrast among substructures inside the branches. A large spread of ages is observed for all branches. The Hercules branch is detected in all subsamples with ages older than ~ 2 Gyr and the set of the other three branches is well established for stars > 400 Myr. The age-metallicity relation of each branch is examined and the relation between kinematics and metallicity is studied. Conclusions. Not all of these observational constraints are successfully explained by the recent models proposed for the formation of such kinematic structures. Simulations incorporating stellar ages and metallicities are essential for future studies. The comparison of the observed and simulated distributions obtained by WD will provide a physical interpretation of the existence of the branches in terms of local or large-scale dynamics. [Abridged]