Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically chosen $\mathrm{[Eu/Fe]} = +1.0$ level previously in use, based on statistical tests of the complete set of RPA data released to date. We also statistically justify the minimum level of [Eu/Fe] for definition of the $r$-I stars, [Eu/Fe] $> +0.3$. Redefining the separation between $r$-I and $r$-II stars will aid in analysis of the possible progenitors of these two classes of stars and whether these signatures arise from separate astrophysical sources at all. Applying this redefinition to previous RPA data, the number of identified $r$-II and $r$-I stars changes to 51 and 121, respectively, from the initial set of data releases published thus far. In this data release, we identify 21 new $r$-II, 111 new $r$-I (plus three re-identified), and 7 new (plus one re-identified) limited-$r$ stars out of a total of 232 target stars, resulting in a total sample of 72 new $r$-II stars, 232 new $r$-I stars, and 42 new limited-$r$ stars identified by the RPA to date.

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.

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.

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.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

The Pristine survey uses narrow-band photometry to derive precise metallicities down to the extremely metal-poor regime ([Fe/H] < -3), and currently consists of over 4 million FGK-type stars over a sky area of $\sim 2~500\, \mathrm{deg}^2$. We focus our analysis on a subsample of $\sim 80~000$ main sequence turnoff stars with heliocentric distances between 6 and 20 kpc, which we take to be a representative sample of the inner halo. The resulting metallicity distribution function (MDF) has a peak at [Fe/H] = -1.6, and a slope of $\Delta$(LogN)/$\Delta[Fe/H] = 1.0 \pm 0.1$ in the metallicity range of -3.4 < [Fe/H] < -2.5. This agrees well with a simple closed-box chemical enrichment model in this range, but is shallower than previous spectroscopic MDFs presented in the literature, suggesting that there may be a larger proportion of metal-poor stars in the inner halo than previously reported. We identify the Monoceros/TriAnd/ACS/EBS/A13 structure in metallicity space in a low latitude field in the anticenter direction, and also discuss the possibility that the inner halo is dominated by a single, large merger event, but cannot strongly support or refute this idea with the current data. Finally, based on the MDF of field stars, we estimate the number of expected metal-poor globular clusters in the Milky Way halo to be 5.4 for [Fe/H] < -2.5 and 1.5 for [Fe/H] < -3, suggesting that the lack of low metallicity globular clusters in the Milky Way is not due simply to statistical undersampling.

This paper focuses on carbon that is one of the most abundant elements in the Universe and is of high importance in the field of nucleosynthesis and galactic and stellar evolution. Even nowadays, the origin of carbon and the relative importance of massive and low- to intermediate-mass stars in producing it is still a matter of debate. In this paper we aim at better understanding the origin of carbon by studying the trends of [C/H], [C/Fe],and [C/Mg] versus [Fe/H], and [Mg/H] for 2133 FGK dwarf stars from the fifth Gaia-ESO Survey internal data release (GES iDR5). The availability of accurate parallaxes and proper motions from Gaia DR2 and radial velocities from GES iDR5 allows us to compute Galactic velocities, orbits and absolute magnitudes and, for 1751 stars, ages via a Bayesian approach. Three different selection methodologies have been adopted to discriminate between thin and thick disk stars. In all the cases, the two stellar groups show different abundance ratios, [C/H], [C/Fe], and [C/Mg], and span different age intervals, with the thick disk stars being, on average, older than those in the thin disk. The behaviours of [C/H], [C/Fe], and [C/Mg] versus [Fe/H], [Mg/H], and age all suggest that C is primarily produced in massive stars like Mg. The increase of [C/Mg] for young thin disk stars indicates a contribution from low-mass stars or the increased C production from massive stars at high metallicities due to the enhanced mass loss. The analysis of the orbital parameters Rmed and |Zmax| support an "inside-out" and "upside-down" formation scenario for the disks of Milky Way.

We find that the combined LF of N- and SC-type stars are consistent with a Gaussian distribution peaking at M_bol~ -5.2 mag. The resulting LF however shows two tails at lower and higher luminosities more extended than those previously found, indicating that AGB carbon stars with Solar metallicity may reach M_bol~-6.0 mag. We find that J-type stars are about half a magnitude fainter on average than N- and SC-type stars, while R-hot stars are half a magnitude brighter than previously found. The Galactic spatial distribution and velocity components of the N-, SC- and J-type stars are very similar, while about 30 % of the R-hot stars in the sample are located at distances larger than ~ 500 pc from the Galactic Plane, and show a significant drift with respect to the local standard of rest. The LF derived for N- and SC-type in the Solar neighbourhood fully agrees with the expected luminosity of stars of 1.5-3 M_o on the AGB. On a theoretical basis, the existence of an extended low luminosity tail would require a contribution of extrinsic low mass carbon stars, while the high luminosity one would imply that stars with mass up to ~5 Mo may become carbon star on the AGB. J-type stars not only differ significantly in their chemical composition with respect to the N- and SC-types but also in their LF, which reinforces the idea that these carbon stars belong to a dvifferent type whose origin is still unknown. The derived luminosities of R-hot stars make these stars unlikely to be in the red-clump as previously claimed. On the other hand, the derived spatial distribution and kinematic properties, together with their metallicity, indicate that most of the N-, SC- and J-type stars belong to the thin disc population, while a significant fraction of R-hot stars show characteristics compatible with the thick disc.

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0\sigma$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

The combination of the final version of the RAVE spectroscopic survey data release 6 with radial velocities and astrometry from Gaia DR2 allows us to identify and create a catalog of single lined binary star candidates (SB1), their inferred orbital parameters, and to inspect possible double lined binary stars (SB2). A probability function for the detection of radial velocity (RV) variations is used for identifying SB1 candidates. The estimation of orbital parameters for main sequence dwarfs is performed by matching the measured RVs with theoretical velocity curves sampling the orbital parameter space. The method is verified by studying a mock sample from the SB9 catalogue. Studying the boxiness and asymmetry of the spectral lines allows us to identify possible SB2 candidates, while matching their spectra to a synthetic library indicates probable properties of their components. From the RAVE catalog we select 37,664 stars with multiple RV measurements and identify 3838 stars as SB1 candidates. Joining RAVE and Gaia DR2 yields 450,646 stars with RVs measured by both surveys and 27,716 of them turn out to be SB1 candidates, which is an increase by an order of magnitude over previous studies. For main sequence dwarf candidates we calculate their most probable orbital parameters: orbital periods are not longer than a few years and primary components have masses similar to the Solar mass. All our results are available via Vizier/CDS.

We present the results of a 3-year long, medium-resolution spectroscopic campaign aimed at identifying very metal-poor stars from candidates selected with the CaHK, metallicity-sensitive Pristine survey. The catalogue consists of a total of 1007 stars, and includes 146 rediscoveries of metal-poor stars already presented in previous surveys, 707 new very metal-poor stars with [Fe/H] < -2.0, and 95 new extremely metal-poor stars with [Fe/H] < -3.0. We provide a spectroscopic [Fe/H] for every star in the catalogue, and [C/Fe] measurements for a subset of the stars (10% with [Fe/H] < -3 and 24% with -3 < [Fe/H] < -2) for which a carbon determination is possible, contingent mainly on the carbon abundance, effective temperature and S/N of the stellar spectra. We find an average carbon enhancement fraction ([C/Fe] >= +0.7) of 41 +- 4% for stars with -3 < [Fe/H] < -2 and 58 +- 14% for stars with [Fe/H] < -3, and report updated success rates for the Pristine survey of 56 % and 23 % to recover stars with [Fe/H] < -2.5 and [Fe/H] < -3, respectively. Finally, we discuss the current status of the survey and its preparation for providing targets to upcoming multi-object spectroscopic surveys such as WEAVE.

In the era of large high-resolution spectroscopic surveys, high-quality spectra can contribute to our understanding of the Galactic chemical evolution, providing chemical abundances belonging to the different nucleosynthesis channels, and also providing constraints to stellar age. Some abundance ratios have been proven to be excellent indicators of stellar ages. We aim at providing an empirical relationship between stellar ages and [C/N] using, as calibrators, open star clusters observed by both the Gaia-ESO and APOGEE surveys. We use stellar parameters and abundances from the Gaia-ESO and APOGEE of the Galactic field and open cluster stars. Ages of star clusters are retrieved from the literature sources and validated using a common set of isochrones. We use the same isochrones to determine, for each age and metallicity, the surface gravity at which the first dredge-up and red giant branch bump occur. We study the effect of extra-mixing processes in our sample of giant stars, and we derive the mean [C/N] in evolved stars, including only stars without evidence of extra-mixing. Combining the Gaia-ESO and APOGEE samples of open clusters, we derive a linear relationship between [C/N] and logarithmic cluster ages. We apply our relationship to selected giant field stars in both Gaia-ESO and APOGEE. We find an age separation between thin and thick disc stars and age trends within their populations, with an increasing age towards lower metallicity populations. With such empirical relationship, we are able to provide an age estimate for giant stars in which C and N abundances are measured. Isochrone fitting is less sensitive for giant than dwarf stars at the turn off. The present method can be thus considered as an additional tool to give an independent estimate of the age of giant stars, with uncertainties in their ages comparable to those obtained using isochrone fitting for dwarf stars.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way.

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 will study the formation history of the Milky Way, and the earliest phases of its chemical enrichment, with a sample of more than 1.5 million stars at high galactic latitude. Elemental abundances of up to 20 elements with a precision of better than 0.2 dex will be derived for these stars. The sample will include members of kinematically coherent substructures, which we will associate with their possible birthplaces by means of their abundance signatures and kinematics, allowing us to test models of galaxy formation. Our target catalogue is also expected to contain 30,000 stars at a metallicity of less than one hundredth that of the Sun. This sample will therefore be almost a factor of 100 larger than currently existing samples of metal-poor stars for which precise elemental abundances are available (determined from high-resolution spectroscopy), enabling us to study the early chemical evolution of the Milky Way in unprecedented detail.

