Search

The typically large distances, extinction, and crowding of Galactic supermassive star clusters have so far hampered the identification of their very low mass members, required to extend our understanding of star and planet formation, and early stellar evolution, to starburst. This situation has now evolved thanks to the James Webb Space Telescope (JWST), and its unmatched resolution and sensitivity in the infrared. In this paper, the third of the series of the Extended Westerlund 1 and 2 Open Clusters Survey (EWOCS), we present JWST/NIRCam and JWST/MIRI observations of the supermassive star cluster Westerlund 1. These observations are specifically designed to unveil the cluster members down to the BD mass regime, and to allow us to select and study the protoplanetary disks and to study the mutual feedback between the cluster members and the surrounding environment. Westerlund 1 was observed as part of JWST GO-1905 for 23.6 hours. The data have been reduced using the JWST calibration pipeline, together with specific tools necessary to remove artifacts. Source identification and photometry were performed with DOLPHOT. The MIRI images show a plethora of different features. Diffuse nebular emission is observed around the cluster, which is typically composed of myriads of droplet-like features pointing toward the cluster center or the group of massive stars surrounding the WR star W72/A. A long pillar is also observed in the NW. The MIRI images also show resolved shells and outflows surrounding the M-type RSG W20, W26, W75, and W237, the sgB[e] star W9 and the YHG W4. The color-magnitude diagrams built using the NIRCam photometry show a clear cluster sequence, which is marked in its upper part by the 1828 NIRCam stars with X-ray counterparts. NIRCam observations using the F115W filter have reached the 23.8 mag limit with 50\% completeness (roughly corresponding to a 0.06 Msol brown dwarf)., Comment: Accepted for publication in Astronomy & Astrophysics

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

Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which aims to investigate the influence of the starburst environment on the formation of stars and planets, and on the evolution of both low and high mass stars. The primary targets of this project are Westerlund 1 and 2, the closest supermassive star clusters to the Sun. Methods. The project is based primarily on recent observations conducted with the Chandra and JWST observatories. Specifically, the Chandra survey of Westerlund 1 consists of 36 new ACIS-I observations, nearly co-pointed, for a total exposure time of 1 Msec. Additionally, we included 8 archival Chandra/ACIS-S observations. This paper presents the resulting catalog of X-ray sources within and around Westerlund 1. Sources were detected by combining various existing methods, and photon extraction and source validation were carried out using the ACIS-Extract software. Results. The EWOCS X-ray catalog comprises 5963 validated sources out of the 9420 initially provided to ACIS-Extract, reaching a photon flux threshold of approximately 2x10^-8 photons/cm^2/s. The X-ray sources exhibit a highly concentrated spatial distribution, with 1075 sources located within the central 1 arcminute. We have successfully detected X-ray emissions from 126 out of the 166 known massive stars of the cluster, and we have collected over 71000 photons from the magnetar CXO J164710.20-455217, Comment: The paper has been accepted for publication by Astronomy and Astrophysics

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

A complete map of the youngest stellar populations of the Milky Way in the era of all-sky surveys, is one of the most challenging goals in modern astrophysics. The characterisation of the youngest stellar component is crucial not only for a global overview of the Milky Way structure, of the Galactic thin disk, and its spiral arms, but also for local studies. In fact, the identification of the star forming regions (SFRs) and the comparison with the environment in which they form are also fundamental to put them in the context of the surrounding giant molecular clouds and to understand still unknown physical mechanisms related to the star and planet formation processes. In 10 yrs of observations, Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will achieve an exquisite photometric depth that will allow us to significantly extend the volume within which we will be able to discover new SFRs and to enlarge the domain of a detailed knowledge of our own Galaxy. We describe here a metrics that estimates the total number of young stars with ages t < 10 Myr and masses >0.3M$_\odot$ that will be detected with the Rubin LSST observations in the gri bands at a 5 {\sigma} magnitude significance. We examine the results of our metrics adopting the most recent simulated Rubin-LSST survey strategies in order to evaluate the impact that different observing strategies might have on our science case., Comment: 11 pages, 5 figures, 1 table; accepted for publication in The Astrophysical Journal Supplement Series

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

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

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

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

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

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

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

