We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the tadpole, a small globule in the Carina Nebula that hosts the HH 900 jet+outflow system. Our data include $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1, $^{13}$CO, C$^{18}$O J=3-2, and serendipitous detections of DCN J=3-2 and CS J=7-6. With angular resolution comparable to the Hubble Space Telescope (HST), our data reveal for the first time the bipolar molecular outflow in CO, seen only inside the globule, that is launched from the previously unseen jet-driving protostar (the HH 900 YSO). The biconical morphology joins smoothly with the externally irradiated outflow seen in ionized gas tracers outside the globule, tracing the overall morphology of a jet-driven molecular outflow. Continuum emission at the location of the HH 900 YSO appears to be slightly flattened perpendicular to outflow axis. Model fits to the continuum have a best-fit spectral index of $\sim 2$, suggesting cold dust and the onset of grain growth. In position-velocity space, $^{13}$CO and C$^{18}$O gas kinematics trace a C-shaped morphology, similar to infall profiles seen in other sources, although the global dynamical behaviour of the gas remains unclear. Line profiles of the CO isotopologues display features consistent with externally heated gas. We estimate a globule mass of $\sim 1.9$ M$_{\odot}$, indicating a remaining lifetime of $\sim 4$ Myr, assuming a constant photoevaporation rate. This long globule lifetime will shield the disk from external irradiation perhaps prolonging its life and enabling planet formation in regions where disks are typically rapidly destroyed.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

Forming high-mass stars have a significant effect on their natal environment. Their feedback pathways, including winds, outflows, and ionising radiation, shape the evolution of their surroundings which impacts the formation of the next generation of stars. They create or reveal dense pillars of gas and dust towards the edges of the cavities they clear. They are modelled in feedback simulations, and the sizes and shapes of the pillars produced are consistent with those observed. However, these models predict measurably different kinematics which provides testable discriminants. Here we present the first ALMA Compact Array (ACA) survey of 13 pillars in Carina, observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=2-1, and the 230 GHz continuum. The pillars in this survey were chosen to cover a wide range in properties relating to the amount and direction of incident radiation, proximity to nearby irradiating clusters and cloud rims, and whether they are detached from the cloud. With these data, we are able to discriminate between models. We generally find pillar velocity dispersions of $<$ 1 km s$^{-1}$ and that the outer few layers of molecular emission in these pillars show no significant offsets from each other, suggesting little bulk internal motions within the pillars. There are instances where the pillars are offset in velocity from their parental cloud rim, and some with no offset, hinting at a stochastic development of these motions.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

We present new MUSE/VLT observations of a small globule in the Carina H II region that hosts the HH 900 jet+outflow system. Data were obtained with the GALACSI ground-layer adaptive optics system in wide-field mode, providing spatially-resolved maps of diagnostic emission lines. These allow us to measure the variation of the physical properties in the globule and jet+outflow system. We find high temperatures ($T_e \approx 10^4$ K), modest extinction ($A_V \approx 2.5$ mag), and modest electron densities ($n_e \approx 200$ cm$^{-3}$) in the ionized gas. Higher excitation lines trace the ionized outflow; both the excitation and ionization in the outflow increase with distance from the opaque globule. In contrast, lower excitation lines that are collisionally de-excited at densities $\gtrsim 10^4$ cm$^{-3}$ trace the highly collimated protostellar jet. Assuming the globule is an isothermal sphere confined by the pressure of the ionization front, we compute a Bonnor-Ebert mass of $\sim 3.7$ M$_{\odot}$. This is two orders of magnitude higher than previous mass estimates, calling into question whether small globules like the Tadpole contribute to the bottom of the IMF. Derived globule properties are consistent with a cloud that has been and/or will be compressed by the ionization front on its surface. At the estimated globule photoevaporation rate of $\sim 5 \times 10^{-7}$ M$_{\odot}$ yr$^{-1}$, the globule will be completely ablated in $\sim 7$ Myr. Stars that form in globules like the Tadpole will emerge into the H II later and may help resolve some of the temporal tension between disk survival and enrichment.

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.

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

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.

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

Context: Simulations suggest that gas heating due to radiative feedback is a key factor in whether or not multiple protostellar systems will form. Chemistry is a good tracer of the physical structure of a protostellar system, since it depends on the temperature structure. Aims: We aim to study the relationship between envelope gas temperature and protostellar multiplicity. Methods: Single dish observations of various molecules that trace the cold, warm and UV-irradiated gas are used to probe the temperature structure of multiple and single protostellar systems on 7000 AU scales. Results: Single, close binary, and wide multiples present similar current envelope gas temperatures, as estimated from H$_2$CO and DCO$^+$ line ratios. The temperature of the outflow cavity, traced by c-C$_3$H$_2$, on the other hand, shows a relation with bolometric luminosity and an anti-correlation with envelope mass. Although the envelope gas temperatures are similar for all objects surveyed, wide multiples tend to exhibit a more massive reservoir of cold gas compared to close binary and single protostars. Conclusions: Although the sample of protostellar systems is small, the results suggest that gas temperature may not have a strong impact on fragmentation. We propose that mass, and density, may instead be key factors in fragmentation.

The JCMT Gould Belt Survey was one of the first Legacy Surveys with the James Clerk Maxwell Telescope in Hawaii, mapping 47 square degrees of nearby (< 500 pc) molecular clouds in both dust continuum emission at 850 $\mu$m and 450 $\mu$m, as well as a more-limited area in lines of various CO isotopologues. While molecular clouds and the material that forms stars have structures on many size scales, their larger-scale structures are difficult to observe reliably in the submillimetre regime using ground-based facilities. In this paper, we quantify the extent to which three subsequent data-reduction methods employed by the JCMT GBS accurately recover emission structures of various size scales, in particular, dense cores which are the focus of many GBS science goals. With our current best data-reduction procedure, we expect to recover 100% of structures with Gaussian sigma sizes of $\le$30" and intensity peaks of at least five times the local noise for isolated peaks of emission. The measured sizes and peak fluxes of these compact structures are reliable (within 15% of the input values), but source recovery and reliability both decrease significantly for larger emission structures and for fainter peaks. Additional factors such as source crowding have not been tested in our analysis. The most recent JCMT GBS data release includes pointing corrections, and we demonstrate that these tend to decrease the sizes and increase the peak intensities of compact sources in our dataset, mostly at a low level (several percent), but occasionally with notable improvement.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

