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.

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.

The kinematic structure of the Cygnus OB2 association is investigated. No evidence of expansion or contraction is found at any scale within the region. Stars that are within $\sim$ 0.5 parsecs of one another are found to have more similar velocities than would be expected by random chance, and so it is concluded that velocity substructure exists on these scales. At larger scales velocity substructure is not found. We suggest that bound substructures exist on scales of $\sim$ 0.5 parsecs, despite the region as a whole being unbound. We further suggest that any velocity substructure that existed on scales > 0.5 parsecs has been erased. The results of this study are then compared to those of other kinematic studies of Cygnus OB2.

The kinematics of low-mass stars in nearby OB associations can provide clues about their origins and evolution. Combining the precise positions, proper motions and parallaxes given in the second Gaia Data Release with radial velocity measurements obtained with the Hermes spectrograph at the Anglo-Australian Telescope, we have an opportunity to study in detail the kinematics of low-mass stars belonging to the nearby $\gamma$ Vel cluster and the Vela OB2 association it is projected against. The presence of lithium is used to confirm the youth of our targets. We separate our sample into the cluster and association populations based on the membership probabilities of Jeffries et al. 2014, their parallaxes, and kinematics. We find strong evidence for expansion in the OB association population with at least 4$\sigma$ significance along all three axes, though the expansion is notably anisotropic. We discuss these results in the context of cluster and association dispersal theories.

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<M/M_{\odot}<1.0$) stars, are kinematically associated with known MGs. Whilst we find some evidence that six of the non-MG stars may be connected with the Lower Centaurus-Crux association, the rest form a kinematically "hotter" population, much more broadly dispersed in velocity, and with no obvious concentrations in space. The mass distributions of the MG members and non-MG stars is similar, but the non-MG stars may be older on average. We briefly discuss several explanations for the origin of the non-MG population.

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.

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.

We examine the mass-dependence of the velocity dispersion of stars in the young cluster NGC 6530 to better understand how it formed. Using a large sample of members we find that the proper motion velocity dispersion increases with stellar mass. While this trend is the opposite to that predicted if the cluster were developing energy equipartition, it is in agreement with recent N-body simulations that find such a trend develops because of the Spitzer instability. In these simulations the massive stars sink to the centre of the cluster and form a self-gravitating system with a higher velocity dispersion. If the cluster has formed by the cool collapse of an initially substructured distribution then this occurs within 1-2 Myr, in agreement with our observations of NGC 6530. We therefore conclude that NGC 6530 formed from much more extended initial conditions and has since collapsed to form the cluster we see now. This cluster formation model is inconsistent with the idea that all stars form in dense, compact clusters and provides the first dynamical evidence that star clusters can form by hierarchical mergers between subclusters.

We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We find that the fraction of unbound stars differs depending on the initial conditions. After 10 Myr, in initially highly subvirial, substructured simulations, the high-mass and lower-mass unbound fractions are similar at $\sim$23 per cent. In initially virialised, substructured simulations, we find only $\sim$16 per cent of all high-mass stars are unbound, whereas $\sim$37 per cent of all lower-mass stars are. The velocity distributions of unbound stars only show differences for extremely different initial conditions. The distributions are dominated by large numbers of lower-mass stars becoming unbound just above the escape velocity of $\sim$3 km s$^{-1}$ with unbound high-mass stars moving faster on average than lower-mass unbound stars. We see no high-mass runaway stars (velocity > 30 km s$^{-1}$) from any of our initial conditions and only an occasional lower-mass runaway star from initially subvirial/substructured simulations. In our simulations, we find a small number of lower-mass walkaway stars (with velocity 5-30 km s$^{-1}$) from all of our initial conditions. These walkaway stars should be observable around many nearby star-forming regions with Gaia.

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.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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

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.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

All cool stars show magnetic activity, and X-ray emission is the hallmark of this activity. Gaining an understanding of activity aids us in answering fundamental questions about stellar astrophysics and in determining the impact of activity on the exoplanets that orbit these stars. Stellar activity is driven by magnetic fields, which are ultimately powered by convection and stellar rotation. However, the resulting dynamo properties heavily depend on the stellar interior structure and are far from being understood. X-ray radiation can evaporate exoplanet atmospheres and damage organic materials on the planetary surface, reducing the probability that life can form or be sustained, but also provides an important source of energy for prebiotic chemical reactions. Over the next two years, the TESS mission will deliver a catalog of the closest exoplanets, along with rotation periods and activity diagnostics for millions of stars, whether or not they have a planet. We propose to include all cool stars that are TESS targets and bright enough for Chandra observations, as determined by their detection and flux in the ROSAT all-sky survey (RASS), to the list of Chandra Cool Attitude Targets (CATs). For each target, the signal will be sufficient to fit the coronal plasma with at least two temperature components, and compare abundances of groups of elements with low, medium, or high first ionization potential. Similar to the known relation between X-ray luminosity $L_X$ and rotation period, we can correlate stellar properties with coronal temperatures and abundances to constrain models for stellar activity, coronal heating, and stellar dynamos. Detailed X-ray characterization for even a subset of planet-hosting systems would dramatically advance our knowledge of what impact these emissions have on orbiting planets.

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.

For the first time, we have explored the spatial substructure of the Cygnus OB2 association using parallaxes from the recent second Gaia data release. We find significant line-of-sight substructure within the association, which we quantify using a parameterised model that reproduces the observed parallax distribution. This inference approach is necessary due to the non-linearity of the parallax-distance transformation and the asymmetry of the resulting probability distribution. By using a Markov Chain Monte Carlo ensemble sampler and an unbinned maximum likelihood test we identify two different stellar groups superposed on the association. We find the main Cygnus OB2 group at ~1760 pc, further away than recent estimates have envisaged, and a foreground group at ~1350 pc. We also calculate individual membership probabilities and identify outliers as possible non-members of the association.

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively.

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.

The formation and evolution of young star clusters and OB associations is fundamental to our understanding of the star formation process, the conditions faced by young binary and planetary systems, and the formation of long-lived open and globular clusters. Despite this our understanding of the physical processes that drive this evolution has been limited by the static nature of most observations. This is all changing thanks to a revolution in kinematic data quality from large-scale radial velocity surveys and new astrometric facilities such as Gaia. Here I summarise recent studies of multiple OB associations from both {\it Gaia} and ground-based astrometric surveys. These observations show that OB associations have considerable kinematic substructure and no evidence for the radial expansion pattern predicted by theories such as residual gas expulsion. This means that, contrary to the standard view of OB associations as expanded star clusters, these systems could never have been dense star clusters in the past and were most likely born as extended and highly substructured groups of stars. This places strong constraints on the primordial clustering of young stars and the conditions faced by young planetary systems.

The Friends of Friends algorithm identifies groups of objects with similar spatial and kinematic properties, and has recently been used extensively to quantify the distributions of gas and stars in young star-forming regions. We apply the Friends of Friends algorithm to $N$-body simulations of the dynamical evolution of subvirial (collapsing) and supervirial (expanding) star-forming regions. We find that the algorithm picks out a wide range of groups (1 -- 25) for statistically identical initial conditions, and cannot distinguish between subvirial and supervirial regions in that we obtain similar mode and median values for the number of groups it identifies. We find no correlation between the number of groups identified initially and either the initial or subsequent spatial and kinematic tracers of the regions' evolution, such as the amount of spatial substructure, dynamical mass segregation, or velocity dispersion. We therefore urge caution in using the Friends of Friends algorithm to quantify the initial conditions of star formation.

The first Gaia Data Release presents an opportunity to characterise the low-mass population of OB associations, providing larger statistical samples and better constraints on the formation and evolution of clusters and associations. Using previously known low mass members in a small region of Vela OB2 we have designed selection criteria that combine Gaia and 2MASS photometry, independently of any astrometric information, to identify low-mass pre-main-sequence (PMS) stars over the wider association area. Our method picks out the known clusters of young stars around $\gamma^2$ Velorum and NGC-2547, but also identifies other over-densities that may represent previously unknown clusters. There are clear differences in the spatial distributions of the low-mass and the high-mass OB populations, suggesting either that the structure and dynamics of these populations has evolved separately or that the initial mass function can vary considerably on small scales within a single association.

