Search

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

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

IC 4665 is one of only a dozen young open clusters with a ``lithium depletion boundary" (LDB) age. Using an astrometrically and spectroscopically filtered sample of cluster members, we show that both the positions of its low mass stars in Gaia absolute colour-magnitude diagrams and the lithium depletion seen among its K- and early M-stars are discordant with the reported LDB age of (32 +4/-5) Myr. Re-analysis of archival spectra suggests that the LDB of IC 4665 has not been detected and that the published LDB age should be interpreted as a lower limit. Empirical comparisons with similar datasets from other young clusters with better-established LDB ages indicate that IC 4665 is bracketed in age by the clusters IC 2602 and IC 2391 at (55 +/- 3) Myr., Comment: For publication in MNRAS, accepted version

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

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

The Gaia-ESO Survey (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

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

The kinematics of stars in OB associations can provide insights into their formation, dynamical evolution, and eventual fate. The low-mass stellar content of OB associations are sufficiently numerous as to provide a detailed sampling of their kinematic properties, however spectroscopy is required to confirm the youth of individual stars and to get 3D kinematics. In this paper we present and analyse results from a large spectroscopic survey of Vela OB2 conducted using 2dF/HERMES on the AAT. This spectroscopy is used to confirm the youth of candidate young stars and determine radial velocities, which are combined with proper motions and parallaxes from Gaia to measure 3-dimensional positions and velocities. We identify multiple separate kinematic groups in the region, for which we measure velocity dispersions and infer their virial states. We measure expansion rates for all these groups and find strong evidence for anisotropic expansion in the Vela OB2 association of at least 11$\sigma$ significance in all three dimensions, as well as some evidence for expansion in the $\gamma$ Vel and P Puppis clusters. We trace back the motions of these groups into the past and find that the open cluster NGC 2547 is an interloper in the Vela OB2 region and actually formed $>$100 pc away from the association. We conclude that Vela OB2 must have formed with considerable spatial and kinematic substructure over a timescale of $\sim$10 Myr, with clear temporal substructure within the association, but no clear evidence for an age gradient., Comment: 19 pages, 13 figures, 4 tables, submitted to MNRAS

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

With an unprecedented astrometric and photometric data precision, Gaia EDR3 gives us, for the first time, the opportunity to systematically detect and map in the optical bands, the low mass populations of the star forming regions (SFRs) in the Milky Way. We provide a catalogue of the Gaia EDR3 data (photometry, proper motions and parallaxes) of the young stellar objects (YSOs) identified in the Galactic Plane (|b|<30 deg) within about 1.5 kpc. The catalogue of the SFRs to which they belong is also provided to study the properties of the very young clusters and put them in the context of the Galaxy structure. We applied the machine learning unsupervised clustering algorithm DBSCAN on a sample of Gaia EDR3 data photometrically selected on the region where very young stars (t<10 Myr) are expected to be found, with the aim to identify co-moving and spatially consistent stellar clusters. A subsample of 52 clusters, selected among the 7323 found with DBSCAN, has been used as template data set, to identify very young clusters from the pattern of the observed color-absolute magnitude diagrams through a pattern match process. We find 124440 candidate YSOs clustered in 354 SFRs and stellar clusters younger than 10 Myr and within about 1.5 Kpc. In addition, 65863 low mass members of 322 stellar clusters located within about 500 pc and with ages 10 Myr

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

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 exploit the unique characteristics of a sample of open clusters (OCs) and field stars for which high-precision 7Li abundances and stellar parameters are homogeneously derived by the Gaia-ESO Survey (GES). We derive possibly undepleted 7Li abundances for 26 OCs and star forming regions with ages from young to old spanning a large range of Galactocentric distances, which allows us to reconstruct the local late Galactic evolution of lithium as well as its current abundance gradient along the disc. Field stars are added to look further back in time and to constrain 7Li evolution in other Galactic components. The data are then compared to theoretical tracks from chemical evolution models that implement different 7Li forges. We find that the upper envelope of the 7Li abundances measured in field stars of nearly solar metallicities traces very well the level of lithium enrichment attained by the ISM as inferred from observations of cluster stars. We confirm previous findings that the abundance of 7Li in the solar neighbourhood does not decrease at supersolar metallicity. The comparison of the data with the chemical evolution model predictions favours a scenario in which the majority of the 7Li abundance in meteorites comes from novae. Current data also seem to suggest that the nova rate flattens out at later times. This requirement might have implications for the masses of the white dwarf nova progenitors and deserves further investigation. Neutrino-induced reactions taking place in core-collapse supernovae also produce some fresh lithium. This likely makes a negligible contribution to the meteoritic abundance, but could be responsible for a mild increase of the 7Li abundance in the ISM of low-metallicity systems that would counterbalance the astration processes., Comment: 17 pages, 10 figures, 3 tables (full tables 1 and 3 only available at the CDS), accepted for publication in Astronomy and Astrophysics. Version 2 includes corrections from language editor

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.

New fibre spectroscopy and radial velocities from the WIYN telescope are used to measure photospheric lithium in 242 high-probability, zero-age-main-sequence (ZAMS) F- to K-type members of the rich cluster M35. Combining these with published rotation periods, the connection between lithium depletion and rotation is studied in unprecedented detail. At $T_{\rm eff}<5500$ K there is a strong relationship between faster rotation and less Li depletion, although with a dispersion larger than measurement uncertainties. Components of photometrically identified binary systems follow the same relationship. A correlation is also established between faster rotation rate (or smaller Rossby number), decreased Li depletion and larger stellar radius at a given $T_{\rm eff}$. These results support models where starspots and interior magnetic fields lead to inflated radii and reduced Li depletion during the pre main sequence (PMS) phase for the fastest rotators. However, the data are also consistent with the idea that all stars suffered lower levels of Li depletion than predicted by standard PMS models, perhaps because of deficiencies in those models or because saturated levels of magnetic activity suppress Li depletion equally in PMS stars of similar $T_{\rm eff}$ regardless of rotation rate, and that slower rotators subsequently experience more mixing and post-PMS Li depletion., Comment: Accepted version, for publication in MNRAS

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

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

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

The kinematics of low-mass stars in nearby OB associations can provide clues about their origins and evolution. Combining the precise positions, proper motions and parallaxes given in the second Gaia Data Release with radial velocity measurements obtained with the Hermes spectrograph at the Anglo-Australian Telescope, we have an opportunity to study in detail the kinematics of low-mass stars belonging to the nearby $\gamma$ Vel cluster and the Vela OB2 association it is projected against. The presence of lithium is used to confirm the youth of our targets. We separate our sample into the cluster and association populations based on the membership probabilities of Jeffries et al. 2014, their parallaxes, and kinematics. We find strong evidence for expansion in the OB association population with at least 4$\sigma$ significance along all three axes, though the expansion is notably anisotropic. We discuss these results in the context of cluster and association dispersal theories., Comment: 9 pages, 8 figures, 1 table, accepted for publication in MNRAS

In the last three decades several hundred nearby members of young stellar moving groups (MGs) have been identified, but there has been less systematic effort to quantify or characterise young stars that do not belong to previously identified MGs. Using a kinematically unbiased sample of 225 lithium-rich stars within 100 pc, we find that only $50 \pm 10$ per cent of young ($\lesssim 125$ Myr), low-mass ($0.5

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

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

In the era of large high-resolution spectroscopic surveys, high-quality spectra can contribute to our understanding of the Galactic chemical evolution, providing chemical abundances belonging to the different nucleosynthesis channels, and also providing constraints to stellar age. Some abundance ratios have been proven to be excellent indicators of stellar ages. We aim at providing an empirical relationship between stellar ages and [C/N] using, as calibrators, open star clusters observed by both the Gaia-ESO and APOGEE surveys. We use stellar parameters and abundances from the Gaia-ESO and APOGEE of the Galactic field and open cluster stars. Ages of star clusters are retrieved from the literature sources and validated using a common set of isochrones. We use the same isochrones to determine, for each age and metallicity, the surface gravity at which the first dredge-up and red giant branch bump occur. We study the effect of extra-mixing processes in our sample of giant stars, and we derive the mean [C/N] in evolved stars, including only stars without evidence of extra-mixing. Combining the Gaia-ESO and APOGEE samples of open clusters, we derive a linear relationship between [C/N] and logarithmic cluster ages. We apply our relationship to selected giant field stars in both Gaia-ESO and APOGEE. We find an age separation between thin and thick disc stars and age trends within their populations, with an increasing age towards lower metallicity populations. With such empirical relationship, we are able to provide an age estimate for giant stars in which C and N abundances are measured. Isochrone fitting is less sensitive for giant than dwarf stars at the turn off. The present method can be thus considered as an additional tool to give an independent estimate of the age of giant stars, with uncertainties in their ages comparable to those obtained using isochrone fitting for dwarf stars., Comment: 28 pages, 24 figures. Accepted by A&A

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

The Vela-Puppis region is known to host the Vela OB2 association as well as several young clusters featuring OB and pre-main sequence stars. Several spatial and kinematic subgroups have been identified in the recent years. By grouping stars based on their positions and velocity we can address the quetion of the dynamical history of the region and the mechanisms that drove stellar formation. The Gaia DR2 astrometry and photometry enables us to characterise the 3D spatial and 3D kinematic distribution of young stars and to estimate the ages of the identified components. We use an unsupervised classification method to group stars based on their proper motions and parallax. We perform a study of the expansion rates of the different identified groups from 3D velocities, and from corrected tangential velocities. We make use of theoretical isochrones to estimate ages. The young stars can be separated into seven main groups of different ages and kinematical distribution. All groups are found to be expanding, although the expansion is mostly not isotropic. The size of the region, the age substructure, and the anistropic expansion rates are compatible with a prolonged period of star formation in a turbulent molecular cloud, and that the dispersion of the stars cannot be explained by gas expulsion alone., Comment: accepted A&A

Rotation periods from Kepler K2 are combined with projected rotation velocities from the WIYN 3.5-m telescope, to determine projected radii for fast-rotating, low-mass ($0.15 \leq M/M_{\odot} \leq 0.6$) members of the Praesepe cluster. A maximum likelihood analysis that accounts for observational uncertainties, binarity and censored data, yields marginal evidence for radius inflation -- the average radius of these stars is $6\pm4$ per cent larger at a given luminosity than predicted by commonly-used evolutionary models. This over-radius is smaller (at 2-sigma confidence) than was found for similar stars in the younger Pleiades using a similar analysis; any decline appears due to changes occurring in higher mass ($>0.25 M_{\odot}$) stars. Models incorporating magnetic inhibition of convection predict an over-radius, but do not reproduce this mass dependence unless super-equipartition surface magnetic fields are present at lower masses. Models incorporating flux-blocking by starspots can explain the mass dependence but there is no evidence that spot coverage diminishes between the Pleiades and Praesepe samples to accompany the decline in over-radius. The fastest rotating stars in both Praesepe and the Pleiades are significantly smaller than the slowest rotators for which a projected radius can be measured. This may be a selection effect caused by more efficient angular momentum loss in larger stars leading to their progressive exclusion from the analysed samples. Our analyses assume random spin-axis orientations; any alignment in Praesepe, as suggested by Kovacs (2018), is strongly disfavoured by the broad distribution of projected radii., Comment: Accepted version for publication in MNRAS, 14 pages

We present a study of the kinematics and structure of the Cep OB3b cluster based on new spectra obtained with the Hectoschelle spectrograph on the MMT and data from Spitzer, Chandr}, and Gaia. At a distance of 819+/-16 pc, Cep OB3b is one of the closest examples of a young (~3 - 5 Myr), large (~3000 total members) cluster at the late stages of gas dispersal. The cluster is broken into two sub-clusters surrounded by a lower density halo. We fit the empirical density law of King (1962) to each sub-cluster to constrain their sizes and structure. The richer eastern sub-cluster has circular symmetry, a modest central density, and lacks molecular gas toward its core suggesting it has undergone expansion due to gas dispersal. In contrast, the western sub-cluster deviates from circular symmetry, has a smaller core size, and contains significant molecular gas near its core, suggesting that it is in an earlier phase of gas dispersal. We present posterior probability distributions for the velocity dispersions from the Hectoschelle spectra. The east will continue to expand and likely form a bound cluster with ~35% of stars remaining. The west is undergoing slower gas dispersal and will potentially form a bound cluster with ~75% of stars remaining. If the halo dissipates, this will leave two independent clusters with ~300 members; proper motions suggest that the two sub-clusters are not bound to each other., Comment: 44 pages, 15 figures, 9 tables, accepted in ApJ

The first Gaia Data Release presents an opportunity to characterise the low-mass population of OB associations, providing larger statistical samples and better constraints on the formation and evolution of clusters and associations. Using previously known low mass members in a small region of Vela OB2 we have designed selection criteria that combine Gaia and 2MASS photometry, independently of any astrometric information, to identify low-mass pre-main-sequence (PMS) stars over the wider association area. Our method picks out the known clusters of young stars around $\gamma^2$ Velorum and NGC-2547, but also identifies other over-densities that may represent previously unknown clusters. There are clear differences in the spatial distributions of the low-mass and the high-mass OB populations, suggesting either that the structure and dynamics of these populations has evolved separately or that the initial mass function can vary considerably on small scales within a single association., Comment: 6 pages, 4 figures, accepted for publication in MNRAS

Gaia-ESO Survey observations of the young Gamma Velorum cluster led to the discovery of two kinematically-distinct populations, Gamma Vel A and B, respectively, with population B extended over several square degrees in the Vela OB2 association. Using the Gaia DR2 data for a sample of high-probability cluster members, we find that the two populations differ not only kinematically, but are also located at different distances along the line of sight, with the main cluster Gamma Vel A being closer. A combined fit of the two populations yields $\varpi_A = 2.895 \pm 0.008$ mas and $\varpi_B = 2.608 \pm 0.017$ mas, with intrinsic dispersions of $0.038 \pm 0.011$ mas and $0.091 \pm 0.016$ mas, respectively. This translates into distances of $345.4^{+1.0+12.4}_{-1.0-11.5}\,$ pc and $383.4^{+2.5+15.3}_{-2.5-14.2}\,$ pc, respectively, showing that Gamma Vel A is closer than Gamma Vel B by $\sim$38 pc. We find that the two clusters are nearly coeval, and that Gamma Vel B is expanding. We suggest that Gamma Vel A and B are two independent clusters located along the same line of sight., Comment: 5 pages, 3 figures. Revised version after referee's comments, accepted for publication in Astronomy & Astrophysics

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

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

Rotation periods obtained with the Kepler satellite have been combined with precise measurements of projected rotation velocity from the WIYN 3.5-m telescope to determine the distribution of projected radii for several hundred low-mass ($0.1 \leq M/M_{\odot} \leq 0.8$), fast-rotating members of the Pleiades cluster. A maximum likelihood modelling technique, that takes account of observational uncertainties, selection effects and censored data, and considers the effects of differential rotation and unresolved binarity, has been used to find that the average radius of these stars is $14 \pm 2$ per cent larger at a given luminosity than predicted by the evolutionary models of Dotter et al. (2008) and Baraffe et al. (2015). The same models are a reasonable match to the interferometric radii of older, magnetically inactive field M-dwarfs, suggesting that the over-radius may be associated with the young, magnetically active nature of the Pleiades objects. No evidence is found for any change in this over-radius above and below the boundary marking the transition to full convection. Published evolutionary models that incorporate either the effects of magnetic inhibition of convection or the blocking of flux by dark starspots do not individually explain the radius inflation, but a combination of the two effects might. The distribution of projected radii is consistent with the adopted hypothesis of a random spatial orientation of spin axes; strong alignments of the spin vectors into cones with an opening semi-angle $<30^{\circ}$ can be ruled out. Any plausible but weaker alignment would increase the inferred over-radius., Comment: Accepted version, to appear in MNRAS

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

We 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

I review progress towards understanding the time-scales of star and cluster formation and of the absolute ages of young stars. I focus in particular on the areas in which Francesco Palla made highly significant contributions - interpretation of the Hertzsprung-Russell diagrams of young clusters and the role of photospheric lithium as an age diagnostic., Comment: To appear in "Francesco's Legacy: Star Formation in Space and Time", Memorie della SAIt, in press. Eds. R. Cesaroni, E. Corbelli and D. Galli. 5pp

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

Recent results using radial-velocity measurements from the Gaia-ESO Survey have led to the discovery of multiple kinematic populations across the Vela OB2 association. We present here a proper-motion study of the same region. Our aim is to test whether the radial-velocity populations have a counterpart in proper-motion space, and if so, how the two sets of kinematical data complement each other. This work is based on parallaxes and proper motions from the TGAS catalog, as part of Gaia DR1. Two distinct proper-motion populations are found dispersed across ~5 degrees (or ~30 pc at their likely distances). Their detailed correspondence to the radial-velocity populations could not be tested, because of the paucity of common objects. However, compelling indications are found that one of the new proper-motion populations consists mostly of members of the young cluster NGC 2547, and the other is related to the gamma Vel cluster. Constraints on the age of the two populations, both of which appear to be only 10-35 Myr old, and their possible mutual interactions within the last 1.5 Myrs, are discussed., Comment: 4 pages, 6 figures, accepted for publication on A&A Letters

We present the first extensive spectroscopic study of the global population in star clusters Trumpler~16, Trumpler~14 and Collinder~232 in the Carina Nebula, using data from the Gaia-ESO Survey, down to solar-mass stars. In addition to the standard homogeneous Survey data reduction, a special processing was applied here because of the bright nebulosity surrounding Carina stars. We find about four hundred good candidate members ranging from OB types down to slightly sub-solar masses. About one-hundred heavily-reddened early-type Carina members found here were previously unrecognized or poorly classified, including two candidate O stars and several candidate Herbig Ae/Be stars. Their large brightness makes them useful tracers of the obscured Carina population. The spectroscopically-derived temperatures for nearly 300 low-mass members allows the inference of individual extinction values, and the study of the relative placement of stars along the line of sight. We find a complex spatial structure, with definite clustering of low-mass members around the most massive stars, and spatially-variable extinction. By combining the new data with existing X-ray data we obtain a more complete picture of the three-dimensional spatial structure of the Carina clusters, and of their connection to bright and dark nebulosity, and UV sources. The identification of tens of background giants enables us also to determine the total optical depth of the Carina nebula along many sightlines. We are also able to put constraints on the star-formation history of the region, with Trumpler~14 stars found to be systematically younger than stars in other sub-clusters. We find a large percentage of fast-rotating stars among Carina solar-mass members, which provide new constraints on the rotational evolution of pre-main-sequence stars in this mass range., Comment: 35 pages, 9 Tables, 23 figures. Accepted for journal Astronomy and Astrophysics

The radial metallicity distribution in the Galactic thin disc represents a crucial constraint for modelling disc formation and evolution. Open clusters allow us to derive both the radial metallicity distribution and its evolution over time. In this paper we perform the first investigation of the present-day radial metallicity distribution based on [Fe/H] determinations in late type members of pre-main-sequence clusters. Because of their youth, these clusters are therefore essential for tracing the current inter-stellar medium metallicity. We used the products of the Gaia-ESO Survey analysis of 12 young regions (age<100 Myr), covering Galactocentric distances from 6.67 to 8.70 kpc. For the first time, we derived the metal content of star forming regions farther than 500 pc from the Sun. Median metallicities were determined through samples of reliable cluster members. For ten clusters the membership analysis is discussed in the present paper, while for other two clusters (Chamaeleon I and Gamma Velorum) we adopted the members identified in our previous works. All the pre-main-sequence clusters considered in this paper have close-to-solar or slightly sub-solar metallicities. The radial metallicity distribution traced by these clusters is almost flat, with the innermost star forming regions having [Fe/H] values that are 0.10-0.15 dex lower than the majority of the older clusters located at similar Galactocentric radii. This homogeneous study of the present-day radial metallicity distribution in the Galactic thin disc favours models that predict a flattening of the radial gradient over time. On the other hand, the decrease of the average [Fe/H] at young ages is not easily explained by the models. Our results reveal a complex interplay of several processes (e.g. star formation activity, initial mass function, supernova yields, gas flows) that controlled the recent evolution of the Milky Way., Comment: 13 pages, 9 figures, accepted for publication

The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured., Comment: Accepted for publication on Astronomy and Astrophysics

Context: Trumpler 23 is a moderately populated, intermediate-age open cluster within the solar circle at a Rgc ~6 kpc. It is in a crowded field very close to the Galactic plane and the color-magnitude diagram shows significant field contamination and possible differential reddening; it is a relatively understudied cluster for these reasons, but its location makes it a key object for determining Galactic abundance distributions. Aims: New data from the Gaia-ESO Survey enable the first ever radial velocity and spectroscopic metallicity measurements for this cluster. We aim to use velocities to isolate cluster members, providing more leverage for determining cluster parameters. Methods: Gaia-ESO Survey data for 167 potential members have yielded radial velocity measurements, which were used to determine the systemic velocity of the cluster and membership of individual stars. Atmospheric parameters were also used as a check on membership when available. Literature photometry was used to re-determine cluster parameters based on radial velocity member stars only; theoretical isochrones are fit in the V, V-I diagram. Cluster abundance measurements of ten radial-velocity member stars with high-resolution spectroscopy are presented for 24 elements. These abundances have been compared to local disk stars, and where possible placed within the context of literature gradient studies. Results: We find Trumpler 23 to have an age of 0.80 +/- 0.10 Gyr, significant differential reddening with an estimated mean cluster E(V-I) of 1.02 +0.14/-0.09, and an apparent distance modulus of 14.15 +/- 0.20. We find an average cluster metallicity of [Fe/H] = 0.14 +/- 0.03 dex, a solar [alpha/Fe] abundance, and notably subsolar [s-process/Fe] abundances., Comment: 18 pages, 10 figures, accepted for publication in Astronomy & Astrophysics

We show that non-magnetic models for the evolution of pre-main-sequence (PMS) stars *cannot* simultaneously describe the colour-magnitude diagram (CMD) and the pattern of lithium depletion seen in the cluster of young, low-mass stars surrounding $\gamma^2$ Velorum. The age of 7.5+/-1 Myr inferred from the CMD is much younger than that implied by the strong Li depletion seen in the cluster M-dwarfs and the Li depletion occurs at much redder colours than predicted. The epoch at which a star of a given mass depletes its Li and the surface temperature of that star are both dependent on its radius. We demonstrate that if the low-mass stars have radii ~10 per cent larger at a given mass and age, then both the CMD and Li depletion pattern of the Gamma Vel cluster are explained at a common age of 18-21 Myr. This radius inflation could be produced by some combination of magnetic suppression of convection and extensive cool starspots. Models that incorporate radius inflation suggest that PMS stars similar to those in the Gamma Vel cluster, in the range 0.230 per cent) than inferred from conventional, non-magnetic models in the Hertzsprung-Russell diagram. Systematic changes of this size may be of great importance in understanding the evolution of young stars, disc lifetimes and the formation of planetary systems., Comment: 10 pages, accepted for publication in MNRAS