The goal of this survey is to study the formation and evolution of the Milky Way halo to deduce its assembly history and the 3D distribution of mass in the Milky Way. The combination of multi-band photometry, Gaia proper motion and parallax data, and radial velocities and the metallicity and elemental abundances obtained from low-resolution spectra of halo giants with 4MOST, will yield an unprecedented characterisation of the Milky Way halo and its interface with the thick disc. The survey will produce a volume- and magnitude-limited complete sample of giant stars in the halo. It will cover at least 10,000 square degrees of high Galactic latitude, and measure line-of-sight velocities with a precision of 1-2 km/s as well as metallicities to within 0.2 dex.

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.

A new moderately r-process-enhanced metal-poor star, RAVE J093730.5-062655, has been identified in the Milky Way halo as part of an ongoing survey by the R-Process Alliance. The temperature and surface gravity indicate that J0937-0626 is likely a horizontal branch star. At [Fe/H] = -1.86, J0937-0626 is found to have subsolar [X/Fe] ratios for nearly every light, alpha, and Fe-peak element. The low [alpha/Fe] ratios can be explained by an ~0.6 dex excess of Fe; J0937-0626 is therefore similar to the subclass of "iron-enhanced" metal-poor stars. A comparison with Milky Way field stars at [Fe/H] = -2.5 suggests that J0937-0626 was enriched in material from an event, possibly a Type Ia supernova, that created a significant amount of Cr, Mn, Fe, and Ni and smaller amounts of Ca, Sc, Ti, and Zn. The r-process enhancement of J0937-0626 is likely due to a separate event, which suggests that its birth environment was highly enriched in r-process elements. The kinematics of J0937-0626, based on Gaia DR2 data, indicate a retrograde orbit in the Milky Way halo; J0937-0626 was therefore likely accreted from a dwarf galaxy that had significant r-process enrichment.

