Search

Recent advancements in stellar evolution modeling offer unprecedented accuracy in predicting the evolution and deaths of stars. We present new stellar evolutionary models computed with the updated PARSEC V2.0 code for a comprehensive and homogeneous grid of metallicities and initial masses. Nuclear reaction networks, mass loss prescriptions, and the treatment of elemental mixing have all been updated in PARSEC V2.0. We computed models for thirteen initial metallicities spanning $Z = 10^{-11}$ to $Z = 0.03$, with masses ranging from 2.0 M$_{\odot}$ to 2000 M$_{\odot}$, consisting of a library of over 1,100 ($\sim 2100$ tracks including pure-He models) full stellar evolution tracks. For each track, the evolution is followed from the pre-main-sequence to the most advanced early-asymptotic-giant-branch or the pre-supernova phases, depending on the stellar mass. Here, we describe the properties of the tracks and their chemical and structural evolution. We computed the final fates and the remnant masses and built the mass spectrum for each metallicity, finding that the combined black hole (BH) pair-instability mass gap spans just between 100 and 130 M$_{\odot}$. Moreover, the remnant masses provide models consistent with observed BH masses, such as those from the primaries of GW190521, Cygnus X-1, and $\textit{Gaia}$ BH3 binary systems. We computed and provided the chemical ejecta from stellar winds and explosive final fates, along with the ionizing photon rates. Our results show strong overall consistency with other tracks computed with different codes. A comparison with a large sample of observed massive stars in the Tarantula Nebula of the Large Magellanic Cloud shows that our tracks nicely reproduce the majority of stars that lie on the main sequence. All the models are publicly available and can be retrieved on the PARSEC database., Comment: Accepted for publication in A&A. 24 pages, 18 figures

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

A fraction of the extreme horizontal branch stars of globular clusters exhibit a periodic light variability that has been attributed to rotational modulation caused by surface spots. These spots are believed to be connected to inhomogeneous surface distribution of elements. However, the presence of such spots has not been tested against spectroscopic data. We analyzed the phase-resolved ESO X-shooter spectroscopy of three extreme horizontal branch stars that are members of the globular cluster $\omega$ Cen and also display periodic light variations. The aim of our study is to understand the nature of the light variability of these stars and to test whether the spots can reproduce the observed variability. Our spectroscopic analysis of these stars did not detect any phase-locked abundance variations that are able to reproduce the light variability. Instead, we revealed the phase variability of effective temperature and surface gravity. In particular, the stars show the highest temperature around the light maximum. This points to pulsations as a possible cause of the observed spectroscopic and photometric variations. However, such an interpretation is in a strong conflict with Ritter's law, which relates the pulsational period to the mean stellar density. The location of the $\omega$ Cen variable extreme horizontal branch stars in HR diagram corresponds to an extension of PG 1716 stars toward lower temperatures or blue, low-gravity, large-amplitude pulsators toward lower luminosities, albeit with much longer periods. Other models of light variability, namely, related to temperature spots, should also be tested further. The estimated masses of these stars in the range of $0.2-0.3\,M_\odot$ are too low for helium-burning objects., Comment: 12 pages, accepted for publication in Astronomy & Astrophysics

The Gaia-ESO Survey (GES) is a large public spectroscopic survey which acquired spectra for more than 100000 stars across all major components of the Milky Way. In addition to atmospheric parameters and stellar abundances that have been derived in previous papers of this series, the GES spectra allow us to detect spectroscopic binaries with one (SB1), two (SB2) or more (SBn $\ge$ 3) components. Cross-correlation functions (CCFs) have been re-computed thanks to a dozen spectral masks probing a range of effective temperatures, surface gravities and metallicities. By optimising the mask choice for a given spectrum, the new computed so-called Nacre (Narrow cross-correlation experiment) CCFs are narrower and allow to unblend more stellar components than standard masks. The Doe (Detection of Extrema) extremum-finding code then selects the individual components and provides their radial velocities. From the sample of HR10 and HR21 spectra corresponding to 37565 objects, the present study leads to the detection of 322 SB2, ten (three of them being tentative) SB3, and two tentative SB4. In particular, compared to our previous study, the Nacre CCFs allow us to multiply the number of SB2 candidates by $\approx$ 1.5. The colour-magnitude diagram reveals, as expected, the shifted location of the SB2 main sequence. A comparison between the SB identified in Gaia DR3 and the ones detected in the present work is performed and the complementarity of the two censuses is discussed. An application to mass-ratio determination is presented, and the mass-ratio distribution of the GES SB2 is discussed. When accounting for the SB2 detection rate, an SB2 frequency of $\approx$ 1.4% is derived within the present stellar sample of mainly FGK-type stars. As primary outliers identified within the GES data, SBn spectra produce a wealth of information and useful constraints for the binary population synthesis studies.

Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Galactic Disc stellar component, which is well-described by power-laws with indices ranging from $1/2$ to $2/3$. The vertically-integrated star formation rate falls from $(1.147 \pm 0.039)\times10^{-8}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at earlier times down to $(6.2 \pm 3.0) \times10^{-9}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at present times, but we find a significant peak of star formation in the 2 to 3 Gyr age bin. The total mass of stars formed per unit area over time is $118.7 \pm 6.2\, \text{M}_{\odot} \text{pc}^{-2}$, which is nearly twice the present stellar mass derived from kinematics within 1~kpc from the Galactic Plane, implying a high degree of matter recycling in successive generations of stars. The method is then modified by adopting an age-dependent correlation between the SFH across the different slices, which results in less noisy and more symmetrical results without significantly changing the previously mentioned quantities. This appears to be a promising way to improve SFH recovery in external galaxies., Comment: Accepted for publication in MNRAS

