Search

We present the most comprehensive and deepest X-ray study to date of the properties of the richest Wolf-Rayet (WR) population observed in a single stellar cluster, Westerlund 1 (Wd1). This work is based on 36 Chandra observations obtained from the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, plus 8 archival Chandra observations. The overall exposure depth (~1.1 Ms) and baseline of the EWOCS observations extending over more than one year enable us to perform a detailed photometric, colour, and spectral analysis, as well as to search for short- and long-term periodicity. In X-rays, we detect 20 out of the 24 known Wolf-Rayet stars in Wd1 down to an observed luminosity of ~7$\times10^{29}$erg s$^{-1}$ (assuming a distance of 4.23 kpc to Wd1), with 8 WR stars being detected in X-rays for the first time. Nine stars show clear evidence of variability over the year-long baseline, with clear signs of periodicity. The X-ray colours and spectral analysis reveal that the vast majority of the WR stars are hard X-ray sources (kT$\geq$2.0keV). The Fe XXV emission line at ~6.7 keV, which commonly originates from the wind-wind collision zone in binary systems, is detected for the first time in the spectra of 17 WR stars in Wd1. In addition the ~6.4 keV fluorescent line is observed in the spectra of three stars, indicating that dense cold material coexists with the hot gas in these systems. Overall, our X-ray results alone suggest a very high binary fraction ($\geq$80%) for the WR star population in Wd1. When combining our results with properties of the WR population from other wavelengths, we estimate a binary fraction of $\geq$92%, which could even reach unity. This suggests that either all the most massive stars are found in binary systems within Wd1, or that binarity is essential for the formation of such a rich population of WR stars., Comment: 20 pages, 11 figures, 6 tables, accepted for publication in A&A

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

The processes that lead from molecular clouds to protostars, stellar associations, and the evolution of both, are active research topics. The target of this study, NGC2264, is a benchmark star-forming region (SFR) in which these issues can be profitably studied. We assembled a new extensive sample of NGC 2264 members. To this end we used new X-ray data obtained with the XMM-Newton telescope, GAIA eDR3 data, and an extensive collection of public and published catalogs. Following a previous suggestion that the SFR might extend significantly beyond the better studied areas, our search covers a wide 2.5x2.5 degrees region in the sky. Our catalog comprises more than 2200 candidate members, a ~100% increase over previous determinations. We analyze their spatial distribution and define new substructures. Using GAIA parallaxes we estimate a new average distance to NGC2264 of 722+/-2pc and suggest that the embedded Spokes subregion is ~20pc farther away within the molecular cloud. A complex dynamics is unveiled by the available proper motions and radial velocities: we observe signs of global expansion and rotation. At the same time, we observe the collapse and coalescence of two substructures in a region where active star formation is taking place. The fraction of stars with disks and of those undergoing circumstellar accretion varies significantly across the field, suggesting that star formation has been occurring for several million years. A particularly low disk fraction around the O VII star S Mon might be attributed to external disk photoevaporation or to an older age of the stars in the region. NGC2264 is not dynamically relaxed and its present configuration is the result of multiple dynamical processes. The cloud has been forming stars for several million years and we identify the process that is likely responsible for the ongoing formation activity., Comment: 24 pages, 18 figures, accepted for publication on Astronomy & Astrophysics

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

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

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

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters., Comment: 26 pages, 14 figures, accepted for publication in Astronomy & Astrophysics; language-edited version; two appendices merged

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

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

We analyze high-resolution (dv=<10km/s) optical and infrared spectra covering the [OI] 6300 angstrom and [NeII] 12.81 micron lines from a sample of 31 disks in different evolutionary stages. Following work at optical wavelengths, we use Gaussian profiles to fit the [NeII] lines and classify them into HVC (LVC) if the line centroid is more (less) blueshifted than 30 km/s with respect to the stellar radial velocity. Unlike for the [OI] where a HVC is often accompanied by a LVC, all 17 sources with a [NeII] detection have either a HVC or a LVC. [NeII] HVCs are preferentially detected toward high accretors (Macc > 10$^{-8}$ Msun/yr) while LVCs are found in sources with low Macc, low [OI] luminosity, and large infrared spectral index (n13-31). Interestingly, the [NeII] and [OI] LVC luminosities display an opposite behaviour with n13-31: as the inner dust disk depletes (higher n13-31) the [NeII] luminosity increases while the [OI] weakens. The [NeII] and [OI] HVC profiles are generally similar with centroids and FWHMs showing the expected behaviour from shocked gas in micro-jets. In contrast, the [NeII] LVC profiles are typically more blueshifted and narrower than the [OI] profiles. The FWHM and centroid vs. disk inclination suggest that the [NeII] LVC predominantly traces unbound gas from a slow, wide-angle wind that has not lost completely the Keplerian signature from its launching region. We sketch an evolutionary scenario that could explain the combined [OI] and [NeII] results and includes screening of hard (~1keV) X-rays in inner, mostly molecular, MHD winds., Comment: Accepted for publication in The Astrophysical Journal

