We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value.

We use the second Gaia data release in combination with the catalog of Sanders & Das (2018) to dissect the Milky Way disc in phase-space and relative ages. We confirm and report the existence of multiple velocity moving groups at low azimuthal velocities and angular momenta, below Arcturus, regularly separated by $\sim~18-20\,\rm{km s^{-1}}$ in azimuthal velocity. Such features were predicted to exist more than ten years ago from the reaction of the Milky Way to a perturbation in the disc undergoing phase-mixing. These structures appear slightly younger than their phase-space surroundings, arguing against an extra-galactic origin. We also identify in relative age many of the classical ridges in the plane of azimuthal velocity vs. Galactocentric radius, traditionally associated with resonance features. These ridges are also younger than their phase-space surroundings in accordance with predictions from recent state-of-the-art cosmological hydrodynamical simulations of Milky Way-like galaxies. We study the response of dynamically young and old stellar disc populations to resonances from an analytic model of a large bar, which, remarkably, qualitatively reproduces the trends seen in the data. Our results re-inforce the idea that the Galactic disc is currently being shaped by both internal and external perturbations, and that, while absolute isochrone ages have to be taken with great care, exploring the dynamical structure of the disc in stellar ages, especially with future asteroseismic data, will provide much stronger constraints than metallicity/abundance trends alone.

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.

Data-driven methods play an increasingly important role in the field of astrophysics. In the context of large spectroscopic surveys of stars, data-driven methods are key in deducing physical parameters for millions of spectra in a short time. Convolutional Neural Networks (CNNs) enable us to connect observables (e.g., spectra, stellar magnitudes) to physical properties (atmospheric parameters, chemical abundances, or labels in general). We trained a CNN, adopting stellar atmospheric parameters and chemical abundances from APOGEE DR16 (resolution R~22500) data as training set labels. As input we used parts of the intermediate-resolution RAVE DR6 spectra (R~7500) overlapping with the APOGEE DR16 data as well as broad-band ALL_WISE and 2MASS photometry together with Gaia DR2 photometry and parallaxes. We derived precise atmospheric parameters Teff, log(g), [M/H] and chemical abundances of [Fe/H], [alpha/M], [Mg/Fe], [Si/Fe], [Al/Fe], and [Ni/Fe] for 420165 RAVE spectra. The precision typically amounts to 60K in Teff, 0.06 in log(g), and 0.02-0.04 dex for individual chemical abundances. Incorporating photometry and astrometry as additional constraints substantially improves the results in terms of accuracy and precision of the derived labels, as long as we operate in parts of the parameter space that is well covered by the training sample. Scientific validation confirms the robustness of the CNN results. We provide a catalog of CNN-trained atmospheric parameters and abundances along with their uncertainties for 420165 stars in the RAVE survey. CNN-based methods provide a powerful way to combine spectroscopic, photometric, and astrometric data, without applying any priors in the form of stellar evolutionary models. The developed procedure can extend the scientific output of RAVE spectra beyond DR6 to on-going and planned surveys such as Gaia RVS and 4MOST.

The density variations in thin stellar streams may encode important information on the nature of the dark matter. For instance, if dark matter aggregates into massive sub-halos, these perturbers are expected to scatter stars out of dynamically cold stellar streams, possibly leading to detectable gaps in those structures. Here we re-examine the density variations in the GD-1 stream, using Gaia DR2 astrometry, Pan-STARRS photometry, together with high precision radial velocities measured with the CFHT/ESPaDOnS and VLT/UVES instruments and complemented with public radial velocity catalogs. We show that after correcting for projection effects, the density profile exhibits high contrast periodic peaks, separated by $2.64\pm0.18$ kpc. An N-body simulation is presented that reproduces this striking morphology with simple epicyclic motion in a smooth Galactic potential. We also discuss the reliability of measuring density variations using ground-based photometric surveys, and for the particular case of GD-1 we highlight some of the artifacts present in the Gaia DR2 catalog along its track. Massive dark subhalos do not appear to be required to explain the density clumping along GD-1.

