Reconstructing the Galactic evolution of lithium (Li) is the main tool used to constrain the source(s) of Li enrichment in the Galaxy. Recent results have suggested a decline in Li at supersolar metallicities, which may indicate reduced production. We exploit the unique characteristics of the Gaia-ESO Survey open star cluster sample to further investigate this issue and to better constrain the evolution of Li at high metallicity. We trace the the upper envelope of Li abundance versus metallicity evolution using 18 clusters and considering members that should not have suffered any Li depletion. At variance with previous claims, we do not find any evidence of a Li decrease at high metallicity. The most metal-rich clusters in the sample ([Fe/H] about 0.3) actually show the highest Li abundances, with A(Li) > 3.4. Our results clearly show that previous findings, which were based on field stars, were affected by selection effects. The metal-rich population in the solar neighbourhood is composed of relatively old and cool stars that have already undergone some Li depletion; hence, their measured Li does not represent the initial interstellar medium abundance, but a lower limit to it.

The Gaia-ESO Survey (GES) observed many open clusters as part of its programme to spectroscopically characterise the various Milky Way populations. GES spectroscopy and Gaia astrometry from its second data release are used here to assign membership probabilities to targets towards 32 open clusters with ages from 1-3800 Myr, based on maximum likelihood modelling of the 3D kinematics of the cluster and field populations. From a parent catalogue of 14398 individual targets, 5033 stars with uniformly determined 3D velocities, $T_{\rm eff}$, $\log g$ and chemistry are assigned cluster membership with probability $>0.9$, and with an average probability of 0.991. The robustness of the membership probabilities is demonstrated using independent membership criteria (lithium and parallax) in two of the youngest clusters. The addition of radial velocities improves membership discrimination over proper motion selection alone, especially in more distant clusters. The kinematically-selected nature of the membership lists, independent of photometry and chemistry, makes the catalogue a valuable resource for testing stellar evolutionary models and investigating the time evolution of various parameters.

In the era of large spectroscopic surveys, massive databases of high-quality spectra provide tools to outline a new picture of our Galaxy. In this framework, an important piece of information is provided by our ability to infer stellar ages. We aim to provide empirical relations between stellar ages and abundance ratios for a sample of solar-like stars. We investigate the dependence on metallicity, and we apply our relations to Gaia-ESO samples of open clusters and field stars. We analyse high-resolution and high-S/N HARPS spectra of a sample of solar-like stars to obtain precise determinations of their atmospheric parameters and abundances through differential spectral analysis and age through isochrone fitting. We investigate the relations between ages and abundance ratios. For the abundance ratios with a steeper dependence on age, we perform multivariate linear regressions, including the dependence on metallicity. We apply our relations to a sample of open clusters located in 4<R$_{GC}$<16 kpc. Using them, we are able to recover the literature ages only for clusters located at R$_{GC}$>7 kpc. In these clusters, the content of s-elements is lower than expected from chemical evolution models, and consequently the [s/$\alpha$] are lower than in clusters of the same age located in the solar neighbourhood. With our chemical evolution model and a set of empirical yields, we suggest that a strong dependence on the star formation history and metallicity-dependent yields of s-elements can substantially modify the slope of the [s/$\alpha$]-[Fe/H]-age relation in different regions of the Galaxy. Our results point towards a non-universal relation [s/$\alpha$]-[Fe/H]-age, indicating the existence of relations at different R$_{GC}$ or for different star formation history. A better understanding of the s-process at high metallicity is necessary to fully understand the origin of these variations.