Using Gaia DR2 data, combined with OmegaCAM ground-based optical photometry from the AD-HOC survey, and detailed Radial Velocity measurements from ESO-Gaia, we analyse in detail a 10x5 deg region around the Wolf-Rayet star $\gamma^2$ Vel, including the previously known clusters Gamma Vel and NGC2547. Using clustering analysis that considers positions, proper motions and parallax, we discover 6 clusters or associations -- 4 of which appear new. Analysis of the colour-magnitude diagram for these clusters show that 4 of them formed coevally from the same molecular clouds 10 Myr ago, while NGC 2547 formed together with a newly discovered cluster 30 Myr ago. This study shows the incredible wealth of data provided by Gaia for the study of young stellar clusters.

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.

The coronal activity-rotation relationship is considered to be a proxy for the underlying stellar dynamo responsible for magnetic activity in solar and late-type stars. While this has been studied in considerable detail for partly-convective stars that are believed to operate an interface dynamo, it is poorly unconstrained in fully-convective stars that lack the necessary shear layer between radiative core and the convective envelope. We present new X-ray observations of 19 slowly-rotating fully-convective stars with rotation periods from the MEarth Project. We use these to calculate X-ray luminosities (or upper limits for undetected sources) and combine these with existing measurements from \citet{wrig16b}. We confirm the existence of fully-convective stars in the X-ray unsaturated regime and find that these objects follow the same rotation-activity relationship seen for partly-convective stars. We measure the power-law slope of the relationship between Rossby number (the ratio of the rotation period to the convective turnover time) and the fractional X-ray luminosity for X-ray unsaturated fully-convective stars for the first time, and find it to be consistent with that found for partly-convective stars. We discuss this implications of this result for our understanding of stellar magnetic dynamos in fully- and partly-convective stars. Finally, we also use this data to improve empirical estimates of the convective turnover time for fully-convective stars.

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.

The morphology and optical spectrum of IPHASXJ210205+471015, a nebula classified as a possible planetary nebula, are however strikingly similar to those of ATCnc, a classical nova shell around a dwarf nova. To investigate its true nature, we have obtained high-resolution narrow-band [O III] and [N II] images and deep GTC OSIRIS optical spectra. The nebula shows an arc of [N II]-bright knots notably enriched in nitrogen, whilst an [O III]-bright bow-shock is progressing throughout the ISM. Diagnostic line ratios indicate that shocks are associated with the arc and bow-shock. The central star of this nebula has been identified by its photometric variability. Time-resolved photometric and spectroscopic data of this source reveal a period of 4.26 hours, which is attributed to a binary system. The optical spectrum is notably similar to that of RWSex, a cataclysmic variable star (CV) of the UXUMa nova-like (NL) type. Based on these results, we propose that IPHASX J210205+471015 is a classical nova shell observed around a CV-NL system in quiescence.

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

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.

We present a kinematic study of the Scorpius-Centaurus (Sco-Cen) OB association (Sco OB2) using Gaia DR1 parallaxes and proper motions. Our goal is to test the classical theory that OB associations are the expanded remnants of dense and compact star clusters disrupted by processes such as residual gas expulsion. Gaia astrometry is available for 258 out of 433 members of the association, with revised Hipparcos astrometry used for the remainder. We use this data to confirm that the three subgroups of Sco-Cen are gravitationally unbound and have non-isotropic velocity dispersions, suggesting they have not had time to dynamically relax. We also explore the internal kinematics of the subgroups to search for evidence of expansion. We test Blaauw's classical linear model of expansion, search for velocity trends along the Galactic axes, compare the expanding and non-expanding convergence points, perform traceback analysis assuming both linear trajectories and using an epicycle approximation, and assess the evidence for expansion in proper motions corrected for virtual expansion / contraction. None of these methods provide coherent evidence for expansion of the subgroups, with no evidence to suggest that the subgroups had a more compact configuration in the past. We find evidence for kinematic substructure within the subgroups that supports the view that they were not formed by the disruption of individual star clusters. We conclude that Sco-Cen was likely born highly substructured, with multiple small-scale star formation events contributing to the overall OB association, and not as single, monolithic bursts of clustered star formation.

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.

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.

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.

In solar-type stars (with radiative cores and convective envelopes), the magnetic field powers star spots, flares and other solar phenomena, as well as chromospheric and coronal emission at ultraviolet to X-ray wavelengths. The dynamo responsible for generating the field depends on the shearing of internal magnetic fields by differential rotation. The shearing has long been thought to take place in a boundary layer known as the tachocline between the radiative core and the convective envelope. Fully convective stars do not have a tachocline and their dynamo mechanism is expected to be very different, although its exact form and physical dependencies are not known. Here we report observations of four fully convective stars whose X-ray emission correlates with their rotation periods in the same way as in Sun-like stars. As the X-ray activity - rotation relationship is a well-established proxy for the behaviour of the magnetic dynamo, these results imply that fully convective stars also operate a solar-type dynamo. The lack of a tachocline in fully convective stars therefore suggests that this is not a critical ingredient in the solar dynamo and supports models in which the dynamo originates throughout the convection zone.

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.

We present a high-precision proper motion study of 873 X-ray and spectroscopically selected stars in the massive OB association Cygnus OB2 as part of the DANCe project. These were calculated from images spanning a 15 year baseline and have typical precisions < 1 mas/yr. We calculate the velocity dispersion in the two axes to be $\sigma_\alpha(c) = 13.0^{+0.8}_{-0.7}$ and $\sigma_\delta(c) = 9.1^{+0.5}_{-0.5}$ km/s, using a 2-component, 2-dimensional model that takes into account the uncertainties on the measurements. This gives a 3-dimensional velocity dispersion of $\sigma_{3D} = 17.8 \pm 0.6$ km/s implying a virial mass significantly larger than the observed stellar mass, confirming that the association is gravitationally unbound. The association appears to be dynamically unevolved, as evidenced by considerable kinematic substructure, non-isotropic velocity dispersions and a lack of energy equipartition. The proper motions show no evidence for a global expansion pattern, with approximately the same amount of kinetic energy in expansion as there is in contraction, which argues against the association being an expanded star cluster disrupted by process such as residual gas expulsion or tidal heating. The kinematic substructures, which appear to be close to virial equilibrium and have typical masses of 40-400 M$_\odot$, also do not appear to have been affected by the expulsion of the residual gas. We conclude that Cyg OB2 was most likely born highly substructured and globally unbound, with the individual subgroups born in (or close to) virial equilibrium, and that the OB association has not experienced significant dynamical evolution since then.

Located in the Small Magellanic Cloud (SMC), SXP214 is an X-ray pulsar in a high mass X-ray binary system with a Be-star companion. A recent survey of the SMC under a Chandra X-ray Visionary program found the source in a transition when the X-ray flux was on a steady rise. The Lomb-Scargle periodogram revealed a pulse period of 211.49 +/- 0.42 s, which is significantly (>5sigma) shorter than the previous measurements with XMM-Newton and RXTE. This implies that the system has gone through sudden spin-up episodes recently. The pulse profile shows a sharp eclipse-like feature with a modulation amplitude of >95%. The linear rise of the observed X-ray luminosity from <~2x to 7x10^35 erg s^-1 is correlated with steady softening of the X-ray spectrum, which can be described by the changes in the local absorption from N_H ~ 10^24 to <~10^20 cm^-2 for an absorbed power-law model. The soft X-ray emission below 2 keV was absent in the early part of the observation when only the pulsating hard X-ray component was observed, whereas at later times both soft and hard X-ray components were observed pulsating. A likely explanation is that the neutron star was initially hidden in the circumstellar disk of the companion, and later came out of the disk with the accreted material that continued fueling the observed pulsation.

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.

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.

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.