The dispersal of protoplanetary disks sets the timescale available for planets to assemble, and thus it is one of the fundamental parameters in theories of planetary formation. Disk dispersal is determined by several properties of the central star, the disk itself, and the surrounding environment. In particular, the metallicity of disks may impact their evolution, even if to date controversial results exist: in low-metallicity clusters disks seem to rapidly disperse, while in the Magellanic Clouds some evidence supports the existence of accreting disks few tens of Myrs old. In this paper we study the dispersal timescale of disks in Dolidze~25, the young cluster in proximity of the Sun with lowest metallicity, with the aim of understanding whether disk evolution is impacted by the low-metallicity of the cluster. We have analyzed Chandra/ACIS-I observations of the cluster and combined the resulting source catalog with existing optical and infrared catalogs of the region. We selected the disk-bearing population and the disk-less population of Dolidze 25. We have derived stellar parameters from isochrones fitted to color-magnitude diagrams. We derived a disk fraction of about 34% and a median age of 1.2 Myrs. By comparing this estimate with existing estimates of the disk fraction of clusters younger than 10 Myrs, our study suggests that the disk fraction of Dolidze 25 is lower than what is expected from its age alone. Even if our results are not conclusive given the intrinsic uncertainty on stellar ages estimated from isochrones fitting to color-magnitude diagrams, we suggest that disk evolution in Dolidze 25 may be impacted by the environment. Given the poor O star population and low stellar density of the cluster, it is more likely that disks dispersal timescale is dictated more by the low metallicity of the cluster rather than external photoevaporation or dynamical encounters., Comment: Accepted for publication by Astronomy & Astrophysics

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

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

Pismis 18 is a moderately populated, intermediate-age open cluster located within the solar circle at a Galactocentric distance of about 7 kpc. Few open clusters have been studied in detail in the inner disc region before the Gaia-ESO Survey. New data from the Gaia-ESO Survey allowed us to conduct an extended radial velocity membership study as well as spectroscopic metallicity and detailed chemical abundance measurements for this cluster. Gaia-ESO Survey data for 142 potential members, lying on the upper MS and on the red clump, yielded radial velocity measurements, which, together with proper motion measurements from the Gaia DR2, were used to determine the systemic velocity of the cluster and membership of individual stars. Photometry from Gaia DR2 was used to re-determine cluster parameters based on high confidence member stars only. Cluster abundance measurements of six radial-velocity member stars with UVES high-resolution spectroscopy are presented for 23 elements. According to the new estimates, based on high confidence members, Pismis 18 has an age of $700^{+40}_{-50}$ Myr, interstellar reddening of E(B-V) = $0.562^{+0.012}_{-0.026}$ mag and a de-reddened distance modulus of $DM_0 = 11.96^{+0.10}_{-0.24}$ mag. The median metallicity of the cluster (using the six UVES stars) is [Fe/H] = $+0.23 \pm 0.05$ dex, with [$\alpha$/Fe]= $0.07 \pm 0.13$ and a slight enhancement of s- and r- neutron-capture elements. With the present work, we fully characterized the open cluster Pismis 18, confirming its present location in the inner disc. We estimated a younger age than the previous literature values and gave, for the first time, its metallicity and its detailed abundances. Its [$\alpha$/Fe] and [s-process/Fe], both slightly super-solar, are in agreement with other inner-disc open clusters observed by the Gaia-ESO survey. [abridged], Comment: 15 pages, 10 figures, 8 tables (Table 2 available at the CDS)

The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts., Comment: Accepted for publication in A&A; 23 pages, 17 figures, 3 tables

(abridged) NGC6530 is a young cluster with a complex morphology and star-formation history. We present a statistical study of its global properties using a new large list of candidate members down to masses of 0.2-0.4 Msun and Gaia DR2 astrometry. We consider a wider sky region compared to previous studies, to investigate the whole cluster and its neighborhood. We study the distribution of extinction and age across the different regions, and obtain constraints on the star-formation history. Cluster membership is determined on the basis of X-ray data, H-alpha emission, near-IR and UV excesses from VPHAS+, UKIDSS and other surveys, and Gaia DR2 astrometry; we also use a novel method for photometric selection of M-type pre-main-sequence cluster members. The list of candidates includes nearly 3700 stars, of which we estimate ~2700 genuine NGC6530 members. Using Gaia parallaxes, the cluster distance is found to be 1325 pc, with errors of 0.5% (statistical) and 8.5% (systematic), in agreement with previous determinations. The cluster morphology and boundaries are established with great confidence, from members selected using different criteria. There is no diffuse population of members around the cluster, but there are minor member subgroups (associated with the stars 7 Sgr (F2II-III) and HD 164536 (O7.5V)) in addition to the two main groups in the cluster core and in the Hourglass nebula. We find a pattern of sequential star formation, between ages of <0.5 Myr and ~5 Myr. Extinction is spatially non-uniform, with part of the population still obscured by thick dust. The precise Gaia proper motion data indicate that the NGC6530 parent cloud collided with the Galactic plane ~4 Myr ago, and we suggest that event as the trigger of the bulk of star formation in NGC6530. The internal cluster dynamics is also partially resolved by the Gaia data, indicating expansion of the main cluster population., Comment: 21 pages, 22 figures, accepted for publication on Astronomy and Astrophysics

