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.

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

Context. Taurus represents an ideal region to study the three-dimensional distribution of the young stellar population and relate it to the associated molecular cloud. Aims. The second Gaia data release (DR2) enables us to investigate the Taurus complex in three dimensions, starting from a previously defined robust membership. The molecular cloud structured in filaments can be traced in emission using the public far-infrared maps from Herschel. Methods. From a compiled catalog of spectroscopically confirmed members, we analyze the 283 sources with reliable parallax and proper motions in the Gaia DR2 archive. We fit the distribution of parallaxes and proper motions with multiple populations described by multivariate Gaussians. We compute the cartesian Galactic coordinates (X,Y,Z) and, for the populations associated with the main cloud, also the galactic space velocity (U,V,W). We discuss the spatial distribution of the populations in relation to the structure of the filamentary molecular cloud traced by Herschel. Results. We discover the presence of six populations which are all well defined in parallax and proper motions, with the only exception being Taurus D. The derived distances range between 130 and 160 pc. We do not find a unique relation between stellar population and the associated molecular cloud: while the stellar population seems to be on the cloud surface, both lying at similar distances, this is not the case when the molecular cloud is structured in filaments. Taurus B is probably moving in the direction of Taurus A, while Taurus E appears to be moving towards them. Conclusions. The Taurus region is the result of a complex star formation history which most probably occurred in clumpy and filamentary structures that are evolving independently.

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<FeH<+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 with 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.

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

The past two decades have seen dramatic progress in our knowledge of the population of young stars of age $< 200\,$Myr that lie within $150\,$pc of the Sun. These nearby, young stars, most of which are found in loose, comoving groups, provide the opportunity to explore (among many other things) the dissolution of stellar clusters and their diffusion into the field star population. Here, we exploit the combination of astrometric and photometric data from {\it Gaia} and photometric data from GALEX (UV) and 2MASS (near-IR) in an attempt to identify additional nearby, young, late-type stars. Specifically, we present a sample of 146 GALEX UV-selected late-type (predominantly K-type) field stars with {\it Gaia}-based distances $< 125\,$pc (based on {\it Gaia} Data Release 1) that have isochronal ages $< 80\,$Myr even if equal-components binaries. We investigate the spectroscopic and kinematic properties of this sample. Despite their young isochronal ages, only $\sim 10$ per cent) of stars among this sample can be confidently associated with established nearby, young moving groups (MGs). These candidate MG members include 5 stars newly identified in this study. The vast majority of our sample of 146 nearby young star candidates have anomalous kinematics relative to the known MGs. These stars may hence represent a previously unrecognised population of young stars that has recently mixed into the older field star population. We discuss the implications and caveats of such a hypothesis---including the intriguing fact that, in addition to their non-young-star-like kinematics, the majority of the UV-selected, isochronally young field stars within $50\,$pc appear surprisingly X-ray faint.

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.

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.

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.

The past two decades have seen dramatic progress in our knowledge of the population of young stars of age <200 Myr that lie within 150 pc of the Sun. These nearby, young stars, most of which are found in loose, comoving groups, provide the opportunity to explore (among many other things) the dissolution of stellar clusters and their diffusion into the field star population. In the age of Gaia, this potential can now be fully exploited. We have identified, and are now investigating, a sample of nearly 400 Galex UV-selected late-type (K and early-M) field stars with Gaia-based distances <120 pc and isochronal ages <=80 Myr (even if binaries). Only a small percentage (<10%) of stars among this (kinematically unbiased) sample can be confidently associated with established nearby, young moving groups (NYMGs). The majority display anomalous kinematics, relative to the known NYMGs. These stars may hence represent a previously unrecognized population of young stars that has recently mixed into the older field star population. We discuss the implications and caveats of such a hypothesis---including the intriguing fact that, in addition to their non-young-star-like kinematics, the majority of the UV-selected, isochronally young field stars within 50 pc appear surprisingly X-ray faint.

Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.

The main aim of this experiment is to provide a complete census of old (t > 10 Gyr, RR Lyrae, type II Cepheids, red horizontal branch), intermediate age (red clump, Miras) and young (classical Cepheids) stellar tracers across the Galactic Bulge. To fully exploit the unique photometric quality of LSST images, we plan to perform a Shallow minisurvey (ugrizy, -20 < l < 20 deg, -15 < b < 10 deg) and a Deep minisurvey (izy, -20 < l < 20 deg, -3 < b < 3 deg). The former one is aimed at constraining the 3D structure of the galactic Bulge across the four quadrants, and in particular, the transition between inner and outer Bulge. The u,g,r,i,z,y LSST bands provide fundamental diagnostics to constrain the evolutionary properties of low and intermediate-mass stars when moving from a metal-poor to a metal-rich regime. The deep minisurvey is aimed at tracing RR Lyrae, Red Clump stars, Miras and classical Cepheids in highly reddened regions of the Galactic center. These images will allow us to investigate the role that baryonic mass and dark matter played in the early formation and evolution of the MW.

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.

The Vela OB2 association is a group of 10 Myr stars exhibiting a complex spatial and kinematic substructure. The all-sky Gaia DR2 catalogue contains proper motions, parallaxes (a proxy for distance) and photometry that allow us to separate the various components of Vela OB2. We characterise the distribution of the Vela OB2 stars on a large spatial scale, and study its internal kinematics and dynamic history. We make use of Gaia DR2 astrometry and published Gaia-ESO Survey data. We apply an unsupervised classification algorithm to determine groups of stars with common proper motions and parallaxes. We find that the association is made up of a number of small groups, with a total current mass over 2330 Msun. The three-dimensional distribution of these young stars trace the edge of the gas and dust structure known as the IRAS Vela Shell across 180 pc and shows clear signs of expansion. We propose a common history for Vela OB2 and the IRAS Vela Shell. The event that caused the expansion of the shell happened before the Vela OB2 stars formed, imprinted the expansion in the gas the stars formed from, and most likely triggered star formation.

Context. Chamaeleon I represents an ideal laboratory to study the cluster formation in a low-mass environment. Recently, two sub clusters spatially located in the northern and southern parts of Chamaeleon I were found with different ages and radial velocities. Aims. In this letter we report new insights into the structural properties, age, and distance of Chamaeleon I based on the astrometric parameters from Gaia data-release 2 (DR2). Methods. We identified 140 sources with a reliable counterpart in the Gaia DR2 archive. We determined the median distance of the cluster using Gaia parallaxes and fitted the distribution of parallaxes and proper motions assuming the presence of two clusters. We derived the probability of each single source of belonging to the northern or southern sub-clusters, and compared the HR diagram of the most probable members to pre-main sequences isochrones. Results. The median distance of Chamaeleon I is ~190 pc. This is about 20 pc larger than the value commonly adopted in the literature. From a Kolmogorov-Smirnov test of the parallaxes and proper-motion distributions we conclude that the northern and southern clusters do not belong to the same parent population. The northern population has a distance dN = 192.7+/-0.4 pc, while the southern one dS = 186.5+/-0.7 pc. The two sub-clusters appear coeval, at variance with literature results, and most of the sources are younger than 3 Myr. The northern cluster is more elongated and extends towards the southern direction partially overlapping with the more compact cluster located in the south. A hint of a relative rotation between the two sub-clusters is also found.

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.

