Search

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

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: Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims: This study aims to investigate correlations between the Mean Interdistance ($\bar{D_{i}}$), Mean Closest Interdistance ($\bar{D_{c}}$) and Median Weighted Central Interdistance ($\bar{D_{cc}}$) with the age of star clusters, examining their evolutionary trends and assessing the robustness of these quantities as possible age indicators. Methods: We selected a sample of open clusters in the solar region and with a representative number of members (e.g. well populated and without outliers). The interdistances are derived from the spatial distribution of member stars within a cluster. Their evolution over time allows us to use them as an age indicators for star clusters. Results: Our investigation reveals a high-significant correlation between the interdistances and cluster age. Considering the full sample of clusters between 7 and 9 kpc, the relationship is very broad. This is due to uncertainties in parallax, which increase with increasing distance. In particular, we must limit the sample to a maximum distance from the Sun of about 200 pc to avoid artificial effects on cluster shape and on the spatial distribution of their stars along the line of sight. Conclusions: By conservatively restraining the distance to a maximum of $\sim$200 pc, we have established a relationship between the interdistances and the age of the clusters. In our sample, the relationship is mainly driven by the internal expansion of the clusters, and is marginally affected by external perturbative effects. Such relation might enhance our comprehension of cluster dynamics and might be used to derive cluster dynamical ages., Comment: This paper was accepted for publication in Astronomy & Astrophysics on June 28, 2024. The document includes 9 main figures and 4 additional figures in the appendix. Additionally, it contains one main table and an extract of another table in the appendix, available at the CDS

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

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

Spectroscopy from the final internal data release of the Gaia-ESO Survey (GES) has been combined with Gaia EDR3 to assign membership probabilities to targets observed towards 63 Galactic open clusters and 7 globular clusters. The membership probabilities are based chiefly on maximum likelihood modelling of the 3D kinematics of the targets, separating them into cluster and field populations. From 43211 observed targets, 13985 are identified as highly probable cluster members ($P>0.9$), with an average membership probability of 0.993. The addition of GES radial velocities successfully drives down the fraction of false positives and we achieve better levels of discrimination in most clusters over the use of astrometric data alone, especially those at larger distances. Since the membership selection is almost purely kinematic, the union of this catalogue with GES and Gaia is ideal for investigating the photometric and chemical properties of clusters as a function of stellar mass, age and Galactic position., Comment: Accepted version for publication in MNRAS. 16 pages + 38 pages of Appendices

In large scale maps of the Galactic disc, the Carina Arm stands out as a clear spiral feature, hosting prominent star clusters and associations rich in massive stars. We study the proper motions of 4199 O and early B most likely in the far Carina Arm, at distances mainly in excess of 4 kpc from the Sun, within the sky region, $282^{\circ} < \ell < 294^{\circ}$ and $-3^{\circ} < b < +1^{\circ}$ (Galactic coordinates). The sample is constructed by extending an existing blue-selected catalogue, and cross-matching with Gaia EDR3 astrometry. The observed pattern of proper motions is modulated into a saw-tooth pattern, with full amplitude approaching 1 mas yr$^{-1}$, recurring roughly every 2--3 degrees of longitude (200--300 pc at the median OB-star distance of 5.8 kpc). Kinematic perturbation of underlying circular rotation is most likely present. The data also reveal a moving group containing $>50$ OB stars at $\ell \sim 286^{\circ}$, $b \sim -1^{\circ}.4$ behind the main run of the far arm. An analysis of relative proper motions is performed that yields an incidence of runaway O stars of at least 10\% (potentially $>20$\% when full space motions become available). To map where runaways have run away from, we set up simulations for the region that assume linear trajectories and test for trajectory impact parameter in order to identify likely ejection hot spots. We find the method currently gives good results for times of flight of up to $\sim$4 Myr. It shows convincingly that only NGC 3603 and Westerlund 2 have ejected OB stars in significant numbers. Indeed, both clusters have experienced intense spells of ejection between 0.6--0.9 and 0.5--0.8 Myr ago, respectively., Comment: 19 pages, 12 figures, (accepted 5th October for publication in MNRAS, in original form 6th July)