The second data release of the Gaia mission has revealed, in stellar velocity and action space, multiple ridges, the exact origin of which is still debated. Recently, we demonstrated that a large Galactic bar with pattern speed 39 km/s/kpc does create most of the observed ridges. Among those ridges, the Hercules moving group would then be associated to orbits trapped at the co-rotation resonance of the bar. Here we show that a distinctive prediction of such a model is that the angular momentum of Hercules at the Sun's radius must significantly decrease with increasing Galactocentric azimuth, i.e. when getting closer to the major axis of the bar. We show that such a dependence of the angular momentum of trapped orbits on the azimuth would on the other hand not happen close to the outer Lindblad resonance of a faster bar, unless the orbital distribution is still far from phase-mixed, namely for a bar perturbation younger than ~ 2 Gyr. Using Gaia DR2 and Bayesian distances from the StarHorse code, and tracing the average Galactocentric radial velocity as a function of angular momentum and azimuth, we show that the Hercules angular momentum changes significantly with azimuth as expected for the co-rotation resonance of a dynamically old large bar.

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ 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.

The second data release of the Gaia mission has revealed a very rich structure in local velocity space. In terms of in-plane motions, this rich structure is also seen as multiple ridges in the actions of the axisymmetric background potential of the Galaxy. These ridges are probably related to a combination of effects from ongoing phase-mixing and resonances from the spiral arms and the bar. We have recently developed a method to capture the behaviour of the stellar phase-space distribution function at a resonance, by re-expressing it in terms of a new set of canonical actions and angles variables valid in the resonant region. Here, by properly treating the distribution function at resonances, and by using a realistic model for a slowly rotating large Galactic bar with pattern speed 39 km/s/kpc, we show that no less than six ridges in local action space can be related to resonances with the bar. Two of these at low angular momentum correspond to the corotation resonance, and can be associated to the Hercules moving group in local velocity space. Another one at high angular momentum corresponds to the outer Lindblad resonance, and can tentatively be associated to the velocity structure seen as an arch at high azimuthal velocities in Gaia data. The other ridges are associated to the 3:1, 4:1 and 6:1 resonances. The latter can be associated to the so-called 'horn' of the local velocity distribution. While it is clear that effects from spiral arms and incomplete phase-mixing related to external perturbations also play a role in shaping the complex kinematics revealed by Gaia data, the present work demonstrates that, contrary to common misconceptions, the bar alone can create multiple prominent ridges in velocity and action space.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We determine the total dynamical density in the solar neighbourhood using the Tycho-Gaia Astrometric Solution (TGAS) catalogue. Astrometric measurements of proper motion and parallax of stars inform us of both the stellar number density distribution and the velocity distribution of stars close to the plane. Assuming equilibrium, these distributions are interrelated through the local dynamical density. For the first time, we do a full joint fit of the velocity and stellar number density distributions while accounting for the astrometric error of individual stars, in the framework of Bayesian Hierarchical Modelling. We use a sample of stars whose distance extends to approximately 160 pc from the Sun. We find a local matter density of $\rho_0 = 0.119^{+0.015}_{-0.012}~M_\odot$pc$^{-3}$, where the result is presented as the median to the posterior distribution, plus/minus its 16th and 84th percentiles. We find the Sun's position above the Galactic plane to be $z_\odot = 15.3^{+2.24}_{-2.16}$ pc, and the Sun's velocity perpendicular to the Galactic plane to be $W_\odot = 7.19^{+0.18}_{-0.18}$ km/s.

We measure the escape speed curve of the Milky Way based on the analysis of the velocity distribution of $\sim 2850$ counter-rotating halo stars from the Gaia DR2. The distances were estimated through the StarHorse code, and only stars with distance errors smaller than 10 per cent were used in the study. The escape speed curve is measured at Galactocentric radii ranging from $\sim 5$ kpc to $\sim 10.5$ kpc. The local Galactic escape at the Sun's position is estimated to be $v_\mathrm{e}(r_\odot)=580 \pm 63~\mathrm{km~s^{-1}}$, and it rises towards the Galactic center. Defined as the minimum speed required to reach three virial radii, our estimate of the escape speed as a function of radius implies, for a Navarro-Frenk-White profile and local circular velocity of $240~\mathrm{km~s^{-1}}$, a dark matter mass $M_{200}=1.28^{+0.68}_{-0.50} \times 10^{12}~M_\odot$ and a high concentration $c_{200}=11.09^{+2.94}_{-1.79}$. Assuming the mass-concentration relation of $\Lambda$CDM, we get $M_{200}=1.55_{-0.51}^{+0.64}\times 10^{12}~M_\odot$, $c_{200}=7.93_{-0.27}^{+0.33}$, for a local circular velocity of $228~\mathrm{km~s^{-1}}$.