Mass segregation in star clusters is often thought to indicate the onset of energy equipartition, where the most massive stars impart kinetic energy to the lower-mass stars and brown dwarfs/free floating planets. The predicted net result of this is that the centrally concentrated massive stars should have significantly lower velocities than fast-moving low-mass objects on the periphery of the cluster. We search for energy equipartition in initially spatially and kinematically substructured N-body simulations of star clusters with N = 1500 stars, evolved for 100 Myr. In clusters that show significant mass segregation we find no differences in the proper motions or radial velocities as a function of mass. The kinetic energies of all stars decrease as the clusters relax, but the kinetic energies of the most massive stars do not decrease faster than those of lower-mass stars. These results suggest that dynamical mass segregation -- which is observed in many star clusters -- is not a signature of energy equipartition from two-body relaxation.

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 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.

We follow the dynamical evolution of young star-forming regions with a wide range of initial conditions and examine how the radial velocity dispersion, $\sigma$, evolves over time. We compare this velocity dispersion to the theoretically expected value for the velocity dispersion if a region were in virial equilibrium, $\sigma_{\rm vir}$ and thus assess the virial state ($\sigma / \sigma_{\rm vir}$) of these systems. We find that in regions that are initially subvirial, or in global virial equilibrium but subvirial on local scales, the system relaxes to virial equilibrium within several million years, or roughly 25 - 50 crossing times, according to the measured virial ratio. However, the measured velocity dispersion, $\sigma$, appears to be a bad diagnostic of the current virial state of these systems as it suggests that they become supervirial when compared to the velocity dispersion estimated from the virial mass, $\sigma_{\rm vir}$. We suggest that this discrepancy is caused by the fact that the regions are never fully relaxed, and that the early non-equilibrium evolution is imprinted in the one-dimensional velocity dispersion at these early epochs. If measured early enough ($<$2 Myr in our simulations, or $\sim$20 crossing times), the velocity dispersion can be used to determine whether a region was highly supervirial at birth without the risk of degeneracy. We show that combining $\sigma$, or the ratio of $\sigma$ to the interquartile range (IQR) dispersion, with measures of spatial structure, places stronger constraints on the dynamical history of a region than using the velocity dispersion in isolation.

The European Space Agency's Gaia space telescope, launched in 2013, aims to measure the positions, parallaxes, and proper motions of a billion stars in our Galaxy and throughout the Local Group. In doing so it will include hundreds of thousands of young stars in star forming regions, star clusters and OB associations. This will provide us with an unprecedented view of the role of dynamics in the star formation process and the evolution of young star clusters. Data from this ambitious mission is expected very soon, with the first data release scheduled for 2016. This review discusses the current state of our understanding of the kinematics of star formation and how the community can best prepare for Gaia data.

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.

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.

The optical Galactic Plane H$\alpha$ surveys IPHAS and VPHAS+ are dramatically improving our understanding of Galactic stellar populations and stellar evolution by providing large samples of stars in short lived, but important, evolutionary phases, and high quality homogeneous photometry and images over the entire Galactic Plane. Here I summarise some of the contributions these surveys have already made to our understanding of a number of key areas of stellar and Galactic astronomy.

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<R_V<4$.

