We present the results of a J, H, and K photometric variability survey of the central 0.78 square degrees of the young OB association Cygnus OB2. We used data observed with the Wide-Field CAMera at the United Kingdom Infrared Telescope in 2007 (spanning 217 days) to investigate the light curves of 5083 low mass candidate members in the association and explore the occurrence and main characteristics of their near-infrared variability. We identified 2529 stars ($\sim$50$\%$ of the sample) with significant variability with time-scales ranging from days to months. We classified the variable stars into the following three groups according to their light curve morphology: periodic variability (1697 stars), occultation variability (124 stars), and other types of variability (726 stars). We verified that the disk-bearing stars in our sample are significantly more variable in the near-infrared than diskless stars, with a steep increase in the disk-fraction among stars with higher variability amplitude. We investigated the trajectories described by variable stars in the color-space and measured slopes for 335 stars describing linear trajectories. Based on the trajectories in the color-space, we inferred that the sample analyzed is composed of a mix of young stars presenting variability due to hot and cold spots, extinction by circumstellar material, and changes in the disk emission in the near-infrared. We contemplated using the use of near-infrared variability to identify disk-bearing stars and verified that 53.4$\%$ of the known disk-bearing stars in our sample could have been identified as such based solely on their variability. We present 18 newly identified disk-bearing stars and 14 eclipsing binary candidates among CygOB2 lower-mass members.

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]

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.

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.

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.

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.

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.

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.

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

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.

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.

In the last decades, the early pre main sequence stellar rotational evolution picture has been constrained by studies targeting different young regions at a variety of ages. Observational studies suggest a mass-rotation dependence, and for some mass ranges a connection between rotation and the presence of a circumstellar disk. Not still fully explored, though, is the role of environmental conditions on the rotational regulation. We investigate the rotational properties of candidate members of the young massive association Cygnus OB2. The Stetson variability index, Lomb-Scargle periodogram, Saunders statistics, string/rope length method, and visual verification of folded light curves were applied to select 1224 periodic variable stars. Completeness and contamination of the periodic sample was derived from Monte Carlo simulations, out of which 894 periods were considered reliable. Our study was considered reasonably complete for periods from 2 to 30 days. The general rotational scenario seen in other young regions is confirmed by Cygnus OB2 period distributions, with disked stars rotating on average slower than non-disked stars. A mass-rotation dependence was also verified, but as in NGC 6530, lower mass stars are rotating on average slower than higher mass stars, with an excess of slow rotators among the lower mass population. The effect of the environment on the rotational properties of the association was investigated by re-analysing the results while taking into account the incident UV radiation arising from O stars in the association. Results compatible with the disk-locking scenario were verified for stars with low UV incidence, but no statistical significant relation between rotation and disk presence was verified for stars with high UV incidence suggesting that massive stars can have an important role on regulating the rotation of nearby low mass stars.

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.

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.

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects.

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.

The Chandra Cygnus OB2 Legacy Survey is a wide and deep X-ray survey of the nearby and massive Cygnus OB2 association. The survey has detected ~8,000 X-ray sources, the majority of which are pre-main sequence X-ray emitting young stars in the association itself. To facilitate quantitative scientific studies of these sources as well as the underlying OB association it is important to understand the sensitivity of the observations and the level of completeness the observations have obtained. Here we describe the use of a hierarchical Monte Carlo simulation to achieve this goal by combining the empirical properties of the observations, analytic estimates of the source verification process, and an extensive set of source detection simulations. We find that our survey reaches a 90% completeness level for a pre-main-sequence population at the distance of Cyg OB2 at an X-ray luminosity of 4 x 10^30 ergs/s and a stellar mass of 1.3 Msun for a randomly distributed population. For a spatially clustered population such as Cyg~OB2 the 90% completeness level is reached at 1.1 Msun instead, as the sources are more concentrated in areas of our survey with a high exposure. These simulations can easily be adapted for use with other X-ray observations and surveys, and we provide X-ray detection efficiency curves for a very wide array of source and background properties to allow these simulations to be easily exploited by other users.

