Search

Context: The light element (anti-)correlations shown by globular clusters (GCs) are the main spectroscopic signature of multiple stellar populations. These internal abundance variations provide us with fundamental constraints on the formation mechanism of stellar clusters. Aims: Using Gaia-ESO, the largest and most homogeneous survey of open clusters (OCs), we intend to check whether these stellar aggregates display the same patterns. Based on previous studies of many GCs, several young and massive clusters in the Magellanic Clouds, as well as a few OCs, we do not expect to find any anti-correlation, given the low mass of Milky Way OCs. Methods: We used the results based on UVES spectra of stars in Gaia-ESO to derive the distribution of Na and O abundances and seevwhether they show an unexplained dispersion or whether they are anti-correlated. By selecting only high-probability members with high-precision stellar parameters, we ended up with more than 700 stars in 74 OCs. We examined the O-Na distribution in 28 OCsvwith at least 4 stars available as well as the Na distribution in 24 OCs, with at least 10 stars available. Results: We find that the distribution of Na abundances is compatible with a single-value population, within the errors. The fewvapparent exceptions can be explained by differences in the evolutionary phase (main sequence and giant post first dredge-up episode) or by difficulties in analysing low gravity giants. We did not find any indication of an Na-O anti-correlation in any of the clusters for which O has been derived. Conclusions: Based on the very small spread we find, OCs maintain the status of single stellar populations. However, a definitive answer requires studying more elements and larger samples covering different evolutionary phases. This will be possible with the next generation of large surveys, Comment: Based on the Gaia-ESO Survey. In press on Astronomy & Astrophysics. Figure A.2 updated

The Li abundance observed in pre-main sequence and main sequence late-type stars is strongly age-dependent, but also shows a complex pattern depending on several parameters, such as rotation, chromospheric activity and metallicity. The best way to calibrate these effects, with the aim of studying Li as an age indicator for FGK stars, is to calibrate coeval groups of stars, such as open clusters (OCs) and associations. We present a considerable target sample of 42 OCs and associations, ranging from 1 Myr to 5 Gyr, observed within the Gaia-ESO survey (GES), and using the latest data provided by GES iDR6 and the most recent release of Gaia that was then available, EDR3. As part of this study, we update and improve the membership analysis for all 20 OCs presented in our previous article. We perform detailed membership analyses for all target clusters to identify likely candidates, using all available parameters provided by GES and based on numerous criteria: from radial velocity distributions, to the astrometry and photometry provided by Gaia, to gravity indicators, [Fe/H] metallicity, and Li content. We obtain updated lists of cluster members for the whole target sample, as well as a selection of Li-rich giant contaminants obtained as an additional result of the membership process. Each selection of cluster candidates was thoroughly contrasted with numerous existing membership studies using data from Gaia to ensure the most robust results. These final cluster selections will be used in the third and last paper of this series, which reports the results of a comparative study characterising the observable Li dispersion in each cluster and analysing its dependence on several parameters, allowing us to calibrate a Li-age relation and obtain a series of empirical Li envelopes for key ages in our sample., Comment: 83 pages, 246 figures. To be published in A&A, accepted January 27, 2024

The Gaia-ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted $10^5$ stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ($R\sim47,000$) and GIRAFFE medium ($R\sim20,000$) resolution spectra. The analysis of the Gaia-ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a bayesian Inference method developed specifically for the Gaia-ESO Survey by WG11. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6,231 with typical uncertainties for Teff, log g and [Fe/H] of 32~K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6,188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76,675 with typical uncertainties for Teff, log g and [Fe/H] of 64~K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64,177 of the stars., Comment: 31 pages, 19 figures

We present the first large-scale 3D kinematic study of ~2000 spectroscopically-confirmed young stars (<20 Myr) in 18 star clusters and OB associations (hereafter groups) from the combination of Gaia astrometry and Gaia-ESO Survey spectroscopy. We measure 3D velocity dispersions for all groups, which range from 0.61 to 7.4 km/s (1D velocity dispersions of 0.35 to 4.3 km/s). We find the majority of groups have anisotropic velocity dispersions, suggesting they are not dynamically relaxed. From the 3D velocity dispersions, measured radii and estimates of total mass we estimate the virial state and find that all systems are super-virial when only the stellar mass is considered, but that some systems are sub-virial when the mass of the molecular cloud is taken into account. We observe an approximately linear correlation between the 3D velocity dispersion and the group mass, which would imply that the virial state of groups scales as the square root of the group mass. However, we do not observe a strong correlation between virial state and group mass. In agreement with their virial state we find that nearly all of the groups studied are in the process of expanding and that the expansion is anisotropic, implying that groups were not spherical prior to expansion. One group, Rho Oph, is found to be contracting and in a sub-virial state (when the mass of the surrounding molecular cloud is considered). This work provides a glimpse of the potential of the combination of Gaia and data from the next generation of spectroscopic surveys., Comment: 26 pages, 15 figures, accepted for publication in MNRAS

We present an empirical model of age-dependent photospheric lithium depletion, calibrated using a large, homogeneously-analysed sample of 6200 stars in 52 open clusters, with ages from 2--6000 Myr and $-0.3<{\rm [Fe/H}]<0.2$, observed in the Gaia-ESO spectroscopic survey. The model is used to obtain age estimates and posterior age probability distributions from measurements of the Li I 6708A equivalent width for individual (pre) main sequence stars with $3000 < T_{\rm eff}/{\rm K} <6500$, a domain where age determination from the HR diagram is either insensitive or highly model-dependent. In the best cases, precisions of 0.1 dex in log age are achievable; even higher precision can be obtained for coeval groups and associations where the individual age probabilities of their members can be combined. The method is validated on a sample of exoplanet-hosting young stars, finding agreement with claimed young ages for some, but not others. We obtain better than 10 per cent precision in age, and excellent agreement with published ages, for seven well-studied young moving groups. The derived ages for young clusters ($<1$ Gyr) in our sample are also in good agreement with their training ages, and consistent with several published, model-insensitive lithium depletion boundary ages. For older clusters there remain systematic age errors that could be as large as a factor of two. There is no evidence to link these errors to any strong systematic metallicity dependence of (pre) main sequence lithium depletion, at least in the range $-0.29 < {\rm [Fe/H]} < 0.18$. Our methods and model are provided as software -- "Empirical AGes from Lithium Equivalent widthS" (EAGLES)., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. Code available at https://github.com/robdjeff/eagles

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results., Comment: A&A accepted, minor revision, 36 pages, 38 figures

