Search

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-South survey achieved a root mean square noise level of $\sim$ 0.11 mJy beam$^{-1}$, using the 6A configuration of the Australia Telescope Compact Array (ATCA). In this paper, we discuss the observations, data processing and measurements of the source properties. Above a 7$\sigma$ detection limit, 4701 sources were detected, and their ensemble properties show similar distributions with their northern counterparts. The catalogue is highly reliable and is complete to 90 per cent at a flux density level of 1.1 mJy. We developed a new way of measuring the integrated flux densities and angular sizes of non-Gaussian sources. The catalogue primarily provides positions, flux density measurements and angular sizes. All sources with IR counterparts at 8$\mu m$ have been visually classified, utilizing additional imaging data from optical, near-IR, mid-IR, far-IR and sub-millimetre galactic plane surveys. This has resulted in the detection of 524 H II regions of which 255 are ultra-compact H II regions, 287 planetary nebulae, 79 radio stars and 6 massive young stellar objects. The rest of the sources are likely to be extra-galactic. These data are particularly important in the characterization and population studies of compact ionized sources such as UCHII regions and PNe towards the Galactic mid-plane.

With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We present the analysis of K-band spectro-interferometric observations toward the Massive Young Stellar Object IRAS 13481-6124, which is known to host an accreting dusty disc. Using GRAVITY on the VLTI, we trace the crucial au-scales of the warm inner interface between the star and the accretion dusty disc. We detect and spatially resolve the Na I doublet and He I transitions towards an object of this class for the first time. The new observations in combination with our geometric models allowed us to probe the smallest au-scales of accretion/ejection around an MYSO. We find that Na I originates in the disc at smaller radii than the dust disc and is more compact than any of the other spatially resolved diagnostics (Br$\gamma$, He I, and CO). Our findings suggest that Na I can be a new powerful diagnostic line in tracing the warm star/disc accreting interface of forming (massive) stars, while the similarities between He I and Br$\gamma$ point towards an accretion/ejection origin of He I, Comment: Accepted for publication in MNRAS Letters; 5 pages, 3 Figures

Young massive stars influence their surroundings from local to galactic scales, but the observational challenges associated with their distance and embedded nature has, until the recent decade, made high-resolution studies of these objects difficult. In particular, comparative analyses of massive young stellar object (MYSO) discs are currently lacking and our understanding of their evolution is limited. Here, we combine the results of two studies with the aim to attribute geometrical features to an evolutionary sequence for a sample of seven MYSOs. The time evolution is based on a near-IR spectral features, while the geometry is determined from a multi size-scale study of MYSOs. We find that MYSO discs with determined geometrical substructure turn out to be also spectroscopically more evolved. This implies that disc evolution and dispersal are occurring within MYSOs, similar to low-mass YSO disc evolution, despite their faster formation timescales.

In conjunction with a previous southern-hemisphere work, we present the largest radio survey of jets from massive protostars to date with high-resolution, ($\sim 0.04^{\prime\prime}$) VLA observations towards two subsamples of massive star-forming regions of different evolutionary statuses: 48 infrared-bright, massive YSOs and 8 IRDCs containing 16 luminous (${\rm L_{bol}}>10^3\,{\rm L_\odot}$) cores. For $94\%$ of the MYSO sample we detect thermal radio ($\alpha \geq -0.1$ whereby $S_\nu \propto \nu^\alpha$) sources coincident with the protostar, of which $84\%$ (13 jets and 25 candidates) are jet-like. Radio luminosity is found to scale with ${\rm L_{bol}}$ similarly to the low-mass case supporting a common mechanism for jet production across all masses. Associated radio lobes tracing shocks are seen towards $52\%$ of jet-like objects and are preferentially detected towards jets of higher radio and bolometric luminosities, resulting from our sensitivity limitations. We find jet mass loss rate scales with bolometric luminosity as $\dot{m}_{\rm jet} \propto {\rm L_{bol}}^{0.9\pm0.2}$, thereby discarding radiative, line-driving mechanisms as the dominant jet-launching process. Calculated momenta show that the majority of jets are mechanically capable of driving the massive, molecular outflow phenomena since $p_{\rm jet}>p_{\rm outflow}$. Finally, from their physical extent we show that the radio emission can not originate from small, optically-thick \textsc{Hii} regions. Towards the IRDC cores, we observe increasing incidence rates/radio fluxes with age using the proxy of increasing luminosity-to-mass $\left( \frac{L}{M} \right)$ and decreasing infrared flux ratios $\left(\frac{S_{70 \rm \mu m}}{S_{24 \rm \mu m}}\right)$. Cores with $\frac{L}{M}<40\,{\rm L_{sol}}{\rm M_{sol}}^{-1}$ are not detected above ($5.8{\rm GHz}$) radio luminosities of $\sim1{\rm mJy\,kpc}^2$., Comment: To be published in MNRAS. 20 pages main text (4 tables, 14 figures), 90 pages appendices (6 tables, 66 figures)

The rarity and deeply embedded nature of stars with masses larger than 8 solar masses has limited our understanding of their formation. Previous work has shown that complementing spectral energy distributions with interferometric and imaging data can probe the circumstellar environments of massive young stellar objects (MYSOs) well. However, complex studies of single objects often use different approaches in their analysis. Therefore the results of these studies cannot be directly compared. This work aims to obtain the physical characteristics of a sample of MYSOs at ~0.01" scales, at 0.1" scales, and as a whole, which enables us to compare the characteristics of the sources. We apply the same multi-scale method and analysis to a sample of MYSOs. High-resolution interferometric data, near-diffraction-limited imaging data, and a multi-wavelength spectral energy distribution are combined. By fitting simulated observables derived from 2.5D radiative transfer models of disk-outflow-envelope systems to our observations, the properties of the MYSOs are constrained. We find that the observables of all the MYSOs can be reproduced by models with disk-outflow-envelope geometries, analogous to the Class I geometry associated with low-mass protostars. The characteristics of the envelopes and the cavities within them are very similar across our sample. On the other hand, the disks seem to differ between the objects, in particular with regards to what we interpret as evidence of complex structures and inner holes. This is comparable to the morphologies observed for low-mass young stellar objects. A strong correlation is found between the luminosity of the central MYSO and the size of the transition disk-like inner hole for the MYSOs, implying that photoevaporation or the presence of binary companions may be the cause.

We report the results of the Jansky Very Large Array (JVLA) observation of five massive protostars at 6 and 22.2 GHz. The aim of the study was to compare their current fluxes and positions with previous observations to search for evidence of variability. Most of the observed sources present the morphologies of a thermal core, hosting the protostar and exhibiting no proper motion, and associated non-thermal radio lobes that are characterised by proper motions and located away from the thermal core. Some of the protostars drive jets whose lobes have dissimilar displacement vectors, implying precession of the jets or the presence of multiple jet drivers. The jets of the protostars were found to have proper motions that lie in the range 170$ \leq v \leq $650 kms$^{-1}$, and precessions of periods 40$\leq p \leq$50 years and angles 2$\leq \alpha \leq$10$^\circ$, assuming that their velocities $v=$500 kms$^{-1}$. The core of one of the sources, S255 NIRS 3, which was in outburst at the time of our observations, showed a significant change in flux compared to the other sources. Its spectral index decreased during the outburst, consistent with the model of an expanding gas bubble. Modelling the emission of the outburst as that of a new non-thermal lobe that is emerging from a thermal core whose emission enshrouds that of the lobe also has the potential to account for the increase in flux and a decrease in the spectral index of the source's outburst., Comment: 16 pages, 15 figures