Pre-main sequence stars are variable sources. In diskless stars this variability is mainly due to the rotational modulation of dark photospheric spots and active regions, as in main sequence stars even if associated with a stronger magnetic activity. Aims. We aim at analyzing the simultaneous optical and X-ray variability in these stars to unveil how the activity in the photosphere is connected with that in the corona, to identify the dominant surface magnetic activity, and to correlate our results with stellar properties, such as rotation and mass. Methods. We analyzed the simultaneous optical and X-ray variability in stars without inner disks (e.g., class III objects and stars with transition disks) in NGC 2264 from observations obtained with Chandra/ACIS-I and CoRoT as part of the Coordinated Synoptic Investigation of NGC 2264. We searched for those stars whose optical and X-ray variability is correlated, anti-correlated, or not correlated by sampling their optical and X-ray light curves in suitable time intervals and studying the correlation between the flux observed in optical and in X-rays. We then studied how this classification is related with stellar properties. Results. Starting from a sample of 74 class III/transition disk (TD) stars observed with CoRoT and detected with Chandra with more than 60 counts, we selected 16 stars whose optical and X-ray variability is anti-correlated, 11 correlated, and 17 where there is no correlation. The remaining stars did not fall in any of these groups. We interpreted the anti-correlated optical and X-ray variability as typical of spot-dominated sources, due to the rotational modulation of photospheric spots spatially coincident to coronal active regions, and correlated variability typical of faculae-dominated sources, where the brightening due to faculae is dominant over the darkening due to spots. [Conclusions not shown in the pre-print], Comment: Accepted for publication on Astronomy & Astrophysics. Appendices are not included in this preprint

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

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

[Abridged] We investigated the X-ray characteristics of the Class I YSO Elias 29 with joint XMM-Newton and NuSTAR observations of 300 ks and 450 ks, respectively. These are the first observations of a very young (<1 Myr) stellar object in a band encompassing simultaneously both soft and hard X-rays. In addition to the hot Fe complex at 6.7 keV, we observed fluorescent emission from Fe at $\sim6.4$ keV, confirming the previous findings. The line at 6.4 keV is detected during quiescent and flaring states and its flux is variable. The equivalent width is found varying in the $\approx 0.15--0.5$ keV range. These values make unrealistic a simple model with a centrally illuminated disk and suggest a role of the cavity containing Elias 29 and possibly reverberation processes that could occur in it. We observed two flares, with duration of 20 ks and 50 ks, respectively. We systematically observed an increase of $N_H$ during the flares of a factor five. This behavior has been observed during flares previously detected in Elias 29 with XMM-Newton and ASCA. The phenomenon hints that the flaring regions could be buried under the accretion streams and at high stellar latitudes, as the X-rays from flares pass through gas denser than the gas along the line of sight of the quiescent corona. In a different scenario, a contribution from scattered soft photons to the primary coronal emission could mimic a shallower $N_H$ in the quiescent spectrum. In the spectrum of the full NuSTAR exposure, we detect hard X-ray emission in the band $\approx20-80$ keV in excess with respect to the thermal emission. The hard X-ray emission could be due to a population of energetic electrons accelerated by the magnetic field along the accretion streams. These particles could concur to pumping up the Fe fluorescence of cold Fe of the disk along with X-ray photons with $E>7.11$ keV., Comment: 15 pages, 14 figures, accepted for publication on A&A