The Gaia-ESO survey sample of massive OB stars in the Carina Nebula consists of 234 stars. The addition of brighter sources from the Galactic O-Star Spectroscopic Survey and additional sources from the literature allows us to create the most complete census of massive OB stars done so far in the region. It contains a total of 316 stars, being 18 of them in the background and four in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of non-supergiant B type and 6 are of WR or non-O supergiant (II to Ia) spectral class. We identify 20 spectroscopic binary systems with an O-star primary, of which 6 are reported for the first time, and another 18 with a B-star primary, of which 13 are new detections. The average observed double-lined binary fraction of O-type stars in the surveyed region is 0.35, which represents a lower limit. We find a good correlation between the spectroscopic n-qualifier and the projected rotational velocity of the stars. The fraction of candidate runaways among the stars with and without the n-qualifier is 4.4% and 2.4%, respectively, although non resolved double-lined binaries can be contaminating the fast rotators sample., Comment: Accepted for publication in A&A, 62 pages, 8 figures

We test a scenario in which radial migration could affect the Li abundance pattern of dwarf stars in the solar neighbourhood. This may confirm that the Li abundance in these stars can not serve as a probe for the Li abundance in the interstellar medium. We use the high-quality data (including Li abundances) from the 6th internal Data Release of the Gaia-ESO survey. In this sample, we group stars by similarity in chemical abundances via hierarchical clustering. Our analysis treats both measured Li abundances and upper limits. The Li envelope of the previously identified radially migrated stars is well below the benchmark meteoritic value (<3.26 dex); the star with the highest detected abundance has A(Li) = 2.76 dex. This confirms the previous trends observed for old dwarf stars (median ages $\sim$ 8 Gyr), where Li decreases for [Fe/H]$\gtrsim$0. This result acts as supporting evidence that the abundance of Li measured in the upper envelope of old dwarf stars should not be considered a proxy for the interstellar medium Li. Our scenario also indicates that the stellar yields for [M/H]>0 should not be decreased, as recently proposed in the literature. Our study backs the recent studies that claimed that old dwarfs on the hot side of the dip are efficient probes of the ISM abundance of Li, provided atomic diffusion does not lower significantly the initial Li abundance in the atmospheres of metal-rich objects., Comment: 8 pages, 3 figures. Abridged abstract to fit ArXiv's requirements. Letter published in A&A Letters. Version after the language proofs