One of the aims of LSST is to perform a systematic survey of star clusters and star forming regions (SFRs) in our Galaxy. In particular, the observations obtained with LSST will make a big difference in Galactic regions that have been poorly studied in the past, such as the anticenter and the disk beyond the Galactic center, and they will have a strong impact in discovering new distant SFRs. These results can be achieved by exploiting the exquisite depth that will be attained if the wide-fast-deep (WFD) observing strategy of the main survey is also adopted for the Galactic plane, in the g, r, and i filters., Comment: Submitted White Paper for the cadence optimization call of LSST, November 30 2018

We aim to test a purely photometric approach to statistically identify a young clustered population embedded in a large population of field stars, with no prior knowledge on the nature of stars in the field. We conducted our blind test study on the NGC 2264 region, which hosts a well-known young cluster. We selected a large (4 sq. deg.) area around the NGC 2264 cluster, and assembled an extensive r,i,J catalog of the field from pre-existing large-scale surveys. We then mapped the stellar color locus on the (i-J, r-i) diagram to select M-type stars, which: i) comprise a significant fraction of the Galactic stellar population; ii) exhibit the strongest luminosity evolution from the PMS to the MS; iii) have r,i,J color properties that provide a direct and empirical estimate of A_V. A comparative analysis of the photometric and spatial properties of M-type stars as a function of A_V enabled us to probe the structure and stellar content of our field. We could identify two distinct populations: a diffuse field population and a clustered population. The presence of occulting material, spatially associated with the clustered population, allowed us to derive an estimate of its distance (800-900 pc) and age (~0.5-5 Myr), consistent with the literature parameters for the NGC 2264 star-forming region. The extracted clustered population exhibits a hierarchical structure, in excellent agreement with the NGC 2264 subregions reported in the literature. Our selection of clustered members is coherent with the literature census of the NGC 2264 cluster for about 95% of the objects in the inner regions of the field, where the contamination rate by field stars in our sample is only 2%. The method tested here can be readily applied to surveys like Pan-STARRS and the future LSST to undertake a first complete census of low-mass, young star populations down to distances of several kpc across the Galactic plane., Comment: 19 pages, 22 figures; accepted for publication in Astronomy & Astrophysics; abstract shortened for the ArXiv listing

Sco OB2 is the nearest OB association, extending over approximately 2000 sq.deg. on the sky. Only its brightest members are already known (from Hipparcos) across its entire size, while studies of its lower-mass population refer only to small portions of its extent. In this work we exploit the capabilities of Gaia DR2 measurements to search for Sco OB2 members across its entire size and down to the lowest stellar masses. We use both Gaia astrometric and photometric data to select association members, using minimal assumptions derived mostly from the Hipparcos studies. Gaia resolves small details in both the kinematics of individual Sco OB2 subgroups and their distances from the Sun. We develop methods to explore the 3D kinematics of stellar populations covering large sky areas. We find ~11000 pre-main sequence (PMS) Sco OB2 members (with <3% contamination), plus ~3600 MS candidate members with a larger (10-30%) field-star contamination. A higher-confidence subsample of ~9200 PMS (and ~1340 MS) members is also selected (<1% contamination for the PMS), affected however by larger (~15%) incompleteness. We classify separately stars in compact and diffuse populations. Most members belong to a few kinematically distinct diffuse populations, whose ensemble outlines the association shape. Upper Sco is the densest part of Sco OB2, with a complex spatial and kinematical structure, and no global pattern of motion. Other dense subclusters are found in Upper Centaurus-Lupus and in Lower Centaurus-Crux. Most clustered stars appear to be younger than the diffuse PMS population, suggesting star formation in small groups which rapidly disperse and dilute, while keeping memory of their original kinematics. We also find that the open cluster IC 2602 has a similar dynamics to Sco OB2, and its PMS members are evaporating and forming a ~10 deg halo around its double-peaked core., Comment: 27 pages, 37 figures. Accepted for publication in Astronomy and Astrophysics

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