The Vera C. Rubin Observatory will undertake the Legacy Survey of Space and Time, providing an unprecedented, volume-limited catalog of star clusters in the Southern Sky, including Galactic and extragalactic star clusters. The Star Clusters subgroup of the Stars, Milky Way and Local Volume Working Group has identified key areas where Rubin Observatory will enable significant progress in star cluster research. This roadmap represents our science cases and preparation for studies of all kinds of star clusters from the Milky Way out to distances of tens of megaparsecs., Comment: accepted to PASP

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

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

We report the identification of a set of old super metal-rich dwarf stars with orbits of low eccentricity that reach a maximum height from the Galactic plane between ~0.5-1.5 kpc. We discuss their properties to understand their origins. We use data from the internal data release 6 of the Gaia-ESO Survey. We selected stars observed at high resolution with abundances of 21 species of 18 individual elements. We apply hierarchical clustering to group the stars with similar chemical abundances within the complete chemical abundance space. According to their chemical properties, this set of super metal-rich stars can be arranged into five subgroups. Four seem to follow a chemical enrichment flow, where nearly all abundances increase in lockstep with Fe. The fifth subgroup shows different chemical characteristics. All subgroups have the following features: median ages of the order of 7-9 Gyr, Solar or sub-Solar [Mg/Fe] ratios, maximum height between 0.5-1.5 kpc, low eccentricities, and a detachment from the expected metallicity gradient with guiding radius. The high metallicity of our stars is incompatible with a formation in the Solar neighbourhood. Their dynamic properties agree with theoretical expectations that these stars travelled from the inner Galaxy due to blurring and, most importantly, to churning. We suggest that most of this population's stars originated in the Milky Way's inner regions (inner disc and/or the bulge) and later migrated to the Solar neighbourhood. The region from where the stars originated had a complex chemical enrichment history, with contributions from supernovae types Ia and II and possibly asymptotic giant branch stars., Comment: 18 pages, 13 figures. Accepted for publication in A&A. Version after language proofs. The catalogue remains to be released. Abridged abstract to fit arxiv's requirements

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 describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Galaxy, to generate mock stellar samples through a population synthesis approach. The main bodies of the Magellanic Clouds are also included. A complete simulation is provided for single stars, down to the $r=27.5$ mag depth of the co-added wide-fast-deep survey images. A second simulation is provided for a fraction of the binaries, including the interacting ones, as derived with the BinaPSE module of TRILEGAL. We illustrate the main properties and numbers derived from these simulations, including: comparisons with real star counts; the expected numbers of Cepheids, long-period variables and eclipsing binaries; the crowding limits as a function of seeing and filter; the star-to-galaxy ratios, etc. Complete catalogs are accessible through the NOIRLab Astro Data Lab, while the stellar density maps are incorporated in the LSST metrics analysis framework (MAF)., Comment: Accepted for publication in the LSST focused ApJS issue