The mass segregation of stellar clusters could be primordial rather than dynamical. Despite the abundance of studies of mass segregation for stellar clusters, those for stellar progenitors are still scarce, so the question on the origin and evolution of mass segregation is still open. Our goal is to characterize the structure of the NGC 2264 molecular cloud and compare the populations of clumps and young stellar objects (YSOs) in this region whose rich YSO population has shown evidence of sequential star formation. We separated the Herschel column density map of NGC 2264 in three subregions and compared their cloud power spectra using a multiscale segmentation technique. We identified in the whole NGC 2264 cloud a population of 256 clumps with typical sizes of ~0.1 pc and masses ranging from 0.08 Msun to 53 Msun. Although clumps have been detected all over the cloud, the central subregion of NGC 2264 concentrates most of the massive, bound clumps. The local surface density and the mass segregation ratio indeed indicate a strong degree of mass segregation for the 15 most massive clumps, with a median $\Sigma_6$ three time that of the whole clumps population and $\Lambda_{MSR}$ about 8. We showed that this cluster of massive clumps is forming within a high-density cloud ridge, itself formed and probably still fed by the high concentration of gas observed on larger scales in the central subregion. The time sequence obtained from the combined study of the clump and YSO populations in NGC 2264 suggests that the star formation started in the northern subregion, that it is now actively developing at the center and will soon start in the southern subregion. Taken together, the cloud structure and the clump and YSO populations in NGC 2264 argue for a dynamical scenario of star formation., Comment: 26 pages, 19 figures. Accepted for publication in A&A, abstract abridged for arXiv

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

Massive young stellar objects (MYSOs) have recently been shown to drive jets whose particles can interact with either the magnetic fields of the jet or ambient medium to emit non-thermal radiation. We report a search for non-thermal radio emission from a sample of 15 MYSOs to establish the prevalence of the emission in the objects. We used their spectra across the L-, C- and Q-bands along with spectral index maps to characterise their emission. We find that about 50% of the sources show evidence for non-thermal emission with 40% showing clear non-thermal lobes, especially sources of higher bolometric luminosity. The common or IRAS names of the sources that manifest non-thermal lobes are; V645Cyg, IRAS 22134+5834, NGC 7538 IRS 9, IRAS 23262+640, AFGL 402d and AFGL 490. All the central cores of the sources are thermal with corresponding mass-loss rates that lie in the range 3X10^{-7} to 7X10^{-6} solar masses per year. Given the presence of non-thermal lobes in some of the sources and the evidence of non-thermal emission from some spectral index maps, it seems that magnetic fields play a significant role in the jets of massive protostars. Also noted is that some of the sources show evidence of binarity and variability., Comment: Accepted for publication in MNRAS

The rarity of young massive stars combined with the fact that they are often deeply embedded has limited the understanding of their formation. Ground based mid-infrared (IR) interferometry is one way of securing the spatial resolution required to study massive young stellar objects (MYSOs) and as the spatial-frequency coverage of such observations is often incomplete, direct-imaging can be supplementary to such a dataset. By consolidating these observations with modelling, the features of a massive protostellar environment can be constrained. This work simultaneously fits the aforementioned observations and a spectral energy distribution (SED) with a 2.5D radiative transfer model, providing an extensive view of the physical characteristics of the accreting regions of the MYSO G305.20+0.21. The high-resolution observations were obtained using the Very Large Telescope's MIDI and VISIR instruments, producing visibilities in the N-band and near-diffraction-limited imaging in the Q-band respectively. A model including a central protostar with a luminosity of $\sim$5$\times$10$^{4}$L$_{\odot}$ surrounded by a low-density bipolar cavity, a flared 1M$_{\odot}$ disk and an envelope provides a sufficient fit all three types of observation. The need to include a disk in the model implies that this MYSO follows a scaled-up version of the low-mass star formation process. The weak silicate absorption feature within the SED requires low-density envelope cavities to be successfully fit and is an atypical characteristic in comparison to previously studied MYSOs. Additionally, the inner radius of the disk must be three times the dust sublimation radius to satisfy the MIDI visibilities. The low density, low extinction environment implies the object is a more evolved MYSO and this combined with large inner radius of the disk suggests that it could be an example of a transitional disk around an MYSO., Comment: 16 pages, 8 figures

We present sensitive ($\sigma<10{\rm \mu Jy \, beam^{-1}}$), radio continuum observations using the Australian telescope compact array (ATCA) at frequencies of $6$ and $9$ GHz towards 4 MYSOs. From a previous, less sensitive work, these objects are known to harbour ionized jets associated with radio lobes, which result from shock processes. In comparison with that work, further emission components are detected towards each MYSO. These include extended, direct, thermal emission from the ionized jet's stream, new radio lobes indicative of shocks close ($<10^5\, {\rm au}$) to the MYSO, 3 radio Herbig-Haro objects separated by up to $3.8 {\rm pc}$ from the jet's launching site and an IR-dark source coincident with CH$_3$OH maser emission. No significant, integrated flux variability is detected towards any jets or shocked lobes, and only one proper motion is observed ($1806\pm596\,{\rm km \, s^{-1}}$ parallel to the jet axis of G310.1420+00.7583A). Evidence for precession is detected in all four MYSOs with precession-periods and angles ranging from $66-15480 \, {\rm yr}$ and $6-36 ^\circ$ respectively. Should precession be the result of the influence from a binary companion, we infer orbital radii of $30-1800\, {\rm au} $.

We present a new radio-selected sample of PNe from the CORNISH survey. This is a radio continuum survey of the inner Galactic plane covering Galactic longitude, $10^\circ 1000 K) and located closer than 7 kpc. Within this sample is a water-maser PN with a spectral index of $-0.55\pm 0.08$, which indicates non-thermal radio emission. Such a radio-selected sample, unaffected by extinction, will be particularly useful to compare with population synthesis models and should contribute to the understanding of the formation and evolution of PNe.