Recent and on-going large ground-based multi-object spectroscopic surveys allow to significantly increase the sample of spectroscopic binaries to get insight into their statistical properties. We investigate the repeated spectral observations of the Gaia-ESO Survey (GES) internal data release 5 to identify and characterize spectroscopic binaries with one visible component (SB1) in fields covering the discs, the bulge, the CoRot fields, and stellar clusters and associations. A statistical chi2-test is performed on spectra of the iDR5 sub-sample of approximately 43500 stars characterized by at least 2 observations and a S/N > 3. Our sample of RV variables is cleaned from contamination by pulsation/convection-induced variables using Gaia DR2 parallaxes and photometry. Monte-Carlo simulations using the SB9 catalogue of spectroscopic orbits allow to estimate our detection efficiency and to correct the SB1 rate to evaluate the GES SB1 binary fraction and its dependence with effective temperature and metallicity. We find 641 (resp., 803) FGK SB1 candidates at the 5 sigma (resp., 3 sigma) level. The orbital-period distribution is estimated from the RV standard-deviation distribution and reveals that the detected SB1 probe binaries with log(P[d]) < 4. We estimate the global GES SB1 fraction to be in the range 7-14% with a typical uncertainty of 4%. The GES SB1 frequency decreases with metallicity at a rate of -9+/-3%/dex in the metallicity range -2.7<FeH<+0.6. This anticorrelation is obtained with a confidence level higher than 93% on a homogeneous sample covering spectral types FGK and a large range of metallicities. When the present-day mass function is accounted for, this rate turns to 4+/-2%/dex with a confidence level higher than 88%. In addition we provide the variation of the SB1 fraction with metallicity separately for F, G, and K spectral types, as well as for dwarf and giant primaries.

We present the second data release (DR2) of the X-Shooter Spectral Library (XSL), which contains all the spectra obtained over the six semesters of that program. This release supersedes our first data release from Chen et al. 2014, with a larger number of spectra (813 observations of 666 stars) and with a more extended wavelength coverage as the data from the near-infrared arm of the X-Shooter spectrograph are now included. The DR2 spectra then consist of three segments that were observed simultaneously and, if combined, cover the range between $\sim$300 nm and $\sim$2.45 $\mu$m at a spectral resolving power close to $R=10\,000$. The spectra were corrected for instrument transmission and telluric absorption, and they were also corrected for wavelength-dependent flux-losses in 85% of the cases. On average, synthesized broad-band colors agree with those of the MILES library and of the combined IRTF and Extended IRTF libraries to within $\sim\!1$%. The scatter in these comparisons indicates typical errors on individual colors in the XSL of 2$-$4 %. The comparison with 2MASS point source photometry shows systematics of up to 5% in some colors, which we attribute mostly to zero-point or transmission curve errors and a scatter that is consistent with the above uncertainty estimates. The final spectra were corrected for radial velocity and are provided in the rest-frame (with wavelengths in air). The spectra cover a large range of spectral types and chemical compositions (with an emphasis on the red giant branch), which makes this library an asset when creating stellar population synthesis models or for the validation of near-ultraviolet to near-infrared theoretical stellar spectra across the Hertzsprung-Russell diagram.

