Gaia-ESO Survey Research Articles

Using neural network models to estimate stellar ages from lithium equivalent widths: an eagles expansion

ABSTRACT We present an artificial neural network (ANN) model of photospheric lithium depletion in cool stars ($3000\lt T_{\rm eff}/{\rm K} \lt 6500$), producing estimates and probability distributions of age from $^7$Li 6708 Å equivalent width (LiEW) and effective temperature data inputs. The model is trained on the same sample of 6200 stars from 52 open clusters, observed in the Gaia-ESO spectroscopic survey, and used to calibrate the previously published analytical eagles model, with ages 2–6000 Myr and $-0.3 \lt $ [Fe/H] $\lt 0.2$. The additional flexibility of the ANN provides some improvements, including better modelling of the ‘lithium dip’ at ages $\lt 50$ Myr and $T_{\rm eff}\sim 3500$ K, and of the intrinsic dispersion in LiEW at all ages. Poor age discrimination is still an issue at ages >1 Gyr, confirming that additional modelling flexibility is not sufficient to fully represent the LiEW–age–T$_{\text{eff}}$ relationship, and suggesting the involvement of further astrophysical parameters. Expansion to include such parameters–rotation, accretion, and surface gravity–is discussed, and the use of an ANN means these can be more easily included in future iterations, alongside more flexible functional forms for the LiEW dispersion. Our methods and ANN model are provided in an updated version 2.0 of the eagles software.

The Gaia-ESO Survey: No sign of multiple stellar populations in open clusters from their sodium and oxygen abundances

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 see whether 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 OCs with 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 few apparent 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.

The Gaia-ESO Survey: Calibrating the lithium–age relation with open clusters and associations

Context. 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, and 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. Aims. We present a considerable target sample of 42 OCs and associations – with an age range 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. Methods. We perform detailed membership analyses for all target clusters to identify likely candidates, using all available parameters provided by GES, complemented with detailed bibliographical searches, and based on numerous criteria: from radial velocity distributions, to the astrometry (proper motions and parallaxes) and photometry provided by Gaia, to gravity indicators (log g and the γ index), [Fe/H] metallicity, and Li content in diagrams of (Li equivalent widths) EW(Li) versus Teff. Results. 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. Conclusions. 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.

The Gaia-ESO Survey: The DR5 analysis of the medium-resolution GIRAFFE and high-resolution UVES spectra of FGK-type stars

