J. C. Mottram

THE RED MSX SOURCE SURVEY: THE MASSIVE YOUNG STELLAR POPULATION OF OUR GALAXY

We present the Red MSX Source survey, the largest statistically selected catalog of young massive protostars and H II regions to date. We outline the construction of the catalog using mid- and near-infrared color selection. We also discuss the detailed follow up work at other wavelengths, including higher spatial resolution data 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 arcsec. We briefly summarize some of the results that can be obtained from studying the properties of the objects in the catalog as a whole; we find evidence that the most massive stars form: (1) preferentially nearer the Galactic center than the anti-center; (2) in the most heavily reddened environments, suggestive of high accretion rates; and (3) from the most massive cloud cores.

ATLASGAL — towards a complete sample of massive star forming clumps ⋆

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 ◦ <l< 350 ◦ and 10 ◦ <l< 60 ◦ with| b|< 1.5 ◦ . 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 sig nificant di fferences 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 s tar 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-regio n 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 asso ciated with the most massive clumps, probably due to very short pre-main sequence lifetimes in the most luminous sources.

The Herschel and JCMT Gould Belt Surveys: Constraining Dust Properties in the Perseus B1 Clump with PACS, SPIRE, and SCUBA-2

We present Herschel observations from the Herschel Gould Belt Survey and SCUBA-2 science verification observations from the JCMT Gould Belt Survey of the B1 clump in the Perseus molecular cloud. We determined the dust emissivity index using four different techniques to combine the Herschel PACS+SPIRE data at 160 - 500 microns with the SCUBA-2 data at 450 microns and 850 microns. Of our four techniques, we found the most robust method was to filter-out the large-scale emission in the Herschel bands to match the spatial scales recovered by the SCUBA-2 reduction pipeline. Using this method, we find beta ~ 2 towards the filament region and moderately dense material and lower beta values (beta > 1.6) towards the dense protostellar cores, possibly due to dust grain growth. We find that beta and temperature are more robust with the inclusion of the SCUBA-2 data, improving estimates from Herschel data alone by factors of ~ 2 for beta and by ~ 40% for temperature. Furthermore, we find core mass differences of < 30% compared to Herschel-only estimates with an adopted beta = 2, highlighting the necessity of long wavelength submillimeter data for deriving accurate masses of prestellar and protostellar cores.

The RMS Survey: Ammonia and water maser analysis of massive star forming regions. ⋆

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 H 2O masers which corresponds to an overall detection rate of∼50 per cent. We find no significant di fference in the detection rate for Hii regions and MYSOs which would suggest that the conditions required to produce maser emission are equally likely in both phases. Comparing the detection rates as a function of luminosity we find the H 2O detection rate has a positive dependence on the source luminosity, with the detection rate increasing with increasing luminosity. We detect ammonia emission towards 479 of these massive young stars, which corresponds to∼80 per cent. 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 −1 and 2× 10 15 cm −2 , respectively. We find that the NH 3 (1,1) line width and kinetic temperature are correlated with luminosity and finding no underlying dep endence 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.

The JCMT Gould Belt Survey: first results from the SCUBA-2 observations of the Ophiuchus molecular cloud and a virial analysis of its prestellar core population

In this paper, we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N2H+ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialized. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N2H+ is a good tracer of the bound material of prestellar cores, although we find some evidence for N2H+ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C18O and that traced by N2H+, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor–Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from south-west to north-east, indicating sequential star formation across the region.

Waterfalls around protostars - Infall motions towards Class 0/I envelopes as probed by water

Context. For stars to form, material must fall inwards from core scales through the envelope towards the central protostar. While theories of how this takes place have been around for some time, the velocity profile around protostars is poorly constrai ned. The combination of observations in multiple transitions of a tr acer which is sensitive to kinematics and radiative transfe r modelling of those lines has the potential to break this deadlock. Aims. Seven protostars observed with the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory as part of the “Water in star-forming regions with Herschel” (WISH) survey show infall signatures in water line observations. We aim to constrain the infall velocity and the radii over which inf all is taking place within the protostellar envelopes of the se sources. We will also use these data to constrain the chemistry of cold water. Methods. We use 1-D non-LTE RATRAN radiative transfer models of the observed water lines to constrain the infall velocity and chemistry in the protostellar envelopes of 6 Class 0 protost ars and one Class I source. We assume a free-fall velocity profile and, having found the best fit, vary the radii over which infall tak es place. Results. In the well-studied Class 0 protostar NGC1333-IRAS4A we find that infall takes place over the whole envelope to which our observations are sensitive (r&1000 AU). For L1527, L1157, BHR71 and IRAS15398 infall takesplace on core to envelope scales (i.e. ∼10000−3000 AU). In Serpens-SMM4 and GSS30 the inverse P-Cygni profiles seen in the ground-state lines are more likely due to larger-scale motions or foreground clouds. Models including a simple consideration of the chemistry are consistent with the observations, while using step abundance profiles are not. T he non-detection of excited water in the inner envelope in 6 out of 7 protostars is further evidence that water must be heavily depleted from the gas-phase at these radii. Conclusions. Infall in four of the sources is supersonic and infall in all s ources must take place at the outer edge of the envelope, which may be evidence that collapse is global or outside-in rather than inside-out. The mass infall rate in NGC1333-IRAS4A is large (&10 −4 M⊙ yr −1 ), higher than the mass outflow rate and expected mass accreti on rates onto the star, suggesting that any flattened disk-like structure on small scales will be gravitationall y unstable, potentially leading to rotational fragmentati on and/or episodic accretion.