Flares are powerful events ignited by a sudden release of magnetic energy. With the aim of studying flares in the 125-Myr-old stars in the Pleiades observed simultaneously in optical and X-ray light, we obtained new XMM-Newton observations of this cluster during the observations of Kepler K2 Campaign 4. Our objective is to characterize the most powerful flares observed in both bands and to constrain the energy released in the optical and X-ray, the geometry of the loops, and their time evolution. We aim to compare our results to existing studies of flares occurring in the Sun and stars at different ages. We selected bright X-ray/optical flares occurred in 12 known members of the Pleiades from their K2 and XMM-Newton light curves. The sample includes ten K-M stars, one F9 star, and one G8 star. Flare average properties were obtained from integrated analysis of the light curves during the flares. The time evolution of the plasma in the magnetic loops is constrained with time-resolved X-ray spectral analysis. Most of the flares studied in this work emitted more energy in optical than in X-rays, as in most solar flares, even if the Pleiades flares output a larger fraction of their total energy in X-rays than typical solar flares do. Additionally, the energy budget in the two bands is weakly correlated. We also found comparable flare duration in optical and X-rays and observed that rapidly rotating stars (e.g., with rotation period shorter than 0.5 days) preferentially host short flares. We estimated the slope of the cooling path of the flares in the log(EM)-versus-log(T) plane. The values we obtained are affected by large uncertainties, but their nominal values suggest that the flares analyzed in this paper are mainly due to single loops with no sustained heating occurring during the cooling phase. We also observed and analyzed oscillations with a period of 500 s during one of the flares., Comment: Accepted for publication by Astronomy & Astrophysics. Some figures are included with low resolution. The abstract is shorter than that of the orginal paper

We analyze the X-ray spectra of the $\sim$8000 sources detected in the Cygnus OB2 Chandra Legacy Survey (Drake et al., this issue), with the goals of characterizing the coronal plasma of the young low-mass stars in the region and estimating their intrinsic X-ray luminosities. We adopt two different strategies for X-ray sources for which more or less than 20 photons were detected. For the brighter sample we fit the spectra with absorbed isothermal models. In order to limit uncertainties, for most of the fainter Cygnus OB2 members in this sample, we constrain the spectral parameters to characteristic ranges defined from the brightest stars. For X-ray sources with $<$20 net photons we adopt a conversion factor from detected photon flux to intrinsic flux. This was defined, building on the results for the previous sample, as a function of the 20% quantile of the detected photon energy distributions, which we prove to also correlate well with extinction. We then use the X-ray extinction from the spectral fits to constrain the ratio between optical and X-ray extinction toward Cygnus OB2, finding it consistent with standard "Galactic"' values, when properly accounting for systematics. Finally we exploit the large number of sources to constrain the average coronal abundances of several elements, through two different ensemble analyses of the X-ray spectra of low-mass Cygnus OB2 members. We find the pattern of abundances to be largely consistent with that derived for the young stellar coronae in the Orion Nebula Cluster., Comment: 28 pages, 16 figures, accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Survey

Flares from the Sun and other stars are most prominently observed in the soft X-ray band. Most of the radiated energy, however, is released at optical/UV wavelengths. In spite of decades of investigation, the physics of flares is not fully understood. Even less is known about the powerful flares routinely observed from pre-main sequence stars, which might significantly influence the evolution of circumstellar disks. Observations of the NGC2264 star forming region were obtained in Dec. 2011, simultaneously with three telescopes, Chandra (X-rays), CoRoT (optical), and Spitzer (mIR), as part of the "Coordinated Synoptic Investigation of NGC2264" (CSI-NGC2264). Shorter Chandra and CoRoT observations were also obtained in March 2008. We analyzed the lightcurves to detect X-ray flares with an optical and/or mIR counterpart. Basic flare properties from the three datasets, such as emitted energies and peak luminosities, were then compared to constrain the spectral energy distribution of the flaring emission and the physical conditions of the emitting regions. Flares from stars with and without circumstellar disks were also compared to establish any difference that might be attributed to the presence of disks. Seventy-eight X-ray flares with an optical and/or mIR counterpart were detected. Their optical emission is found to correlate well with, and to be significantly larger than, the X-ray emission. The slopes of the correlations suggest that the difference becomes smaller for the most powerful flares. The mIR flare emission seems to be strongly affected by the presence of a circumstellar disk: flares from stars with disks have a stronger mIR emission with respect to stars without disks. This might be attributed to the reprocessing of the optical (and X-ray) flare emission by the inner circumstellar disk, providing evidence for flare-induced disk heating., Comment: 16 pages (36 including appendixes), 8 figures (main text), accepted for publication by Astronomy & Astrophysics (section 8)