This paper presents the detailed abundances and r-process classifications of 126 newly identified metal-poor stars as part of an ongoing collaboration, the R-Process Alliance. The stars were identified as metal-poor candidates from the RAdial Velocity Experiment (RAVE) and were followed-up at high spectral resolution (R~31,500) with the 3.5~m telescope at Apache Point Observatory. The atmospheric parameters were determined spectroscopically from Fe I lines, taking into account <3D> non-LTE corrections and using differential abundances with respect to a set of standards. Of the 126 new stars, 124 have [Fe/H]<-1.5, 105 have [Fe/H]<-2.0, and 4 have [Fe/H]<-3.0. Nine new carbon-enhanced metal-poor stars have been discovered, 3 of which are enhanced in r-process elements. Abundances of neutron-capture elements reveal 60 new r-I stars (with +0.3<=[Eu/Fe]<=+1.0 and [Ba/Eu]<0) and 4 new r-II stars (with [Eu/Fe]>+1.0). Nineteen stars are found to exhibit a `limited-r' signature ([Sr/Ba]>+0.5, [Ba/Eu]<0). For the r-II stars, the second- and third-peak main r-process patterns are consistent with the r-process signature in other metal-poor stars and the Sun. The abundances of the light, alpha, and Fe-peak elements match those of typical Milky Way halo stars, except for one r-I star which has high Na and low Mg, characteristic of globular cluster stars. Parallaxes and proper motions from the second Gaia data release yield UVW space velocities for these stars which are consistent with membership in the Milky Way halo. Intriguingly, all r-II and the majority of r-I stars have retrograde orbits, which may indicate an accretion origin.

Finding solar siblings, that is, stars that formed in the same cluster as the Sun, will yield information about the conditions at the Sun's birthplace. We search for solar sibling candidates in AMBRE, the very large spectra database of solar vicinity stars. Since the ages and chemical abundances of solar siblings are very similar to those of the Sun, we carried out a chemistry- and age-based search for solar sibling candidates. We used high-resolution spectra to derive precise stellar parameters and chemical abundances of the stars. We used these spectroscopic parameters together with Gaia DR2 astrometric data to derive stellar isochronal ages. Gaia data were also used to study the kinematics of the sibling candidates. From the about 17000 stars that are characterized within the AMBRE project, we first selected 55 stars whose metallicities are closest to the solar value (-0.1 < [Fe/H] < 0.1 dex). For these stars we derived precise chemical abundances of several iron-peak, alpha- and neutron-capture elements, based on which we selected 12 solar sibling candidates with average abundances and metallicities between -0.03 to 0.03 dex. Our further selection left us with 4 candidates with stellar ages that are compatible with the solar age within observational uncertainties. For the 2 of the hottest candidates, we derived the carbon isotopic ratios, which are compatible with the solar value. HD186302 is the most precisely characterized and probably the most probable candidate of our 4 best candidates. Very precise chemical characterization and age estimation is necessary to identify solar siblings. We propose that in addition to typical chemical tagging, the study of isotopic ratios can give further important information about the relation of sibling candidates with the Sun. Ideally, asteroseismic age determinations of the candidates could solve the problem of imprecise isochronal ages.

We test the performance of the semi-analytic self-consistent Just-Jahrei{\ss} disc model (JJ model) with the astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) sub-catalogue of the first Gaia data release (Gaia DR1), as well as the radial velocities from the fifth data release of the Radial Velocity Experiment survey (RAVE DR5). We use a sample of 19,746 thin disc stars from the TGAS$\times$RAVE cross-match selected in the local solar cylinder of 300 pc radius and 1 kpc height below the Galactic plane and simulate this sample via the forward modelling technique. First, we convert the predicted vertical density laws of the thin disc populations into a mock sample. Then the obtained mock populations are reddened with a 3D dust map and are subjected to the selection criteria corresponding to the RAVE and TGAS observational limitations as well as to additional cuts applied to the data sample. We calculate the quantities of interest separately at different heights above the Galactic plane taking into account the distance error effects separately in horizontal and vertical directions. We investigate the simulated sample in terms of the vertical number density profiles, Hess diagrams and velocity distribution functions. Basing on a good agreement of our simulations with the data, we conclude that our fiducial disc model confidently reproduces the vertical trends in the thin disc stellar population properties. Thus, it can serve as a starting point for the future extension of the JJ model to other Galactocentric distances.

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.

We present a semi-empirical, largely model-independent approach for estimating Galactic birth radii, r_birth, for Milky Way disk stars. The technique relies on the justifiable assumption that a negative radial metallicity gradient in the interstellar medium (ISM) existed for most of the disk lifetime. Stars are projected back to their birth positions according to the observationally derived age and [Fe/H] with no kinematical information required. Applying our approach to the AMBRE:HARPS and HARPS-GTO local samples, we show that we can constrain the ISM metallicity evolution with Galactic radius and cosmic time, [Fe/H]_ISM(r, t), by requiring a physically meaningful r_birth distribution. We find that the data are consistent with an ISM radial metallicity gradient that flattens with time from ~-0.15 dex/kpc at the beginning of disk formation, to its measured present-day value (-0.07 dex/kpc). We present several chemo-kinematical relations in terms of mono-r_birth populations. One remarkable result is that the kinematically hottest stars would have been born locally or in the outer disk, consistent with thick disk formation from the nested flares of mono-age populations and predictions from cosmological simulations. This phenomenon can be also seen in the observed age-velocity dispersion relation, in that its upper boundary is dominated by stars born at larger radii. We also find that the flatness of the local age-metallicity relation (AMR) is the result of the superposition of the AMRs of mono-r_birth populations, each with a well-defined negative slope. The solar birth radius is estimated to be 7.3+-0.6 kpc, for a current Galactocentric radius of 8 kpc.

The abundance ratio N/O is a useful tool to study the interplay of galactic processes, e.g. star formation efficiency, time-scale of infall and outflow loading factor We aim to trace log(N/O) versus [Fe/H] in the Milky Way and to compare it with a set of chemical evolution models to understand the role of infall, outflow and star formation efficiency in the building-up of the Galactic disc. We use the abundances from idr2-3, idr4, idr5 data releases of the Gaia-ESO Survey both for Galactic field and open cluster stars.We determine membership and average composition of open clusters and we separate thin and thick disc field stars.We consider the effect of mixing in the abundance of N in giant stars. We compute a grid of chemical evolution models, suited to reproduce the main features of our Galaxy, exploring the effects of the star formation efficiency, the infall time-scale and the differential outflow. With our samples, we map the metallicity range -0.6<[Fe/H]<0.3 with a corresponding -1.2<log(N/O)<-0.2, where the secondary production of N dominates. Thanks to the wide range of Galactocentric distances covered by our samples, we can distinguish the behaviour of log(N/O) in different parts of the Galaxy. Our spatially resolved results allow us to distinguish differences in the evolution of N/O with Galactocentric radius. Comparing the data with our models, we can characterise the radial regions of our Galaxy. A shorter infall time-scale is needed in the inner regions, while the outer regions need a longer infall time-scale, coupled with a higher star formation efficiency. We compare our results with nebular abundances obtained in MaNGA galaxies, finding in our Galaxy a much wider range of log(N/O) than in integrated observations of external galaxies of similar stellar mass, but similar to the ranges found in studies of individual H ii regions.

Several works have found an increase of the abundances of the s-process neutron-capture elements in the youngest Galactic stellar populations, giving important constraints to stellar and Galactic evolution. We aim to trace the abundance patterns and the time-evolution of five s-process elements in the first peak, Y and Zr, and in the second peak, Ba, La and Ce using the Gaia-ESO idr5 results. From the UVES spectra of cluster member stars, we determined the average composition of clusters with ages >0.1 Gyr. We derived statistical ages and distances of field stars, and we separated them in thin and thick disc populations. We studied the time evolution and dependence on metallicity of abundance ratios using open clusters and field stars. Using our large and homogeneous sample of open clusters, thin and thick disc stars, spanning an age range larger than 10 Gyr, we confirm an increase towards young ages of s-process abundances in the Solar neighbourhood. These trends are well defined for open clusters and stars located nearby the solar position and they may be explained by a late enrichment due to significant contribution to the production of these elements from long-living low-mass stars. At the same time, we found a strong dependence of the s-process abundance ratios with the Galactocentric distance and with the metallicity of the clusters and field stars. Our results, derived from the largest and homogeneous sample of s-process abundances in the literature, confirm the growth with decreasing stellar ages of the s-process abundances in both field and open cluster stars. At the same time, taking advantage of the abundances of open clusters located in a wide Galactocentric range, they open a new view on the dependence of the s-process evolution on the metallicity and star formation history, pointing to different behaviours at various Galactocentric distances.

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 present results from a medium-resolution (R ~ 2, 000) spectroscopic follow-up campaign of 1,694 bright (V < 13.5), very metal-poor star candidates from the RAdial Velocity Experiment (RAVE). Initial selection of the low-metallicity targets was based on the stellar parameters published in RAVE Data Releases 4 and 5. Follow-up was accomplished with the Gemini-N and Gemini-S, the ESO/NTT, the KPNO/Mayall, and the SOAR telescopes. The wavelength coverage for most of the observed spectra allows for the determination of carbon and {\alpha}-element abundances, which are crucial for con- sidering the nature and frequency of the carbon-enhanced metal-poor (CEMP) stars in this sample. We find that 88% of the observed stars have [Fe/H] <= -1.0, 61% have [Fe/H] <= -2.0, and 3% have [Fe/H] <= -3.0 (with four stars at [Fe/H] <= -3.5). There are 306 CEMP star candidates in this sample, and we identify 169 CEMP Group I, 131 CEMP Group II, and 6 CEMP Group III stars from the A(C) vs. [Fe/H] diagram. Inspection of the [alpha/C] abundance ratios reveals that five of the CEMP Group II stars can be classified as "mono-enriched second-generation" stars. Gaia DR1 matches were found for 734 stars, and we show that transverse velocities can be used as a confirmatory selection criteria for low-metallicity candidates. Selected stars from our validated list are being followed-up with high-resolution spectroscopy, to reveal their full chemical abundance patterns for further studies.

The second Gaia data release (Gaia-DR2) contains, beyond the astrometry, three-band photometry for 1.38 billion sources. We have used these three broad bands to infer stellar effective temperatures, Teff, for all sources brighter than G=17 mag with Teff in the range 3000-10 000 K (161 million sources). Using in addition the parallaxes, we infer the line-of-sight extinction, A_G, and the reddening, E[BP-RP], for 88 million sources. Together with a bolometric correction we derive luminosity and radius for 77 million sources. These quantities as well as their estimated uncertainties are part of Gaia-DR2. Here we describe the procedures by which these quantities were obtained, including the underlying assumptions, comparison with literature estimates, and the limitations of our results. Typical accuracies are of order 324 K (Teff), 0.46 mag (A_G), 0.23 mag (E[BP-RP]), 15% (luminosity), and 10% (radius). Being based on only a small number of observable quantities and limited training data, our results are necessarily subject to some extreme assumptions that can lead to strong systematics in some cases (not included in the aforementioned accuracy estimates). One aspect is the non-negativity contraint of our estimates, in particular extinction. Yet in several regions of parameter space our results show very good performance, for example for red clump stars and solar analogues. Large uncertainties render the extinctions less useful at the individual star level, but they show good performance for ensemble estimates. We identify regimes in which our parameters should and should not be used and we define a "clean" sample. Despite the limitations, this is the largest catalogue of uniformly-inferred stellar parameters to date. More precise and detailed astrophysical parameters based on the full BP/RP spectrophotometry are planned as part of the third Gaia data release.

We combine parallaxes from the first Gaia data release with the spectrophotometric distance estimation framework for stars in the fifth RAVE survey data release. The combined distance estimates are more accurate than either determination in isolation - uncertainties are on average two times smaller than for RAVE-only distances (three times smaller for dwarfs), and 1.4 times smaller than TGAS parallax uncertainties (two times smaller for giants). We are also able to compare the estimates from spectrophotometry to those from Gaia, and use this to assess the reliability of both catalogues and improve our distance estimates. We find that the distances to the lowest log g stars are, on average, overestimated and caution that they may not be reliable. We also find that it is likely that the Gaia random uncertainties are smaller than the reported values. As a byproduct we derive ages for the RAVE stars, many with relative uncertainties less than 20 percent. These results for 219566 RAVE sources have been made publicly available, and we encourage their use for studies that combine the radial velocities provided by RAVE with the proper motions provided by Gaia. A sample that we believe to be reliable can be found by taking only the stars with the flag notification flag_any=0

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.

Lithium abundance in most of the warm metal-poor main sequence stars shows a constant plateau (A(Li)~2.2 dex) and then the upper envelope of the lithium vs. metallicity distribution increases as we approach solar metallicity. Meteorites, which carry information about the chemical composition of the interstellar medium at the solar system formation time, show a lithium abundance A(Li)~3.26 dex. This pattern reflects the Li enrichment history of the interstellar medium during the Galaxy lifetime. After the initial Li production in Big Bang Nucleosynthesis, the sources of the enrichment include AGB stars, low-mass red giants, novae, type II supernovae, and Galactic cosmic rays. The total amount of enriched Li is sensitive to the relative contribution of these sources. Thus different Li enrichment histories are expected in the Galactic thick and thin disc. We investigate the main sequence stars observed with UVES in Gaia-ESO Survey iDR4 catalog and find a Li-[alpha/Fe] anticorrelation independent of [Fe/H], Teff, and log(g). Since in stellar evolution different {\alpha} enhancements at the same metallicity do not lead to a measurable Li abundance change, the anticorrelation indicates that more Li is produced during the Galactic thin disc phase than during the Galactic thick disc phase. We also find a correlation between the abundance of Li and s-process elements Ba and Y, and they both decrease above the solar metallicity, which can be explained in the framework of the adopted Galactic chemical evolution models.

We construct the rotation curve of the Milky Way in the extended solar neighbourhood using a sample of SEGUE (Sloan Extension for Galactic Understanding and Exploration) G-dwarfs. We investigate the rotation curve shape for the presence of any peculiarities just outside the solar radius as has been reported by some authors. We approach the problem in a framework of classical Jeans analysis. Using the most recent data from RAVE (RAdial Velocity Experiment), we determine the solar peculiar velocity and the radial scalelengths for the three populations of different metallicities representing the Galactic thin disc. Then with the same binning in metallicity for the SEGUE G-dwarfs, we construct the rotation curve in the range of Galactocentric distances 7-10 kpc. We derive the circular velocity by correcting the mean tangential velocity for the asymmetric drift in each distance bin. With SEGUE data we also calculate the radial scalelength of the thick disc taking as known the derived peculiar motion of the Sun and the slope of the rotation curve. The rotation curve constructed through SEGUE G-dwarfs appears to be smooth in the selected radial range. The local kinematics of the thin disc rotation as determined in the framework of our new careful analysis does not favour the presence of a massive overdensity ring just outside the solar radius.

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.

We attempt to determine the relative fraction of stars that have undergone significant radial migration by studying the orbital properties of metal-rich ([Fe/H]$>0.1$) stars within 2 kpc of the Sun using a sample of more than 3,000 stars selected from iDR4 of the Gaia-ESO Survey. We investigate the kinematic properties, such as velocity dispersion and orbital parameters, of stellar populations near the sun as a function of [Mg/Fe] and [Fe/H], which could show evidence of a major merger in the past history of the Milky Way. This was done using the stellar parameters from the Gaia-ESO Survey along with proper motions from PPMXL to determine distances, kinematics, and orbital properties for these stars to analyze the chemodynamic properties of stellar populations near the Sun. Analyzing the kinematics of the most metal-rich stars ([Fe/H]$>0.1$), we find that more than half have small eccentricities ($e<0.2$) or are on nearly circular orbits. Slightly more than 20\% of the metal-rich stars have perigalacticons $R_p>7$ kpc. We find that the highest [Mg/Fe], metal-poor populations have lower vertical and radial velocity dispersions compared to lower [Mg/Fe] populations of similar metallicity by $\sim10$ km s$^{-1}$. The median eccentricity increases linearly with [Mg/Fe] across all metallicities, while the perigalacticon decreases with increasing [Mg/Fe] for all metallicities. Finally, the most [Mg/Fe]-rich stars are found to have significant asymmetric drift and rotate more than 40 km s$^{-1}$ slower than stars with lower [Mg/Fe] ratios. While our results cannot constrain how far stars have migrated, we propose that migration processes are likely to have played an important role in the evolution of the Milky Way, with metal-rich stars migrating from the inner disk toward to solar neighborhood and past mergers potentially driving enhanced migration of older stellar populations in the disk.

We investigate here the effect of the selection function on the metallicity distribution function (MDF) and on the vertical metallicity gradient by studying similar lines of sight using four different spectroscopic surveys (APOGEE, LAMOST, RAVE, and Gaia-ESO), which have different targeting strategies and therefore different selection functions. We use common fields between the spectroscopic surveys of APOGEE, LAMOST, RAVE (ALR) and APOGEE, RAVE, Gaia-ESO (AGR) and use two stellar population synthesis models, GALAXIA and TRILEGAL, to create mock fields for each survey. We apply the selection function in the form of colour and magnitude cuts of the respective survey to the mock fields to replicate the observed source sample. We make a basic comparison between the models to check which best reproduces the observed sample distribution. We carry out a quantitative comparison between the synthetic MDF from the mock catalogues using both models to understand the effect of the selection function on the MDF and on the vertical metallicity gradient. Using both models, we find a negligible effect of the selection function on the MDF for APOGEE, LAMOST, and RAVE. We find a negligible selection function effect on the vertical metallicity gradients as well, though GALAXIA and TRILEGAL have steeper and shallower slopes, respectively, than the observed gradient. After applying correction terms on the metallicities of RAVE and LAMOST with respect to our reference APOGEE sample, our observed vertical metallicity gradients between the four surveys are consistent within 1-sigma. We also find consistent gradient for the combined sample of all surveys in ALR and AGR. We estimated a mean vertical metallicity gradient of -0.241+/-0.028 dex kpc-1. There is a significant scatter in the estimated gradients in the literature, but our estimates are within their ranges.

We apply the twin method to determine parallaxes to 232,545 stars of the RAVE survey using the parallaxes of Gaia DR1 as a reference. To search for twins in this large dataset, we apply the t-stochastic neighbour embedding t-SNE projection which distributes the data according to their spectral morphology on a two dimensional map. From this map we choose the twin candidates for which we calculate a chi^2 to select the best sets of twins. Our results show a competitive performance when compared to other model-dependent methods relying on stellar parameters and isochrones. The power of the method is shown by finding that the accuracy of our results is not significantly affected if the stars are normal or peculiar since the method is model free. We find twins for 60% of the RAVE sample which is not contained in TGAS or that have TGAS uncertainties which are larger than 20%. We could determine parallaxes with typical errors of 28%. We provide a complementary dataset for the RAVE stars not covered by TGAS, or that have TGAS uncertainties which are larger than 20%, with model-free parallaxes scaled to the Gaia measurements.

The paucity of hypervelocity stars (HVSs) known to date has severely hampered their potential to investigate the stellar population of the Galactic Centre and the Galactic Potential. The first Gaia data release gives an opportunity to increase the current sample. The challenge is the disparity between the expected number of hypervelocity stars and that of bound background stars. We have applied a novel data mining algorithm based on machine learning techniques, an artificial neural network, to the Tycho-Gaia astrometric solution (TGAS) catalogue. With no pre-selection of data, we could exclude immediately $\sim 99 \%$ of the stars in the catalogue and find 80 candidates with more than $90\%$ predicted probability to be HVSs, based only on their position, proper motions, and parallax. We have cross-checked our findings with other spectroscopic surveys, determining radial velocities for 30 and spectroscopic distances for 5 candidates. In addition, follow-up observations have been carried out at the Isaac Newton Telescope for 22 stars, for which we obtained radial velocities and distance estimates. We discover 14 stars with a total velocity in the Galactic rest frame > 400 km/s, and 5 of these have a probability $>50\%$ of being unbound from the Milky Way. Tracing back their orbits in different Galactic potential models we find one possible unbound HVS with velocity $\sim$ 520 km/s, 5 bound HVSs, and, notably, 5 runaway stars with median velocity between 400 and 780 km/s. At the moment, uncertainties in the distance estimates and ages are too large to confirm the nature of our candidates by narrowing down their ejection location, and we wait for future Gaia releases to validate the quality of our sample. This test successfully demonstrates the feasibility of our new data mining routine.

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 present a novel analysis of the metal-poor star sample in the complete Radial Velocity Experiment (RAVE) Data Release 5 catalog with the goal of identifying and characterizing all very metal-poor stars observed by the survey. Using a three-stage method, we first identified the candidate stars using only their spectra as input information. We employed an algorithm called t-SNE to construct a low-dimensional projection of the spectrum space and isolate the region containing metal-poor stars. Following this step, we measured the equivalent widths of the near-infrared CaII triplet lines with a method based on flexible Gaussian processes to model the correlated noise present in the spectra. In the last step, we constructed a calibration relation that converts the measured equivalent widths and the color information coming from the 2MASS and WISE surveys into metallicity and temperature estimates. We identified 877 stars with at least a 50% probability of being very metal-poor $(\rm [Fe/H] < -2\,\rm dex)$, out of which 43 are likely extremely metal-poor $(\rm [Fe/H] < -3\,\rm dex )$. The comparison of the derived values to a small subsample of stars with literature metallicity values shows that our method works reliably and correctly estimates the uncertainties, which typically have values $\sigma_{\rm [Fe/H]} \approx 0.2\,\mathrm{dex}$. In addition, when compared to the metallicity results derived using the RAVE DR5 pipeline, it is evident that we achieve better accuracy than the pipeline and therefore more reliably evaluate the very metal-poor subsample. Based on the repeated observations of the same stars, our method gives very consistent results. The method used in this work can also easily be extended to other large-scale data sets, including to the data from the Gaia mission and the upcoming 4MOST survey.

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

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.

The spatial distribution of elemental abundances in the disc of our Galaxy gives insights both on its assembly process and subsequent evolution, and on the stellar nucleogenesis of the different elements. Gradients can be traced using several types of objects as, for instance, (young and old) stars, open clusters, HII regions, planetary nebulae. We aim at tracing the radial distributions of abundances of elements produced through different nucleosynthetic channels -the alpha-elements O, Mg, Si, Ca and Ti, and the iron-peak elements Fe, Cr, Ni and Sc - by using the Gaia-ESO idr4 results of open clusters and young field stars. From the UVES spectra of member stars, we determine the average composition of clusters with ages >0.1 Gyr. We derive statistical ages and distances of field stars. We trace the abundance gradients using the cluster and field populations and we compare them with a chemo-dynamical Galactic evolutionary model. Results. The adopted chemo-dynamical model, with the new generation of metallicity-dependent stellar yields for massive stars, is able to reproduce the observed spatial distributions of abundance ratios, in particular the abundance ratios of [O/Fe] and [Mg/Fe] in the inner disc (5 kpc<RGC <7 kpc), with their differences, that were usually poorly explained by chemical evolution models. Often, oxygen and magnesium are considered as equivalent in tracing alpha-element abundances and in deducing, e.g., the formation time-scales of different Galactic stellar populations. In addition, often [alpha/Fe] is computed combining several alpha-elements. Our results indicate, as expected, a complex and diverse nucleosynthesis of the various alpha-elements, in particular in the high metallicity regimes, pointing towards a different origin of these elements and highlighting the risk of considering them as a single class with common features.

Correlations between stellar chemistry and kinematics have long been used to gain insight into the evolution of the Milky Way Galaxy. Orbital angular momentum is a key physical parameter and it is often estimated from three-dimensional space motions. We here demonstrate the lower uncertainties that can be achieved in the estimation of one component of velocity through selection of stars in key directions and use of line-of-sight velocity alone (i.e. without incorporation of proper motion data). In this first paper we apply our technique to stars observed in the direction of Galactic rotation in the APOGEE survey. We first derive the distribution of azimuthal velocities, $v_\phi$, then from these and observed radial coordinates, estimate the stellar guiding centre radii, $R_g$, within $6.9\leq R \leq 10$ kpc with uncertainties smaller than (or of the order of) 1kpc. We show that there is no simple way to select a clean stellar sample based on low errors on proper motions and distances to obtain high-quality 3D velocities and hence one should pay particular attention when trying to identify kinematically peculiar stars based on velocities derived using the proper motions. Using our azimuthal velocity estimations, we investigate the joint distribution of elemental abundances and rotational kinematics free from the blurring effects of epicyclic motions, and we derive the $\partial v_\phi / \partial R_g$ trends for the thin and thick discs as a function of radius. Our analysis provides further evidence for radial migration within the thin disc and hints against radial migration playing a significant role in the evolution of the thick disc.

We present a set of 87 RAVE stars with detected solar like oscillations, observed during Campaign 1 of the K2 mission (RAVE K2-C1 sample). This dataset provides a useful benchmark for testing the gravities provided in RAVE Data Release 4 (DR4), and is key for the calibration of the RAVE Data Release 5 (DR5). In the present work, we use two different pipelines, GAUFRE (Valentini et al. 2013) and Sp_Ace (Boeche & Grebel 2016), to determine atmospheric parameters and abundances by fixing log(g) to the seismic one. Our strategy ensures highly consistent values among all stellar parameters, leading to more accurate chemical abundances. A comparison of the chemical abundances obtained here with and without the use of seismic log(g) information has shown that an underestimated (overestimated) gravity leads to an underestimated (overestimated) elemental abundance (e.g. [Mg/H] is underestimated by ~0.25 dex when the gravity is underestimated by 0.5 dex). We then perform a comparison between the seismic gravities and the spectroscopic gravities presented in the RAVE DR4 catalogue, extracting a calibration for log(g) of RAVE giants in the colour interval 0.50<(J - Ks)<0.85. Finally, we show a comparison of the distances, temperatures, extinctions (and ages) derived here for our RAVE K2-C1 sample with those derived in RAVE DR4 and DR5.DR5 performs better than DR4 thanks to the seismic calibration, although discrepancies can still be important for objects for which the difference between DR4/DR5 and seismic gravities differ by more than ~0.5 dex. The method illustrated in this work will be used for analysing RAVE targets present in the other K2 campaigns, in the framework of Galactic Archaeology investigations.

A large sample of over 38,000 chromospherically active candidate solar-like stars and cooler dwarfs from the RAVE survey is addressed in this paper. An improved activity identification with respect to the previous study was introduced to build a catalog of field stars in the Solar neighborhood with an excess emission flux in the calcium infrared triplet wavelength region. The central result of this work is the calibration of the age--activity relation for the main sequence dwarfs in a range from a few $10 \; \mathrm{Myr}$ up to a few Gyr. It enabled an order of magnitude age estimation of the entire active sample. Almost 15,000 stars are shown to be younger than $1\;\mathrm{Gyr}$ and $\sim$2000 younger than $100\;\mathrm{Myr}$. The young age of the most active stars is confirmed by their position off the main sequence in the $J-K$ versus $N_{UV}-V$ diagram showing strong ultraviolet excess, mid-infrared excess in the $J-K$ versus $W_1-W_2$ diagram and very cool temperatures ($J-K>0.7$). They overlap with the reference pre-main sequence RAVE stars often displaying X-ray emission. The activity level increasing with the color reveals their different nature from the solar-like stars and probably represents an underlying dynamo generating magnetic fields in cool stars. 50\% of the RAVE objects from DR5 are found in the TGAS catalog and supplemented with accurate parallaxes and proper motions by Gaia. This makes the database of a large number of young stars in a combination with RAVE's radial velocities directly useful as a tracer of the very recent large-scale star formation history in the Solar neighborhood. The data are available online in the Vizier database.

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.

The orbits, atmospheric parameters, chemical abundances, and ages of individual stars in the Milky Way provide the most comprehensive illustration of galaxy formation available. The Tycho-Gaia Astrometric Solution (TGAS) will deliver astrometric parameters for the largest ever sample of Milky Way stars, though its full potential cannot be realized without the addition of complementary spectroscopy. Among existing spectroscopic surveys, the RAdial Velocity Experiment (RAVE) has the largest overlap with TGAS ($\gtrsim$200,000 stars). We present a data-driven re-analysis of 520,781 RAVE spectra using The Cannon. For red giants, we build our model using high-fidelity APOGEE stellar parameters and abundances for stars that overlap with RAVE. For main-sequence and sub-giant stars, our model uses stellar parameters from the K2/EPIC. We derive and validate effective temperature $T_{\rm eff}$, surface gravity $\log{g}$, and chemical abundances of up to seven elements (O, Mg, Al, Si, Ca, Fe, Ni). We report a total of 1,685,851 elemental abundances with a typical precision of 0.07 dex, a substantial improvement over previous RAVE data releases. The synthesis of RAVE-on and TGAS is the most powerful data set for chemo-dynamic analyses of the Milky Way ever produced.

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $\alpha$-low (`thin disc') and $\alpha$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $\alpha$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrm{\phi}}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrm{\phi}}$/$\partial$[Fe/H] for the $\alpha$-high disc, consistent with results from high-resolution surveys. We also find differences between the $\alpha$-low and $\alpha$-high discs in all three components of velocity dispersion. We discuss the implications of an $\alpha$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$\alpha$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy-peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal rich stars also exhibit peaky features in their line-of-sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy-peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The 'banana' or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data.

