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

We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar parameter pipeline on the final resulting catalogue, which in principle limits the parameter space over which our selection function is valid. Our final selection function is the product of the completeness fraction and the selection function of the pipeline. We then test if the application of the selection function introduces biases in the derived parameters. To do this, we compare a parent mock catalogue generated using Galaxia with a mock-RAVE catalogue where the selection function of RAVE has been applied. We conclude that for stars brighter than I = 12, between $4000 \rm K < T_{\rm eff} < 8000 \rm K$ and $0.5 < \rm{log}\,g < 5.0$, RAVE is kinematically and chemically unbiased with respect to expectations from Galaxia.

In a companion paper, we have presented the combined asteroseismic-spectroscopic dataset obtained from CoRoT lightcurves and APOGEE infra-red spectra for 678 solar-like oscillating red giants in two fields of the Galactic disc (CoRoGEE). We have measured chemical abundance patterns, distances, and ages of these field stars which are spread over a large radial range of the Milky Way's disc. Here we show how to simulate this dataset using a chemodynamical Galaxy model. We also demonstrate how the observation procedure influences the accuracy of our estimated ages.

Starting from the hypothesis that the Galaxy's dark halo responded adiabatically to the infall of baryons, we have constructed a self-consistent dynamical model of the Galaxy that satisfies a large number of observations, including measurements of gas terminal velocities and masers, the kinematics of a 180,000 giant stars from the RAVE survey, and star count data from the SDSS. The stellar disc and the dark halo are both specified by distribution functions (DFs) of the action integrals. The model is obtained by extending the work of Piffl Penoyre & Binney (2015} from the construction of a single model to a systematic search of model space. Whereas the model of Piffl et al violated constraints on the terminal-velocity curve, our model respects these constraints by adopting a long scale length R_d=3.66 kpc for the thin and thick discs. The model is, however, inconsistent with the measured optical depth for microlensing of bulge stars because it attributes too large a fraction of the density at R <~ 3 kpc to dark matter rather than stars. Moreover, it now seems likely that the thick disc's scale-length is significantly shorter than the model implies. Shortening this scale-length would cause the constraints from the rotation curve to be violated anew. We conclude that we can now rule out adiabatic compression of our Galaxy's dark halo.

We present a new method to construct fully self-consistent equilibrium models of multi-component disc galaxies similar to the Milky Way. We define distribution functions for the stellar disc and dark halo that depend on phase space position only through action coordinates. We then use an iterative approach to find the corresponding gravitational potential. We study the adiabatic response of the initially spherical dark halo to the introduction of the baryonic component and find that the halo flattens in its inner regions with final minor-major axis ratios $q$ = 0.75 - 0.95. The extent of the flattening depends on the velocity structure of the halo particles with radially biased models exhibiting a stronger response. In this latter case, which is according to cosmological simulations the most likely one, the new density structure resembles a "dark disc" superimposed on a spherical halo. We discuss the implications of these results for our recent estimate of the local dark matter density. The velocity distribution of the dark-matter particles near the Sun is very non-Gaussian. All three principal velocity dispersions are boosted as the halo contracts, and at low velocities a plateau develops in the distribution of $v_z$. For models similar to a state-of-the-art Galaxy model we find velocity dispersions around 155 km s$^{-1}$ for $v_z$ and the tangential velocity, $v_\varphi$, and 140 - 175 km s$^{-1}$ for the in-plane radial velocity, $v_R$, depending on the anisotropy of the model.