A catalogue of 239 ultra-compact HII regions (UCHIIs) found in the CORNISH survey at 5 GHz and 1.5" resolution in the region $10^{\circ} < l < 65^{\circ}, ~|b| < 1^{\circ}$ is presented. This is the largest complete and well-selected sample of UCHIIs to date and provides the opportunity to explore the global and individual properties of this key state in massive star formation at multiple wavelengths. The nature of the candidates was validated, based on observational properties and calculated spectral indices, and the analysis is presented in this work. The physical sizes, luminosities and other physical properties were computed by utilising literature distances or calculating the distances whenever a value was not available. The near- and mid-infrared extended source fluxes were measured and the extinctions towards the UCHIIs were computed. The new results were combined with available data at longer wavelengths and the spectral energy distributions (SEDs) were reconstructed for 177 UCHIIs. The bolometric luminosities obtained from SED fitting are presented. By comparing the radio flux densities to previous observational epochs, we find about 5% of the sources appear to be time variable. This first high-resolution area survey of the Galactic plane shows that the total number of UCHIIs in the Galaxy is ~ 750 - a factor of 3-4 fewer than found in previous large area radio surveys. It will form the basis for future tests of models of massive star formation., Comment: 30 pages; accepted for publication in A&A on 21.03.2018

HD 100546 is one of the few known pre-main-sequence stars that may host a planetary system in its disk. We analyze new VLT/SPHERE/ZIMPOL polarimetric images of HD 100546 with filters in H$\alpha$ and the adjacent continuum. We have probed the disk gap and the surface layers of the outer disk, covering a region < 500 mas (< 55 au at 109 pc) from the star, at an angular resolution of ~ 20 mas. Our data show an asymmetry: the SE and NW regions of the outer disk are more polarized than the SW and NE. This can be explained from a preferential scattering angle close to 90$^o$, consistent with previous polarization images. The outer disk extends from 13 $\pm$ 2 to 45 $\pm$ 9 au, with a position angle and inclination of 137 $\pm$ 5$^o$ and 44 $\pm$ 8$^o$. The comparison with previous estimates suggests that the disk inclination could increase with the stellocentric distance, although the different measurements are still consistent within the error bars. In addition, no direct signature of the innermost candidate companion is detected from polarimetry, confirming recent results based on intensity imagery. We set an upper limit to its mass accretion rate < 10$^{-8}$ M$_{\odot}$/yr for a sub-stellar mass of 15M$_{Jup}$. Finally, we report the first detection (> 3$\sigma$) of a ~ 20 au bar-like structure that crosses the gap through the central region of HD 100546. It is tentatively suggested that the bar could be dust dragged by infalling gas that radially flows from the outer disk. This could represent an exceptional case in which a small-scale radial inflow is observed in a single system. If this scenario is confirmed, it could explain the presence of atomic gas in the inner disk that would otherwise accrete on to the central star on a timescale of a few months/years, as previously indicated from spectro-interferometric data, and could be related with additional (undetected) planets., Comment: 8 pages, 7 figures. Accepted in A&A

We present medium-resolution (R~7000) near-infrared echelle spectroscopic data for 36 MYSOs drawn from the Red MSX Source (RMS) survey. This is the largest sample observed at this resolution at these wavelengths of MYSOs to date. The spectra are characterized mostly by emission from hydrogen recombination lines and accretion diagnostic lines. One MYSO shows photospheric HI absorption, a comparison with spectral standards indicates that the star is an A type star with a low surface gravity, implying that the MYSOs are probably swollen, as also suggested by evolutionary calculations. An investigation of the Br$\gamma$ line profiles shows that most are in pure emission, while 13$\pm$5 % display P Cygni profiles, indicative of outflow, while less than 8$\pm$4 % have inverse P Cygni profiles, indicative of infall. These values are comparable with investigations into the optically bright Herbig Be stars, but not with those of Herbig Ae and T Tauri stars, consistent with the notion that the more massive stars undergo accretion in a different fashion than lower mass objects which are undergoing magnetospheric accretion. Accretion luminosities and rates as derived from the Br$\gamma$ line luminosities agree with results for lower mass sources, providing tentative evidence for massive star formation theories based on scaling of low-mass scenarios. We present Br$\gamma$/Br12 line profile ratios exploiting the fact that optical depth effects can be traced as a function of Doppler shift across the lines. These show that the winds of MYSOs in this sample are nearly equally split between constant, accelerating, and decelerating velocity structures. There are no trends between the types of features we see and bolometric luminosities or near-infrared colours., Comment: 30 pages, 21 figures (with appendices). Accepted for publication in Monthly Notices of the Royal Astronomical Society

Until recently, there have been few studies of the protostellar evolution of intermediate-mass (IM) stars, which may bridge the low- and high-mass regimes. This paper aims to investigate whether the properties of an IM protostar within the IRAS 05373+2349 embedded cluster are similar to that of low- and/or high-mass protostars. We carried out Very Large Array as well as Combined Array for Research in Millimeter Astronomy continuum and 12CO(J=1-0) observations, which uncover seven radio continuum sources (VLA 1-7). The spectral index of VLA 2, associated with the IM protostar is consistent with an ionised stellar wind or jet. The source VLA 3 is coincident with previously observed H2 emission line objects aligned in the north-south direction (P.A. -20 to -12 deg), which may be either an ionised jet emanating from VLA 2 or (shock-)ionised cavity walls in the large-scale outflow from VLA 2. The position angle between VLA 2 and 3 is slightly misaligned with the large-scale outflow we map at ~5-arcsec resolution in 12CO (P.A. ~30 deg), which in the case of a jet suggests precession. The emission from the mm core associated with VLA 2 is also detected; we estimate its mass to be 12-23 Msun, depending on the contribution from ionised gas. Furthermore, the large-scale outflow has properties intermediate between outflows from low- and high-mass young stars. Therefore, we conclude that the IM protostar within IRAS 05373+2349 is phenomenologically as well as quantitatively intermediate between the low- and high-mass domains., Comment: 11 pages, 7 figures and 4 tables. Accepted by MNRAS

We present the results of the application of locally linear embedding (LLE) to reduce the dimensionality of dereddened and continuum subtracted near-infrared spectra using a combination of models and real spectra of massive protostars selected from the Red MSX Source survey database. A brief comparison is also made with two other dimension reduction techniques; Principal Component Analysis (PCA) and Isomap using the same set of spectra as well as a more advanced form of LLE, Hessian locally linear embedding. We find that whilst LLE certainly has its limitations, it significantly outperforms both PCA and Isomap in classification of spectra based on the presence/absence of emission lines and provides a valuable tool for classification and analysis of large spectral data sets., Comment: 8 pages, 7 figures. Accepted for publication in MNRAS 2016 June 21