We have compiled a significantly updated and comprehensive census of massive stars in the nearby Cygnus OB2 association by gathering and homogenising data from across the literature. The census contains 169 primary OB stars, including 52 O-type stars and 3 Wolf-Rayet stars. Spectral types and photometry are used to place the stars in a Hertzprung-Russell diagram, which is compared to both non-rotating and rotating stellar evolution models, from which stellar masses and ages are calculated. The star formation history and mass function of the association are assessed, and both are found to be heavily influenced by the evolution of the most massive stars to their end states. We find that the mass function of the most massive stars is consistent with a `universal' power-law slope of Gamma = 1.3. The age distribution inferred from stellar evolutionary models with rotation and the mass function suggest the majority of star formation occurred more or less continuously between 1 and 7 Myr ago, in agreement with studies of low- and intermediate mass stars in the association. We identify a nearby young pulsar and runaway O-type star that may have originated in Cyg OB2 and suggest that the association has already seen its first supernova. Finally we use the census and mass function to calculate the total mass of the association of 16500^+3800_-2800 Msun, at the low end, but consistent with, previous estimates of the total mass of Cyg OB2. Despite this Cyg OB2 is still one of the most massive groups of young stars known in our Galaxy making it a prime target for studies of star formation on the largest scales.

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.

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.

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.

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.

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

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 .

The nature of the mechanisms apparently driving X-rays from intermediate mass stars lacking strong convection zones or massive winds remains poorly understood, and the possible role of hidden, lower mass close companions is still unclear. A 20ks Chandra HRC-I observation of HR 4796A, an 8 Myr old main sequence A0 star devoid of close stellar companions, has been used to search for a signature or remnant of magnetic activity from the Herbig Ae phase. X-rays were not detected and the X-ray luminosity upper limit was L_X =< 1.3x10^27 erg/s. The result is discussed in the context of various scenarios for generating magnetic activity, including rotational shear and subsurface convection. A dynamo driven by natal differential rotation is unlikely to produce observable X rays, chiefly because of the difficulty in getting the dissipated energy up to the surface of the star. A subsurface convection layer produced by the ionisation of helium could host a dynamo that should be effective throughout the main-sequence but can only produce X-ray luminosities of order 10^25 erg/s. This luminosity lies only moderately below the current detection limit for Vega. Our study supports the idea that X-ray production in Herbig Ae/Be stars is linked largely to the accretion process rather than the properties of the underlying star, and that early A stars generally decline in X-ray luminosity at least 100,000 fold in only a few million years.

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.

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.

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.

We examine substructure and mass segregation in the massive OB association Cygnus OB2 to better understand its initial conditions. Using a well understood Chandra X-ray selected sample of young stars we find that Cyg OB2 exhibits considerable physical substructure and has no evidence for mass segregation, both indications that the association is not dynamically evolved. Combined with previous kinematical studies we conclude that Cyg OB2 is dynamically very young, and what we observe now is very close to its initial conditions: Cyg OB2 formed as a highly substructured, unbound association with a low volume density (< 100 stars/pc^3). This is inconsistent with the idea that all stars form in dense, compact clusters. The massive stars in Cyg OB2 show no evidence for having formed particularly close to one another, nor in regions of higher than average density. Since Cyg OB2 contains stars as massive as ~100 Mo this result suggests that very massive stars can be born in relatively low-density environments. This would imply that the massive stars in Cyg OB2 did not form by competitive accretion, or by mergers.

We model the dynamical evolution of star forming regions with a wide range of initial properties. We follow the evolution of the regions' substructure using the Q-parameter, we search for dynamical mass segregation using the Lambda_MSR technique, and we also quantify the evolution of local density around stars as a function of mass using the Sigma_LDR method. The amount of dynamical mass segregation measured by Lambda_MSR is generally only significant for subvirial and virialised, substructured regions - which usually evolve to form bound clusters. The Sigma_LDR method shows that massive stars attain higher local densities than the median value in all regions, even those that are supervirial and evolve to form (unbound) associations. We also introduce the Q-Sigma_LDR plot, which describes the evolution of spatial structure as a function of mass-weighted local density in a star forming region. Initially dense (>1000 stars pc^{-2}), bound regions always have Q >1, Sigma_LDR > 2 after 5Myr, whereas dense unbound regions always have Q < 1, Sigma_LDR > 2 after 5Myr. Less dense regions (<100 stars pc^{-2}) do not usually exhibit Sigma_LDR > 2 values, and if relatively high local density around massive stars arises purely from dynamics, then the Q-Sigma_LDR plot can be used to estimate the initial density of a star forming region.

We present WISE All-Sky mid-infrared (IR) survey detections of 55% (17/31) of the known QSOs at z>6 from a range of surveys: the SDSS, the CFHT-LS, FIRST, Spitzer and UKIDSS. The WISE catalog thus provides a substantial increase in the quantity of IR data available for these sources: 17 are detected in the WISE W1 (3.4-micron) band, 16 in W2 (4.6-micron), 3 in W3 (12-micron) and 0 in W4 (22-micron). This is particularly important with Spitzer in its warm-mission phase and no faint follow-up capability at wavelengths longwards of 5 microns until the launch of JWST. WISE thus provides a useful tool for understanding QSOs found in forthcoming large-area optical/IR sky surveys, using PanSTARRS, SkyMapper, VISTA, DES and LSST. The rest-UV properties of the WISE-detected and the WISE-non-detected samples differ: the detections have brighter i/z-band magnitudes and redder rest-UV colors. This suggests that a more aggressive hunt for very-high-redshift QSOs, by combining WISE W1 and W2 data with red observed optical colors could be effective at least for a subset of dusty candidate QSOs. Stacking the WISE images of the WISE-non-detected QSOs indicates that they are on average significantly fainter than the WISE-detected examples, and are thus not narrowly missing detection in the WISE catalog. The WISE-catalog detection of three of our sample in the W3 band indicates that their mid-IR flux can be detected individually, although there is no stacked W3 detection of sources detected in W1 but not W3. Stacking analyses of WISE data for large AGN samples will be a useful tool, and high-redshift QSOs of all types will be easy targets for JWST.

Using a new uniform sample of 824 solar and late-type stars with measured X-ray luminosities and rotation periods we have studied the relationship between rotation and stellar activity that is believed to be a probe of the underlying stellar dynamo. Using an unbiased subset of the sample we calculate the power law slope of the unsaturated regime of the activity -- rotation relationship as $L_X/L_{bol}\propto Ro^\beta$, where $\beta=-2.70\pm0.13$. This is inconsistent with the canonical $\beta = -2$ slope to a confidence of 5$\sigma$ and argues for an interface-type dynamo. We map out three regimes of coronal emission as a function of stellar mass and age, using the empirical saturation threshold and theoretical super-saturation thresholds. We find that the empirical saturation timescale is well correlated with the time at which stars transition from the rapidly rotating convective sequence to the slowly rotating interface sequence in stellar spin-down models. This may be hinting at fundamental changes in the underlying stellar dynamo or internal structure. We also present the first discovery of an X-ray unsaturated, fully convective M star, which may be hinting at an underlying rotation - activity relationship in fully convective stars hitherto not observed. Finally we present early results from a blind search for stellar X-ray cycles that can place valuable constraints on the underlying ubiquity of solar-like activity cycles.

We present H\alpha images of an ionized nebula surrounding the M2-5Ia red supergiant (RSG) W26 in the massive star cluster Westerlund 1. The nebula consists of a circumstellar shell or ring ~0.1pc in diameter and a triangular nebula ~0.2pc from the star that in high-resolution Hubble Space Telescope images shows a complex filamentary structure. The excitation mechanism of both regions is unclear since RSGs are too cool to produce ionizing photons and we consider various possibilities. The presence of the nebula, high stellar luminosity and spectral variability suggest that W26 is a highly evolved RSG experiencing extreme levels of mass-loss. As the only known example of an ionized nebula surrounding a RSG W26 deserves further attention to improve our understanding of the final evolutionary stages of massive stars.

We present the result of the cross-matching between UV-excess sources selected from the UV-excess survey of the Northern Galactic Plane (UVEX) and several infrared surveys (2MASS, UKIDSS and WISE). From the position in the (J-H) vs. (H-K) colour-colour diagram we select UV-excess candidate white dwarfs with an M-dwarf type companion, candidates that might have a lower mass, brown-dwarf type companion, and candidates showing an infrared-excess only in the K-band, which might be due to a debris disk. Grids of reddened DA+dM and sdO+MS/sdB+MS model spectra are fitted to the U,g,r,i,z,J,H,K photometry in order to determine spectral types and estimate temperatures and reddening. From a sample of 964 hot candidate white dwarfs with (g-r)<0.2, the spectral energy distribution fitting shows that ~2-4% of the white dwarfs have an M-dwarf companion, ~2% have a lower-mass companion, and no clear candidates for having a debris disk are found. Additionally, from WISE 6 UV-excess sources are selected as candidate Quasi-Stellar Objects (QSOs). Two UV-excess sources have a WISE IR-excess showing up only in the mid-IR W3 band of WISE, making them candidate Luminous InfraRed Galaxies (LIRGs) or Sbc star-burst galaxies.

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.

We present a determination of the average space density and birth rate of hydrogen-line (DA) white dwarfs within a radius of 1 kpc around the Sun, based on an observational sample of 360 candidate white dwarfs with g<19.5 and (g-r)<0.4, selected from the UV-excess Survey of the Northern Galactic Plane (UVEX), in combination with a theoretical white dwarf population that has been constructed to simulate the observations, including the effects of reddening and observational selection effects. The main uncertainty in the derivation of the white dwarf space density and current birth rate lies in the absolute photometric calibration and the photometric scatter of the observational data, which influences the classification method on colours, the completeness and the pollution. Corrections for these effects are applied. We derive an average space density of hydrogen-line (DA) white dwarfs with T_eff > 10,000K (M_V<12.2) of (3.8 +/- 1.1) x 1e-4 pc^-3, and an average DA white dwarf birth rate over the last 7e7 years of (5.4 + 1.5) x 1e-13 pc^-3 yr^-1. Additionally, we show that many estimates of the white dwarf space density from different studies are consistent with each other, and with our determination here.

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.

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.

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.

Stellar activity cycles are known to be a widespread phenomenon amongst moderately active solar- and late-type stars from long-term periodic variations in chromospheric Ca II H and K emission lines, yet to date only a handful of coronal X-ray cycles are known. We have surveyed serendipitously observed stellar sources in fields observed multiple times in the last decade by XMM-Newton and present our analysis of 9 stars from 6 fields. Since our sample is flux-limited, it is strongly biased towards higher levels of X-ray activity. We fit a single temperature APEC spectrum to each source and search for significant periodicities using a Lomb-Scargle Periodogram (LSP). We use a Monte Carlo (MC) algorithm to yield robust analysis of the statistical significance of cycle detections and non-detections. None of the 9 stellar lightcurves show any convincing indications of periodicity. From MC simulations we simulate the detection capabilities of our methodology and, assuming a uniform distribution of cycle periods and strengths over the domain searched, we conclude with 95% confidence that less than 72% of the stars represented by our sample of active stars have 5-13 year coronal X-ray cycles.

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.

Using a new catalog of 824 solar and late-type stars with X-ray luminosities and rotation periods we have studied the relationship between rotation and stellar activity. From an unbiased subset of this sample the power law slope of the unsaturated regime, $L_X/L_{bol}\propto Ro^\beta$, is fit as $\beta=-2.70\pm0.13$. This is inconsistent with the canonical $\beta=-2$ slope to a confidence of 5$\sigma$ and argues for an interface-type dynamo. Super-saturation is observed for the fastest rotators in our sample and its parametric dependencies are explored. Significant correlations are found with both the corotation radius and the excess polar updraft, the latter theory being supported by other observations. We also present a new X-ray population synthesis model of the mature stellar component of our Galaxy and use it to reproduce deep observations of a high Galactic latitude field. The model, XStar, can be used to test models of stellar spin-down and dynamo decay, as well as for estimating stellar X-ray contamination rates for non-stellar studies.

IRAS 20050+2720 is young star forming region at a distance of 700 pc without apparent high mass stars. We present results of our multiwavelength study of IRAS 20050+2720 which includes observations by Chandra and Spitzer, and 2MASS and UBVRI photometry. In total, about 300 YSOs in different evolutionary stages are found. We characterize the distribution of young stellar objects (YSOs) in this region using a minimum spanning tree (MST) analysis. We newly identify a second cluster core, which consists mostly of class II objects, about 10 arcmin from the center of the cloud. YSOs of earlier evolutionary stages are more clustered than more evolved objects. The X-ray luminosity function (XLF) of IRAS 20050+2720 is roughly lognormal, but steeper than the XLF of the more massive Orion nebula complex. IRAS 20050+2720 shows a lower N_H/A_K ratio compared with the diffuse ISM.

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.

We report the spectral and temporal variability properties of 18 candidate transient and potential transient (TC and PTC) sources detected in deep multi-epoch Chandra observation of the nearby elliptical galaxies, NGC 3379, NGC 4278 and NGC 4697. Only one source can be identified with a background counterpart, leaving 17 TCs + PTCs in the galaxies. Of these, 14 are in the galaxy field, supporting the theoretical picture that the majority of field X-ray binaries (XRBs) will exhibit transient accretion for >75% of their lifetime. Three sources are coincident with globular clusters (GCs), including two high-luminosity candidate black hole (BH) XRBs, with Lx=5.4E38 erg/s, and Lx=2.8E39 erg/s, respectively. The spectra, luminosities and temporal behavior of these 17 sources suggest that the transient population is heterogeneous, including neutron star (NS) and BH XRBs in both normal and high-rate accretion modes, and super soft sources containing white dwarf binaries. Our TC and PTC detections are noticeably fewer that the number expected from the populations synthesis (PS) models of Fragos et al. (2009), tailored to our new Chandra pointings of NGC 4278. We attribute this discrepancy to the PS assumption that the transient population is composed of NS XRBs, as well as differences between the statistical analysis and error estimates used in the model and our observations.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra program that has imaged the central 0.9 deg^2 of the COSMOS field down to limiting depths of 1.9 10^-16 erg cm^-2 s-1 in the 0.5-2 keV band, 7.3 10^-16 erg cm^-2 s^-1 in the 2-10 keV band, and 5.7 10^-16 erg cm^-2 s-1 in the 0.5-10 keV band. In this paper we report the i, K and 3.6micron identifications of the 1761 X-ray point sources. We use the likelihood ratio technique to derive the association of optical/infrared counterparts for 97% of the X-ray sources. For most of the remaining 3%, the presence of multiple counterparts or the faintness of the possible counterpart prevented a unique association. For only 10 X-ray sources we were not able to associate a counterpart, mostly due to the presence of a very bright field source close by. Only 2 sources are truly empty fields. Making use of the large number of X-ray sources, we update the "classic locus" of AGN and define a new locus containing 90% of the AGN in the survey with full band luminosity >10^42 erg/s. We present the linear fit between the total i band magnitude and the X-ray flux in the soft and hard band, drawn over 2 orders of magnitude in X-ray flux, obtained using the combined C-COSMOS and XMM-COSMOS samples. We focus on the X-ray to optical flux ratio (X/O) and we test its known correlation with redshift and luminosity, and a recently introduced anti-correlation with the concentration index (C). We find a strong anti-correlation (though the dispersion is of the order of 0.5 dex) between C and X/O, computed in the hard band, and that 90% of the obscured AGN in the sample with morphological information live in galaxies with regular morphology (bulgy and disky/spiral), suggesting that secular processes govern a significant fraction of the BH growth at X-ray luminosities of 10^43- 10^44.5 erg/s.

This paper describes the first data release of the Kepler-INT Survey (KIS), that covers a 116 deg2 region of the Cygnus and Lyra constellations. The Kepler field is the target of the most intensive search for transiting planets to date. Despite the fact that the Kepler mission provides superior time series photometry, with an enormous impact on all areas of stellar variability, its field lacks optical photometry complete to the confusion limit of the Kepler instrument necessary for selecting various classes of targets. For this reason, we follow the observing strategy and data reduction method used in the IPHAS and UVEX galactic plane surveys in order to produce a deep optical survey of the Kepler field. This initial release concerns data taken between May and August 2011, using the Isaac Newton Telescope on the island of La Palma. Four broadband filters were used, U, g, r, i, as well as one narrowband one, Halpha, reaching down to a 10-sigma limit of around 20th mag in the Vega system. Observations covering around 50 deg2, thus about half of the field, passed our quality control thresholds and constitute this first data release. We derive a global photometric calibration by placing the KIS magnitudes as close as possible to the Kepler Input Catalog (KIC) photometry. The initial data release catalogue containing around 6 million sources from all the good photometric fields is available for download from the KIS webpage, as well as via MAST.

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 sample of 824 solar and late-type stars with X-ray luminosities and rotation periods. This is used to study the relationship between rotation and stellar activity and derive a new estimate of the convective turnover time. From an unbiased subset of this sample the power law slope of the unsaturated regime, L_X / L_bol = Ro^\beta, is fit as \beta = -2.70 +/- 0.13. This is inconsistent with the canonical \beta=-2 slope to a confidence of 5 sigma, and argues for an additional term in the dynamo number equation. From a simple scaling analysis this implies \Delta\Omega / \Omega = \Omega^0.7, i.e. the differential rotation of solar-type stars gradually declines as they spin down. Super-saturation is observed for the fastest rotators in our sample and its parametric dependencies are explored. Significant correlations are found with both the corotation radius and the excess polar updraft, the latter theory providing a stronger dependence and being supported by other observations. We estimate mass-dependent empirical thresholds for saturation and super-saturation and map out three regimes of coronal emission. Late F-type stars are shown never to pass through the saturated regime, passing straight from super-saturated to unsaturated X-ray emission. The theoretical threshold for coronal stripping is shown to be significantly different from the empirical saturation threshold (Ro < 0.13), suggesting it is not responsible. Instead we suggest that a different dynamo configuration is at work in stars with saturated coronal emission. This is supported by a correlation between the empirical saturation threshold and the time when stars transition between convective and interface sequences in rotational spin-down models.

V1647 Ori is a young eruptive variable star, illuminating a reflection nebula (McNeil's Nebula). It underwent an outburst in 2003 before fading back to its pre-outburst brightness in 2006. In 2008, V1647 Ori underwent a new outburst. The observed properties of the 2003-2006 event are different in several respects from both the EXor and FUor type outbursts, and suggest that this star might represent a new class of eruptive young stars, younger and more deeply embedded than EXors, and exhibiting variations on shorter time scales than FUors. In outburst, the star lights up the otherwise invisible McNeil's nebular - a conical cloud likely accumulated from previous outbursts. We present follow-up photometric as well as optical and near-IR spectroscopy of the nebula obtainted during the 2008-2009 outburst. We will also present results from contemporaneous X-ray observations. These multi-wavelength observations of V1647 Ori, obtained at this key early stage of the outburst, provide a snapshot of the "lighting up" of the nebula, probe its evolution through the event, and enable comparison with the 2003-2006 outburst.

We present a 3D photoionization model of the PN NGC 6302, one of the most complex objects of its kind. Our Mocassin model is composed of an extremely dense circumstellar disk and a large pair of diffuse bipolar lobes, a combination necessary to reproduce the observed spectrum. The masses of these components gives a total nebular mass of 4.7Mo. Discrepancies between our model fit and the observations are attributed to complex density inhomogeneities in the nebula. The potential to resolve such discrepancies with more complex models is confirmed by a range of models introducing small-scale structures. Compared to solar abundances He is enhanced by 50%, C is slightly subsolar, O is solar, and N is enhanced by a factor of 6. These imply a significant 3rd dredge-up coupled with hot-bottom burning CN-cycle conversion of dredged-up C to N. The central star is partly obscured by the edge-on circumstellar disk and its properties are not well constrained. Emission from a number of high-ionization `coronal' lines provides constraints on the form of the high-energy ionizing flux. Using a solar abundance stellar atmosphere we are unable to fit all of the observed line fluxes, but a substantially better fit was obtained using a 220,000K H-deficient stellar atmosphere with L*=14,300 Lo. The H-deficient nature of the central star suggests it has undergone a late thermal pulse, and fits to evolutionary tracks imply a central star mass of 0.73-0.82Mo. Timescales for these tracks suggest the object left the top of the AGB ~2100 years ago, in agreement with studies of the recent mass-loss event that formed the bipolar lobes. Based on the modelled nebular and central star masses we estimate the initial mass of the central star to be 5.5Mo, in agreement with that derived from evolutionary tracks. (Abstract truncated)