We use the kinematics of $\sim200\,000$ giant stars that lie within $\sim 1.5$ kpc of the plane to measure the vertical profile of mass density near the Sun. We find that the dark mass contained within the isodensity surface of the dark halo that passes through the Sun ($(6\pm0.9)\times10^{10}\,\mathrm{M_\odot}$), and the surface density within $0.9$ kpc of the plane ($(69\pm10)\,\mathrm{M_\odot\,pc^{-2}}$) are almost independent of the (oblate) halo's axis ratio $q$. If the halo is spherical, 46 per cent of the radial force on the Sun is provided by baryons, and only 4.3 per cent of the Galaxy's mass is baryonic. If the halo is flattened, the baryons contribute even less strongly to the local radial force and to the Galaxy's mass. The dark-matter density at the location of the Sun is $0.0126\,q^{-0.89}\,\mathrm{M_\odot\,pc^{-3}}=0.48\,q^{-0.89}\,\mathrm{GeV\,cm^{-3}}$. When combined with other literature results we find hints for a mildly oblate dark halo with $q \simeq 0.8$. Our value for the dark mass within the solar radius is larger than that predicted by cosmological dark-matter-only simulations but in good agreement with simulations once the effects of baryonic infall are taken into account. Our mass models consist of three double-exponential discs, an oblate bulge and a Navarro-Frenk-White dark-matter halo, and we model the dynamics of the RAVE stars in the corresponding gravitational fields by finding distribution functions $f(\mathbf{J})$ that depend on three action integrals. Statistical errors are completely swamped by systematic uncertainties, the most important of which are the distance to the stars in the photometric and spectroscopic samples and the solar distance to the Galactic centre. Systematics other than the flattening of the dark halo yield overall uncertainties $\sim 15$ per cent.

Stellar population studies of globular clusters have suggested that the brightest clusters in the Galaxy might actually be the remnant nuclei of dwarf spheroidal galaxies. If the present Galactic globular clusters formed within larger stellar systems, they are likely surrounded by extra-tidal halos and/or tails made up of stars that were tidally stripped from their parent systems. The stellar surroundings around globular clusters are therefore one of the best places to look for the remnants of an ancient dwarf galaxy. Here an attempt is made to search for tidal debris around the supernovae enriched globular clusters M22 and NGC 1851 as well as the kinematically unique cluster NGC 3201. The stellar parameters from the Radial Velocity Experiment (RAVE) are used to identify stars with RAVE metallicities, radial velocities and elemental-abundances consistent with the abundance patterns and properties of the stars in M22, NGC 1851 and NGC 3201. The discovery of RAVE stars that may be associated with M22 and NGC 1851 are reported, some of which are at projected distances of ~10 degrees away from the core of these clusters. Numerous RAVE stars associated with NGC 3201 suggest that either the tidal radius of this cluster is underestimated, or that there are some unbound stars extending a few arc minutes from the edge of the cluster's radius. No further extra-tidal stars associated with NGC 3201 could be identified. The bright magnitudes of the RAVE stars make them easy targets for high resolution follow-up observations, allowing an eventual further chemical tagging to solidify (or exclude) stars outside the tidal radius of the cluster as tidal debris. In both our radial velocity histograms of the regions surrounding NGC 1851 and NGC 3201, a peak of stars at 230 km/s is seen, consistent with extended tidal debris from omega Centauri.

We provide new constraints on the chemo-dynamical models of the Milky Way by measuring the radial and vertical chemical gradients for the elements Mg, Al, Si, Ti, and Fe in the Galactic disc and the gradient variations as a function of the distance from the Galactic plane ($Z$). We selected a sample of giant stars from the RAVE database using the gravity criterium 1.7$<$log g$<$2.8. We created a RAVE mock sample with the Galaxia code based on the Besan\c con model and selected a corresponding mock sample to compare the model with the observed data. We measured the radial gradients and the vertical gradients as a function of the distance from the Galactic plane $Z$ to study their variation across the Galactic disc. The RAVE sample exhibits a negative radial gradient of $d[Fe/H]/dR=-0.054$ dex kpc$^{-1}$ close to the Galactic plane ($|Z|<0.4$ kpc) that becomes flatter for larger $|Z|$. Other elements follow the same trend although with some variations from element to element. The mock sample has radial gradients in fair agreement with the observed data. The variation of the gradients with $Z$ shows that the Fe radial gradient of the RAVE sample has little change in the range $|Z|\lesssim0.6$ kpc and then flattens. The iron vertical gradient of the RAVE sample is slightly negative close to the Galactic plane and steepens with $|Z|$. The mock sample exhibits an iron vertical gradient that is always steeper than the RAVE sample. The mock sample also shows an excess of metal-poor stars in the [Fe/H] distributions with respect to the observed data. These discrepancies can be reduced by decreasing the number of thick disc stars and increasing their average metallicity in the Besan\c con model.