Radio continuum observations using the Australia telescope compact array at 5.5, 9.0, 17.0 and 22.8 GHz have detected free-free emission associated with 45 of 49 massive young stellar objects and HII regions. Of these, 26 sources are classified as ionized jets (12 of which are candidates), 2 as ambiguous jets or disc winds, 1 as a disc-wind, 14 as HII regions and 2 were unable to be categorised. Classification as ionized jets is based upon morphology, radio flux and spectral index, in conjunction with previous observational results at other wavelengths. Radio-luminosity and momentum are found to scale with bolometric luminosity in the same way as low-mass jets, indicating a common mechanism for jet production across all masses. In 13 of the jets, we see associated non-thermal/optically-thin lobes resulting from shocks either internal to the jet and/or at working surfaces. Ten jets display non-thermal (synchrotron emission) spectra in their lobes, with an average spectral index of -0.55 consistent with Fermi acceleration in shocks. This shows that magnetic fields are present, in agreement with models of jet formation incorporating magnetic fields. Since the production of collimated radio jets is associated with accretion processes, the results presented in this paper support the picture of disc-mediated accretion for the formation of massive stars with an upper-limit on the jet phase lasting approximately $6.5 \times 10^4 yr$. Typical mass loss rates in the jet are found to be $1.4 \times 10^{-5} M_\odot yr^{-1}$ with associated momentum rates of the order $(1-2) \times 10^{-2} M_\odot km s^{-1} yr^{-1}$.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average., Comment: Accepted by MNRAS

We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molecular outflows associated with massive star formation. Of these, 59, 17 and 13 sources (66, 19 and 15 percent) are found to have outflows, show some evidence of outflow, and have no evidence of outflow, respectively. The time-dependent parameters of the high-velocity molecular flows are calculated using a spatially variable dynamic timescale. The canonical correlations between the outflow parameters and source luminosity are recovered and shown to scale with those of low-mass sources. For coeval star formation we find the scaling is consistent with all the protostars in an embedded cluster providing the outflow force, with massive stars up to $\sim$30 M$_{\odot}$ generating outflows. Taken at face value, the results support the model of a scaled-up version of the accretion-related outflow-generation mechanism associated with discs and jets in low-mass objects with time-averaged accretion rates of $\sim$10$^{-3}$ M$_{\odot}$ yr$^{-1}$ onto the cores. However, we also suggest an alternative model, in which the molecular outflow dynamics are dominated by the entrained mass and are unrelated to the details of the acceleration mechanism. We find no evidence that outflows contribute significantly to the turbulent kinetic energy of the surrounding dense cores.

We analyse C$^{18}$O ($J=3-$2) data from a sample of 99 infrared-bright massive young stellar objects (MYSOs) and compact HII regions that were identified as potential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a distance limited (D $<$ 6 kpc) sample shown to be representative of star formation covering the transition between the source types. At the spatial resolution probed, Larson-like relationships are found for these cores, though the alternative explanation, that Larson's relations arise where surface-density-limited samples are considered, is also consistent with our data. There are no significant differences found between source properties for the MYSOs and HII regions, suggesting that the core properties are established prior to the formation of massive stars, which subsequently have little impact at the later evolutionary stages investigated. There is a strong correlation between dust-continuum and C$^{18}$O-gas masses, supporting the interpretation that both trace the same material in these IR-bright sources. A clear linear relationship is seen between the independently established core masses and luminosities. The position of MYSOs and compact HII regions in the mass-luminosity plane is consistent with the luminosity expected a cluster of protostars when using a $\sim$40 percent star-formation efficiency and indicates that they are at a similar evolutionary stage, near the end of the accretion phase.

Correlations between the accretion luminosity and emission line luminosities (L_acc and L_line) of pre-main sequence (PMS) stars have been published for many different spectral lines, which are used to estimate accretion rates. Despite the origin of those correlations is unknown, this could be attributed to direct or indirect physical relations between the emission line formation and the accretion mechanism. This work shows that all (near-UV/optical/near-IR) L_acc-L_line correlations are the result of the fact that the accretion luminosity and the stellar luminosity (L_star) are correlated, and are not necessarily related with the physical origin of the line. Synthetic and observational data are used to illustrate how the L_acc-L_line correlations depend on the L_acc-L_star relationship. We conclude that because PMS stars show the L_acc-L_star correlation immediately implies that L_acc also correlates with the luminosity of all emission lines, for which the L_acc-L_line correlations alone do not prove any physical connection with accretion but can only be used with practical purposes to roughly estimate accretion rates. When looking for correlations with possible physical meaning, we suggest that L_acc/L_star and L_line/L_star should be used instead of L_acc and L_line. Finally, the finding that L_acc has a steeper dependence on L_star for T-Tauri stars than for intermediate-mass Herbig Ae/Be stars is also discussed. That is explained from the magnetospheric accretion scenario and the different photospheric properties in the near-UV., Comment: 9 pages, 5 figures, 2 tables. Accepted in MNRAS

We present the results of ammonia observations towards 66 massive star forming regions identified by the Red MSX source survey. We have used the Green Bank Telescope and the K-band focal plane array to map the ammonia NH3 (1,1) and (2,2) inversion emission at a resolution of 30 arcsec in 8 arcmin regions towards the positions of embedded massive star formation. We have identified a total of 115 distinct clumps, approximately two-thirds of which are associated with an embedded massive young stellar object or compact HII region, while the others are classified as quiescent. There is a strong spatial correlation between the peak NH3 emission and the presence of embedded objects. We derive the spatial distribution of the kinetic gas temperatures, line widths, and NH$_3$ column densities from these maps, and by combining these data with dust emission maps we estimate clump masses, H$_2$ column densities and ammonia abundances. The clumps have typical masses of ~1000 Msun and radii ~0.5 pc, line widths of ~2 km/s and kinetic temperatures of ~16-20 K. We find no significant difference between the sizes and masses of the star forming and quiescent subsamples; however, the distribution maps reveal the presence of temperature and line width gradients peaking towards the centre for the star forming clumps while the quiescent clumps show relatively uniform temperatures and line widths throughout. Virial analysis suggests that the vast majority of clumps are gravitationally bound and are likely to be in a state of global free fall in the absence of strong magnetic fields. The similarities between the properties of the two subsamples suggest that the quiescent clumps are also likely to form massive stars in the future, and therefore provide a excellent opportunity to study the initial conditions of massive pre-stellar and protostellar clumps., Comment: Accepted for publication by MNRAS. 25 pages, 21 figures and 6 tables. This paper also contains 100 pages of additional online images and tables that are not included here. A full resolution version of the paper including all of the online material can be found here http://rms.leeds.ac.uk/cgi-bin/public/RMS_PUBLICATIONS.cgi

By matching infrared-selected, massive young stellar objects (MYSOs) and compact HII regions in the RMS survey to massive clumps found in the submillimetre ATLASGAL survey, we have identified ~1000 embedded young massive stars between 280\degr < $\ell$ < 350\degr and 10degr < $\ell$ < 60\degr with |b|<1.5degr. Combined with an existing sample of radio-selected methanol masers and compact HII regions, the result is a catalogue of ~1700 massive stars embedded within ~1300 clumps located across the inner Galaxy, containing three observationally distinct subsamples, methanol-maser, MYSO and HII-region associations, covering the most important tracers of massive star formation, thought to represent key stages of evolution. We find that massive star formation is strongly correlated with the regions of highest column density in spherical, centrally condensed clumps. We find no significant differences between the three samples in clump structure or the relative location of the embedded stars, which suggests that the structure of a clump is set before the onset of star formation, and changes little as the embedded object evolves towards the main sequence. There is a strong linear correlation between clump mass and bolometric luminosity, with the most massive stars forming in the most massive clumps. We find that the MYSO and HII-region subsamples are likely to cover a similar range of evolutionary stages and that the majority are near the end of their main accretion phase. We find few infrared-bright MYSOs associated with the most massive clumps, probably due to very short pre-main sequence lifetimes in the most luminous sources., Comment: Accepted for publication in MNRAS: 33 pages, 24 figures and 9 tables. Full versions of Fig. 3 and Tables 1, 3, 6 and 7 will only be available in electronic form via CDS. The resolution of Figs. 1, 5, 7 and 21 have been reduced but a higher resolution version of the paper is available upon request