The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and young low-mass stars and the feedback provided by massive stars on star and planet formation process. Cygnus OB2 has been recently observed with Chandra/ACIS-I as part of the 1.08Msec Chandra Cygnus OB2 Legacy Project. This survey detected 7924 X-ray sources in a square degree area centered on Cyg OB2. Since a proper classification and study of the observed X-ray sources also requires the analysis of their optical and infrared counterparts, we combined a large and deep set of optical and infrared catalogs available for this region with our new X-ray catalog. In this paper we describe the matching procedure and present the combined catalog containing 5703 sources. We also briefly discuss the nature of the X-ray sources with optical and infrared counterparts using their position in the color-magnitude and color-color diagrams.

Cygnus OB2 is the most massive association within 2kpc from the Sun, hosting hundreds of massive stars, thousands of young low mass members, and some sights of active star formation in the surrounding cloud. Recently, 10 photoevaporating proplyd-like objects with tadpole-shaped morphology were discovered in the outskirts of the OB association, approximately 6-14pc away from its center. The classification of these objects is ambiguous, being either evaporating residuals of the parental cloud which are hosting a protostar inside, or disk-bearing stars with an evaporating disk, such as the evaporating proplyds observed in the Trapezium Cluster in Orion. In this paper we present a study based on low resolution optical spectroscopic observations made with the Optical System for Imaging and low Resolution Integrated Spectroscopy (OSIRIS), mounted on the 10.4m Gran Telescopio CANARIAS (GTC), of two of these protostars. The spectrum of one of the objects shows evidence of accretion but not of outflows. In the latter object, the spectra show several emission lines indicating the presence of an actively accreting disk with outflow. We present estimates of the mass loss rate and the accretion rate from the disk, showing that the former exceeds the latter as observed in other known objects with evaporating disks. We also show evidence of a strong variability in the integrated flux observed in these objects, as well as in the accretion and outflow diagnostics.

The Cygnus OB2 association is the largest concentration of young and massive stars within 2 kpc of the Sun, including an estimated 65 O-type stars and hundreds of OB stars. The Chandra Cygnus OB2 Legacy Survey is a large imaging program undertaken with the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. The survey has imaged the central 0.5 deg^2 of the Cyg OB2 association with an effective exposure of 120ks and an outer 0.35 deg^2 area with an exposure of 60ks. Here we describe the survey design and observations, the data reduction and source detection, and present a catalog of 8,000 X-ray point sources. The survey design employs a grid of 36 heavily (~50%) overlapping pointings, a method that overcomes Chandra's low off-axis sensitivity and produces a highly uniform exposure over the inner 0.5 deg^2. The full X-ray catalog is described here and is made available online.

We report on the analysis of the Chandra-ACIS data of O, B and WR stars in the young association Cyg OB2. X-ray spectra of 49 O-stars, 54 B-stars and 3 WR-stars are analyzed and for the brighter sources, the epoch dependence of the X-ray fluxes is investigated. The O-stars in Cyg\,OB2 follow a well-defined scaling relation between their X-ray and bolometric luminosities: log(Lx/Lbol) = -7.2 +/- 0.2. This relation is in excellent agreement with the one previously derived for the Carina OB1 association. Except for the brightest O-star binaries, there is no general X-ray overluminosity due to colliding winds in O-star binaries. Roughly half of the known B-stars in the surveyed field are detected, but they fail to display a clear relationship between Lx and Lbol. Out of the three WR stars in Cyg OB2, probably only WR144 is itself responsible for the observed level of X-ray emission, at a very low log(Lx/Lbol) = -8.8 +/- 0.2. The X-ray emission of the other two WR-stars (WR145 and 146) is most probably due to their O-type companion along with a moderate contribution from a wind-wind interaction zone.

