{'fraction_non_alphanumeric': 0.04605803925939342, 'fraction_numerical': 0.02399030357068223, 'mean_word_length': 4.404306580334287, 'pattern_counts': {'":': 0, '<': 37, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Gaia-ESO Survey: the Galactic Thick to Thin Disc transition

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


Introduction
Understanding the chemodynamical evolution of the Milky Way disc stellar populations is a crucial step to reconstruct the history of our Galaxy. This approach can shed new light on the detailed processes that led the Milky Way to form and evolve as a disc-type galaxy, in the more general context of galaxy evolution. In particular, since the first works highlighting the existence of the Milky Way's thin/thick disc dichotomy (Yoshii, 1982;Gilmore & Reid, 1983), the study of the structural and chemodynamical transition between the two disc populations has opened new pathways to constraint Galactic evolution models. External mechanisms, supported by the hierarchical formation of galaxies in the ΛCDM paradigm, have been invoked to explain the chemodynamical characteristics of thick disc stars. Among them, the accretion of dwarf galaxies (e.g. Statler, 1988;Abadi et al., 2003) and minor mergers of satellites (e.g. Quinn et al., 1993;Villalobos & Helmi, 2008) heating dynamically the disc, have been proposed. Similarly, Jones & Wyse (1983) and Brook et al. (2004Brook et al. ( , 2007 have proposed the accretion of a gasrich merger from the collapse of which the thick disc would have formed. On the other hand, purely internal formation mechanisms have also been proposed, like the early turbulent phase of the primordial disc (e.g. Bournaud et al., 2009) and the stars radial migration due to resonances with the spiral structure or the bar of the Milky Way (e.g. Schönrich & Binney, 2009c;Minchev & Famaey, 2010;Loebman et al., 2011). The relative importance of the different proposed physical processes of evolution can now be tested with increasing robustness, thanks to the Galactic spectroscopic surveys targeting the disc stars.
Spectroscopic studies of the disc populations were initially based on kinematically selected samples of solar neighbourhood stars (e.g. Prochaska et al., 2000;Fuhrmann, 2004;Bensby et al., 2005;Reddy et al., 2006), with some exceptions like Edvardsson et al. (1993) or Fuhrmann (1998). In the classical picture that emerged from those studies, the thick disc population is characterized, compared to the thin disc one, as kinematically hotter (Casetti-Dinescu et al., 2011), and therefore, composed by stars moving in Galactic orbits with a higher scale height (∼900 pc vs. ∼300 pc, e.g. Jurić et al., 2008). In addition, the thick disc metallicity distribution peaks at lower values (∼-0.5 dex compared to ∼-0.2 dex for the thin disc, e.g Wyse & Gilmore, 1995;Kordopatis et al., 2011;Lee et al., 2011b). With respect to the thin disc, the thick disc stars are older and have an enhanced ratio of α-element abundances over iron ([α/Fe], e.g. Bensby et al., 2005Bensby et al., , 2007Fuhrmann, 2008). Nevertheless, it progressively appeared that the distributions in the above mentioned physical parameters displayed significant overlaps between the thin and the thick disc populations (e.g. Bensby & Feltzing, 2010). This blurred our comprehension of the interplay between the discs and questionned the classical kinematically-based definitions (Bovy et al., 2012c).
On the other hand, as discussed in Bovy et al. (2012b), defining stellar populations by abundances patterns is a better approach than the traditional kinematical criteria, as chemical abundances can correlate with disc structure, but are formally independent of it. Moreover, the identification of possible chemical evolutionary paths and/or gaps in the abundance ratios distributions can be crucial to disentangle the otherwise overlapping populations in the kinematical space.
Based on observations collected with the FLAMES spectrograph at the VLT/UT2 telescope (Paranal Observatory, ESO, Chile), for the  In the recent years, the number of stars analysed with high enough spectroscopic resolution to provide detailed chemical diagnostics has increased from a few hundreds to several tens of thousands. The advent of Milky Way spectroscopic surveys and of automatised chemical analysis techniques have improved both the statistical robustness and the homogeneity of the data as a consequence. From the theoretical side, new approaches based, at least partially, on the chemical identification of disc sub-populations, as the chemical-tagging (Freeman & Bland-Hawthorn, 2002) or the mono-abundance populations  methods, are opening promising pathways for constraining the Milky Way's evolutionary processes. This is also challenging our previous conception of the Milky Way disc and the way in which stellar populations can be defined.
In this context, the low resolution (R ∼ 2 000) SEGUE spectroscopic survey addressed the question of the thin to thick disc transition with a robust statistical approach, inside and outside the solar neighbourhood. Lee et al. (2011b) analysed ∼ 17 300 G-type dwarfs with errors of about 0.23 dex in [Fe/H] (Smolinski et al., 2011) and 0.1 dex for [α/Fe] (Lee et al., 2011a). This was complemented by Schlesinger et al. (2012) who analysed 24 270 G and 16 847 K dwarfs at distances from 0.2 to 2.3 kpc from the Galactic plane. The Lee et al. (2011b) [α/Fe] versus [Fe/H] distribution of number densities (see their Fig. 2) already motivated the authors to chemically separate the thin and the thick disc populations, although no clear gap is observed between the two. The SEGUE data interpretations differ from authors proposing no thin-thick disc distinction (Bovy et al., 2012b) to those allowing the existence of a distinct thick-disc component formed through an external mechanism (Liu & van de Ven, 2012). Nevertheless, selection effects have been invoked to explain some of the distinctions between the Lee et al. (2011a) and the Bovy et al. (2012b) chemical distributions. On the other hand, the data interpretation in the sense of no thin-thick disc distinction agrees with the results of the Schönrich & Binney (2009a) theoretical study.
More recently, Boeche et al. (2013) analysed the chemokinematical characteristics of 9 131 giants included in the last RAVE survey data release, with available [M/H] and [α/Fe]. The RAVE survey data have a higher resolution (R=7 500) than SEGUE spectra, but with a wavelength range limited to the infrared Ca II triplet. The Boeche et al. (2013) analysis, based on the eccentricity-Zmax plane combined with additional orbital parameters, allowed them to identify three stellar populations that could be associated with the Galactic thin disc, a dissipative component composed mostly of thick-disc stars, and the halo. Nevertheless, their thin and thick disc populations, defined with the above mentioned dynamical criteria, do not clearly separate in their [α/Fe] v.s. [M/H] plane (c.f. their Fig. 10).
On the other hand, recent studies of solar neighbourhood samples, analysed with high resolution spectroscopic data, have revealed the existence of a gap in the [α/Fe] versus [Fe/H] plane (e.g Fuhrmann, 2004;Reddy et al., 2006;Bensby et al., 2007). Recently, Adibekyan et al. (2013) used the stellar sample of 1 111 long-lived FGK dwarf stars from Adibekyan et al. (2012) to separate and characterise the different Galactic stellar subsystems, taking into account the existence of the above mentioned gap. Their [α/Fe] measurements are based on Mg, Si and Ti averaged abundances and the typical relative uncertainties in the metallicity and [α/Fe] are of about 0.03 dex. Haywood et al. (2013) derived ages of the Adibekyan et al. (2012) sample, concluding that two regimes appear in the age-[α/Fe] plane, that they identify as the epochs of the thick and thin disc for-mation. In the Haywood et al. (2013) scenario, the thick disc "formed from a well mixed interstellar medium, probably first in starburst, then in a more quiescent mode, over a time scale of 4-5 Gyr". They also suggest that "the youngest thick disk set the initial conditions from which the inner thin disk started to form 8 Gyr ago, at [Fe/H] in the range of (-0.1,+0.1) dex and [α/Fe]=0.1 dex". Finally, they interpret the low-metallicity tail of the thin disc as having an outer disc origin, somewhat disconnected from the inner disc evolution. In addition, similar separations in the [α/Fe] versus [Fe/H] plane are distinguishable in the Ramírez et al. (2013) oxygen-abundance versus [Fe/H] measurements. These authors also report similar age-[α/Fe] regimes as in Haywood et al. (2013).
In summary, the increasing number of Galactic disc stars with known kinematical and chemical characteristics has emphasized the question of the continuity or not between the thick and the thin disc components, in structural properties, kinematics, abundances and ages. Different authors like Norris (1999) and Bovy et al. (2012b) have even raised the possibility of considering the Galactic disc as a single component, without a thick/thin disc distinction, analysing for instance the massweighted scale-height distribution of disc stars.
In this context, the GIRAFFE observations of the Gaia-ESO Survey (GES, Gilmore et al., 2012) offer a unique opportunity to extend the previously high-resolution studies of the solar neighbourhood to larger and more radially extended samples. In its first data release, GES already includes ∼10 000 spectra (R∼20 000 and R∼16 000) of stars in the Milky Way field, from the halo, the thick disc, the thin disc and the bulge. The expected final number of targets, at the end of the five years of observations, is ∼100 000.
The present paper addresses the thick-thin disc transition as seen in the GES kinematical and chemical ([α/Fe] and global metallicity) first release data of FGK-type stars. Our point of view, exploiting the high resolution of the GES data, is that of the chemical characterisation and definition of the disc populations (essentially identified as thin and thick discs). As a consequence, the subsequent analysis of the structural properties, kinematics and orbital parameters will be done in the framework of the chemically defined disc populations. Section 2 describes the data sample, including the stellar parameters and [α/Fe] measurements. Section 3 presents the derivation of the stellar distances and kinematics. Section 4 focuses on the analysis of the number density distribution in the [α/Fe] versus [M/H] plane. Then, the distributions of the distances to the Galactic plane (Section 5), of the azimuthal velocities (Section 6), the azimuthal and vertical velocity dispersions (Section 7) and the stellar orbital parameters (Section 8) are discussed, in the context of a proposed thick-thin disc separation based on the chemical criteria. Section 9 presents the rotational and abundance gradients with galactocentric radius and, finally, we discuss our results in Section 10.

The GES data sample, the stellar parameters and the [α/Fe] measurements.
The present work is based on the data collected by the Gaia-ESO Survey during the first nine months of observations. As explained in Gilmore et al. (2012), the GES consortium has a work package structure that manages the data flows, from target selection, through data reduction, spectrum analysis, astrophysical parameters determination, calibration and homogenisation, to delivery of science data for verification analysis. The general data processing is described in Gilmore et al. (2014, in preparation). In this paper, we analyse the iDR1 GES results for ∼5 000 Milky Way disc field stars observed with the GIRAFFE spectrograph: 4 534 stars observed with both the HR10 (R∼19 800) and the HR21 (R∼16 200) setup modes, and 394 stars with only HR10 spectra. Fig. 1 shows the l and b coordinates for the lines of sight of the GES iDR1 (bulge and star clusters fields are excluded). All the targets were selected from Visible and Infrared Survey Telescope for Astronomy (VISTA) photometry. As explained in Gilmore et al. (2012), the target selection for the Halo and disc targets is based on the Disc/Halo transition seen in the Sloan Digital Sky Survey photometry at 17 < r < 18 and 0.2 < g−r <0.4. The equivalent selection from VISTA near-infrared photometry is used by GES with two main selection boxes: -Blue box: 0.0< J-K <0.45 and 14.0 < J < 17.5 mag -Red box: 0.4< J-K <0.70 and 12.5 < J < 15.0 mag The colour boxes are shifted according to the Schlegel extinction maps (Schlegel et al., 1998a). In addition, when the stellar density was not enough to fill in all the FLAMES fibers, as for high latitude fields, the red box was extended to around J-K <0.85 and J<17 mag to allow for second priority fibers. The spectrum analysis has shown that the blue box has mainly selected dwarf stars with T eff between 6500 K and 5250 K. The corresponding metallicity distribution peaks at around -0.6 dex for high latitude fields (|b| >30) and at -0.3 dex for low latitude fields. On the other hand, the red box has targeted dwarfs with T eff between 6000 K and 4500 K and a small number of giant stars. The metallicity distribution of the red box targets is quite independent on the latitude and it peaks around -0.4 dex. Finally, the extention of the selection boxes in the low density fields has added main sequence dwarfs with T eff between 6000 K and 4000 K with a doubled peaked metallicity distribution (peaks at ∼-0.8 dex and -0.3 dex). The signal to noise value of this second priority targets is in general lower than the main targets data. As a consequence of the selection criteria, the thin disc is preferentially targeted in the low latitude fields with respect to the thick disc, in a proportion of about 3:1, against about 1:1 for the high latitude fields.
The radial velocity measurements are based on a spectral fitting technique, that includes a first guess determination of atmospheric parameters. The method is described in detail by Koposov (2014, in preparation). Typical errors in radial velocity are of the order of 0.3 km/s. The stellar atmospheric parameters and [α/Fe] abundances were determined by the GES work package in charge of the GIRAFFE spectrum analysis for FGK-type stars. The final recommended effective temperatures (T eff ), surface gravities (log g), global metallicities ([M/H]) and [α/Fe] abundances ([α/Fe]) were estimated thanks to the combination of three different procedures: MATISSE (Recio-Blanco et al., 2006), FERRE (Allende , and further developments) and SME (Valenti & Piskunov, 1996, and further developments). The relative error distributions peak at 70 K for T eff , 0.10 dex for log g, 0.08 dex for [M/H] and 0.03 dex for [α/Fe]. The homogeneity of the parameters between the stars observed with HR10 and HR21 and those observed with HR10 only was verified during the GES parameters validation process. The adopted reference solar abundances are those of Grevesse et al. (2007) and the Gustafsson et al. (2008) MARCS model atmospheres were used for the analysis. For more details about the related GES parametrisation pipeline we refer the reader to Recio-Blanco et al. (2014). In addition, 95% of the targets have proper motions available from the PPMXL catalogue (Roeser et al., 2010) with typical errors of 8 mas/yr.  To enhance the best quality data, the points size is inversely proportional to the associated errors in the atmospheric parameters.

The global [α/Fe] estimations
The GES [α/Fe] estimations used in this paper were determined, at the same time as the stellar atmospheric parameters, by the three above mentioned groups with different methodologies . More particularly, the final GES iDR1 [α/Fe] values are the mean of the three independent determinations: the FERRE and the MATISSE estimations, based on a four dimensional (T eff , log g, [M/H], and [α/Fe]) grid of synthetic spectra (see de Laverny et al., 2012, for a description of the methodologies adopted for this grid); and the SME determination, based on a mean of the preliminary individual abundance estimations of Mg, Ca, Ti and Si. This global [α/Fe] is estimated in the first analysis phase, together with T eff , log g and [M/H] , and before the final GES recommended individual abundances are determined.
As explained in Recio-Blanco et al. (2014), the final global [α/Fe] reflects the information on the α-element abundances present in the analysed wavelength ranges for each type of star and spectrum quality. This approach, similar to the SDSS/SEGUE one described by Lee et al. (2011c), estimates the overall α-elements behaviour. The relation between the GES iDR1 global [α/Fe] and the individual α-element abundances, also included in the GES iDR1 (Mg, Ca, Ti and Si), is presented in Recio-Blanco et al. (2014). This analysis shows that the global [α/Fe] is dominated by the Mg abundance, with a standard deviation of the difference equal to 0.069 dex.
Finally, GES iDR1 measurements of the Fe abundance are also available, although they are not used in this paper.

Distance estimations and kinematics
The determination of the stellar distances is based on the method presented in Kordopatis et al. (2011) and successfully applied in Gazzano et al. (2013) and Kordopatis et al. (2013b). This method projects the atmospheric parameters T eff , logg, [M/H] and their errors on a set of isochrones, obtaining the most likely absolute magnitude of the star, given an a priori knowledge on the lifetime spent by a star on each region of the HR diagram (i.e. mainsequence stars are, statistically, more likely to be observed).
The isochrones on which the atmospheric parameters are projected are the ones of Yonsei-Yale (Demarque et al., 2004), combined with the colour tables of Lejeune et al. (1998). Using the provided interpolation code, we generated a set of isochrones with a constant step in age of 1 Gyr, starting from 1 Gyr to 14 Gyr. As far as the metallicities are concerned, the isochrones are within a range of −3.0 <[Fe/H]< 0.8 dex, constantly spaced by 0.1 dex. The α−enhancements of the isochrones have been selected in the following way, according to the typical [α/Fe] of disc and halo stars: On the other hand, Zwitter et al. (2010), in a similar approach to the one used here, shown that adopting isochrones of different alpha abundances does not affect the distance determination significantly.
In addition to the atmospheric parameters, the (J − K s ) colours are also used in order to obtain the most likely absolute magnitude of the star. Nevertheless, in order to do so, one must correct by the effect of the reddening, which is a function of the stellar distance and of the dust distribution. For that reason, we proceed in two steps. First, we compute the distances without using the colours, and we estimate iteratively the extinction at the distance of the star. Then, we apply the colour correction to the (J − K s ) of the star, and re-compute the distances one final time using the de-reddened colour.

Extinction correction
The band in which we compute the distance modulus is the J V HS band. For the brightest targets, no V HS (VISTA Hemisphere Survey) photometry is available, due to the saturation limits of the survey, and hence we adopt the 2MASS photometry, in which system we have J 2MAS S = J V HS and K 2MAS S = K V HS . In addition, the (JK) ES O provided by the Yonsei-Yale isochrones match very well the (JK s ) 2MAS S and (JK s ) V HS , so no colour transformation is needed when manipulating the magnitudes from the different photometric systems (Carpenter, 2001).
Once the absolute magnitude of the star is obtained using the Kordopatis et al. (2011) pipeline, a first estimation of the distance of the star is obtained, assuming the Schlegel et al. (1998b) reddening. For the lines-of-sight with the largest reddening (E(B − V) higher than 0.10 mag), we have applied the correction suggested by Bonifacio et al. (2000): where E(B − V) sch is the Schlegel et al. (1998b) reddening, and E(B − V) a is the adopted one. The extinction in the (JK s ) bands is then estimated using the following equations from McCall (2004): By applying the above extinction on the apparent magnitude and deriving the line-of-sight distance, we over-estimate the distance of the stars because the Schlegel et al. (1998b) reddening is a value integrated along the entire length of the line-of-sight. We hence compute the following correction factor to apply for the Schlegel reddening, assuming an exponential disc of dust: where h dust = 125 pc is the scale-height of the disc of dust causing the reddening (Misiriotis et al., 2006), D is the line-of-sight distance and b is the Galactic latitude of the target. Once the correction applied, we recompute the distance to the star, and following Ruchti et al. (2011), we repeat iteratively the above step until the difference between two consecutive distance estimations are smaller than 2%. For the few stars that lie within the dust lane, no more than four iterations are needed until convergence. The adopted reddening correction, as a function of the distance from the Galactic plane is shown in Fig. 3.
Once the final extinctions have been computed, the dereddened (J − K s ) colour of the star is estimated, and the pipeline in order to estimate the absolute magnitude of the star is ran, using this time the (J − K s ) colours. The effect on the derived lineof-sight distances is plotted in Fig. 4. The latter distances are the ones that are used in the remainder (see for instance Fig. 14).

Fig. 3.
Reddening correction computed for each star depending on its line-of-sight distance, as a function of its distance from the Galactic plane Z. The gray-scale colours represent the starcounts in log 10 (N). The contour lines contain 33%, 66%, 90% and 99% of the stars. As expected the corrections made to E(B − V) a are only for the stars located closer than 400 pc from the Galactic plane.

Errors in the line-of-sight distances
In order to estimate the errors in the derived distances, we performed 5 000 Monte-Carlo realisations on the errors on the distance modulus. Those errors include the error on the derived absolute magnitude J (which takes into account the errors in T eff , log g, [M/H]), and the error on the apparent J and K s magnitudes. The uncertainty in the atmospheric parameters is taken into account when computing the absolute magnitude, in the same fashion as described in Kordopatis et al. (2011), i.e. by estimating the dispersion of the weighted absolute magnitude of the stars on the isochrones. The parameter dominating the error in the distances is T eff for main sequence stars and log g for giants. Fig. 5 shows the cumulative histogram of the derived distance errors, for the entire sample and for a sub-sample with signal-to-noise ratios higher than 15 and errors in T eff and log g lower than 400 K and 0.5 dex, respectively. σ D is the inferred error bar when the above factors are taken into account. As shown in Fig. 5, 45% of the total sample has σ D /D < 20% and 80% has σ D /D < 50%. These values increase to 55% and 90% for the selected sub-sample of SNR>15 (2107 stars). In addition, Fig. 6 and Fig. 7 show the dependences of the distance errors on the metallicity and the distance value itself, for stars with SNR>15. The final errors in positions for the individual stars are included in Table 1.

Validation of the distances
We verified the accuracy of our results by evaluating the distances at which were found the globular clusters (GC) stars, observed and labelled as such by GES. The reference distances are the ones from Harris (1996), in the updated version (December Table 1. Derived stellar distances(D), positions (in the R and Z coordinates) and kinematics (radial, azimuthal and vertical velocities), plus their associated errors (∆D,∆R,∆Z,∆V R ,∆V p hi and ∆V Z ).  Table 2. Adopted globular cluster parameters (from Harris, 1996, in the catalogue updated version of December 2010). The two last columns show the derived distances with their associated errors (quadratic sum of the standard deviation of the members distances and their median individual errors).
Cluster  Fig. 4. Change in the line-of-sight distance estimates when using only the effective temperature, surface gravity and iron abundance (y−axis) and when using the previously mentioned atmospheric parameters in combination with the (J − K s ) colour (x−axis). The differences are mainly noticeable for the most distant stars. The colour-coding and the contour-lines are the same as in Fig. 3 2010) of the catalogue 1 , summarised in Table 2. In order to remove foreground contamination, in addition to using the GES labels, we have rejected all the targets that were outside the 3σ of the velocity dispersion of the GC (an arbitrary value of σ =10 km s −1 was adopted when this information was not available) and had iron abundances lower or greater than 0.5 dex of the mean cluster value (see Table 2). Figure 8 shows the histogram of the difference between the derived distance, D and the reference one, d Harris , normalised by the estimated individual error on the distances σ D (obtained from our pipeline). The error bars indicate the Poisson noise, i.e. the square root of the stars for each bin. Three stars were removed from the figure, as they were lying respectively at 1 http://www.physics.mcmaster.ca/ harris/mwgc.dat In the case where our analysis is bias free and the errors correctly estimated, the histogram of Fig. 8 is supposed to be Gaussian centred on zero and of unit dispersion. We can see that it is the case, validating in that way our approach. More precisely, once the three contaminators have been removed, we find a distribution centered at 0.20 and a standard deviation of 1.14. The mean value of the (D − d Harris )/σ D distribution corresponds to a bias in the estimated distances ((D − d Harris )/D) of less than 5%. The value of the sigma indicates an under-estimation of the errors of ∼ 14% with respect to the differences between the literature values and the ones derived in this work. However, we would like to stress that the Harris errors on the distance of the GCs could probably explain most of this difference. On the other hand, the slight asymmetry of the histogram comes from the  measurements of NGC 4 372, which has also the lowest metallicity and hence the atmospheric parameters with the largest uncertainties. In addition, it is worth noticing that the GES iDR1 globular cluster targets are giant stars, while the majority of the observed field stars are main sequence objects. Although GES is a magnitude limited sample, 80% of the targets are placed at smaller distances than the above analysed GCs. The red Gaussian is not a fit of the histogram but represents the unit Gaussian which should follow perfect bias-free measurements with correct error estimations. A very satisfactory result is obtained, with an estimated bias of ∼5% and a standard deviation of ∼ 14%.
Finally, it is worth mentioning that for dwarf stars, as explained in Recio-Blanco et al. (2014), a degeneracy between T eff and log g appears progressively, as T eff decreases, starting around T eff ∼5000 K. This degeneracy gives rise to a negative bias (stars may appear cooler and at lower gravities). We have explored the consequence on the derived distances on the above mentioned degeneracy, taking into account that the distances are calculated after the projection of the atmospheric parameters on the isochrones. In the corresponding cool part of the main sequence, T eff is the main parameter influencing the absolute magnitude determination, as the isochrones points of different ages and metallicities are all tightly located in the HR diagram. In addition, as the slope of the isochrones is very low, the determined distances are not very sensitive to errors in T eff . Moreover, the projection of the determined parameters on the isochrones, for the distance determination, corrects both T eff and log g. This correction is of the same order of the bias for stars with T eff >4500 K. For cooler stars (4500>T eff >4200 K), the residual bias in T eff , after the isochrone projection, is around 100 K. This error in T eff implies around 6−7% of distance error, depending on the metallicity. For this work, the selected stellar subsample of iDR1 (c.f. Section 4) contains only 3% of stars cooler than 4500 K in the thin disc sequence, and 2% in the thick disc one. Therefore, in the worst case, a bias of only ∼7% in distance, well inside the typical errors of the sample, can appear for about 2-3% of our final sample.

Derivation of Galactocentric velocities
From the derived line-of-sight distances, the three-dimensional Galactic positions are obtained for all the analysed stars (assuming R = 8 kpc and Z = 0 pc, see Reid, 1993). The derived values together with their associated errors are reported in Table 1. Once the three-dimensional Galactic positions are obtained, the galactocentric radial, azimuthal, and vertical velocities are computed using the PPMXL proper motions, and the equations from the appendix of Williams et al. (2013). We have used for (U , V , W ) = (11.1, 12.24, 7.25) km s −1 (Schönrich, Binney, & Dehnen, 2010) and the Local Standard of Rest (LSR) is at V LS R = 220 km s −1 . The obtained galactocentric velocities are reported in Table 1.

Error estimation in velocities
In order to evaluate the errors on the stellar velocities, we performed 5 000 Monte-Carlo realisations, considering that the proper motions, the radial velocities, and the line-of-sight distances are independent. Depending on the galactic coordinates (l,b) of a star, the uncertainties on the distances, the proper motions and the radial velocities will affect the 3D-velocity estimations in a different manner. Typically, the dominant source for the uncertainty comes from the errors in proper motions (typical errors of 8 mas/yr, see Sect.2), whereas the errors on the radial velocities (typically of 0.3 km/s) have in general a negligible effect. As an example, for a given star at 1 kpc in 290<l<340 and 10<b<30, the contribution to the 3D-velocity uncertainties will be of 5−10 km/s for an error of 10−20% on the distance and of 15−30 km/s for an error of 5−10 mas/yr on the proper motions. The final errors in kinematics for the individual stars are included in Table 1.

Chemical characterisation of the thin and the thick discs
Thanks to their high-resolution, the GES spectra offer an exceptional opportunity to characterise the distribution of disc stars in the [α/Fe] versus [M/H] plane, with a robust enough statistical analysis and rather low abundance errors. To this purpose, we have first selected sub-samples of stars with progressively higher errors in the abundance determination. This allows, on one hand, to find the balance between number statistics and abundance accuracy, and on the other hand, to test the effects that low numbers of stars and high measurement errors could have in the characterisation of the stellar populations. Fig. 9, shows the distribution of stars in the [α/Fe] versus [M/H] plane for four progressively larger sub-samples with increasingly higher error limits in the abundance determination. Table 3 shows the maximum errors and signal-to-noise limits considered for each of the analysed sub-samples. In addition, Fig. 10  The well-known decreasing evolution of [α/Fe] with metallicity is present in all the panels, within the expected values for disc stars. In addition, two sequences separated by a low density region can be identified: the thick disc one in the high-α regime and the thin disc one in the low-α regime (see below for the selection criteria). Both sequences overlap in metallicity from about -0.8 dex to -0.3 dex. For the last sub-sample of 1 952 objects, the increase in the scatter of the measurements blurs the chemical separation, although not completely.
First, as shown in Fig. 11, we have studied the position of the low density region along the [α/Fe] vs. [M/H] plane, thanks to the analysis of the [α/Fe] distribution in four intervals of metallicity. The largest sub-sample of stars in Fig. 11 (panel d of Fig. 9) was used for this purpose. The identified position of the gap in the [α/Fe] vs. [M/H] plane will be used in the following as the separation criterion between the thick and thin disc populations for further analysis. This low density region is found, after a visual examination of the histograms, at [α/Fe]=0.13, 0.19, 0.21 and 0.22 dex for the four central considered [M/H] intervals in Fig. 11.
Then, in order to more precisely characterise the thick and thin disc populations previously defined, and to take into account of the natural increase of the abundance measurement errors for metal-poor star's spectra , we have defined a new sub-sample of objects, hereafter called e (c.f. Table 3 The final sample is shown in Fig. 12 and contains 1 016 stars. A black line (based on the low density regions defined in Fig. 11) is plotted in Fig. 12 to identify the separation between the thick and thin disc sequences. Fig. 13 shows the Hertzsprung-Russel (HR) diagram of this selected sub-sample e. The red and blue points correspond to stars in the defined thin and thick disc sequences, respectively.  illustrates the values for 1 008 stars with errors smaller than 0.15 dex and 0.05 dex, respectively. Finally, panel d) shows all the stars with errors in T eff lower than 400 K, errors in log g lower than 0.5 dex and a spectral signal-to-noise ratio higher than 15 for the HR10 configuration (1 952 stars).  In addition, a high-α metal-poor population is also visible in our data, well above the Galactic plane (see Sect. 5), and occupying a flat sequence from [M/H]∼-1.0 dex down to the lowest sampled metallicities (filled triangles in Fig. 12). The mean [α/Fe] value of these halo stars is around 0.35 dex. This highα metal-poor sequence could be identified with the classical high-α halo sequence or, at least partially, with the metal-weak thick disc population (e.g. Ruchti et al., 2011;Kordopatis et al., 2013a). Moreover, in the same metal-poor metallicity regime ([M/H] values lower than ∼ -0.9 dex), a metal-poor low-α sequence seems to exist (asterisks in Fig. 12). After an individual inspection of these low-α metal-poor stars spectra, most of them have signal-to-noise values higher than 25 in the HR10 setup, and have thus reliable [M/H] and [α/Fe] measurements. When only the lower error measurements are considered (stars with errors smaller than 0.05 dex in [α/Fe] and 0.15 dex in [M/H] and signal-to-noise ratio for the HR10 spectra higher than 20), the low-α metal-poor stars represent about 10% of the observed stars with [M/H] values lower than ∼ -0.9 dex. These stars could confirm the low-α metal-poor halo stars first observed by Nissen & Schuster (2010), and also reported by Nissen & Schuster (2011), Schuster et al. (2012, Adibekyan et al. (2012) and Ishigaki et al. (2013). On the other hand, the possible link between those objects and the metal-weak thick disc population (e.g. Morrison et al., 1990;Carollo et al., 2010;Kordopatis et al., 2013b) seems to be ruled out by their low [α/Fe] ratios, as metal-weak thick disc stars are reported to have [α/Fe]∼ 0.3 dex (Ruchti et al., 2011). As a complementary description of these stars, the individual element abundances of these stars are analysed in Mikolaitis et al. (2014, in preparation).

Comparison with literature studies
First of all, we have compared the spatial distribution of the analysed targets with that of previous spectroscopic studies of large numbers of disc stars, with available measurements of the [α/Fe] ratio, as it is the case of the Lee et al. (2011b) and Adibekyan et al. (2012) analysis. We recall that the former is a low-resolution study of stars with Galactic radius in the range 7−10 kpc, and the latter is a high-resolution analysis of solarneighbourhood stars, located at a maximum distance of about   600 pc from the Sun. Fig. 14 shows the location of the GES targets analysed in this paper (e sub-sample of Fig. 12), in vertical distance to the Galactic plane and radial cylindrical Galactic coordinate R. To understand how many of the targets are outside the solar neighbourhood, the points are colour coded by distance. All the objects with distances smaller than D=600 pc from the Sun are plotted in red and those with larger distances in black. As expected from the photometric target selection, 90% of the stars are at distances larger than D=600 pc from the Sun (i.e. much further than the Adibekyan et al. (2012) sample, comprised mainly of stars within 50 pc).
Regarding the thin to thick disc separation in the [α/Fe] vs. [M/H] plane, the above identified gap is in agreement with the separation seen in the Adibekyan et al. (2012) data. Fig. 12 shows the Adibekyan et al. (2012) division, as a red line with three segments, that can be compared to our proposed separation, in black. Both lines are compatible, with very small differences probably due to possible offsets in the abundance measurements between the GES analysis and the Adibekyan et al. ones could exist. This could come mainly from the identification of the metal-rich part of thick disc, that would be more difficult to separate from the thin disc from low resolution data, due to the small differences in the [α/Fe] values. Furthermore, the difference of slope observed in our sample between the thin and the thick disc sequences seems in agreement with the Haywood et al. (2013) analysis of Adibekyan et al. (2012) data, that already points out an iron enrichment of the thin disc lower than the thick disc one.
Finally, concerning the metal-rich end of the thin disc, no αrich metal-rich stars are observed in our sample, nor in the Lee et al. (2011b) one, contrary to the results of Adibekyan et al. (2011) and Gazzano et al. (2013) and Boeche et al. (2013), restricted to, or predominantly in, smaller distances from the Sun. Furthermore, we note here that the number of stars with [M/H] higher than -0.1 dex is rather small in our sample, compared to the previously mentioned nearby studies.

Distribution of distances to the Galactic plane
As already mentioned in the Introduction, the thin and the thick discs are characterised by different scale heights. In order to explore the transition between these two populations, keeping in mind their chemical separation in the [α/Fe] vs. [M/H] plane as described in the above section, we have analysed the distribution of the derived distances to the Galactic plane (from the absolute values of the Galactic Z coordinate) for the target stars.
The left panel of Fig. 15 shows a three-dimensional view of the distance to the Galactic plane, in the Z-axis, as a function of [M/H] and [α/Fe], for the analysed stars (sub-sample e of Fig. 12). This illustration is complemented, in the right panel, with a grey scale-coded two-dimensional plot of the distribution of Z values in the [α/Fe] vs. [M/H] plane. The red line shows the separation between the thick and the thin discs sequences, defined in Sec. 4 (c.f. Fig. 12). We point out that the possible halo stars discussed in the previous section are not shown in Fig. 15, because of the low number of objects per bin.
Regarding the thick disc population, the mean distances to the Galactic plane seem to decrease along the sequence, from the metal-poor end to the metal-rich one: Therefore, as the metallicity increases, the thick disc seems to lay in thinner and thinner layers. On the other hand, the distribution of vertical distances to the plane with respect to the Galactocentric radial distance is not perfectly homogeneous (c.f. Fig. 14). More particularly, our sample is dominated by stars further away from the galactic plane in the inner parts of the disc with respect to the typical heights probed in regions outside of the solar radius. Thus, different Z heights can be representative of different ranges in R, and, if there were significant metallicity variations with respect to R, this would manifest as a variation in metallicity with height. However, as later explained in Sec. On the other hand, the thin disc sequence shows a quite constant value of the mean distances of its stars to the Galactic plane of Z∼0.5 kpc, from the metal-poor side ([M/H] ∼-0.8 dex) to the metal-rich one ([M/H] ∼0.3 dex). Again, this conclusion should not be influenced by the spatial coverage of the studied sample, despite the fact that the thin disc has a metallicity gradient with respect to R (c.f. Sec. 9.2). This is because the thin disc Z range is fairly homogeneously covered at the full range of Galactocentric radius, and, in particular, the mean Z value of the stars in the thin disc sequence does not present any clear trend with R.
As a consequence of the above conclusions, the transition from the thick disc sequence to the thin disc one, in the [α/Fe] vs. [M/H]plane, implies a steeper change in the mean Z distances at [M/H] around -0.8 dex than around -0.3 dex. This coincides with the fact that the chemical separation between both sequences (by a gap or a lower density region) seems clearer in the metal poor regime than in the metal rich one.
Finally, we point out that, as expected, there is a correlation between the Z distances and the associated errors, with larger errors for stars that are further from the plane. More particularly, the stars in the metal-poor high-alpha end ([M/H] around -1.0 dex) have errors around 0.45 kpc, while at [M/H] around -0.7 dex the error is about 0.20 kpc for the thick disc stars and around 0.075 kpc for the thin disc. Nevertheless, while the stars that are further from the plane have larger absolute distance errors, the percent error is actually smaller, as indicated in Figures 6 and 7.

Distribution of rotational velocities
The distribution of the cylindrical rotational velocities, V Φ , in the [α/Fe] vs. [M/H] plane, is shown in Fig. 16. The rotational velocities of the stars located in the thick disc sequence have a mean value of V Φ =176±16 km/s, and those in the thin disc have a mean V Φ =208±6 km/s. In addition, for the thick disc, the V Φ values seem to progressively increase from the metalpoor α-rich end (with V Φ ∼150 km/s) to the metal-rich α-poor one (V Φ ∼190 km/s), reaching values close to the thin disc ones.
On the other hand, a correlation between rotational velocity and metallicity is observed for both the thin and the thick disc stars, although with opposite signs. In order to exploit the fact that both [α/Fe] and [M/H] are available for our sample, we have decided to analyse the correlation of rotational veloc-  First, all the bins along the thin and thick disc sequences were considered. For the thin disc stars, we find a gradient of Galactic rotation with metallicity of -17±6 km s −1 dex −1 . We have checked the possible influence of the large range of sounded galactocentric radius on the derived gradient, by splitting the sample in two: stars inwards (R< 8kpc) and outwards (R> 8kpc) of the solar radius. The resulting values of the gradient are -20±13 and -22±9 km s −1 dex −1 , suggesting that the derived slopes are not indirectly created by radial gradients in metallicity. For the thick disc population, we find a positive gradient of 43±13 km s −1 dex −1 . The derived values of the gradient for the thin and the thick disc are in very good agreement with the Lee et al. (2011b) analysis of SEGUE data (who find -22.6±1.6 km s −1 dex −1 for the thin disc and 45.8±1.8 km s −1 dex −1 for the thick disc) and with the Adibekyan et al. (2013) measurements in the solar neighbourhood (-16.8±3.7 km s −1 dex −1 for the thin disc and 41.9±18.1 km s −1 dex −1 for the thick disc). Our gradient for the thick disc is also in agreement with the measurements of Kordopatis et al. (2011, 45±12 km s −1 dex −1 ) and Spagna et al. (2010), derived, respectively, from low resolution FLAMES spectroscopy and SDSS spectro-photometry, with no estimations of [α/Fe]. It also confirms the disagreement with the SDSS result of Ivezić et al. (2008) Fig. 17. The gradient of rotation with metallicity seems to be steeper than the mean (64±9 km s −1 dex −1 instead of 43±13 km s −1 dex −1 ) when the high-α bins are considered. On the contrary, the value of ∆V Φ /∆[M/H] is lower than the mean one when it is determined from the low-α bins (30±10 km s −1 dex −1 ). Given that the [α/Fe] dispersion is rather comparable to the measurement errors, the existence of a real intrinsic dispersion can not be concluded. Indeed, the described dependence of the ∆V Φ /∆[M/H] gradient with the considered [α/Fe] bins could be the sign of a contamination of the thick disc stars from the thin disc ones in the lower part of the thick disc sequence. If this is correct, the steeper gradients derived from the high-α bins could be more characteristics of the thick disc population than the previously published ones of about 43±13 km s −1 dex −1 . To test this possibility, we have also derived the value of the slope, using all the bins along the thick disc sequence with a mean error in [α/Fe] lower than 0.04 dex. The resulting slope is 56±19 km s −1 dex −1 instead of 43±13 km s −1 dex −1 , confirming that the slope seems steeper for higher quality data.
Finally, to analyse the direct influence of [α/Fe] in the stellar rotation for thick disc stars, we have also derived the correlation between rotation and [α/Fe] separate from [M/H]. As shown in Fig. 18, V Φ and [α/Fe] are highly anticorrelated (the Spearman's coefficient value is -0.85, corresponding to 0.72% of the V Φ variance). In addition, the correlation rank between V Φ and [α/Fe] is even more significant than the one between V Φ and [M/H] (having a Spearman's coefficient value of 0.73, that accounts for 0.53% of the V Φ variance). Furthermore, we have checked how azimuthal velocity varies with vertical height for the thick disc sequence stars. Fig. 19 shows the mean V Φ value as a function of Z, in four bins of 0.5 kpc. The error bars indicate the standard error of the mean, in each bin. Clearly, the azimuthal velocity seems to decrease with height, in particular for Z values higher than about 1 kpc. As a consequence, this could explain why in situ samples of the Fig. 19. Mean V Φ value as a function of vertical height for the thick disc sequence stars. thick disc can show a higher value of the thick disc lag (hence lower V Φ values) than the local surveys (see for instance Wyse et al., 2006).

Distribution of azimuthal and vertical velocity dispersions
The stellar velocity dispersions can also be studied for the two chemically separated thin and thick disc sequences with our data. To this purpose, as low quality velocity determinations increase artificially the velocity dispersions, we have further restricted our analysis sub-sample, allowing only stars with errors in azimuthal and vertical velocity lower than 70 km/s. This value was found imposing a good balance between the statistical robustness and the allowed maximum error in velocity. This new restriction, in addition to those of sub-sample e presented in Table  3, leaves a data set of 645 stars, with a mean error in V Φ of 31 km/s. Fig. 20 shows the distribution of the measured (not error free) azimuthal and vertical velocity dispersions (σ Φ and σ Z ) as a function of metallicity and distance to the Galactic plane. The error bars come from the standard error of the derived standard deviations. As expected, both velocity dispersions are higher for the thick disc sequence than for the thin disc one. The mean value of σ Φ is 54±5 km/s for the thick disc and 40±3 km/s for the thin disc. Regarding σ Z , the mean values are 48±5 km/s for the thick disc sequence and 32±3 km/s for the thin disc one. As shown in the upper panels of Fig. 20, no significant dependencies of either σ Φ or σ Z with metallicity or Z seem to be present in the data. Nevertheless, for the thick disc, the measurements suggest a possible increase of σ Φ with distance to the Galactic plane that is not visible for σ Z . This conclusion seems strengthened by the increase of the uncertainty on the derived mean V Φ values as Z increases, as shown in Fig. 19. A possible decrease of σ Φ with increasing metallicity could also be present for the thick disc. In addition, around [M/H]∼-0.3 dex, the dispersion in the azimuthal velocity of the thick and the thin discs seem roughly the same, with a further decrease of σ Φ for the metal rich end of the thin disc. This seems to be in agreement with the behaviour of the azimuthal velocity as a function of the metallicity, shown in Fig. 17 (see also Section 6). Moreover, the metal-poor side of the thin disc sequence, that presents higher rotation values, shows also a slightly lower azimuthal velocity dispersion.
On the other hand, we point out that our determinations of the velocity dispersions could be affected by the lack of stars at very low Z values (only 37% of the thin disc stars and 17% of the thick disc ones, included the first considered bin of vertical distance, are at abs(z)<250 pc). This could explain the slightly higher values of dispersion found close to the Galactic plane with respect to other surveys like RAVE (e.g. Burnett, 2010).
Finally, for the thick disc sequence, the azimuthal velocity dispersion seems to be of the order of the vertical velocity dispersion, within the error bars. In any case, the analysed data do not suggest, as predicted by Binney (2012), a smaller σ Φ than the corresponding σ Z for the thick disc.

Orbital parameters
To calculate the orbital parameters of the stars in the sample we integrated orbits in a Galactic potential consisting of a Miyamoto & Nagai (1975) disc, a Hernquist (1990 bulge and a spherical logarithmic halo. Their characteristic parameters took the values M disk = 9.3 × 10 10 M and a = 6.5 kpc, b = 0.26 kpc for the disk, M bulge = 3 × 10 10 M and c = 0.7 kpc for the bulge, and v h = 134.5 km/s, and d = 12 kpc for the halo. In particular, the eccentricities of the stars were defined as (r apo − r peri )/(r apo + r peri ) where r apo (r peri ) is the maximum (minimum) distance reached by the star in its orbit in the past 2 Gyr. It is today widely accepted that the shape of the eccentricity distribution of the stellar orbits can provide important constraints on the formation mechanisms of the Galactic disc (Sales et al., 2009). Figure 21 shows the distribution of orbital eccentricities for the thin (red line) and the thick disc (blue line) chemically Fig. 20. σ Φ (left panels) and σ Z (right panels) as a function of metallicity (upper panels) and distance to the Galactic plane (lower panels). Filled circles correspond to the chemically selected thick disc sequence and open circles to the thin disc sequence. selected stars. As expected, the thin disc distribution is strongly peaked at a low eccentricity value of ∼ 0.15. On the other hand, the thick disc is significantly broader and its median eccentricity is around 0.3. Compared to previous determinations, such as those of Wilson et al. (2011), Lee et al. (2011b) and Kordopatis et al. (2011), the shape is flatter and does not decline as fast at high eccentricities, although because of the small number of stars, the significance of this may be questionable.

Eccentricity distribution
In order to evaluate the errors in the computed eccentricities, we have calculated the median, 25% and 75% percentile values of the eccentricity from all orbits starting from 100 realisations of the stars position and velocity, convolved with the errors. Assuming a Gaussian distribution, we have derived a value of the standard deviation for each stellar eccentricity. The mean value of the standard deviation depends slightly on the eccentricity itself, and its equal to 0.14 at eccentricity=0.20 and 0.13 at eccentricity=0.80. Furthermore, Figure 22 presents the stellar orbital eccentricities as a function of metallicity for the thin and the thick disc sequences (lower and upper panels, respectively). On one hand, the thin disc sequence presents a very similar eccentricity distribution through its whole metallicity range, with a majority of low eccentric orbits. In addition, when only thin disc stars with a distance to the plane lower than 400 pc are considered (blue open circles in Fig. 22), most of the thin disc stars with higher eccentric orbits disappear. In addition, this high eccentricity tail of the thin disc suffers from a very low statistics and the relatively large errors in the derived eccentricities. On the other hand, the thick disc sequence seems to present highly eccentric orbits mainly in the metal-poor regime ([M/H]< −0.6 dex). Additionally, the metal-rich part of the thick disc sequence shows very similar eccentricities to the thin disc one. A similar behaviour of the mean eccentricity distribution with metallicity has also been observed by Adibekyan et al. (2013) in the solar neighbourhood. Finally, the derived slopes of the eccentricity dependence on metallicity, ∆e/∆[M/H], are -0.158±0.081 dex −1 for the thick disc and -0.027±0.045 dex −1 for the thin disc. These values are also in agreement with those determined by Adibekyan et al. (2013) for the solar neighbourhood (-0.184±0.078 dex −1 and -0.023±0.015 dex −1 for the thick and the thin disc respectively).

Maximum height above the plane
We have also studied the distribution of maximum heights above the Galactic plane, Zmax, for the stars in the thick disc sequence. To this purpose, as shown in Fig. 23, we have plotted the normalized Zmax distribution for three different metallicity intervals. The associated error to Zmax was calculated through the above described procedure, already used for the eccentricity error, and it is of the order of 0.6 kpc.
As expected from the analysis of the mean distances to the Galactic plane in Sect. 5, the Zmax distribution seems to peak at lower values as the metallicity increases. Since we are working with a colour/magnitude selected sample, the quantification of a gradient would need correction for the volume sampled as a function of metallicity. Nevertheless, the observed Z and Zmax distributions seem to suggest a rather steep vertical metallicity gradient in the thick disk. On the other hand, Katz et al. (2011) found a low gradient of -0.068±0.009 dex kpc −1 , derived without a chemical selection of the thick disc in the [α/Fe] vs.
[M/H] plane. On the other hand, Bilir et al. (2012) derived a value of the vertical metallicity gradient in the thick disc of -0.034±0.003 dex kpc −1 , from a kinematically selected sample of RAVE stars. Finally, Ruchti et al. (2011) found an iron abundance vertical gradient of -0.09±0.05 dex kpc −1 for metal poor thick disc stars.

Radial gradients in rotation, metallicity and [α/Fe]
9.1. The thin disc rotational velocity behaviour with Galactocentric radius.
First of all, as illustrated in Fig. 24, we examined a possible dependence of V Φ with the Galactocentric radius for the thin disc. To this purpose, we have selected from the thin disc stars of sub-sample e those within a distance to the Galactic plane lower than 700 pc. We have also selected stars having errors in V Φ smaller than 40 km/s for 7<R<9 kpc, up to 60 km/s for 6<R<7 kpc and 9<R<10 kpc and 50% elsewhere. In addition, we have performed a 2-σ clipping around the mean values of V Φ in three metallicity intervals ( Demers & Battinelli (2007), from gas clouds using the terminal velocities of the HI and CO lines (Sofue et al., 2009) and from APOGEE data (Bovy et al., 2012a), find an approximately flat behaviour, as what is found for the rotational velocity with our data.
In addition, we have tested the possible effect of the V Φ correlation with metallicity (c.f. Sect. 6) on the derived V Φ vs. R relation. To this purpose, we have divided the above described sample in three sub-samples for different metallicity intervals:  Fig 24). Then, three values of the slope were derived for each of the sub-samples. For all these metallicity intervals, the rotational velocity stays flat within the error bar. Therefore, no measurable influence of the stellar rotation correlation with metallicity seems to exist in the thin disc, within the studied galactocentric distance range. In agreement with Fig. 17, the mean value of V Φ decreases with metallicity, in the three considered metallicity intervals. Finally, we have checked a possible trend of the V Φ dispersion with R in the examined sample, finding a flat σ Φ with R in the range 7<R<10 kpc, for which the number of stars allows a robust enough measurement of the V Φ dispersion.

Radial abundance gradients
In order to derive the radial abundance gradients in the thin disc, we have used the sample of stars described in Section 7.1 and performed least square fits to the [M/H] v.s. R distribution and the [α/Fe] v.s. R distribution (c.f. Fig. 25). We find a negative gradient of -0.058±0.008 dex kpc −1 for [M/H] and a very small positive gradient of 0.012±0.002 dex kpc −1 for the [α/Fe]. As in the previous section, the uncertainties come from the standard error of the slope, derived by the least square linear regression.
On the other hand, the very small positive gradient in the [α/Fe] abundance ratio seems not to be influenced by the cut in vertical distance to the Galactic plane, as its value does not change if the cut is done at |z| <400 pc instead of 700 pc. This could be the case if a contamination by thick disc stars was affecting the radial gradient measurement. However, the presence of numerous metal-poor thin disc stars with [M/H]< −0.5 dex in our sample, influences the measurement, steepening the gradient, as those stars present high values of [α/Fe] and dominate the outer regions. If real, this slightly positive gradient in [α/Fe] could be in agreement with the measurements of Yong et al. (2005)  Moreover, we have checked the influence of the data spatial distribution in the derived gradients. For that purpose, it is necessary to test the robustness of the two derived gradients on the inhomogeneity of the survey in the R and Z distribution. As al-ready explained in Section 5, our total sample is dominated by stars further away from the galactic plane in the inner parts of the disc with respect to the typical heights probed in regions outside of the solar radius. However, thanks to the restriction of the thin disc sample to Z values smaller than 700 pc, the above mentioned correlation between Z and R practically dissapears: the mean distance to the Galactic plane is ∼450 pc at R∼6 kpc and ∼300 pc at R∼12 kpc. This small trend can not be responsible for the measured gradients in metallicity and [α/Fe], especially, after the chemical selection of the stars from the thin disc sequence of the [M/H] vs. [α/Fe] plane. In addition, in the case of the [α/Fe] gradient, the small negative trend of Z with R excludes a higher spurious thick disc contamination in the outer regions with respect to the inner ones, as we would expect if the [α/Fe] gradient was artificially created by a bad thin-thick disc separation in the outer regions. Therefore, as a consequence of the two gradients in [M/H] and [α/Fe], the outer parts of the thin disc sounded by the present sample of GES data (up to R ∼12 kpc) appear more metal-poor and, possibly, slightly richer in [α/Fe] than the inner regions. Finally, we have also analysed the radial abundance gradients in the thick disc population (c.f. Fig 26)

Summary and discussion
We have used the atmospheric parameters, [α/Fe] abundance ratios and radial velocities, determined from the Gaia-ESO Survey GIRAFFE spectra of FGK-type stars (first nine months of observations), to present the kinematical and chemical characterisation of the observed stellar populations. First, we have chosen to identify the thin to thick disc separation in the [α/Fe] v.s. [M/H] plane, thanks to the presence of a low-density region in the number density distribution. This separation is in agreement with the recently proposed one by Adibekyan et al. (2013) from a sample of solar neighbourhood stars, and the ones observed before by Fuhrmann (2004), Reddy et al. (2006) and Bensby et al. (2007), among others.
Furthermore, we have studied the distribution of the distances to the Galactic plane, in the [α/Fe] v.s. [M/H] plane, with particular attention to the chemically defined thin and thick disc populations. Regarding the thick disc, our results show that the stars lay in progressively thinner layers along the sequence (as metallicity increases and [α/Fe] decreases). This finding could be in agreement with the Haywood et al. (2013) proposal of a thick disc forming stars in thinner and thinner layers for 4-5 Gyrs, also predicted by the dissipational settling formation scenario of Burkert et al. (1992). Similarly, Bovy et al. (2012c) already pointed out the scenario of an evolution in gradually thinner layers (as suggested by the dependence of their scale height distribution on [α/Fe]), although no clear separation between the thin and the thick disc sequences appears in their data. In addition, the recent hydrodynamic simulations of Bird et al. (2013) and Stinson et al. (2013) also show discs in which the younger stars have progressively shorter vertical scale heights. On the other hand, the thin disc sequence presents a constant value of the mean distance to the Galactic plane at all metallicities. As a consequence, the thick to thin disc transition in Z values appears more abrupt in the low metallicity regime.
Concerning the kinematics, the mean rotational velocities found for the thick and the thin disc sequences are in agreement with the canonical ones (V Φ =176±16 km/s and V Φ =208±6 km/s). In addition, the mean Galactic rotation also seems to increase progressively along the thick disc sequence, as the metallicity increases. Our data also confirm the already observed correlations between V Φ and [M/H] for the two discs. In the case of the thick disc sequence, the mean value of the ∆V Φ /∆[M/H] is 43±13 km s −1 dex −1 , is in agreement with the previous findings of Lee et al. (2011b), Kordopatis et al. (2011) and Spagna et al. (2010). This would confirm the disagreement with the SDSS result of Ivezić et al. (2008), based on photometric metallicities. Moreover, we have explored the possible influence of contaminating stars from the thin disc in the analysis of the correlation. If only the stars with the higher [α/Fe] values at each metallicity bin are considered, the derived ∆V Φ /∆[M/H] is 64±9 km s −1 dex −1 . As a consequence, this steeper value of the gradient, cleaned from the influence of possible thin disc contaminants could be more characteristic of the real correlation for the thick disc. For the thin disc, a negative gradient of ∆V Φ /∆[M/H]=-17±6 km s −1 dex −1 is found. Consequently, the difference in the characteristic rotational velocities between the thin and the thick disc are again accentuated in the low metallicity regime, as for the mean distances to the Galactic plane. Furthermore, the analysis of the velocity dispersions for the chemically separated thin and thick disc stars confirms a higher dispersion for the thick disc population than for the thin disc one, as expected, for both V Φ and V z . Further, the stars in the thick and the thin disc sequences present also different distributions of their stellar orbit eccentricities, with the thick disc stars presenting a broader distribution peaked at higher eccentricities than the thin disc one.
Finally, we have analysed the properties of the thin and thick discs as a function of the Galactocentric radius. For the thin disc, we estimate a flat behaviour of the thin disc rotational velocity, a metallicity gradient equal to -0.058±0.008 dex kpc −1 and, a very small positive [α/Fe] gradient. For the thick disc, no gradients in [M/H] and [α/Fe] are found.
The general picture emerging from our analysis is that of a thick disc that lays in progressively thinner and thinner layers as the metallicity increases and the [α/Fe] decreases with time (until about [M/H] = −0.25 dex and [α/Fe] = 0.1 dex). During this settling process, the thick disc rotation increases progressively and, possibly, the azimuthal velocity dispersion decreases. No measurable metallicity or [α/Fe] gradients with galactocentric distance seem to remain for the thick disc. At about [M/H] = −0.25 dex, the mean characteristics of the thick disc in [α/Fe], vertical distance to the plane, rotation and rotational dispersion are in agreement with (or just slightly different than) that of the thin disc stars of the same metallicity, suggesting a possible connection between the two populations at a certain moment of the disc evolution. This is in agreement with the Haywood et al. (2013) scenario, in which the inner thin disk (R<10kpc) inherited from the chemical conditions left at the end of the thick disk phase (the observed metal-rich end of the thick disc sequence in the [α/Fe] vs. [M/H] plane), some 8 Gyr ago. In their view, there may be an age gap in the star formation possibly causing the gap in [α/Fe] enrichment, after which, star formation proceeds in a thin disk.
For the thin disc, no clear gradient of the mean vertical distance to the plane with the stellar metallicity is observed. On the contrary, the stars in the metal-poor thin disc sequence tail seem to have a mean Galactic rotation velocity higher than those in the metal-rich part. Moreover, if confirmed, the existence of a small positive gradient in [α/Fe] with R leads to a thin disc that would be slightly more metal-poor and [α/Fe]-rich in the outer parts. These two characteristics, similar to those of the metalpoor end of the thin disc sequence in the [α/Fe] v.s. [M/H] plane, would favour the picture, suggested by Haywood et al. (2013) and Haywood (2006), of this metal-poor thin disc being formed in the outer parts.
Contrary to the similarities between the thin and the thick discs around [M/H]∼ −0.3dex, the stars in the metal-poor overlapping regime are more clearly separated in their mean physical properties (chemistry, kinematics, possibly orbital eccentricities). This separation seems to be confirmed by the gap (or lowdensity region) in the [α/Fe] vs. [M/H].
On the other hand, the tight bimodal distribution of [α/Fe] at every value of [M/H] poses severe constraints to models predictions. Recent attempts to define an age-metallicity relation for stars currently in the solar neighbourhood (e.g. Haywood et al., 2013;Mishenina et al., 2013) hint at an evolutionary path in which a very wide range of ages corresponds to a single [M/H] value. Haywood et al. (2013) results (c.f. their Fig. 8) show that at a value of about [M/H]=-0.5 dex stars have an age range of around 10 Gyr. On the other hand, the present work and other studies focussed in the solar neighbourhood (e.g. Adibekyan et al., 2013) show that the distribution of [α/Fe] at [M/H]=-0.5 dex appears narrowly bimodal, with two essentially unresolved but quite discrete values of [α/Fe] being found. As a consequence, although stars with [M/H]≈-0.5 dex formed over 10Gyr, they seem to be formed with only two values of [α/Fe]Ṫhis is a quite remarkably severe constraint that can be analysed in the light of Galactic models in which radial migration (e.g. Sellwood & Binney, 2002;Schönrich & Binney, 2009c,b;Roškar et al., 2008), and churning, is a major factor. Schönrich & Binney (2009b) claim that a bimodal, although continuous, distribution of [α/Fe] can naturally come out as a consequence of the standard assumptions about star formation rates and metal enrichment, with no need to invoke breaks in the Galaxy's star formation history or accretion events. However, the tightness of the well populated observed thick and thin disc sequences in the [α/Fe] vs. [M/H] plane, with possibly only two values of [α/Fe] for a given metallicity, poses severe constraints.
If radial migration occurs, one can suppose an evolutionary path in which enrichment in the high-α branch happened rapidly, up to metallicities close to [M/H]= −0.25 dex. This should occur across a large enough radial range of the proto-galactic disk. In fact, to preserve the tightness of the high-α sequence, the stars should had formed in a narrow time interval from radially wellmixed interstellar medium (ISM) gas, through all the galactocentric radii which generated stars that are today at a given R (ex. in the solar neighbourhod). Then, to explain the low-α sequence, one can postulate a significant gas accretion, reducing the disc ISM metallicity to, for instance, -0.6 dex. This should be followed by disc evolution in which large scale gas flows (as in fountains, c.f. Fraternali et al., 2013, and references therein) retain a tight [α/Fe]-[M/H] correlation, over the radial region which today populates the solar neighbourhood. In addition, this scenario should, somehow, also conciliate the higher rotation values observed in the metal-poor thin disc regime. Othewise, other possibility to explain the observed [α/Fe] vs. [M/H] distributions, is to suppose that strong radial migration in the sense of churning did not occure (see also the discussion in Haywood et al., 2013).
Finally, this study also shows the importance of precise chemical abundance measurements to disentangle the stellar population puzzle of the Galactic disc. The stellar chemical patterns can guide the definition of useful stellar sub-samples to unveil the evolutionary paths in the disc formation history and to give clearer constraints to the models. Moreover, in the near future, the Gaia mission of the European Space Agency will allow precise distances and ages estimations for all the analysed stars. The results presented here, based only on the first months of the GES observations, confirm how crucial are today large high-resolution surveys outside the solar neighbourhood for our understanding of the Milky Way history.