The INT Photometric H-Alpha Survey (IPHAS) is a 1800 deg^2 survey of the Northern Galactic Plane, reaching down to r'~21. We demonstrate how the survey can be used to (1) reliably select classical T Tauri star candidates and (2) constrain the mass accretion rates with an estimated relative uncertainty of 0.6 dex. IPHAS is a necessary addition to spectroscopic surveys because it allows large and uniform samples of accretion rates to be obtained with a precise handle on the selection effects. We apply the method on a region of 7 deg^2 towards the HII region IC1396 in Cepheus OB2 and identify 158 pre-main sequence candidates with masses between 0.2 and 2.0 Msun and accretion rates between 10^-9.2 and 10^-7.0 Msun yr^-1. We find a power-law dependency between the stellar mass and the accretion rates with a slope of \alpha = 1.1 +/- 0.2, which is less steep than indicated by previous studies. We discuss the influence of method-dependent systematic effects on the determination of this relationship. The majority of our sample consists of faint, previously unknown, low-mass T Tauri candidates (56 per cent between 0.2 and 0.5 Msun). Many candidates are clustered in front of three bright-rimmed molecular clouds, which are being ionized by the massive star HD206267 (O6.5V). We discover a spatio-temporal gradient of increasing accretion rates, increasing Spitzer infrared excess, and younger ages away from the ionizing star, providing a strong indication that the formation of these clusters has been sequentially triggered by HD206267 during the last ~1 Myr.