Attention is given to a population of 110 stars in the NGC 6611 cluster of the Eagle Nebula that have prominent near-infrared (NIR) excess and optical colours typical of pre-main sequence (PMS) stars older than 8 Myr. At least half of those for which spectroscopy exists have a Halpha emission line profile revealing active accretion. In principle, the V-I colours of all these stars would be consistent with those of young PMS objects (< 1 Myr) whose radiation is heavily obscured by a circumstellar disc seen at high inclination and in small part scattered towards the observer by the back side of the disc. However, using theoretical models it is shown here that objects of this type can only account for a few percent of this population. In fact, the spatial distribution of these objects, their X-ray luminosities, their optical brightness, their positions in the colour-magnitude diagram and the weak Li absorption lines of the stars studied spectroscopically suggest that most of them are at least 8 times older than the ~1 Myr-old PMS stars already known in this cluster and could be as old as ~30 Myr. This is the largest homogeneous sample to date of Galactic PMS stars considerably older than 8 Myr that are still actively accreting from a circumstellar disc and it allows us to set a lower limit of 7% to the disc frequency at ~16 Myr in NGC 6611. These values imply a characteristic exponential lifetime of ~6 Myr for disc dissipation.

The formation of stars in massive clusters is one of the main modes of the star formation process. However, the study of massive star forming regions is hampered by their typically large distances to the Sun. One exception to this is the massive star forming region Cygnus OB2 in the Cygnus X region, at the distance of about 1400 pc. Cygnus OB2 hosts very rich populations of massive and low-mass stars, being the best target in our Galaxy to study the formation of stars, circumstellar disks, and planets in presence of massive stars. In this paper we combine a wide and deep set of photometric data, from the r band to 24 micron, in order to select the disk bearing population of stars in Cygnus OB2 and identify the class I, class II, and stars with transition and pre-transition disks. We selected 1843 sources with infrared excesses in an area of 1 degree x 1 degree centered on Cyg OB2 in several evolutionary stages: 8.4% class I, 13.1% flat-spectrum sources, 72.9% class II, 2.3% pre-transition disks, and 3.3% transition disks. The spatial distribution of these sources shows a central cluster surrounded by a annular overdensity and some clumps of recent star formation in the outer region. Several candidate subclusters are identified, both along the overdensity and in the rest of the association.

The young open cluster NGC 6611 includes among its candidate members a class of peculiar objects with interesting properties: blue stars with infrared IR excesses. These stars show excesses in IR bands, signature of the presence of a circumstellar disk, but optical colors typical of older field stars. In order to confirm their membership to the cluster, it is therefore important to use new spectroscopic observations, together with previous photometric data. We aim at confirming the membership of these objects and at investigating their physical properties to verify whether the observed colors are intrinsic or altered by the presence of the disk or by the accretion processes. We analyze the intermediate resolution spectroscopic data obtained for a subsample of blue stars in NGC 6611 with FLAMES. In particular, we focus on the study of: 1) the profile of the Ha emission line, to select stars with accretion and outflow activity; 2) the Li absorption line, used as a youth indicator; 3) the radial velocity. Using the spectroscopic analysis, it has been possible to investigate the presence of the Li absorption line as well as to discriminate between stars with inert or active disk. In particular, from the analysis of the Ha emission line we were able to infer the activity due to the accretion and outflow processes and the variability of the emission. We also investigated the binarity of the blue stars and their membership to NGC 6611.

In order to fully understand the gravitational collapse of molecular clouds, the star formation process and the evolution of circumstellar disks, these phenomena must be studied in different Galactic environments with a range of stellar contents and positions in the Galaxy. The young massive association Cygnus OB2, in the Cygnus-X region, is an unique target to study how star formation and the evolution of circumstellar disks proceed in the presence of a large number of massive stars. We present a catalog obtained with recent optical observations in r,i,z filters with OSIRIS, mounted on the $10.4\,m$ GTC telescope, which is the deepest optical catalog of Cyg OB2 to date. The catalog consist of 64157 sources down to M=0.15 solar masses at the adopted distance and age of Cyg OB2. A total of 38300 sources have good photometry in all three bands. We combined the optical catalog with existing X-ray data of this region, in order to define the cluster locus in the optical diagrams. The cluster locus in the r-i vs. i-z diagram is compatible with an extinction of the optically selected cluster members in the 2.64<AV<5.57 range. We derive an extinction map of the region, finding a median value of AV=4.33 in the center of the association, decreasing toward the north-west. In the color-magnitude diagrams, the shape of the distribution of main sequence stars is compatible with the presence of an obscuring cloud in the foreground at about 850+/-25 pc from the Sun.