Abridged: We have used the well-selected sample of ~1750 embedded, young, massive stars identified by the RMS survey to investigate the Galactic distribution of recent massive star formation. We describe the various methods used to assign distances extracted from the literature, and solve the distance ambiguities towards ~200 sources located within the Solar circle using archival HI data. These distances are used to calculate bolometric luminosities and estimate the survey completeness (~2x10^4 lsun). In total, we calculate the distance and luminosity of ~1650 sources, one third of which are above the survey's completeness threshold. Examination of the sample's longitude, latitude, radial velocities and mid-infrared images has identified ~120 small groups of sources, many of which are associated with well known star formation complexes, such as W43, W49 and W51. We compare the positional distribution of the sample with the expected locations of the spiral arms, assuming a model of the Galaxy consisting of four gaseous arms. The distribution of young massive stars in the Milky Way is spatially correlated with the spiral arms, with strong peaks in the source position and luminosity distributions at the arms' Galactocentric radii. The overall source and luminosity surface densities are both well correlated with the surface density of the molecular gas, which suggests that the massive star formation rate per unit molecular mass is approximately constant across the Galaxy. We estimate the total luminosity of the embedded massive star population to be ~0.76x10^8 lsun, 30% of which is associated with the ten most active star forming complexes. We measure the scale height as a function of Galactocentric distance and find that it increases only modestly from ~20-30 pc between 4 and 8 kpc, but much more rapidly at larger distances., Comment: Accepted for publication in MNRAS. Paper consists of 17 pages, four table and 11 figures. The full versions of Tables 1,2 and 3 and Figures 1 and 4 will only be available from CDS

We present medium resolution ($R\approx5300$) $K$-band integral field spectroscopy of six MYSOs. The targets are selected from the RMS survey, and we used the NIFS on the Gemini North telescope. The data show various spectral line features including Br$\gamma$, CO, H$_2$, and \mbox{He\,{\sc i}}. The Br$\gamma$ line is detected in emission in all objects with $v_\mathrm{FWHM}\sim100$ -- 200 kms$^{-1}$. V645 Cyg shows a high-velocity P-Cygni profile between -800 kms$^{-1}$ and -300 kms$^{-1}$. We performed three-dimensional spectroastrometry to diagnose the circumstellar environment in the vicinity of the central stars using the Br$\gamma$ line. We measured the centroids of the velocity components with sub-mas precision. The centroids allow us to discriminate the blueshifted and redshifted components in a roughly east--west direction in both IRAS 18151--1208 and S106 in Br$\gamma$. This lies almost perpendicular to observed larger scale outflows. We conclude, given the widths of the lines and the orientation of the spectroastrometric signature, that our results trace a disc wind in both IRAS 18151--1208 and S106. The CO $\nu=2-0$ absorption lines at low $J$ transitions are detected in IRAS 18151--1208 and AFGL 2136. We analysed the velocity structure of the neutral gas discs. In IRAS 18151--1208, the absorption centroids of the blueshifted and redshifted components are separated in a direction of north-east to south-west, nearly perpendicular to that of the larger scale $H_2$ jet. The position-velocity relations of these objects can be reproduced with central masses of 30 $M_{\sun}$ for IRAS 18151--1208 and 20 $M_{\sun}$ for AFGL 2136. We also detect CO $\nu=2-0$ bandhead emission in IRAS 18151--1208, S106 and V645 Cyg. The results can be fitted reasonably with a Keplerian rotation model, with masses of 15, 20 and 20 $M_{\sun}$ respectively., Comment: 17 pages, 10 figures, accepted by MNRAS

We present the Red MSX Source (RMS) Survey, the largest statistically selected catalog of young massive protostars and HII regions to date. We outline the construction of the catalog using mid and near infrared color selection, as well as the detailed follow up work at other wavelengths, and at higher spatial resolution in the infrared. We show that within the adopted selection bounds we are more than 90% complete for the massive protostellar population, with a positional accuracy of the exciting source of better than 2 arcseconds. We briefly summarize some of the results that can be obtained from studying the properties of the objects in the catalog as a whole, and find evidence that the most massive stars form: (i) preferentially nearer the Galactic centre than the anti-centre; (ii) in the most heavily reddened environments, suggestive of high accretion rates; and (iii) from the most massive cloud cores., Comment: Accepted for publication in ApJS. The full version of Table 1 will be made available electronically, but the most up-to-date version of our database will always be that accessible at http://rms.leeds.ac.uk

We present a complete sample of molecular clumps containing compact and ultra-compact (UC) HII regions between \ell=10\degr and 60\degr\ and $|b|<1\degr, identified by combining the the ATLASGAL submm and CORNISH radio continuum surveys with visual examination of archival infrared data. Our sample is complete to optically thin, compact and UCHII regions driven by a zero age main sequence star of spectral type B0 or earlier embedded within a 1,000 Msun clump. In total we identify 213 compact and UCHII regions, associated with 170 clumps. Unambiguous kinematic distances are derived for these clumps and used to estimate their masses and physical sizes, as well as the Lyman continuum fluxes and sizes of their embedded HII regions. We find a clear lower envelope for the surface density of molecular clumps hosting massive star formation of 0.05 g cm^{-2}, which is consistent with a similar sample of clumps associated with 6.7 GHz masers. The mass of the most massive embedded stars is closely correlated with the mass of their natal clump. Young B stars appear to be significantly more luminous in the ultraviolet than predicted by current stellar atmosphere models. The properties of clumps associated with compact and UCHII regions are very similar to those associated with 6.7 GHz methanol masers and we speculate that there is little evolution in the structure of the molecular clumps between these two phases. Finally, we identify a significant peak in the surface density of compact and UCHII regions associated with the W49A star-forming complex, noting that this complex is truly one of the most massive and intense regions of star formation in the Galaxy., Comment: Accepted for publication in MNRAS. Paper consists of 29 pages, 26 figures and 4 tables. The resolution of Figures 21 and 25 has been reduced to comply with file size restrictions. Full versions of Figures 11 and 12 and Tables 3 and 4 will only be available in the electronic version of the journal and through CDS