Most stars born in clusters and recent results suggest that star formation (SF) preferentially occurs in subclusters. Studying the morphology and SF history of young clusters is crucial to understanding early SF. We identify the embedded clusters of young stellar objects (YSOs) down to M stars, in the HII regions RCW33, RCW32 and RCW27 of the Vela Molecular Ridge. Our aim is to characterise their properties, such as morphology and extent of the clusters in the three HII regions, derive stellar ages and the connection of the SF history with the environment. Through public photometric surveys such as Gaia, VPHAS, 2MASS and Spitzer/GLIMPSE, we identify YSOs with IR, Halpha and UV excesses, as signature of circumstellar disks and accretion. In addition, we implement a method to distinguish M dwarfs and giants, by comparing the reddening derived in several optical/IR color-color diagrams, assuming suitable theoretical models. Since this diagnostic is sensitive to stellar gravity, the procedure allows us to identify pre-main sequence stars. We find a large population of YSOs showing signatures of circumstellar disks with or without accretion. In addition, with the new technique of M-type star selection, we find a rich population of young M stars with a spatial distribution strongly correlated to the more massive population. We find evidence of three young clusters, with different morphology. In addition, we identify field stars falling in the same region, by securely classifying them as giants and foreground MS stars. We identify the embedded population of YSOs, down to about 0.1 Msun, associated with the HII regions RCW33, RCW32 and RCW27 and the clusters Vela T2, Cr197 and Vela T1, respectively, showing very different morphologies. Our results suggest a decreasing SF rate in Vela T2 and triggered SF in Cr197 and Vela T1., Comment: Accepted for publication in A&A; 20 pages, 22 figures, 6 tables

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

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

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

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

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

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

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

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

Abridged: We used the fourth internal data release of the Gaia-ESO survey to characterize the bulge chemistry, spatial distribution, kinematics, and to compare it chemically with the thin and thick disks. The sample consist on ~2500 red clump stars in 11 bulge fields ($-10^\circ\leq l\leq+8^\circ$ and $-10^\circ\leq b\leq-4^\circ$), and a set of ~6300 disk stars selected for comparison. The bulge MDF is confirmed to be bimodal across the whole sampled area, with metal-poor stars dominating at high latitudes. The metal-rich stars exhibit bar-like kinematics and display a bimodality in their magnitude distribution, a feature which is tightly associated with the X-shape bulge. They overlap with the metal-rich end of the thin disk sequence in the [Mg/Fe] vs. [Fe/H] plane. Metal-poor bulge stars have a more isotropic hot kinematics and do not participate in the X-shape bulge. With similar Mg-enhancement levels, the position of the metal-poor bulge sequence "knee" is observed at [Fe/H]$_{knee}=-0.37\pm0.09$, being 0.06 dex higher than that of the thick disk. It suggests a higher SFR for the bulge than for the thick disk. Finally, we present a chemical evolution model that suitably fits the whole bulge sequence by assuming a fast ($<1$ Gyr) intense burst of stellar formation at early epochs. We associate metal-rich stars with the B/P bulge formed from the secular evolution of the early thin disk. On the other hand, the metal-poor subpopulation might be the product of an early prompt dissipative collapse dominated by massive stars. Nevertheless, our results do not allow us to firmly rule out the possibility that these stars come from the secular evolution of the early thick disk. This is the first time that an analysis of the bulge MDF and $\alpha$-abundances has been performed in a large area on the basis of a homogeneous, fully spectroscopic analysis of high-resolution, high S/N data., Comment: Accepted in A&A 30 March 2017

