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.

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

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.

The Gaia-ESO Survey is a public spectroscopic survey that targeted ≳105 stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its final public release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O (Teff ~ 30 000–52 000 K) all the way to K-M (≳3500 K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different types. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products, such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, for example, activity and accretion indicators in young stars and membership probability in star clusters for more than 114 000 stars. The spectral analysis is coordinated by a number of working groups (WGs) within the survey, each specialised in one or more of the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all survey results are placed on a common scale in order to arrive at a single set of recommended results for use by all survey collaborators. We also present some general quality and consistency checks performed on the entirety of the survey results.

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

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.

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.

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.

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

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]

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.

ABSTRACT Astrometry and photometry from Gaia and spectroscopic data from the Gaia-ESO Survey (GES) are used to identify the lithium depletion boundary (LDB) in the young cluster NGC 2232. A specialized spectral line analysis procedure was used to recover the signature of undepleted lithium in very low luminosity cluster members. An age of 38 ± 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 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-eff < 5500K 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-eff. These results support models where star-spots 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-eff regardless of rotation rate, and that slower rotators subsequently experience more mixing and post-PMS Li depletion.

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.

Context. 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. Aims. 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. Methods. We trace the upper envelope of Li abundance versus metallicity evolution using 18 clusters and considering members that should not have suffered any Li depletion. Results. 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] = ∼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.

Context.Multiple stellar systems play a fundamental role in the formation and evolution of stellar populations in galaxies. Recent and ongoing large ground-based multi-object spectroscopic surveys significantly increase the sample of spectroscopic binaries (SBs) allowing analyses of their statistical properties.Aims.We investigate the repeated spectral observations of theGaia-ESO Survey internal data release 5 (GES iDR5) to identify and characterise SBs with one visible component (SB1s) in fields covering mainly the discs, the bulge, the CoRot fields, and some stellar clusters and associations.Methods.A statisticalχ2-test is performed on spectra of the iDR5 subsample of approximately 43 500 stars characterised by at least two observations and a signal-to-noise ratio larger than three. In the GES iDR5, most stars have four observations generally split into two epochs. A careful estimation of the radial velocity (RV) uncertainties is performed. Our sample of RV variables is cleaned from contamination by pulsation- and/or convection-induced variables usingGaiaDR2 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 relation to effective temperature and metallicity.Results.We find 641 (resp., 803) FGK SB1 candidates at the 5σ(resp., 3σ) level. The maximum RV differences range from 2.2 km s−1at the 5σconfidence level (1.6 km s−1at 3σ) to 133 km s−1(in both cases). Among them a quarter of the primaries are giant stars and can be located as far as 10 kpc. The orbital-period distribution is estimated from the RV standard-deviation distribution and reveals that the detected SB1s probe binaries with logP[d] ⪅ 4. We show that SB1s with dwarf primaries tend to have shorter orbital periods than SB1s with giant primaries. This is consistent with binary interactions removing shorter period systems as the primary ascends the red giant branch. For two systems, tentative orbital solutions with periods of 4 and 6 d are provided. After correcting for detection efficiency, selection biases, and the present-day mass function, we estimate the global GES SB1 fraction to be in the range 7–14% with a typical uncertainty of 4%. A small increase of the SB1 frequency is observed from K- towards F-type stars, in agreement with previous studies. The GES SB1 frequency decreases with metallicity at a rate of (−9 ± 3)% dex−1in the metallicity range −2.7 ≤ [Fe/H] ≤ +0.6. This anticorrelation is obtained with a confidence level higher than 93% on a homogeneous sample covering spectral types FGK and a large range of metallicities. When the present-day mass function is accounted for, this rate turns to (−4 ± 2)% dex−1with a confidence level higher than 88%. In addition we provide the variation of the SB1 fraction with metallicity separately for F, G, and K spectral types, as well as for dwarf and giant primaries.

ABSTRACT 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 γ 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 Gaia-ESO Survey membership probabilities their parallaxes, and kinematics. We find strong evidence for expansion in the OB association population with at least 4σ significance along all three axes, though the expansion is notably anisotropic. We discuss these results in the context of cluster and association dispersal theories.