Near-infrared H- and K-band spectra are presented for 247 objects, selected from the Red MSX Source (RMS) survey as potential young stellar objects (YSOs). 195 (~80%) of the targets are YSOs, of which 131 are massive YSOs (L_BOL > 5x10^3 L_solar), M > 8M_solar. This is the largest spectroscopic study of massive YSOs to date, providing a valuable resource for the study of massive star formation. In this paper we present our exploratory analysis of the data. The YSOs observed have a wide range of embeddedness (2.7 < A_V < 114), demonstrating that this study covers minimally obscured objects right through to very red, dusty sources. Almost all YSOs show some evidence for emission lines, though there is a wide variety of observed properties. The most commonly detected lines are Brgamma, H_2, fluorescent FeII, CO bandhead, [FeII] and HeI 2-1 2^1S-2^1P, in order of frequency of occurrence. In total, ~40% of the YSOs display either fluorescent FeII 1.6878um or CO bandhead emission (or both), indicative of a circumstellar disc; however, no correlation of the strength of these lines with bolometric luminosity was found. We also find that ~60% of the sources exhibit [FeII] or H_2 emission, indicating the presence of an outflow. Three quarters of all sources have Brgamma in emission. A good correlation with bolometric luminosity was observed for both the Brgamma and H_2 emission line strengths, covering 1 L_solar< L_BOL < 3.5x10^5 L_solar. This suggests that the emission mechanism for these lines is the same for low-, intermediate-, and high-mass YSOs, i.e. high-mass YSOs appear to resemble scaled-up versions of low-mass YSOs., Comment: Accepted for publication in MNRAS

Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 deg < l < 65 deg), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.5" resolution Stokes I map with a root-mean-squared noise level better than 0.4 mJy/beam. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalogue of compact radio-emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimised to detect emission on size scales up to 14" and for unresolved sources the catalogue is more than 90 percent complete at a flux density of 3.9 mJy. We have detected 3,062 sources above a 7-sigma detection limit and present their ensemble properties. The catalogue is highly reliable away from regions containing poorly-sampled extended emission, which comprise less than two percent of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artefacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: regions, planetary nebulae and radio-galaxies. The CORNISH data and catalogue are available online at http://cornish.leeds.ac.uk, Comment: 25 pages, 22 figures. To appear in ApJ. Sup

We describe the motivation, design and implementation of the CORNISH survey, an arcsecond resolution radio continuum survey of the inner Galactic plane at 5 GHz using the Karl G. Jansky Very Large Array (VLA). It is a blind survey co-ordinated with the northern Spitzer GLIMPSE I region covering 10 deg < l < 65 deg and |b| < 1 deg at similar resolution. We discuss in detail the strategy that we employed to control the shape of the synthesised beam across this survey that covers a wide range of fairly low declinations. Two snapshots separated by 4 hours in hour angle kept the beam elongation to less that 1.5 over 75% of the survey area and less than 2 over 98% of the survey. The prime scientific motivation is to provide an unbiased survey for ultra-compact HII regions to study this key phase in massive star formation. A sensitivity around 2 mJy will allow the automatic distinction between radio loud and quiet mid-IR sources found in the Spitzer surveys. This survey has many legacy applications beyond star formation including evolved stars, active stars and binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other Galactic plane surveys that target diffuse and non-thermal sources as well as atomic and molecular phases to build up a complete picture of the ISM in the Galaxy., Comment: 17 pages, 14 figures. Accepted for publication in PASP (16th-August-2012)

We present a study of the kinematical properties of a small sample of nearby near-infrared bright massive and intermediate mass young stellar objects using emission lines sensitive to discs and winds. We show for the first time that the broad ($\sim500$kms$^{-1}$) symmetric line wings on the HI Brackett series lines are due to Stark broadening or electron scattering, rather than pure Doppler broadening due to high speed motion. The results are consistent with the presence of a very dense circumstellar environment. In addition, many of these lines show evidence for weak line self-absorption, suggestive of a wind or disc-wind origin for that part of the absorbing material. The weakness of the self-absorption suggests a large opening angle for such an outflow. We also study the fluorescent 1.688$\mu$m FeII line, which is sensitive to dense material. We fitted a Keplerian disc model to this line, and find reasonable fits in all bar one case, in agreement with previous finding for classical Be stars that fluorescent iron transitions are reasonable disc tracers. Overall the picture is one in which these stars still have accretion discs, with a very dense inner circumstellar environment which may be tracing either the inner regions of a disc, or of a stellar wind, and in which ionised outflow is also present. The similarity with lower mass stars is striking, suggesting that at least in this mass range they form in a similar fashion., Comment: 20 pages, 7 figures, accepted for publication in MNRAS

Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accretion in massive star formation are not well understood. This study evaluates the source of the Q-band, lambda=19.5 microns, emission of massive young stellar objects (MYSOs). This allows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 1000 AU scales. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6", corresponds to approximately 1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis-symmetric dust radiative transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact HII regions and are thus more evolved. We find that, in general, the diverse 20 micron morphology of MYSOs can be attributed to warm dust in the walls of outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation., Comment: Accepted for publication in A&A. The images in this version have been compressed. A high resolution version is available on request

We present high-resolution HI data obtained using the Australia Telescope Compact Array to resolve the near/far distance ambiguities towards a sample of compact HII regions from the Red MSX Source (RMS) survey. The high resolution data are complemented with lower resolution archival HI data extracted from the Southern and VLA Galactic Plane surveys. We resolve the distance ambiguity for nearly all of the 105 sources where the continuum was strong enough to allow analysis of the HI absorption line structure. This represents another step in the determination of distances to the total RMS sample, which with over 1,000 massive young stellar objects and compact HII regions, is the largest and most complete sample of its kind. The full sample will allow the distribution of massive star formation in the Galaxy to be examined., Comment: Accepted by MNRAS. This paper consists of 15 pages and contains 10 figures and 5 tables

The Red MSX Source (RMS) survey has identified a sample of ~1200 massive young stellar objects (MYSOs), compact and ultra compact HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. We have used the 100 m Green Bank telescope to search for 22-24 GHz water maser and ammonia (1,1), (2,2) and (3,3) emission towards ~600 RMS sources located within the northern Galactic plane. We have identified 308 H2O masers which corresponds to an overall detection rate of ~50%. Abridged: We detect ammonia emission towards 479 of these massive young stars, which corresponds to ~80%. Ammonia is an excellent probe of high density gas allowing us to measure key parameters such as gas temperatures, opacities, and column densities, as well as providing an insight into the gas kinematics. The average kinetic temperature, FWHM line width and total NH3 column density for the sample are approximately 22 K, 2 km/s and 2x10^{15} cm^{-2}, respectively. We find that the NH3 (1,1) line width and kinetic temperature are correlated with luminosity and finding no underlying dependence of these parameters on the evolutionary phase of the embedded sources, we conclude that the observed trends in the derived parameters are more likely to be due to the energy output of the central source and/or the line width-clump mass relationship. The velocities of the peak H2O masers and the NH3 emission are in excellent agreement with each other, which would strongly suggest an association between the dense gas and the maser emission. Moreover, we find the bolometric luminosity of the embedded source and the isotropic luminosity of the H2O maser are also correlated. We conclude from the correlations of the cloud and water maser velocities and the bolometric and maser luminosity that there is a strong dynamical relationship between the embedded young massive star and the H2O maser., Comment: 17 pages and 17 figures and 8 tables. Tables\,2 and 5 and full versions of Figs. 3 and 7 are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/