Observations of the main sequence F3 V star KIC 8462852 (also known as Boyajian's star) revealed extreme aperiodic dips in flux up to 20% during the four years of the Kepler mission. Smaller dips (< 2%) were also observed with ground-based telescopes between May and September 2017. We investigated possible correlation between recent dips and the major dips in the last 100 days of the Kepler mission. We compared Kepler light curve data, 2017 data from two observatories (TFN, OGG) which are part of the Las Cumbres Observatory (LCO) network and Sternberg observatory archival data, and determined that observations are consistent with a 1,574-day (4.31 year) periodicity of a transit (or group of transits) orbiting Boyajian's star within the habitable zone. It is unknown if transits that have produced other major dips as observed during the Kepler mission (e.g. D792) share the same orbital period. Nevertheless, the proposed periodicity is a step forward in guiding future observation efforts.

The abundance ratio N/O is a useful tool to study the interplay of galactic processes, e.g. star formation efficiency, time-scale of infall and outflow loading factor We aim to trace log(N/O) versus [Fe/H] in the Milky Way and to compare it with a set of chemical evolution models to understand the role of infall, outflow and star formation efficiency in the building-up of the Galactic disc. We use the abundances from idr2-3, idr4, idr5 data releases of the Gaia-ESO Survey both for Galactic field and open cluster stars.We determine membership and average composition of open clusters and we separate thin and thick disc field stars.We consider the effect of mixing in the abundance of N in giant stars. We compute a grid of chemical evolution models, suited to reproduce the main features of our Galaxy, exploring the effects of the star formation efficiency, the infall time-scale and the differential outflow. With our samples, we map the metallicity range -0.6<[Fe/H]<0.3 with a corresponding -1.2<log(N/O)<-0.2, where the secondary production of N dominates. Thanks to the wide range of Galactocentric distances covered by our samples, we can distinguish the behaviour of log(N/O) in different parts of the Galaxy. Our spatially resolved results allow us to distinguish differences in the evolution of N/O with Galactocentric radius. Comparing the data with our models, we can characterise the radial regions of our Galaxy. A shorter infall time-scale is needed in the inner regions, while the outer regions need a longer infall time-scale, coupled with a higher star formation efficiency. We compare our results with nebular abundances obtained in MaNGA galaxies, finding in our Galaxy a much wider range of log(N/O) than in integrated observations of external galaxies of similar stellar mass, but similar to the ranges found in studies of individual H ii regions.

We present a high resolution synthetic spectral library, INTRIGOSS, designed for studying FGK stars. The library is based on atmosphere models computed with specified individual element abundances via ATLAS12 code. Normalized SPectra (NSP) and surface Flux SPectra (FSP), in the 4830-5400 A, wavelength range, were computed with the SPECTRUM code. INTRIGOSS uses the solar composition by Grevesse et al. 2007 and four [alpha/Fe] abundance ratios and consists of 15,232 spectra. The synthetic spectra are computed with astrophysical gf-values derived by comparing synthetic predictions with a very high SNR solar spectrum and the UVES-U580 spectra of five cool giants. The validity of the NSPs is assessed by using the UVES-U580 spectra of 2212 stars observed in the framework of the Gaia-ESO Survey and characterized by homogeneous and accurate atmospheric parameter values and by detailed chemical compositions. The greater accuracy of NSPs with respect to spectra from the AMBRE, GES_Grid, PHOENIX, C14, and B17 synthetic spectral libraries is demonstrated by evaluating the consistency of the predictions of the different libraries for the UVES-U580 sample stars. The validity of the FSPs is checked by comparing their prediction with both observed spectral energy distribution and spectral indices. The comparison of FSPs with SEDs derived from ELODIE, INDO--U.S., and MILES libraries indicates that the former reproduce the observed flux distributions within a few percent and without any systematic trend. The good agreement between observational and synthetic Lick/SDSS indices shows that the predicted blanketing of FSPs well reproduces the observed one, thus confirming the reliability of INTRIGOSS FSPs.

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.

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.

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.

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

We attempt to determine the relative fraction of stars that have undergone significant radial migration by studying the orbital properties of metal-rich ([Fe/H]$>0.1$) stars within 2 kpc of the Sun using a sample of more than 3,000 stars selected from iDR4 of the Gaia-ESO Survey. We investigate the kinematic properties, such as velocity dispersion and orbital parameters, of stellar populations near the sun as a function of [Mg/Fe] and [Fe/H], which could show evidence of a major merger in the past history of the Milky Way. This was done using the stellar parameters from the Gaia-ESO Survey along with proper motions from PPMXL to determine distances, kinematics, and orbital properties for these stars to analyze the chemodynamic properties of stellar populations near the Sun. Analyzing the kinematics of the most metal-rich stars ([Fe/H]$>0.1$), we find that more than half have small eccentricities ($e<0.2$) or are on nearly circular orbits. Slightly more than 20\% of the metal-rich stars have perigalacticons $R_p>7$ kpc. We find that the highest [Mg/Fe], metal-poor populations have lower vertical and radial velocity dispersions compared to lower [Mg/Fe] populations of similar metallicity by $\sim10$ km s$^{-1}$. The median eccentricity increases linearly with [Mg/Fe] across all metallicities, while the perigalacticon decreases with increasing [Mg/Fe] for all metallicities. Finally, the most [Mg/Fe]-rich stars are found to have significant asymmetric drift and rotate more than 40 km s$^{-1}$ slower than stars with lower [Mg/Fe] ratios. While our results cannot constrain how far stars have migrated, we propose that migration processes are likely to have played an important role in the evolution of the Milky Way, with metal-rich stars migrating from the inner disk toward to solar neighborhood and past mergers potentially driving enhanced migration of older stellar populations in the disk.

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.

The typical methodology for comparing simulated galaxies with observational surveys is usually to apply a spatial selection to the simulation to mimic the region of interest covered by a comparable observational survey sample. In this work we compare this approach with a more sophisticated post-processing in which the observational uncertainties and selection effects (photometric, surface gravity and effective temperature) are taken into account. We compare a `solar neighbourhood analogue' region in a model Milky Way-like galaxy simulated with RAMSES-CH with fourth release Gaia-ESO survey data. We find that a simple spatial cut alone is insufficient and that observational uncertainties must be accounted for in the comparison. This is particularly true when the scale of uncertainty is large compared to the dynamic range of the data, e.g. in our comparison, the [Mg/Fe] distribution is affected much more than the more accurately determined [Fe/H] distribution. Despite clear differences in the underlying distributions of elemental abundances between simulation and observation, incorporating scatter to our simulation results to mimic observational uncertainty produces reasonable agreement. The quite complete nature of the Gaia-ESO survey means that the selection function has minimal impact on the distribution of observed age and metal abundances but this would become increasingly more important for surveys with narrower selection functions.

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.

The Galex Nearby Young Star Survey (GALNYSS) has yielded a sample of $\sim$2000 UV-selected objects that are candidate nearby ($D \stackrel{<}{\sim}$150 pc), young (age $\sim$10--100 Myr), late-type stars. Here, we evaluate the distances and ages of the subsample of (19) GALNYSS stars with Gaia Data Release 1 (DR1) parallax distances $D \le 120$ pc. The overall youth of these 19 mid-K to early-M stars is readily apparent from their positions relative to the loci of main sequence stars and giants in Gaia-based color-magnitude and color-color diagrams constructed for all Galex- and WISE-detected stars with parallax measurements included in DR1. The isochronal ages of all 19 stars lie in the range $\sim$10--100 Myr. Comparison with Li-based age estimates indicates a handful of these stars may be young main-sequence binaries rather than pre-main sequence stars. Nine of the 19 objects have not previously been considered as nearby, young stars, and all but one of these are found at declinations north of $+$30$^\circ$. The Gaia DR1 results presented here indicate that the GALNYSS sample includes several hundred nearby, young stars, a substantial fraction of which have not been previously recognized as having ages $\stackrel{<}{\sim}$100 Myr.

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.

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.

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<RGC <7 kpc), with their differences, that were usually poorly explained by chemical evolution models. Often, oxygen and magnesium are considered as equivalent in tracing alpha-element abundances and in deducing, e.g., the formation time-scales of different Galactic stellar populations. In addition, often [alpha/Fe] is computed combining several alpha-elements. Our results indicate, as expected, a complex and diverse nucleosynthesis of the various alpha-elements, in particular in the high metallicity regimes, pointing towards a different origin of these elements and highlighting the risk of considering them as a single class with common features.

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.

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.

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.

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.

We present Chandra X-ray Observatory and Submillimeter Array (SMA) imaging of HBC 515, a system consisting of multiple young stellar objects (YSOs). The five members of HBC 515 represent a remarkably diverse array of YSOs, ranging from the low-mass Class I/II protostar HBC 515B, through Class II and transition disk objects (HBC 515D and C, respectively), to the "diskless", intermediate- mass, pre-main sequence binary HBC 515A. Our Chandra/ACIS imaging establishes that all five components are X-ray sources, with HBC 515A - a subarcsecond-separation binary that is partially resolved by Chandra - being the dominant X-ray source. We detect an X-ray flare associated with HBC 515B. In the SMA imaging, HBC 515B is detected as a strong 1.3 mm continuum emission source; a second, weaker mm continuum source is coincident with the position of the transition disk object HBC 515C. These results strongly support the protostellar nature of HBC 515B, and firmly establish HBC 515A as a member of the rare class of relatively massive, X-ray luminous "weak-lined T Tauri stars" that are binaries and have shed their disks at very early stages of pre-MS evolution. The coexistence of two such disparate objects within a single, presumably coeval multiple YSO system highlights the influence of pre- MS star mass, binarity, and X-ray luminosity in regulating the lifetimes of circumstellar, planet-forming disks and the timescales of star-disk interactions.

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.2<M/Msun<0.7, are at least a factor of two older and ~7 per cent cooler than previously thought and that their masses are much larger (by >30 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.

The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

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.

Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known. Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.

We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy-peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal rich stars also exhibit peaky features in their line-of-sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy-peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The 'banana' or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data.

The nature of the metallicity gradient inside the solar circle (R_GC < 8 kpc) is poorly understood, but studies of Cepheids and a small sample of open clusters suggest that it steepens in the inner disk. We investigate the metallicity gradient of the inner disk using a sample of inner disk open clusters that is three times larger than has previously been studied in the literature to better characterize the gradient in this part of the disk. We used the Gaia-ESO Survey (GES) [Fe/H] values and stellar parameters for stars in 12 open clusters in the inner disk from GES-UVES data. Cluster mean [Fe/H] values were determined based on a membership analysis for each cluster. Where necessary, distances and ages to clusters were determined via comparison to theoretical isochrones. The GES open clusters exhibit a radial metallicity gradient of -0.10+-0.02 dex/kpc, consistent with the gradient measured by other literature studies of field red giant stars and open clusters in the range R_GC ~ 6-12 kpc. We also measure a trend of increasing [Fe/H] with increasing cluster age, as has also been found in the literature. We find no evidence for a steepening of the inner disk metallicity gradient inside the solar circle as earlier studies indicated. The age-metallicity relation shown by the clusters is consistent with that predicted by chemical evolution models that include the effects of radial migration, but a more detailed comparison between cluster observations and models would be premature.

We report the discovery of two Li-rich giants, with A(Li) ~ 1.50, in an analysis of a sample of 40 giants of the open cluster Trumpler 20 (with turnoff mass ~ 1.8 Msun). The cluster was observed in the context of the Gaia-ESO Survey. The atmospheric parameters and Li abundances were derived using high-resolution UVES spectra. The Li abundances were corrected for nonlocal thermodynamical equilibrium (non-LTE) effects. Only upper limits of the Li abundance could be determined for the majority of the sample. Two giants with detected Li turned out to be Li rich: star MG 340 has A(Li) non-LTE = 1.54 \pm 0.21 dex and star MG 591 has A(Li) non-LTE = 1.60 \pm 0.21 dex. Star MG 340 is on average ~ 0.30 dex more rich in Li than stars of similar temperature, while for star MG 591 this difference is on average ~ 0.80 dex. Carbon and nitrogen abundances indicate that all stars in the sample have completed the first dredge-up. The Li abundances in this unique sample of 40 giants in one open cluster clearly show that extra mixing is the norm in this mass range. Giants with Li abundances in agreement with the predictions of standard models are the exception. To explain the two Li-rich giants, we suggest that all events of extra mixing have been inhibited. This includes rotation-induced mixing during the main sequence and the extra mixing at the red giant branch luminosity bump. Such inhibition has been suggested in the literature to occur because of fossil magnetic fields in red giants that are descendants of main-sequence Ap-type stars.

The Gaia-ESO Survey was designed to target all major Galactic components (i.e., bulge, thin and thick discs, halo and clusters), with the goal of constraining the chemical and dynamical evolution of the Milky Way. This paper presents the methodology and considerations that drive the selection of the targeted, allocated and successfully observed Milky Way field stars. The detailed understanding of the survey construction, specifically the influence of target selection criteria on observed Milky Way field stars is required in order to analyse and interpret the survey data correctly. We present the target selection process for the Milky Way field stars observed with VLT/FLAMES and provide the weights that characterise the survey target selection. The weights can be used to account for the selection effects in the Gaia-ESO Survey data for scientific studies. We provide a couple of simple examples to highlight the necessity of including such information in studies of the stellar populations in the Milky Way.

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.

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.

We present contemporaneous XMM-Newton X-ray and ground-based optical/near-IR spectroscopic observations of the nearby ($D \approx 42$ pc), low-mass (mid-M) binary system TWA 30A and 30B. The components of this wide (separation $\sim$3400 AU) binary are notable for their nearly edge-on disk viewing geometries, high levels of variability, and evidence for collimated stellar outflows. We obtained XMM-Newton X-ray observations of TWA 30A and 30B in 2011 June and July, accompanied (respectively) by IRTF SpeX (near-IR) and VLT XSHOOTER (visible/near-IR) spectroscopy obtained within $\sim$20 hours of the X-ray observations. TWA 30A was detected in both XMM-Newton observations at relatively faint intrinsic X-ray luminosities ($L_{X}$$\sim$$8\times10^{27}$ $erg$ $s^{-1}$) compared to stars of similar mass and age . The intrinsic (0.15-2.0 keV) X-ray luminosities measured in 2011 had decreased by a factor 20-100 relative to a 1990 (ROSAT) X-ray detection. TWA 30B was not detected, and we infer an upper limit of ($L_{X}$ $\lesssim$ 3.0 $\times$ $10^{27}$ erg s$^{-1}$). We measured a large change in visual extinction toward TWA 30A (from $A_V \approx 14.9$ to $A_V \approx 4.7$) between the two 2011 observing epochs, and we find evidence for a corresponding significant decrease in X-ray absorbing column ($N_H$). The apparent correlated change in $A_V$ and $N_H$ is suggestive of variable obscuration of the stellar photosphere by disk material composed of both gas and dust. However, in both observations, the inferred $N_{H}$ to $A_{V}$ ratio is lower than that typical of the ISM, suggesting that the disk is either depleted of gas or is deficient in metals in the gas phase.

Stellar evolution models predict that internal mixing should cause some sodium overabundance at the surface of red giants more massive than ~ 1.5--2.0 Msun. The surface aluminium abundance should not be affected. Nevertheless, observational results disagree about the presence and/or the degree of the Na and Al overabundances. In addition, Galactic chemical evolution models adopting different stellar yields lead to quite different predictions for the behavior of [Na/Fe] and [Al/Fe] versus [Fe/H]. Overall, the observed trends of these abundances with metallicity are not well reproduced. We readdress both issues, using new Na and Al abundances determined within the Gaia-ESO Survey, using two samples: i) more than 600 dwarfs of the solar neighborhood and of open clusters and ii) low- and intermediate-mass clump giants in six open clusters. Abundances of Na in giants with mass below ~2.0 Msun, and of Al in giants below ~3.0 Msun, seem to be unaffected by internal mixing processes. For more massive giants, the Na overabundance increases with stellar mass. This trend agrees well with predictions of stellar evolutionary models. Chemical evolution models that are able to fit well the observed [Na/Fe] vs. [Fe/H] trend in solar neighborhood dwarfs can not simultaneously explain the run of [Al/Fe] with [Fe/H], and viceversa. The comparison with stellar ages is hampered by severe uncertainties. Indeed, reliable age estimates are available for only a half of the stars of the sample. We conclude that Al is underproduced by the models, except for stellar ages younger than about 7 Gyr. In addition, some significant source of late Na production seems to be missing in the models. Either current Na and Al yields are affected by large uncertainties, and/or some important Galactic source(s) of these elements has not been taken into account up to now. [abridged]

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<M/M_\odot <1.3$ and $\alpha=1.1\pm0.4$ for $0.16<M/M_\odot <0.5$ and is consistent with a standard IMF. The similarity of the IMF of the young population around $\gamma^2 $Vel to that in other star forming regions and the field suggests it may have formed through very similar processes.

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

(Abridged) Recent spectroscopic surveys have begun to explore the Galactic disk system outside the solar neighborhood on the basis of large data samples. In this way, they provide valuable information for testing spatial and temporal variations of disk structure kinematics and chemical evolution. We used a Gaussian mixture model algorithm, as a rigurous mathematical approach, to separate in the [Mg/Fe] vs. [Fe/H] plane a clean disk star subsample from the Gaia-ESO survey internal data release 2. We find that the sample is separated into five groups associated with major Galactic components; the metal-rich end of the halo, the thick disk, and three subgroups for the thin disk sequence. This is confirmed with a sample of red clump stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The two metal-intermediate and metal-rich groups of the thin disk decomposition ([Fe/H]>-0.25 dex) highlight a change in the slope at solar metallicity. This holds true at different radial regions. The distribution of Galactocentric radial distances of the metal-poor part of the thin disk ([Fe/H]<-0.25 dex) is shifted to larger distances than those of the more metal-rich parts. Moreover, the metal-poor part of the thin disk presents indications of a scale height intermediate between those of the thick and the rest of the thin disk, and it displays higher azimuthal velocities than the latter. These stars might have formed and evolved in parallel and/or dissociated from the inside-out formation taking place in the internal thin disk. Their enhancement levels might be due to their origin from gas pre-enriched by outflows from the thick disk or the inner halo. The smooth trends of their properties (their spatial distribution with respect to the plane, in particular) with [Fe/H] and [Mg/Fe] suggested by the data indicates a quiet dynamical evolution, with no relevant merger events.

Understanding the history and the evolution of the Milky Way disc is one of the main goals of modern astrophysics. We study the velocity dispersion behaviour of Galactic disc stars as a function of the [Mg/Fe] ratio, which can be used as a proxy of relative age. This key relation is essential to constrain the formation mechanisms of the disc stellar populations as well as the cooling processes. We used the recommended parameters and chemical abundances of 7800 FGK Milky Way field stars from the second internal data release of the Gaia-ESO Survey. These stars were observed with the GIRAFFE spectrograph, and cover a large spatial volume (6<R<10kpc and |Z|<2kpc). Based on the [Mg/Fe] and [Fe/H] ratios, we separated the thin- from the thick-disc sequence. From analysing the Galactocentric velocity of the stars for the thin disc, we find a weak positive correlation between Vphi and [Fe/H], due to a slowly rotating Fe-poor tail. For the thick disc, a strong correlation with [Fe/H] and [Mg/Fe] is established. We have detected an inversion of the radial velocity dispersion with [Mg/Fe] for thick-disc stars with [Fe/H]<-0.1dex and [Mg/Fe]>+0.2dex. First, the velocity dispersion increases with [Mg/Fe] at all [Fe/H] ratios for the thin-disc stars, and then it decreases for the thick-disc at the highest [Mg/Fe] abundances. Similar trends are observed within the errors for the azimuthal velocity dispersion, while a continuous increase with [Mg/Fe] is observed for the vertical velocity dispersion. The velocity dispersion decrease agrees with previous measurements of the RAVE survey, although it is observed here for a greater metallicity interval and a larger spatial volume. We confirm the existence of [Mg/Fe]-rich thick-disc stars with cool kinematics in the generally turbulent context of the primitive Galactic disc. This is discussed in the framework of the different disc formation scenarios.

It is well known that newly formed planetary systems undergo processes of orbital reconfiguration and planetary migration. As a result, planets or protoplanetary objects may accrete onto the central star, being fused and mixed into its external layers. If the accreted mass is sufficiently high and the star has a sufficiently thin convective envelope, such events may result in a modification of the chemical composition of the stellar photosphere in an observable way, enhancing it with elements that were abundant in the accreted mass. The recent Gaia-ESO Survey observations of the 10-20 Myr old Gamma Velorum cluster have enabled identifying a star that is significantly enriched in iron with respect to other cluster members. In this Letter we further investigate the abundance pattern of this star, showing that its abundance anomaly is not limited to iron, but is also present in the refractory elements, whose overabundances are correlated with the condensation temperature. This finding strongly supports the hypothesis of a recent accretion of rocky material.

We investigate, using the Gaia-ESO Survey internal Data-Release 2, the properties of the double sequence of the Milky Way discs (defined chemically as the high-alpha and low-alpha populations), and discuss their compatibility with discs defined by other means such as metallicity, kinematics or positions. This investigation uses two different approaches: in velocity space for stars located in the extended Solar neighbourhood, and in chemical space for stars at different ranges of Galactocentric radii and heights from the plane. The separation we find in velocity space allows us to investigate, in a novel manner, the extent in metallicity of each of the two sequences, identifying them with the two discs, without making any assumption about the shape of their metallicity distribution functions. Then, using the separation in chemical space, we characterise the spatial variation of the slopes of the [alpha/Fe] - [Fe/H] sequences for the thick and thin discs and the way in which the relative proportions of the two discs change across the Galaxy. We find that the thick disc (high-alpha sequence), extends up to [Fe/H]~ +0.2 and the thin disc (low-alpha sequence), at least down to [Fe/H]~ -0.8. Radial and vertical gradients in alpha-abundances are found for the thin disc, with mild spatial variations in its [alpha/Fe] - [Fe/H] paths, whereas for the thick disc we do not detect any such spatial variations. The small variations in the spatial [alpha/Fe] - [Fe/H] paths of the thin disc do not allow us to distinguish between formation models of this structure. On the other hand, the lack of radial gradients and [alpha/Fe] - [Fe/H] variations for the thick disc indicate that the mechanism responsible for the mixing of the metals in the young Galaxy (e.g. radial stellar migration or turbulent gaseous disc) was more efficient before the (present) thin disc started forming.

We present results from a spectroscopic Spitzer and Herschel mid-to-far-infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines of [Ne II], H_2 S(1), CO_2 and HCN, while Herschel PACS and SPIRE spectra reveal emission from [O I], OH, and tentative detections of H_2O and high-J transitions of CO. We measure [Ne III]/[Ne II] < 0.13, which is comparable to other X-ray/EUV luminous T Tauri stars that lack jets. We use the H_2 S(1) line luminosity to estimate the gas mass in the relatively warm surface layers of the inner disk. The presence of [O I] emission suggests that CO, H_2O, and/or OH is being photodissociated, and the lack of [C I] emission suggests any excess C may be locked up in HCN, CN and other organic molecules. Modeling of silicate dust grain emission features in the mid-infrared indicates that the inner disk is composed mainly of large (r~5 um) amorphous pyroxene and olivine grains (~86% by mass) with a relatively large proportion of crystalline silicates. These results are consistent with other lines of evidence indicating that planet building is ongoing in regions of the disk within ~30 AU of the central, close binary.

We discuss the properties of H$\alpha$ emission stars across the sample of 22035 spectra from the Gaia-ESO Survey internal data release, observed with the GIRAFFE instrument and largely belonging to stars in young open clusters. Automated fits using two independent Gaussian profiles and a third component that accounts for the nebular emission allow us to discern distinct morphological types of H$\alpha$ line profiles with the introduction of a simplified classification scheme. All in all we find 3765 stars with intrinsic emission and sort their spectra into eight distinct morphological categories: single--component emission, emission blend, sharp emission peaks, double emission, P-Cygni, inverted P-Cygni, self--absorption, and emission in absorption. We have more than one observation for 1430 stars in our sample, thus allowing a quantitative discussion of the degree of variability of H$\alpha$ emission profiles, which is expected for young, active objects. We present a catalogue of stars with properties of their H$\alpha$ emission line profiles, morphological classification, analysis of variability with time and the supplementary information from the SIMBAD, VizieR, and ADS databases. The records in SIMBAD indicate the presence of H$\alpha$ emission for roughly 25% of all stars in our catalogue, while at least 305 of them have already been more thoroughly investigated according to the references in ADS. The most frequently identified morphological categories in our sample of spectra are emission blend (23%), emission in absorption (22%), and self--absorption (16%). Objects with repeated observations demonstrate that our classification into discrete categories is generally stable through time, but categories P-Cygni and self--absorption seem less stable, which is the consequence of discrete classification rules, as well as of the fundamental change in profile shape.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

Context. The inner disc, linking the thin disc with the bulge, has been somehow neglected in the past because of intrinsic difficulties in its study, due, e.g., to crowding and high extinction. Open clusters located in the inner disc are among the best tracers of its chemistry at different ages and distances. Aims. We analyse the chemical patterns of four open clusters located within 7 kpc of the Galactic Centre and of field stars to infer the properties of the inner disc with the Gaia-ESO survey idr2/3 data release. Methods. We derive the parameters of the newly observed cluster, Berkeley 81, finding an age of about 1 Gyr and a Galactocentric distance of 5.4 kpc. We construct the chemical patterns of clusters and we compare them with those of field stars in the Solar neighbourhood and in the inner-disc samples. Results. Comparing the three populations we observe that inner-disc clusters and field stars are both, on average, enhanced in [O/Fe], [Mg/Fe] and [Si/Fe]. Using the idr2/3 results of M67, we estimate the non-local thermodynamic equilibrium (NLTE) effect on the abundances of Mg and Si in giant stars. After empirically correcting for NLTE effects, we note that NGC 6705 and Be 81 still have a high [{\alpha}/Fe]. Conclusions. The location of the four open clusters and of the field population reveals that the evolution of the metallicity [Fe/H] and of [alpha/Fe] can be explained within the framework of a simple chemical evolution model: both [Fe/H] and [{\alpha}/Fe] of Trumpler 20 and of NGC 4815 are in agreement with expectations from a simple chemical evolution model. On the other hand, NGC 6705, and at a lower level Berkeley 81, have higher [{\alpha}/Fe] than expected for their ages, location in the disc, and metallicity. These differences might originate from local enrichment processes as explained in the inhomogeneous evolution framework.

The Gaia-ESO Survey has recently unveiled the complex kinematic signature of the Gamma Velorum cluster: this cluster is composed of two kinematically distinct populations (hereafter, population A and B), showing two different velocity dispersions and a relative ~2 km s^-1 radial velocity (RV) shift. In this paper, we propose that the two populations of the Gamma Velorum cluster originate from two different sub-clusters, born from the same parent molecular cloud. We investigate this possibility by means of direct-summation N-body simulations. Our scenario is able to reproduce not only the RV shift and the different velocity dispersions, but also the different centroid (~0.5 pc), the different spatial concentration and the different line-of-sight distance (~5 pc) of the two populations. The observed 1-2 Myr age difference between the two populations is also naturally explained by our scenario, in which the two sub-clusters formed in two slightly different star formation episodes. Our simulations suggest that population B is strongly supervirial, while population A is close to virial equilibrium. We discuss the implications of our models for the formation of young star clusters and OB associations in the Milky Way.

The outer (>30 AU) regions of the dusty circumstellar disk orbiting the ~2-5 Myr-old, actively accreting solar analog LkCa 15 are known to be chemically rich, and the inner disk may host a young protoplanet within its central cavity. To obtain a complete census of the brightest molecular line emission emanating from the LkCa 15 disk over the 210-270 GHz (1.4 - 1.1 mm) range, we have conducted an unbiased radio spectroscopic survey with the Institute de Radioastronomie Millimetrique (IRAM) 30 meter telescope. The survey demonstrates that, in this spectral region, the most readily detectable lines are those of CO and its isotopologues 13CO and C18O, as well as HCO+, HCN, CN, C2H, CS, and H2CO. All of these species had been previously detected in the LkCa 15 disk; however, the present survey includes the first complete coverage of the CN (2-1) and C2H (3-2) hyperfine complexes. Modeling of these emission complexes indicates that the CN and C2H either reside in the coldest regions of the disk or are subthermally excited, and that their abundances are enhanced relative to molecular clouds and young stellar object environments. These results highlight the value of unbiased single-dish line surveys in guiding future high resolution interferometric imaging of disks.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

The nearby (distance~350-400 pc), rich Vela OB2 association, includes $\gamma^2$ Velorum, one of the most massive binaries in the solar neighbourhood and an excellent laboratory for investigating the formation and early evolution of young clusters. Recent Gaia-ESO survey observations have led to the discovery of two kinematically distinct populations in the young (10-15 Myr) cluster immediately surrounding $\gamma^2$ Velorum. Here we analyse the results of Gaia-ESO survey observations of NGC 2547, a 35 Myr cluster located two degrees south of $\gamma^2$ Velorum. The radial velocity distribution of lithium-rich pre-main sequence stars shows a secondary population that is kinematically distinct from and younger than NGC 2547. The radial velocities, lithium absorption lines, and the positions in a colour-magnitude diagram of this secondary population are consistent with those of one of the components discovered around $\gamma^2$ Velorum. This result shows that there is a young, low-mass stellar population spread over at least several square degrees in the Vela OB2 association. This population could have originally been part of a cluster around $\gamma^2$ Velorum that expanded after gas expulsion or formed in a less dense environment that is spread over the whole Vela OB2 region.

We use the fundamental parameters delivered by the GES consortium in the first internal data release to select the members of $\gamma$ Vel and Cha I among the UVES and GIRAFFE spectroscopic observations. A total of 140 $\gamma$ Vel members and 74 Cha I members were studied. We calculated stellar luminosities through spectral energy distributions, while stellar masses were derived by comparison with evolutionary tracks. The spectral subtraction of low-activity and slowly rotating templates, which are rotationally broadened to match the $v\sin i$ of the targets, enabled us to measure the equivalent widths (EWs) and the fluxes in the H$\alpha$ and H$\beta$ lines. The H$\alpha$ line was also used for identifying accreting objects and for evaluating the mass accretion rate ($\dot M_{\rm acc}$). The distribution of $v\sin i$ for the members of $\gamma$ Vel displays a peak at about 10 km s$^{-1}$ with a tail toward faster rotators. There is also some indication of a different $v\sin i$ distribution for the members of its two kinematical populations. Only a handful of stars in $\gamma$ Vel display signatures of accretion, while many more accretors were detected in the younger Cha~I. Accreting and active stars occupy two different regions in a $T_{\rm eff}$-flux diagram and we propose a criterion for distinguishing them. We derive $\dot M_{\rm acc}$ in the ranges $10^{-11}$-$10^{-9} M_\odot$yr$^{-1}$ and $10^{-10}$-$10^{-7} M_\odot$yr$^{-1}$ for $\gamma$ Vel and Cha I accretors, respectively. We find less scatter in the $\dot M_{\rm acc}-M_\star$ relation derived through the H$\alpha$ EWs, when compared to the H$\alpha$ $10\%W$ diagnostics, in agreement with other authors.

The Gaia-ESO survey is a large public spectroscopic survey aimed at investigating the origin and formation history of our Galaxy by collecting spectroscopy of representative samples (about 10^5 Milky Way stars) of all Galactic stellar populations, in the field and in clusters. The survey uses globular clusters as intra- and inter-survey calibrators, deriving stellar atmospheric parameters and abundances of a significant number of stars in clusters, along with radial velocity determinations. We used precise radial velocities of a large number of stars in seven globular clusters (NGC 1851, NGC 2808, NGC 4372, NGC 4833, NGC 5927, NGC 6752, and NGC 7078) to validate pipeline results and to preliminarily investigate the cluster internal kinematics. Radial velocity measurements were extracted from FLAMES/GIRAFFE spectra processed by the survey pipeline as part of the second internal data release of data products to ESO. We complemented our sample with ESO archival data obtained with different instrument configurations. Reliable radial velocity measurements for 1513 bona fide cluster star members were obtained in total. We measured systemic rotation, estimated central velocity dispersions, and present velocity dispersion profiles of all the selected clusters, providing the first velocity dispersion curve and the first estimate of the central velocity dispersion for the cluster NGC~5927. Finally, we explore the possible link between cluster kinematics and other physical parameters. The analysis we present here demonstrates that Gaia-ESO survey data are sufficiently accurate to be used in studies of kinematics of stellar systems and stellar populations in the Milky Way.

We present the first results of the EMBLA survey (Extremely Metal-poor BuLge stars with AAOmega), aimed at finding metal-poor stars in the Milky Way bulge, where the oldest stars should now preferentially reside. EMBLA utilises SkyMapper photometry to pre-select metal-poor candidates, which are subsequently confirmed using AAOmega spectroscopy. We describe the discovery and analysis of four bulge giants with -2.72<=[Fe/H]<=-2.48, the lowest metallicity bulge stars studied with high-resolution spectroscopy to date. Using FLAMES/UVES spectra through the Gaia-ESO Survey we have derived abundances of twelve elements. Given the uncertainties, we find a chemical similarity between these bulge stars and halo stars of the same metallicity, although the abundance scatter may be larger, with some of the stars showing unusual [{\alpha}/Fe] ratios.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

Chemically inhomogeneous populations are observed in most globular clusters, but not in open clusters. Cluster mass seems to play a key role in the existence of multiple populations. Studying the chemical homogeneity of the most massive open clusters is necessary to better understand the mechanism of their formation and determine the mass limit under which clusters cannot host multiple populations. Here we studied NGC6705, that is a young and massive open cluster located towards the inner region of the Milky Way. This cluster is located inside the solar circle. This makes it an important tracer of the inner disk abundance gradient. This study makes use of BVI and ri photometry and comparisons with theoretical isochrones to derive the age of NGC6705. We study the density profile of the cluster and the mass function to infer the cluster mass. Based on abundances of the chemical elements distributed in the first internal data release of the Gaia-ESO Survey, we study elemental ratios and the chemical homogeneity of the red clump stars. Radial velocities enable us to study the rotation and internal kinematics of the cluster. The estimated ages range from 250 to 316Myr, depending on the adopted stellar model. Luminosity profiles and mass functions show strong signs of mass segregation. We derive the mass of the cluster from its luminosity function and from the kinematics, finding values between 3700 M$_{\odot}$ and 11 000 M$_{\odot}$. After selecting the cluster members from their radial velocities, we obtain a metallicity of [Fe/H]=0.10$\pm$0.06 based on 21 candidate members. Moreover, NGC6705 shows no sign of the typical correlations or anti-correlations between Al, Mg, Si, and Na, that are expected in multiple populations. This is consistent with our cluster mass estimate, which is lower than the required mass limit proposed in literature to develop multiple populations.

Context. Recent metallicity determinations in young open clusters and star-forming regions suggest that the latter may be characterized by a slightly lower metallicity than the Sun and older clusters in the solar vicinity. However, these results are based on small statistics and inhomogeneous analyses. The Gaia-ESO Survey is observing and homogeneously analyzing large samples of stars in several young clusters and star-forming regions, hence allowing us to further investigate this issue. Aims. We present a new metallicity determination of the Chamaeleon I star-forming region, based on the products distributed in the first internal release of the Gaia-ESO Survey. Methods. 48 candidate members of Chamaeleon I have been observed with the high-resolution spectrograph UVES. We use the surface gravity, lithium line equivalent width and position in the Hertzsprung-Russell diagram to confirm the cluster members and we use the iron abundance to derive the mean metallicity of the region. Results. Out of the 48 targets, we confirm 15 high probability members. Considering the metallicity measurements for 9 of them, we find that the iron abundance of Chamaeleon I is slightly subsolar with a mean value [Fe/H]=-0.08+/-0.04 dex. This result is in agreement with the metallicity determination of other nearby star-forming regions and suggests that the chemical pattern of the youngest stars in the solar neighborhood is indeed more metal-poor than the Sun. We argue that this evidence may be related to the chemical distribution of the Gould Belt that contains most of the nearby star-forming regions and young clusters.

Knowledge of the abundance distribution of star forming regions and young clusters is critical to investigate a variety of issues, from triggered star formation and chemical enrichment by nearby supernova explosions to the ability to form planetary systems.In spite of this, detailed abundance studies are currently available for relatively few regions. In this context, we present the analysis of the metallicity of the Gamma Velorum cluster, based on the products distributed in the first internal release of the Gaia-ESO Survey. The Gamma Velorum candidate members have been observed with FLAMES, using both UVES and Giraffe, depending on the target brightness and spectral type. In order to derive a solid metallicity determination for the cluster, membership of the observed stars must be first assessed. To this aim, we use several membership criteria including radial velocities, surface gravity estimates, and the detection of the photospheric lithium line. Out of the 80 targets observed with UVES, we identify 14 high-probability members. We find that the metallicity of the cluster is slightly subsolar, with a mean [Fe/H]=-0.057+/-0.018 dex. Although J08095427-4721419 is one of the high-probability members, its metallicity is significantly larger than the cluster average. We speculate about its origin as the result of recent accretion episodes of rocky bodies of ~60 M_Sun hydrogen-depleted material from the circumstellar disk.

We study spectral diagnostics available from optical spectra with R=17000 obtained with the VLT/Giraffe HR15n setup, using observations from the Gaia-ESO Survey, on the gamma Vel young cluster, in order to determine the fundamental parameters of these stars. We define a set of spectroscopic indices, sampling TiO bands, H-alpha core and wings, and many temperature- and gravity-sensitive lines. Combined indices tau (gamma) are also defined as Teff (log g) indicators over a wide spectral-type range. H-alpha emission-line indices are also chromospheric activity or accretion indicators. A metallicity-sensitive index is also defined. These indices enable us to find a clear difference between gravities of main-sequence and pre-main-sequence stars (as well as giant stars): the (gamma,tau) diagram is thus argued to be a promising distance-independent age measurement tool for young clusters. Our indices were quantitatively calibrated by means of photometry and literature reference spectra (from UVES-POP and ELODIE 3.1 Libraries, and other datasets), over a wide range of stellar parameters. The ability of our indices to select peculiar, or other rare classes of stars is also established. Finally, our gravity determinations support the existence of an older pre-main-sequence population in the gamma Vel sky region, in agreement with the lithium depletion pattern of those stars (abridged).

We study the relationship between age, metallicity, and alpha-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpc to 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages > 9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more alpha-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

(Abridged) We have used the atmospheric parameters, [alpha/Fe] abundances and radial velocities, determined from the Gaia-ESO Survey GIRAFFE spectra of FGK-type stars (iDR1), to provide a chemo-kinematical characterisation of the disc stellar populations. We focuss on a subsample of 1016 stars with high quality parameters, covering the volume |Z|<4.5kpc and R in the range 2-13kpc. We have identified a thin to thick disc separation in the [alpha/Fe] vs [M/H] plane, thanks to the presence of a low-density region in the number density distribution. The thick disc stars seem to lie in progressively thinner layers above the Galactic plane, as metallicity increases and [alpha/Fe] decreases. The thin disc population presents a constant value of the mean distance to the plane at all metallicities. Our data confirm the already known correlations between V_phi and [M/H] for the two discs. For the thick disc sequence, a study of the possible contamination by thin disc stars suggests a gradient up to 64km/s/dex. The distributions of V_phi, V_Z, and orbital parameters are analysed for the chemically separated samples. Concerning the gradients with galactocentric radius, we find for the thin disc a flat behaviour of V_phi, a [M/H] gradient of -0.058dex/kpc and a small positive [alpha/Fe] gradient. For the thick disc, flat gradients in [M/H] and [alpha/Fe] are derived. Our chemo-kinematical analysis suggests a picture in which the thick disc seems to have experienced a settling process, during which its rotation increased progressively, and, possibly, the V_phi dispersion decreased. At [M/H]-0.25dex and [alpha/Fe]0.1dex, the mean characteristics of the thick disc in distance to the Galactic plane, V_phi, V_phi dispersion and eccentricity agree with those of the thin disc stars, suggesting a possible connection between these populations at a certain epoch of the disc evolution.

NGC 4815 is a populous ~500 Myr open cluster at a Galactocentric radius of 7 kpc observed in the first six months of the Gaia-ESO Survey. Located in the inner Galactic disk, NGC 4815 is an important tracer of the abundance gradient, where few intermediate age open clusters are found. We use the survey derived radial velocities, stellar atmospheric parameters, metallicity, and elemental abundances for stars targeted as potential members of this cluster to carry out an analysis of cluster properties. The radial velocity distribution of stars in the cluster field is used to define the cluster systemic velocity and derive likely cluster membership for stars observed by the Gaia-ESO Survey. We investigate the distributions of Fe and Fe-peak elements, alpha-elements, and the light elements Na and Al and characterize the cluster's internal chemical homogeneity comparing it to the properties of radial velocity non-member stars. Utilizing these cluster properties, the cluster color-magnitude diagram is analyzed and theoretical isochrones are fit to derive cluster reddening, distance, and age. NGC 4815 is found to have a mean [Fe/H]=+0.03. Elemental abundances of cluster members show typically small internal variation, with internal dispersions of ~0.05 dex. [Ca/Fe] and [Si/Fe] show solar ratios, but [Mg/Fe] is moderately enhanced, while [Ti/Fe] appears slightly deficient. As with many open clusters, [Na/Fe] and [Al/Fe] are enhanced, [Na/Fe] significantly so, although the role of internal mixing and the assumption of LTE in the analysis remain to be investigated. From isochrone fits to color-magnitude diagrams, we find a cluster age of 0.5 to 0.63 Gyr, a reddening of E(B-V) = 0.59 to 0.65, and a true distance modulus (m-M) = 11.95 to 12.20, depending on the choice of theoretical models, leading to a Galactocentric distance of 6.9 kpc.

The Gaia-ESO Survey is a large public spectroscopic survey that aims to derive radial velocities and fundamental parameters of about 10^5 Milky Way stars in the field and in clusters. Observations are carried out with the multi-object optical spectrograph FLAMES, using simultaneously the medium resolution (R~20,000) GIRAFFE spectrograph and the high resolution (R~47,000) UVES spectrograph. In this paper, we describe the methods and the software used for the data reduction, the derivation of the radial velocities, and the quality control of the FLAMES-UVES spectra. Data reduction has been performed using a workflow specifically developed for this project. This workflow runs the ESO public pipeline optimizing the data reduction for the Gaia-ESO Survey, performs automatically sky subtraction, barycentric correction and normalisation, and calculates radial velocities and a first guess of the rotational velocities. The quality control is performed using the output parameters from the ESO pipeline, by a visual inspection of the spectra and by the analysis of the signal-to-noise ratio of the spectra. Using the observations of the first 18 months, specifically targets observed multiple times at different epochs, stars observed with both GIRAFFE and UVES, and observations of radial velocity standards, we estimated the precision and the accuracy of the radial velocities. The statistical error on the radial velocities is sigma~0.4 km s^-1 and is mainly due to uncertainties in the zero point of the wavelength calibration. However, we found a systematic bias with respect to the GIRAFFE spectra (~0.9 km s^-1) and to the radial velocities of the standard stars (~0.5 kms^-1) retrieved from the literature. This bias will be corrected in the future data releases, when a common zero point for all the setups and instruments used for the survey will be established.

Context: A key science goal of the Gaia-ESO survey (GES) is to use the kinematics of low-mass stars in young clusters to probe their dynamical histories and how they populate the field as they become unbound. The clustering of low-mass stars around the massive W-R binary gamma^2 Velorum was one of the first GES targets. Aims: To empirically determine the radial velocity (RV) precision of GES data, construct a kinematically unbiased sample of cluster members and characterise their dynamical state. Methods: Targets were selected from colour-magnitude diagrams and intermediate resolution spectroscopy used to derive RVs and assess membership from the strength of the Li6708A line. The RV distribution was analysed using a maximum likelihood technique that accounts for unresolved binaries. Results: The GES RV precision is about 0.25km/s and sufficient to resolve velocity structure in the low-mass population around gamma^2 Vel. The structure is well fitted by two kinematic components with roughly equal numbers of stars; the first has an intrinsic dispersion of 0.34+/-0.16km/s, consistent with virial equilibrium. The second has a broader dispersion of 1.60+/-0.37km/s and is offset from the first by ~2km/s. The first population is older by 1-2Myr based on a greater level of Li depletion seen among its M-stars and is probably more centrally concentrated around gamma^2 Vel. Conclusions: We consider several formation scenarios, concluding that the two kinematic components are a bound remnant of the original, denser cluster that formed gamma^2 Vel, and a dispersed population from the wider Vela OB2 association, of which gamma^2 Vel is the most massive member. The apparent youth of gamma^2 Vel compared to the older (>=10Myr) low-mass population surrounding it suggests a scenario where the massive binary formed in a clustered environment after the formation of the bulk of the low-mass stars.[abridged]

Open clusters are key tools to study the spatial distribution of abundances in the disk and their evolution with time. Using the first release of stellar parameters and abundances of the Gaia-ESO Survey, we analyse the chemical properties of stars in three old/intermediate-age open clusters, namely NGC 6705, NGC 4815, and Trumpler 20, all located in the inner part of the Galactic disk at Galactocentric radius R$_{GC}\sim$7 kpc, aiming at proving their homogeneity and at comparing them with the field population. We study the abundance ratios of elements belonging to two different nucleosynthetic channels: $\alpha$-elements and iron-peak elements. The main results can be summarised as follows: i) cluster members are chemically homogeneous within 3-$\sigma$ in all analysed elements; ii) the three clusters have comparable [El/Fe] patters within $\sim$1-$\sigma$, but they differ in their global metal content [El/H], with NGC 4815 having the lowest metallicity. Their [El/Fe] ratios show differences and analogies with those of the field population, both in the solar neighbourhood and in the bulge/inner disk; iii) comparing the abundance ratios with the results of two chemical evolution models and with field star abundance distributions, we find that the abundance ratios of Mg, Ni, Ca in NGC 6705 might require an inner birthplace, implying a subsequent variation of its R$_{GC}$ during its lifetime, consistent with previous orbit determination. The full dataset of the Gaia-ESO Survey will be a superlative tool to constrain the chemical evolution of our Galaxy by disentangling different formation and evolution scenarios.

Trumpler 20 is an old open cluster (OC) located toward the Galactic centre, at about 3 kpc from the Sun and $\sim$7 kpc from the Galactic centre. Its position makes this cluster particularly interesting in the framework of the chemical properties of the Galactic disc because very few old OCs reside in the inner part of the disc. For this reason it has been selected as a cluster target of the Gaia-ESO Survey, and spectra of many stars in the main sequence and red clump phases are now available. Moreover, although it has been studied by several authors in the past, no consensus on the evolutionary status of Tr 20 has been reached. The heavy contamination of field stars (the line of sight of Tr 20 crosses the Carina spiral arm) complicates a correct interpretation. Another interesting aspect of the cluster is that it shows a broadened main-sequence turn-off and a prominent and extended red-clump, characteristics that are not easily explained by classical evolutionary models. Exploiting both spectroscopic information from the Gaia-ESO Survey (and the ESO archive) and literature photometry, we obtain a detailed and accurate analysis of the properties of the cluster. We make use of the first accurate metallicity measurement ever obtained from several spectra of red clump stars, and of cluster membership determination using radial velocities. According to the evolutionary models adopted, we find that Tr 20 has an age in the range 1.35-1.66 Gyr, an average reddening $E(B-V)$ in the range 0.31-0.35 mag, and a distance modulus $(m-M)_0$ between 12.64 and 12.72 mag. The spectroscopic metallicity is [Fe/H]=+0.17 dex. We discuss the structural properties of the object and constrain possible hypotheses for its broadened upper main sequence by estimating the effect of differential reddening and its extended red clump.

We seek to establish the presence and properties of gas in the circumstellar disk orbiting T Cha, a nearby (d~110 pc), relatively evolved (age ~5-7 Myr) yet actively accreting 1.5 Msun T Tauri star. We used the APEX 12 m radiotelescope to search for submillimeter molecular emission from the T Cha disk, and we reanalyzed archival XMM-Newton spectroscopy of T Cha to ascertain the intervening absorption due to disk gas along the line of sight to the star (N_H). We detected submillimeter rotational transitions of 12CO, 13CO, HCN, CN and HCO+ from the T Cha disk. The 12CO line appears to display a double-peaked line profile indicative of Keplerian rotation. Analysis of the CO emission line data indicates that the disk around T Cha has a mass (M_disk,H_2 = 80 M_earth) similar to, but more compact (R_disk, CO~80 AU) than, other nearby, evolved molecular disks (e.g. V4046 Sgr, TW Hya, MP Mus) in which cold molecular gas has been previously detected. The HCO+/13CO and HCN/13CO, line ratios measured for T Cha appear similar to those of other evolved circumstellar disks (i.e. TW Hya and V4046 Sgr), while the CN/13CO ratio appears somewhat weaker. Analysis of the XMM-Newton data shows that the atomic absorption $N_H$ toward T Cha is 1-2 orders of magnitude larger than toward the other nearby T Tauri with evolved disks. Furthermore, the ratio between atomic absorption and optical extinction N_H/A_V toward T Cha is higher than the typical value observed for the interstellar medium and young stellar objects in the Orion Nebula Cluster. This may suggest that the fraction of metals in the disk gas is higher than in the interstellar medium. Our results confirm that pre-main sequence stars older than ~5 Myr, when accreting, retain cold molecular disks, and that those relatively evolved disks display similar physical and chemical properties.

We present a preliminary 3D potential field extrapolation model of the joint magnetosphere of the close accreting PMS binary V4046 Sgr. The model is derived from magnetic maps obtained as part of a coordinated optical and X-ray observing program.

V4046 Sagittarii AB is a close short-period classical T Tauri binary. It is a circularised and synchronised system accreting from a circumbinary disk. In 2009 it was observed as part of a coordinated program involving near-simultaneous spectropolarimetric observations with ESPaDOnS at the Canada-France-Hawai'i Telescope and high-resolution X-ray observations with XMM-Newton. Magnetic maps of each star were derived from Zeeman-Doppler imaging. After briefly highlighting the most significant observational findings, we present a preliminary 3D model of the binary magnetosphere constructed from the magnetic maps using a newly developed binary magnetic field extrapolation code. The large-scale fields (the dipole components) of both stars are highly tilted with respect to their rotation axes, and their magnetic fields are linked.

We report initial results from a quasi-simultaneous X-ray/optical observing campaign targeting V4046 Sgr, a close, synchronous-rotating classical T Tauri star (CTTS) binary in which both components are actively accreting. V4046 Sgr is a strong X-ray source, with the X-rays mainly arising from high-density (n_e ~ 10^(11-12) cm^(-3)) plasma at temperatures of 3-4 MK. Our multiwavelength campaign aims to simultaneously constrain the properties of this X-ray emitting plasma, the large scale magnetic field, and the accretion geometry. In this paper, we present key results obtained via time-resolved X-ray grating spectra, gathered in a 360 ks XMM-Newton observation that covered 2.2 system rotations. We find that the emission lines produced by this high-density plasma display periodic flux variations with a measured period, 1.22+/-0.01 d, that is precisely half that of the binary star system (2.42 d). The observed rotational modulation can be explained assuming that the high-density plasma occupies small portions of the stellar surfaces, corotating with the stars, and that the high-density plasma is not azimuthally symmetrically distributed with respect to the rotational axis of each star. These results strongly support models in which high-density, X-ray-emitting CTTS plasma is material heated in accretion shocks, located at the base of accretion flows tied to the system by magnetic field lines.

The early-K star T Cha, a member of the relatively nearby (D ~ 100 pc) epsilon Cha Association, is a relatively "old" (age ~7 Myr) T Tauri star that is still sporadically accreting from an orbiting disk whose inner regions are evidently now being cleared by a close, substellar companion. We report the identification, via analysis of proper motions, serendipitous X-ray imaging spectroscopy, and followup optical spectroscopy, of a new member of the epsilon Cha Association that is very likely a low-mass companion to T Cha at a projected separation of ~38 kAU. The combined X-ray and optical spectroscopy data indicate that the companion, T Cha B (= 2M1155-79), is a weak-lined T Tauri star (wTTS) of spectral type M3 and age ~<10 Myr. The serendipitous X-ray (XMM-Newton) observation of T Cha B, which targeted T Cha, also yields serendipitous detections of two background wTTS in the Chamaeleon cloud complex, including one newly discovered, low-mass member of the Cha cloud pre-MS population. T Cha becomes the third prominent example of a nearby, "old" yet still actively accreting, K-type pre-MS star/disk system (the others being TW Hya and V4046 Sgr) to feature a low-mass companion at very large (12-40 kAU) separation, suggesting that such wide-separation companions may affect the conditions and timescales for planet formation around solar-mass stars.

Constraining the spatial and thermal structure of the gaseous component of circumstellar disks is crucial to understand star and planet formation. Models predict that the [Ne II] line at 12.81 {\mu}m detected in young stellar objects with Spitzer traces disk gas and its response to high energy radiation, but such [Ne II] emission may also originate in shocks within powerful outflows. To distinguish between these potential origins for mid-infrared [Ne II] emission and to constrain disk models, we observed 32 young stellar objects using the high resolution (R~30000) mid-infrared spectrograph VISIR at the VLT. We detected the 12.81 {\mu}m [Ne II] line in 12 objects, tripling the number of detections of this line in young stellar objects with high spatial and spectral resolution spectrographs. We obtain the following main results: a) In Class I objects the [Ne II] emission observed from Spitzer is mainly due to gas at a distance of more than 20-40 AU from the star, where neon is, most likely, ionized by shocks due to protostellar outflows. b) In transition and pre-transition disks, most of the emission is confined to the inner disk, within 20-40 AU from the central star. c) Detailed analysis of line profiles indicates that, in transition and pre-transition disks, the line is slightly blue-shifted (2-12 km s{^-1}) with respect to the stellar velocity, and the line width is directly correlated with the disk inclination, as expected if the emission is due to a disk wind. d) Models of EUV/X-ray irradiated disks reproduce well the observed relation between the line width and the disk inclination, but underestimate the blue-shift of the line.

X-shooter is the first 2nd generation instrument of the ESO Very Large Telescope(VLT). It is a very efficient, single-target, intermediate-resolution spectrograph that was installed at the Cassegrain focus of UT2 in 2009. The instrument covers, in a single exposure, the spectral range from 300 to 2500 nm. It is designed to maximize the sensitivity in this spectral range through dichroic splitting in three arms with optimized optics, coatings, dispersive elements and detectors. It operates at intermediate spectral resolution (R~4,000 - 17,000, depending on wavelength and slit width) with fixed echelle spectral format (prism cross-dispersers) in the three arms. It includes a 1.8"x4" Integral Field Unit as an alternative to the 11" long slits. A dedicated data reduction package delivers fully calibrated two-dimensional and extracted spectra over the full wavelength range. We describe the main characteristics of the instrument and present its performance as measured during commissioning, science verification and the first months of science operations.

We explore the possibility that GSC 07396-00759 (spectral type M1e) is a widely separated (~2.82', or projected separation ~12,350 AU) companion to the "old" (age ~12 Myr) classical T Tauri binary system V4046 Sgr AB, as suggested by the proximity and similar space motions of the two systems. If the two systems are equidistant and coeval, then GSC 07396--00759, like V4046 Sgr AB, must be a spectroscopic binary with nearly equal-mass components, and V4046 Sgr must be at least ~8 Myr old. Analysis of a serendipitous Chandra X-ray gratings spectrum and light curve as well as XMM-Newton light curves and CCD spectra of GSC 07396-00759 obtained during long exposures targeting V4046 Sgr AB reveals a relatively hard (T_X ~ 10^7 K) X-ray spectrum, strong flaring, and relatively low-density plasma. These X-ray characteristics of GCS 07396--00759 are indicative of a high level of coronal activity, consistent with its apparent weak-lined T Tauri star status. Interactions between V4046 Sgr AB and GCS 07396-00759 when the two systems were more closely bound may be responsible for (a) their dissolution ~10^6 yr ago, (b) the present tight, circular orbit of V4046 Sgr AB, and (c) the persistence of the gaseous circumbinary disk still orbiting V4046 Sgr AB.

We report here the first results of a multi-wavelength campaign focussing on magnetospheric accretion processes within the close binary system V4046 Sgr, hosting two partly-convective classical T Tauri stars of masses ~0.9 Msun and age ~12 Myr. In this paper, we present time-resolved spectropolarimetric observations collected in 2009 September with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) and covering a full span of 7d or ~2.5 orbital/rotational cycles of V4046 Sgr. Small circularly polarised Zeeman signatures are detected in the photospheric absorption lines but not in the accretion-powered emission lines of V4046 Sgr, thereby demonstrating that both system components host large-scale magnetic fields weaker and more complex than those of younger, fully-convective cTTSs of only a few Myr and similar masses. Applying our tomographic imaging tools to the collected data set, we reconstruct maps of the large-scale magnetic field, photospheric brightness and accretion-powered emission at the surfaces of both stars of V4046 Sgr. We find that these fields include significant toroidal components, and that their poloidal components are mostly non-axisymmetric with a dipolar component of 50-100G strongly tilted with respect to the rotation axis; given the similarity with fields of partly-convective main-sequence stars of similar masses and rotation periods, we conclude that these fields are most likely generated by dynamo processes. We also find that both stars in the system show cool spots close to the pole and extended regions of low-contrast, accretion-powered emission; it suggests that mass accretion is likely distributed rather than confined in well defined high-contrast accretion spots, in agreement with the derived magnetic field complexity.

We studied the X-ray properties of the young (~1-8 Myr) open cluster around the O8III star lambda Ori using XMM-Newton. We detected 167 X-ray sources, of which 58 are identified with known cluster members, from massive stars down to low-mass stars with spectral types M5.5. Another 23 sources were identified with new possible photometric candidates. Late-type stars have a median log LX/Lbol ~ -3.3, close to the saturation limit. Variability was observed in 35% of late-type members, including six flaring sources. The emission from the central hot star lambda Ori is dominated by plasma at 0.2-0.3 keV, with a weaker component at 0.7 keV, consistently with a wind origin. The coronae of late-type stars can be described by two plasma components with temperatures T1~0.3-0.8 keV and T2~0.8-3 keV, and subsolar abundances Z~0.1-0.3 Zsun, similar to what is found in other star-forming regions and associations. No significant difference was observed between stars with and without circumstellar discs, although the smallness of the sample of stars with discs and accretion does not allow us to draw definitive conclusions. The X-ray properties of lambda Ori late-type stars are comparable to those of the coeval sigma Ori cluster, suggesting that stellar activity in lambda Ori has not been significantly affected by the different ambient environment.

High resolution X-ray spectroscopy has revealed soft X-rays from high density plasma in Classical T-Tauri stars (CTTSs), probably arising from the accretion shock region. However, the mass accretion rates derived from the X-ray observations are consistently lower than those derived from UV/optical/NIR studies. We aim to test the hypothesis that the high density soft X-ray emission is from accretion by analysing optical accretion tracers from an X-ray selected sample of CTTSs in a homogeneous manner. We analyse optical spectra of a sample of CTTSs and calculate the accretion rates based on measuring optical emission lines. These are then compared to the accretion rates derived from the X-ray spectroscopy. We find that, for each CTTS in our sample, the different optical tracers predict mass accretion rates that agree within the errors, albeit with a spread of ~1 order of magnitude. Typically, mass accretion rates derived from Halpha and HeI 5876 Ang are larger than those derived from Hbeta, Hgamma and OI. When comparisons of the optical mass accretion rates are made to the X-ray derived mass accretion rates, we find that: a) the latter are always lower (but by varying amounts); b) the latter range within a factor of ~2 around 2x10^{-10} M_odot yr^{-1}, despite the fact that the former span a range of ~3 orders of magnitude. We suggest that the systematic underestimation of the X-ray derived mass accretion rates could depend on the density distribution inside the accretion streams, where the densest part of the stream is not visible in the X-ray band because of the absorption by the stellar atmosphere. We also suggest that a non-negligible optical depth of X-ray emission lines produced by post-shock accreting plasma may explain the almost constant mass accretion rates derived in X-rays if the effect is larger in stars with larger optical mass accretion rates.

We have used the Atacama Pathfinder Experiment 12 m telescope to detect circumstellar CO emission from MP Mus (K1 IVe), a nearby (D ~ 100 pc), actively accreting, ~7 Myr-old pre-main sequence (pre-MS) star. The CO emission line profile measured for MP Mus is indicative of an orbiting disk with radius ~120 AU, assuming the central star mass is 1.2 solar masses and the disk inclination is ~30 degrees, and the inferred disk molecular gas mass is ~3 Earth masses. MP Mus thereby joins TW Hya and V4046 Sgr as the only late-type (low-mass), pre-MS star systems within ~100 pc of Earth that are known to retain orbiting, molecular disks. We also report the nondetection (with the Institut de Radio Astronomie Millimetrique 30 m telescope) of CO emission from another ten nearby (D ~ 100 pc or less), dusty, young (age ~10-100 Myr) field stars of spectral type A-G. We discuss the implications of these results for the timescales for stellar and Jovian planet accretion from, and dissipation of, molecular disks around young stars.

Context. High resolution X-ray observations of classical T Tauri stars (CTTSs) show a soft X-ray excess due to high density plasma (n_e=10^11-10^13 cm^-3). This emission has been attributed to shock-heated accreting material impacting onto the stellar surface. Aims. We investigate the observability of the shock-heated accreting material in the X-ray band as a function of the accretion stream properties (velocity, density, and metal abundance) in the case of plasma-beta<<1 in the post-shock zone. Methods. We use a 1-D hydrodynamic model describing the impact of an accretion stream onto the chromosphere, including the effects of radiative cooling, gravity and thermal conduction. We explore the space of relevant parameters and synthesize from the model results the X-ray emission in the [0.5-8.0] keV band and in the resonance lines of O VII (21.60 Ang) and Ne IX (13.45 Ang), taking into account the absorption from the chromosphere. Results. The accretion stream properties influence the temperature and the stand-off height of the shocked slab and its sinking in the chromosphere, determining the observability of the shocked plasma. Our model predicts that X-ray observations preferentially detect emission from low density and high velocity shocked accretion streams due to the large absorption of dense post-shock plasma. In all the cases examined, the post-shock zone exhibits quasi-periodic oscillations due to thermal instabilities, but in the case of inhomogeneous streams and beta<<1, the shock oscillations are hardly detectable. Conclusions. We suggest that, if accretion streams are inhomogeneous, the selection effect introduced by the absorption on observable plasma components may explain the discrepancy between the accretion rate measured by optical and X-ray data as well as the different densities measured using different He-like triplets in the X-ray band.

AIMS. We investigate the stability and dynamics of accretion shocks in CTTSs, considering the case of beta >= 1 in the post-shock region. In these cases the 1D approximation is not valid and a multi-dimensional MHD approach is necessary. METHODS. We model an accretion stream propagating through the atmosphere of a CTTS and impacting onto its chromosphere, by performing 2D axisymmetric MHD simulations. The model takes into account the stellar magnetic field, the gravity, the radiative cooling, and the thermal conduction (including the effects of heat flux saturation). RESULTS. The dynamics and stability of the accretion shock strongly depends on the plasma beta. In the case of shocks with beta > 10, violent outflows of shock-heated material (and possibly MHD waves) are generated at the base of the accretion column and strongly perturb the surrounding stellar atmosphere and the accretion column itself (modifying, therefore, the dynamics of the shock). In shocks with beta ~ 1, the post-shock region is efficiently confined by the magnetic field. The shock oscillations induced by cooling instability are strongly influenced by beta: for beta > 10, the oscillations may be rapidly dumped by the magnetic field, approaching a quasi-stationary state, or may be chaotic with no obvious periodicity due to perturbation of the stream induced by the post-shock plasma itself; for beta ~ 1 the oscillations are quasi-periodic, although their amplitude is smaller and the frequency higher than those predicted by 1D models.

Determining the metal content of low-mass members of young associations provides a tool that addresses different issues, such as triggered star formation or the link between the metal-rich nature of planet-host stars and the early phases of planet formation. The Orion complex is a well known example of possible triggered star formation and is known to host a rich variety of proto-planetary disks around its low-mass stars. Available metallicity measurements yield discrepant results. We analyzed FLAMES/UVES and Giraffe spectra of low-mass members of three groups/clusters belonging to the Orion association. Our goal is the homogeneous determination of the metallicity of the sample stars, which allows us to look for [Fe/H] differences between the three regions and for the possible presence of metal-rich stars. Nine members of the ONC and one star each in the $\lambda$ Ori cluster and OB1b subgroup were analyzed. After the veiling determination, we retrieved the metallicity by means of equivalent widths and/or spectral synthesis using MOOG. We obtain an average metallicity for the ONC [Fe/H]=-0.01\pm 0.04. No metal-rich stars were detected and the dispersion within our sample is consistent with measurement uncertainties. The metallicity of the $\lambda$ Ori member is also solar, while the OB1b star has an [Fe/H] significantly below the ONC average. If confirmed by additional [Fe/H] determinations in the OB1b subgroup, this result would support the triggered star formation and the self-enrichment scenario for the Orion complex.

The photoevaporation of circumstellar disks is a powerful process in the disk dissipation at the origin of the Orion proplyds. This Letter reports the first detection of a photoevaporating disk in the final but long-lasting phase of its evolution. The disk is associated to a low-mass T Tauri member of the sigma Orionis Cluster. It is characterized by a very low (if any) accretion rate and by a tenuous (Mloss ~ 10^{-9} Msun/yr) photoevaporation wind, which is unambiguously detected in the optical spectrum of the object. The wind emits strong forbidden lines of [SII] and [NII] because the low-mass star is close to a powerful source of ionizing photons, the O9.5 star sigma Ori.

High spectral resolution X-ray observations of CTTSs demonstrate the presence of plasma at T~2-3X10^6 K and n_e~10^11-10^13 cm^-3, unobserved in non-accreting stars. Stationary models suggest that this emission is due to shock-heated accreting material, but they do not allow to analyze the stability of such material and its position in the stellar atmosphere. We investigate the dynamics and the stability of shock-heated accreting material in CTTSs and the role of the stellar chromosphere in determining the position and the thickness of the shocked region. We perform 1-D HD simulations of the impact of the accretion flow onto chromosphere of a CTTS, including the effects of gravity, radiative losses from optically thin plasma, thermal conduction and a well tested detailed model of the stellar chromosphere. Here we present the results of a simulation based on the parameters of the CTTS MP Mus. We find that the accretion shock generates an hot slab of material above the chromosphere with a maximum thickness of 1.8X10^9 cm, density n_e~10^11-10^2 cm^-3, temperature T~3X10^6 K and uniform pressure equal to the ram pressure of the accretion flow (~450 dyn cm^-2). The base of the shocked region penetrates the chromosphere and stays where the ram pressure is equal to the thermal pressure. The system evolves with quasi-periodic instabilities of the material in the slab leading to cyclic disappearance and re-formation of the slab. For an accretion rate of ~10^-10 M_sun yr^-1, the shocked region emits a time-averaged X-ray luminosity L_X~7X10^29 erg s^-1, which is comparable to the X-ray luminosity observed in CTTSs of the same mass. Furthermore, the X-ray spectrum synthesized from the simulation matches in detail all the main features of the O VIII and O VII lines of the star MP Mus.

We have analyzed optical spectra of 25 X-ray sources identified as potential new members of the Taurus molecular cloud (TMC), in order to confirm their membership in this SFR. Fifty-seven candidates were previously selected among the X-ray sources in the XEST survey, having a 2MASS counterpart compatible with a PMS star based on color-magnitude and color-color diagrams. We obtained high-resolution optical spectra for 7 of these candidates with the SARG spectrograph at the TNG telescope, which were used to search for Li absorption and to measure the Ha line and the radial and rotational velocities; 18 low-resolution optical spectra obtained with DOLORES for other candidate members were used for spectral classification, for Ha measurements, and to assess membership together with IR color-color and color-magnitude diagrams and additional information from the X-ray data. We found that 3 sources show Li absorption, with equivalent widths of ~500 mA, broad spectral line profiles, indicating v sin i ~20-40 km/s, radial velocities consistent with those for known members, and Ha emission. Two of them are classified as new WTTSs, while the EW (~ -9 Ang) of the Ha line and its broad asymmetric profile clearly indicate that the third star (XEST-26-062) is a CTTS. Fourteen sources observed with DOLORES are M-type stars. Fifteen sources show Ha emission; 6 of them have spectra that indicate surface gravity lower than in MS stars, and their de-reddened positions in IR color-magnitude diagrams are consistent with their derived spectral type and with PMS models at the distance of the TMC. The K-type star XEST-11-078 is confirmed as a new member from the strength of its Ha emission line. Overall, we confirm membership to the TMC for 10 out of 25 X-ray sources observed in the optical. Three sources remain uncertain.

Aims. We performed a detailed membership selection and studied the accretion properties of low-mass stars in the two apparently very similar young (1-10 Myr) clusters sigma Ori and lambda Ori. Methods. We observed 98 and 49 low-mass (0.2-1.0 M_sun) stars in sigma Ori and lambda Ori respectively, using the multi-object optical spectrograph FLAMES at the VLT, with the high-resolution (R=17,000) HR15N grating (6470-6790 AA). We used radial velocities, Li and Halpha to establish cluster membership and Halpha and other optical emission lines to analyze the accretion properties of members. Results. We identified 65 and 45 members of the sigma Ori and lambda Ori clusters, respectively and discovered 16 new candidate binary systems. We also measured rotational broadening for 20 stars and estimated the mass accretion rates in 25 stars of the sigma Ori cluster, finding values between 10^-11 and 10^-7.7 M_sun yr^-1 and in 4 stars of the lambda Ori cluster, finding values between 10^-11 and 10^-10.1 M_sun yr-1. Comparing our results with the infrared photometry obtained by the Spitzer satellite, we find that the fraction of stars with disks and the fraction of active disks is larger in the sigma Ori cluster (52+-9% and 78+-16%) than in lambda Ori (28+-8% and 40+-20%) Conclusions. The different disk and accretion properties of the two clusters could be due either to the effect of the high-mass stars and the supernova explosion in the lambda Ori cluster or to different ages of the cluster populations. Further observations are required to draw a definitive conclusion.

Accretion disks around young stars evolve in time with time scales of few million years. We present here a study of the accretion properties of a sample of 35 stars in the ~3 million year old star-forming region sigma Ori. Of these, 31 are objects with evidence of disks, based on their IR excess emission. We use near-IR hydrogen recombination lines (Pa_gamma) to measure their mass accretion rate. We find that the accretion rates are significantly lower in sigma Ori than in younger regions, such as rho-Oph, consistently with viscous disk evolution. The He I 1.083 micron line is detected (either in absorption or in emission) in 72% of the stars with disks, providing evidence of accretion-powered activity also in very low accretors, where other accretion indicators dissapear.

Context: Knowledge of the age distribution of stars in young clusters and associations is crucial to constrain models of star formation. HR diagrams of different young clusters and associations suggest the presence of age spreads, but the influence of errors on the derived ages is still largely debated. Determination of lithium abundances in low-mass stars represents an alternative and robust way to infer stellar ages. Aims: We measured lithium in a sample of low mass members of the young (4-5 Myr) sigma Ori cluster with the main goal of investigating its star formation history. Methods: Using the FLAMES multi-object spectrograph on VLT/UT2, we obtained spectra of 98 candidate cluster members. The spectra were used to determine radial velocities, to infer the presence of Halpha emission, and to measure the strength of the Li 670.8 nm absorption line. Results: Using radial velocities, Halpha and Li, together with information on X-ray emission, we identified 59 high probability cluster members. Three of them show severe Li depletion. The nuclear ages inferred for these highly depleted stars exceed 10-15 Myr; for two of them these values are in good agreement with the isochronal age, while for the third star the nuclear age exceeds the isochronal one.