We present 3-dimensional kinematical observations of the massive OB association Cygnus OB2 to identify the mechanisms responsible for disrupting young star clusters. The picture revealed by these observations is of a highly-substructured, dynamically unmixed OB association that does not exhibit the position-velocity correlations predicted by the theories of infant mortality or tidal stripping. These observations would appear to support a picture of hierarchical star formation.

Mechanisms regulating the origin of X-rays in YSOs and the correlation with their evolutionary stage are under debate. Studies of the X-ray properties in young clusters allow to understand these mechanisms. One ideal target for this analysis is the Eagle Nebula (M16), with its central cluster NGC6611. At 1750 pc from the Sun, it harbors 93 OB stars, together with a population of low-mass stars from embedded protostars to disk-less Class III objects, with age <= 3Myrs. We study an archival 78 ksec Chandra/ACIS-I observation of NGC6611, and two new 80ksec observations of the outer region of M16, one centered on the Column V, and one on a region of the molecular cloud with ongoing star-formation. We detect 1755 point sources, with 1183 candidate cluster members (219 disk-bearing and 964 disk-less). We study the global X-ray properties of M16 and compare them with those of the Orion Nebula Cluster. We also compare the level of X-ray emission of Class II and Class III stars, and analyze the X-ray spectral properties of OB stars. Our study supports the lower level of X-ray activity for the disk-bearing stars with respect to the disk-less members. The X-ray Luminosity Function (XLF) of M16 is similar to that of Orion, supporting the universality of the XLF in young clusters. 85% of the O stars of NGC6611 have been detected in X-rays. With only one possible exception, they show soft spectra with no hard component, indicating that mechanisms for the production of hard X-ray emission in O stars are not operating in NGC 6611.