A renewed interest about the origin of \emph{r}-process elements has been stimulated by the multi-messenger observation of the gravitational event GW170817, with the detection of both gravitational waves and electromagnetic waves corresponding to the merger of two neutron stars. Such phenomenon has been proposed as one of the main sources of the \emph{r}-process. However, the origin of the \emph{r}-process elements at different metallicities is still under debate. We aim at investigating the origin of the \emph{r}-process elements in the Galactic thin disc population. From the sixth internal data release of the \emph{Gaia}-ESO we have collected a large sample of Milky Way thin- and thick-disc stars for which abundances of Eu, O, and Mg are available. The sample consists of members of 62 open clusters, located at a Galactocentric radius from $\sim 5$ kpc to $\sim 20$ kpc in the disc, in the metallicity range $[-0.5, 0.4]$ and covering an age interval from 0.1 to 7 Gy, and about 1300 Milky Way disc field stars in the metallicity range $[-1.5, 0.5]$. We compare the observations with the results of a chemical evolution model, in which we varied the nucleosynthesis sources for the three considered elements. Our main result is that Eu in the thin disc is predominantly produced by sources with short lifetimes, such as magneto-rotationally driven SNe. There is no strong evidence for additional sources at delayed times. Our findings do not imply that there cannot be a contribution from mergers of neutron stars in other environments, as in the halo or in dwarf spheroidal galaxies, but such a contribution is not needed to explain Eu abundances at thin disc metallicities., Comment: 20 pages, accepted for publication in A&A

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

The last decade has seen a revolution in our knowledge of the Galaxy thanks to the Gaia and asteroseismic space missions and the ground-based spectroscopic surveys. To complete this picture, it is necessary to map the ages of its stellar populations. During recent years, the dependence on time of abundance ratios involving slow (s) neutron-capture and $\alpha$ elements (called chemical-clocks) has been used to provide estimates of stellar ages, usually in a limited volume close to the Sun. We aim to analyse the relations of chemical clocks in the Galactic disc extending the range to R$_{\rm GC}\sim$6-20~kpc. Using the sixth internal data release of the Gaia-ESO survey, we calibrated several relations between stellar ages and abundance ratios [s/$\alpha$] using a sample of open clusters, the largest one so far used with this aim. Thanks to their wide galactocentric coverage, we investigated the radial variations of the shape of these relations, confirming their non-universality. We estimated our accuracy and precision in recovering the global ages of open clusters, and the ages of their individual members. We applied the multi-variate relations with the highest correlation coefficients to the field star population. We confirm that there is no single age-chemical clock relationship valid for the whole disc, but that there is a dependence on the galactocentric position, which is related to the radial variation of the star formation history combined with the non-monotonic dependence on metallicity of the yields of the s-process elements from low- and intermediate-mass stars. Finally, the abundance ratios [Ba/$\alpha$] are more sensitive to age than those with [Y/$\alpha$] for young disc stars, and their slopes vary less with galactocentric distance., Comment: 14 pages, 10 figures + Appendix (3 tables and 2 figures)

The Gaia-ESO Survey (GES) is a public, high-resolution spectroscopic survey with FLAMES@VLT. GES targeted in particular a large sample of open clusters (OCs) of all ages. The different kinds of OCs are useful to reach the main science goals, which are the study of the OC structure and dynamics, the use of OCs to constrain and improve stellar evolution models, and the definition of Galactic disc properties (e.g. metallicity distribution). GES is organised in 19 working groups (WGs). We describe here the work of three of them, WG4 in charge of the selection of the targets within each cluster), WG1 responsible for defining the most probable candidate members, and WG6 in charge of the preparation of the observations. As GES has been conducted before Gaia DR2, we could not make use of the Gaia astrometry to define cluster members. We made use of public and private photometry to select the stars to be observed with FLAMES. Candidate target selection was based on ground-based proper motions, radial velocities, and X-ray properties when appropriate, and it was mostly used to define the position of the clusters' evolutionary sequences in the colour-magnitude diagrams. Targets for GIRAFFE were selected near the sequences in an unbiased way. We used available information on membership only for the few UVES stars. We collected spectra for 62 confirmed OCs (a few more were taken from the ESO archive). Among them are very young clusters, where the main targets are pre-main sequence stars, clusters with very hot and massive stars currently on the main sequence, intermediate-age and old clusters where evolved stars are the main targets. The selection of targets was as inclusive and unbiased as possible and we observed a representative fraction of all possible targets, thus collecting the largest, most accurate, and most homogeneous spectroscopic data set on ever achieved. [abridged], Comment: Accepted for publication on A&A. Abstract abridged from the accepted version

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

