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 detaile
d 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.
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 ($sim50
0$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.
Using near infrared UKIDSS Galactic Plane Survey data, we make extinction measurements to individual stars along the same line of sight as molecular clouds. Using an existing 3D extinction map of the inner Galaxy, that provides line of sight specific
extinction-distance relationships, we convert the measured extinction of molecular clouds to a corresponding distance. These distances are derived independently from kinematic methods, typically used to derive distances to molecular clouds, and as such they have no near/far ambiguity. The near/far distance ambiguity has been resolved for 27 clouds and distances have been derived to 20 clouds. The results are found to be in good agreement with kinematic measurements to molecular clouds where the ambiguity has already been resolved, using HI self-absorption techniques.
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 envelopes 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.
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.
We determine the slope of the near infrared extinction power law (A$_{lambda} propto lambda^{-alpha}$) for 8 regions of the Galaxy between l$sim27^{circ}$ and $sim100^{circ}$. UKIDSS Galactic Plane Survey data are compared, in colour-colour space, wi
th Galactic population synthesis model data reddened using a series of power laws and convolved through the UKIDSS filter profiles. Monte Carlo simulations allow us to determine the best fit value of $alpha$ and evaluate the uncertainty. All values are consistent with each other giving an average extinction power law of $alpha$=2.14$^{+0.04}_{-0.05}$. This is much steeper than most laws previously derived in the literature from colour excess ratios, which are typically between 1.6 and 1.8. We show that this discrepancy is due to an inappropriate choice of filter wavelength in conversion from colour excess ratios to $alpha$ and that effective rather than isophotal wavelengths are more appropriate. In addition, curved reddening tracks, which depend on spectral type and filter system, should be used instead of straight vectors.
Massive young stellar objects (MYSO) are surrounded by massive dusty envelopes. Our aim is to establish their density structure on scales of ~1000 AU, i.e. a factor 10 increase in angular resolution compared to similar studies performed in the (sub)m
m. We have obtained diffraction-limited (0.6) 24.5 micron images of 14 well-known massive star formation regions with Subaru/COMICS. The images reveal the presence of discrete MYSO sources which are resolved on arcsecond scales. For many sources, radiative transfer models are capable of satisfactorily reproducing the observations. They are described by density powerlaw distributions (n(r) ~ r^(-p)) with p = 1.0 +/-0.25. Such distributions are shallower than those found on larger scales probed with single-dish (sub)mm studies. Other sources have density laws that are shallower/steeper than p = 1.0 and there is evidence that these MYSOs are viewed near edge-on or near face-on, respectively. The images also reveal a diffuse component tracing somewhat larger scale structures, particularly visible in the regions S140, AFGL 2136, IRAS 20126+4104, Mon R2, and Cep A. We thus find a flattening of the MYSO envelope density law going from ~10 000 AU down to scales of ~1000 AU. We propose that this may be evidence of rotational support of the envelope (abridged).
We present interferometric and single-dish mid-infrared observations of a sample of massive young stellar objects (BN-type objects), using VLTI-MIDI (10 micron) and Subaru-COMICS (24.5 micron). We discuss the regions S140, Mon R2, M8E-IR, and W33A. T
he observations probe the inner regions of the dusty envelope at scales of 50 milli arcsecond and 0.6 arcsec (100-1000 AU), respectively. Simultaneous model fits to spectral energy distributions and spatial data are achieved using self-consistent spherical envelope modelling. We conclude that those MYSO envelopes that are best described by a spherical geometry, the commensurate density distribution is a powerlaw with index -1.0. Such a powerlaw is predicted if the envelope is supported by turbulence on the 100-1000AU scales probed with MIDI and COMICS, but the role of rotation at these spatial scales need testing.
We present high-resolution observations made with the Very Large Array (VLA) in its A configuration at frequencies between 5 and 43 GHz of a sample of five massive young stellar objects (YSOs): LkHa101, NGC2024-IRS2, S106-IR, W75N and S140-IRS1. The
resolution varied from 0.04 arcsec (at 43 GHz) to 0.5 arcsec (at 5 GHz), corresponding to a linear resolution as high as 17 AU for our nearest source. A MERLIN observation of S106-IR at 23 GHz with 0.03-arcsec resolution is also presented. S106-IR and S140-IRS1 are elongated at 43 GHz perpendicular to their large scale bipolar outflows. This confirms the equatorial wind picture for these sources seen previously in MERLIN 5 GHz observations. The other sources are marginally resolved at 43 GHz. The spectral indices we derive for the sources in our sample range from +0.2 to +0.8, generally consistent with ionized stellar winds. We have modelled our sources as uniform, isothermal spherical winds, with LkHa101 and NGC2024-IRS2 yielding the best fits. However, in all cases our fits give wind temperatures of only 2000 to 5000 K, much less than the effective temperatures of main-sequence stars of the same luminosity, a result which is likely due to the clumpy nature of the winds.