We construct new estimates on the Galactic escape speed at various Galactocentric radii using the latest data release of the Radial Velocity Experiment (RAVE DR4). Compared to previous studies we have a database larger by a factor of 10 as well as reliable distance estimates for almost all stars. Our analysis is based on the statistical analysis of a rigorously selected sample of 90 high-velocity halo stars from RAVE and a previously published data set. We calibrate and extensively test our method using a suite of cosmological simulations of the formation of Milky Way-sized galaxies. Our best estimate of the local Galactic escape speed, which we define as the minimum speed required to reach three virial radii $R_{340}$, is $533^{+54}_{-41}$ km/s (90% confidence) with an additional 5% systematic uncertainty, where $R_{340}$ is the Galactocentric radius encompassing a mean over-density of 340 times the critical density for closure in the Universe. From the escape speed we further derive estimates of the mass of the Galaxy using a simple mass model with two options for the mass profile of the dark matter halo: an unaltered and an adiabatically contracted Navarro, Frenk & White (NFW) sphere. If we fix the local circular velocity the latter profile yields a significantly higher mass than the un-contracted halo, but if we instead use the statistics on halo concentration parameters in large cosmological simulations as a constraint we find very similar masses for both models. Our best estimate for $M_{340}$, the mass interior to $R_{340}$ (dark matter and baryons), is $1.3^{+0.4}_{-0.3} \times 10^{12}$ M$_\odot$ (corresponding to $M_{200} = 1.6^{+0.5}_{-0.4} \times 10^{12}$ M$_\odot$). This estimate is in good agreement with recently published independent mass estimates based on the kinematics of more distant halo stars and the satellite galaxy Leo I.

We use multi-epoch near-infrared observations from the VISTA survey of the Magellanic Cloud system (VMC) to measure the proper motion of stars of the LMC, in one tile of 1.5 deg^2 centred at (alpha, delta) = (05:59:23.136, -66:20:28.68) and including the South Ecliptic Pole, with respect to their 2MASS position over a time baseline of ~10 years and from VMC observations only, spanning a time range of ~1 year. Stars of different ages are selected from the (J-Ks) vs. Ks diagram and their average coordinate displacement is computed from the difference between Ks band observations for stars as faint as Ks=19 mag. Proper motions are derived by averaging up to seven 2MASS-VMC combinations in the first case and from the slope of the best fit line among the seven VMC epochs in the second case. Separate proper motion values are obtained for variable stars in the field. The proper motion of ~40,000 LMC stars in the tile, with respect to ~8000 background galaxies, obtained from VMC data alone, is mu_alpha cos(delta) = +2.20+/-0.06 (stat) +/-0.29 (sys) and mu_delta = +1.70+/-0.06 (stat) +/-0.30 (sys) mas/yr. This value agrees with recent ground-based determinations in a nearby field but is larger than studies with the HST; the cause of this discrepancy may be due to additional systematic errors in the data. The LMC proper motion is also clearly distinct from the proper motion derived for stars of the Milky Way. The relative proper motion between the foreground stars and the LMC stars is ~5 mas/yr. Furthermore, we measure a decrease of the proper motion with increasing logarithm of stellar age for LMC stars. This study, based on just one VMC tile, shows the potential of the 2MASS-VMC and VMC-VMC combinations for a comprehensive investigation of the proper motion across the Magellanic system. (Abridged)

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database.

Aim: We aim at measuring the chemical gradients of the elements Mg, Al, Si, and Fe along the Galactic radius to provide new constraints on the chemical evolution models of the Galaxy and Galaxy models such as the Besancon model. Methods: We analysed three different samples selected from three independent datasets: a sample of 19,962 dwarf stars selected from the RAVE database, a sample of 10,616 dwarf stars selected from the Geneva-Copenhagen Survey (GCS) dataset, and a mock sample (equivalent to the RAVE sample) created by using the GALAXIA code, which is based on the Besancon model. We measured the chemical gradients as functions of the guiding radius (Rg) at different distances from the Galactic plane reached by the stars along their orbit (Zmax). Results: The chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at Zmax<0.4 kpc and 4.5<Rg(kpc)<9.5, the iron gradient for the RAVE sample is d[Fe/H]/dRg=-0.065 dex kpc^{-1}, whereas for the GCS sample it is d[Fe/H]/dRg=-0.043 dex kpc^{-1} with internal errors +-0.002 and +-0.004 dex kpc^{-1}, respectively. The gradients of the RAVE and GCS samples become flatter at larger Zmax. Conversely, the mock sample has a positive iron gradient of d[Fe/H]/dRg=+0.053+-0.003 dex kpc^{-1} at Zmax<0.4 kpc and remains positive at any Zmax. These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besancon model. The discrepancies between the observational samples and the mock sample can be reduced by i) decreasing the density, ii) decreasing the vertical velocity, and iii) increasing the metallicity of the thick disc in the Besancon model.