We derive the spatially-resolved star formation history (SFH) for a $96$ deg$^2$ area across the main body of the Large Magellanic Cloud (LMC), using the near-infrared photometry from the VISTA survey of the Magellanic Clouds (VMC). The data and analyses are characterised by a great degree of homogeneity and a low sensitivity to the interstellar extinction. 756 subregions of size $0.125$ deg$^2$ -- corresponding to projected sizes of about $296\times322\,\mathrm{pc}^{2}$ in the LMC -- are analysed. The resulting SFH maps, with typical resolution of $0.2$--$0.3$ dex in logarithm of age, reveal main features in the LMC disc at different ages: the patchy star formation at recent ages, the concentration of star formation on three spiral arms and on the Bar up to ages of $\sim\!1.6$ Gyr, and the wider and smoother distribution of older populations. The period of most intense star formation occurred roughly between 4 and 0.5 Gyr ago, at rates of $\sim\!0.3\,\mathrm{M}_{\odot}\mathrm{yr}^{-1}$. We compare young and old star formation rates with the observed numbers of RR Lyrae and Cepheids. We also derive a mean extinction and mean distance for every subregion, and the plane that best describes the spatial distribution of the mean distances. Our results cover an area about 50 per cent larger than the classical SFH maps derived from optical data by Harris & Zaritsky (2009). Main differences with respect to those maps are lower star formation rates at young ages, and a main peak of star formation being identified at ages slightly younger than $1$ Gyr., Comment: Accepted for publication in MNRAS

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

Context. The TOPoS project has the goal to find and analyse Turn-Off (TO) stars of extremely low metallicity. To select the targets for spectroscopic follow-up at high spectral resolution, we have relied on low-resolution spectra from the Sloan Digital Sky Survey. Aims. In this paper we use the metallicity estimates we have obtained from our analysis of the SDSS spectra to construct the metallicity distribution function (MDF) of the Milky Way, with special emphasis on its metal-weak tail. The goal is to provide the underlying distribution out of which the TOPoS sample was extracted. Methods. We make use of SDSS photometry, Gaia photometry and distance estimates derived from the Gaia parallaxes to derive a metallicity estimate for a large sample of over 24 million TO stars. This sample is used to derive the metallicity bias of the sample for which SDSS spectra are available. Results. We determined that the spectroscopic sample is strongly biased in favour of metal-poor stars, as intended. A comparison with the unbiased photometric sample allows to correct for the selection bias. We select a sub-sample of stars with reliable parallaxes for which we combine the SDSS radial velocities with Gaia proper motions and parallaxes to compute actions and orbital parameters in the Galactic potential. This allows us to characterize the stars dynamically, and in particular to select a sub-sample that belongs to the Gaia-Sausage-Enceladus (GSE) accretion event. We are thus able to provide also the MDF of GSE. Conclusions. The metal-weak tail derived in our study is very similar to that derived in the H3 survey and in the Hamburg/ESO Survey. This allows us to average the three MDFs and provide an error bar for each metallicity bin. Inasmuch the GSE structure is representative of the progenitor galaxy that collided with the Milky Way, that galaxy appears to be strongly deficient in metal-poor stars compared to the Milky Way, suggesting that the metal-weak tail of the latter has been largely formed by accretion of low mass galaxies rather than massive galaxies, such as the GSE progenitor., Comment: Astronomy and Astrophysics - A\&A, EDP Sciences, in press

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.

We analyze the oxygen abundances of a stellar sample representative of the two major Galactic populations: the thin and thick disks. The aim is to investigate the differences between members of the Galactic disks and to contribute to the understanding on the origin of oxygen chemical enrichment in the Galaxy. The analysis is based on the [O\,{\sc i}]=6300.30\,\AA~ oxygen line in HR spectra ($R\sim$52,500) obtained from the GES Survey. By comparing the observed spectra with a theoretical dataset, computed in LTE with the SPECTRUM synthesis and ATLAS12 codes, we derive the oxygen abundances of 516 FGK dwarfs for which we have previously measured carbon abundances. Based on kinematic, chemical and dynamical considerations we identify 20 thin and 365 thick disk members. We study potential trends of both subsamples in terms of their chemistry ([O/H], [O/Fe], [O/Mg], and [C/O] versus [Fe/H] and [Mg/H]), age, and position in the Galaxy. Main results are: (a) [O/H] and [O/Fe] ratios versus [Fe/H] show systematic differences between thin and thick disk stars with enhanced O abundance of thick disk stars with respect to thin disk members and a monotonic decrement of [O/Fe] with increasing metallicity, even at metal-rich regime; (b) a smooth correlation of [O/Mg] with age in both populations, suggesting that this abundance ratio can be a good proxy of stellar ages within the Milky Way; (c) thin disk members with [Fe/H]$\simeq0$ display a [C/O] ratio smaller than the solar value, suggesting a possibly outward migration of the Sun from lower Galactocentric radii., Comment: 23 pages 9 figures Accepted on 30-October 2020; to be published in AJ

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

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

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

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