The GES survey using FLAMES at the VLT has obtained high-resolution UVES spectra for a large number of giant stars, allowing a determination of the abundances of the key chemical elements C and N at their surface. The surface abundances of these chemical species are well-known to change in stars during their evolution on the red giant branch after the first dredge-up episod, as a result of extra-mixing phenomena. We investigate the effects of thermohaline mixing on C and N abundances using the first comparison between the GES [C/N] determinations with simulations of the observed fields using a model of stellar population synthesis. We explore the effects of thermohaline mixing on the chemical properties of giants through stellar evolutionary models computed with the stellar evolution code STAREVOL. We include these stellar evolution models in the Besan\c{c}on Galaxy model to simulate the [C/N] distributions determined from the UVES spectra of the GES and compare them with the observations. Theoretical predictions including the effect of thermohaline mixing are in good agreement with the observations. However, the field stars in the GES with C and N-abundance measurements have a metallicity close to solar, where the efficiency of thermohaline mixing is not very large. The C and N abundances derived by the GES in open and globular clusters clearly show the impact of thermohaline mixing at low-metallicity, allowing to explain the [C/N] ratio observed in lower-mass and older giant stars. Using independent observations of carbon isotopic ratio in clump field stars and open clusters, we also confirm that thermohaline mixing should be taken into account to explain the behavior of 12C/13C ratio as a function of stellar age. Overall the current model including thermohaline mixing is able to reproduce very well the C- and N-abundances over the whole metallicity range investigated by the GES data., Comment: Accepted for publication in A&A ; 12 pages, 12 figures

We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission. Diffuse emission was detected in the sub-bands Soft [0.5 : 1.2] and Medium [1.2 : 2.5], and marginally in the Hard [2.5 : 7.0] keV band. From X-ray spectral analysis of stacked spectra we compute a total [0.5 : 7.0 keV] diffuse X-ray luminosity of L$_{\rm x}^{\rm diff}\approx$4.2$\times$10$^{\rm 34}$ erg s$^{-1}$, characterized with plasma temperature components at kT$\approx$ 0.11, 0.40 and 1.18 keV, respectively. The HI absorption column density corresponding to these temperatures has a distribution consistent with N$_{\rm H}$ = 0.43, 0.80 and 1.39 $\times$10$^{22}$ cm$^{-2}$. The extended medium band energy emission likely arises from O-type stellar winds thermalized by wind-wind collisions in the most populated regions of the association, while the soft band emission probably arises from less energetic termination shocks against the surrounding Interstellar-Medium. Super-soft and Soft diffuse emission appears more widely dispersed and intense than the medium band emission. The diffuse X-ray emission is generally spatially coincident with low-extinction regions that we attribute to the ubiquitous influence of powerful stellar winds from massive stars and their interaction with the local Interstellar-Medium. Diffuse X-ray emission is volume-filling, rather than edge-brightened, oppositely to other star-forming regions. We reveal the first observational evidence of X-ray haloes around some evolved massive stars., Comment: Accepted for publication on ApJSS, Cygnus OB2 Special Issue (2018). 20 pages, 19 figures

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

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

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

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

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

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

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

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

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

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

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

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

Pre-main sequence stars are variable sources. In stars with disks, this variability is related to the morphology of the inner circumstellar region (<0.1 AU) and that of the photosphere and corona, all impossible to be spatially resolved with present day techniques. This has been the main motivation for the Coordinated Synoptic Investigation of NGC 2264. In this paper, we focus on the stars with disks. We analyze the X-ray spectral properties extracted during optical bursts and dips in order to unveil the nature of these phenomena. We analyze simultaneous CoRoT and Chandra/ACIS-I observations to search for coherent optical and X-ray flux variability in stars with disks. Then, stars are analyzed in two different samples. In stars with variable extinction, we look for a simultaneous increase of optical extinction and X-ray absorption during the optical dips; in stars with accretion bursts, we search for soft X-ray emission and increasing X-ray absorption during the bursts. Results. We find evidence for coherent optical and X-ray flux variability among the stars with variable extinction. In 9/24 stars with optical dips, we observe a simultaneous increase of X-ray absorption and optical extinction. In seven dips, it is possible to calculate the NH/AV ratio in order to infer the composition of the obscuring material. In 5/20 stars with optical accretion bursts, we observe increasing soft X-ray emission during the bursts that we associate to the emission of accreting gas. It is not surprising that these properties are not observed in all the stars with dips and bursts, since favorable geometric configurations are required. The observed variable absorption during the dips is mainly due to dust-free material in accretion streams. In stars with accretion bursts, we observe on average a larger soft X-ray spectral component not observed in non accreting stars., Comment: Accepted for publication by Astronomy & Astrophysics

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

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