The discovery of lithium-rich giants contradicts expectations from canonical stellar evolution. Here we report on the serendipitous discovery of 20 Li-rich giants observed during the Gaia-ESO Survey, which includes the first nine Li-rich giant stars known towards the CoRoT fields. Most of our Li-rich giants have near-solar metallicities, and stellar parameters consistent with being before the luminosity bump. This is difficult to reconcile with deep mixing models proposed to explain lithium enrichment, because these models can only operate at later evolutionary stages: at or past the luminosity bump. In an effort to shed light on the Li-rich phenomenon, we highlight recent evidence of the tidal destruction of close-in hot Jupiters at the sub-giant phase. We note that when coupled with models of planet accretion, the observed destruction of hot Jupiters actually predicts the existence of Li-rich giant stars, and suggests Li-rich stars should be found early on the giant branch and occur more frequently with increasing metallicity. A comprehensive review of all known Li-rich giant stars reveals that this scenario is consistent with the data. However more evolved or metal-poor stars are less likely to host close-in giant planets, implying that their Li-rich origin requires an alternative explanation, likely related to mixing scenarios rather than external phenomena., Comment: Accepted to MNRAS

Understanding the properties of young open clusters, such as the Initial Mass Function (IMF), star formation history and dynamic evolution, is crucial to obtain reliable theoretical predictions of the mechanisms involved in the star formation process. We want to obtain a list, as complete as possible, of confirmed members of the young open cluster Gamma Velorum, with the aim of deriving general cluster properties such as the IMF. We used all available spectroscopic membership indicators within the Gaia-ESO public archive together with literature photometry and X-ray data and, for each method, we derived the most complete list of candidate cluster members. Then, we considered photometry, gravity and radial velocities as necessary conditions to select a subsample of candidates whose membership was confirmed by using the lithium and H$\alpha$ lines and X-rays as youth indicators. We found 242 confirmed and 4 possible cluster members for which we derived masses using very recent stellar evolutionary models. The cluster IMF in the mass range investigated in this study shows a slope of $\alpha=2.6\pm0.5$ for $0.5