The spatial distribution of elemental abundances and their time evolution are among the major constraints to disentangle the scenarios of formation and evolution of the Galaxy. We used the sample of open clusters available in the final release of the Gaia-ESO survey to trace the Galactic radial abundance and abundance to iron ratio gradients, and their time evolution. We selected member stars in 62 open clusters, with ages from 0.1 to about 7~Gyr, located in the Galactic thin disc at Galactocentric radii from about 6 to 21~kpc. We analysed the shape of the resulting [Fe/H] gradient, the average gradients [El/H] and [El/Fe] combining elements belonging to four different nucleosynthesis channels, and their individual abundance and abundance ratio gradients. We also investigated the time evolution of the gradients dividing open clusters in three age bins. The[Fe/H] gradient has a slope of -0.054 dex~kpc-1. We saw different behaviours for elements belonging to different channels. We found that the youngest clusters in the inner disc have lower metallicity than their older counterpart and they outline a flatter gradient. We considered some possible explanations, including the effects of gas inflow and migration. We suggested that it might be a bias introduced by the standard spectroscopic analysis producing lower metallicities in low gravity stars. To delineate the shape of the `true' gradient, we should limit our analysis to stars with low surface gravity logg>2.5 and xi<1.8 km~s-1. Based on this reduced sample, we can conclude that the gradient has minimally evolved over the time-frame outlined by the open clusters, indicating a slow and stationary formation of the thin disc in the latest Gyr. We found a secondary role of clusters' migration in shaping the gradient, with a more prominent role of migration for the oldest clusters., Comment: 25 pages, 14 figures and 4 tables in the main text, 3 figures and 7 tables in the Appendix. Accepted for publication in A&A

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that was carried out using the multi-object FLAMES spectrograph at the Very Large Telescope. The survey provides accurate radial velocities, stellar parameters, and elemental abundances for ~115,000 stars in all Milky Way components. In this paper we describe the method adopted in the final data release to derive lithium equivalent widths (EWs) and abundances. Lithium EWs were measured using two different approaches for FGK and M-type stars, to account for the intrinsic differences in the spectra. For FGK stars, we fitted the lithium line using Gaussian components, while direct integration over a predefined interval was adopted for M-type stars. Care was taken to ensure continuity between the two regimes. Abundances were derived using a new set of homogeneous curves of growth that were derived specifically for GES, and which were measured on a synthetic spectral grid consistently with the way the EWs were measured. The derived abundances were validated by comparison with those measured by other analysis groups using different methods. Lithium EWs were measured for ~40,000 stars, and abundances could be derived for ~38,000 of them. The vast majority of the measures (80%) have been obtained for stars in open cluster fields. The remaining objects are stars in globular clusters, or field stars in the Milky Way disc, bulge, and halo. The GES dataset of homogeneous lithium abundances described here will be valuable for our understanding of several processes, from stellar evolution and internal mixing in stars at different evolutionary stages to Galactic evolution., Comment: 15 pages, 15 figures. Accepted by Astronomy & Astrophysics

With its origin coming from several sources (Big Bang, stars, cosmic rays) and given its strong depletion during its stellar lifetime, the lithium element is of great interest as its chemical evolution in the Milky Way is not well understood at present. To help constrain stellar and galactic chemical evolution models, numerous and precise lithium abundances are necessary for a large range of evolutionary stages, metallicities, and Galactic volume. In the age of stellar parametrization on industrial scales, spectroscopic surveys such as APOGEE, GALAH, RAVE, and LAMOST have used data-driven methods to rapidly and precisely infer stellar labels (atmospheric parameters and abundances). To prepare the ground for future spectroscopic surveys such as 4MOST and WEAVE, we aim to apply machine learning techniques to lithium measurements and analyses. We trained a convolution neural network (CNN), coupling Gaia-ESO Survey iDR6 stellar labels (Teff, log(g), [Fe/H], and A(Li)) and GIRAFFE HR15N spectra, to infer the atm parameters and lithium abundances for ~40,000 stars. We show that the CNN properly learns the physics of the stellar labels, from relevant spectral features through a broad range of evolutionary stages and stellar parameters. The Li feature at 6707.8 A is successfully singled out by our CNN, among the thousands of lines. Rare objects such as Li-rich giants are found in our sample. This level of performance is achieved thanks to a meticulously built, high-quality, and homogeneous training sample. The CNN approach is very well adapted for the next generations of spectroscopic surveys aimed at studying (among other elements) lithium, such as the 4MIDABLE-LR/HR (4MOST Milky Way disk and bulge low- and high-resolution) surveys. In this context, the caveats of ML applications should be appropriately investigated, along with the realistic label uncertainties and upper limits for abundances., Comment: A&A, accepted 05 December 2022, 23 pages, 23 figures. Codes and trained models available at https://github.com/SamirNepal/Li\_CNN\_2022

We aim to prepare the machine-learning ground for the next generation of spectroscopic surveys, such as 4MOST and WEAVE. Our goal is to show that convolutional neural networks can predict accurate stellar labels from relevant spectral features in a physically meaningful way. We built a neural network and trained it on GIRAFFE spectra with associated stellar labels from the sixth internal Gaia-ESO data release. Our neural network predicts the atmospheric parameters Teff and log(g) as well as the chemical abundances [Mg/Fe], [Al/Fe], and [Fe/H] for 30115 stellar spectra. The scatter of predictions from eight slightly different network models shows a high internal precision of the network results: 24 K for Teff, 0.03 for log(g), 0.02 dex for [Mg/Fe], 0.03 dex for [Al/Fe], and 0.02 dex for [Fe/H]. The network gradients reveal that the network is inferring the labels in a physically meaningful way from spectral features. Validation with benchmark stars and several scientific applications confirm that our network predictions are accurate for individual stars and recover the properties of different stellar populations in the Milky Way galaxy. Such a study provides very good insights into the application of machine-learning for the spectral analysis of large-scale spectroscopic surveys, such as WEAVE and 4MIDABLE-LR and -HR (4MOST Milky Way disk and bulge low- and high-resolution). The community will have to put a substantial effort into building proactive training sets for machine-learning methods to minimize the possible systematics., Comment: Submitted to A&A

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

We present a spectroscopic analysis of the GIRAFFE and UVES data collected by the Gaia-ESO survey for the young open cluster NGC 3293. Archive spectra from the same instruments obtained in the framework of the `VLT-FLAMES survey of massive stars' are also analysed. Atmospheric parameters, non-LTE chemical abundances for six elements, or variability information are reported for a total of about 160 B stars spanning a wide range in terms of spectral types (B1 to B9.5) and rotation rate (up to 350 km/s). We take advantage of the multi-epoch observations to detect several binary systems or intrinsically line-profile variables. A deconvolution algorithm is used to infer the current, true (deprojected) rotational velocity distribution. We find a broad, Gaussian-like distribution peaking around 200-250 km/s. Although some stars populate the high-velocity tail, most stars in the cluster appear to rotate far from critical. We discuss the chemical properties of the cluster, including the low occurrence of abundance peculiarities in the late B stars and the paucity of objects showing CN-cycle burning products at their surface. We argue that the former result can largely be explained by the inhibition of diffusion effects because of fast rotation, while the latter is generally in accord with the predictions of single-star evolutionary models under the assumption of a wide range of initial spin rates at the onset of main-sequence evolution. However, we find some evidence for a less efficient mixing in two quite rapidly rotating stars that are among the most massive objects in our sample. Finally, we obtain a cluster age of ~20 Myrs through a detailed, star-to-star correction of our results for the effect of stellar rotation. This is significantly older than previous estimates from turn-off fitting that fully relied on classical, non-rotating isochrones. [abridged], Comment: 29 pages, 24 figures, accepted for publication in A&A

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

A growing disquiet has emerged in recent years that standard stellar models are at odds with observations of the colour-magnitude diagrams (CMDs) and lithium depletion patterns of pre main sequence (PMS) stars in clusters. In this work we select 1,246 high probability K/M-type constituent members of 5 young open clusters (5--125\,Myr) in the Gaia-ESO Survey to test a series of models that use standard input physics and others that incorporate surface magnetic fields or cool starspots. We find that: standard models provide systematically under-luminous isochrones for low-mass stars in the CMD and fail to predict Li-depletion of the right strength at the right colour; magnetic models provide better CMD fits with isochrones that are $\sim 1.5-2$ times older, and provide better matches to Li depletion patterns. We investigate how rotation periods, most of which are determined here for the first time from Transiting Exoplanet Survey Satellite data, correlate with CMD position and Li. Among the K-stars in the older clusters we find the brightest and least Li-depleted are the fastest rotators, demonstrating the classic "Li-rotation connection" for the first time at $\sim 35$ Myr in NGC 2547, and finding some evidence that it exists in the early M-stars of NGC 2264 at $<10\,$Myr. However, the wide dispersion in Li depletion observed in fully-convective M-dwarfs in the $\gamma$ Vel cluster at $\sim 20$ Myr appears not to be correlated with rotation and is challenging to explain without a very large ($>10$ Myr) age spread., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society: 28 pages, 14 figures