The Gaia-ESO Survey was designed to target all major Galactic components (i.e., bulge, thin and thick discs, halo and clusters), with the goal of constraining the chemical and dynamical evolution of the Milky Way. This paper presents the methodology and considerations that drive the selection of the targeted, allocated and successfully observed Milky Way field stars. The detailed understanding of the survey construction, specifically the influence of target selection criteria on observed Milky Way field stars is required in order to analyse and interpret the survey data correctly. We present the target selection process for the Milky Way field stars observed with VLT/FLAMES and provide the weights that characterise the survey target selection. The weights can be used to account for the selection effects in the Gaia-ESO Survey data for scientific studies. We provide a couple of simple examples to highlight the necessity of including such information in studies of the stellar populations in the Milky Way.

Multi-object spectrographs have opened a new window on the analyses of the chemo-dynamical properties of old Milky Way stars. These analyses allow us to trace back the internal mechanisms and the external factors that have influenced the evolution of our Galaxy, and therefore understand fundamental aspects of galaxy evolution in general. Here, we present recent results from the RAdial Velocity Experiment (RAVE) and the Gaia-ESO survey. These surveys explore the Milky Way properties in different ways, in terms of sample size and selection, magnitude range, and spectral resolution. We focus here on (i) the first direct detection of evidence for radial migration within the thin disc, providing insight into the history of spiral structure of the Milky Way, and (ii) the chemo-dynamical characterisation of the metal-weak thick and thin discs, for which chemo-dynamical models still have difficulties in reproducing.