The Gaia ESO Public Spectroscopic Survey (GES) is providing the astronomical community with high-precision measurements of many stellar parameters including radial velocities (RVs) of stars belonging to several young clusters and star-forming regions. One of the main goals of the young cluster observations is to study of their dynamical evolution and provide insight into their future, revealing if they will eventually disperse to populate the field, rather than evolve into bound open clusters. In this paper we report the analysis of the dynamical state of L1688 in the $\rho$~Ophiuchi molecular cloud using the dataset provided by the GES consortium. We performed the membership selection of the more than 300 objects observed. Using the presence of the lithium absorption and the location in the Hertzspung-Russell diagram, we identify 45 already known members and two new association members. We provide accurate RVs for all 47 confirmed members.A dynamical analysis, after accounting for unresolved binaries and errors, shows that the stellar surface population of L1688 has a velocity dispersion $\sigma \sim$1.14$\pm$0.35 km s$^{-1}$ that is consistent with being in virial equilibrium and is bound with a $\sim$80% probability. We also find a velocity gradient in the stellar surface population of $\sim$1.0 km s$^{-1}$pc$^{-1}$ in the northwest/southeast direction, which is consistent with that found for the pre-stellar dense cores, and we discuss the possibility of sequential and triggered star formation in L1688., Comment: 12 pages, 7 figures, 1 online material table