The chemical makeup of a star provides the fossil information of the environment where it formed. Under this premise, it should be possible to use chemical abundances to tag stars that formed within the same stellar association. This idea - known as chemical tagging - has not produced the expected results, especially within the thin disk where open stellar clusters have chemical patterns that are difficult to disentangle. The ultimate goal of this study is to probe the feasibility of chemical tagging within the thin disk population using high-quality data from a controlled sample of stars. We also aim at improving the existing techniques of chemical tagging and giving guidance on different strategies of clustering analysis in the elemental abundance space. Here we develop the first blind search of open clusters' members through clustering analysis in the elemental abundance space using the OPTICS algorithm applied to data from the Gaia-ESO survey. First, we evaluate different strategies of analysis, determining which ones are more performing. Second, we apply these methods to a data set including both field stars and open clusters attempting a blind recover of as many open clusters as possible. We show how specific strategies of data analysis can improve the final results. Specifically, we demonstrate that open clusters can be more efficaciously recovered with the Manhattan metric and on a space whose dimensions are carefully selected. Using these (and other) prescriptions we are able to recover open clusters hidden in our data set and find new members of these stellar associations. Our results indicate that there are chances of recovering open clusters' members via clustering analysis in the elemental abundance space. Presumably, the performances of chemical tagging will further increase with higher quality data and more sophisticated clustering algorithms., Comment: Accepted for publication by A&A. 21 pages, 11 figures

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters., Comment: 26 pages, 14 figures, accepted for publication in Astronomy & Astrophysics; language-edited version; two appendices merged

Pre-main sequence models with inflated radii are needed to simultaneously reproduce the colour-magnitude diagram and the lithium depletion pattern in young open clusters. We tested a new set of PMS models including radius inflation due to starspots or magnetic inhibition of convection, using five clusters observed by the Gaia-ESO Survey, spanning the age range ~10-100 Myr where such effects could be important. Gaia-ESO radial velocities were combined with Gaia EDR3 astrometry to obtain clean lists of high-probability members for the five clusters. A Bayesian maximum likelihood method was adopted to derive the best model parameters and the cluster reddening and age. Models were calculated for different values of the mixing length parameter ($\alpha_{ML}=2.0$, 1.5 and 1.0), without spots or with effective spot coverage $\beta_{spot}=0.2$ and 0.4. To reproduce the CMD and the Li depletion pattern in Gamma Vel A and B and in 25 Ori we need both a reduced convection efficiency $\alpha_{ML}=1.0$ and an effective spot coverage of ~20%. We obtained ages of 18 Myr and 21 Myr for Gamma Vel A and B, respectively, and 19 Myr for 25 Ori. However, a single isochrone is not sufficient to account for the Li dispersion, and an increasing level of spot coverage as mass decreases seems to be required. The older clusters (NGC2451B at 30 Myr, NGC2547 at 35 Myr, and NGC2516 at 138 Myr) are consistent with standard models ($\alpha_{ML}=2.0$ and no spots) except at low masses: a 20% spot coverage seems to reproduce the sequence of M-type stars better and might explain the observed abundance spread. The quality of Gaia-ESO data combined with Gaia allows us to gain important insights on PMS evolution. Models including starspots can provide a consistent explanation of the cluster sequences and Li abundances of young clusters, although a range of starspot coverage is required to fully reproduce the data., Comment: Accepted for publication in A&A, 18 pages, 10 figures. Revised version after language editing, revised acknowledgements

Context. NGC 1851 is one of several globular clusters for which multiple stellar populations of the subgiant branch have been clearly identified and a difference in metallicity detected. A crucial piece of information on the formation history of this cluster can be provided by the sum of A(C+N+O) abundances. However, these values have lacked a general consensus thus far. The separation of the subgiant branch can be based on age and/or A(C+N+O) abundance differences. Aims. Our main aim was to determine carbon, nitrogen, and oxygen abundances for evolved giants in the globular cluster NGC1851 in order to check whether or not the double populations of stars are coeval. Methods. High-resolution spectra, observed with the FLAMES-UVES spectrograph on the ESO VLT telescope, were analysed using a differential model atmosphere method. Results. We provide abundances of up to 29 chemical elements for a sample of 45 giants in NGC 1851. The investigated stars can be separated into two populations with a difference of 0.07 dex in the mean metallicity, 0.3 dex in the mean C/N, and 0.35 dex in the mean s-process dominated element-to-iron abundance ratios [s/Fe]. No significant difference was determined in the mean values of A(C+N+O) as well as in abundance to iron ratios of carbon, alpha- and iron-peak-elements, and of europium. Conclusions. As the averaged A(C+N+O) values between the two populations do not differ, additional evidence is given that NGC 1851 is composed of two clusters, the metal-rich cluster being by about 0.6 Gyr older than the metal-poor one. A global overview of NGC 1851 properties and the detailed abundances of chemical elements favour its formation in a dwarf spheroidal galaxy that was accreted by the Milky Way., Comment: 13 pages, 17 figures, accepted to Astronomy and Astrophysics

Spectroscopy from the final internal data release of the Gaia-ESO Survey (GES) has been combined with Gaia EDR3 to assign membership probabilities to targets observed towards 63 Galactic open clusters and 7 globular clusters. The membership probabilities are based chiefly on maximum likelihood modelling of the 3D kinematics of the targets, separating them into cluster and field populations. From 43211 observed targets, 13985 are identified as highly probable cluster members ($P>0.9$), with an average membership probability of 0.993. The addition of GES radial velocities successfully drives down the fraction of false positives and we achieve better levels of discrimination in most clusters over the use of astrometric data alone, especially those at larger distances. Since the membership selection is almost purely kinematic, the union of this catalogue with GES and Gaia is ideal for investigating the photometric and chemical properties of clusters as a function of stellar mass, age and Galactic position., Comment: Accepted version for publication in MNRAS. 16 pages + 38 pages of Appendices