We present full spectral scans from 200-670$\mu$m of 26 Class 0+I protostellar sources, obtained with $Herschel$-SPIRE, as part of the "COPS-SPIRE" Open Time program, complementary to the DIGIT and WISH Key programs. Based on our nearly continuous, line-free spectra from 200-670 $\mu$m, the calculated bolometric luminosities ($L_{\rm bol}$) increase by 50% on average, and the bolometric temperatures ($T_{\rm bol}$) decrease by 10% on average, in comparison with the measurements without Herschel. Fifteen protostars have the same Class using $T_{\rm bol}$ and $L_{\rm bol}$/$L_{\rm submm}$. We identify rotational transitions of CO lines from J=4-3 to J=13-12, along with emission lines of $^{13}$CO, HCO$^+$, H$_{2}$O, and [CI]. The ratios of $^{12}$CO to $^{13}$CO indicate that $^{12}$CO emission remains optically thick for $J_{\rm up}$ < 13. We fit up to four components of temperature from the rotational diagram with flexible break points to separate the components. The distribution of rotational temperatures shows a primary population around 100 K with a secondary population at $\sim$350 K. We quantify the correlations of each line pair found in our dataset, and find the strength of correlation of CO lines decreases as the difference between $J$-level between two CO lines increases. The multiple origins of CO emission previously revealed by velocity-resolved profiles are consistent with this smooth distribution if each physical component contributes to a wide range of CO lines with significant overlap in the CO ladder. We investigate the spatial extent of CO emission and find that the morphology is more centrally peaked and less bipolar at high-$J$ lines. We find the CO emission observed with SPIRE related to outflows, which consists two components, the entrained gas and shocked gas, as revealed by our rotational diagram analysis as well as the studies with velocity-resolved CO emission.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars (Class 0/I). The Herschel-PACS spectral maps covering $\sim55-210$ $\mu$m with a field-of-view of $\sim$50'' are used to quantify the gas excitation conditions and spatial extent using rotational transitions of H$_{2}$O, high-$J$ CO, and OH, as well as [O I] and [C II]. We confirm that a warm ($\sim$300 K) CO reservoir is ubiquitous and that a hotter component ($760\pm170$ K) is frequently detected around protostars. The line emission is extended beyond $\sim$1000 AU spatial scales in 40/90 objects, typically in molecular tracers in Class 0 and atomic tracers in Class I objects. High-velocity emission ($\gtrsim90$ km s$^{-1}$) is detected in only 10 sources in the [O I] line, suggesting that the bulk of [O I] arises from gas that is moving slower than typical jets. Line flux ratios show an excellent agreement with models of $C$-shocks illuminated by UV photons for pre-shock densities of $\sim$$10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-IR molecular and atomic lines are a unique diagnostic of feedback from UV emission and shocks in envelopes of deeply embedded protostars.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$\alpha$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$\alpha$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$\alpha$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region.

We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

(Abridged) Through spectrally unresolved observations of high-J CO transitions, Herschel-PACS has revealed large reservoirs of warm (300 K) and hot (700 K) molecular gas around low-mass protostars. We aim to shed light on the excitation and origin of the CO ladder observed toward protostars, and on the water abundance in different physical components using spectrally resolved Herschel-HIFI data. Observations are presented of the highly excited CO line J=16-15 with Herschel-HIFI toward 24 low-mass protostellar objects. The spectrally resolved profiles show two distinct velocity components: a broad component with an average FWHM of 20 km/s, and a narrower component with a FWHM of 5 km/s that is often offset from the source velocity. The average rotational temperature over the entire profile, as measured from comparison between CO J=16-15 and 10-9 emission, is ~300 K. A radiative-transfer analysis shows that the average H2O/CO column-density ratio is ~0.02, suggesting a total H2O abundance of ~2x10^-6. Two distinct velocity profiles observed in the HIFI line profiles suggest that the CO ladder observed with PACS consists of two excitation components. The warm component (300 K) is associated with the broad HIFI component, and the hot component (700 K) is associated with the offset HIFI component. The former originates in either outflow cavity shocks or the disk wind, and the latter in irradiated shocks. The ubiquity of the warm and hot CO components suggests that fundamental mechanisms govern the excitation of these components; we hypothesize that the warm component arises when H2 stops being the dominant coolant. In this scenario, the hot component arises in cooling molecular H2-poor gas just prior to the onset of H2 formation. High spectral resolution observations of highly excited CO transitions uniquely shed light on the origin of warm and hot gas in low-mass protostellar objects.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

In molecular outflows from forming low-mass protostars, most oxygen is expected to be locked up in water. However, Herschel observations have shown that typically an order of magnitude or more of the oxygen is still unaccounted for. To test if the oxygen is instead in atomic form, SOFIA-GREAT observed the R1 position of the bright molecular outflow from NGC1333-IRAS4A. The [OI] 63 um line is detected and spectrally resolved. From an intensity peak at +15 km/s, the intensity decreases until +50 km/s. The profile is similar to that of high-velocity (HV) H2O and CO 16-15, the latter observed simultaneously with [OI]. A radiative transfer analysis suggests that ~15% of the oxygen is in atomic form toward this shock position. The CO abundance is inferred to be ~10^-4 by a similar analysis, suggesting that this is the dominant oxygen carrier in the HV component. These results demonstrate that a large portion of the observed [OI] emission is part of the outflow. Further observations are required to verify whether this is a general trend.

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.

[Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms{}) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,{\sc i}] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,{\sc i}] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,{\sc i}] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

We present the physical and evolutionary properties of prestellar and protostellar clumps in the Herschel Infrared GALactic plane survey (Hi-GAL) in two large areas centered in the Galactic plane and covering the tips of the long Galactic bar at the intersection with the spiral arms. The areas fall in the longitude ranges 19 < l < 33 and 340 < l < 350, while latitude is -1 < b < 1. Newly formed high mass stars and prestellar objects are identified and their properties derived and compared. A study is also presented on five giant molecular complexes at the further edge of the bar. The star-formation rate was estimated from the quantity of proto-stars expected to form during the collapse of massive turbulent clumps into star clusters. This new method was developed by applying a Monte Carlo procedure to an evolutionary model of turbulent cores and takes into account the wide multiplicity of sources produced during the collapse. The star-formation rate density values at the tips are 1.2 +- 0.3 10-3 Msol/yr/kpc2 and 1.5+-0.3 10-3 Msol/yr/kpc2 in the first and fourth quadrant, respectively. The same values estimated on the entire field of view, that is including the tips of the bar and background and foreground regions, are 0.9+-0.2 10-3 Msol/yr/kpc2 and 0.8+-0.2 10-3 Msol/yr/kpc2. The conversion efficiency is approximately 0.8% in the first quadrant and 0.5% in the fourth quadrant, and does not show a significant difference in proximity of the bar. The star forming regions identified through CO contours at the further edge of the bar show star-formation rate densities larger than the surrounding regions but their conversion efficiencies are comparable. Our results suggest that the star-formation activity at the bar is due to a large amount of dust and molecular material rather than being due to a triggering process.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