The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5

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

Reconstructing the Galactic evolution of lithium (Li) is the main tool used to constrain the source(s) of Li enrichment in the Galaxy. Recent results have suggested a decline in Li at supersolar metallicities, which may indicate reduced production. We exploit the unique characteristics of the Gaia-ESO Survey open star cluster sample to further investigate this issue and to better constrain the evolution of Li at high metallicity. We trace the the upper envelope of Li abundance versus metallicity evolution using 18 clusters and considering members that should not have suffered any Li depletion. At variance with previous claims, we do not find any evidence of a Li decrease at high metallicity. The most metal-rich clusters in the sample ([Fe/H] about 0.3) actually show the highest Li abundances, with A(Li) > 3.4. Our results clearly show that previous findings, which were based on field stars, were affected by selection effects. The metal-rich population in the solar neighbourhood is composed of relatively old and cool stars that have already undergone some Li depletion; hence, their measured Li does not represent the initial interstellar medium abundance, but a lower limit to it., Comment: 5 pages, 3 figures, accepted for publication in Astronomy & Astrophysics Letters

The Gaia-ESO Survey (GES) observed many open clusters as part of its programme to spectroscopically characterise the various Milky Way populations. GES spectroscopy and Gaia astrometry from its second data release are used here to assign membership probabilities to targets towards 32 open clusters with ages from 1-3800 Myr, based on maximum likelihood modelling of the 3D kinematics of the cluster and field populations. From a parent catalogue of 14398 individual targets, 5033 stars with uniformly determined 3D velocities, $T_{\rm eff}$, $\log g$ and chemistry are assigned cluster membership with probability $>0.9$, and with an average probability of 0.991. The robustness of the membership probabilities is demonstrated using independent membership criteria (lithium and parallax) in two of the youngest clusters. The addition of radial velocities improves membership discrimination over proper motion selection alone, especially in more distant clusters. The kinematically-selected nature of the membership lists, independent of photometry and chemistry, makes the catalogue a valuable resource for testing stellar evolutionary models and investigating the time evolution of various parameters., Comment: Accepted for publication in MNRAS