It has recently been suggested that all giant stars with mass below 2 $M_{\odot}$ suffer an episode of surface lithium enrichment between the tip of the red giant branch (RGB) and the red clump (RC). We test if the above result can be confirmed in a sample of RC and RGB stars that are members of open clusters. We discuss Li abundances in six open clusters with ages between 1.5 and 4.9 Gyr (turn-off masses between 1.1 and 1.7 $M_{\odot}$). These observations are compared with the predictions of different models that include rotation-induced mixing, thermohaline instability, mixing induced by the first He flash, and energy losses by neutrino magnetic moment. In six clusters, we find about 35\% RC stars with Li abundances that are similar or higher than those of upper RGB stars. This can be a sign of fresh Li production. Because of the extra-mixing episode connected to the luminosity bump, the expectation was for RC stars to have systematically lower surface Li abundances. However, we cannot confirm that the possible Li production is ubiquitous. For about 65\% RC giants we can only determine abundance upper limits that could be hiding very low Li abundances. Our results indicate a possible production of Li during the RC, at levels that would not classify the stars as Li rich. Determination of their carbon isotopic ratio would help to confirm that the RC giants have suffered extra mixing followed by Li enrichment. The Li abundances of the RC stars can be qualitatively explained by the models with an additional mixing episode close to the He flash., Comment: Accepted for publication in A&A, 8 pages, 5 figures, tables available online and under request

Young open clusters (t<200 Myr) have been observed to exhibit several peculiarities in their chemical compositions, from a slightly sub-solar iron content, super-solar abundances of some atomic species (e.g. ionised chromium), and atypical enhancements of [Ba/Fe], with values up to +0.7 dex. Regarding the behaviour of the other $s$-process elements like yttrium, zirconium, lanthanum, and cerium, there is general disagreement in the literature. In this work we expand upon our previous analysis of a sample of five young open clusters (IC2391, IC2602, IC4665, NGC2516, and NGC2547) and one star-forming region (NGC2264), with the aim of determining abundances of different neutron-capture elements, mainly CuI, SrI, SrII, YII, ZrII, BaII, LaII, and CeII. We analysed high-resolution, high signal-to-noise spectra of 23 solar-type stars observed within the \textit{Gaia}-ESO survey. We find that our clusters have solar [Cu/Fe] within the uncertainties, while we confirm the super-solar [Ba/Fe] values (from +0.22 to +0.64 dex). Our analysis also points to mildly enhanced [Y/Fe] values (from 0 and +0.3 dex). For the other $s$-process elements we find that [X/Fe] ratios are solar at all ages. It is not possible to reconcile the anomalous behaviour of Ba and Y at young ages with standard stellar yields and Galactic chemical evolution model predictions. Thus, we explore different possible scenarios related to the behaviour of spectral lines, from the sensitivity to the presence of magnetic fields to the first ionisation potential effect. We also investigate the possibility that they may arise from alterations of the structure of the stellar photosphere due to higher levels of activity in such young stars. We are still unable to explain these enhancements, but we suggest that other elements (i.e. La) might be more reliable tracer of the $s$-process at young ages and encourage further observations., Comment: 26 pages, 14 figures, accepted for publication in A&A

We aim to constrain the mixing processes in low-mass stars by investigating the behaviour of the Li surface abundance after the main sequence. We take advantage of the data from the sixth internal data release of Gaia-ESO, idr6, and from the Gaia Early Data Release 3, edr3. We select a sample of main sequence, sub-giant, and giant stars in which Li abundance is measured by the Gaia-ESO survey, belonging to 57 open clusters with ages from 120~Myr to about 7 Gyr and to Milky Way fields, covering a range in [Fe/H] between -1.0 and +0.5dex. We study the behaviour of the Li abundances as a function of stellar parameters. We compare the observed Li behaviour in field giant stars and in giant stars belonging to individual clusters with the predictions of a set of classical models and of models with mixing induced by rotation and thermohaline instability. The comparison with stellar evolution models confirms that classical models cannot reproduce the lithium abundances observed in the metallicity and mass regimes covered by the data. The models that include the effects of both rotation-induced mixing and thermohaline instability account for the Li abundance trends observed in our sample, in all metallicity and mass ranges. The differences between the results of the classical models and of the rotation models largely differ (up to ~2 dex), making lithium the best element to constrain stellar mixing processes in low-mass stars. For stars with well-determined masses, we find a better agreement between observed surface abundances and models with rotation-induced and thermohaline mixings, the former dominating during the main sequence and the first phases of the post-main sequence evolution and the latter after the bump in the luminosity function., Comment: Accepted for publication in A&A, 19 pages, 16 figures

Astrometry and photometry from {\it Gaia} and spectroscopic data from the {\it Gaia}-ESO Survey (GES) are used to identify the lithium depletion boundary (LDB) in the young cluster NGC 2232. A specialised spectral line analysis procedure was used to recover the signature of undepleted lithium in very low luminosity cluster members. An age of $38\pm 3$ Myr is inferred by comparing the LDB location in absolute colour-magnitude diagrams (CMDs) with the predictions of standard models. This is more than twice the age derived from fitting isochrones to low-mass stars in the CMD with the same models. Much closer agreement between LDB and CMD ages is obtained from models that incorporate magnetically suppressed convection or flux-blocking by dark, magnetic starspots. The best agreement is found at ages of $45-50$\,Myr for models with high levels of magnetic activity and starspot coverage fractions $>50$ per cent, although a uniformly high spot coverage does not match the CMD well across the full luminosity range considered.

Previous studies of open clusters have shown that lithium depletion is not only strongly age dependent but also shows a complex pattern with other parameters that is not yet understood. For pre- and main-sequence late-type stars, these parameters include metallicity, mixing mechanisms, convection structure, rotation, and magnetic activity. We perform a thorough membership analysis for a large number of stars observed within the Gaia-ESO survey (GES) in the field of 20 open clusters, ranging in age from young clusters and associations, to intermediate-age and old open clusters. Based on the parameters derived from the GES spectroscopic observations, we obtained lists of candidate members for each of the clusters in the sample by deriving RV distributions and studying the position of the kinematic selections in the EW(Li) versus Teff plane to obtain lithium members. We used gravity indicators to discard field contaminants and studied [Fe/H] metallicity to further confirm the membership of the candidates. We also made use of studies using recent data from the Gaia DR1 and DR2 releases to assess our member selections. We identified likely member candidates for the sample of 20 clusters observed in GES (iDR4) with UVES and GIRAFFE, and conducted a comparative study that allowed us to characterize the properties of these members, as well as identify field contaminant stars, both lithium-rich giants and non-giant outliers. This work is the first step towards the calibration of the lithium-age relation and its dependence on other GES parameters. During this project we aim to use this relation to infer the ages of GES field stars, and identify their potential membership to young associations and stellar kinematic groups of different ages., Comment: 264 pages, 105 figures. To be published in A&A, accepted 29th July 2020