Open clusters (OCs) are recognised as excellent tracers of Galactic thin-disc properties. At variance with intermediate-age and old OCs, for which a significant number of studies is now available, clusters younger than 150 Myr have been mostly overlooked in terms of their chemical composition, with few exceptions. On the other hand, previous investigations seem to indicate an anomalous behaviour of young clusters, which includes slightly sub-solar iron (Fe) abundances and extreme, unexpectedly high barium (Ba) enhancements. In a series of papers, we plan to expand our understanding of this topic and investigate whether these chemical peculiarities are instead related to abundance analysis techniques. We present a new determination of the atmospheric parameters for 23 dwarf stars observed by the Gaia-ESO survey in five young OCs (younger than 150 Myr) and one star-forming region (NGC 2264). We exploit a new method based on titanium (Ti) lines to derive the spectroscopic surface gravity, and most importantly, the microturbulence parameter. A combination of Ti I and Fe I lines is used to obtain effective temperatures. We also infer the abundances of Fe II, Ti II, Na I, Mg I, Al I, Si I, Ca I, Cr I and Ni I. Our findings are in fair agreement with Gaia-ESO iDR5 results for effective temperatures and surface gravities, but suggest that for very young stars, the microturbulence parameter is over-estimated when Fe lines are employed. This affects the derived chemical composition and causes the metal content of very young clusters to be under-estimated. Our clusters display a metallicity [Fe/H] between +0.04 and +0.12; they are not more metal poor than the Sun. Although based on a relatively small sample size, our explorative study suggests that we may not need to call for ad hoc explanations to reconcile the chemical composition of young OCs with Galactic chemical evolution models.

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.

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 X-shooter Spectral Library (XSL) is an empirical stellar library at medium spectral resolution covering the wavelength range from 3000 \AA to 24 800 \AA. This library aims to provide a benchmark for stellar population studies. In this work, we present a uniform set of stellar atmospheric parameters, effective temperatures, surface gravities, and iron abundances for 754 spectra of 616 XSL stars. We used the full-spectrum fitting package ULySS with the empirical MILES library as reference to fit the ultraviolet-blue (UVB) and visible (VIS) spectra. We tested the internal consistency and we compared our results with compilations from the literature. The stars cover a range of effective temperature 2900 < Teff < 38 000 K, surface gravity 0 < log g < 5.7, and iron abundance -2.5 < [Fe/H] < +1.0, with a couple of stars extending down to [Fe/H] = -3.9. The precisions of the measurements for the G- and K-type stars are 0.9%, 0.14, and 0.06 in Teff, log g, and [Fe/H], respectively. For the cool giants with log g < 1, the precisions are 2.1%, 0.21, and 0.22, and for the other cool stars these values are 1%, 0.14, and 0.10. For the hotter stars (Teff > 6500 K), these values are 2.6%, 0.20, and 0.10 for the three parameters.

Pismis 18 is a moderately populated, intermediate-age open cluster located within the solar circle at a Galactocentric distance of about 7 kpc. Few open clusters have been studied in detail in the inner disc region before the Gaia-ESO Survey. New data from the Gaia-ESO Survey allowed us to conduct an extended radial velocity membership study as well as spectroscopic metallicity and detailed chemical abundance measurements for this cluster. Gaia-ESO Survey data for 142 potential members, lying on the upper MS and on the red clump, yielded radial velocity measurements, which, together with proper motion measurements from the Gaia DR2, were used to determine the systemic velocity of the cluster and membership of individual stars. Photometry from Gaia DR2 was used to re-determine cluster parameters based on high confidence member stars only. Cluster abundance measurements of six radial-velocity member stars with UVES high-resolution spectroscopy are presented for 23 elements. According to the new estimates, based on high confidence members, Pismis 18 has an age of $700^{+40}_{-50}$ Myr, interstellar reddening of E(B-V) = $0.562^{+0.012}_{-0.026}$ mag and a de-reddened distance modulus of $DM_0 = 11.96^{+0.10}_{-0.24}$ mag. The median metallicity of the cluster (using the six UVES stars) is [Fe/H] = $+0.23 \pm 0.05$ dex, with [$\alpha$/Fe]= $0.07 \pm 0.13$ and a slight enhancement of s- and r- neutron-capture elements. With the present work, we fully characterized the open cluster Pismis 18, confirming its present location in the inner disc. We estimated a younger age than the previous literature values and gave, for the first time, its metallicity and its detailed abundances. Its [$\alpha$/Fe] and [s-process/Fe], both slightly super-solar, are in agreement with other inner-disc open clusters observed by the Gaia-ESO survey. [abridged]

The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories.

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts.

The X-Shooter Spectral Library (XSL) contains more than 800 spectra of stars across the color-magnitude diagram, that extend from near-UV to near-IR wavelengths (320-2450 nm). We summarize properties of the spectra of O-rich Long Period Variables in XSL, such as phase-related features, and we confront the data with synthetic spectra based on static and dynamical stellar atmosphere models. We discuss successes and remaining discrepancies, keeping in mind the applications to population synthesis modeling that XSL is designed for.