We apply mid-infrared spectro-interferometry to the massive young stellar object CRL2136. The observations were performed with the Very Large Telescope Interferometer and the MIDI instrument at a 42m baseline probing angular scales of 50 milli-arcseconds. We model the observed visibilities in parallel with diffraction-limited images at both 24.5micron and in the N-band (with resolutions of 0.6" and 0.3", respectively), as well as the spectral energy distribution. The arcsec-scale spatial information reveals the well-resolved emission from the dusty envelope. By simultaneously modelling the spatial and spectral data, we find that the bulk of the dust emission occurs at several dust sublimation radii (approximately 170 AU). This reproduces the high mid-infrared fluxes and at the same time the low visibilities observed in the MIDI data for wavelengths longward of 8.5micron. However, shortward of this wavelength the visibility data show a sharp up-turn indicative of compact emission. We discuss various potential sources of this emission. We exclude a dust disk being responsible for the observed spectral imprint on the visibilities. A cool supergiant star and an accretion disk are considered and both shown to be viable origins of the compact mid-infrared emission. We propose that CRL2136 is embedded in a dusty envelope, which truncates at several times the dust sublimation radius. A dust torus is manifest in the equatorial region. We find that the spectro-interferometric N-band signal can be reproduced by either a gaseous disk or a bloated central star. If the disk extends to the stellar surface, it accretes at a rate of 3.0 10^(-3) Msun/yr., Comment: Accepted for publication in A&A

Context: The Red MSX Source (RMS) survey is returning a large sample of massive young stellar objects (MYSOs) and ultra-compact (UC) \HII{} regions using follow-up observations of colour-selected candidates from the MSX point source catalogue. Aims: To obtain the bolometric fluxes and, using kinematic distance information, the luminosities for young RMS sources with far-infrared fluxes. Methods: We use a model spectral energy distribution (SED) fitter to obtain the bolometric flux for our sources, given flux data from our work and the literature. The inputs to the model fitter were optimised by a series of investigations designed to reveal the effect varying these inputs had on the resulting bolometric flux. Kinematic distances derived from molecular line observations were then used to calculate the luminosity of each source. Results: Bolometric fluxes are obtained for 1173 young RMS sources, of which 1069 have uniquely constrained kinematic distances and good SED fits. A comparison of the bolometric fluxes obtained using SED fitting with trapezium rule integration and two component greybody fits was also undertaken, and showed that both produce considerable scatter compared to the method used here. Conclusions: The bolometric flux results allowed us to obtain the luminosity distributions of YSOs and UC\HII{} regions in the RMS sample, which we find to be different. We also find that there are few MYSOs with L $\geq$ 10$^{5}$\lsol{}, despite finding many MYSOs with 10$^{4}$\lsol{} $\geq$ L $\geq$ 10$^{5}$\lsol{}., Comment: 12 pages, 12 figures, 3 tables, accepted to A&A. The full versions of tables 1 and 2 will be available via the CDS upon publication

The Red MSX Source (RMS) survey has identified a large sample of massive young stellar objects (MYSOs) and ultra compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. Using a recent catalogue of molecular clouds derived from the Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey (GRS), and by applying a Galactic scaleheight cut off of 120 pc, we solve the distance ambiguity for RMS sources located within 18\degr < |l| > 54\degr. These two steps yield kinematic distances to 291 sources out of a possible 326 located within the GRS longitude range. Combining distances and integrated fluxes derived from spectral energy distributions, we estimate luminosities to these sources and find that > 90% are indicative of the presence of a massive star. We find the completeness limit of our sample is ~10^4 Lsun, which corresponds to a zero age main sequence (ZAMS) star with a mass of ~12 Msun. Selecting only these sources, we construct a complete sample of 196 sources. Comparing the properties of the sample of young massive stars with the general population, we find the RMS-clouds are generally larger, more massive, and more turbulent. We examine the distribution of this sub-sample with respect to the location of the spiral arms and the Galactic bar and find them to be spatially correlated. We identify three significant peaks in the source surface density at Galactocentric radii of approximately 4, 6 and 8 kpc, which correspond to the proposed positions of the Scutum, Sagittarius and Perseus spiral arms, respectively. Fitting a scale height to the data we obtain an average value of ~29+-0.5 pc, which agrees well with other reported values in the literature, however, we note a dependence of the scale height on galactocentric radius with it increases from 30 pc to 45 pc between 2.5 and 8.5 kpc., Comment: Accepted for publication by MNRAS. Paper consists of 15 pages including 12 figures and four tables. Full versions of Tables 2 and 3 will only be available online. The resolution of Figure 9 has been reduced - a full resolution version of the paper can be download from here: http://www.ast.leeds.ac.uk/cgi-bin/RMS/RMS_PUBLICATIONS.cgi

We present high resolution (R~50,000) spectroastrometry over the CO 1st overtone bandhead of a sample of seven intermediate/massive young stellar objects. These are primarily drawn from the red MSX source (RMS) survey, a systematic search for young massive stars which has returned a large, well selected sample of such objects. The mean luminosity of the sample is approximately 5 times 10^4 L_\odot, indicating the objects typically have a mass of ~15 solar masses. We fit the observed bandhead profiles with a model of a circumstellar disc, and find good agreement between the models and observations for all but one object. We compare the high angular precision (0.2-0.8 mas) spectroastrometric data to the spatial distribution of the emitting material in the best-fitting models. No spatial signatures of discs are detected, which is entirely consistent with the properties of the best-fitting models. Therefore, the observations suggest that the CO bandhead emission of massive young stellar objects originates in small-scale disks, in agreement with previous work. This provides further evidence that massive stars form via disc accretion, as suggested by recent simulations., Comment: Accepted for publication in MNRAS

The circumstellar structure on 100 AU scales of the massive young stellar object W33A is probed using the VLTI and the MIDI instrument. N-band visibilities on 4 baselines are presented which are inconsistent with a spherically symmetric geometry. The visibility spectra and SED are simultaneously compared to 2D axi-symmetric dust radiative transfer models with a geometry including a rotationally flattened envelope and outflow cavities. We assume an O7.5 ZAMS star as the central source, consistent with the observed bolometric luminosity. The observations are also compared to models with and without (dusty and gaseous) accretion disks. A satisfactory model is constructed which reproduces the visibility spectra for each (u,v) point. It fits the silicate absorption, the mid-IR slope, the far-infrared peak, and the (sub)mm of the SED. It produces a 350 micron morphology consistent with observations. The 10 micron emission on 100 AU scales is dominated by the irradiated walls of the cavity sculpted by the outflow. The visibilities rule out the presence of dust disks with total (gas and dust) masses more than 0.01 Msun. However, optically thick accretion disks, interior to the dust sublimation radius, are allowed to accrete at rates equalling the envelope's mass infall rate (up to 10^(-3) Msun/yr) without substantially affecting the visibilities due to the extinction by the extremely massive envelope of W33A., Comment: Accepted for publication in A&A