Cygnus OB2 is the nearest example of a massive star forming region, containing over 50 O-type stars and hundreds of B-type stars. We have analyzed two Chandra pointings in Cyg OB2, detecting ~1700 X-ray sources, of which ~1450 are thought to be members of the association. Optical and near-IR photometry has been obtained for ~90% of these sources from recent deep Galactic plane surveys. We have performed isochrone fits to the near-IR color-magnitude diagram, deriving ages of 3.5(+0.75,-1.0) and 5.25(+1.5,-1.0) Myrs for sources in the two fields, both with considerable spreads around the pre-MS isochrones. The presence of a second population in the region, somewhat older than the present-day O-type stars, has been suggested by other authors and fits with the ages derived here. The fraction of sources with inner circumstellar disks (as traced by the K-band excess) is found to be very low, but appropriate for a population of age ~5 Myrs. We measure the stellar mass functions and find a power-law slope of Gamma = -1.09 +/- 0.13, in good agreement with the global mean value estimated by Kroupa. A steepening of the mass function at high masses is observed and we suggest this is due to the presence of the previous generation of stars that have lost their most massive members. Finally, combining our mass function and an estimate of the radial density profile of the association suggests a total mass of Cyg OB2 of ~30,000 Msun, similar to that of many of our Galaxy's most massive star forming regions.

Existing stellar X-ray surveys suggest major problems in our understanding of the evolution of stellar magnetic activity in solar and late-type stars, reaching conflicting conclusions over the rate of decay of X-ray activity and the spectral types responsible. We are confronting these discrepancies with a new model of the expected stellar X-ray luminosity distribution, combining a Galactic population synthesis model with current theories for rotational spin- down and the rotation - activity relation for the stellar magnetic dynamo. Here we test our model using new observations of the stellar content of the Chandra COSMOS survey, for which 60 new stellar X-ray sources are identified from the thin disk and Galactic halo populations. Our model is in approximate agreement with the observed X-ray luminosity distribution and the distribution of spectral types responsible. However, slight differences in the form of the X-ray luminosity distribution exist that may hint at problems in our understanding of stellar X-ray emission.

The supersonic stellar and disk winds possessed by massive young stellar objects will produce shocks when they collide against the interior of a pre-existing bipolar cavity (resulting from an earlier phase of jet activity). The shock heated gas emits thermal X-rays which may be observable by spaceborne observa- tories such as the Chandra X-ray Observatory. Hydrodynamical models are used to explore the wind-cavity interaction. Radiative transfer calculations are performed on the simulation output to produce synthetic X-ray observations, allowing constraints to be placed on model parameters through comparisons with observations. The model reveals an intricate interplay between the inflowing and outflowing material and is successful in reproducing the observed X-ray count rates from massive young stellar objects.

We present an analysis of the X-ray properties of a sample of solar- and late-type field stars identified in the Chandra Cosmic Evolution Survey (COSMOS), a deep (160ks) and wide (0.9 deg2) extragalactic survey. The sample of 60 sources was identified using both morphological and photometric star/galaxy separation methods. We determine X-ray count rates, extract spectra and light curves and perform spectral fits to determine fluxes and plasma temperatures. Complementary optical and near-IR photometry is also presented and combined with spectroscopy for 48 of the sources to determine spectral types and distances for the sample. We find distances ranging from 30pc to ~12kpc, including a number of the most distant and highly active stellar X-ray sources ever detected. This stellar sample extends the known coverage of the L_X-distance plane to greater distances and higher luminosities, but we do not detect as many intrinsically faint X-ray sources compared to previous surveys. Overall the sample is typically more luminous than the active Sun, representing the high-luminosity end of the disk and halo X-ray luminosity functions. The halo population appears to include both low-activity spectrally hard sources that may be emitting through thermal bremsstrahlung, as well as a number of highly active sources in close binaries.

Cygnus OB2 is the nearest example of a massive star forming region, containing over 50 O-type stars and hundreds of B-type stars. We have analysed the properties of young stars in two fields in Cyg OB2 using the recently published deep catalogue of Chandra X-ray point sources with complementary optical and near-IR photometry. Our sample is complete to 1 Msun (excluding A and B-type stars that do not emit X-rays), making this the deepest study of the stellar properties and star formation history in Cyg OB2 to date. From Siess et al. (2000) isochrone fits to the near-IR color-magnitude diagram, we derive ages of 3.5 (+0.75/-1.0) and 5.25 (+1.5/-1.0) Myrs for sources in the two fields, both with considerable spreads around the pre-MS isochrones. The presence of a stellar population somewhat older than the present-day O-type stars, also fits in with the low fraction of sources with inner circumstellar disks (as traced by the K-band excess) that we find to be very low, but appropriate for a population of age ~5 Myrs. We also find that the region lacks a population of highly embedded sources that is often observed in young star forming regions, suggesting star formation in the vicinity has declined. We measure the stellar mass functions in this limit and find a power-law slope of Gamma = -1.09 +/- 0.13, in good agreement with the global mean value estimated by Kroupa (2002). A steepening of the slope at higher masses is observed and suggested as due to the presence of the previous generation of stars that have lost their most massive members. Finally, combining our mass function and an estimate of the radial density profile of the association suggests a total mass of Cyg OB2 of ~3 x 10^4 Msun, similar to that of many of our Galaxy's most massive star forming regions.

We present the results of the search for candidate Planetary Nebulae interacting with the interstellar medium (PN-ISM) in the framework of the INT Photometric H$\alpha$ Survey (IPHAS) and located in the right ascension range 18h-20h. The detection capability of this new Northern survey, in terms of depth and imaging resolution, has allowed us to overcome the detection problem generally associated to the low surface brightness inherent to PNe-ISM. We discuss the detection of 21 IPHAS PN-ISM candidates. Thus, different stages of interaction were observed, implying various morphologies i.e. from the unaffected to totally disrupted shapes. The majority of the sources belong to the so-called WZO2 stage which main characteristic is a brightening of the nebula's shell in the direction of motion. The new findings are encouraging as they would be a first step into the reduction of the scarcity of observational data and they would provide new insights into the physical processes occurring in the rather evolved PNe.

The INT Photometric H-Alpha Survey (IPHAS) is particular effective in tracing asymptotic giant branch (AGB) stars across the Galactic disk due to its use of broad-band r' and i' filters, and the clear separation of the dwarf and giant branches in the colour-colour diagram. Here we discuss some results from our studies of the Galactic AGB population that has been revealed by IPHAS. This includes a photometric study, our spectroscopic follow-up campaign, and the photometric technique that allows us to separate O-rich, S-type and carbon stars across the entire Galactic disk based on their IPHAS photometry alone. We also discuss the potential impact of this technique on our understanding of the dredge-up process in AGB stars, Galactic chemical evolution, and the structure and properties of the Galactic disk.