The solar neighborhood contains disc stars that have recently crossed spiral arms in the Galaxy. We propose that boundaries in local velocity distributions separate stars that have recently crossed or been perturbed by a particular arm from those that haven't. Ridges in the stellar velocity distributions constructed from the second Gaia data release trace orbits that could have touched nearby spiral arms at apocentre or pericentre. The multiple ridges and arcs seen in local velocity distributions are consistent with the presence of multiple spiral features and different pattern speeds and imply that the outer Galaxy is flocculent rather than grand design.

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.

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.

We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.

The present-day response of a Galactic disc stellar population to a non-axisymmetric perturbation of the potential has previously been computed through perturbation theory within the phase-space coordinates of the unperturbed axisymmetric system. Such an Eulerian linearized treatment however leads to singularities at resonances, which prevent quantitative comparisons with data. Here, we manage to capture the behaviour of the distribution function (DF) at a resonance in a Lagrangian approach, by averaging the Hamiltonian over fast angle variables and re-expressing the DF in terms of a new set of canonical actions and angles variables valid in the resonant region. We then follow the prescription of Binney (2016), assigning to the resonant DF the time average along the orbits of the axisymmetric DF expressed in the new set of actions and angles. This boils down to phase-mixing the DF in terms of the new angles, such that the DF for trapped orbits only depends on the new set of actions. This opens the way to quantitatively fitting the effects of the bar and spirals to Gaia data in terms of distribution functions in action space.

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.

We report on the detection of a small overdensity of stars in velocity space with systematically higher Galactocentric rotation velocity than the Sun by about 20 km s$^{-1}$ in the $Gaia$ Data Release 1 Tycho-Gaia astrometric solution (TGAS) data. We find these fast Galactic rotators more clearly outside of the Solar radius, compared to inside of the Solar radius. In addition, the velocity of the fast Galactic rotators is independent of the Galactocentric distance up to $R-R_{\odot}\sim0.6$ kpc. Comparing with numerical models, we qualitatively discuss that a possible cause of this feature is the co-rotation resonance of the Perseus spiral arm, where the stars in peri-centre phase in the trailing side of the Perseus spiral arm experience an extended period of acceleration owing to the torque from the Perseus arm.

The length and pattern speed of the Milky Way bar are still controversial. Photometric and spectroscopic surveys of the inner Galaxy, as well as gas kinematics, favour a long and slowly rotating bar, with corotation around a Galactocentric radius of 6 kpc. On the other hand, the existence of the Hercules stream in local velocity space favours a short and fast bar with corotation around 4 kpc. This follows from the fact that the Hercules stream looks like a typical signature of the outer Lindblad resonance of the bar. As we showed recently, reconciling this local stream with a slow bar would need to find a yet unknown alternative explanation, based for instance on the effect of spiral arms. Here, by combining the TGAS catalogue of the Gaia DR1 with LAMOST radial velocities, we show that the position of Hercules in velocity space as a function of radius in the outer Galaxy indeed varies exactly as predicted by fast bar models with a pattern speed no less than 1.8 times the circular frequency at the Sun's position.

Both the three-dimensional density of red clump giants and the gas kinematics in the inner Galaxy indicate that the pattern speed of the Galactic bar could be much lower than previously estimated. Here, we show that such slow bar models are unable to reproduce the bimodality observed in local stellar velocity space. We do so by computing the response of stars in the Solar neighbourhood to the gravitational potential of slow and fast bars, in terms of their perturbed distribution function in action-angle space up to second order, as well as by identifying resonantly trapped orbits. We also check that the bimodality is unlikely to be produced through perturbations from spiral arms, and conclude that, contrary to gas kinematics, local stellar kinematics still favour a fast bar in the Milky Way, with a pattern speed of the order of almost twice (and no less than 1.8 times) the circular frequency at the Sun's position. This leaves open the question of the nature of the long flat extension of the bar in the Milky Way.