In a previous paper, we assembled a collection of medium-resolution spectra of 35 carbon stars, covering optical and near-infrared wavelengths from 400 to 2400 nm. The sample includes stars from the Milky Way and the Magellanic Clouds, with a variety of $(J-K_s)$ colors and pulsation properties. In the present paper, we compare these observations to a new set of high-resolution synthetic spectra, based on hydrostatic model atmospheres. We find that the broad-band colors and the molecular-band strengths measured by spectrophotometric indices match those of the models when $(J-K_s)$ is bluer than about 1.6, while the redder stars require either additional reddening or dust emission or both. Using a grid of models to fit the full observed spectra, we estimate the most likely atmospheric parameters $T_\mathrm{eff}$, $\log(g)$, $[\mathrm{Fe/H}]$ and C/O. These parameters derived independently in the optical and near-infrared are generally consistent when $(J-K_s)<1.6$. The temperatures found based on either wavelength range are typically within $\pm$100K of each other, and $\log(g)$ and $[\mathrm{Fe/H}]$ are consistent with the values expected for this sample. The reddest stars ($(J-K_s)$ $>$ 1.6) are divided into two families, characterized by the presence or absence of an absorption feature at 1.53\,$\mu$m, generally associated with HCN and C$_2$H$_2$. Stars from the first family begin to be more affected by circumstellar extinction. The parameters found using optical or near-infrared wavelengths are still compatible with each other, but the error bars become larger. In stars showing the 1.53\,$\mu$m feature, which are all large-amplitude variables, the effects of pulsation are strong and the spectra are poorly matched with hydrostatic models. For these, atmospheric parameters could not be derived reliably, and dynamical models are needed for proper interpretation.

We provide a new collection of spectra of 35 carbon stars obtained with the ESO/VLT X-shooter instrument as part of the X-shooter Spectral Library project. The spectra extend from 0.3$\mu$m to 2.4$\mu$m with a resolving power above $\sim$ 8000. The sample contains stars with a broad range of (J-K) color and pulsation properties located in the Milky Way and the Magellanic Clouds. We show that the distribution of spectral properties of carbon stars at a given (J-K) color becomes bimodal (in our sample) when (J-K) is larger than about 1.5. We describe the two families of spectra that emerge, characterized by the presence or absence of the absorption feature at 1.53$\mu$m, generally associated with HCN and C$_2$H$_2$. This feature appears essentially only in large-amplitude variables, though not in all observations. Associated spectral signatures that we interpret as the result of veiling by circumstellar matter, indicate that the 1.53$\mu$m feature might point to episodes of dust production in carbon-rich Miras.

While the near-infrared wavelength regime is becoming more and more important for astrophysics there is a marked lack of spectrophotometric standard star data that would allow the flux calibration of such data. Furthermore, flux calibrating medium- to high-resolution \'echelle spectroscopy data is challenging even in the optical wavelength range, because the available flux standard data are often too coarsely sampled. We will provide standard star reference data that allow users to derive response curves from 300nm to 2500nm for spectroscopic data of medium to high resolution, including those taken with \'echelle spectrographs. In addition we describe a method to correct for moderate telluric absorption without the need of observing telluric standard stars. As reference data for the flux standard stars we use theoretical spectra derived from stellar model atmospheres. We verify that they provide an appropriate description of the observed standard star spectra by checking for residuals in line cores and line overlap regions in the ratios of observed (X-shooter) spectra to model spectra. The finally selected model spectra are then corrected for remaining mismatches and photometrically calibrated using independent observations. The correction of telluric absorption is performed with the help of telluric model spectra.We provide new, finely sampled reference spectra without telluric absorption for six southern flux standard stars that allow the users to flux calibrate their data from 300 nm to 2500 nm, and a method to correct for telluric absorption using atmospheric models.

We present the first release of XSL, the X-Shooter Spectral Library. This release contains 237 stars spanning the wavelengths 3000--10200 \AA\ observed at a resolving power $R \equiv \lambda / \Delta\lambda \sim 10000$. The spectra were obtained at ESO's 8-m Very Large Telescope (VLT). The sample contains O -- M, long-period variable (LPV), C and S stars. The spectra are flux-calibrated and telluric-corrected. We describe a new technique for the telluric correction. The wavelength coverage, spectral resolution and spectral type of this library make it well suited to stellar population synthesis of galaxies and clusters, kinematical investigation of stellar systems and studying the physics of cool stars.