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.

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.

Recent observations suggest a "double-branch" behaviour of Li/H versus metallicity in the local thick and thin discs. This is reminiscent of the corresponding O/Fe versus Fe/H behaviour, which has been explained as resulting from radial migration in the Milky Way disc. We use a semi-analytical model of disc evolution with updated chemical yields and parameterised radial migration. We explore the cases of long-lived (red giants of a few Gy lifetime) and shorter-lived (Asymptotic Giant Branch stars of several 10$^8$ yr) stellar sources of Li, as well as those of low and high primordial Li. We show that both factors play a key role in the overall Li evolution. We find that the observed "two-branch" Li behaviour is only directly obtained in the case of long-lived stellar Li sources and low primordial Li. In all other cases, the data imply systematic Li depletion in stellar envelopes, thus no simple picture of the Li evolution can be obtained. This concerns also the reported Li/H decrease at supersolar metallicities.

The chemical evolution of lithium in the Milky Way represents a major problem in modern astrophysics. Indeed, lithium is, on the one hand, easily destroyed in stellar interiors, and, on the other hand, produced at some specific stellar evolutionary stages that are still not well constrained. The goal of this paper is to investigate the lithium stellar content of Milky Way stars in order to put constraints on the lithium chemical enrichment in our Galaxy, in particular in both the thin and thick discs. Thanks to high-resolution spectra from the ESO archive and high quality atmospheric parameters, we were able to build a massive and homogeneous catalogue of lithium abundances for 7300 stars derived with an automatic method coupling, a synthetic spectra grid, and a Gauss-Newton algorithm. We validated these lithium abundances with literature values, including those of the Gaia benchmark stars. In terms of lithium galactic evolution, we show that the interstellar lithium abundance increases with metallicity by 1 dex from [M/H]=-1 dex to +0.0 dex. Moreover, we find that this lithium ISM abundance decreases by about 0.5 dex at super-solar metalllicity. Based on a chemical separation, we also observed that the stellar lithium content in the thick disc increases rather slightly with metallicity, while the thin disc shows a steeper increase. The lithium abundance distribution of alpha-rich, metal-rich stars has a peak at A(Li)~3 dex. We conclude that the thick disc stars suffered of a low lithium chemical enrichment, showing lithium abundances rather close to the Spite plateau while the thin disc stars clearly show an increasing lithium chemical enrichment with the metallicity, probably thanks to the contribution of low-mass stars.

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

In the current era of large spectroscopic surveys of the Milky Way, reference stars for calibrating astrophysical parameters and chemical abundances are of paramount importance. We determine elemental abundances of Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Co and Ni for our predefined set of Gaia FGK benchmark stars. By analysing high-resolution and high-signal to noise spectra taken from several archive datasets, we combined results of eight different methods to determine abundances on a line-by-line basis. We perform a detailed homogeneous analysis of the systematic uncertainties, such as differential versus absolute abundance analysis, as well as we assess errors due to NLTE and the stellar parameters in our final abundances. Our results are provided by listing final abundances and the different sources of uncertainties, as well as line-by-line and method-by-method abundances. The Gaia FGK benchmark stars atmospheric parameters are already being widely used for calibration of several pipelines applied to different surveys. With the added reference abundances of 10 elements this set is very suitable to calibrate the chemical abundances obtained by these pipelines.

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.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

High-resolution spectroscopy is a powerful tool to study the dynamical structure of pulsating stars atmosphere. We aim at comparing the line asymmetry and velocity of the two delta Sct stars rho Pup and DX Cet with previous spectroscopic data obtained on classical Cepheids and beta Cep stars. We obtained, analysed and discuss HARPS high-resolution spectra of rho Pup and DX Cet. We derived the same physical quantities as used in previous studies, which are the first-moment radial velocities and the bi-Gaussian spectral line asymmetries. The identification of f=7.098 (1/d) as a fundamental radial mode and the very accurate Hipparcos parallax promote rho Pup as the best standard candle to test the period-luminosity relations of delta Sct stars. The action of small-amplitude nonradial modes can be seen as well-defined cycle-to-cycle variations in the radial velocity measurements of rho Pup. Using the spectral-line asymmetry method, we also found the centre-of-mass velocities of rho Pup and DX Cet, V_gamma = 47.49 +/- 0.07 km/s and V_gamma = 25.75 +/- 0.06 km/s, respectively. By comparing our results with previous HARPS observations of classical Cepheids and beta Cep stars, we confirm the linear relation between the atmospheric velocity gradient and the amplitude of the radial velocity curve, but only for amplitudes larger than 22.5 km/s. For lower values of the velocity amplitude (i.e., < 22.5 km/s), our data on rho Pup seem to indicate that the velocity gradient is null, but this result needs to be confirmed with additional data. We derived the Baade-Wesselink projection factor p = 1.36 +/- 0.02 for rho Pup and p = 1.39 +/- 0.02 for DX Cet. We successfully extended the period-projection factor relation from classical Cepheids to delta Scuti stars.

Aims. The Baade-Wesselink method of distance determination is based on the oscillations of pulsating stars. The key parameter of this method is the projection factor used to convert the radial velocity into the pulsation velocity. Our analysis was aimed at deriving for the first time the projection factor of delta Scuti stars, using high-resolution spectra of the high-amplitude pulsator AI Vel and of the fast rotator beta Cas. Methods. The geometric component of the projection factor (i.e. p0) was calculated using a limb-darkening model of the intensity distribution for AI Vel, and a fast-rotator model for beta Cas. Then, using SOPHIE/OHP data for beta Cas and HARPS/ESO data for AI Vel, we compared the radial velocity curves of several spectral lines forming at different levels in the atmosphere and derived the velocity gradient associated to the spectral-line-forming regions in the atmosphere of the star. This velocity gradient was used to derive a dynamical projection factor p. Results. We find a flat velocity gradient for both stars and finally p = p0 = 1.44 for AI Vel and p = p0 = 1.41 for beta Cas. By comparing Cepheids and delta Scuti stars, these results bring valuable insights into the dynamical structure of pulsating star atmospheres. They suggest that the period-projection factor relation derived for Cepheids is also applicable to delta Scuti stars pulsating in a dominant radial mode.