In the last three decades several hundred nearby members of young stellar moving groups (MGs) have been identified, but there has been less systematic effort to quantify or characterise young stars that do not belong to previously identified MGs. Using a kinematically unbiased sample of 225 lithium-rich stars within 100 pc, we find that only $50 \pm 10$ per cent of young ($\lesssim 125$ Myr), low-mass ($0.5

The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. Here, we present the IGAPS point source catalogue. It contains 295.4 million rows providing photometry in the filters, i, r, narrow-band Halpha, g and U_RGO. The IGAPS footprint fills the Galactic coordinate range, |b| < 5deg and 30deg < l < 215deg. A uniform calibration, referred to the Pan-STARRS system, is applied to g, r and i, while the Halpha calibration is linked to r and then is reconciled via field overlaps. The astrometry in all 5 bands has been recalculated on the Gaia DR2 frame. Down to i ~ 20 mag (Vega system), most stars are also detected in g, r and Halpha. As exposures in the r band were obtained within the IPHAS and UVEX surveys a few years apart, typically, the catalogue includes two distinct r measures, r_I and r_U. The r 10sigma limiting magnitude is ~21, with median seeing 1.1 arcsec. Between ~13th and ~19th magnitudes in all bands, the photometry is internally reproducible to within 0.02 magnitudes. Stars brighter than r=19.5 have been tested for narrow-band Halpha excess signalling line emission, and for variation exceeding |r_I-r_U| = 0.2 mag. We find and flag 8292 candidate emission line stars and over 53000 variables (both at >5sigma confidence). The 174-column catalogue will be available via CDS Strasbourg., Comment: 28 pages, 22 figures

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2., Comment: 41 pages, 12 figures, accepted in A&A

The very bright and compact massive young cluster, NGC 3603, has been cited as an example of a starburst in the Milky Way and compared with the much-studied R136/30 Doradus region in the Large Magellanic Cloud. Here we build on the discovery by Mohr-Smith et al. (2017) of a large number of reddened O stars around this cluster. We construct a list of 288 candidate O stars with proper motions, in a region of sky spanning 1.5x1.5 square degrees centered on NGC 3603, by cross-matching the Mohr-Smith et al. (2017) catalogue with Gaia DR2 (Gaia Collaboration et al. 2018). This provides the basis for a first comprehensive examination of the proper motions of these massive stars in the halo of NGC 3603, relative to the much better studied central region. We identify up to 11 likely O star ejections -- 8 of which would have been ejected between 0.60 and 0.95 Myr ago (supporting the age of c.1 Myr that has been attributed to the bright cluster centre). Seven candidate ejections are arranged in a partial ring to the south of the cluster core spanning radii of 9-18 arcmin (18-36 pc if the cluster is 7 kpc away). We also show that the cluster has a halo of a further 100 O stars extending to a radius of at least 5 arcmin, adding to the picture of NGC 3603 as a scaled down version of the R136/30 Dor region., Comment: 11 pages, with a 6-page table (Appendix B). Accepted by MNRAS, 21st March 2019

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

The Vela-Puppis region is known to host the Vela OB2 association as well as several young clusters featuring OB and pre-main sequence stars. Several spatial and kinematic subgroups have been identified in the recent years. By grouping stars based on their positions and velocity we can address the quetion of the dynamical history of the region and the mechanisms that drove stellar formation. The Gaia DR2 astrometry and photometry enables us to characterise the 3D spatial and 3D kinematic distribution of young stars and to estimate the ages of the identified components. We use an unsupervised classification method to group stars based on their proper motions and parallax. We perform a study of the expansion rates of the different identified groups from 3D velocities, and from corrected tangential velocities. We make use of theoretical isochrones to estimate ages. The young stars can be separated into seven main groups of different ages and kinematical distribution. All groups are found to be expanding, although the expansion is mostly not isotropic. The size of the region, the age substructure, and the anistropic expansion rates are compatible with a prolonged period of star formation in a turbulent molecular cloud, and that the dispersion of the stars cannot be explained by gas expulsion alone., Comment: accepted A&A

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

An unsettled question concerning the formation and distribution of massive stars is whether they must be born in massive clusters and, if found in less dense environments, whether they must have migrated there. With the advent of wide-area digital photometric surveys, it is now possible to identify massive stars away from prominent Galactic clusters without bias. In this study we consider 40 candidate OB stars found in the field around the young massive cluster, Westerlund 2, by Mohr-Smith et al (2017): these are located inside a box of 1.5x1.5 square degrees and are selected on the basis of their extinctions and K magnitudes. We present VLT/X-shooter spectra of two of the hottest O stars, respectively 11 and 22 arcmin from the centre of Westerlund 2. They are confirmed as O4V stars, with stellar masses likely to be in excess of 40 Msun. Their radial velocities relative to the non-binary reference object, MSP 182, in Westerlund 2 are -29.4 +/- 1.7 and -14.4 +/- 2.2 km/s, respectively. Using Gaia DR2 proper motions we find that between 8 and 11 early O/WR stars in the studied region (including the two VLT targets, plus WR 20c and WR 20aa) could have been ejected from Westerlund 2 in the last one million years. This represents an efficiency of massive-star ejection of up to 25%. On sky, the positions of these stars and their proper motions show a near N--S alignment. We discuss the possibility that these results are a consequence of prior sub-cluster merging combining with dynamical ejection., Comment: Accepted for publication in MNRAS, 13 July 2018. 16 pages, plus one-page table in an appendix

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

The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD46223, from NGC 2244 occurred 1.73 (+0.34,-0.25)Myr (1$\sigma$ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD46150, the age is set for the famous ionised region at more than ten times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40M$_{\odot}$ star in a thin sheet-like molecular cloud. We form the 135,000M$_{\odot}$ cloud from thermally-unstable diffuse interstellar medium under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5Myr, including cavity size, stellar age, magnetic field and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback., Comment: This is the Version of Record, published RAS Open Access CC BY online in MNRAS, 2018 Feb 13. 15 Pages, 9 figures

Previous studies of the rotation law in the outer Galactic disc have mainly used gas tracers or clump giants. Here, we explore A and F stars as alternatives: these provide a much denser sampling in the outer disc than gas tracers and have experienced significantly less velocity scattering than older clump giants. This first investigation confirms the suitability of A stars in this role. Our work is based on spectroscopy of $\sim$ 1300 photometrically-selected stars in the red calcium-triplet region, chosen to mitigate against the effects of interstellar extinction. The stars are located in two low Galactic latitude sightlines, at longitudes $\ell = 118^{\circ}$, sampling strong Galactic rotation shear, and $\ell = 178^{\circ}$, near the Anticentre. With the use of Markov Chain Monte Carlo parameter fitting, stellar parameters and radial velocities are measured, and distances computed. The obtained trend of radial velocity with distance is inconsistent with existing flat or slowly rising rotation laws from gas tracers (Brand & Blitz 1993; Reid et al. 2014). Instead, our results fit in with those obtained by Huang et al. (2016) from disc clump giants that favoured rising circular speeds. An alternative interpretation in terms of spiral arm perturbation is not straight forward. We assess the role that undetected binaries in the sample and distance error may have in introducing bias, and show that the former is a minor factor. The random errors in our trend of circular velocity are within $\pm 5$ km s$^{-1}$.

As part of the Accretion Discs in H$\alpha$ with OmegaCAM (ADHOC) survey, we imaged in r, i and H-alpha a region of 12x8 square degrees around the Orion Nebula Cluster. Thanks to the high-quality photometry obtained, we discovered three well-separated pre-main sequences in the color-magnitude diagram. The populations are all concentrated towards the cluster's center. Although several explanations can be invoked to explain these sequences we are left with two competitive, but intriguing, scenarios: a population of unresolved binaries with an exotic mass ratio distribution or three populations with different ages. Independent high-resolution spectroscopy supports the presence of discrete episodes of star formation, each separated by about a million years. The stars from the two putative youngest populations rotate faster than the older ones, in agreement with the evolution of stellar rotation observed in pre-main sequence stars younger than 4 Myr in several star forming regions. Whatever the final explanation, our results prompt for a revised look at the formation mode and early evolution of stars in clusters., Comment: A&A in print

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

Massive OB stars are critical to the ecology of galaxies, and yet our knowledge of OB stars in the Milky Way, fainter than $V \sim 12$, remains patchy. Data from the VST Photometric H$\alpha$ Survey (VPHAS+) permit the construction of the first deep catalogues of blue excess-selected OB stars, without neglecting the stellar field. A total of 14900 candidates with 2MASS cross-matches are blue-selected from a 42 square-degree region in the Galactic Plane, capturing the Carina Arm over the Galactic longitude range $282^{\circ} \lesssim \ell \lesssim 293^{\circ}$. Spectral energy distribution fitting is performed on these candidates' combined VPHAS+ $u,g,r,i$ and 2MASS $J,H,K$ magnitudes. This delivers: effective temperature constraints, statistically separating O from early-B stars; high-quality extinction parameters, $A_0$ and $R_V$ (random errors typically $< 0.1$). The high-confidence O-B2 candidates number 5915 and a further 5170 fit to later B spectral type. Spectroscopy of 276 of the former confirms 97% of them. The fraction of emission line stars among all candidate B stars is 7--8% . Greyer ($R_V > 3.5$) extinction laws are ubiquitous in the region, over the distance range 2.5--3 kpc to $\sim$10~kpc. Near prominent massive clusters, $R_V$ tends to rise, with particularly large and chaotic excursions to $R_V \sim 5$ seen in the Carina Nebula. The data reveal a hitherto unnoticed association of 108 O-B2 stars around the O5If$+$ star LSS 2063 ($\ell = 289.77^{\circ}$, $b = -1.22^{\circ}$). Treating the OB star scale-height as a constant within the thin disk, we find an orderly mean relation between extinction ($A_0$) and distance in the Galactic longitude range, $287.6^{\circ} < \ell < 293.5^{\circ}$, and infer the subtle onset of thin-disk warping. A halo around NGC 3603, roughly a degree in diameter, of $\sim$500 O-B2 stars with $4 < A_0 (\rm{mag}) < 7$ is noted., Comment: Accepted for publication in MNRAS. 26 pages, 20 figures, 4 tables, 1 appendix

We report the detection of an infrared selected transient which has lasted at least 5 years, first identified by a large mid-infrared and optical outburst from a faint X-ray source detected with the Chandra X-ray Observatory. In this paper we rule out several scenarios for the cause of this outburst, including a classical nova, a luminous red nova, AGN flaring, a stellar merger, and intermediate luminosity optical transients, and interpret this transient as the result of a Young Stellar Object (YSO) of at least solar mass accreting material from the remains of the dusty envelope from which it formed, in isolation from either a dense complex of cold gas or massive star formation. This object does not fit neatly into other existing categories of large outbursts of YSOs (FU Orionis types) which may be a result of the object's mass, age, and environment. It is also possible that this object is a new type of transient unrelated to YSOs., Comment: 13 pages, 6 figures; Accepted to Monthly Notices of the Royal Astronomical Society

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 present a semi-automatic procedure to obtain fundamental physical parameters and distances of classical Be (CBe) stars, based on the Barbier-Chalonge-Divan (BCD) spectrophotometric system. Our aim is to apply this procedure to a large sample of CBe stars detected by the IPHAS photometric survey, to determine their fundamental physical parameters and to explore their suitability as galactic structure tracers. In this paper we describe the methodology used and the validation of the procedure by comparing our results with those obtained from different independent astrophysical techniques for subsamples of stars in common with other studies. We also present a test case study of the galactic structure in the direction of the Perseus Galactic Arm, in order to compare our results with others recently obtained with different techniques and the same sample of stars. We did not find any significant clustering of stars at the expected positions of the Perseus and Outer Galactic Arms, in agreement with previous studies in the same area that we used for verification., Comment: 27 pages,15 figures, Accepted in A&A

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., Comment: Received February 18, 2016; accepted April 5, 2016, by A&A

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

We estimate the accretion rates of 235 Classical T Tauri star (CTTS) candidates in the Lagoon Nebula using $ugri$H$\alpha$ photometry from the VPHAS+ survey. Our sample consists of stars displaying H$\alpha$-excess, the intensity of which is used to derive accretion rates. For a subset of 87 stars, the intensity of the $u$-band excess is also used to estimate accretion rates. We find the mean variation in accretion rates measured using H$\alpha$ and $u$-band intensities to be $\sim$ 0.17 dex, agreeing with previous estimates (0.04-0.4 dex) but for a much larger sample. The spatial distribution of CTTS align with the location of protostars and molecular gas suggesting that they retain an imprint of the natal gas fragmentation process. Strong accretors are concentrated spatially, while weak accretors are more distributed. Our results do not support the sequential star forming processes suggested in the literature., Comment: Accepted for publication in MNRAS. 23 pages and figs. Tables 2 and 4 in full can be obtained by contacting the author, and will be given in the online version of the journal

O and early B stars are at the apex of galactic ecology, but in the Milky Way, only a minority of them may yet have been identified. We present the results of a pilot study to select and parametrise OB star candidates in the Southern Galactic plane, down to a limiting magnitude of $g=20$. A 2 square-degree field capturing the Carina Arm around the young massive star cluster, Westerlund 2, is examined. The confirmed OB stars in this cluster are used to validate our identification method, based on selection from the $(u-g, g-r)$ diagram for the region. Our Markov Chain Monte Carlo fitting method combines VPHAS+ $u, g, r, i$ with published $J, H, K$ photometry in order to derive posterior probability distributions of the stellar parameters $\log(\rm T_{\rm eff})$ and distance modulus, together with the reddening parameters $A_0$ and $R_V$. The stellar parameters are sufficient to confirm OB status while the reddening parameters are determined to a precision of $\sigma(A_0)\sim0.09$ and $\sigma(R_V)\sim0.08$. There are 489 objects that fit well as new OB candidates, earlier than $\sim$B2. This total includes 74 probable massive O stars, 5 likely blue supergiants and 32 reddened subdwarfs. This increases the number of previously known and candidate OB stars in the region by nearly a factor of 10. Most of the new objects are likely to be at distances between 3 and 6 kpc. We have confirmed the results of previous studies that, at these longer distances, these sight lines require non-standard reddening laws with $3.5

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., Comment: Matched catalog can be provided by the authors; Accepted for publication by ApJ

We present a catalogue of 247 photometrically and spectroscopically confirmed fainter classical Be stars (13 < r < 16) in the direction of the Perseus Arm of the Milky Way (-1 < b < +4, 120 < l < 140). The catalogue consists of 181 IPHAS-selected new classical Be stars, in addition to 66 objects that were studied by Raddi et al. (2013) more closely, and 3 stars identified as classical Be stars in earlier work. This study more than doubles the number known in the region. Photometry spanning 0.6 to 5 micron, spectral types, and interstellar reddenings are given for each object. The spectral types were determined from low-resolution spectra (lambda / Delta-lambda ~ 800-2000), to a precision of 1-3 subtypes. The interstellar reddenings are derived from the (r - i) colour, using a method that corrects for circumstellar disc emission. The colour excesses obtained range from E(B-V) = 0.3 up to 1.6 - a distribution that modestly extends the range reported in the literature for Perseus-Arm open clusters. For around half the sample, the reddenings obtained are compatible with measures of the total sightline Galactic extinction. Many of these are likely to lie well beyond the Perseus Arm., Comment: 20 pages, 13 figures, 5 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

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., Comment: Accepted in ApJ

We present the first results of our search for new, extended Planetary Nebulae (PNe) based on careful, systematic, visual scrutiny of the imaging data from the INT Photometric H-alpha Survey of the Northern Galactic Plane (IPHAS). The newly uncovered PNe will help to improve the census of this important population of Galactic objects that serve as key windows into the late stage evolution of low to intermediate mass stars. They will also facilitate study of the faint end of the ensemble Galactic PN luminosity function. The sensitivity and coverage of IPHAS allows PNe to be found in regions of greater extinction in the Galactic Plane and/or those PNe in a more advanced evolutionary state and at larger distances compared to the general Galactic PN population. Using a set of newly revised optical diagnostic diagrams in combination with access to a powerful, new, multi-wavelength imaging database, we have identified 159 true, likely and possible PNe for this first catalogue release. The ability of IPHAS to unveil PNe at low Galactic latitudes and towards the Galactic Anticenter, compared to previous surveys, makes this survey an ideal tool to contribute to the improvement of our knowledge of the whole Galactic PN population, Comment: 58 pages including the main article, data and image catalogues, Submitted to MNRAS

The INT/WFC Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS) is a 1800 square degrees imaging survey covering Galactic latitudes |b| < 5 deg and longitudes l = 30 to 215 deg in the r, i and H-alpha filters using the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) in La Palma. We present the first quality-controlled and globally-calibrated source catalogue derived from the survey, providing single-epoch photometry for 219 million unique sources across 92% of the footprint. The observations were carried out between 2003 and 2012 at a median seeing of 1.1 arcsec (sampled at 0.33 arcsec/pixel) and to a mean 5\sigma-depth of 21.2 (r), 20.0 (i) and 20.3 (H-alpha) in the Vega magnitude system. We explain the data reduction and quality control procedures, describe and test the global re-calibration, and detail the construction of the new catalogue. We show that the new calibration is accurate to 0.03 mag (rms) and recommend a series of quality criteria to select the most reliable data from the catalogue. Finally, we demonstrate the ability of the catalogue's unique (r-Halpha, r-i) diagram to (1) characterise stellar populations and extinction regimes towards different Galactic sightlines and (2) select H-alpha emission-line objects. IPHAS is the first survey to offer comprehensive CCD photometry of point sources across the Galactic Plane at visible wavelengths, providing the much-needed counterpart to recent infrared surveys., Comment: Accepted for publication in MNRAS. The catalogue is available in full from the survey website at http://www.iphas.org and has been submitted to Vizier