The extensive stellar spectroscopic datasets that are available for studies in Galactic Archeaology thanks to, for example, the Gaia-ESO Survey, now benefit from having a significant number of targets that overlap with asteroseismology projects such as Kepler, K2 and CoRoT. Combining the measurements from spectroscopy and asteroseismology allows us to attain greater accuracy with regard to the stellar parameters needed to characterise the stellar populations of the Milky Way. The aim of this Gaia-ESO Survey special project is to produce a catalogue of self-consistent stellar parameters by combining measurements from high-resolution spectroscopy and precision asteroseismology. We carried out an iterative analysis of 90 K2@Gaia-ESO red giants. The spectroscopic values of Teff were used as input in the seismic analysis to obtain log(g) values. The seismic estimates of log(g) were then used to re-determine the spectroscopic values of Teff and [Fe/H]. Only one iteration was required to obtain parameters that are in good agreement for both methods and thus, to obtain the final stellar parameters. A detailed analysis of outliers was carried out to ensure a robust determination of the parameters. The results were then combined with Gaia DR2 data to compare the seismic log(g) with a parallax-based log(g) and to investigate instances of variations in the velocity and possible binaries within the dataset. This analysis produced a high-quality catalogue of stellar parameters for 90 red giant stars observed by both K2 and Gaia-ESO that were determined through iterations between spectroscopy and asteroseismology. We compared the seismic gravities with those based on Gaia parallaxes to find an offset which is similar to other studies that have used asteroseismology. Our catalogue also includes spectroscopic chemical abundances and radial velocities, as well as indicators for possible binary detections., Comment: 33 pages, 26 figures

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., Comment: 5 pages, 3 figures, accepted for publication in Astronomy & Astrophysics Letters

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., Comment: Accepted for publication in MNRAS

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 47 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., Comment: 18 pages, 15 figures, to be published in A&A

Lithium rich stars in globular clusters are rare. In fact, only 14 have been found so far, in different evolutionary phases from dwarfs to giants. Different mechanisms have been proposed to explain this enhancement, but it is still an open problem. Using spectra collected within the Gaia-ESO Survey, obtained with the GIRAFFE spectrograph at the ESO Very Large Telescope, we present the discovery of the first Li-rich star in the cluster NGC 1261, the second star known in the red giant branch bump phase. The star shows an extreme Li overabundance of A(Li)_LTE=3.92\pm0.14, corresponding to A(Li)_NLTE=3.40 dex. We propose that the Li enhancement is caused by fresh Li production through an extra mixing process (sometimes referred to as {\em cool bottom burning}) or could be a pre-existing Li overabundance resulting from binary mass transfer, likely from a red giant branch star, because of the low barium abundance. To unambiguously explain the Li enhancement in globular cluster stars, however, a reliable determination of the abundance of key species like Be, 6Li, 12C/13C, and several s-process elements is required, as well as detailed modeling of chromospheric activity indicators., Comment: 5 pages, 3 figures, accepted for publication on A&ALetter

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

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., Comment: 16 pages, 11 figures, accepted by A&A

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., Comment: 28 pages, 24 figures. Accepted by A&A

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], Comment: 15 pages, 10 figures, 8 tables (Table 2 available at the CDS)

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., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in A&A; 23 pages, 17 figures, 3 tables

The abundance ratio N/O is a useful tool to study the interplay of galactic processes, e.g. star formation efficiency, time-scale of infall and outflow loading factor We aim to trace log(N/O) versus [Fe/H] in the Milky Way and to compare it with a set of chemical evolution models to understand the role of infall, outflow and star formation efficiency in the building-up of the Galactic disc. We use the abundances from idr2-3, idr4, idr5 data releases of the Gaia-ESO Survey both for Galactic field and open cluster stars.We determine membership and average composition of open clusters and we separate thin and thick disc field stars.We consider the effect of mixing in the abundance of N in giant stars. We compute a grid of chemical evolution models, suited to reproduce the main features of our Galaxy, exploring the effects of the star formation efficiency, the infall time-scale and the differential outflow. With our samples, we map the metallicity range -0.6<[Fe/H]<0.3 with a corresponding -1.2

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

The GES survey using FLAMES at the VLT has obtained high-resolution UVES spectra for a large number of giant stars, allowing a determination of the abundances of the key chemical elements C and N at their surface. The surface abundances of these chemical species are well-known to change in stars during their evolution on the red giant branch after the first dredge-up episod, as a result of extra-mixing phenomena. We investigate the effects of thermohaline mixing on C and N abundances using the first comparison between the GES [C/N] determinations with simulations of the observed fields using a model of stellar population synthesis. We explore the effects of thermohaline mixing on the chemical properties of giants through stellar evolutionary models computed with the stellar evolution code STAREVOL. We include these stellar evolution models in the Besan\c{c}on Galaxy model to simulate the [C/N] distributions determined from the UVES spectra of the GES and compare them with the observations. Theoretical predictions including the effect of thermohaline mixing are in good agreement with the observations. However, the field stars in the GES with C and N-abundance measurements have a metallicity close to solar, where the efficiency of thermohaline mixing is not very large. The C and N abundances derived by the GES in open and globular clusters clearly show the impact of thermohaline mixing at low-metallicity, allowing to explain the [C/N] ratio observed in lower-mass and older giant stars. Using independent observations of carbon isotopic ratio in clump field stars and open clusters, we also confirm that thermohaline mixing should be taken into account to explain the behavior of 12C/13C ratio as a function of stellar age. Overall the current model including thermohaline mixing is able to reproduce very well the C- and N-abundances over the whole metallicity range investigated by the GES data., Comment: Accepted for publication in A&A ; 12 pages, 12 figures