We recover the spatially resolved star formation history across the entire main body and Wing of the Small Magellanic Cloud (SMC), using fourteen deep tile images from the VISTA survey of the Magellanic Clouds (VMC), in the YJKs filters. The analysis is performed on 168 subregions of size 0.143 deg2, covering a total contiguous area of 23.57 deg2. We apply a colour-magnitude diagram (CMD) reconstruction method that returns the best-fitting star formation rate SFR(t), age--metallicity relation, distance and mean reddening, together with their confidence intervals, for each subregion. With respect to previous analyses, we use a far larger set of VMC data, updated stellar models, and fit the two available CMDs (Y-Ks versus Ks and J-Ks versus Ks) independently. The results allow us to derive a more complete and more reliable picture of how the mean distances, extinction values, star formation rate, and metallicities vary across the SMC, and provide a better description of the populations that form its Bar and Wing. We conclude that the SMC has formed a total mass of (5.31+-0.05)x10^8 Msun in stars over its lifetime. About two thirds of this mass is expected to be still locked in stars and stellar remnants. 50 per cent of the mass was formed prior to an age of 6.3 Gyr, and 80 per cent was formed between 8 and 3.5 Gyr ago. We also illustrate the likely distribution of stellar ages and metallicities in different parts of the CMD, to aid the interpretation of data from future astrometric and spectroscopic surveys of the SMC., Comment: To appear in MNRAS

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

We present the first spatially resolved map of stellar proper motions within the central ($\sim$3.1 $\times$ 2.4 kpc) regions of the Small Magellanic Cloud (SMC). The data used for this study encompasses four tiles from the ongoing near-infrared VISTA survey of the Magellanic Clouds system and covers a total contiguous area on the sky of $\sim$6.81 deg$^2$. Proper motions have been calculated independently in two dimensions from the spatial offsets in the $K_s$ filter over time baselines between 22 and 27 months. The reflex motions of $\sim$33~000 background galaxies are used to calibrate the stellar motions to an absolute scale. The resulting catalog is composed of more than 690 000 stars which have been selected based on their position in the $(J-K_s, K_s)$ color-magnitude diagram. For the median absolute proper motion of the SMC, we find ($\mu_{\alpha}\mathrm{cos}(\delta)$, $\mu_{\delta}$) = (1.087 $\pm$ 0.192 (sys.) $\pm$ 0.003 (stat.), $-$1.187 $\pm$ 0.008 (sys.) $\pm$ 0.003 (stat.)) mas yr$^{-1}$, consistent with previous studies. Mapping the proper motions as a function of position within the SMC reveals a non uniform velocity pattern indicative of a tidal feature behind the main body of the SMC and a flow of stars in the South-East moving predominantly along the line-of-sight., Comment: 4 pages, 2 figures, accepted for publication in A&A

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

Extremely metal-poor stars provide us with indirect information on the first generations of massive stars. The TOPoS survey has been designed to increase the census of these stars and to provide a chemical inventory that is as detailed as possible. Seven of the most iron-poor stars have been observed with the UVES spectrograph at the ESO VLT Kueyen 8.2m telescope to refine their chemical composition. We analysed the spectra based on 1D LTE model atmospheres, but also used 3D hydrodynamical simulations of stellar atmospheres. We measured carbon in six of the seven stars: all are carbon-enhanced and belong to the low-carbon band, defined in the TOPoS II paper. We measured lithium (A(Li)=1.9) in the most iron-poor star (SDSS J1035+0641, [Fe/H] < -5.2 ). We were also able to measure Li in three stars at [Fe/H]~ -4.0, two of which lie on the Spite plateau. We confirm that SDSS J1349+1407 is extremely rich in Mg, but not in Ca. It is also very rich in Na. Several of our stars are characterised by low alpha-to-iron ratios. The lack of high-carbon band stars at low metallicity can be understood in terms of evolutionary timescales of binary systems. The detection of Li in SDSS J1035+0641 places a strong constraint on theories that aim at solving the cosmological lithium problem. The Li abundance of the two warmer stars at [Fe/H]~ -4.0 places them on the Spite plateau, while the third, cooler star, lies below. We argue that this suggests that the temperature at which Li depletion begins increases with decreasing [Fe/H]. SDSS J1349+1407 may belong to a class of Mg-rich EMP stars. We cannot assess if there is a scatter in alpha-to-iron ratios among the EMP stars or if there are several discrete populations. However, the existence of stars with low alpha-to-iron ratios is supported by our observations., Comment: A&A, in press

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

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