We analyzed data accumulated during 2005 and 2006 by the Taiwan-American Occultation Survey (TAOS) in order to detect short-period variable stars (periods of <~ 1 hour) such as delta Scuti. TAOS is designed for the detection of stellar occultation by small-size Kuiper Belt Objects (KBOs) and is operating four 50cm telescopes at an effective cadence of 5Hz. The four telescopes simultaneously monitor the same patch of the sky in order to reduce false positives. To detect short-period variables, we used the Fast Fourier Transform algorithm (FFT) inasmuch as the data points in TAOS light-curves are evenly spaced. Using FFT, we found 41 short-period variables with amplitudes smaller than a few hundredths of a magnitude and periods of about an hour, which suggest that they are low-amplitude delta Scuti stars (LADS). The light-curves of TAOS delta Scuti stars are accessible online at the Time Series Center website (http://timemachine.iic.harvard.edu)

We present results from a catalogue of 1696 X-ray point sources detected in the massive star forming region (SFR) Cygnus OB2, the majority of which have optical or near-IR associations. We derive ages of 3.5 and 5.25 Myrs for the stellar populations in our two fields, in agreement with recent studies that suggest that the central 1-3 Myr OB association is surrounded and contaminated by an older population with an age of 5-10 Myrs. The fraction of sources with proto-planetary disks, as traced by K-band excesses, are unusually low. Though this has previously been interpreted as due to the influence of the large number of OB stars in Cyg OB2, contamination from an older population of stars in the region could also be responsible. An initial mass function is derived and found to have a slope of Gamma = -1.27, in agreement with the canonical value. Finally we introduce the recently approved Chandra Cygnus OB2 Legacy Survey that will image a 1 square degree area of the Cygnus OB2 association to a depth of 120 ks, likely detecting ~10,000 stellar X-ray sources.

We present a library of 139 near-IR spectra of cool asymptotic giant branch stars that will be useful for comparison with theoretical model atmosphere calculations and for modeling the integrated emission from intermediate-age stellar populations. The source list was selected from the `extremely red' region of the INT Photometric H Alpha Survey (IPHAS) colour-colour plane that is overwhelmingly dominated by very late-type stars. The spectral library also includes a large fraction of S-type and carbon stars. We present a number of spectral classification sequences highlighting the various molecular features identified and discuss a number of rare features with uncertain identifications in the literature. With its focus on particularly cool photospheres this catalogue serves as a companion to recent spectroscopic atlases of MK standards in the near-IR. Finally the relationship between IPHAS (r'-i') and (r'-H Alpha) colours and spectroscopically determined properties is discussed and a strong correlation between (r'-H Alpha) colour and the C/O abundance index for S-type and carbon stars is noted. This relation has the potential to separate O-rich, S-type and carbon stars in the Galaxy based on their photometry alone.

We present a catalogue of 1696 X-ray sources detected in the massive star forming region (SFR) Cygnus OB2 and extracted from two archival Chandra observations of the center of the region. A deep source extraction routine, exploiting the low background rates of Chandra observations was employed to maximize the number of sources extracted. Observations at other wavelengths were used to identify low count-rate sources and remove likely spurious sources. Monte Carlo simulations were also used to assess the authenticity of these sources. We used a Bayesian technique to identify optical or near-IR counterparts for 1501 (89%) of our sources, using deep observations from the INT Photometric H-alpha Survey, the Two Micron All Sky Survey, and the UKIRT Infrared Deep Sky Survey-Galactic plane Survey. 755 (45%) of these objects have six-band r', H-alpha, i', J, H, and K optical and near-IR photometry. From an analysis of the Poisson false-source probabilities for each source we estimate that our X-ray catalogue includes < 1% of false sources, and an even lower fraction when only sources with optical or near-IR associations are considered. A Monte Carlo simulation of the Bayesian matching scheme allows this method to be compared to more simplified matching techniques and enables the various sources of error to be quantified. The catalogue of 1696 objects presented here includes X-ray broad band fluxes, model fits, and optical and near-IR photometry in what is one of the largest X-ray catalogue of a single SFR to date. The high number of stellar X-ray sources detected from relatively shallow observations confirms the status and importance of Cygnus OB2 as one of our Galaxy's most massive SFRs.

2D axis-symmetric hydrodynamical simulations are presented which explore the interaction of stellar and disk winds with surrounding infalling cloud material. The star, and its accompanying disk, blow winds inside a cavity cleared out by an earlier jet. The collision of the winds with their surroundings generates shock heated plasma which reaches temperatures up to ~10^8 K. Attenuated X-ray spectra are calculated from solving the equation of radiative transfer along lines-of-sight. This process is repeated at various epochs throughout the simulations to examine the evolution of the intrinsic and attenuated flux. We find that the dynamic nature of the wind-cavity interaction fuels intrinsic variability in the observed emission on timescales of several hundred years. This is principally due to variations in the position of the reverse shock which is influenced by changes in the shape of the cavity wall. The collision of the winds with the cavity wall can cause clumps of cloud material to be stripped away. Mixing of these clumps into the winds mass-loads the flow and enhances the X-ray emission measure. The position and shape of the reverse shock plays a key role in determining the strength and hardness of the X-ray emission. In some models the reverse shock is oblique to much of the stellar and disk outflows, whereas in others it is closely normal over a wide range of polar angles. For reasonable stellar and disk wind parameters the integrated count rate and spatial extent of the intensity peak for X-ray emission agree with \textit{Chandra} observations of the deeply embedded MYSOs S106 IRS4, Mon R2 IRS3 A, and AFGL 2591.(abridged)

We present a deep and wide field-of-view (4'x 7') image of the planetary nebula (PN) NGC 7293 (the Helix Nebula) in the 2.12 micron H2 v=1-0 S(1) line. The excellent seeing (0.4'') at the Subaru Telescope, allows the details of cometary knots to be examined. The knots are found at distances of 2.2'-6.4' from the central star (CS). At the inner edge and in the inner ring (up to 4.5' fromthe CS), the knot often show a `tadpole' shape, an elliptical head with a bright crescent inside and a long tail opposite to the CS. In detail, there are variations in the tadpole shapes, such as narrowing tails, widening tails, meandering tails, or multi-peaks within a tail. In the outer ring (4.5'-6.4' from the CS), the shapes are more fractured, and the tails do not collimate into a single direction. The transition in knot morphology from the inner edge to the outer ring is clearly seen. The number density of knots governs the H2 surface brightness in the inner ring: H2 exists only within the knots. Possible mechanisms which contribute to the shaping of the knots are discussed, including photo-ionization and streaming motions. A plausible interpretation of our images is that inner knots are being overrun by a faster wind, but that this has not (yet) reached the outer knots. Based on H2 formation and destruction rates, H2 gas can survive in knots from formation during the late asymptotic giant branch (AGB) phase throughout the PN phase. These observations provide new constraints on the formation and evolution of knots, and on the physics of molecular gas embedded within ionized gas.

We present photometric analysis and follow-up spectroscopy for a population of extremely red stellar objects extracted from the point-source catalogue of the INT Photometric H-Alpha Survey (IPHAS) of the northern galactic plane. The vast majority of these objects have no previous identification. Analysis of optical, near- and mid-infrared photometry reveals that they are mostly highly-reddened asymptotic giant branch stars, with significant levels of circumstellar material. We show that the distribution of these objects traces galactic extinction, their highly reddened colours being a product of both interstellar and circumstellar reddening. This is the first time that such a large sample of evolved low-mass stars has been detected in the visual and allows optical counterparts to be associated with sources from recent infrared surveys. Follow-up spectroscopy on some of the most interesting objects in the sample has found significant numbers of S-type stars which can be clearly separated from oxygen-rich objects in the IPHAS colour-colour diagram. We show that this is due to the positions of different molecular bands relative to the narrow-band H-alpha filter used for IPHAS observations. The IPHAS (r' - H-alpha) colour offers a valuable diagnostic for identifying S-type stars. A selection method for identifying S-type stars in the galactic plane is briefly discussed and we estimate that over a thousand new objects of this type may be discovered, potentially doubling the number of known objects in this short but important evolutionary phase.

Nova V458 Vul erupted on 2007 August 8th and reached a visual magnitude of 8.1 a few days later. H$\alpha$ images obtained six weeks before the outburst as part of the IPHAS galactic plane survey reveal an 18th magnitude progenitor surrounded by an extended nebula. Subsequent images and spectroscopy of the nebula reveal an inner nebular knot increasing rapidly in brightness due to flash ionization by the nova event. We derive a distance of 13 kpc based on light travel time considerations, which is supported by two other distance estimation methods. The nebula has an ionized mass of 0.2 M$_{\odot}$ and a low expansion velocity: this rules it out as ejecta from a previous nova eruption, and is consistent with it being a ~14,000 year old planetary nebula, probably the product of a prior common envelope (CE) phase of evolution of the binary system. The large derived distance means that the mass of the erupting WD component of the binary is high. We identify two possible evolutionary scenarios, in at least one of which the system is massive enough to produce a Type Ia supernova on merging.

The INT/WFC Photometric H-alpha Survey of the Northern Galactic Plane (IPHAS) is an imaging survey being carried out in H-alpha, r' and i' filters, with the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) to a depth of r'=20 (10 sigma). The survey is aimed at revealing large scale structure in our local galaxy, and also the properties of key early and late populations making up the Milky Way. Mapping emission line objects enables a particular focus on objects in the young and old stages of stellar evolution ranging from early T-Tauri stars to late planetary nebulae. In this paper we present the IPHAS Initial Data Release, primarily a photometric catalogue of about 200 million unique objects, coupled with associated image data covering about 1600 square degrees in three passbands. We note how access to the primary data products has been implemented through use of standard virtual observatory publishing interfaces. Simple traditional web access is provided to the main IPHAS photometric catalogue, in addition to a number of common catalogues (such as 2MASS) which are of immediate relevance. Access through the AstroGrid VO Desktop opens up the full range of analysis options, and allows full integration with the wider range of data and services available through the Virtual Observatory. The IDR represents the largest dataset published primarily through VO interfaces to date, and so stands as an examplar of the future of survey data mining. Examples of data access are given, including a cross-matching of IPHAS photometry with sources in the UKIDSS Galactic Plane Survey that validates the existing calibration of the best data.

We report on the discovery of over 50 strong Halpha emitting objects towards the large OB association Cyg OB2 and the HII region DR 15 on its southern periphery. This was achieved using the INT Photometric Halpha Survey of the Northern Galactic Plane (IPHAS), combined with follow-up spectroscopy using the MMT multi-object spectrometer HectoSpec. We present optical spectra, supplemented with optical r', i' and Halpha photometry from IPHAS, and near-infrared J, H, and K photometry from 2MASS. The position of the objects in the (J - H) versus (H - K) diagram strongly suggests most of them are young. Many show CaII IR triplet emission indicating that they are in a pre-main sequence phase of evolution of T Tauri and Herbig Ae nature. Among these, we have uncovered pronounced clustering of T Tauri stars roughly a degree south of the centre of Cyg OB2, in an arc close to the HII region DR 15, and the radio ring nebula G79.29+0.46, for which we discuss its candidacy as a luminous blue variable (LBV). The emission line objects toward Cyg OB2 itself could be the brightest most prominent component of a population of lower mass pre-main sequence stars that has yet to be uncovered. Finally, we discuss the nature of the ongoing star formation in Cyg OB2 and the possibility that the central OB stars have triggered star formation in the periphery.

We present a three-dimensional photoionisation and dust radiative transfer model of NGC 6302, an extreme, high-excitation planetary nebula. We use the 3D photoionisation code Mocassin} to model the emission from the gas and dust. We have produced a good fit to the optical emission-line spectrum, from which we derived a density distribution for the nebula. A fit to the infrared coronal lines places strong constraints on the properties of the unseen ionising source. We find the best fit comes from using a 220,000 K hydrogen-deficient central star model atmosphere, indicating that the central star of this PN may have undergone a late thermal pulse. We have also fitted the overall shape of the ISO spectrum of NGC 6302 using a dust model with a shallow power-law size distribution and grains up to 1.0 micron in size. To obtain a good fit to the infrared SED the dust must be sufficiently recessed within the circumstellar disk to prevent large amounts of hot dust at short wavelengths, a region where the ISO spectrum is particularly lacking. These and other discoveries are helping to unveil many properties of this extreme object and trace it's evolutionary history.