Context: The Red MSX Source (RMS) survey is a multi-wavelength campaign of follow-up observations of a colour-selected sample of candidate massive young stellar objects (MYSOs) in the galactic plane. This survey is returning the largest well-selected sample of MYSOs to date, while identifying other dust contaminant sources with similar mid-infrared colours including a large number of new ultra-compact (UC)HII regions. Aims:To measure the far-infrared (IR) flux, which lies near the peak of the spectral energy distribution (SED) of MYSOs and UCHII regions, so that, together with distance information, the luminosity of these sources can be obtained. Methods:Less than 50% of RMS sources are associated with IRAS point sources with detections at 60 micron and 100 micron, though the vast majority are visible in Spitzer MIPSGAL or IRAS Galaxy Atlas (IGA) images. However, standard aperture photometry is not appropriate for these data due to crowding of sources and strong spatially variable far-IR background emission in the galactic plane. A new technique using a 2-dimensional fit to the background in an annulus around each source is therefore used to obtain far-IR photometry for young RMS sources. Results:Far-IR fluxes are obtained for a total of 1113 RMS candidates identified as young sources. Of these 734 have flux measurements using IGA 60 micron and 100 micron images and 724 using MIPSGAL 70 micron images, with 345 having measurements in both data sets., Comment: 10 pages, 10 figures, 2 Tables, accepted to A&A. A full version of table 1 is available from the lead author or at the CDS upon publication

Context: The Red MSX Source (RMS) survey has identified a large sample of candidate massive young stellar objects (MYSOs) and ultra compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. Aims: To search for H2O masers towards a large sample of young high mass stars and to investigate the statistical correlation of H2O masers with the earliest stages of massive star formation. Methods: We have used the Mopra Radio telescope to make position-switched observations towards ~500 UCHII regions and MYSOs candidates identified from the RMS survey and located between 190\degr < l < 30\degr. These observations have a 4$\sigma$ sensitivity of ~1 Jy and a velocity resolution of ~0.4 km/s.} Results: We have detected 163 H2O masers, approximately 75% of which were previously unknown. Comparing the maser velocities with the velocities of the RMS sources, determined from 13CO observations, we have identified 135 RMS-H2O maser associations, which corresponds to a detection rate of ~27%. Taking into account the differences in sensitivity and source selection we find our detection rate is in general agreement with previously reported surveys. Conclusions: We find similar detection rates for UCHII regions and MYSOs candidates, suggesting that the conditions needed for maser activity are equally likely in these two stages of the star formation process. Looking at the detection rate as a function of distance from the Galactic centre we find it significantly enhanced within the solar circle, peaking at ~37% between 6-7 kpc, which is consistent with previous surveys of UC HII regions, possibly indicating the presence of a high proportion of more luminous YSOs and HII regions., Comment: Accepted for publication in Astronomy and Astrophysics. 10 pages, 7 figures. Full versions of Figure 2 and Tables 1 and 2 will only be available at the CDS or from the lead author

(Abridged) Context: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, well-selected sample of massive young stellar objects (MYSOs). We have identified $\sim$2000 MYSO candidates located throughout the Galaxy by comparing the colours of MSX and 2MASS point sources to those of known MYSOs. Aims: To identify the populations of UCHII regions and PNe within the sample and examine their Galactic distribution. Method: We have conducted high resolution radio continuum observations at 6 cm towards 659 MYSO candidates in the northern hemisphere ($10\degr< l < 250\degr$) using the VLA. In addition to these targeted observations we present archival data towards a further 315 RMS sources extracted from a previous VLA survey of the inner Galaxy. Results: We find radio emission towards 272 ($\sim$27% of the observed sample). Using results from other parts of our multi-wavelength survey we separate these RMS-radio associations into two distinct types of objects, classifying 51 as PNe and a further 208 as either compact or UC HII regions. Using this well selected sample of HII regions we estimate their Galactic scale height to be 0.6\degr. Conclusions: Using radio continuum and archival data we have identified 79 PNe and 391 HII regions within the northern RMS catalogue. We estimate the total fraction of contamination by PNe in the RMS sample is of order 10%. The sample of HII regions is probably the best representation to date of the Galactic population of HII regions as a whole., Comment: Accepted for publication in Astronomy and Astrophysics. 15 pages, 9 figures and 5 tables. Full versions of Tables 3, 4 and 5 and Figs. 2, 4 and 7 will only be available via CDS or the RMS website at http:/www.ast.leeds.ac.uk/cgi-bin/RMS/RMS_VLA_IMAGES.cgi

Context: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, high-resolution mid-infrared colour-selected sample of massive young stellar objects. Aims: A critical part of our follow-up programme is to conduct 13CO molecular line observations in order to determine kinematic distances to all of our MYSO candidates. These distances will allow us to identify and remove nearby low-mass YSOs and help in identifying evolved stars which are weak CO emitters. Method: We have used the 15 m James Clerk Maxwell Telescope (JCMT), the 13.7 m telescope of the Purple Mountain Observatory (PMO), the 20 m Onsala telescope and the 22m Mopra telescope to conduct molecular line observations towards 508 MYSOs candidates located in the 1st and 2nd Quadrants. Results: We detect 13CO emission towards 780 RMS sources which corresponds to approximately 84% of those observed (911). A total of 2595 emission components are detected above 3sigma level (typically T^*_{\rm{A}} > 0.3K), with multiple components being observed towards the majority of these sources -- 520 sources (~56%) -- with an average of ~4 molecular clouds detected along each line of sight. We have used archival CS (J=2-1) and maser velocities to resolve the component multiplicity towards 175 sources (~20%) and have derived a criterion which is used to identify the most likely component for a further 191 multiple component sources. Combined with the single component detections we have obtained unambiguous kinematic velocities for 638 of the 780 MYSOs candidates towards which CO is detected (~80% of the detections). Using the Galactic rotation curve we calculate kinematic distances for all detected components., Comment: 13 pages and 9 figures, accepted for publication in A&A. Includes a large table and figure to be made available via CDS

Here we describe the Red MSX Source (RMS) survey which is the largest, systematic, galaxy-wide search for massive young stellar objects (MYSOs) yet undertaken. Mid-IR bright point sources from the MSX satellite survey have been followed-up with ground-based radio, millimetre, and infrared observations to identify the contaminating sources and characterise the MYSOs and UCHII regions. With the initial classification now complete the distribution of sources in the galaxy will be discussed, as well as some programmes being developed to exploit our sample., Comment: 8 pages, 3 figures, to appear in the proceedings for "Massive Star Formation: Observations confront Theory 2007" edited by H. Beuther