We report 20 new lithium-rich giants discovered within the Gaia-ESO Survey, including the first Li-rich giant with evolutionary stage confirmed by CoRoT data. Atmospheric parameters and abundances were derived in model atmosphere analyses using medium-resolution GIRAFFE or high-resolution UVES spectra. These results are part of the fifth internal data release of Gaia-ESO. The Li abundances were corrected for non-LTE effects. We used Gaia DR2 parallaxes to estimate distances and luminosities. The giants have A(Li) > 2.2 dex. The majority of them (14 out of 20 stars) are in the CoRoT fields. Four giants are located in the field of three open clusters but are not members. Two giants were observed in fields towards the Galactic bulge but are likely in the inner disk. One of the bulge field giants is super Li-rich with A(Li) = 4.0 dex. We identified one giant with infrared excess at 22 microns. Two other giants, with large vsin i, might be Li-rich because of planet engulfment. Another giant is found to be barium enhanced and thus could have accreted material from a former AGB companion. Otherwise, besides the Li enrichment, the evolutionary stages are the only other connection between these new Li-rich giants. The CoRoT data confirm that one Li-rich giant is at the core-He burning stage. The other giants are concentrated in close proximity to the RGB luminosity bump, the core-He burning stages, or the early-AGB. This is very clear when looking at the Gaia-based luminosities of the Li-rich giants. This is also seen when the CoRoT Li-rich giants are compared to a larger sample of 2252 giants observed in the CoRoT fields by the Gaia-ESO Survey, which are distributed all over the RGB in the Teff-logg diagram. These observations show that evolutionary stage is a major factor behind the Li enrichment in giants. Other processes, like planet accretion, contribute to a smaller scale. [abridged], Comment: accepted for publication in A&A, new version after the language corrections

Investigating the chemical homogeneity of stars born from the same molecular cloud at virtually the same time is very important for our understanding of the chemical enrichment of the interstellar medium and with it the chemical evolution of the Galaxy. One major cause of inhomogeneities in the abundances of open clusters is stellar evolution of the cluster members. In this work, we investigate variations in the surface chemical composition of member stars of the old open cluster M67 as a possible consequence of atomic diffusion effects taking place during the main-sequence phase. The abundances used are obtained from high-resolution UVES/FLAMES spectra within the framework of the Gaia-ESO Survey. We find that the surface abundances of stars on the main sequence decrease with increasing mass reaching a minimum at the turn-off. After deepening of the convective envelope in sub-giant branch stars, the initial surface abundances are restored. We found the measured abundances to be consistent with the predictions of stellar evolutionary models for a cluster with the age and metallicity of M67. Our findings indicate that atomic diffusion poses a non-negligible constraint on the achievable precision of chemical tagging methods., Comment: 14 pages, 7 figures, accepted by MNRAS

Context. The origin and dynamical evolution of star clusters is an important topic in stellar astrophysics. Several models have been proposed to understand the formation of bound and unbound clusters and their evolution, and these can be tested by examining the kinematical and dynamical properties of clusters over a wide range of ages and masses. Aims. We use the Gaia-ESO Survey products to study four open clusters (IC 2602, IC 2391, IC 4665, and NGC 2547) that lie in the age range between 20 and 50 Myr. Methods. We employ the gravity index $\gamma$ and the equivalent width of the lithium line at 6708 $\AA$, together with effective temperature $\rm{T_{eff}}$, and the metallicity of the stars in order to discard observed contaminant stars. Then, we derive the cluster radial velocity dispersions $\sigma_c$, the total cluster mass $\rm{M}_{tot}$, and the half mass radius $r_{hm}$. Using the $Gaia$-DR1 TGAS catalogue, we independently derive the intrinsic velocity dispersion of the clusters from the astrometric parameters of cluster members. Results. The intrinsic radial velocity dispersions derived by the spectroscopic data are larger than those derived from the TGAS data, possibly due to the different masses of the considered stars. Using $\rm{M}_{tot}$ and $r_{hm}$ we derive the virial velocity dispersion $\sigma_{vir}$ and we find that three out of four clusters are supervirial. This result is in agreement with the hypothesis that these clusters are dispersing, as predicted by the "residual gas expulsion" scenario. However, recent simulations show that the virial ratio of young star clusters may be overestimated if it is determined using the global velocity dispersion, since the clusters are not fully relaxed., Comment: 14 pages, 6 figures

We describe the methodologies that, taking advantage of Gaia-DR1 and the Gaia-ESO Survey data, enable the comparison of observed open star cluster sequences with stellar evolutionary models. The final, long-term goal is the exploitation of open clusters as age calibrators. We perform a homogeneous analysis of eight open clusters using the Gaia-DR1 TGAS catalogue for bright members, and information from the Gaia-ESO Survey for fainter stars. Cluster membership probabilities for the Gaia-ESO Survey targets are derived based on several spectroscopic tracers. The Gaia-ESO Survey also provides the cluster chemical composition. We obtain cluster parallaxes using two methods. The first one relies on the astrometric selection of a sample of bona fide members, while the other one fits the parallax distribution of a larger sample of TGAS sources. Ages and reddening values are recovered through a Bayesian analysis using the 2MASS magnitudes and three sets of standard models. Lithium depletion boundary (LDB) ages are also determined using literature observations and the same models employed for the Bayesian analysis. For all but one cluster, parallaxes derived by us agree with those presented in Gaia Collaboration et al. (2017), while a discrepancy is found for NGC 2516; we provide evidence supporting our own determination. Inferred cluster ages are robust against models and are generally consistent with literature values. The systematic parallax errors inherent in the Gaia DR1 data presently limit the precision of our results. Nevertheless, we have been able to place these eight clusters onto the same age scale for the first time, with good agreement between isochronal and LDB ages where there is overlap. Our approach appears promising and demonstrates the potential of combining Gaia and ground-based spectroscopic datasets., Comment: Accepted for publication by A&A. 22 pages main text plus 8 pages appendices, 12 figures main text, 7 figures appendices

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

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