The RMS survey: near-IR spectroscopy of massive young stellar objects

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 per cent) of the targets are YSOs, of which 131 are massive YSOs (L BOL > 5×10 3 L ⊙ , M > 8 M ⊙ ). 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 Brγ, H 2 , fluorescent Fe ii, CO bandhead, [Fe ii] and He i 2-1 1 S- 1 P, in order of frequency of occurrence. In total, ∼40 per cent of the YSOs display either fluorescent Fe ii 1.6878 μm 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 per cent of the sources exhibit [Fe ii] or H 2 emission, indicating the presence of an outflow. Three quarters of all sources have Brγ in emission. A good correlation with bolometric luminosity was observed for both the Brγ and H 2 emission line strengths, covering 1 < L BOL < 3.5 × 10 5 L ⊙ . 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.

CO bandhead emission of massive young stellar objects: determining disc properties

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 μm spectra of 20 massive young stellar objects (MYSOs) from the Red MSX Source (RMS) data base, 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 data base, 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.

Cores in infrared dark clouds (IRDCs) seen in the Hi-GAL survey between l=300 degrees and 330 degrees

We have used data taken as part of the Herschel infrared Galactic Plane survey (Hi-GAL) to study 3171 infrared dark cloud (IRDC) candidates that were identified in the mid-IR (8 μm) by Spitzer (we refer to these as ‘Spitzer-dark’ regions). They all lie in the range l= 300–330° and |b|≤ 1°. Of these, only 1205 were seen in emission in the far-IR (250–500 μm) by Herschel (we call these ‘Herschel-bright’ clouds). It is predicted that a dense cloud will not only be seen in absorption in the mid-IR, but will also be seen in emission in the far-IR at the longest Herschel wavebands (250–500 μm). If a region is dark at all wavelengths throughout the mid-IR and far-IR, then it is most likely to be simply a region of lower background IR emission (a ‘hole in the sky’). Hence, it appears that previous surveys, based on Spitzer and other mid-IR data alone, may have overestimated the total IRDC population by a factor of ∼2. This has implications for estimates of the star formation rate in IRDCs in the Galaxy. We studied the 1205 Herschel-bright IRDCs at 250 μm and found that 972 of them had at least one clearly defined 250-μm peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-μm point source somewhere within the cloud, 149 (15 per cent) contained a 24-μm point source but no 8-μm source and 170 (18 per cent) contained no 24- or 8-μm point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

Deuterated water in the solar-type protostars NGC 1333 IRAS 4A and IRAS 4B

Aims. The aim of this paper is to study deuterated water in the solar-type protostars NGC1333 IRAS4A and IRAS4B, compare their HDO abundance distribution with other star-forming regions and constrain their HDO/H2O ratios. Methods. Using the Herschel/HIFI instrument as well as ground-based telescopes, we observed several HDO lines covering a large excitation range (Eup/k=22-168 K) towards these protostars and an outflow position. Non-LTE radiative transfer codes were then used to determine the HDO abundance profiles in these sources. Results. The HDO fundamental line profiles show a very broad component, tracing the molecular outflows, in addition to a narrower emission component as well as a narrow absorbing component. In the protostellar envelope of NGC1333 IRAS4A, the HDO inner (T>100 K) and outer (T<100 K) abundances with respect to H2 are estimated at 7.5x10^{-9} and 1.2x10^{-11} respectively, whereas, in NGC1333 IRAS4B, they are 1.0x10^{-8} and 1.2x10^{-10} respectively. Similarly to the low-mass protostar IRAS16293-2422, an absorbing outer layer with an enhanced abundance of deuterated water is required to reproduce the absorbing components seen in the fundamental lines at 465 and 894 GHz in both sources. This water-rich layer is probably extended enough to encompass the two sources as well as parts of the outflows. In the outflows emanating from NGC1333 IRAS4A, the HDO column density is estimated at about (2-4)x10^{13} cm^{-2}, leading to an abundance of about (0.7-1.9)x10^{-9}. An HDO/H2O ratio between 7x10^{-4} and 9x10^{-2} is derived in the outflows. In the warm inner regions of these two sources, we estimate the HDO/H2O ratios at about 1x10^{-4}-4x10^{-3}. This ratio seems higher (a few %) in the cold envelope of IRAS4A, whose possible origin is discussed in relation to formation processes of HDO and H2O.