The low spin rates measured for solar-type stars at an age of a few Myr (~10% of the break-up velocity) indicate that some mechanism of angular momentum regulation must be at play in the early pre-main sequence. We characterize the rotation properties for members of the region NGC 2264 (~3 Myr), and investigate the accretion-rotation connection at an age where about 50% of the stars have already lost their disks. We examined a sample of 500 cluster members whose photometric variations were monitored in the optical for 38 consecutive days with CoRoT. Light curves were analyzed for periodicity using the Lomb-Scargle periodogram, the autocorrelation function and the string-length method. The period distribution obtained for the cluster consists of a smooth distribution centered around P=5.2 d with two peaks at P=1-2 d and 3-4 d. A separate analysis of CTTS and WTTS indicates that the P=1-2 d peak is associated with the latter, while both groups contribute to the P=3-4 d peak. The comparison between CTTS and WTTS supports the idea of a rotation-accretion connection: their respective rotational properties are statistically different, and CTTS rotate on average more slowly than WTTS. We also observe a clear dearth of fast rotators with strong accretion signatures (large UV flux excess). This is consistent with earlier findings that fast rotators in young star clusters are typically devoid of dusty disks. Our sample shows some evidence of a mass dependence in the rotation properties of NGC 2264 members, lower-mass stars spinning on average faster. This study confirms that disks influence the rotational evolution of young stars. The idea of disk-locking may be consistent with the picture of rotation and rotation-accretion connection that we observe for the NGC 2264 cluster. However, the origin of the several substructures that we observe in the period distribution deserves further investigation., Comment: 45 pages (27 pages for the main text, plus data table and 5 Appendices); 4 tables, 24 figures. Accepted for publication in Astronomy & Astrophysics. Abstract shortened for the arXiv listing

The space telescope CoRoT has provided light curves of T Tauri stars belonging to the star-forming region of NGC 2264 with unprecedented continuity and precision in the framework of a coordinated multi-wavelength observational project. We perform spot modelling of the optical light curves of five weak-line T Tauri stars whose variability is dominated by starspots. We apply two-spot and evolving single spot models in the framework of a Bayesian Monte Carlo Markov Chain approach to derive the a posteriori distribution of the starspot parameters and the inclination of the star rotation axis. We focus on the rotation periods of the spots that can provide evidence for differential rotation in those stars. We find meaningful results in the case of three stars with an inclination higher than 50 deg and a slow variation of the light curve amplitude. The relative difference of the spot rotation periods ranges from 0.02 to 0.05 that is 3-10 times larger than the amplitude of the differential rotation found in similar stars with Doppler imaging techniques. We conclude that the intrinsic starspot evolution, although very slow, has a significant impact on the determination of the differential rotation by means of our spot modelling approach. We estimate typical timescales for the evolution of the starspot pattern between about 20 and 50 rotation periods in our stars., Comment: 12 pages, 13 figures, 11 tables; research note accepted by Astronomy and Astrophysics

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., Comment: Accepted to A&A; 9 pages, 5 figures; Full Table 2 in source

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy., Comment: Accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Project

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

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