Due to their volatile nature, when sulfur and zinc are observed in external galaxies, their determined abundances represent the gas-phase abundances in the interstellar medium. This implies that they can be used as tracers of the chemical enrichment of matter in the Universe at high redshift. Comparable observations in stars are more difficult and, until recently, plagued by small number statistics. We wish to exploit the Gaia ESO Survey (GES) data to study the behaviour of sulfur and zinc abundances of a large number of Galactic stars, in a homogeneous way. By using the UVES spectra of the GES sample, we are able to assemble a sample of 1301 Galactic stars, including stars in open and globular clusters in which both sulfur and zinc were measured. We confirm the results from the literature that sulfur behaves as an alpha-element. We find a large scatter in [Zn/Fe] ratios among giant stars around solar metallicity. The lower ratios are observed in giant stars at Galactocentric distances less than 7.5 kpc. No such effect is observed among dwarf stars, since they do not extend to that radius. Given the sample selection, giants and dwarfs are observed at different Galactic locations, and it is plausible, and compatible with simple calculations, that Zn-poor giants trace a younger population more polluted by SN Ia yields. It is necessary to extend observations in order to observe both giants and dwarfs at the same Galactic location. Further theoretical work on the evolution of zinc is also necessary., Comment: 15 pages, 17 figures, accepted for publication in A&A