The Gaia-ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted 105 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 ~ 47 000) and GIRAFFE medium (R ~ 20 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 method was also used for the chemical abundance homogenisation within WG11, however, the WG10 chemical abundance data set was too sparsely populated so basic corrections for each node analysis were employed for the homogenisation instead. 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 64177 of the stars.

Observational constraints on the origin of the elements – VII. NLTE analysis of Y ii lines in spectra of cool stars and implications for Y as a Galactic chemical clock

ABSTRACT Yttrium (Y), a key s-process element, is commonly used in nucleosynthesis studies and as a Galactic chemical clock when combined with magnesium (Mg). We study the applicability of the previously assumed local thermal equilibrium (LTE) line formation assumption in Y abundance studies of main-sequence and red giant stars, and probe the impact of NLTE (non-LTE) effects on the [Y/Mg] ratio, a proposed stellar age indicator. We derive stellar parameters, ages, and NLTE abundances of Fe, Mg, and Y for 48 solar analogue stars from high-resolution spectra acquired within the Gaia-ESO survey. For Y, we present a new NLTE atomic model. We determine a solar NLTE abundance of A(Y)NLTE = 2.12 ± 0.04 dex, 0.04 dex higher than LTE. NLTE effects on Y abundance are modest for optical Y ii lines, which are frequently used in Sun-like stars diagnostics. NLTE has a small impact on the [Y/Mg] ratio in such stars. For metal-poor red giants, NLTE effects on Y ii lines are substantial, potentially exceeding +0.5 dex. For the Gaia/4MOST/WEAVE benchmark star, HD 122563, we find the NLTE abundance ratio of [Y/Fe]NLTE = −0.55 ± 0.04 dex with consistent abundances obtained from different Y ii lines. NLTE has a differential effect on Y abundance diagnostics in late-type stars. They notably affect Y ii lines in red giants and very metal-poor stars, which are typical Galactic enrichment tracers of neutron-capture elements. For main-sequence stars, NLTE effects on optical diagnostic Y ii lines remain minimal across metallicities. This affirms the [Y/Mg] ratio’s reliability as a cosmochronometer for Sun-like stars.

The Gaia-ESO Survey: Homogenisation of stellar parameters and elemental abundances

The Gaia-ESO Survey is a public spectroscopic survey that targeted ≳105 stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its final public 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 (Teff ~ 30 000–52 000 K) all the way to K-M (≳3500 K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different types. 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, for example, 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, each specialised in one or more of 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 in order to arrive at a single set of recommended results for use by all survey collaborators. We also present some general quality and consistency checks performed on the entirety of the survey results.

Origin of neutron-capture elements with the Gaia-ESO survey: the evolution of s- and r-process elements across the Milky Way

ABSTRACT We investigate the origin of neutron-capture elements by analysing their abundance patterns and radial gradients in the Galactic thin disc. We adopt a detailed two-infall chemical evolution model for the Milky Way, including state-of-the-art nucleosynthesis prescriptions for neutron-capture elements. We consider r-process nucleosynthesis from merging neutron stars (MNS) and magneto-rotational supernovae (MR-SNe), and s-process synthesis from low- and intermediate-mass stars (LIMS) and rotating massive stars. The predictions of our model are compared with data from the sixth data release of the Gaia-ESO survey, from which we consider 62 open clusters with age ≳ 0.1 Gyr and ∼1300 Milky Way disc field stars. We conclude that: (i) the [Eu/Fe] versus [Fe/H] diagram is reproduced by both prompt and delayed sources, with the prompt source dominating Eu production; (ii) rotation in massive stars significantly contributes to the first peak s-process elements, but MNS and MR-SNe are necessary to match the observations; and (iii) our model slightly underpredicts Mo and Nd, while accurately reproducing the [Pr/Fe] versus [Fe/H] trend. Regarding the radial gradients, we find that: (i) our predicted [Fe/H] gradient slope agrees with observations from Gaia-ESO and other high-resolution spectroscopic surveys; (ii) the predicted [Eu/H] radial gradient slope is steeper than the observed one, regardless of how quick the production of Eu is, prompting discussion on different Galaxy-formation scenarios and stellar radial migration effects; and (iii) elements in the second s-process peak as well as Nd and Pr exhibit a plateau at low-Galactocentric distances, likely due to enhanced enrichment from LIMS in the inner regions.

StarHorse results for spectroscopic surveys and Gaia DR3: Chrono-chemical populations in the solar vicinity, the genuine thick disk, and young alpha-rich stars

The Gaia mission has provided an invaluable wealth of astrometric data for more than a billion stars in our Galaxy. The synergy between Gaia astrometry, photometry, and spectroscopic surveys gives us comprehensive information about the Milky Way. Using the Bayesian isochrone-fitting code StarHorse, we derive distances and extinctions for more than 10 million unique stars listed in both Gaia Data Release 3 and public spectroscopic surveys: 557 559 in GALAH+ DR3, 4 531 028 in LAMOST DR7 LRS, 347 535 in LAMOST DR7 MRS, 562 424 in APOGEE DR17, 471 490 in RAVE DR6, 249 991 in SDSS DR12 (optical spectra from BOSS and SEGUE), 67 562 in the Gaia-ESO DR5 survey, and 4 211 087 in the Gaia RVS part of the Gaia DR3 release. StarHorse can increase the precision of distance and extinction measurements where Gaia parallaxes alone would be uncertain. We used StarHorse for the first time to derive stellar ages for main-sequence turnoff and subgiant branch stars, around 2.5 million stars, with age uncertainties typically around 30%; the uncertainties drop to 15% for subgiant-branch-only stars, depending on the resolution of the survey. With the derived ages in hand, we investigated the chemical-age relations. In particular, the α and neutron-capture element ratios versus age in the solar neighbourhood show trends similar to previous works, validating our ages. We used the chemical abundances from local subgiant samples of GALAH DR3, APOGEE DR17, and LAMOST MRS DR7 to map groups with similar chemical compositions and StarHorse ages, using the dimensionality reduction technique t-SNE and the clustering algorithm HDBSCAN. We identify three distinct groups in all three samples, confirmed by their kinematic properties: the genuine chemical thick disk, the thin disk, and a considerable number of young alpha-rich stars (427) that are also a part of the delivered catalogues. We confirm that the genuine thick disk’s kinematics and age properties are radically different from those of the thin disk and compatible with high-redshift (z ≈ 2) star-forming disks with high dispersion velocities. We also find a few extra chemical populations in GALAH DR3 thanks to the availability of neutron-capture element information.

The Gaia-ESO Survey: empirical estimates of stellar ages from lithium equivalent widths (eagles)

ABSTRACT 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 to 6000 Myr and −0.3 < [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 6708 Å equivalent width for individual (pre) main-sequence stars with 3000 < Teff/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 2. 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 < [Fe/H] < 0.18. Our methods and model are provided as software – ‘Empirical AGes from Lithium Equivalent widthS’ (eagles).

LOTUS: A (Non-) LTE Optimization Tool for Uniform Derivation of Stellar Atmospheric Parameters

Precise fundamental atmospheric stellar parameters and abundance determination of individual elements in stars are important for all stellar population studies. Non–local thermodynamic equilibrium (non-LTE; hereafter NLTE) models are often important for such high precision, however, can be computationally complex and expensive, which renders the models less utilized in spectroscopic analyses. To alleviate the computational burden of such models, we developed a robust 1D, NLTE fundamental atmospheric stellar parameter derivation tool, LOTUS, to determine the effective temperature T eff, surface gravity , metallicity [Fe/H], and microturbulent velocity v mic for FGK-type stars, from equivalent width (EW) measurements of Fe i and Fe ii lines. We utilize a generalized curve of growth method to take into account the EW dependencies of each Fe i and Fe ii line on the corresponding atmospheric stellar parameters. A global differential evolution optimization algorithm is then used to derive the fundamental parameters. Additionally, LOTUS can determine precise uncertainties for each stellar parameter using a Markov Chain Monte Carlo algorithm. We test and apply LOTUS on a sample of benchmark stars, as well as stars with available asteroseismic surface gravities from the K2 survey, and metal-poor stars from the Gaia-ESO and R-Process Alliance surveys. We find very good agreement between our NLTE-derived parameters in LOTUS to nonspectroscopic values on average within T eff = ±30 K, and = ±0.10 dex for benchmark stars. We provide open access of our code, as well as of the interpolated precomputed NLTE EW grids available on Github (the software is available on GitHub 3 3 https://github.com/Li-Yangyang/LOTUS under an MIT License, and version 0.1.1 (as the persistent version) is archived in Zenodo) and documentation with working examples on the Readthedocs book.

Gaia-ESO survey: Massive stars in the Carina Nebula

Context. The Carina Nebula is one of the major massive star-forming regions in the Galaxy. Its relatively nearby distance (2.35 kpc) makes it an ideal laboratory for the study of massive star formation, structure, and evolution, both for individual stars and stellar systems. Thanks to the high-quality spectra provided by the Gaia-ESO survey and the LiLiMaRlin library, as well as Gaia EDR3 astrometry, a detailed and homogeneous spectroscopic characterization of its massive stellar content can be carried out. Aims. Our main objective is to spectroscopically characterize all massive members of the Carina Nebula in the Gaia-ESO survey footprint to provide an updated census of massive stars in the region and an updated estimate of the binary fraction of O stars. Methods. We performed accurate spectral classification using an interactive code that compares spectra with spectral libraries of OB standard stars, as well as line-based classic methods. We calculated membership using our own algorithm based on Gaia EDR3 astrometry. To check the correlation between the spectroscopic n-qualifier and the rotational velocity, we used a semi-automated tool for the line-broadening characterization of OB stars based on a combined Fourier transform and goodness-of-fit methodology. Results. 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 so far in the region. It contains a total of 316 stars, with 18 of them in the background and 4 in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of nonsupergiant B type, and 6 are of WR or nonO 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 nonresolved double-lined binaries could be contaminating the sample of fast rotators.

TheGaia-ESO Survey: Preparing the ground for 4MOST and WEAVE galactic surveys

Context.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.Aims.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.Methods.We trained a convolution neural network (CNN), couplingGaia-ESO Survey iDR6 stellar labels (Teff, log(g), [Fe/H], and A(Li)) and GIRAFFE HR15N spectra, to infer the atmospheric parameters and lithium abundances for ∼40 000 stars. The CNN architecture and accompanying notebooks are available online via GitHub.Results.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 lithium feature at 6707.8 Å is successfully singled out by our CNN, among the thousands of lines in the GIRAFFE HR15N setup. Rare objects such as lithium-rich giants are found in our sample. This level of performance is achieved thanks to a meticulously built, high-quality, and homogeneous training sample.Conclusions.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 machine-learning applications should be appropriately investigated, along with the realistic label uncertainties and upper limits for abundances.

The Gaia-ESO survey: Mapping the shape and evolution of the radial abundance gradients with open clusters

Context. The spatial distribution of elemental abundances and their time evolution are among the major constraints to disentangling the scenarios of formation and evolution of the Galaxy. Aims. In this paper 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. Methods. 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 (RGC) 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. Results. The [Fe/H] gradient has a slope of −0.054 dex kpc−1. It can be better approximated with a two-slope shape, steeper for RGC ≤ 11.2 kpc and flatter in the outer regions. We saw different behaviours for elements belonging to different channels. For the time evolution of the gradient, we found that the youngest clusters (age < 1 Gyr) in the inner disc have lower metallicity than their older counterparts and that they outline a flatter gradient. We considered some possible explanations, including the effects of gas inflow and migration. We suggest that the most likely one may be related to a bias introduced by the standard spectroscopic analysis producing lower metallicities in the analysis of low-gravity stars. Conclusions. To delineate the shape of the ‘true’ gradient, we should most likely limit our analysis to stars with low surface gravity log g > 2.5 and microturbulent parameter ξ < 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 over the last 3 Gyr. We found a secondary role of cluster migration in shaping the gradient, with a more prominent role of migration for the oldest clusters.

The Gaia-ESO Survey: Lithium measurements and new curves of growth

Context. 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. Aims. In this paper, we describe the method adopted in the final data release to derive lithium equivalent widths (EWs) and abundances. Methods. 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. Results. 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. Conclusions. 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.

The Gaia-ESO Survey: Chemical tagging in the thin disk

Context. 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. Aims. 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 some kind of guidance on different strategies of clustering analysis in the elemental abundance space. Methods. 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 (e.g., choice of the algorithm, data preprocessing techniques, metric, space of data clustering), 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. Results. 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 (i.e., escapers, binaries). Conclusions. Our results indicate that there are chances of recovering open clusters’ members via clustering analysis in the elemental abundance space, albeit in a data set that has a very high fraction of cluster members compared to an average field star sample. Presumably, the performances of chemical tagging will further increase with higher quality data and more sophisticated clustering algorithms, which will likely became available in the near future.

The Gaia-ESO Public Spectroscopic Survey: Implementation, data products, open cluster survey, science, and legacy

Context. 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, 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. Aims. 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 reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. Methods. We made use of the information recorded and archived in the observing blocks; during the observing runs; in a number of relevant documents; in the spectra and master catalogue of spectra; in the parameters delivered by the analysis nodes and the working groups; in the final catalogue; and in the science papers. Based on these sources, we critically analyse and discuss the output and products of the Survey, including science highlights. We also determined the average metallicities of the open clusters observed as science targets and of a sample of clusters whose spectra were retrieved from the ESO archive. Results. The Gaia-ESO Survey 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. Conclusions. The final catalogue will be released through the ESO archive in the first half of 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.

The Gaia-ESO Public Spectroscopic Survey: Motivation, implementation, GIRAFFE data processing, analysis, and final data products

Context. 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 the homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. Aims. 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 introduces the survey results. Methods. 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. Results. 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. Conclusions. 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. A companion article reviews the survey implementation, scientific highlights, the open cluster survey, and data products.

The Gaia-ESO survey: A spectroscopic study of the young open cluster NGC 3293

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-local thermodynamic equilibrium (LTE) chemical abundances for six elements (He, C, N, Ne, Mg, and Si), 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−1). Our analysis leads to about a five-fold increase in the number of cluster members with an abundance determination and it characterises the late B-star population in detail for the first time. We take advantage of the multi-epoch observations on various timescales and a temporal baseline, sometimes spanning ∼15 years, 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−1. 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 Myr through a detailed, star-to-star correction of our results for the effect of stellar rotation (e.g., gravity darkening). This is significantly older than previous estimates from turn-off fitting that fully relied on classical, non-rotating isochrones.

Assessment of [Fe/H] determinations for FGK stars in spectroscopic surveys

Context. The iron abundance [Fe/H] in the atmosphere of FGK-type stars is crucial in stellar and galactic physics. The number of stars with a measured value of [Fe/H] is increasing considerably thanks to spectroscopic surveys. However, different methodologies, inputs, and assumptions used in spectral analyses lead to different precisions in [Fe/H] and possibly to systematic differences, which need to be evaluated. It is essential to understand the characteristics of each survey to fully exploit their potential, in particular if the surveys are combined to probe a larger galactic volume and to improve statistics. Aims. The purpose of this study is to compare [Fe/H] determinations from the largest spectroscopic surveys to other catalogues taken as reference. Offsets and dispersions of the residuals are examined, as are their trends with other parameters. The investigated surveys are the latest public releases of APOGEE, GALAH, RAVE, LAMOST, SEGUE, and the Gaia-ESO Survey. Methods. We use reference samples that provide independent determinations of [Fe/H], which are compared to those from the surveys for common stars. The distribution of the residuals is assessed through simple statistics that measure the offset between two catalogues and the dispersion representative of the precision of both catalogues. When relevant, linear fits are performed. A large sample of FGK-type stars with [Fe/H] based on high-resolution, high-signal-to-noise spectroscopy was built from the PASTEL catalogue to provide a reference sample. We also use FGK members in open and globular clusters to assess the internal consistency of [Fe/H] of each survey. The agreement of median [Fe/H] values for clusters observed by different surveys is discussed. Results. All the surveys overestimate the low metallicities, and some of them also underestimate the high metallicities. They perform well in the most populated intermediate metallicity range whatever the resolution. In most cases, the typical precision that we deduce from the comparisons is in good agreement with the uncertainties quoted in the catalogues. Some exceptions to this general behaviour are discussed.

The Gaia-ESO survey: Constraining evolutionary models and ages for young low mass stars with measurements of lithium depletion and rotation

Abstract 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 ∼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 ∼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 γ Vel cluster at ∼20 Myr appears not to be correlated with rotation and is challenging to explain without a very large (>10 Myr) age spread.