We investigate the uncertainties of fitted X-ray model parameters and fluxes for relatively faint Chandra ACIS-I source spectra. Monte-Carlo (MC) simulations are employed to construct a large set of 150,000 fake X-ray spectra in the low photon count statistics regime (from 20 to 350 net counts) using the XSPEC spectral model fitting package. The simulations employed both absorbed thermal (APEC) and non-thermal (power-law) models, in concert with the Chandra ACIS-I instrument response and interstellar absorption. Simulated X-ray spectra were fit assuming a wide set of different input parameters and C-statistic minimization criteria to avoid numerical artifacts in the accepted solutions. Results provide an error estimate for each parameter (absorption, NH, plasma temperature, kT, or power-law slope, Gamma, and flux, and for different background contamination levels. The distributions of these errors are studied as a function of the 1 sigma quantiles and we show how these correlate with different model parameters, net counts in the spectra and relative background level. Maps of uncertainty in terms of the 1 sigma quantiles for parameters and flux are computed as a function of spectrum net counts. We find very good agreement between our estimated X-ray spectral parameter and flux uncertainties and those recovered from spectral fitting for a subset of the X-ray sources detected in the Chandra Cygnus OB2 Legacy Survey diagnosed to be Association members and that have between 20 and 350 net counts. Our method can provide uncertainties for spectral parameters whenever formal X-ray spectral fits cannot be well-constrained, or are unavailable, and predictions useful for computing Chandra ACIS-I exposure times for observation planning., Comment: Accepted for publication to the Astrophysical Journal Supplement Series on the Cygnus OB2 Special Issue (2016). 18 pages, 13 figures

The extension of the corona of classical T Tauri stars (CTTS) is under discussion. The standard model of magnetic configuration of CTTS predicts that coronal magnetic flux tubes connect the stellar atmosphere to the inner region of the disk. However, differential rotation may disrupt these long loops. The results from Hydrodynamic modeling of X-ray flares observed in CTTS confirming the star-disk connection hypothesis are still controversial. Some authors suggest the presence of the accretion disk prevent the stellar corona to extent beyond the co-rotation radius, while others simply are not confident with the methods used to derive loop lengths. We use independent procedures to determine the length of flaring loops in stars of the Orion Nebula Cluster previously analyzed using Hydrodynamic models. Our aim is to disentangle between the two scenarios proposed. We present a different approach to determine the length of flaring loops based on the oscillatory nature of the loops after strong flares. We use wavelet tools to reveal oscillations during several flares. The subsequent analysis of such oscillations is settle on the Physics of coronal seismology. Our results likely confirm the large extension of the corona of CTTS and the hypothesis of star-disk magnetic interaction in at least three CTTS of the Orion Nebula Cluster. Analyzing oscillations in flaring events is a powerful tool to determine the physical characteristics of magnetic loops in coronae in stars other than the Sun. The results presented in this work confirm the star-disk magnetic connection in CTTS., Comment: Accepted for publication in A&A

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

The magnetic activity of late-type MS stars is characterized by different regimes, and their activity levels are well described by Ro, the ratio between P_rot and the convective turnover time. Very young PMS stars show, similarly to MS stars, intense magnetic activity. However they do not show clear activity-rotation trends, and it still debated which stellar parameters determine their magnetic activity levels. To bridge the gap between MS and PMS stars, we studied the activity-rotation relation in the young cluster h Per, a ~13 Myr old cluster, that contains both fast and slow rotators, whose members have ended their accretion phase and have already developed a radiative core. It offers us the opportunity to study the activity level of intermediate-age PMS stars with different rotational velocities, excluding any interactions with the circumstellar environment. We constrained the magnetic activity levels of h Per members measuring their X-ray emission from a Chandra observation, while P_rot were obtained by Moraux et al. (2013). We collected a final catalog of 414 h Per members with known P_rot, T_eff, M_star, with 169 of them having also detected X-ray emission. We found that h Per members, with 1.0 M_sun < M_star < 1.4 M_sun, display different activity regimes: fast rotators show supersaturation, while slower rotators have activity levels compatible to the non-saturated regime. At 13 Myr h Per is therefore the youngest cluster showing activity-rotation regimes analogous to that of MS stars, indicating that, at this age, magnetic field production is likely regulated by the alpha-Omega type dynamo. Moreover we observed that supersaturation is better described by P_rot than Ro, and that the observed patterns are compatible with the hypothesis of centrifugal stripping. In this scenario we inferred that coronae can produce structures as large as ~2 R_star above the stellar surface., Comment: 36 pages (the paper is 14 pages; two tables to be published online are 22 pages), 13 figures, accepted for publication in A&A

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

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

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+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., Comment: 14 pages, 13 figures, A&A accepted

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., Comment: Accepted for publication

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., Comment: 21 pages, 18 figures, accepted for publication in Astronomy & Astrophysics