In preparation for future, large-scale, multi-object, high-resolution spectroscopic surveys of the Galaxy, we present a series of tests of the precision in radial velocity and chemical abundances that any such project can achieve at a 4m class telescope. We briefly discuss a number of science cases that aim at studying the chemo-dynamical history of the major Galactic components (bulge, thin and thick disks, and halo) - either as a follow-up to the Gaia mission or on their own merits. Based on a large grid of synthetic spectra that cover the full range in stellar parameters of typical survey targets, we devise an optimal wavelength range and argue for a moderately high-resolution spectrograph. As a result, the kinematic precision is not limited by any of these factors, but will practically only suffer from systematic effects, easily reaching uncertainties <1 km/s. Under realistic survey conditions (namely, considering stars brighter than r=16 mag with reasonable exposure times) we prefer an ideal resolving power of R~20000 on average, for an overall wavelength range (with a common two-arm spectrograph design) of [395;456.5] nm and [587;673] nm. We show for the first time on a general basis that it is possible to measure chemical abundance ratios to better than 0.1 dex for many species (Fe, Mg, Si, Ca, Ti, Na, Al, V, Cr, Mn, Co, Ni, Y, Ba, Nd, Eu) and to an accuracy of about 0.2 dex for other species such as Zr, La, and Sr. While our feasibility study was explicitly carried out for the 4MOST facility, the results can be readily applied to and used for any other conceptual design study for high-resolution spectrographs.

The 4MOST consortium is currently halfway through a Conceptual Design study for ESO with the aim to develop a wide-field (>3 square degree, goal >5 square degree), high-multiplex (>1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope (VISTA). 4MOST will run permanently on the telescope to perform a 5 year public survey yielding more than 20 million spectra at resolution R~5000 ({\lambda}=390-1000 nm) and more than 2 million spectra at R~20,000 (395-456.5 nm & 587-673 nm). The 4MOST design is especially intended to complement three key all-sky, space-based observatories of prime European interest: Gaia, eROSITA and Euclid. Initial design and performance estimates for the wide-field corrector concepts are presented. We consider two fibre positioner concepts, a well-known Phi-Theta system and a new R-Theta concept with a large patrol area. The spectrographs are fixed configuration two-arm spectrographs, with dedicated spectrographs for the high- and low-resolution. A full facility simulator is being developed to guide trade-off decisions regarding the optimal field-of-view, number of fibres needed, and the relative fraction of high-to-low resolution fibres. Mock catalogues with template spectra from seven Design Reference Surveys are simulated to verify the science requirements of 4MOST. The 4MOST consortium aims to deliver the full 4MOST facility by the end of 2018 and start delivering high-level data products for both consortium and ESO community targets a year later with yearly increments.

[abridged] In this work we study in detail the kinematics of tidal debris stars to investigate the implications of the new scenario that the observed sample of Hypervelocity stars could partly originate from a dwarf-host galaxy collision. We use a suite of N-body simulations following the encounter of a satellite galaxy with its Milky Way-type host galaxy to gather statistics on the properties of stripped-off stars. We study especially the orbital energy distribution of this population. We quantify the typical pattern in angular and phase space formed by the debris stars. We further develop a simple stripping model predicting the kinematics of stripped-off stars. We show that the distribution of orbital energies in the tidal debris has a typical form which can be described quite accurately by a simple function. Based on this we develop a method to predict the energy distribution which allows us to evaluate the significance and the implications of high velocity stars in satellite tidal debris. Generally tidal collisions of satellite galaxy produce stars which escape into the intragalactic space even if the satellite itself is on a bound orbit. The main parameters determining the maximum energy kick a tidal debris star can get is the initial mass of the satellite and only to a lower extent its orbit. Main contributors to an unbound stellar population created in this way are massive satellites (M_sat > 10^9 M_sun). We thus expect intragalactic stars to have a metallicity higher than the surviving satellite population of the Milky Way. However, the probability that the observed HVS population is significantly contaminated by tidal debris stars appears small in the light of our results.