We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Three stars, which are among the 4 most likely candidates for NGC 3201 tidal debris, are separated by 80 degrees on the sky yet are well matched by the 12 Gyr, [Fe/H] = -1.5 isochrone appropriate for the cluster. This is the first time tidal debris around this cluster has been reported over such a large spatial extent, with implications for the cluster$'$s origin and dynamical evolution.

We study the chemodynamical properties of the Carina dwarf spheroidal by combining an intermediate spectroscopic resolution dataset of more than 900 red giant and red clump stars, with high-precision photometry to derive the atmospheric parameters, metallicities and age estimates for our targets. Within the red giant branch population, we find evidence for the presence of three distinct stellar sub-populations with different metallicities, spatial distributions, kinematics and ages. As in the Fornax and Sculptor dwarf spheroidals, the subpopulation with the lowest average metallicity is more extended and kinematically hotter than all other populations. However, we identify an inversion in the parallel ordering of metallicity, kinematics and characteristic length scale in the two most metal rich subpopulations, which therefore do not contribute to a global negative chemical gradient. Contrary to common trends in the chemical properties with radius, the metal richest population is more extended and mildly kinematically hotter than the main component of intermediate metallicity. More investigations are required to ascertain the nature of this inversion, but we comment on the mechanisms that might have caused it.

(Abridged) Recent spectroscopic surveys have begun to explore the Galactic disk system outside the solar neighborhood on the basis of large data samples. In this way, they provide valuable information for testing spatial and temporal variations of disk structure kinematics and chemical evolution. We used a Gaussian mixture model algorithm, as a rigurous mathematical approach, to separate in the [Mg/Fe] vs. [Fe/H] plane a clean disk star subsample from the Gaia-ESO survey internal data release 2. We find that the sample is separated into five groups associated with major Galactic components; the metal-rich end of the halo, the thick disk, and three subgroups for the thin disk sequence. This is confirmed with a sample of red clump stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The two metal-intermediate and metal-rich groups of the thin disk decomposition ([Fe/H]>-0.25 dex) highlight a change in the slope at solar metallicity. This holds true at different radial regions. The distribution of Galactocentric radial distances of the metal-poor part of the thin disk ([Fe/H]<-0.25 dex) is shifted to larger distances than those of the more metal-rich parts. Moreover, the metal-poor part of the thin disk presents indications of a scale height intermediate between those of the thick and the rest of the thin disk, and it displays higher azimuthal velocities than the latter. These stars might have formed and evolved in parallel and/or dissociated from the inside-out formation taking place in the internal thin disk. Their enhancement levels might be due to their origin from gas pre-enriched by outflows from the thick disk or the inner halo. The smooth trends of their properties (their spatial distribution with respect to the plane, in particular) with [Fe/H] and [Mg/Fe] suggested by the data indicates a quiet dynamical evolution, with no relevant merger events.

Understanding the history and the evolution of the Milky Way disc is one of the main goals of modern astrophysics. We study the velocity dispersion behaviour of Galactic disc stars as a function of the [Mg/Fe] ratio, which can be used as a proxy of relative age. This key relation is essential to constrain the formation mechanisms of the disc stellar populations as well as the cooling processes. We used the recommended parameters and chemical abundances of 7800 FGK Milky Way field stars from the second internal data release of the Gaia-ESO Survey. These stars were observed with the GIRAFFE spectrograph, and cover a large spatial volume (6<R<10kpc and |Z|<2kpc). Based on the [Mg/Fe] and [Fe/H] ratios, we separated the thin- from the thick-disc sequence. From analysing the Galactocentric velocity of the stars for the thin disc, we find a weak positive correlation between Vphi and [Fe/H], due to a slowly rotating Fe-poor tail. For the thick disc, a strong correlation with [Fe/H] and [Mg/Fe] is established. We have detected an inversion of the radial velocity dispersion with [Mg/Fe] for thick-disc stars with [Fe/H]<-0.1dex and [Mg/Fe]>+0.2dex. First, the velocity dispersion increases with [Mg/Fe] at all [Fe/H] ratios for the thin-disc stars, and then it decreases for the thick-disc at the highest [Mg/Fe] abundances. Similar trends are observed within the errors for the azimuthal velocity dispersion, while a continuous increase with [Mg/Fe] is observed for the vertical velocity dispersion. The velocity dispersion decrease agrees with previous measurements of the RAVE survey, although it is observed here for a greater metallicity interval and a larger spatial volume. We confirm the existence of [Mg/Fe]-rich thick-disc stars with cool kinematics in the generally turbulent context of the primitive Galactic disc. This is discussed in the framework of the different disc formation scenarios.

We investigate, using the Gaia-ESO Survey internal Data-Release 2, the properties of the double sequence of the Milky Way discs (defined chemically as the high-alpha and low-alpha populations), and discuss their compatibility with discs defined by other means such as metallicity, kinematics or positions. This investigation uses two different approaches: in velocity space for stars located in the extended Solar neighbourhood, and in chemical space for stars at different ranges of Galactocentric radii and heights from the plane. The separation we find in velocity space allows us to investigate, in a novel manner, the extent in metallicity of each of the two sequences, identifying them with the two discs, without making any assumption about the shape of their metallicity distribution functions. Then, using the separation in chemical space, we characterise the spatial variation of the slopes of the [alpha/Fe] - [Fe/H] sequences for the thick and thin discs and the way in which the relative proportions of the two discs change across the Galaxy. We find that the thick disc (high-alpha sequence), extends up to [Fe/H]~ +0.2 and the thin disc (low-alpha sequence), at least down to [Fe/H]~ -0.8. Radial and vertical gradients in alpha-abundances are found for the thin disc, with mild spatial variations in its [alpha/Fe] - [Fe/H] paths, whereas for the thick disc we do not detect any such spatial variations. The small variations in the spatial [alpha/Fe] - [Fe/H] paths of the thin disc do not allow us to distinguish between formation models of this structure. On the other hand, the lack of radial gradients and [alpha/Fe] - [Fe/H] variations for the thick disc indicate that the mechanism responsible for the mixing of the metals in the young Galaxy (e.g. radial stellar migration or turbulent gaseous disc) was more efficient before the (present) thin disc started forming.