We present a three dimensional map of extinction in the Northern Galactic Plane derived using photometry from the IPHAS survey. The map has fine angular ($\sim 10$ arcmin) and distance (100 pc) sampling allied to a significant depth ($\gtrsim 5$ kpc). We construct the map using a method based on a hierarchical Bayesian model as previously described by Sale (2012). In addition to mean extinction, we also measure differential extinction, which arises from the fractal nature of the ISM, and show that it will be the dominant source of uncertainty in estimates of extinction to some arbitrary position. The method applied also furnishes us with photometric estimates of the distance, extinction, effective temperature, surface gravity, and mass for $\sim 38$ million stars. Both the extinction map and the catalogue of stellar parameters are made publicly available via http://www.iphas.org/extinction ., Comment: 17 pages, 13 figures, accepted for publication in MNRAS. Data can be downloaded from http://www.iphas.org/extinction . References updated

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is surveying the southern Milky Way in u, g, r, i and Halpha at 1 arcsec angular resolution. Its footprint spans the Galactic latitude range -5 < b < +5 at all longitudes south of the celestial equator. Extensions around the Galactic Centre to Galactic latitudes +/-10 bring in much of the Galactic Bulge. This ESO public survey, begun on 28th December 2011, reaches down to 20th magnitude (10-sigma) and will provide single-epoch digital optical photometry for around 300 million stars. The observing strategy and data pipelining is described, and an appraisal of the segmented narrowband Halpha filter in use is presented. Using model atmospheres and library spectra, we compute main-sequence (u - g), (g - r), (r - i) and (r - Halpha) stellar colours in the Vega system. We report on a preliminary validation of the photometry using test data obtained from two pointings overlapping the Sloan Digital Sky Survey. An example of the (u - g, g - r) and (r - Halpha, r - i) diagrams for a full VPHAS+ survey field is given. Attention is drawn to the opportunities for studies of compact nebulae and nebular morphologies that arise from the image quality being achieved. The value of the u band as the means to identify planetary-nebula central stars is demonstrated by the discovery of the central star of NGC 2899 in survey data. Thanks to its excellent imaging performance, the VST/OmegaCam combination used by this survey is a perfect vehicle for automated searches for reddened early-type stars, and will allow the discovery and analysis of compact binaries, white dwarfs and transient sources., Comment: 24 pages, 23 figures (supplementary tables not included, see www.vphas.eu for publisher's version). Accepted for publication by MNRAS, 24th February 2014. Replacing original arxiv upload to correct corrupt figures

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., Comment: Accepted for an ApJS Special Issue devoted to the Chandra Cygnus OB2 Legacy Survey

We present initial results of the first panoramic search for high-amplitude near-infrared variability in the Galactic Plane. We analyse the widely separated two-epoch K-band photometry in the 5th and 7th data releases of the UKIDSS Galactic Plane Survey. We find 45 stars with \Delta K > 1 mag, including 2 previously known OH/IR stars and a Nova. Even though the mid-plane is not yet included in the dataset, we find the majority (66%) of our sample to be within known star forming regions (SFRs), with two large concentrations in the Serpens OB2 association (11 stars) and the Cygnus-X complex (12 stars). Sources in SFRs show spectral energy distributions (SEDs) that support classification as Young Stellar Objects (YSOs). This indicates that YSOs dominate the Galactic population of high amplitude infrared variable stars at low luminosities and therefore likely dominate the total high amplitude population. Spectroscopic follow up of the DR5 sample shows at least four stars with clear characteristics of eruptive pre-main-sequence variables, two of which are deeply embedded. Our results support the recent concept of eruptive variability comprising a continuum of outburst events with different timescales and luminosities, but triggered by a similar physical mechanism involving unsteady accretion. Also, we find what appears to be one of the most variable classical Be stars., Comment: 26 pages, 16 Figures, accepted for publication in MNRAS

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., Comment: Accepted by ApJ

As part of a systematic search programme of a 10-degree wide strip of the Northern Galactic plane we present preliminary evidence for the discovery of four (and possibly five) new supernova remnants (SNRs). The pilot search area covered the 19-20 hour right ascension zone sampling from +20 to +55 degrees in declination using binned mosaic images from the INT Photometric H-alpha Survey (IPHAS). The optical identification of the candidate SNRs was based mainly on their filamentary and arc-like emission morphologies, their apparently coherent, even if fractured structure and clear disconnection from any diffuse neighbouring HII region type nebulosity. Follow-up optical spectroscopy was undertaken, sampling carefully across prominent features of these faint sources. The resulting spectra revealed typical emission line ratios for shock excited nebulae which are characteristic of SNRs, which, along with the latest diagnostic diagrams, strongly support the likely SNR nature of these sources: G038.7-1.3 (IPHASX J190640.5+042819); G067.6+0.9 (IPHASX J195744.9+305306); G066.0-0.0 (IPHASX J195749.2+290259) and G065.8-0.5 (IPHASX J195920.4+283740). A fifth possible younger, higher density nebula SNR candidate, G067.8+0.5 (IPHASX J200002.4+305035) was discovered about 5 arcmins to the west of IPHASX J195744.9+305306, and warrants further study. A multi-wavelength cross-check from available archived data in the regions of these candidates was also performed with a focus on possible radio counterparts., Comment: 14 pages, 8 figures. Accepted for publication in MNRAS

We investigate a region of the Galactic plane, between 120 <= l <= 140 and -1 <= b <= +4, and uncover a population of moderately reddened (E(B-V) \sim 1) classical Be stars within and beyond the Perseus and Outer Arms. 370 candidate emission line stars (13 <= r <= 16) selected from the INT Photometric H-alpha Survey of the Northern Galactic plane (IPHAS) have been followed up spectroscopically. A subset of these, 67 stars with properties consistent with those of classical Be stars, have been observed at sufficient spectral resolution (Delta_lambda \sim 2 - 4 Angstrom) at blue wavelengths to narrow down their spectral types. We determine these to a precision estimated to be +/- 1 sub-type and then we measure reddenings via SED fitting with reference to appropriate model atmospheres. Corrections for contribution to colour excess from circumstellar discs are made using an established scaling to H-alpha emission equivalent width. Spectroscopic parallaxes are obtained after luminosity class has been constrained via estimates of distances to neighbouring A/F stars with similar reddenings. Overwhelmingly, the stars in the sample are confirmed as luminous classical Be stars at heliocentric distances ranging from 2 kpc up to \sim 12 kpc. However, the errors are presently too large to enable the cumulative distribution function with respect to distance to distinguish between models placing the stars exclusively in spiral arms, or in a smooth exponentially-declining distribution., Comment: 21 pages, 13 figures. Accepted for publication in MNRAS

We identify 69 X-ray sources discovered by the Galactic Bulge Survey (GBS) that are coincident with, or very close to bright stars in the Tycho-2 catalog. Additionally, two other GBS sources are resolved binary companions to Tycho-2 stars where both components are separately detected in X-rays. Most of these are likely to be real matches, but we identify nine objects with large and significant X-ray to optical offsets as either detections of resolved binary companions or chance alignments. We collate known spectral types for these objects, and also examine 2MASS colors, variability information from the All-Sky Automated Survey (ASAS), and X-ray hardness ratios for the brightest objects. Nearly a third of the stars are found to be optically variable, divided roughly evenly between irregular variations and periodic modulations. All fall among the softest objects identified by the GBS. The sample forms a very mixed selection, ranging in spectral class from O9 to M3. In some cases the X-ray emission appears consistent with normal coronal emission from late-type stars, or wind emission from early-types, but the sample also includes one known Algol, one W UMa system, two Be stars, and several X-ray bright objects likely to be coronally active stars or binaries. Surprisingly, a substantial fraction of the spectroscopically classified, non-coincidental sample (12 out of 38 objects) have late B or A type counterparts. Many of these exhibit redder near-IR colors than expected for their spectral type and/or variability, and it is likely that the X-rays originate from a late-type companion star in most or all of these objects., Comment: 24 pages. Accepted for publication in the Astrophysical Journal Supplement

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