We explore the structure and star formation history of the open cluster NGC 2264 (~3 Myr). We combined spectroscopic data from the Gaia-ESO Survey (GES) with multi-wavelength photometry from the Coordinated Synoptic Investigation of NGC 2264 (CSI 2264). We examined a sample of 655 cluster members, including both disk-bearing and disk-free young stars. We find a significant age spread of 4-5 Myr among cluster members. Disk-bearing objects are statistically associated with younger isochronal ages than disk-free sources. The cluster has a hierarchical structure, with two main blocks. The northern half develops around the O-type binary star S Mon; the southern half, close to the tip of the Cone Nebula, contains the most embedded regions of NGC 2264, populated mainly by objects with disks and ongoing accretion. The median ages of objects at different locations within the cluster, and the spatial distribution of disked and non-disked sources, suggest that star formation began in the north of the cluster, over 5 Myr ago, and was ignited in its southern region a few Myr later. Star formation is likely still ongoing in the most embedded regions of the cluster, while the outer regions host a widespread population of more evolved objects. We find a detectable lag between the typical age of disk-bearing objects and that of accreting objects in the inner regions of NGC 2264: the first tend to be older than the second, but younger than disk-free sources at similar locations within the cluster. This supports earlier findings that the characteristic timescales of disk accretion are shorter than those of disk dispersal, and smaller than the average age of NGC 2264 (i.e., < 3 Myr). At the same time, disks in the north of the cluster tend to be shorter-lived (~2.5 Myr) than elsewhere; this may reflect the impact of massive stars within the region (notably S Mon), that trigger rapid disk dispersal., Comment: 24 pages, 14 figures, 7 tables; accepted for publication in A&A. Abstract shortened for arXiv listing

The typical methodology for comparing simulated galaxies with observational surveys is usually to apply a spatial selection to the simulation to mimic the region of interest covered by a comparable observational survey sample. In this work we compare this approach with a more sophisticated post-processing in which the observational uncertainties and selection effects (photometric, surface gravity and effective temperature) are taken into account. We compare a `solar neighbourhood analogue' region in a model Milky Way-like galaxy simulated with RAMSES-CH with fourth release Gaia-ESO survey data. We find that a simple spatial cut alone is insufficient and that observational uncertainties must be accounted for in the comparison. This is particularly true when the scale of uncertainty is large compared to the dynamic range of the data, e.g. in our comparison, the [Mg/Fe] distribution is affected much more than the more accurately determined [Fe/H] distribution. Despite clear differences in the underlying distributions of elemental abundances between simulation and observation, incorporating scatter to our simulation results to mimic observational uncertainty produces reasonable agreement. The quite complete nature of the Gaia-ESO survey means that the selection function has minimal impact on the distribution of observed age and metal abundances but this would become increasingly more important for surveys with narrower selection functions., Comment: Accepted for publication in MNRAS; 13 pages, 7 figures, 1 table

The Gaia-ESO Survey (GES) is a large spectroscopic survey that provides a unique opportunity to study the distribution of spectroscopic multiple systems among different populations of the Galaxy. We aim at detecting binarity/multiplicity for stars targeted by the GES from the analysis of the cross-correlation functions (CCFs) of the GES spectra with spectral templates. We develop a method based on the computation of the CCF successive derivatives to detect multiple peaks and determine their radial velocities, even when the peaks are strongly blended. The parameters of the detection of extrema (DOE) code have been optimized for each GES GIRAFFE and UVES setup to maximize detection. This code therefore allows to automatically detect multiple line spectroscopic binaries (SBn, n>1). We apply this method on the fourth GES internal data release and detect 354 SBn candidates (342 SB2, 11 SB3 and even one SB4), including only 9 SB2 known in the literature. This implies that about 98% of these SBn candidates are new (because of their faint visual magnitude that can reach V=19). Visual inspection of the SBn candidate spectra reveals that the most probable candidates have indeed a composite spectrum. Among SB2 candidates, an orbital solution could be computed for two previously unknown binaries: 06404608+0949173 (known as V642 Mon) in NGC 2264 and 19013257-0027338 in Berkeley 81. A detailed analysis of the unique SB4 (four peaks in the CCF) reveals that HD 74438 in the open cluster IC 2391 is a physically bound stellar quadruple system. The SB candidates belonging to stellar clusters are reviewed in detail to discard false detections. We warn against the use of atmospheric parameters for these system components rather than by SB-specific pipelines. Our implementation of an automatic detection of spectroscopic binaries within the GES has allowed an efficient discovery of many new multiple systems., Comment: 36 pages, 30 figures, 12 tables. Accepted for publication by A&A on 2017-07-05

We present a spectroscopic study of the dynamics of the ionized and neutral gas throughout the Lagoon nebula (M8), using VLT/FLAMES data from the Gaia-ESO Survey. We explore the connections between the nebular gas and the stellar population of the associated star cluster NGC6530. We characterize through spectral fitting emission lines of H-alpha, [N II] and [S II] doublets, [O III], and absorption lines of sodium D doublet, using data from the FLAMES/Giraffe and UVES spectrographs, on more than 1000 sightlines towards the entire face of the Lagoon nebula. Gas temperatures are derived from line-width comparisons, densities from the [S II] doublet ratio, and ionization parameter from H-alpha/[N II] ratio. Although doubly-peaked emission profiles are rarely found, line asymmetries often imply multiple velocity components along the line of sight. This is especially true for the sodium absorption, and for the [O III] lines. Spatial maps for density and ionization are derived, and compared to other known properties of the nebula and of its massive stars 9 Sgr, Herschel 36 and HD 165052 which are confirmed to provide most of the ionizing flux. The detailed velocity fields across the nebula show several expanding shells, related to the cluster NGC6530, the O stars 9 Sgr and Herschel 36, and the massive protostar M8East-IR. The origins of kinematical expansion and ionization of the NGC6530 shell appear to be different. We are able to put constrains on the line-of-sight (relative or absolute) distances between some of these objects and the molecular cloud. The large obscuring band running through the middle of the nebula is being compressed by both sides, which might explain its enhanced density. We also find an unexplained large-scale velocity gradient across the entire nebula. At larger distances, the transition from ionized to neutral gas is studied using the sodium lines., Comment: 26 pages, 31 figures, accepted on Astronomy and Astrophysics journal