We investigate models of the Milky Way disc taking into account simultaneously the bar and a two-armed quasi-static spiral pattern. Away from major resonance overlaps, the mean stellar radial motions in the plane are essentially a linear superposition of the isolated effects of the bar and spirals. Thus, provided the bar is strong enough, even in the presence of spiral arms, these mean radial motions are predominantly affected by the Galactic bar for large scale velocity fluctuations. This is evident when comparing the peculiar line-of-sight velocity power spectrum of our coupled models with bar-only models. However, we show how forthcoming spectroscopic surveys could disentangle bar-only non-axisymmetric models of the Galaxy from models in which spiral arms have a significant amplitude. We also point out that overlaps of low-order resonances are sufficient to enhance stellar churning within the disc, even when the spirals amplitude is kept constant. Nevertheless, for churning to be truly non-local, stronger or (more likely) transient amplitudes would be needed: otherwise the disc is actually mostly unaffected by churning in the present models. Finally, regarding vertical breathing modes, the combined effect of the bar and spirals on vertical motions is a clear non-linear superposition of the isolated effects of both components, significantly superseding the linear superposition of modes produced by each perturber separately, thereby providing an additional effect to consider when analysing the observed breathing mode of the Galactic disc in the extended Solar neighbourhood.

Starting from an axisymmetric equilibrium distribution function (DF) in action space, representing a Milky Way thin disc stellar population, we use the linearized Boltzmann equation to explicitly compute the response to a three-dimensional spiral potential in terms of the perturbed DF. This DF, valid away from the main resonances, allows us to investigate a snapshot of the velocity distribution at any given point in three-dimensional configuration space. Moreover, the first order moments of the DF give rise to non-zero radial and vertical bulk flows -- namely breathing modes -- qualitatively similar to those recently observed in the extended Solar neighbourhood. We show that these analytically predicted mean stellar motions are in agreement with the outcome of test-particle simulations. Moreover, we estimate for the first time the reduction factor for the vertical bulk motions of a stellar population compared to the case of a cold fluid. Such an explicit expression for the full perturbed DF of a thin disc stellar population in the presence of spiral arms will be helpful in order to dynamically interpret the detailed information on the Milky Way disc stellar kinematics that will be provided by upcoming large astrometric and spectroscopic surveys of the Galaxy.

We investigate the kinematic signatures induced by spiral and bar structure in a set of simulations of Milky Way-sized spiral disc galaxies. The set includes test particle simulations that follow a quasi-stationary density wave-like scenario with rigidly rotating spiral arms, and $N$-body simulations that host a bar and transient, co-rotating spiral arms. From a location similar to that of the Sun, we calculate the radial, tangential and line-of-sight peculiar velocity fields of a patch of the disc and quantify the fluctuations by computing the power spectrum from a two-dimensional Fourier transform. We find that the peculiar velocity power spectrum of the simulation with a bar and transient, co-rotating spiral arms fits very well to that of APOGEE red clump star data, while the quasi-stationary density wave spiral model without a bar does not. We determine that the power spectrum is sensitive to the number of spiral arms, spiral arm pitch angle and position with respect to the spiral arm. However, it is necessary to go beyond the line of sight velocity field in order to distinguish fully between the various spiral models with this method. We compute the power spectrum for different regions of the spiral discs, and discuss the application of this analysis technique to external galaxies.

Evidence for non-zero mean stellar velocities in the direction perpendicular to the Galactic plane has been accumulating from various recent large spectroscopic surveys. Previous analytical and numerical work has shown that a "breathing mode" of the Galactic disc, similar to what is observed in the Solar vicinity, can be the natural consequence of a non-axisymmetric internal perturbation of the disc. Here we provide a general analytical framework, in the context of perturbation theory, allowing us to compute the vertical bulk motions generated by a single internal perturber (bar or spiral pattern). In the case of the Galactic bar, we show that these analytically predicted bulk motions are well in line with the outcome of a numerical simulation. The mean vertical motions induced by the Milky Way bar are small (mean velocity of less than 1 km/sec) and cannot be responsible alone for the observed breathing mode, but they are existing. Our analytical treatment is valid close to the plane for all the non-axisymmetric perturbations of the disc that can be described by small-amplitude Fourier modes. Further work should study how the coupling of multiple internal perturbers and external perturbers is affecting the present analytical results.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

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.

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.

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