We present the identification of potential members of nearby Galactic globular clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 (RAVE-DR4) survey database. Our identifications are based on three globular clusters -- NGC 3201, NGC 5139 ($\omega$ Cen) and NGC 362 -- all of which are shown to have |RV|>100 km/s. The identification of globular cluster stars in RAVE DR4 data offers a unique opportunity to test the precision and accuracy of the stellar parameters determined with the currently available Stellar Parameter Pipelines (SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 and $\omega$ Cen, there is compelling evidence for numerous members (> 10) in the RAVE database; in the case of NGC 362 the evidence is more ambiguous, and there may be significant foreground and/or background contamination in our kinematically-selected sample. A comparison of the RAVE-derived stellar parameters and abundances with published values for each cluster and with BASTI isochrones for ages and metallicities from the literature reveals overall good agreement, with the exception of the apparent underestimation of surface gravities for giants, in particular for the most metal-poor stars. Moreover, if the selected members are part of the main body of each cluster our results would also suggest that the distances from Binney et al. 2013, where only isochrones more metal-rich than -0.9 dex were used, are typically underestimated by ~ 40% with respect to the published distances for the clusters, while the distances from Zwitter et al. 2010 show stars ranging from 1 to ~ 6.5 kpc -- with indications of a trend toward higher distances at lower metallicities -- for the three clusters analysed in this study.

According to our current cosmological model, galaxies like the Milky Way are expected to experience many mergers over their lifetimes. The most massive of the merging galaxies will be dragged towards the disc-plane, depositing stars and dark matter into an accreted disc structure. In this work, we utilize the chemo-dynamical template developed in Ruchti et al. to hunt for accreted stars. We apply the template to a sample of 4,675 stars in the third internal data release from the Gaia-ESO Spectroscopic Survey. We find a significant component of accreted halo stars, but find no evidence of an accreted disc component. This suggests that the Milky Way has had a rather quiescent merger history since its disc formed some 8-10 billion years ago and therefore possesses no significant dark matter disc.

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.

Large spectroscopic surveys have enabled in the recent years the computation of three-dimensional interstellar extinction maps thanks to accurate stellar atmospheric parameters and line-of-sight distances. Such maps are complementary to 3D maps extracted from photometry, allowing a more thorough study of the dust properties. Our goal is to use the high-resolution spectroscopic survey Gaia-ESO in order to obtain with a good distance resolution the interstellar extinction and its dependency as a function of the environment and the Galactocentric position. We use the stellar atmospheric parameters of more than 5000 stars, obtained from the Gaia-ESO survey second internal data release, and combine them with optical (SDSS) and near-infrared (VISTA) photometry as well as different sets of theoretical stellar isochrones, in order to calculate line-of-sight extinction and distances. The extinction coefficients are then compared with the literature to discuss their dependancy on the stellar parameters and position in the Galaxy. Within the errors of our method, our work does not show that there is any dependence of the interstellar extinction coefficient on the atmospheric parameters of the stars. We do not find any evidence of the variation of E(J-H)/E(J-K) with the angle from the Galactic centre nor with Galactocentric distance. This suggests that we are dealing with a uniform extinction law in the SDSS ugriz bands and the near-IR JHKs bands. Therefore, extinction maps using mean colour-excesses and assuming a constant extinction coefficient can be used without introducing any systematic errors.

Using the RAdial Velocity Experiment fourth data release (RAVE DR4), and a new metallicity calibration that will be also taken into account in the future RAVE DR5, we investigate the existence and the properties of super-solar metallicity stars ([M/H] > +0.1 dex) in the sample, and in particular in the Solar neighbourhood. We find that RAVE is rich in super-solar metallicity stars, and that the local metallicity distribution function declines remarkably slowly up to +0.4 dex. Our results show that the kinematics and height distributions of the super-solar metallicity stars are identical to those of the [M/H] < 0 thin-disc giants that we presume were locally manufactured. The eccentricities of the super-solar metallicity stars indicate that half of them are on a roughly circular orbit (e < 0.15), so under the assumption that the metallicity of the interstellar medium at a given radius never decreases with time, they must have increased their angular momenta by scattering at corotation resonances of spiral arms from regions far inside the Solar annulus. The likelihood that a star will migrate radially does not seem to decrease significantly with increasing amplitude of vertical oscillations within range of oscillation amplitudes encountered in the disc.

The RAdial Velocity Experiment (RAVE) published in November 2013 its fourth data release with the stellar atmospheric parameters, abundances, distances, radial velocities, proper motions and spectral morphological flags of more than 480,000 targets. With that, a plethora of papers ranging from the mass of the Milky Way to the mapping of the Diffuse Interstellar Band, and from the chemo-dynamical history and properties of the disc to the Galaxy's bar pattern speed have also been published. Being one of the largest spectroscopic surveys in the magnitude range of Gaia, RAVE has helped to pave the way for the exploitation of the Gaia catalogs. Here, we review some of these results and present some perspectives about future RAVE data releases.

We aim to characterize high-velocity (HiVel) stars in the solar vicinity both chemically and kinematically using the fourth data release of the RAdial Velocity Experiment (RAVE). We used a sample of 57 HiVel stars with Galactic rest-frame velocities larger than 275 km s$^{-1}$. With 6D position and velocity information, we integrated the orbits of the HiVel stars and found that, on average, they reach out to 13 kpc from the Galactic plane and have relatively eccentric orbits consistent with the Galactic halo. Using the stellar parameters and [$\alpha$/Fe] estimates from RAVE, we found the metallicity distribution of the HiVel stars peak at [M/H] = -1.2 dex and is chemically consistent with the inner halo. There are a few notable exceptions that include a hypervelocity star (HVS) candidate, an extremely high-velocity bound halo star, and one star that is kinematically consistent with the halo but chemically consistent with the disk. High-resolution spectra were obtained for the metal-rich HiVel star candidate and the second highest velocity star in the sample. Using these high-resolution data, we report the discovery of a metal-rich halo star that has likely been dynamically ejected into the halo from the Galactic thick disk. This discovery could aid in explaining the assembly of the most metal-rich component of the Galactic halo.

We performed a detailed study of the ratio of low-{\alpha} to high-{\alpha} stars in the Galactic halo as observed by the Gaia-ESO Survey. Using a sample of 381 metal-poor stars from the second internal data release, we found that the value of this ratio did not show evidence of systematic trends as a function of metallicity, surface gravity, Galactic latitude, Galactic longitude, height above the Galactic plane, and Galactocentric radius. We conclude that the {\alpha}-poor/{\alpha}-rich value of 0.28 $\pm$ 0.08 suggests that in the inner halo, the larger portion of stars were formed in a high star formation rate environment, and about 15% of the metal-poor stars originated from much lower star formation rate environments.

We determine the Galactic potential in the solar neigbourhood from RAVE observations. We select red clump stars for which accurate distances, radial velocities, and metallicities have been measured. Combined with data from the 2MASS and UCAC catalogues, we build a sample of 4600 red clump stars within a cylinder of 500 pc radius oriented in the direction of the South Galactic Pole, in the range of 200 pc to 2000 pc distances. We deduce the vertical force and the total mass density distribution up to 2 kpc away from the Galactic plane by fitting a distribution function depending explicitly on three isolating integrals of the motion in a separable potential locally representing the Galactic one with four free parameters. Because of the deep extension of our sample, we can determine nearly independently the dark matter mass density and the baryonic disc surface mass density. We find (i) at 1kpc Kz/(2piG) = 68.5 pm 1.0 Msun/pc2, and (ii) at 2 kpc Kz/(2piG) = 96.9 pm 2.2 Msun/pc2. Assuming the solar Galactic radius at R0 = 8.5 kpc, we deduce the local dark matter density rhoDM (z=0) = 0.0143 pm 0.0011Msun pc3 = 0.542 pm 0.042 Gev/cm3 and the baryonic surface mass density Sigma = 44.4 pm 4.1 Msun/pc2 . Our results are in agreement with previously published Kz determinations up to 1 kpc, while the extension to 2 kpc shows some evidence for an unexpectedly large amount of dark matter. A flattening of the dark halo of order 0.8 can produce such a high local density in combination with a circular velocity of 240 km/s . Another explanation, allowing for a lower circular velocity, could be the presence of a secondary dark component, a very thick disc resulting either from the deposit of dark matter from the accretion of multiple small dwarf galaxies, or from the presence of an effective phantom thick disc in the context of effective galactic-scale modifications of gravity.

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.

Most high-resolution spectroscopic studies of the Galactic discs were mostly confined to objects in the solar vicinity. Here we aim at enlarging the volume in which individual chemical abundances are used to characterise both discs, using the first internal data release of the Gaia-ESO survey. We derive and discuss the abundances of eight elements (Mg, Al, Si, Ca, Ti, Fe, Cr, Ni, and Y). The trends of these elemental abundances with iron are very similar to those in the solar neighbourhood. We find a natural division between alpha-rich and alpha-poor stars, best seen in the bimodality of the [Mg/M] distributions in bins of metallicity, which we attribute to thick- and thin-disc sequences, respectively. With the possible exception of Al, the observed dispersion around the trends is well described by the expected errors, leaving little room for astrophysical dispersion. Using previously derived distances from Recio-Blanco et al. (2014b), we further find that the thick-disc is more extended vertically and is more centrally concentrated towards the inner Galaxy than the thin-disc, which indicates a shorter scale-length. We derive the radial and vertical gradients in metallicity, iron, four alpha-element abundances, and Al for the two populations, taking into account the identified correlation between R_GC and |Z|. Radial metallicity gradient is found in the thin disc. The positive radial individual [alpha/M] gradients found are at variance from the gradients observed in the RAVE survey. The thin disc also hosts a negative vertical metallicity gradient, accompanied by positive individual [alpha/M] and [Al/M] gradients. The thick-disc, presents no radial metallicity gradient, a shallower vertical metallicity gradient than the thin-disc, an alpha-elements-to-iron radial gradient in the opposite sense than that of the thin disc, and positive vertical individual [alpha/M] and [Al/M] gradients.

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 APASS photometry in the Landolt BV and Sloan g'r'i' bands for all the 425,743 stars included in the latest 4th RAVE Data Release. The internal accuracy of the APASS photometry of RAVE stars, expressed as error of the mean of data obtained and separately calibrated over a median of 4 distinct observing epochs and distributed between 2009 and 2013, is 0.013, 0.012, 0.012, 0.014 and 0.021 mag for B, V, g', r' and i' band, respectively. The equally high external accuracy of APASS photometry has been verified on secondary Landolt and Sloan photometric standard stars not involved in the APASS calibration process, and on a large body of literature data on field and cluster stars, confirming the absence of offsets and trends. Compared with the Carlsberg Meridian Catalog (CMC-15), APASS astrometry of RAVE stars is accurate to a median value of 0.098 arcsec. Brightness distribution functions for the RAVE stars have been derived in all bands. APASS photometry of RAVE stars, augmented by 2MASS JHK infrared data, has been chi2 fitted to a densely populated synthetic photometric library designed to widely explore in temperature, surface gravity, metallicity and reddening. Resulting Teff and E(B-V), computed over a range of options, are provided and discussed, and will be kept updated in response to future APASS and RAVE data releases. In the process it is found that the reddening caused by an homogeneous slab of dust, extending for 140 pc on either side of the Galactic plane and responsible for E(B-V,poles)=0.036 +/- 0.002 at the galactic poles, is a suitable approximation of the actual reddening encountered at Galactic latitudes |b|>=25 deg.

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~{\AA}ngstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

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 investigate the kinematic parameters of the Milky Way disc using the RAVE and GCS stellar surveys. We do this by fitting a kinematic model to the data taking the selection function of the data into account. For stars in the GCS we use all phase-space coordinates, but for RAVE stars we use only $(l,b,v_{\rm los})$. Using MCMC technique, we investigate the full posterior distributions of the parameters given the data. We investigate the `age-velocity dispersion' relation for the three kinematic components ($\sigma_R,\sigma_{\phi},\sigma_z$), the radial dependence of the velocity dispersions, the Solar peculiar motion ($U_{\odot},V_{\odot}, W_{\odot} $), the circular speed $\Theta_0$ at the Sun and the fall of mean azimuthal motion with height above the mid-plane. We confirm that the Besan\c{c}on-style Gaussian model accurately fits the GCS data, but fails to match the details of the more spatially extended RAVE survey. In particular, the Shu distribution function (DF) handles non-circular orbits more accurately and provides a better fit to the kinematic data. The Gaussian distribution function not only fits the data poorly but systematically underestimates the fall of velocity dispersion with radius. We find that correlations exist between a number of parameters, which highlights the importance of doing joint fits. The large size of the RAVE survey, allows us to get precise values for most parameters. However, large systematic uncertainties remain, especially in $V_{\odot}$ and $\Theta_0$. We find that, for an extended sample of stars, $\Theta_0$ is underestimated by as much as $10\%$ if the vertical dependence of the mean azimuthal motion is neglected. Using a simple model for vertical dependence of kinematics, we find that it is possible to match the Sgr A* proper motion without any need for $V_{\odot}$ being larger than that estimated locally by surveys like GCS.

We study the relationship between age, metallicity, and alpha-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpc to 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages > 9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more alpha-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

(Abridged) We have used the atmospheric parameters, [alpha/Fe] abundances and radial velocities, determined from the Gaia-ESO Survey GIRAFFE spectra of FGK-type stars (iDR1), to provide a chemo-kinematical characterisation of the disc stellar populations. We focuss on a subsample of 1016 stars with high quality parameters, covering the volume |Z|<4.5kpc and R in the range 2-13kpc. We have identified a thin to thick disc separation in the [alpha/Fe] vs [M/H] plane, thanks to the presence of a low-density region in the number density distribution. The thick disc stars seem to lie in progressively thinner layers above the Galactic plane, as metallicity increases and [alpha/Fe] decreases. The thin disc population presents a constant value of the mean distance to the plane at all metallicities. Our data confirm the already known correlations between V_phi and [M/H] for the two discs. For the thick disc sequence, a study of the possible contamination by thin disc stars suggests a gradient up to 64km/s/dex. The distributions of V_phi, V_Z, and orbital parameters are analysed for the chemically separated samples. Concerning the gradients with galactocentric radius, we find for the thin disc a flat behaviour of V_phi, a [M/H] gradient of -0.058dex/kpc and a small positive [alpha/Fe] gradient. For the thick disc, flat gradients in [M/H] and [alpha/Fe] are derived. Our chemo-kinematical analysis suggests a picture in which the thick disc seems to have experienced a settling process, during which its rotation increased progressively, and, possibly, the V_phi dispersion decreased. At [M/H]-0.25dex and [alpha/Fe]0.1dex, the mean characteristics of the thick disc in distance to the Galactic plane, V_phi, V_phi dispersion and eccentricity agree with those of the thin disc stars, suggesting a possible connection between these populations at a certain epoch of the disc evolution.

We analyse the kinematics of ~400000 RAVE stars. We split the sample into hot and cold dwarfs, red-clump and non-clump giants. The kinematics of the clump giants are consistent with being identical with those of non-clump giants. We fit Gaussian velocity ellipsoids to the meridional-plane components of velocity of each star class and give formulae from which the shape and orientation of the velocity ellipsoid can be determined at any location. The data are consistent with the giants and the cool dwarfs sharing the same velocity ellipsoids; sigma_z rises from 21 kms in the plane to sim 55 kms at |z|=2 kpc, while sigma_r rises from 37 kms to 82 kms. At (R,z) the longest axis of one of these velocity ellipsoids is inclined to the Galactic plane by an angle ~0.8 arctan(z/R). We use a novel formula to obtain precise fits to the highly non-Gaussian distributions of v_phi components. We compare the observed velocity distributions with the predictions of a dynamical model fitted to the velocities of stars that lie within ~150 pc of the Sun and star counts towards the Galactic pole. The model accurately reproduces the non-Gaussian nature of the v_r and v_z distributions and provides excellent fits to the data for v_z at all locations. The model v_phi distributions for the cool dwarfs fit the data extremely well, while those for the hot dwarfs have displacements to low v_phi that grow with |z| from very small values near the plane. At |z|>0.5 kpc, the theoretical v_phi distributions for giants show a deficit of stars with large v_phi and the model v_r distributions are too narrow. Systematically over-estimating distances by 20 per cent introduces asymmetry into the model v_r and v_z distributions near the plane and but significantly improves the fits to the data at |z|>0.5 kpc. The quality of the fits lends credence to the assumed, disc-dominated, gravitational potential.

Probability density functions are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N>=10. Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature absolute-magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schoenrich, Binney & Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants <pi> tends to be larger than the true distance by up to 30 percent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are over-estimating distances to giants, especially in young clusters.

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.

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

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.

By selecting in the RAVE-DR4 survey the stars located between 1 and 2 kpc above the Galactic plane, we question the consistency of the simplest three-component model (thin disc, thick disc, halo) for the Milky Way. We confirm that the metallicity and azimuthal velocity distribution functions of the thick disc are not Gaussian. In particular, we find that the thick disc has an extended metallicity tail going at least down to [M/H]=-2 dex, contributing roughly at 3% of the entire thick disc population and having a shorter scale-length compared to the canonical thick disc. The mean azimuthal velocity of these metal-poor stars allows us to estimate the correlation between the metallicity and the orbital velocity, which is an important constraint on the formation mechanisms of the Galactic thick disc. Given our simple approach, we find dVphi/d[M/H] ~ 50 km/s/dex, which is in a very good agreement with previous literature values. We complete the study with a brief discussion on the implications of the formation scenarios for the thick disc, and suggest that given the above mentioned characteristics, a thick disc mainly formed by radial migration mechanisms seems unlikely.

Context. Galactic open clusters (OCs) mainly belong to the young stellar population in the Milky Way disk, but are there groups and complexes of OCs that possibly define an additional level in hierarchical star formation? Current compilations are too incomplete to address this question, especially regarding radial velocities (RVs) and metallicities ($[M/H]$). Aims. Here we provide and discuss newly obtained RV and $[M/H]$ data, which will enable us to reinvestigate potential groupings of open clusters and associations. Methods.We extracted additional RVs and $[M/H]$ from the RAdial Velocity Experiment (RAVE) via a cross-match with the Catalogue of Stars in Open Cluster Areas (CSOCA). For the identified OCs in RAVE we derived RV and $[M/H]$ from a cleaned working sample and compared the results with previous findings. Results. Although our RAVE sample does not show the same accuracy as the entire survey, we were able to derive reliable RV for 110 Galactic open clusters. For 37 OCs we publish RV for the first time. Moreover, we determined $[M/H]$ for 81 open clusters, extending the number of OCs with $[M/H]$ by 69.

RAVE, the unbiased magnitude limited survey of the southern sky stars, contained 456,676 medium-resolution spectra at the time of our analysis. Spectra cover the CaII IRT range which is a known indicator of chromospheric activity. Our previous work (Matijevi\v{c} et al. 2012) classified all spectra using locally linear embedding. It identified 53,347 cases with a suggested emission component in calcium lines. Here we use a spectral subtraction technique to measure the properties of this emission. Synthetic templates are replaced by the observed spectra of non-active stars to bypass the difficult computations of non-LTE profiles of the line cores and stellar parameter dependence. We derive both the equivalent width of the excess emission for each calcium line on a 5\AA\ wide interval and their sum EW_IRT for ~44,000 candidate active dwarf stars with S/N>20 and with no respect to the source of their emission flux. From these ~14,000 show a detectable chromospheric flux with at least 2\sigma\ confidence level. Our set of active stars vastly enlarges previously known samples. Atmospheric parameters and in some cases radial velocities of active stars derived from automatic pipeline suffer from systematic shifts due to their shallower calcium lines. We re-estimate the effective temperature, metallicity and radial velocities for candidate active stars. The overall distribution of activity levels shows a bimodal shape, with the first peak coinciding with non-active stars and the second with the pre main-sequence cases. The catalogue will be publicly available with the next RAVE public data releases.

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.

The RAVE survey, combined with proper motions and distance estimates, can be used to study in detail stellar kinematics in the extended solar neighbourhood (solar suburb). Using the red clump, we examine the mean velocity components in 3D between an R of 6 and 10 kpc and a Z of -2 to 2 kpc, concentrating on North-South differences. Simple parametric fits to the R, Z trends for VPHI and the velocity dispersions are presented. We confirm the recently discovered gradient in mean Galactocentric radial velocity, VR, finding that the gradient is more marked below the plane, with a Z gradient also present. The vertical velocity, VZ, also shows clear structure, with indications of a rarefaction-compression pattern, suggestive of wave-like behaviour. We perform a rigorous error analysis, tracing sources of both systematic and random errors. We confirm the North-South differences in VR and VZ along the line-of-sight, with the VR estimated independent of the proper motions. The complex three-dimensional structure of velocity space presents challenges for future modelling of the Galactic disk, with the Galactic bar, spiral arms and excitation of wave-like structures all probably playing a role.

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.

Diffuse interstellar bands are usually observed in spectra of hot stars, where interstellar lines are rarely blended with stellar ones. The need for hot stars is a strong limitation in the number of sightlines we can observe and the distribution of sightlines in the Galaxy, as hot stars are rare and concentrated in the Galactic plane. We are introducing a new method, where interstellar lines can be observed in spectra of cool stars in large spectroscopic surveys. The method is completely automated and does not require prior knowledge of the stellar parameters. If known, the stellar parameters only reduce the computational time and are not involved in the extraction of the interstellar spectrum. The main step in extracting interstellar lines is a construction of the stellar spectrum, which is in our method done by finding other observed spectra that lack interstellar features and are otherwise very similar to the spectrum in question. Such spectra are then combined into a single stellar spectrum template, which matches the stellar component in an observed spectrum. We demonstrate the performance of this new method on a sample of 482,430 spectra observed in RAVE survey. However, many spectra have to be combined (48 on average) in order to achieve a S/N ratio high enough to measure the DIB's profile, hence limiting the spatial information about the ISM. Only one strong interstellar line is included in the RAVE spectral range, a diffuse interstellar band at 8620 \AA. We compare its equivalent width with extinction maps and with Bayesian reddening, calculated for individual stars, and provide a linear relation between the equivalent width and reddening. Separately from the introduced method, we calculate equivalent widths of the diffuse interstellar band in spectra of hot stars with known extinction and compare all three linear relations with each other and with relations from the literature.

We aim to understand the accretion history of the Milky Way by exploring the vertical and radial properties of the Galactic thick disc. We study the chemical and kinematic properties of roughly a thousand spectra of faint magnitude foreground Galactic stars observed serendipitously during extra-galactic surveys in four lines-of-sight: three in the southern Galactic hemisphere (surveys of the Carina, Fornax and Sculptor dwarf spheroidal galaxies) and one in the northern Galactic hemisphere (a survey of the Sextans dwarf spheroidal galaxy). The foreground stars span distances up to 3kpc from the Galactic plane and Galactocentric radii up to 11kpc. Only three lines-of-sight have a sufficient number of foreground stars for a robust analysis. Towards Sextans in the Northern Galactic hemisphere and Sculptor in the South, we measure a consistent decrease in mean metallicity with height from the Galactic plane, suggesting a chemically symmetric thick disc. This decrease can either be due to an intrinsic thick disc metallicity gradient, or simply due to a change in the thin disc/thick disc population ratio and no intrinsic metallicity gradients for the thick disc. We favour the latter explanation. In contrast, we find evidence of an unpredicted metal-poor population in the direction of Carina. This population was earlier detected by Wyse et al. (2006), but our more detailed analysis provides robust estimates of its location (|Z|<1kpc), metallicity (-2<[M/H]<-1 dex) and azimuthal orbital velocity (V_phi ~120 km/s). Given the chemo-dynamical properties of the over-density towards the Carina line-of-sight, we suggest that it represents the metal-poor tail of the canonical thick disc. In spite of the small number of stars available, we suggest that this metal-weak thick disc follows the often suggested thick disc velocity-metallicity correlation of dV_phi/d[M/H]~40-50km/s/dex.

Aims: We study the relations between stellar kinematics and chemical abundances of a large sample of RAVE giants in search for selection criteria needed for disentangling different Galactic stellar populations. Methods: We select a sample of 2167 giant stars with signal-to-noise per spectral measurements above 75 from the RAVE chemical catalogue and follow the analysis performed by Gratton and colleagues on 150 subdwarf stars spectroscopically observed at high-resolution. We then use a larger sample of 9131 giants (with signal-to-noise above 60) to investigate the chemo-kinematical characteristics of our stars by grouping them into nine subsamples with common eccentricity ($e$) and maximum distance achieved above the Galactic plane ($Z_max$). Results: The RAVE kinematical and chemical data proved to be reliable by reproducing the results by Gratton et al. obtained with high-resolution spectroscopic data. Our analysis, based on the $e$-$Z_max$ plane combined with additional orbital parameters and chemical information, provides an alternative way of identifying different populations of stars. In addition to extracting canonical thick- and thin-disc samples, we find a group of stars in the Galactic plane ($Z_max<1$ kpc and 0.4 $< e < $0.6), which show homogeneous kinematics but differ in their chemical properties. We interpret this as a clear sign that some of these stars have experienced the effects of heating and/or radial migration, which have modified their original orbits. The accretion origin of such stars cannot be excluded.

We use kinematical and chemical properties of 754 Corot stars to characterise the stellar populations of the Milky Way disc in three beams close the Galactic plane. From the atmospheric parameters derived in Gazzano et al. (2010) with the Matisse algorithm, we derived stellar distances using isochrones. Combining these data with proper motions, we provide the complete kinematical description of stars in three Corot fields. Finally, we used kinematical criteria to identify the Galactic populations in our sample and study their characteristics, particularly their chemistry. Comparing our kinematics with the Besancon Galactic model, we show that, within 3-sigma, simulated and observed kinematical distributions are in good agreement. We study the characteristics of the thin disc, finding a correlation that is significant at a value of 2-sigma between the V-velocity component and the metallicity for two different radial distance bins (8-9kpc and 9-10kpc; but not for the most inner bin 7-8kpc, probably because of the uncertainties in the abundances) which could be interpreted as radial migration evidence. We also measured a radial metallicity gradient value of -0.097+/-0.015dex/kpc with giant stars, and -0.053+/-0.015dex/kpc with dwarfs. Finally, we identified metal-rich stars with peculiar high [alpha/Fe] values in the directions pointing to the inner part of the Galaxy. Applying the same methodology to the planet-hosting stars detected by Corot shows that they mainly belong to the thin disc population with normal chemical and kinematical properties.

The Gaia mission, with its unprecedented astrometric and photometric precision, combined with its Radial Velocity Spectrometer, will provide to the astronomical community a wealth of necessary constraints to disentangle between the different formation scenarios of the Galactic thick disc. The aim of this review is to present some of the recent results obtained spectroscopically concerning this Galactic structure, and highlight the open questions that still remain to be answered under the Gaia era. These concern mainly the measurement of the chemo-dynamical properties of the Milky Way at the inner and outer parts, which allow us to determine the total accreted mass from the mergers with satellite galaxies, and will give us an estimate of the strength of the radial migration phenomena to form such a structure.

Our ability to extract information from the spectra of stars depends on reliable models of stellar atmospheres and appropriate techniques for spectral synthesis. Various model codes and strategies for the analysis of stellar spectra are available today. We aim to compare the results of deriving stellar parameters using different atmosphere models and different analysis strategies. The focus is set on high-resolution spectroscopy of cool giant stars. Spectra representing four cool giant stars were made available to various groups and individuals working in the area of spectral synthesis, asking them to derive stellar parameters from the data provided. The results were discussed at a workshop in Vienna in 2010. Most of the major codes currently used in the astronomical community for analyses of stellar spectra were included in this experiment. We present the results from the different groups, as well as an additional experiment comparing the synthetic spectra produced by various codes for a given set of stellar parameters. Similarities and differences of the results are discussed. Several valid approaches to analyze a given spectrum of a star result in quite a wide range of solutions. The main causes for the differences in parameters derived by different groups seem to lie in the physical input data and in the details of the analysis method. This clearly shows how far from a definitive abundance analysis we still are.

We performed a spectroscopic survey of nearly 700 stars probing the galactic thick disc far from the solar neighbourhood towards the galactic coordinates (l~277, b~47). The derived effective temperatures, surface gravities and overall metallicities were then combined with stellar evolution isochrones, radial velocities and proper motions to derive the distances, kinematics and orbital parameters of the sample stars. The targets belonging to each galactic component (thin disc, thick disc, halo) were selected either on their kinematics or according to their position above the galactic plane, and the vertical gradients were also estimated. We present here atmospheric parameters, distances and kinematics for this sample, and a comparison of our kinematic and metallicity distributions with the Besancon model of the Milky Way. The thick disc far from the solar neighbourhood is found to differ only slightly from the thick disc properties as derived in the solar vicinity. For regions where the thick disc dominates, we measured vertical velocity and metallicity trends of d(V_phi)/dZ = 19 +/- 8 km/s/kpc and d[M/H]/dZ = -0.14 +/- 0.05 dex/kpc, respectively. These trends can be explained as a smooth transition between the different galactic components, although intrinsic gradients could not be excluded. In addition, a correlation d(V_phi)/d[M/H] = -45 +/- 12 km/s/dex between the orbital velocity and the metallicity of the thick disc is detected. This gradient is inconsistent with the SDSS photometric survey analysis, which did not detect any such trend, and challenges radial migration models of thick disc formation. Estimations of the scale heights and scale lengths for different metallicity bins of the thick disc result in consistent values, with hR~3.4 \pm 0.7 kpc, and hZ~694 \pm 45 pc, showing no evidence of relics of destroyed massive satellites.

(Abridged) Galactic archaeology aims to determine the evolution of the Galaxy from the chemical and kinematical properties of its stars. The analysis of current large spectroscopic surveys (thousands of stars) and future ones (millions of stars) require automated analysis techniques to obtain robust estimates of the stellar parameters. Several on-going and planned spectroscopic surveys have selected their wavelength region to contain the IR CaII triplet and this paper focuses on the automatic analysis of such spectra. We investigated two algorithms, MATISSE and DEGAS, both of which compare the observed spectrum to a grid of synthetic spectra, but each uses a different mathematical approach for finding the optimum match and hence the best stellar parameters. We identified degeneracies in different regions of the HR diagram: hot dwarfs and giants share the same spectral signatures. Furthermore, the surface gravity of cooler dwarfs is difficult to determine accurately. These effects are intensified when the information decreases (e.g. metal-poor stars or low SNR spectra). Our results show that the local projection method MATISSE is preferred for high SNR spectra, whereas the decision-tree method DEGAS is preferred for noisier spectra. We therefore propose a hybrid approach of both methods and show that sufficiently accurate results for the purposes of galactic archaeology are retrieved down to SNR~20 for typical thin or thick disc stars, and down to SNR~50 for the more metal-poor halo giants. If unappreciated, degeneracies in stellar parameters can introduce biases in derived quantities for target stars such as distances and full space motions. These biases can be minimised using the knowledge gained by thorough testing of the proposed algorithm, which in turn lead to robust automated methods for the coming extensive stellar spectroscopic surveys in the Local Group.