We present an analysis of the effect of feedback from O- and B-type stars with data from the integral field spectrograph MUSE mounted on the Very Large Telescope of pillar-like structures in the Carina Nebular Complex, one of the most massive star-forming regions in the Galaxy. For the observed pillars, we compute gas electron densities and temperatures maps, produce integrated line and velocity maps of the ionised gas, study the ionisation fronts at the pillar tips, analyse the properties of the single regions, and detect two ionised jets originating from two distinct pillar tips. For each pillar tip we determine the incident ionising photon flux $Q_\mathrm{0,pil}$ originating from the nearby massive O- and B-type stars and compute the mass-loss rate $\dot{M}$ of the pillar tips due to photo-evaporation caused by the incident ionising radiation. We combine the results of the Carina data set with archival MUSE data of a pillar in NGC 3603 and with previously published MUSE data of the Pillars of Creation in M16, and with a total of 10 analysed pillars, find tight correlations between the ionising photon flux and the electron density, the electron density and the distance from the ionising sources, and the ionising photon flux and the mass-loss rate. The combined MUSE data sets of pillars in regions with different physical conditions and stellar content therefore yield an empirical quantification of the feedback effects of ionising radiation. In agreement with models, we find that $\dot{M}\propto Q_\mathrm{0,pil}^{1/2}$.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself. Using the ALMA Science Verification data of HL Tau, we investigate the high velocity molecular gas being removed from the system as a result of the star formation process. We aim to place these motions in context with the optically detected jet, and the disk. With these high resolution ($\sim 1"$) ALMA observations of CO (J=1-0), we quantify the outwards motions of the molecular gas. We find evidence for a bipolar outwards flow, with an opening angle, as measured in the red-shifted lobe, starting off at 90$^\circ$, and narrowing to 60$^\circ$ further from the disk, likely because of magnetic collimation. Its outwards velocity, corrected for inclination angle is of order 2.4 km s$^{-1}$.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar.

The JCMT Gould Belt Legacy Survey obtained SCUBA-2 observations of dense cores within three sub-regions of Orion B: LDN 1622, NGC 2023/2024, and NGC 2068/2071, all of which contain clusters of cores. We present an analysis of the clustering properties of these cores, including the two-point correlation function and Cartwright's Q parameter. We identify individual clusters of dense cores across all three regions using a minimal spanning tree technique, and find that in each cluster, the most massive cores tend to be centrally located. We also apply the independent M-Sigma technique and find a strong correlation between core mass and the local surface density of cores. These two lines of evidence jointly suggest that some amount of mass segregation in clusters has happened already at the dense core stage.

We present the $^{13}$CO/C$^{18}$O (J=3-2) Heterodyne Inner Milky Way Plane Survey (CHIMPS) which has been carried out using the Heterodyne Array Receiver Program on the 15 m James Clerk Maxwell Telescope (JCMT) in Hawaii. The high-resolution spectral survey currently covers |b| < 0.5 deg and 28 < l < 46 deg, with an angular resolution of 15 arcsec in 0.5 km/s velocity channels. The spectra have a median rms of $\sim$ 0.6 K at this resolution, and for optically thin gas at an excitation temperature of 10 K, this sensitivity corresponds to column densities of $N_{\mathrm{H}_{2}} \sim 3 \times 10^{20}\,$cm$^{-2}$ and $N_{\mathrm{H}_{2}} \sim 4 \times 10^{21}\,$cm$^{-2}$ for $^{13}$CO and C$^{18}$O, respectively. The molecular gas that CHIMPS traces is at higher column densities and is also more optically thin than in other publicly available CO surveys due to its rarer isotopologues, and thus more representative of the three-dimensional structure of the clouds. The critical density of the J=3-2 transition of CO is $\gtrsim 10^{4}$ cm$^{-3}$ at temperatures of $\leq 20$ K, and so the higher density gas associated with star formation is well traced. These data complement other existing Galactic plane surveys, especially the JCMT Galactic Plane Survey which has similar spatial resolution and column density sensitivity, and the Herschel infrared Galactic Plane Survey. In this paper, we discuss the observations, data reduction and characteristics of the survey, presenting integrated emission maps for the region covered. Position-velocity diagrams allow comparison with Galactic structure models of the Milky Way, and while we find good agreement with a particular four arm model, there are some significant deviations.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

Context. The classification of young stellar objects (YSOs) is typically done using the infrared spectral slope or bolometric temperature, but either can result in contamination of samples. More accurate methods to determine the evolutionary stage of YSOs will improve the reliability of statistics for the embedded YSO population and provide more robust stage lifetimes. Aims. We aim to separate the truly embedded YSOs from more evolved sources. Methods. Maps of HCO+ J=4-3 and C18O J=3-2 were observed with HARP on the James Clerk Maxwell Telescope (JCMT) for a sample of 56 candidate YSOs in Perseus and Taurus in order to characterize emission from high (column) density gas. These are supplemented with archival dust continuum maps observed with SCUBA on the JCMT and Herschel PACS to compare the morphology of the gas and dust in the protostellar envelopes. The spatial concentration of HCO+ J=4-3 and 850 micron dust emission are used to classify the embedded nature of YSOs. Results. Approximately 30% of Class 0+I sources in Perseus and Taurus are not Stage I, but are likely to be more evolved Stage II pre-main sequence (PMS) stars with disks. An additional 16% are confused sources with an uncertain evolutionary stage. Conclusions. Separating classifications by cloud reveals that a high percentage of the Class 0+I sources in the Perseus star forming region are truly embedded Stage I sources (71%), while the Taurus cloud hosts a majority of evolved PMS stars with disks (68%). The concentration factor method is useful to correct misidentified embedded YSOs, yielding higher accuracy for YSO population statistics and Stage timescales. Current estimates (0.54 Myr) may overpredict the Stage I lifetime on the order of 30%, resulting in timescales of 0.38 Myr for the embedded phase.

Water probes the dynamics in young stellar objects (YSOs) effectively, especially shocks in molecular outflows. It is a key molecule for exploring whether the physical properties of low-mass protostars can be extrapolated to massive YSOs. As part of the WISH key programme, we investigate the dynamics and the excitation conditions of shocks along the outflow cavity wall as function of source luminosity. Velocity-resolved Herschel-HIFI spectra of the H2O 988, 752, 1097 GHz and 12CO J=10-9, 16-15 lines were analysed for 52 YSOs with bolometric luminosities (L_bol) ranging from <1 to >10^5 L_sun. The profiles of the H2O lines are similar, indicating that they probe the same gas. We see two main Gaussian emission components in all YSOs: a broad component associated with non-dissociative shocks in the outflow cavity wall (cavity shocks) and a narrow component associated with quiescent envelope material. More than 60% of the total integrated intensity of the H2O lines (L_H2O) comes from the cavity shock component. The H2O line widths are similar for all YSOs, whereas those of 12CO 10-9 increase slightly with L_bol. The excitation analysis of the cavity shock component, performed with the non-LTE radiative transfer code RADEX, shows stronger 752 GHz emission for high-mass YSOs, likely due to pumping by an infrared radiation field. As previously found for CO, a strong correlation with slope unity is measured between log(L_H2O) and log(L_bol), which can be extrapolated to extragalactic sources. We conclude that the broad component of H2O and high-J CO lines originate in shocks in the outflow cavity walls for all YSOs, whereas lower-J CO transitions mostly trace entrained outflow gas. The higher UV field and turbulent motions in high-mass objects compared to their low-mass counterparts may explain the slightly different kinematical properties of 12CO 10-9 and H2O lines from low- to high-mass YSOs.

Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artifacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software, Starlink, but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only 3 times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.

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.

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.

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 present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

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 N$_{2}$H$^{+}$ 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 virialised. 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 N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, 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 southwest to northeast, indicating sequential star formation across the region.

Context: Outflows are an important part of the star formation process as both the result of ongoing active accretion and one of the main sources of mechanical feedback on small scales. Water is the ideal tracer of these effects because it is present in high abundance in various parts of the protostar. Method: We present \textit{Herschel} HIFI spectra of multiple water-transitions towards 29 nearby Class 0/I protostars as part of the WISH Survey. These are decomposed into different Gaussian components, with each related to one of three parts of the protostellar system; quiescent envelope, cavity shock and spot shocks in the jet and at the base of the outflow. We then constrain the excitation conditions present in the two outflow-related components. Results: Water emission is optically thick but effectively thin, with line ratios that do not vary with velocity, in contrast to CO. The physical conditions of the cavity and spot shocks are similar, with post-shock H$_{2}$ densities of order 10$^{5}-$10$^{8}$\,cm$^{-3}$ and H$_{2}$O column densities of order 10$^{16}-$10$^{18}$\,cm$^{-2}$. H$_{2}$O emission originates in compact emitting regions: for the spot shocks these correspond to point sources with radii of order 10-200\,AU, while for the cavity shocks these come from a thin layer along the outflow cavity wall with thickness of order 1-30\,AU. Conclusions: Water emission at the source position traces two distinct kinematic components in the outflow; J shocks at the base of the outflow or in the jet, and C shocks in a thin layer in the cavity wall. Class I sources have similar excitation conditions to Class 0 sources, but generally smaller line-widths and emitting region sizes. We suggest that it is the velocity of the wind driving the outflow, rather than the decrease in envelope density or mass, that is the cause of the decrease in H$_{2}$O intensity between Class 0 and I.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

The bandwith, sensitivity and sheer survey speed of the SKA offers unique potential for deep spectroscopic surveys of the Milky Way. Within the frequency bands available to the SKA lie many transitions that trace the ionised, radical and molecular components of the interstellar medium and which will revolutionise our understanding of many physical processes. In this chapter we describe the impact on our understanding of the Milky Way that can be achieved by spectroscopic SKA surveys, including "out of the box" early science with radio recombination lines, Phase 1 surveys of the molecular ISM using anomalous formaldehyde absorption, and full SKA surveys of ammonia inversion lines.

We present CO, H2, HI and HISA distributions from a set of simulations of grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate the emission of the 2nd Galactic Quadrant by placing the observer inside the modelled galaxies and post process the simulations using a radiative transfer code, so as to create synthetic observations. We compare the synthetic datacubes to observations of the 2nd Quadrant of the Milky Way to test the ability of the current models to reproduce the basic chemistry of the Galactic ISM, as well as to test how sensitive such galaxy models are to different recipes of chemistry and/or feedback. We find that models which include feedback and self-gravity can reproduce the production of CO with respect to H2 as observed in our Galaxy, as well as the distribution of the material perpendicular to the Galactic plane. While changes in the chemistry/feedback recipes do not have a huge impact on the statistical properties of the chemistry in the simulated galaxies, we find that the inclusion of both feedback and self-gravity are crucial ingredients, as our test without feedback failed to reproduce all of the observables. Finally, even though the transition from H2 to CO seems to be robust, we find that all models seem to underproduce molecular gas, and have a lower molecular to atomic gas fraction than is observed. Nevertheless, our fiducial model with feedback and self-gravity has shown to be robust in reproducing the statistical properties of the basic molecular gas components of the ISM in our Galaxy.

Aims: Our aim is to determine the critical parameters in water chemistry and the contribution of water to the oxygen budget by observing and modelling water gas and ice for a sample of eleven low-mass protostars, for which both forms of water have been observed. Methods: A simplified chemistry network, which is benchmarked against more sophisticated chemical networks, is developed that includes the necessary ingredients to determine the water vapour and ice abundance profiles in the cold, outer envelope in which the temperature increases towards the protostar. Comparing the results from this chemical network to observations of water emission lines and previously published water ice column densities, allows us to probe the influence of various agents (e.g., FUV field, initial abundances, timescales, and kinematics). Results: The observed water ice abundances with respect to hydrogen nuclei in our sample are 30-80ppm, and therefore contain only 10-30% of the volatile oxygen budget of 320ppm. The keys to reproduce this result are a low initial water ice abundance after the pre-collapse phase together with the fact that atomic oxygen cannot freeze-out and form water ice in regions with T(dust)>15 K. This requires short prestellar core lifetimes of less than about 0.1Myr. The water vapour profile is shaped through the interplay of FUV photodesorption, photodissociation, and freeze-out. The water vapour line profiles are an invaluable tracer for the FUV photon flux and envelope kinematics. Conclusions: The finding that only a fraction of the oxygen budget is locked in water ice can be explained either by a short pre-collapse time of less than 0.1 Myr at densities of n(H)~1e4 cm-3, or by some other process that resets the initial water ice abundance for the post-collapse phase. A key for the understanding of the water ice abundance is the binding energy of atomic oxygen on ice.

Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario.

We present the study of the H2O spatial distribution at two bright shocked regions along IRAS4A, one of the strongest H2O emitters among the Class 0 outflows. We obtained Herschel-PACS maps of the IRAS4A outflow and HIFI observations of two shocked positions. The largest HIFI beam of 38 arcsec at 557 GHz was mapped in several key water lines with different upper energy levels, to reveal possible spatial variations of the line profiles. We detect four H2O lines and CO (16-15) at the two selected positions. In addition, transitions from related outflow and envelope tracers are detected. Different gas components associated with the shock are identified in the H2O emission. In particular, at the head of the red lobe of the outflow, two distinct gas components with different excitation conditions are distinguished in the HIFI emission maps: a compact component, detected in the ground-state water lines, and a more extended one. Assuming that these two components correspond to two different temperature components observed in previous H2O and CO studies, the excitation analysis of the H2O emission suggests that the compact (about 3 arcsec) component is associated with a hot (T~1000 K) gas with densities ~(1-4)x10^5 cm^{-3}, whereas the extended one (10-17 arcsec) traces a warm (T~300-500 K) and dense gas (~(3-5)x10^7 cm^{-3}). Finally, using the CO (16-15) emission observed at R2, we estimate the H2O/H2 abundance of the warm and hot components to be (7-10)x10^{-7} and (3-7)x10^{-5}. Our data allowed us, for the first time, to resolve spatially the two temperature components previously observed with HIFI and PACS. We propose that the compact hot component may be associated with the jet that impacts the surrounding material, whereas the warm, dense, and extended component originates from the compression of the ambient gas by the propagating flow.

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.

The major physical processes responsible for shaping and sculpting pillars in the clouds surrounding massive stars (i.e. the `Pillars of Creation') are now being robustly incorporated into models quantifying the ionizing radiation from massive stars. The detailed gas dynamics within these pillars can now be compared with observations. Our goal is to quantify the gas dynamics in a pillar being sculpted by a nearby massive star. To do this, we use the CO, $^{13}$CO, and C$^{18}$O J=1-0 emission towards a pillar in the Vulpecula Rift. These data are a combination of CARMA and FCRAO observations providing high resolution ($\sim5''$) imaging of large scale pillar structures ($>100''$). We find that this cold ($\sim18$ K), low density material ($8\times10^3$ cm$^{-3}$) material is fragmenting on Jeans scales, has very low velocity dispersions ($\sim0.5$ km s$^{-1}$), and appears to be moving away from the ionizing source. We are able to draw direct comparisons with three models from the literature, and find that those with lower velocity dispersions best fit our data, although the dynamics of any one model do not completely agree with our observations. We do however, find that our observed pillar exhibits many of the characteristics expected from simulations.

Water in outflows from protostars originates either as a result of gas-phase synthesis from atomic oxygen at T > 200 K, or from sputtered ice mantles containing water ice. We aim to quantify the contribution of the two mechanisms that lead to water in outflows, by comparing observations of gas-phase water to methanol (a grain surface product) towards three low-mass protostars in NGC1333. In doing so, we also quantify the amount of methanol destroyed in outflows. To do this, we make use of JCMT and Herschel-HIFI data of H2O, CH3OH and CO emission lines and compare them to RADEX non-LTE excitation simulations. We find up to one order of magnitude decrease in the column density ratio of CH3OH over H2O as the velocity increases in the line wings up to ~15 km/s. An independent decrease in X(CH3OH) with respect to CO of up to one order of magnitude is also found in these objects. We conclude that gas-phase formation of H2O must be active at high velocities (above 10 km/s, relative to the source velocity) to re-form the water destroyed during sputtering. In addition, the transition from sputtered water at low velocities to formed water at high velocities must be gradual. We place an upper limit of two orders of magnitude on the destruction of methanol by sputtering effects.

(Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are <~ 100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii > 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.

Aims. The aim of this paper is to study deuterated water in the solar-type protostars NGC1333 IRAS4A and IRAS4B, to compare their HDO abundance distribution with other star-forming regions, and to 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 and 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.

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.

Abridged abstract: For stars to form, material must fall inwards from core scales through the envelope towards the central protostar. The velocity profile around protostars is poorly constrained. 6 Class 0 protostars and one Class I protostars observed with HIFI on board Herschel as part of the "Water in star-forming regions with Herschel" (WISH) survey show infall signatures in water line observations. 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 these sources. We assume a free-fall velocity profile and, having found the best fit, vary the radii over which infall takes place. 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 4 sources infall takes place on core to envelope scales (i.e. ~10000-3000 AU). In 2 sources 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. The non-detection of excited water in the inner envelope in 6/7 protostars is further evidence that water must be heavily depleted from the gas-phase at these radii. Infall in four of the sources is supersonic and infall in all sources 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 IRAS4A is large (>~10^-4 msol\yr), higher than the mass outflow rate and expected mass accretion rates onto the star, suggesting that any flattened disk-like structure on small scales will be gravitationally unstable, potentially leading to rotational fragmentation and/or episodic accretion.

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.

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

The outflow force of molecular bipolar outflows is a key parameter in theories of young stellar feedback on their surroundings. The focus of many outflow studies is the correlation between the outflow force, bolometric luminosity and envelope mass. However, it is difficult to combine the results of different studies in large evolutionary plots over many orders of magnitude due to the range of data quality, analysis methods and corrections for observational effects such as opacity and inclination. We aim to determine the outflow force for a sample of low luminosity embedded sources. We will quantify the influence of the analysis method and the assumptions entering the calculation of the outflow force. We use the James Clerk Maxwell Telescope to map 12CO J=3-2 over 2'x2' regions around 16 Class I sources of a well-defined sample in Ophiuchus at 15" resolution. The outflow force is then calculated using seven different methods differing e.g. in the use of intensity-weighted emission and correction factors for inclination. The results from the analysis methods differ from each other by up to a factor of 6, whereas observational properties and choices in the analysis procedure affect the outflow force by up to a factor of 4. For the sample of Class I objects, bipolar outflows are detected around 13 sources including 5 new detections, where the three non-detections are confused by nearby outflows from other sources. When combining outflow forces from different studies, a scatter by up to a factor of 5 can be expected. Although the true outflow force remains unknown, the separation method (separate calculation of dynamical time and momentum) is least affected by the uncertain observational parameters. The correlations between outflow force, bolometric luminosity and envelope mass are further confirmed down to low luminosity sources.

Disk winds have been postulated as a mechanism for angular momentum release in protostellar systems for decades. HD 163296 is a Herbig Ae star surrounded by a disk and has been shown to host a series of HH knots (HH 409) with bow shocks associated with the farthest knots. Here we present ALMA Science Verification data of CO J=2-1 and J=3-2 emission which are spatially coincident with the blue shifted jet of HH knots, and offset from the disk by -18.6 km/s. The emission has a double corkscrew morphology and extends more than 10" from the disk with embedded emission clumps coincident with jet knots. We interpret this double corkscrew as emission from material in a molecular disk wind, and that the compact emission near the jet knots is being heated by the jet which is moving at much higher velocities. We show that the J=3-2 emission is likely heavily filtered by the interferometer, but the J=2-1 emission suffers less due to the larger beam and measurable angular scales. Excitation analysis suggests temperatures exceeding 900 K in these compact features, with the wind mass, momentum and energy being of order 10^{-5} M_sun, 10^{-4} M_sun km/s and 10^{40} erg respectively. The high mass loss rate suggests that this star is dispersing the disk faster than it is funneling mass onto the star.

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.

Context. G29.96-0.02 is a high-mass star-forming cloud observed at 70, 160, 250, 350, and 500 microns as part of the Herschel survey of the Galactic Plane during the Science Demonstration Phase. Aims. We wish to conduct a far-infrared study of the sources associated with this star-forming region by estimating their physical properties and evolutionary stage, and investigating the clump mass function, the star formation efficiency and rate in the cloud. Methods. We have identified the Hi-GAL sources associated with the cloud, searched for possible counterparts at centimeter and infrared wavelengths, fitted their spectral energy distribution and estimated their physical parameters. Results. A total of 198 sources have been detected in all 5 Hi-GAL bands, 117 of which are associated with 24 microns emission and 87 of which are not associated with 24 microns emission. We called the former sources 24 microns-bright and the latter ones 24 microns-dark. The [70-160] color of the 24 microns-dark sources is smaller than that of the 24 microns-bright ones. The 24 microns-dark sources have lower L_bol and L_bol/M_env than the 24 microns-bright ones for similar M_env, which suggests that they are in an earlier evolutionary phase. The G29-SFR cloud is associated with 10 NVSS sources and with extended centimeter continuum emission well correlated with the 70 microns emission. Most of the NVSS sources appear to be early B or late O-type stars. The most massive and luminous Hi-GAL sources in the cloud are located close to the G29-UC region, which suggests that there is a privileged area for massive star formation towards the center of the G29-SFR cloud. Almost all the Hi-GAL sources have masses well above the Jeans mass but only 5% have masses above the virial mass, which indicates that most of the sources are stable against gravitational collapse. The sources with M_env > M_virial and that ...

Context; Our understanding of the star formation process has traditionally been confined to certain mass or luminosity boundaries because most studies focus only on low-, intermediate- or high-mass star-forming regions. As part of the "Water In Star-forming regions with Herschel" (WISH) key program, water and other important molecules, such as CO and OH, have been observed in 51 embedded young stellar objects (YSOs). The studied sample covers a range of luminosities from <1 to >10^5 L_sol. Aims; We analyse the CO line emission towards a large sample of protostars in terms of both line intensities and profiles. Methods; Herschel-HIFI spectra of the 12CO 10-9, 13CO 10-9 and C18O 5-4, 9-8 and 10-9 lines are analysed for a sample of 51 YSOs. In addition, JCMT spectra of 12CO 3-2 and C18O 3-2 extend this analysis to cooler gas components. Results; All observed CO and isotopologue spectra show a strong linear correlation between the logarithms of the line and bolometric luminosities across six orders of magnitude on both axes. This suggests that the high-J CO lines primarily trace the amount of dense gas associated with YSOs. This relation can be extended to larger (extragalactic) scales. The majority of the detected 12CO line profiles can be decomposed into a broad and a narrow Gaussian component, while the C18O spectra are mainly fitted with a single Gaussian. A broadening of the line profile is also observed from pre-stellar cores to embedded protostars, which is due mostly to non-thermal motions (turbulence/infall). The widths of the broad 12CO 3-2 and 10-9 velocity components correlate with those of the narrow C18O 9-8 profiles, suggesting that the entrained outflowing gas and envelope motions are related independent of the mass of the protostar. These results indicate that physical processes in protostellar envelopes have similar characteristics across the studied luminosity range.

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.

Observations of higher-excited transitions of abundant molecules such as CO are important for determining where energy in the form of shocks is fed back into the parental envelope of forming stars. The nearby prototypical and protobinary low-mass hot core, IRAS16293-2422 (I16293) is ideal for such a study. The source was targeted with ALMA for science verification purposes in band 9, which includes CO J=6-5 (E_up/k_B ~ 116 K), at an unprecedented spatial resolution (~0.2", 25 AU). I16293 itself is composed of two sources, A and B, with a projected distance of 5". CO J=6-5 emission is detected throughout the region, particularly in small, arcsecond-sized hotspots, where the outflow interacts with the envelope. The observations only recover a fraction of the emission in the line wings when compared to data from single-dish telescopes, with a higher fraction of emission recovered at higher velocities. The very high angular resolution of these new data reveal that a bow shock from source A coincides, in the plane of the sky, with the position of source B. Source B, on the other hand, does not show current outflow activity. In this region, outflow entrainment takes place over large spatial scales, >~ 100 AU, and in small discrete knots. This unique dataset shows that the combination of a high-temperature tracer (e.g., CO J=6-5) and very high angular resolution observations is crucial for interpreting the structure of the warm inner environment of low-mass protostars.

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.

In a previous paper we identified cores within infrared dark clouds (IRDCs). We regarded those without embedded sources as the least evolved, and labelled them starless. Here we identify the most isolated starless cores and model them using a three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive the cores' physical parameters and discuss the relation between the mass, temperature, density, size and the surrounding interstellar radiation field (ISRF) for the cores. The masses of the cores were found not to correlate with their radial size or central density. The temperature at the surface of a core was seen to depend almost entirely on the level of the ISRF surrounding the core. No correlation was found between the temperature at the centre of a core and its local ISRF. This was seen to depend, instead, on the density and mass of the core.

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-infrared (8 {\mu}m) by Spitzer (we refer to these as 'Spitzer-dark' regions). They all lie in the range l=300 - 330 \circ and |b| 6 1 \circ. Of these, only 1205 were seen in emission in the far-infrared (250-500 {\mu}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-infrared, but will also be seen in emission in the far-infrared at the longest Herschel wavebands (250-500 {\mu}m). If a region is dark at all wavelengths throughout the mid-infrared and far-infrared, then it is most likely to be simply a region of lower background infrared emission (a 'hole in the sky'). Hence, it appears that previous surveys, based on Spitzer and other mid-infrared data alone, may have over-estimated 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 {\mu}m, and found that 972 of them had at least one clearly defined 250-{\mu}m peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-{\mu}m point source somewhere within the cloud, 149 (15 per cent) contained a 24-{\mu}m point source but no 8-{\mu}m source, and 170 (18 per cent) contained no 24-{\mu}m or 8-{\mu}m point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

We present $^{12}$CO, $^{13}$CO and C$^{18}$O (J=3$-$2) observations of a new cluster of outflows in the Vulpecula Rift with HARP-B on the JCMT. The mass associated with the outflows, measured using the $^{12}$CO HARP-B observations and assuming a distance to the region of 2.3 kpc, is 129 \msol{}, while the mass associated with the dense gas from C$^{18}$O observations is 458 \msol{} and the associated sub-millimeter core has a mass of 327 $\pm$ 112 \msol{} independently determined from Bolocam 1.1mm data. The outflow-to-core mass ratio is therefore $\sim$0.4, making this region one of the most efficient observed thus far with more than an order of magnitude more mass in the outflow than would be expected based on previous results. The kinetic energy associated with the flows, 94$\times10^{45}$ ergs, is enough to drive the turbulence in the local clump, and potentially unbind the local region altogether. The detection of SiO (J=8$-$7) emission toward the outflows indicates that the flow is still active, and not simply a fossil flow. We also model the SEDs of the four YSOs associated with the molecular material, finding them all to be of mid to early B spectral type. The energetic nature of the outflows and significant reservoir of cold dust detected in the sub-mm suggest that these intermediate mass YSOs will continue to accrete and become massive, rather than reach the main sequence at their current mass.

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.

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.

There is currently no accepted theoretical framework for the formation of the most massive stars, and the manner in which protostars continue to accrete and grow in mass beyond \sim10Msun is still a controversial topic. In this study we use several prescriptions of stellar accretion and a description of the Galactic gas distribution to simulate the luminosities and spatial distribution of massive protostellar population of the Galaxy. We then compare the observables of each simulation to the results of the Red MSX Source (RMS) survey, a recently compiled database of massive young stellar objects. We find that the observations are best matched by accretion rates which increase as the protostar grows in mass, such as those predicted by the turbulent core and competitive accretion (i.e. Bondi-Hoyle) models. These 'accelerating accretion' models provide very good qualitative and quantitative fits to the data, though we are unable to distinguish between these two models on our simulations alone. We rule out models with accretion rates which are constant with time, and those which are initially very high and which fall away with time, as these produce results which are quantitatively and/or qualitatively incompatible with the observations. To simultaneously match the low- and high-luminosity YSO distribution we require the inclusion of a 'swollen-star' pre-main-sequence phase, the length of which is well-described by the Kelvin-Helmholz timescale. Our results suggest that the lifetime of the YSO phase is \sim 10^5yrs, whereas the compact Hii-region phase lasts between \sim 2 - 4 \times 10^5yrs depending on the final mass of the star. Finally, the absolute numbers of YSOs are best matched by a globally averaged star-formation rate for the Galaxy of 1.5-2Msun/yr.

We use a Minimum Spanning Tree algorithm to characterize the spatial distribution of Galactic Far-IR sources and derive their clustering properties. We aim to reveal the spatial imprint of different types of star forming processes, e.g. isolated spontaneous fragmentation of dense molecular clouds, or events of triggered star formation around HII regions, and highlight global properties of star formation in the Galaxy. We plan to exploit the entire Hi-GAL survey of the inner Galactic plane to gather significant statistics on the clustering properties of star forming regions, and to look for possible correlations with source properties such as mass, temperature or evolutionary stage. In this paper we present a pilot study based on the two 2x2 square degree fields centered at longitudes l=30 and l=59 obtained during the Science Demonstration Phase (SDP) of the Herschel mission. We find that over half of the clustered sources are associated with HII regions and infrared dark clouds. Our analysis also reveals a smooth chromatic evolution of the spatial distribution where sources detected at short-wavelengths, likely proto-stars surrounded by warm circumstellar material emitting in the far-infrared, tend to be clustered in dense and compact groups around HII regions while sources detected at long-wavelengths, presumably cold and dusty density enhancements of the ISM emitting in the sub-millimeter, are distributed in larger and looser groups.

We present a determination of the luminosity functions of massive young stellar objects (MYSOs) and compact (C)HII regions within the Milky Way Galaxy using the large, well-selected sample of these sources identified by the Red MSX Source (RMS) survey. The MYSO luminosity function decreases monotonically such that there are few with $L\gtrsim 10^{5}$Lsol, whilst the CHII regions are detected up to ~10$^{6}Lsol. The lifetimes of these phases are also calculated as a function of luminosity by comparison with the luminosity function for local main-sequence OB stars. These indicate that the MYSO phase has a duration ranging from 4x10$^{5}$ yrs for 10$^{4}$Lsol to ~7x10$^{4}$ yrs at 10$^{5}$Lsol, whilst the CHII region phase lasts of order 3x10$^{5}$ yrs or ~3-10% of the exciting star's main-sequence lifetime. MYSOs between 10$^{4} Lsol and ~10$^{5}$ Lsol are massive but do not display the radio continuum or near-IR \HI{} recombination line emission indicative of an HII region, consistent with being swollen due to high ongoing or recent accretion rates. Above ~10$^{5}$ Lsol the MYSO phase lifetime becomes comparable to the main-sequence Kelvin-Helmholtz timescale, at which point the central star can rapidly contract onto the main-sequence even if still accreting, and ionise a CHII region, thus explaining why few highly luminous MYSOs are observed.

Context: Despite their profound effect on the universe, the formation of massive stars and stellar clusters remains elusive. In the past decade, compelling evidence has emerged that suggests Infrared Dark Clouds (IRDCs) may be precursors to stellar clusters. However, the usual method for identifying IRDCs is biased by the requirement that they are seen in absorption against bright mid-IR emission, whereas dust continuum observations allow cold, dense pre-stellar-clusters to be identified anywhere. Aims: We aim to understand what physical properties characterize IRDCs, to explore the population of dust continuum sources that are not IRDCs, and to roughly characterize the star formation activity in dust continuum sources. Results: We present temperature and column density maps in the Hi-GAL l=30 and l=59 SDP fields, as well as a robust algorithm for cirrus subtraction and source identification using Hi-GAL. We report on the fraction of Hi-GAL sources which are mid-IR-dark, -neutral, or -bright in both fields. We find significant trends in column density and temperature between mid-IR-dark and -bright pixels; mid-IR-dark pixels are about 10 K colder and have a factor of 2 higher column density on average than mid-IR-bright pixels. We find that Hi-GAL dust continuum sources span a range of evolutionary states from pre- to star-forming, and that warmer sources are associated with more star formation tracers. There is a trend of increasing temperature with tracer type from mid-IR-dark at the coldest, to outflow/maser sources in the middle, and finally to 8 and 24 micron bright sources at the warmest. Finally, we identify five candidate IRDC-like sources on the far-side of the Galaxy. These are cold (~ 20 K), high column density (N(H2) > 10^22 cm^-2) clouds identified with Hi-GAL which, despite bright surrounding mid-IR emission, show little to no absorption at 8 micron. (abridged)

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi-wavelength, spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G031.03+00.26 and G031.03+00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We found masses ranging from 90 to 290 solar masses with temperatures from 8 to 11K at the centre of each core and 18 to 28K at the surface. The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high mass star formation having yet occurred in three of our cores.

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

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.

The Exeter FCRAO CO Galactic Plane Survey consists of 12CO and 13CO (J=1-0) observations over the galactic plane covering 55 degrees <= l <= 102 degrees, |b| >= 1 degree and 141 degrees <= l <= 195 degrees, -3.5 degrees <= b <= 5.5 degrees with a spatial resolution of ~45" and a spectral resolution of ~0.15km/s. We will present the methodology of a threshold-based cloud and clump determination method which retains hierarchical information, then discuss associating sources with clouds in the catalogue. Once complete, this catalogue of clouds and clumps will encompass the majority of the Northern Galactic Plane, providing knowledge of the molecular structure of the galaxy and the starting point for studies of the variation in star formation efficiency. In addition, it will allow us to identify clouds that have no or little star formation taking place inside them, which are often overlooked in the study of the conditions required for star formation to take place.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.

We present a first study of the star-forming compact dust condensations revealed by Herschel in the two 2 \times 2 \degr Galactic Plane fields centered at [l;b] = [30\degr; 0 \degr] and [l;b] = [59\degr; 0 \degr], respectively, and observed during the Science Demonstration Phase for the Herschel infrared Galactic Plane survey (Hi-GAL) Key-Project. Compact source catalogs extracted for the two fields in the five Hi-GAL bands (70, 160, 250, 350 and 500 $\mu$m) were merged based on simple criteria of positional association and spectral energy distribution (SED) consistency into a final catalog which contains only coherent SEDs with counterparts in at least three adjacent Herschel bands. These final source lists contain 528 entries for the l = 30\degr field, and 444 entries for the l = 59\degr field. The SED coverage has been augmented with ancillary data at 24 $\mu$m and 1.1 mm. SED modeling for the subset of 318 and 101 sources (in the two fields, respectively) for which the distance is known was carried out using both a structured star/disk/envelope radiative transfer model and a simple isothermal grey-body. Global parameters like mass, luminosity, temperature and dust properties have been estimated. The Lbol/Menv ratio spans four orders of magnitudes from values compatible with the pre-protostellar phase to embedded massive zero-age main sequence stars. Sources in the l = 59\degr field have on average lower L/M, possibly outlining an overall earlier evolutionary stage with respect to the sources in the l = 30\degr field. Many of these cores are actively forming high-mass stars, although the estimated core surface densities appear to be an order of magnitude below the 1 g cm$^{-2}$ critical threshold for high-mass star formation.

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

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.

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.

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

We present medium resolution spectropolarimetry and long term photo-polarimetry of two massive post-red supergiants, IRC +10420 and HD 179821. The data provide new information on their circumstellar material as well as their evolution. In IRC +10420, the polarization of the Halpha line is different to that of the continuum, which indicates that the electron-scattering region is not spherically symmetric. The observed long term changes in the polarimetry can be associated with an axi-symmetric structure, along the short axis of the extended reflection nebulosity. Long term photometry reveals that the star increased in temperature until the mid-nineties, after which the photospheric flux in the optical levelled off. As the photometric changes are mostly probed in the red, they do not trace high stellar temperatures sensitively. And so, it is not obvious whether the star has halted its increase in temperature or not. For HD 179821 we find no polarization effects across any absorption or emission lines, but observe very large polarization changes of order 5% over 15 years. Unexpectedly, during the same period, the optical photometry displayed modest variability at the 0.2 magnitude level. Several explanations for this puzzling fact are discussed. Most of which, involving asymmetries in the circumstellar material, seem to fail as there is no evidence for the presence of hot, dusty material close to the star. Alternatively, the variations can be explained by the presence of a non-radially pulsating photosphere. Changes in the photometry hint at an increase in temperature corresponding to a change through two spectral subclasses over the past ten years.

We report on resolved interferometric observations with VLTI/MIDI of the massive young stellar object (MYSO) W33A. The MIDI observations deliver spectrally dispersed visibilities with values between 0.03 and 0.06, for a baseline of 45m over the wavelength range 8-13 micron. The visibilities indicate that W33A has a FWHM size of approximately 120AU (0.030'') at 8 micron which increases to 240AU at 13 micron, scales previously unexplored among MYSOs. This observed trend is consistent with the temperature falling off with distance. 1D dust radiative transfer models are simultaneously fit to the visibility spectrum, the strong silicate feature and the shape of the mid infrared spectral energy distribution (SED). For any powerlaw density distribution, we find that the sizes (as implied by the visibilities) and the stellar luminosity are incompatible. A reduction to a third of W33A's previously adopted luminosity is required to match the visibilities; such a reduction is consistent with new high resolution 70 micron data from Spitzer's MIPSGAL survey. We obtain best fits for models with shallow dust density distributions of r^(-0.5) and r^(-1.0) and for increased optical depth in the silicate feature produced by decreasing the ISM ratio of graphite to silicates and using optical grain properties by Ossenkopf et al. (1992).

Abridged abstract: The Red MSX Source (RMS) survey is an ongoing effort to return a large, well-selected sample of massive young stellar objects (MYSOs) within our Galaxy. A series of ground-based follow-up observations are being undertaken in order to remove contaminant objects from our list of 2000 candidates, and to begin characterising these MYSOs. As a part of these follow-up observations, high resolution (~1") mid-IR imaging aids the identification of contaminant objects which are resolved (UCHII regions, PN) as opposed to those which are unresolved (YSOs, evolved stars) as well as identifying YSOs near UCHII regions and other multiple sources. We present 10.4 micron imaging observations for 346 candidate MYSOs in the RMS survey in the Southern Hemisphere, primarily outside the region covered by the GLIMPSE Spitzer Legacy Survey. These were obtained using TIMMI2 on the ESO 3.6m telescope in La Silla, Chile. Our photometric accuracy is of order 0.05Jy, and our astrometric accuracy is 0.8", which is an improvement over the nominal 2" accuracy of the MSX PSC.

We present linear spectropolarimetric data for eight Herbig Be and four Herbig Ae stars at H alpha, H beta and H gamma. Changes in the linear polarisation are detected across all Balmer lines for a large fraction of the observed objects, confirming that the small-scale regions surrounding these objects are flattened (i.e. disk-like). Furthermore, all objects with detections show similar characteristics at the three spectral lines, despite differences in transition probability and optical depth going from H alpha to H gamma. A large fraction of early Herbig Be stars (B0-B3) observed show line depolarisation effects. However the early Herbig Ae stars (A0-A2), observed for comparison, show intrinsic line polarisation signatures. Our data suggest that the popular magnetic accretion scenario for T Tauri objects may be extended to Herbig Ae stars, but that it may not be extended to early Herbig Be stars, for which the available data are consistent with disc accretion.

The Red MSX Source (RMS) survey is a multi-wavelength programme of follow-up observations designed to distinguish between genuine massive young stellar objects (MYSOs) and other embedded or dusty objects, such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe), from a sample of ~2000 MYSOs candidates. These were identified by comparing the colours of objects from the MSX and 2MASS point source catalogues to those of known MYSOs, in order to develop colour selection criteria which have been used to produced the RMS sample of MYSOs candidates. Our ultimate aim is to produce a large unbiased sample of MYSOs (~500) with complementary multi-wavelength data with which to study their properties. Here we report the results of 826 radio continuum and 346 mid-IR imaging observations carried out using the Australia Telescope Compact Array and TIMMI2 on the ESO 3.6 m telescope respectively. These observations are aimed at identifying and removing contaminating sources.

We present new medium-resolution spectropolarimetric observations of HD 45677 in the B and R-bands. A change in polarisation is detected across Halpha, Hbeta and Hgamma confirming that the ionised region around the star is aspherical. (Q,U) points associated with these emission lines occur away from the continuum, defining a polarisation vector which points in the same direction for each of the lines at an average intrinsic polarisation angle of 164 +/- 3 degrees. These data were combined with past photometric and polarimetric data from the literature to investigate any variability. We find that HD 45677 is both photometrically and polarimetrically variable and that these changes are linked. We suggest that these variations may be caused by an aspherical blowout and by deriving a least-squares fit to the B-band polarimetric data in Q-U space, we show that the blowout occurs at an intrinsic polarisation angle of 175 +/- 1 degrees, along the same angle as the proposed geometry of the ionised region.

The Red MSX Source (RMS) survey is a multi-wavelength program of follow-up observations designed to distinguish between genuine massive young stellar objects (MYSOs) and other embedded or dusty objects, such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe). We have identified nearly 2000 massive YSOs candidates by comparing the colours of MSX and 2MASS point sources to those of known MYSOs. Unfortunately, there are several other types of embedded or dust enshrouded objects that have similar colours as MYSOs and contaminate our sample. Two sources of contamination are from UCHII regions and PNe, both of which can be identified from the radio emission emitted by their ionised nebulae. In order to identify UCHII regions and PNe that contaminate our sample we have conducted high resolution radio continuum observations at 3.6 and 6 cm of all southern MYSOs candidates ($235< l < 350$) using the Australia Telescope Compact Array (ATCA). These observations have a spatial resolution of $\sim$1--2\arcsec and typical image rms noise values of $\sim$0.3 mJy -- sensitive enough to detect a B0.5 star at the far side of the Galaxy. Of the 826 RMS sources observed we found 199 to be associated with radio emission, $\sim$25% of the sample. More significant than those sources associated with radio emission are the 627 RMS sources for which no radio emission was detected as these sources are still potential MYSOs. The Galactic distribution, morphologies and spectral indices of the radio sources associated with the RMS sources are consistent with these sources being UCHII regions. In addition to the 826 RMS sources observed we present observations of a further 190 fields. These observations were made towards MSX sources that passed cuts in earlier versions of the survey, but were later excluded.