Rapidly rotating, low-mass members of eclipsing binary systems have measured radii significantly larger than predicted by standard models. It has been proposed that magnetic activity is responsible for radius inflation. By estimating the radii of low-mass stars in three young clusters (NGC 2264, NGC 2547, NGC 2516, with ages of 5, 35 and 140 Myr respectively), we aim to establish whether similar radius inflation is seen in single, magnetically active stars. We use radial velocities from the Gaia-ESO Survey (GES) and published photometry to establish cluster membership and combine GES measurements of vsini with published rotation periods to estimate average radii for groups of fast-rotating cluster members as a function of their luminosity and age. The average radii are compared with the predictions of both standard evolutionary models and variants that include magnetic inhibition of convection and starspots. At a given luminosity, the stellar radii in NGC 2516 and NGC 2547 are larger than predicted by standard evolutionary models at the ages of these clusters. The discrepancy is least pronounced and not significant ~10 percent) in ZAMS stars with radiative cores, but more significant in lower-mass, fully convective pre main-sequence cluster members, reaching 30+/-10 percent. The uncertain age and distance of NGC 2264 preclude a reliable determination of any discrepancy for its members. The median radii we have estimated for low-mass fully convective stars in the older clusters are inconsistent (at the 2-3 sigma level) with non-magnetic evolutionary models and more consistent with models that incorporate the effects of magnetic fields or dark starspots. The available models suggest this requires either surface magnetic fields exceeding 2.5 kG, spots that block about 30 per cent of the photospheric flux, or a more moderate combination of both. [Abridged], Comment: 12 pages. Accepted for publication in Astronomy & Astrophysics