Sh 2-188 is an example of strong interaction between a planetary nebula (PN) and the interstellar medium (ISM). Its structure is postulated to be the result of motion through the ISM. We present new H$\alpha$ images from the Isaac Newton Telescope Photometric H$\alpha$ Survey of the Northern Galactic Plane which reveal new structure. The nebula extends 15 arcmin on the sky in total. We have developed a `triple-wind' hydrodynamical model, comprising of the initial `slow' asymptotic giant branch (AGB) wind and the later `fast' stellar wind plus a third wind reflecting the motion through the ISM. Simulations at various velocities of the central star relative to the ISM indicate that a high velocity of 125 kms is required to reproduce the observed structure. We find most of the structure already forms during the AGB phase. The closure of the ring arises from the slow--fast wind interaction. Most of the mass lost on the AGB has been swept downstream, providing a potential explanation of the missing mass problem in PNe. We report a proper motion for the central star of 30 mas yr^-1 in the direction of the bright limb, implying a distance to the nebula of 850^{+500}_{-420} pc, consistent with a spectroscopic distances. Expansion velocities measured from spectroscopic data are consistent with velocities measured from the simulation. The model shows that the size of the PN was already set during the AGB phase.

JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision. JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each. JEDI will measure the cosmic expansion history H(z) as a free function to < 2% accuracy in redshift bins of 0.2-0.3. Assuming a flat universe and \sigma(\Omega_m)=0.01 (0.03), JEDI could measure a departure from a vanilla \Lambda CDM model (w_0=-1, w'=0) with \sigma(w_0)=0.013 (0.031) and \sigma(w')=0.046 (0.063). JEDI will obtain the well-sampled lightcurves in Z, J, H, K, L bands and spectra of ~ 14,000 type Ia supernovae with redshifts ranging from 0 to 1.7; the redshifts of ~ 10-100 million galaxies to H ~ 23 and z ~ 4 over 1000-10,000 square degrees; and measurements of the shapes of galaxies over 1000-10,000 square degrees in Z,J,H,K,L for > 10^9 galaxies to H \~ 25.

The blue supergiant Sher25 in the massive Galactic cluster NGC3603 is surrounded by a striking emission line nebula. The nebula contains an equatorial ring and probable bi-polar outflows, and is similar in morphology, mass and kinematics to the structure visible around SN1987A. It has been suggested that both nebulae were ejected while Sher25 and the progenitor of SN1987A were in previous red supergiant phases. We present optical high-resolution spectra of Sher25 and a model photosphere and unified stellar wind analysis which determines atmospheric parameters, mass-loss rate and photospheric abundances. We compare CNO abundances to other Galactic B-type supergiants and find that Sher25 does not appear extreme or abnormal in terms of its photospheric nitrogen and helium abundances. The C/N and N/O ratios are compared to surface abundances predicted by stellar evolutionary calculations and are incompatible with the star having a previous red-supergiant phase. The nebula is likely to have been ejected while the star was a blue supergiant. The results are compatible with some degree of rotationally induced mixing having occurred while the star was on or near the main-sequence. Our analsysis suggests the star has a relatively normal stellar wind and mass-loss rate, and sits comfortably within the wind momentum-luminosity relationship. In light of the evidence regarding massive evolved early-type stars in the Galaxy we suggest there is no object which shows clear evidence of having had a previous red supergiant phase and hence of undergoing blue loops in the HR diagram. [ABRIDGED]