We report the discovery of ten proplyd-like objects in the vicinity of the massive OB association Cygnus OB2. They were discovered in IPHAS H-Alpha images and are clearly resolved in broad-band HST/ACS, near-IR and Spitzer mid-IR images. All exhibit the familiar tadpole shape seen in photoevaporating objects such as the Orion proplyds, with a bright ionization front at the head facing the central cluster of massive stars, and a tail stretching in the opposite direction. Many also show secondary ionization fronts, complex tail morphologies or multiple heads. We consider the evidence that these are either proplyds or `evaporating gaseous globules' (EGGs) left over from a fragmenting molecular cloud, but find that neither scenario fully explains the observations. Typical sizes are 50,000--100,000 AU, larger than the Orion proplyds, but in agreement with the theoretical scaling of proplyd size with distance from the ionizing source. These objects are located at projected separations of 6-14pc from the OB association, compared to 0.1pc for the Orion proplyds, but are clearly being photoionized by the 65 O-type stars in Cyg OB2. Central star candidates are identified in near- and mid-IR images, supporting the proplyd scenario, though their large sizes and notable asymmetries is more consistent with the EGG scenario. A third possibility is therefore considered, that these are a unique class of photoevaporating partially-embedded young stellar objects that have survived the destruction of their natal molecular cloud. This has implications for the properties of stars that form in the vicinity of massive stars.

We present a study of the radial distribution of dust species in young brown dwarf disks. Our work is based on a compositional analysis of the 10 and 20 micron silicate emission features for brown dwarfs in the Taurus-Auriga star-forming region. A fundamental finding of our work is that brown dwarfs exhibit stronger signs of dust processing in the cold component of the disk, compared to the higher mass T Tauri stars in Taurus. For nearly all of our targets, we find a flat disk structure, which is consistent with the stronger signs of dust processing observed in these disks. For the case of one brown dwarf, 2M04230607, we find the forsterite mass fraction to be a factor of ~3 higher in the outer disk compared to the inner disk region. Simple large-scale radial mixing cannot account for this gradient in the dust chemical composition, and some local crystalline formation mechanism may be effective in this disk. The relatively high abundance of crystalline silicates in the outer cold regions of brown dwarf disks provides an interesting analogy to comets. In this context, we have discussed the applicability of the various mechanisms that have been proposed for comets on the formation and the outward transport of high-temperature material. We also present Chandra X-ray observations for two Taurus brown dwarfs, 2M04414825 and CFHT-BD-Tau 9. We find 2M04414825, which has a ~12% crystalline mass fraction, to be more than an order of magnitude brighter in X-ray than CFHT-BD-Tau 9, which has a ~35% crystalline mass fraction. Combining with previous X-ray data, we find the inner disk crystalline mass fractions to be anti-correlated with the X-ray strength.

Different environmental conditions can play a crucial role in determining final products of the star formation process and in this context, less favorable activities of star formation are expected in the external regions of our Galaxy. We studied the properties of the young open cluster NGC 1893 located about 12 Kpc from the galactic center, to investigate how different physical conditions can affect the process of star formation. By adopting a multiwavelength approach, we compiled a catalog extending from X-rays to NIR data to derive the cluster membership. In addition, optical and NIR photometric properties are used to evaluate the cluster parameters. We find 415 diskless candidate members plus 1061 young stellar objects with a circumstellar disk or class II candidate members, 125 of which are also Halpha emitters. Considering the diskless candidate members, we find that the cluster distance is 3.6$\pm$0.2 kpc and the mean interstellar reddening is E(B-V)=0.6$\pm$0.1 with evidence of differential reddening in the whole surveyed region. NGC 1893 contains a conspicuous population of pre-main sequence stars together with the well studied main sequence cluster population; we found a disk fraction of about 70% similar to that found in clusters of similar age in the solar neighbour and then, despite expected unfavorable conditions for star formation, we conclude that very rich young clusters can form also in the outer regions of our Galaxy.

Massive SFR are characterized by intense ionizing fluxes, strong stellar winds and supernovae explosions, all of which have important effects on the surrounding media, on the star-formation (SF) process and on the evolution of YSOs and their disks. We present a multiband study of the massive young cluster NGC6611 and M16, to study how OB stars affect the early stellar evolution and the SF. We search for evidence of triggered SF by OB stars in NGC6611 on a large spatial scale (~10 pc) and how the efficiency of disks photoevaporation depends on the central stars mass. We assemble a multiband catalog with photometric data, from B band to 8.0micron, and X-ray data obtained with 2 new and 1 archival ACIS-I observation. We select the stars with disks from IR photometry and disk-less from X-ray emission, both in NGC6611 and the outer region of M16. We study induced photoevaporation searching for the spatial variation of disk frequency for distinct stellar mass ranges. The triggering of SF by OB stars has been investigated by deriving the history of SF across the nebula. We find evidence of sequential SF in the Eagle Nebula going from the SE (2.6 Myrs) to the NW (0.3 Myrs), with the median age of ~1 Myear. We observe a drop of the disk frequency close to OB stars (up to an average distance of 1 pc), without effects at larger distances. Furthermore, disks are more frequent around low-mass stars (<1 M(solar)) than in high-mass stars, regardless of the distance from OB stars. The SF chronology in M16 does not support the hypothesis of a large-scale SF triggered by OB stars in NGC6611. Instead, we speculate that it was triggered by the encounter (~3 Myrs ago) with a giant molecular shell created ~6 Myrs ago.

NGC6611 and its parental cloud, the Eagle Nebula (M16), are well-studied star-forming regions, thanks to their large content of both OB stars and stars with disks and the observed ongoing star formation. We identified 834 disk-bearing stars associated with the cloud, after detecting their excesses in NIR bands from J band to 8.0 micron. In this paper, we study in detail the nature of a subsample of disk-bearing stars that show peculiar characteristics. They appear older than the other members in the V vs. V-I diagram, and/or they have one or more IRAC colors at pure photospheric values, despite showing NIR excesses, when optical and infrared colors are compared. We confirm the membership of these stars to M16 by a spectroscopic analysis. The physical properties of these stars with disks are studied by comparing their spectral energy distributions (SEDs) with the SEDs predicted by models of T-Tauri stars with disks and envelopes. We show that the age of these stars estimated from the V vs. V-I diagram is unreliable since their V-I colors are altered by the light scattered by the disk into the line of sight. Only in a few cases their SEDs are compatible with models with excesses in V band caused by optical veiling. Candidate members with disks and photospheric IRAC colors are selected by the used NIR disk diagnostic, which is sensitive to moderate excesses, such as those produced by disks with low masses. In 1/3 of these cases, scattering of stellar flux by the disks can also be invoked. The photospheric light scattered by the disk grains into the line of sight can affect the derivation of physical parameters of ClassII stars from photometric optical and NIR data. Besides, the disks diagnostic we defined are useful for selecting stars with disks, even those with moderate excesses or whose optical colors are altered by veiling or photospheric scattered light.

Context: the observations of the proplyds in the Orion Nebula Cluster, showing clear evidence of ongoing photoevaporation, have provided a clear proof about the role of the externally induced photoevaporation in the evolution of circumstellar disks. NGC 6611 is an open cluster suitable to study disk photoevaporation, thanks to its large population of massive members and of stars with disk. In a previous work, we obtained evidence of the influence of the strong UV field generated by the massive cluster members on the evolution of disks around low-mass Pre-Main Sequence members. That work was based on a multi-band BVIJHK and X-ray catalog purposely compiled to select the cluster members with and without disk. Aims: in this paper we complete the list of candidate cluster members, using data at longer wavelengths obtained with Spitzer/IRAC, and we revisit the issue of the effects of UV radiation on the evolution of disks in NGC 6611. Methods: we select the candidate members with disks of NGC 6611, in a field of view of 33'x34' centered on the cluster, using IRAC color-color diagrams and suitable reddening-free color indices. Besides, using the X-ray data to select Class III cluster members, we estimate the disks frequency vs. the intensity of the incident radiation emitted by massive members. Results: we identify 458 candidate members with circumstellar disks, among which 146 had not been revealed in our previous work. Comparing of the various color indices we used to select the cluster members with disk, we claim that they detect the excesses due to the emission of the same physical region of the disk: the inner rim at the dust sublimation radius. Our new results confirm that UV radiation from massive stars affects the evolution of nearby circumstellar disks.

Context: the observation of young stars with circumstellar disks suggests that the disks are dissipated, starting from the inner region, by the radiation of the central star and eventually by the formation of rocky planetesimals, over a time scale of several million years. It was also shown that strong UV radiation emitted by nearby massive stars can heat a circumstellar disk up to some thousand degrees, inducing the photoevaporation of the gas. This process strongly reduces the dissipation time scale. Aims: we study whether there exists a correlation between the spatial distribution of stars with circumstellar disks and the position of massive stars with spectral class earlier than B5, in the open cluster NGC 6611. Methods: we created a multiband catalog of the cluster, down to V~23, using optical data from a WFI observation at 2.2m of ESO in the BVI bands, the 2MASS public point source catalog and an archival X-ray observation made with CHANDRA/ACIS. We selected the stars with infrared excess (due to the emission of a circumstellar disk) using suitable color indices independent of extinction, and studied their spatial distribution. Results: the spatial distribution of the stars with K band excess (due to the presence of a circumstellar disk) is anti correlated with that of the massive stars: the disks are more frequent at large distances from these stars. We argue that this is in agreement with the hypothesis that the circumstellar disks are heated by the UV radiation from the massive stars and photoevaporated.