The spatial distribution of elemental abundances in the disc of our Galaxy gives insights both on its assembly process and subsequent evolution, and on the stellar nucleogenesis of the different elements. Gradients can be traced using several types of objects as, for instance, (young and old) stars, open clusters, HII regions, planetary nebulae. We aim at tracing the radial distributions of abundances of elements produced through different nucleosynthetic channels -the alpha-elements O, Mg, Si, Ca and Ti, and the iron-peak elements Fe, Cr, Ni and Sc - by using the Gaia-ESO idr4 results of open clusters and young field stars. From the UVES spectra of member stars, we determine the average composition of clusters with ages >0.1 Gyr. We derive statistical ages and distances of field stars. We trace the abundance gradients using the cluster and field populations and we compare them with a chemo-dynamical Galactic evolutionary model. Results. The adopted chemo-dynamical model, with the new generation of metallicity-dependent stellar yields for massive stars, is able to reproduce the observed spatial distributions of abundance ratios, in particular the abundance ratios of [O/Fe] and [Mg/Fe] in the inner disc (5 kpc

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

Milky Way open clusters are very diverse in terms of age, chemical composition, and kinematic properties. Intermediate-age and old open clusters are less common, and it is even harder to find them inside the solar Galactocentric radius, due to the high mortality rate and strong extinction inside this region. NGC 6802 is one of the inner disk open clusters (IOCs) observed by the $Gaia$-ESO survey (GES). This cluster is an important target for calibrating the abundances derived in the survey due to the kinematic and chemical homogeneity of the members in open clusters. Using the measurements from $Gaia$-ESO internal data release 4 (iDR4), we identify 95 main-sequence dwarfs as cluster members from the GIRAFFE target list, and eight giants as cluster members from the UVES target list. The dwarf cluster members have a median radial velocity of $13.6\pm1.9$ km s$^{-1}$, while the giant cluster members have a median radial velocity of $12.0\pm0.9$ km s$^{-1}$ and a median [Fe/H] of $0.10\pm0.02$ dex. The color-magnitude diagram of these cluster members suggests an age of $0.9\pm0.1$ Gyr, with $(m-M)_0=11.4$ and $E(B-V)=0.86$. We perform the first detailed chemical abundance analysis of NGC 6802, including 27 elemental species. To gain a more general picture about IOCs, the measurements of NGC 6802 are compared with those of other IOCs previously studied by GES, that is, NGC 4815, Trumpler 20, NGC 6705, and Berkeley 81. NGC 6802 shows similar C, N, Na, and Al abundances as other IOCs. These elements are compared with nucleosynthetic models as a function of cluster turn-off mass. The $\alpha$, iron-peak, and neutron-capture elements are also explored in a self-consistent way., Comment: 18 pages, 12 Figures, A&A accepted

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

NGC6231 is a massive young star cluster, near the center of the Sco OB1 association. While its OB members are well studied, its low-mass population has received little attention. We present high-spatial resolution Chandra ACIS-I X-ray data, where we detect 1613 point X-ray sources. Our main aim is to clarify global properties of NGC6231 down to low masses through a detailed membership assessment, and to study the cluster stars' spatial distribution, the origin of their X-ray emission, the cluster age and formation history, and initial mass function. We use X-ray data, complemented by optical/IR data, to establish cluster membership. The spatial distribution of different stellar subgroups also provides highly significant constraints on cluster membership, as does the distribution of X-ray hardness. We perform spectral modeling of group-stacked X-ray source spectra. We find a large cluster population down to ~0.3 Msun (complete to ~1 Msun), with minimal non-member contamination, with a definite age spread (1-8 Myrs) for the low-mass PMS stars. We argue that low-mass cluster stars also constitute the majority of the few hundreds unidentified X-ray sources. We find mass segregation for the most massive stars. The fraction of circumstellar-disk bearing members is found to be ~5%. Photoevaporation of disks under the action of massive stars is suggested by the spatial distribution of the IR-excess stars. We also find strong Halpha emission in 9% of cluster PMS stars. The dependence of X-ray properties on mass, stellar structure, and age agrees with extrapolations based on other young clusters. The cluster initial mass function, computed over ~2 dex in mass, has a slope Gamma~-1.14. The total mass of cluster members above 1 Msun is 2280 Msun, and the inferred total mass is 4380 Msun. We also study the peculiar, hard X-ray spectrum of the Wolf-Rayet star WR79., Comment: 25 pages, 36 figures, accepted for publication on Astronomy and Astrophysics

We investigate here the lithium content of young low-mass stars in the 5~Myr-old star forming region NGC~2264 and its relationship with rotation. We combine lithium equivalent width measurements, EW(Li), from the Gaia-ESO Survey with the determination of rotational periods from the CSI~2264 survey. We consider only bona fide non accreting cluster members in order to minimize uncertainties on EW(Li). We report the existence of a relationship between lithium content and rotation in NGC~2264 at an age of 5~Myr. The Li-rotation connection is seen over a restricted temperature range, Teff=3800--4400K, where fast rotators are Li-rich compared to slow ones. This correlation is similar to, albeit of lower amplitude than, the Li-rotation connection previously reported for K dwarfs in the 125 Myr-old Pleiades cluster. We investigate whether the non-standard pre-main sequence models developed so far to explain the Pleiades results, which are based on episodic accretion, pre-main sequence core-envelope decoupling, and/or radius inflation due to enhanced magnetic activity, can account for an early development of the Li-rotation connection. While radius inflation appears to be the most promising possibility, each of these models has issues. We therefore also discuss external causes that might operate during the first few Myr of pre-main sequence evolution, such as planet engulfment and/or steady disk accretion, as possible candidates for the common origin for Li-excess and fast rotation in young low-mass pre-main sequence stars. The emergence of a connection between lithium content and rotation rate at such an early age as 5~Myr suggests a complex link between accretion processes, early angular momentum evolution, and possibly planet formation, which likely impacts early stellar evolution and still is to be fully deciphered., Comment: 9 pages, 5 figures, Astronomy & Astrophysics, in press

We present observations from the Gaia-ESO Survey in the lines of H$\alpha$, [N II], [S II] and He I of nebular emission in the central part of the Carina Nebula. We investigate the properties of the two already known kinematic components (approaching and receding, respectively), which account for the bulk of emission. Moreover, we investigate the features of the much less known low-intensity high-velocity (absolute RV $>$50 km/s) gas emission. We show that gas giving rise to H$\alpha$ and He I emission is dynamically well correlated, but not identical, to gas seen through forbidden-line emission. Gas temperatures are derived from line-width ratios, and densities from [S II] doublet ratios. The spatial variation of N ionization is also studied, and found to differ between the approaching and receding components. The main result is that the bulk of the emission lines in the central part of Carina arises from several distinct shell-like expanding regions, the most evident found around $\eta$ Car, the Trumpler 14 core, and the star WR25. Such "shells" are non-spherical, and show distortions probably caused by collisions with other shells or colder, higher-density gas. Part of them is also obscured by foreground dust lanes, while only very little dust is found in their interior. Preferential directions, parallel to the dark dust lanes, are found in the shell geometries and physical properties, probably related to strong density gradients in the studied region. We also find evidence that the ionizing flux emerging from $\eta$ Car and the surrounding Homunculus nebula varies with polar angle. The high-velocity components in the wings of H$\alpha$ are found to arise from expanding dust reflecting the $\eta$ Car spectrum., Comment: 25 pages, 42 figures, accepted for publication in Astronomy and Astrophysics

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