The mass growth of protostars is a central element to the determination of fundamental stellar population properties such as the initial mass function. Constraining the accretion history of individual protostars is therefore an important aspect of star formation research. The goal of the study presented here is to determine whether high-mass (proto)stars gain their mass from a compact (<0.1pc) fixed-mass reservoir of gas, often referred to as dense cores, in which they are embedded, or whether the mass growth of high-mass stars is governed by the dynamical evolution of the parsec-scale clump that typically surrounds them. To achieve this goal, we performed a 350micron continuum mapping of 11 infrared dark clouds, along side some of their neighbouring clumps, with the ArT\'eMiS camera on APEX. By identifying about 200 compact ArT\'eMiS sources, and matching them with Herschel Hi-GAL 70micron sources, we have been able to produce mass vs. temperature diagrams. We compare the nature (i.e. starless or protostellar) and location of the ArT\'eMiS sources in these diagrams with modelled evolutionary tracks of both core-fed and clump-fed accretion scenarios. We argue that the latter provide a better agreement with the observed distribution of high-mass star-forming cores. However, a robust and definitive conclusion on the question of the accretion history of high-mass stars requires larger number statistics.

The initial stage of star formation is a complex area study because of its high density and low temperature. Under such conditions, many molecules become depleted from the gas phase by freezing out onto dust grains. However, the deuterated species could remain gaseous and are thus ideal tracers. We investigate the gas dynamics and NH$_2$D chemistry in eight massive pre/protocluster clumps. We present NH$_2$D 1$_{11}$-1$_{01}$ (at 85.926 GHz), NH$_3$ (1, 1) and (2, 2) observations in the eight clumps using the PdBI and the VLA, respectively. We find that the distribution between deuterium fractionation and kinetic temperature shows a number density peak at around $T_{\rm kin}=16.1$ K, and the NH$_2$D cores are mainly located at a temperature range of 13.0 to 22.0 K. We detect seven instances of extremely high deuterium fractionation of $1.0 \leqslant D_{\rm frac} \leqslant 1.41$. We find that the NH$_2$D emission does not appear to coincide exactly with either dust continuum or NH$_3$ peak positions, but often surrounds the star-formation active regions. This suggests that the NH$_{2}$D has been destroyed by the central young stellar object (YSO) due to its heating. The detected NH$_2$D lines are very narrow with a median width of $\rm 0.98\pm0.02 km/s$. The extracted NH$_2$D cores are gravitationally bound ($\alpha_{\rm vir} < 1$), are likely prestellar or starless, and can potentially form intermediate-mass or high-mass stars. Using NH$_3$ (1, 1) as a dynamical tracer, we find very complicated dynamical movement, which can be explained by a combined process with outflow, rotation, convergent flow, collision, large velocity gradient, and rotating toroids. High deuterium fractionation strongly depends on the temperature condition. NH$_2$D is a poor evolutionary indicator of high-mass star formation in evolved stages, but a useful tracer in the starless and prestellar cores.

Massive stars and their associated ionized (HII) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H regions are still poorly known. To investigate the impact of HII regions on the formation of filaments, we imaged the Galactic HII region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. We used the ArT\'eMiS camera on the APEX telescope and combined the ArT\'eMiS data at 350 and 450 microns with Herschel-SPIRE/HOBYS. We studied the dense gas distribution around RCW 120 with a resolution of 8 arcsec (0.05 pc at a distance of 1.34 kpc). Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the HII region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer like profile with a deconvolved Gaussian FWHM of 0.09 pc, as observed for filaments in low-mass star-forming regions. This study suggests that compression exerted by HII regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArT\'eMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure

Classical hot cores are rich in molecular emission, and they show a high abundance of complex organic molecules (COMs). The emergence of molecular complexity is poorly constrained in the early evolution of hot cores. Using the Atacama Large Millimeter Array we put observational constraints on the physical location of COMs in a high-mass protostellar envelope associated with the G328.2551-0.5321 clump. The protostar is single down to ~400au scales and we resolve the emission region of COMs. Using thermodynamic equilibrium modelling of the available 7.5 GHz bandwidth around ~345 GHz, we detect emission from 10 COMs, and identify a line of deuterated water (HDO). The most extended emission originates from methanol, methyl formate and formamide. Together with HDO, these molecules are found to be associated with both the accretion shocks and the inner envelope, for which we estimate a moderate temperature of $T_{\rm kin}\sim$110 K. Our findings reveal a significant difference in the distribution of COMs. O-bearing COMs, such as ethanol, acetone, and ethylene glycol are almost exclusively found and show a higher abundance towards the accretion shocks with $T_{\rm kin}\sim$180 K. Whereas N-bearing COMs with a CN group, such as vinyl and ethyl cyanide peak on the central position, thus the protostar and the accretion disk. This is the first observational evidence for a large column density of COMs seen towards accretion shocks at the centrifugal barrier at the inner envelope. Since the molecular composition is dominated by that of the accretion shocks and the radiatively heated hot inner region is very compact, we propose this source to be a precursor to a classical, radiatively heated hot core.

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 ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

The structure formation of the dense interstellar material and the fragmentation of clumps into cores is a fundamental step to understand how stars and stellar clusters form. We aim to establish a statistical view of clump fragmentation at sub-parsec scales based on a large sample of massive clumps selected from the ATLASGAL survey. We used the APEX/SABOCA camera at 350 $\mu$m to image clumps at a resolution of 8.$''$5. The majority of the sample consists of massive clumps that are weak or in absorption at 24 $\mu$m. We resolve rich filamentary structures and identify the population of compact sources. We use association with mid-infrared 22-24 $\mu$m and 70 $\mu$m point sources to pin down the star formation activity of the cores. We then statistically assess their physical properties, and the fragmentation characteristics of massive clumps. We find a moderate correlation between the clump fragmentation levels with the clump gas density and the predicted number of fragments with pure Jeans fragmentation scenario; we find a strong correlation between the mass of the most massive fragment and the total clump mass, suggesting that the self-gravity may play an important role in the clumps' small scale structure formation. We identify 27 massive quiescent cores with $M_{\rm core}>100$ M$_{\odot}$ within 5 kpc; these are massive enough to be self-gravitating but do not yet show any sign of star-formation. This sample comprises, therefore, promising candidates of massive pre-stellar cores, or deeply embedded high-mass protostars.

Context The Vela Molecular Ridge is one of the nearest (700 pc) giant molecular cloud (GMC) complexes hosting intermediate-mass (up to early B, late O stars) star formation, and is located in the outer Galaxy, inside the Galactic plane. Vela C is one of the GMCs making up the Vela Molecular Ridge, and exhibits both sub-regions of robust and sub-regions of more quiescent star formation activity, with both low- and intermediate(high)-mass star formation in progress. Aims We aim to study the individual and global properties of dense dust cores in Vela C, and aim to search for spatial variations in these properties which could be related to different environmental properties and/or evolutionary stages in the various sub-regions of Vela C. Methods We mapped the submillimetre (345 GHz) emission from vela C with LABOCA (beam size 19.2", spatial resolution ~0.07 pc at 700 pc) at the APEX telescope. We used the clump-finding algorithm CuTEx to identify the compact submillimetre sources. We also used SIMBA (250 GHz) observations, and Herschel and WISE ancillary data. The association with WISE red sources allowed the protostellar and starless cores to be separated, whereas the Herschel dataset allowed the dust temperature to be derived for a fraction of cores. The protostellar and starless core mass functions (CMFs) were constructed following two different approaches, achieving a mass completeness limit of 3.7 Msun. Results We retrieved 549 submillimetre cores, 316 of which are starless and mostly gravitationally bound (therefore prestellar in nature). Both the protostellar and the starless CMFs are consistent with the shape of a Salpeter initial mass function in the high-mass part of the distribution. Clustering of cores at scales of 1--6 pc is also found, hinting at fractionation of magnetised, turbulent gas.

Massive clumps tend to fragment into clusters of cores and condensations, some of which form high-mass stars. In this work, we study the structure of massive clumps at different scales, analyze the fragmentation process, and investigate the possibility that star formation is triggered by nearby HII regions. We present a high angular resolution study of a sample of 8 massive proto-cluster clumps. Combining infrared data, we use few-arcsecond resolution radio- and millimeter interferometric data to study their fragmentation and evolution. Our sample is unique in the sense that all the clumps have neighboring HII regions. Taking advantage of that, we test triggered star formation using a novel method where we study the alignment of the centres of mass traced by dust emission at multiple scales. The eight massive clumps have masses ranging from 228 to 2279 $M_\odot$. The brightest compact structures within infrared bright clumps are typically associated with embedded compact radio continuum sources. The smaller scale structures of $R_{\rm eff}$ $\sim$ 0.02 pc observed within each clump are mostly gravitationally bound and massive enough to form at least a B3-B0 type star. Many condensations have masses larger than 8 $M_\odot$ at small scale of $R_{\rm eff}$ $\sim$ 0.02 pc. Although the clumps are mostly infrared quiet, the dynamical movements are active at clump scale ($\sim$ 1 pc). We studied the spatial distribution of the gas conditions detected at different scales. For some sources we find hints of external triggering, whereas for others we find no significant pattern that indicates triggering is dynamically unimportant. This probably indicates that the different clumps go through different evolutionary paths. In this respect, studies with larger samples are highly desired.

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.

The GLOSTAR survey will study the star formation in the Galactic plane between $-2^{\circ}< \textit{l}< 85^{\circ}$ and $|b| < 1^{\circ}$ with unprecedented sensitivity in both, flux density ($\sim$ 40 $\mu Jy$ beam$^{-1}$) and range of angular scales ($\sim$ 1."5 to the largest radio structures in the Galaxy). In this paper we present the first results obtained from a radio continuum map of a 16 square degree sized region of the Galactic plane centered on $\textit{l} = 32^{\circ}$ and $b = 0^{\circ}$ ($28^{\circ} < \textit{l} < 36^{\circ}$ and $|b| < 1^{\circ}$). This map has a resolution of 18" and sensitivity of $\sim$ 60-150 $\mu Jy$ beam$^{-1}$. We present data acquired with the VLA in D-configuration. Two 1 GHz wide sub-bands were observed simultaneously and centred at 4.7 and 6.9 GHz. These data were calibrated and imaged using the $\textit{Obit}$ software package. The source extraction has been performed using the BLOBCAT software package and verified through a combination of visual inspection and cross-matching with other radio and mid-infrared surveys. The final catalog consists of 1575 discrete radio sources and 27 large scale structures (including W43 and W44). By cross-matching with other catalogs and calculating the spectral indices ($S(\nu) \propto \nu^\alpha$), we have classified 231 continuum sources as HII regions, 37 as ionization fronts, and 46 as planetary nebulae. The longitude and latitude distribution and negative spectral indices are all consistent with the vast majority of the unclassified sources being extragalactic background sources. We present a catalog of 1575 radio continuum sources and discuss their physical properties, emission nature and relation with previously reported. These first GLOSTAR results have increased the number of reliable HII regions in this part of the Galaxy by a factor of four.

The accretion-driven outbursts of young FU Orionis-type stars may be a common stage of pre-main sequence evolution and can have a significant impact on the circumstellar environment as it pertains to the growth of solids and eventually planets. This episodic accretion is thought to be sustained by additional gas in-falling from the circumstellar envelope and disk. We present APEX observations of the CO gas in the envelope around V883 Orionis, a young outbursting star. The observations mapped the $^{12}$CO(4-3), $^{12}$CO(3-2), and $^{13}$CO(3-2) lines with the FLASH$^{+}$ instrument and the $^{12}$CO(6-5) line with the SEPIA instrument. We detected high signal-to-noise emission extending out to radii $>10000$ au and calculated integrated fluxes of $1100~\rm Jy~km~s^{-1}$ for $^{12}$CO(6-5), $2400~\rm Jy~km~s^{-1}$ for $^{12}$CO(4-3), $1600~\rm Jy~km~s^{-1}$ for $^{12}$CO(3-2), and $450~\rm Jy~km~s^{-1}$ for $^{13}$CO(3-2). We used the thermo-chemical code P{\small RO}D{\small I}M{\small O} to test several models and find the data are best described by an envelope structure with $\rm M_{env}\approx 0.2-0.4\,M_{\odot}$ and a mass-infall rate of $\rm \dot{M}_{inf}=1-2\times10^{-6}\,M_{\odot}\,yr^{-1}$. We infer that the observed envelope and outflow structure around V883 Ori could be caused by multiple outbursts, consistent with episodic accretion.

We have conducted a 3-mm molecular-line survey towards 570 high-mass star-forming clumps, using the Mopra telescope. The sample is selected from the 10,000 clumps identified by the ATLASGAL survey and includes all of the most important embedded evolutionary stages associated with massive star formation, classified into five distinct categories (quiescent, protostellar, young stellar objects, \hii\ regions and photo-dominated regions). The observations were performed in broadband mode with frequency coverage of 85.2 to 93.4\,GHz and a velocity resolution of $\sim$0.9\,\kms, detecting emission from 26 different transitions. We find significant evolutionary trends in the detection rates, integrated line intensities, and abundances of many of the transitions and also identify a couple of molecules that appear to be invariant to changes in the dust temperature and evolutionary stage (N$_2$H$^+$\,(1-0) and HN$^{13}$C\,(1-0)). We use the K-ladders for CH$_3$C$_2$H\,(5-4) and CH$_3$CH\,(5-4) to calculate the rotation temperatures and find $\sim$1/3 of the quiescent clumps have rotation temperatures that suggest the presence of an internal heating source. These sources may constitute a population of very young protostellar objects that are still dark at 70\,\mum\ and suggest that the fraction of truly quiescent clumps may only be a few per cent. We also identify a number of line ratios that show a strong correlation with the evolutionary stage of the embedded objects and discuss their utility as diagnostic probes of evolution.

High-mass stars are formed within massive molecular clumps, where a large number of stars form close together. The evolution of the clumps with different masses and luminosities is mainly regulated by its high-mass stellar content and the formation of such objects is still not well understood. In this work, we characterise the mid-J CO emission in a statistical sample of 99 clumps (Top100) selected from the ATLASGAL survey that are representative of the Galactic proto-cluster population. High-spatial resolution APEX-CHAMP+ maps of the CO(6-5) and CO(7-6) transitions were obtained and combined with additional single-pointing APEX-FLASH+ spectra of the CO(4-3) line. We study the correlations of the CO line luminosities and profiles for the three CO transitions with the clump properties and investigate if and how they change as a function of the evolution. All sources were detected above 3-$\sigma$ in all three CO transitions and most of the sources exhibit broad CO emission likely associated with molecular outflows. We found that the extension of the mid-J CO emission is correlated with the size of the dust emission traced by the Herschel-PACS 70 $\mu$m maps. The CO line luminosity is correlated with the luminosity and mass of the clumps. However, it does not correlate with the L/M ratio. The dependency of the CO luminosity with the properties of the clumps is steeper for higher-J transitions. Our data seem to exclude that this trend is biased by self-absorption features in the CO emission, but rather suggest that different J transitions arise from different regions of the inner envelope. Moreover, high-mass clumps show similar trends in CO luminosity as lower mass clumps, but are systematically offset towards larger values, suggesting that higher column density and/or temperature (of unresolved) CO emitters are found inside high-mass clumps.

We present the performance of the upGREAT heterodyne array receivers on the SOFIA telescope after several years of operations. This instrument is a multi-pixel high resolution (R > 10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receivers use 7-pixel subarrays configured in a hexagonal layout around a central pixel. The low frequency array receiver (LFA) has 2x7 pixels (dual polarization), and presently covers the 1.83-2.06 THz frequency range, which allows to observe the [CII] and [OI] lines at 158 um and 145 um wavelengths. The high frequency array (HFA) covers the [OI] line at 63 um and is equipped with one polarization at the moment (7 pixels, which can be upgraded in the near future with a second polarization array). The 4.7 THz array has successfully flown using two separate quantum-cascade laser local oscillators from two different groups. NASA completed the development, integration and testing of a dual-channel closed-cycle cryocooler system, with two independently operable He compressors, aboard SOFIA in early 2017 and since then, both arrays can be operated in parallel using a frequency separating dichroic mirror. This configuration is now the prime GREAT configuration and has been added to SOFIA's instrument suite since observing cycle 6.

G351.776-0.527 is among the most massive, closest, and youngest filaments in the inner Galactic plane and therefore it is an ideal laboratory to study the kinematics of dense gas and mass replenishment on a large scale. In this paper, we present far-infrared (FIR) and submillimetre wavelength continuum observations combined with spectroscopic C$^{18}$O (2-1) data of the entire region to study its temperature, mass distribution, and kinematics. The structure is composed of a main elongated region with an aspect ratio of $\sim 23$, which is associated with a network of filamentary structures. The main filament has a remarkably constant width of 0.2 pc. The total mass of the network (including the main filament) is $\geq 2600$ M$_\odot$, while we estimate a mass of $\sim 2000$ M$_\odot$ for the main structure. Therefore, the network harbours a large reservoir of gas and dust that could still be accreted onto the main structure. From the analysis of the gas kinematics, we detect two velocity components in the northern part of the main filament. The data also reveal velocity oscillations in C$^{18}$O along the spine in the main filament and in at least one of the branches. Considering the region as a single structure, we find that it is globally close to virial equilibrium indicating that the entire structure is approximately in a stable state.

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.

Aims. To constrain the physical processes that lead to the birth of high-mass stars it is mandatory to study the very first stages of their formation. We search for high-mass analogs of low-mass prestellar cores in W43-MM1. Methods. We conducted a 1.3 mm ALMA mosaic of the complete W43-MM1 cloud, which has revealed numerous cores with ~ 2000 au FWHM sizes. We investigated the nature of cores located at the tip of the main filament, where the clustering is minimum. We used the continuum emission to measure the core masses and the $^{13}$CS(5-4) line emission to estimate their turbulence level. We also investigated the prestellar or protostellar nature of these cores by searching for outflow signatures traced by CO(2-1) and SiO(5-4) line emission, and for molecular complexity typical of embedded hot cores. Results. Two high-mass cores of ~ 1300 au diameter and ~ $60~M_\odot$ mass are observed to be turbulent but gravitationally bound. One drives outflows and is associated with a hot core. The other core, W43-MM1\#6, does not yet reveal any star formation activity and thus is an excellent high-mass prestellar core candidate.

Analysing the kinematics of filamentary molecular clouds is a crucial step towards understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of the Inner Galactic Inter Stellar Medium) allows us to analyse the kinematics of these targets and to determine their physical properties at a resolution of 30 arcsec and 0.25 km/s. To do so, we developed an automated algorithm to identify all velocity components along the line-of-sight correlated with the ATLASGAL dust emission, and derive size, mass, and kinematic properties for all velocity components. We find two-third of the filament candidates are coherent structures in position-position-velocity space. The remaining candidates appear to be the result of a superposition of two or three filamentary structures along the line-of-sight. At the resolution of the data, on average the filaments are in agreement with Plummer-like radial density profiles with a power-law exponent of p = 1.5 +- 0.5, indicating that they are typically embedded in a molecular cloud and do not have a well-defined outer radius. Also, we find a correlation between the observed mass per unit length and the velocity dispersion of the filament of $m \sim \sigma_v^2$. We show that this relation can be explained by a virial balance between self-gravity and pressure. Another possible explanation could be radial collapse of the filament, where we can exclude infall motions close to the free-fall velocity.

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.

The central area (40"x40") of the bipolar nebula S106 was mapped in the OI line at 63.2 micron with high angular (6") and spectral resolution, using GREAT on board SOFIA. The OI emission distribution is compared to the CO 16-15, CII 158 micron, and CO 11-10 lines, mm-molecular lines, and continuum. It is composed of several velocity components in the range from -30 km/s to 25 km/s. The high-velocity blue- and redshifted emission can be explained as arising from accelerated photodissociated (PDR) gas associated with a dark lane close to the massive binary system S106 IR, and from shocks caused by the stellar wind and/or a disk--envelope interaction. At velocities from -9 to -4 km/s and 0.5 to 8 km/s line wings are observed that we attribute to cooling in PDRs created by the ionizing radiation impinging on the cavity walls. The bulk velocity range is dominated by PDR emission from the clumpy molecular cloud. Modelling the emission in the different velocity ranges with the KOSMA-tau code constrains a radiation field chi of a few times 10^4 and densities n of a few times 10^4 cm^-3. Considering self-absorption of the OI line results in higher densities (up to 10^6 cm^-3) only for the gas component seen at high blue- and red velocities. The dark lane has a mass of 275 Msun and shows a velocity difference of 1.4 km/s along its projected length of 1 pc, determined from H13CO+ 1-0 mapping. It can be interpreted as a massive accretion flow, or the remains of it, linked to S106 IR/FIR. The most likely explanation is that the binary system is at a stage of its evolution where gas accretion is counteracted by the stellar winds and radiation, leading to the very complex observed spatial and kinematic emission distribution of the various tracers.

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.

Previous radio recombination line (RRL) observations of dust clumps identified in the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) have led to the detection of a large number of RRLs in the 3mm range. Here, we aim to study their excitation with shorter wavelength (sub)millimeter radio recombination line (submm-RRL) observations. We made observations of submm-RRLs with low principal quantum numbers ($n$ $\leq$ 30) using the APEX 12 m telescope, toward 104 HII regions associated with massive dust clumps from ATLASGAL. The observations covered the H25$\alpha$, H28$\alpha$, and H35$\beta$ transitions. Toward a small subsample the H26$\alpha$, H27$\alpha$, H29$\alpha$, and H30$\alpha$ lines were observed to avoid contamination by molecular lines at adjacent frequencies. We have detected submm-RRLs (signal-to-noise $\geq$ 3 $\sigma$) from compact HII regions embedded within 93 clumps. The submm-RRLs are approximately a factor of two brighter than the mm-RRLs and consistent with optically thin emission in local thermodynamic equilibrium (LTE). The average ratio (0.31) of the measured H35$\beta$/H28$\alpha$ fluxes is close to the LTE value of 0.28. No indication of RRL maser emission has been found. The Lyman photon flux, bolometric, and submm-RRL luminosities toward the submm-RRL detected sources present significant correlations. The trends of dust temperature and the ratio of bolometric luminosity to clump mass, $L_{\rm bol}/M_{\rm clump}$, indicate that the HII regions are related to the most massive and luminous clumps. By estimating the production rate of ionizing photons, $Q$, from the submm-RRL flux, we find that the $Q$(H28$\alpha$) measurements provide estimates of the Lyman continuum photon flux consistent with those determined from 5 GHz radio continuum emission.

The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a physical scale of $\sim$400 au in dust continuum and molecular lines. The dust continuum emission reveals a single high-mass protostellar envelope and shows evidence for a marginally resolved continuum source. We detect a rotational line of CH$_3$OH within its $v_{\rm t}$=1 torsionally excited state revealing two bright peaks of emission spatially offset from the dust continuum peak, and exhibiting a distinct velocity component $\pm$4.5 km s$^{-1}$ offset compared to the source $v_{\rm lsr}$. Local thermodynamic equilibrium analysis suggests N(CH$_3$OH)=1.2$-$2$\times$10$^{19}$ cm$^{-2}$, and kinetic temperatures of 160$-$200 K at the position of these peaks. Their velocity shifts correspond well to the expected Keplerian velocity around a central object with 15M$_{\odot}$ consistent with the mass estimate based on the source's bolometric luminosity. We propose a picture where the CH$_3$OH emission peaks trace the accretion shocks around the centrifugal barrier, pinpointing the interaction region between the collapsing envelope and an accretion disk. Because the HC$_3$N $v_{\rm 7}$=1e ($J$=38-37) line shows compact emission, and a velocity pattern consistent with models of Keplerian rotation, we suggest that this could be a new tracer for compact accretion disks around high-mass protostars. The estimated physical properties of the accretion disk suggest a specific angular momentum several times larger than typically observed towards low-mass protostars.

The methylidyne radical CH is commonly used as a proxy for H$_2$ in the cold, neutral phase of the interstellar medium. The optical spectroscopy of CH is limited by interstellar extinction, whereas far-infrared observations provide an integral view through the Galaxy. While the HF ground state absorption, another H$_2$ proxy in diffuse gas, frequently suffers from saturation, CH remains transparent both in spiral-arm crossings and high-mass star forming regions, turning this light hydride into a universal surrogate for H$_2$. However, in slow shocks and in regions dissipating turbulence its abundance is expected to be enhanced by an endothermic production path, and the idea of a "canonical" CH abundance needs to be addressed. The $N=2\leftarrow 1$ ground state transition of CH at $\lambda 149\,\mu$m has become accessible to high-resolution spectroscopy thanks to GREAT aboard SOFIA. Its unsaturated absorption and the absence of emission makes it an ideal candidate for the determination of column densities with a minimum of assumptions. Here we present an analysis of four sightlines towards distant, far-infrared bright Galactic star forming regions. If combined with the sub-millimeter line of CH at $\lambda 560\,\mu$m, environments forming massive stars can be analyzed. For this we present a case study on the "proto-Trapezium" cluster W3 IRS5, and demonstrate that the sub-millimeter/far-infrared lines of CH reliably trace not only diffuse but also dense, molecular gas. While we confirm the global correlation between the column densities of HF and those of CH, clear signposts of an over-abundance of CH are observed towards lower densities. A quiescent ion-neutral chemistry alone cannot account for this over-abundance. Vortices forming in turbulent, diffuse gas may be the setting for an enhanced production path.

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC and its immediate surroundings. This strongest emission component has a velocity similar to that of the $^{13}$CO(2-1) emission of the IRDC, but offset by $\sim3\:{\rm km\:s}^{-1}$ and with a larger velocity width of $\sim9\:{\rm km\:s}^{-1}$. The spatial distribution of the [CII] emission of this component is also offset predominantly to one side of the dense filamentary structure of the IRDC. The CII column density is estimated to be of the order of $\sim10^{17}-10^{18}\,{\rm cm}^{-2}$. We compare these results to the [CII] emission from numerical simulations of magnetized, dense gas filaments formed from giant molecular cloud (GMC) collisions, finding similar spatial and kinematic offsets. These observations and modeling of [CII] add further to the evidence that IRDC G035.39-00.33 has been formed by a process of GMC-GMC collision, which may thus be an important mechanism for initiating star cluster formation.

We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $\mu$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we derive the gas kinetic temperature and spatial density using the measured p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$, 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$, and 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ ratios. The gas kinetic temperatures derived from the p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$ and 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$ line ratios are high, ranging from 43 to $>$300 K with an unweighted average of 91 $\pm$ 4 K. Deduced $T_{\rm kin}$ values from the $J$ = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the p-H$_2$CO 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ line ratios yield 0.6-8.3 $\times$ 10$^6$ cm$^{-3}$ with an unweighted average of 1.5 ($\pm$0.1) $\times$ 10$^6$ cm$^{-3}$. A comparison of kinetic temperatures derived from p-H$_2$CO, NH$_3$, and the dust emission indicates that p-H$_2$CO traces a distinctly higher temperature than the NH$_3$ (2,2)/(1,1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H$_2$CO linewidths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H$_2$CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H$_2$CO increase with time through the evolution of the clumps.

According to the current paradigm of circumstellar disk evolution, gas-rich primordial disks evolve into gas-poor debris disks compose of second-generation dust. To explore the transition between these phases, we searched for $^{12}$CO, $^{13}$CO, and C$^{18}$O emission in seven dust-rich debris disks around young A-type stars, using ALMA in Band 6. We discovered molecular gas in three debris disks. In all these disks, the $^{12}$CO line was optically thick, highlighting the importance of less abundant molecules in reliable mass estimates. Supplementing our target list by literature data, we compiled a volume-limited sample of dust-rich debris disks around young A-type stars within 150 pc. We obtained a CO detection rate of 11/16 above a $^{12}$CO J=2$-$1 line luminosity threshold of $\sim 1.4 \times 10 ^4$ Jykms$^{-1}$pc$^2$ in the sample. This high incidence implies that the presence of CO gas in bright debris disks around young A-type stars is likely more the rule than the exception. Interestingly, dust-rich debris disks around young FG-type stars exhibit, with the same detectability threshold as for A-type stars, significantly lower gas incidence. While the transition from protoplanetary to debris phase is associated with a drop of dust content, our results exhibit a large spread in the CO mass in our debris sample, with peak values comparable to those in protoplanetary Herbig Ae disks. In the particularly CO-rich debris systems the gas may have primordial origin, characteristic of a hybrid disk.

Abridged: ATLASGAL is an unbiased 870 micron submillimetre survey of the inner Galactic plane. It provides a large and systematic inventory of all massive, dense clumps in the Galaxy (>1000 Msun) and includes representative samples of all embedded stages of high-mass star formation. Here we present the first detailed census of the properties (velocities, distances, luminosities and masses) and spatial distribution of a complete sample of ~8000 dense clumps located in the Galactic disk. We derive highly reliable velocities and distances to ~97% of the sample and use mid- and far-infrared survey data to develop an evolutionary classification scheme that we apply to the whole sample. Comparing the evolutionary subsamples reveals trends for increasing dust temperatures, luminosities and line-widths as a function of evolution indicating that the feedback from the embedded proto-clusters is having a significant impact on the structure and dynamics of their natal clumps. We find 88\,per\,cent are already associated with star formation at some level. We also find the clump mass to be independent of evolution suggesting that the clumps form with the majority of their mass in-situ. We estimate the statistical lifetime of the quiescent stage to be ~5 x 10^4 yr for clump masses ~1000 Msun decreasing to ~1 x 10^4 yr for clump masses >10000 Msun. We find a strong correlation between the fraction of clumps associated with massive stars and peak column density. The fraction is initially small at low column densities but reaching 100\,per\,cent for column densities above 10^{23} cm^{-2}; there are no clumps with column density clumps above this value that are not already associated with massive star formation. All of the evidence is consistent with a dynamic view of star formation wherein the clumps form rapidly and are initially very unstable so that star formation quickly ensues.

The initial conditions of molecular clumps in which high-mass stars form are poorly understood. In particular, a more detailed study of the earliest evolutionary phases is needed. The APEX Telescope Large Area Survey of the whole inner Galactic disk at 870 micron, ATLASGAL, has been conducted to discover high-mass star-forming regions at different evolutionary phases. Using the Parkes telescope, we observed the NH3 (1,1) to (3,3) inversion transitions towards 354 ATLASGAL clumps in the fourth quadrant. For a subsample of 289 sources, the N2H+ (1-0) line was measured with the Mopra telescope. We measured a median NH3(1,1) line width of about 2 km/s and rotational temperatures from 12 to 28 K with a mean of 18 K. For a subsample with detected NH3 (2,2) hyperfine components, we found that the commonly used method to compute the (2,2) optical depth from the (1,1) optical depth and the (2,2) to (1,1) main beam brightness temperature ratio leads to an underestimation of the rotational temperature and column density. A larger median virial parameter of about 1 is determined using the broader N2H+ line width than is estimated from the NH3 line width of about 0.5 with a general trend of a decreasing virial parameter with increasing gas mass. We found a warmer surrounding of ATLASGAL clumps than the surrounding of low-mass cores and smaller velocity dispersions in low-mass than high-mass star-forming regions. The NH3 (1,1) inversion transition of 49% of the sources shows hyperfine structure anomalies. The intensity ratio of the outer hyperfine structure lines with a median of 1.27+/-0.03 and a standard deviation of 0.45 is significantly higher than 1, while the intensity ratios of the inner satellites with a median of 0.9+/-0.02 and standard deviation of 0.3 and the sum of the inner and outer hyperfine components with a median of 1.06+/-0.02 and standard deviation of 0.37 are closer to 1.

A long-standing open issue of the paradigm of low-mass star formation is the luminosity problem: most protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Ori-type stars (FUors) are thought to be the visible examples for objects in the high accretion state. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. However, we have insufficient information on the envelope dynamics in FUors, about where and how mass transfer from the envelope to the disk happens. Here we present ALMA observations of the FUor-type star V346 Nor at 1.3 mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyze the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We found that V346 Nor is surrounded by gaseous material on 10000 au scale in which a prominent outflow cavity is carved. Within the central $\sim$700 au, the circumstellar matter forms a flattened pseudo-disk where material is infalling with conserved angular momentum. Within $\sim$350 au, the velocity profile is more consistent with a disk in Keplerian rotation around a central star of 0.1 $M_{\odot}$. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$yr$^{-1}$, a factor of few higher than the quiescent accretion rate from the disk onto the star, hinting for a mismatch between the infall and accretion rates as the cause of the eruption.

(Abridged) Aims: We aim to use the progressive heating of the gas caused by the feedback of high-mass young stellar objects (YSOs) to prove the statistical validity of the most common schemes used to define an evolutionary sequence for high-mass clumps, and characterise the sensitivity of different tracers to this process. Methods: From the spectroscopic follow-ups of the ATLASGAL TOP100 sample, we selected several multiplets of CH3CN, CH3CCH, and CH3OH emission lines to derive and compare the physical properties of the gas in the clumps along the evolutionary sequence. Our findings are compared with results obtained from CO isotopologues, dust, and NH3 from previous studies on the same sample. Results: The chemical properties of each species have a major role on the measured physical properties. Low temperatures are traced by NH3, CH3OH, and CO (in the early phases), the warm and dense envelope can be probed with CH3CN, CH3CCH, and, in evolved sources via CO isotopologues. CH3OH and CH3CN are also abundant in the hot cores, and their high-excitation transitions may be good tools to study the kinematics in the hot gas surrounding the YSOs that these clumps are hosting. All tracers show, to different degrees, progressive warming with evolution. The relation between gas temperature and L/M is reproduced by a toy model of a spherical, internally heated clump. Conclusions: The evolutionary sequence defined for the clumps is statistically valid and we could identify the processes dominating in different intervals of L/M. For L/M<2Lsun/Msun a large quantity of gas is still being accumulated and compressed at the bottom of the potential well. Between 2Lsun/Msun<L/M<40Lsun/Msun the YSOs gain mass and increase in L; the first hot cores appear around L/M=10Lsun/Msun. Finally, for L/M>40Lsun/Msun HII regions become common, showing that dissipation of the parental clump dominates.

The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales.

Observations of millimeter wavelength radio recombination lines (mm-RRLs) are used to search for HII regions in an unbiased way that is complementary to many of the more traditional methods previously used. We carried out targeted mm-RRL observations (39$\leq$ principal quantum number ($n$) $\leq$65 and $\Delta n$ = 1, 2, 3 and 4, i.e. H$n\alpha$, H$n\beta$, H$n\gamma$ and H$n\delta$) using the IRAM 30m and Mopra 22m telescopes. We observed 976 compact dust clumps selected from a catalog of $\sim$10,000 ATLASGAL sources. We detected H$n\alpha$ mm-RRL emission toward 178 clumps; H$n\beta$, H$n\gamma$ and H$n\delta$ were also detected toward 65, 23, and 22 clumps, respectively. This is the largest sample of mm-RRLs detections published to date. Comparing the positions of these clumps with radio continuum surveys we identified compact radio counterparts for 134 clumps. The nature of the other 44 detections is unclear, but 8 detections are thought to be potentially new HII regions. Broad linewidths are seen toward nine clumps (linewidth > 40 km s$^{-1}$) revealing significant turbulent motions within the ionized gas. We find that the linewidth of the H$^{13}$CO$^{+}$ (1$-$0) emission is significantly wider than those without mm-RRL detection. We also find a correlation between the integrated fluxes of the mm-RRLs and the 6 cm continuum flux densities of their radio counterparts. We find that the mm-RRL emission is associated with HII regions with $n_{e}<10^{5}$ cm$^{-3}$ and HII region diameter >0.03 pc. We detected mm-RRLs toward 178 clumps and identified eight new HII region candidates. The broad mm-RRL from nine clumps may indicate that they arise in very young hyper-compact HII regions. The mm-RRLs trace the radio continuum sources detected by high-resolution observations and their line parameters show associations with the embedded radio sources and their parental molecular clumps.

FUors are young stellar objects experiencing large optical outbursts due to highly enhanced accretion from the circumstellar disk onto the star. FUors are often surrounded by massive envelopes, which play a significant role in the outburst mechanism. Conversely, the subsequent eruptions might gradually clear up the obscuring envelope material and drive the protostar on its way to become a disk-only T Tauri star. Here we present an APEX $^{12}$CO and $^{13}$CO survey of eight southern and equatorial FUors. We measure the mass of the gaseous material surrounding our targets. We locate the source of the CO emission and derive physical parameters for the envelopes and outflows, where detected. Our results support the evolutionary scenario where FUors represent a transition phase from envelope-surrounded protostars to classical T Tauri stars.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

The progenitors of high-mass stars and clusters are still challenging to recognise. Only unbiased surveys, sensitive to compact regions of high dust column density, can unambiguously reveal such a small population of particularly massive and cold clumps. Here we study a flux limited sample of compact sources from the ATLASGAL survey to identify a sample of candidate progenitors of massive clusters in the inner Galaxy. Sensitive mid-infrared data at 21-24 $\mu$m from the WISE and MIPSGAL surveys were explored to search for embedded objects, and complementary spectroscopic data were used to investigate their stability and star formation activity. Based on such ancillary data we identify an unbiased sample of infrared-quiet massive clumps in the Galaxy that potentially represent the earliest stages of massive cluster formation. An important fraction of this sample consists of sources that have not been studied in detail before. Comparing their properties to clumps hosting more evolved embedded objects, we find that they exhibit similar physical properties in terms of mass and size, suggesting that infrared-quiet massive clumps are not only capable of forming high-mass stars, but likely also follow a single evolutionary track leading to the formation of massive clusters. The majority of the sources are not in virial-equilibrium, suggesting collapse on the clump scale. This is in line with the low number of infrared-quiet massive clumps and earlier findings that star formation, in particular for high-mass objects is a fast, dynamic process. We propose a scenario in which massive clumps start to fragment and collapse before their final mass is accumulated indicating that strong self-gravity and global collapse is needed to build up rich clusters and the most massive stars.

The [CII] 158 um fine structure line is one of the dominant cooling lines in the interstellar medium (ISM) and is an important tracer of star formation. Recent velocity-resolved studies with Herschel/HIFI and SOFIA/GREAT showed that the [CII] line can constrain the properties of the ISM phases in star-forming regions. The [CII] line as a tracer of star formation is particularly important in low-metallicity environments where CO emission is weak because of the presence of large amounts of CO-dark gas. The nearby irregular dwarf galaxy NGC 4214 offers an excellent opportunity to study an actively star-forming ISM at low metallicity. We analyzed the spectrally resolved [CII] line profiles in three distinct regions at different evolutionary stages of NGC 4214 with respect to ancillary HI and CO data in order to study the origin of the [CII] line. We used SOFIA/GREAT [CII] 158 um observations, HI data from THINGS, and CO(2-1) data from HERACLES to decompose the spectrally resolved [CII] line profiles into components associated with neutral atomic and molecular gas. We use this decomposition to infer gas masses traced by [CII] under different ISM conditions. Averaged over all regions, we associate about 46% of the [CII] emission with the HI emission. However, we can assign only around 9% of the total [CII] emission to the cold neutral medium (CNM). We found that about 79% of the total molecular hydrogen mass is not traced by CO emission. On average, the fraction of CO-dark gas dominates the molecular gas mass budget. The fraction seems to depend on the evolutionary stage of the regions: it is highest in the region covering a super star cluster in NGC 4214, while it is lower in a more compact, more metal-rich region.

The ATLASGAL survey provides an ideal basis for detailed studies of large numbers of massive star forming clumps covering the whole range of evolutionary stages. The ATLASGAL Top100 is a sample of clumps selected from their infrared and radio properties to be representative for the whole range of evolutionary stages. The ATLASGAL Top100 sources are the focus of a number of detailed follow-up studies that will be presented in a series of papers. In the present work we use the dust continuum emission to constrain the physical properties of this sample and identify trends as a function of source evolution. We determine flux densities from mid-infrared to submm wavelength (8-870 micron) images and use these values to fit their spectral energy distributions (SEDs) and determine their dust temperature and flux. Combining these with recent distances from the literature including maser parallax measurements we determine clump masses, luminosities and column densities. We find trends for increasing temperature, luminosity and column density with the proposed evolution sequence, confirming that this sample is representative of different evolutionary stages of massive star formation. We show that most of the sample has the ability to form massive stars (including the most massive O-type stars) and that the majority is gravitationally unstable and hence likely to be collapsing. The highest column density ATLASGAL sources presented cover the whole range of evolutionary stages from the youngest to the most evolved high-mass star forming clumps. Their study provides a unique starting point for more in-depth research on massive star formation in four distinct evolutionary stages whose well defined physical parameters afford more detailed studies. As most of the sample is closer than 5 kpc, these sources are also ideal for follow-up observations with high spatial resolution.

We present a new multi-pixel high resolution (R >10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receiver uses 2 x 7-pixel subarrays in orthogonal polarization, each in an hexagonal array around a central pixel. We present the first results for this new instrument after commissioning campaigns in May and December 2015 and after science observations performed in May 2016 . The receiver is designed to ultimately cover the full 1.8-2.5 THz frequency range but in its first implementation, the observing range was limited to observations of the [CII] line at 1.9 THz in 2015 and extended to 1.83-2.07 THz in 2016. The instrument sensitivities are state-of-the-art and the first scientific observations performed shortly after the commissioning confirm that the time efficiency for large scale imaging is improved by more than an order of magnitude as compared to single pixel receivers. An example of large scale mapping around the Horsehead Nebula is presented here illustrating this improvement. The array has been added to SOFIA's instrument suite already for ongoing observing cycle 4.

A significant fraction of nearby young moving group members harbor circumstellar debris dust disks. Due to their proximity and youth, these disks are attractive targets for studying the early evolution of debris dust and planetesimal belts. Here we present 70 and 160$\mu$m observations of 31 systems in the $\beta$ Pic moving group, and in the Tucana-Horologium, Columba, Carina and Argus associations, using the Herschel Space Observatory. None of these stars were observed at far-infrared wavelengths before. Our Herschel measurements were complemented by photometry from the WISE satellite for the whole sample, and by submillimeter/millimeter continuum data for one source, HD 48370. We identified six stars with infrared excess, four of them are new discoveries. By combining our new findings with results from the literature, we examined the incidence and general characteristics of debris disks around Sun-like members of the selected groups. With their dust temperatures of <45 K the newly identified disks around HD 38397, HD 48370, HD 160305, and BD-20 951 represent the coldest population within this sample. For HD 38397 and HD 48370, the emission is resolved in the 70$\mu$m PACS images, the estimated radius of these disks is ~90 au. Together with the well-known disk around HD 61005, these three systems represent the highest mass end of the known debris disk population around young G-type members of the selected groups. In terms of dust content, they resemble the hypothesized debris disk of the ancient Solar System.

[Abridged] Aims. We study the properties of filamentary structures from the ATLASGAL survey. Methods. We use the DisPerSE algorithm to identify spatially coherent structures located across the inner-Galaxy (300 < l < 60 and |b| < 1.5). Results. We have determined distances, masses and physical sizes for 241 of the filamentary structures. We find a median distance of 3.8 kpc, a mean mass of a few 10^3 m_sun, a mean length of ~6pc and a mass-to-length ratio of (M/L) ~200-2000M_sun/ pc. We also find that these filamentary structures are tightly correlated with the spiral arms in longitude and velocity, and that their semi-major axis is preferentially aligned parallel to the Galactic mid-plane and therefore with the direction of large-scale Galactic magnetic field. We find many examples where the dense filaments identified in ATLASGAL are associated with larger scale filamentary structures (~100 pc), and argue that this is likely to be common, and as such these may indicate a connection between large-scale Galactic dynamics and star formation. Conclusions. We have produced a large and Galaxy-wide catalogue of dense filamentary structures that are representative of a particular size and mass range not previously well studied in the literature. Analyses of the properties and distribution of these filaments reveals that they are correlated with the spiral arms and make a significant contribution to star formation in the Galaxy. Massive star formation is ongoing within ~20% of the filaments and is strongly correlated with the filaments with the largest mass-to- length ratios. The luminosity of the embedded sources has a similar distribution to the Galactic-wide samples of young massive stars and can therefore be considered to be representative.

Probability distribution functions (PDFs) of column densities are an established tool to characterize the evolutionary state of interstellar clouds. Using simulations, we show to what degree their determination is affected by noise, line-of-sight contamination, field selection, and the incomplete sampling in interferometric measurements. We solve the integrals that describe the convolution of a cloud PDF with contaminating sources and study the impact of missing information on the measured column density PDF. The effect of observational noise can be easily estimated and corrected for if the root mean square (rms) of the noise is known. For $\sigma_{noise}$ values below 40% of the typical cloud column density, $N_{peak}$, this involves almost no degradation of the accuracy of the PDF parameters. For higher noise levels and narrow cloud PDFs the width of the PDF becomes increasingly uncertain. A contamination by turbulent foreground or background clouds can be removed as a constant shield if the PDF of the contamination peaks at a lower column or is narrower than that of the observed cloud. Uncertainties in the definition of the cloud boundaries mainly affect the low-column density part of the PDF and the mean density. As long as more than 50% of a cloud are covered, the impact on the PDF parameters is negligible. In contrast, the incomplete sampling of the uv plane in interferometric observations leads to uncorrectable distortions of the PDF of the produced maps. An extension of ALMA's capabilities would allow us to recover the high-column density tail of the PDF but we found no way to measure the intermediate and low column density part of the underlying cloud PDF in interferometric observations.

EX Lupi-type objects (EXors) form a sub-class of T Tauri stars, defined by sudden sporadic flare-ups of 1-5 magnitudes at optical wavelengths. These eruptions are attributed to enhanced mass accretion from the circumstellar disk to the star, and may constitute important events in shaping the structure of the inner disk and the forming planetary system. Although disk properties must play a fundamental role in driving the outbursts, they are surprisingly poorly known. In order to characterize the dust and gas components of EXor disks, here we report on observations of the $^{12}$CO J=3-2 and 4-3 lines, and the $^{13}$CO 3-2 line in EX Lup, the prototype of the EXor class. We reproduce the observed line fluxes and profiles with a line radiative transfer model, and compare the obtained parameters with corresponding ones of other T Tauri disks.

We present a statistical study of the distribution and physical properties of cold dense material in and around the inner Galactic Plane near infrared bubbles as catalogued by the Milky Way Project citizen scientists. Using data from the ATLASGAL 870 um survey, we show that 48 +/- 2% of all cold clumps in the studied survey region (|l| <= 65 degrees, |b| <= 1 degree) are found in close proximity to a bubble, and 25 +/- 2% appear directly projected towards a bubble rim. A two-point correlation analysis confirms the strong correlation of massive cold clumps with expanding bubbles. It shows an overdensity of clumps along bubble rims that grows with increasing bubble size, which shows how interstellar medium material is reordered on large scales by bubble expansion around regions of massive star formation. The highest column density clumps appear resistent to the expansion, remaining overdense towards the bubbles' interior rather than being swept up by the expanding edge. Spectroscopic observations in ammonia show that cold dust clumps near bubbles appear to be denser, hotter and more turbulent than those in the field, offering circumstantial evidence that bubble-associated clumps are more likely to be forming stars. These observed differences in physical conditions persist for beyond the region of the bubble rims.

We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed.

Newborn stars form within the localized, high density regions of molecular clouds. The sequence and rate at which stars form in dense clumps and the dependence on local and global environments are key factors in developing descriptions of stellar production in galaxies. We seek to observationally constrain the rate and latency of star formation in dense massive clumps that are distributed throughout the Galaxy and to compare these results to proposed prescriptions for stellar production. A sample of 24 micron-based Class~I protostars are linked to dust clumps that are embedded within molecular clouds selected from the APEX Telescope Large Area Survey of the Galaxy. We determine the fraction of star-forming clumps, f*, that imposes a constraint on the latency of star formation in units of a clump's lifetime. Protostellar masses are estimated from models of circumstellar environments of young stellar objects from which star formation rates are derived. Physical properties of the clumps are calculated from 870 micron dust continuum emission and NH_3 line emission. Linear correlations are identified between the star formation rate surface density, Sigma_{SFR}, and the quantities Sigma_{H2}/tau_{ff} and Sigma_{H2}/tau_{cross}, suggesting that star formation is regulated at the local scales of molecular clouds. The measured fraction of star forming clumps is 23%. Accounting for star formation within clumps that are excluded from our sample due to 24 micron saturation, this fraction can be as high as 31%. Dense, massive clumps form primarily low mass (< 1-2 msun) stars with emergent 24 micron fluxes below our sensitivity limit or are incapable of forming any stars for the initial 70% of their lifetimes. The low fraction of star forming clumps in the Galactic center relative to those located in the disk of the Milky Way is verified.

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

The origin of massive stars is a fundamental open issue in modern astrophysics. Pre-ALMA interferometric studies reveal precursors to early B to late O type stars with collapsing envelopes of 15-20 M$_\odot$ on 1000-3000 AU size-scales. To search for more massive envelopes we selected the most massive nearby young clumps from the ATLASGAL survey to study their protostellar content with ALMA. Our first results using the intermediate scales revealed by the ALMA ACA array providing 3-5" angular resolution, corresponding to $\sim$0.05-0.1 pc size-scales, reveals a sample of compact objects. These massive dense cores are on average two-times more massive than previous studies of similar types of objects. We expect that once the full survey is completed, it will provide a comprehensive view on the origin of the most massive stars.

Sensitive ground-based submillimeter surveys, such as ATLASGAL, provide a global view on the distribution of cold dense gas in the Galactic plane. Here we use the 353 GHz maps from the Planck/HFI instrument to complement the ground-based APEX/LABOCA observations with information on larger angular scales. The resulting maps reveal the distribution of cold dust in the inner Galaxy with a larger spatial dynamic range. We find examples of elongated structures extending over angular scales of 0.5 degree. Corresponding to >30 pc structures in projection at a distance of 3 kpc, these dust lanes are very extended and show large aspect ratios. Furthermore, we assess the fraction of dense gas ($f_{\rm DG}$), and estimate 2-5% (above A$_{\rm{v}}>$7 mag) on average in the Galactic plane. PDFs of the column density reveal the typically observed log-normal distribution for low- and exhibit an excess at high column densities. As a reference for extragalactic studies, we show the line-of-sight integrated N-PDF of the inner Galaxy, and derive a contribution of this excess to the total column density of $\sim2.2$%, above $N_{\rm H_2} = 2.92\times10^{22}$ cm$^{-2}$. Taking the total flux density, we provide an independent estimate of the mass of molecular gas in the inner Galaxy of $\sim1\times10^9\,M_{\odot}$, which is consistent with previous estimates using CO emission. From the mass and $f_{\rm DG}$ we estimate a Galactic SFR of $\dot M = 1.3\,M_{\odot}$ yr$^{-1}$. While the distribution of diffuse gas is homogenous in the inner Galaxy, the CMZ stands out with a higher dense gas fraction. The low star formation efficiency of the Milky Way is well explained by the low $f_{\rm DG}$ in the Galactic ISM, while the high $f_{\rm DG}$ towards the CMZ, despite its low star formation activity, suggests that, in that particular region of our Galaxy, high-density gas is not the bottleneck for star formation.

The processes leading to the birth of high-mass stars are poorly understood. We characterise here a sample of 430 massive clumps from the ATLASGAL survey, which are representative of different evolutionary stages. To establish a census of molecular tracers of their evolution we performed an unbiased spectral line survey covering the 3-mm atmospheric window between 84-117 GHz with the IRAM 30m. A smaller sample of 128 clumps has been observed in the SiO (5-4) transition with the APEX telescope to complement the SiO (2-1) line and probe the excitation conditions of the emitting gas, which is the main focus of the current study. We report a high detection rate of >75% of the SiO (2-1) line and a >90% detection rate from the dedicated follow-ups in the (5-4) transition. The SiO (2-1) line with broad line profiles and high detection rates, is a powerful probe of star formation activity, while the ubiquitous detection of SiO in all evolutionary stages suggests a continuous star formation process in massive clumps. We find a large fraction of infrared-quiet clumps to exhibit SiO emission, the majority of them only showing a low-velocity component (FWHM~5-6 km/s) centred at the rest velocity of the clump. In the current picture, where this is attributed to low-velocity shocks from cloud-cloud collisions, this can be used to pinpoint the youngest, thus, likely prestellar massive structures. Based on the line ratio of the (5-4) to the (2-1) line, our study reveals a trend of changing excitation conditions that lead to brighter emission in the (5-4) line towards more evolved sources. Our analysis delivers a more robust estimate of SiO column density and abundance than previous studies and questions the decrease of jet activity in massive clumps as a function of age.

With the GREAT receiver at the Stratospheric Observatory for Infrared Astronomy (SOFIA), nine massive molecular clumps have been observed in the ammonia $3_{2+}- 2_{2-}$ line at 1.8~THz in a search for signatures of infall. The sources were selected from the ATLASGAL submillimeter dust continuum survey of our Galaxy. Clumps with high masses covering a range of evolutionary stages based on their infrared properties were chosen. The ammonia line was detected in all sources, leading to five new detections and one confirmation of a previous detection of redshifted absorption in front of their strong THz continuum as a probe of infall in the clumps. These detections include two clumps embedded in infrared dark clouds. The measured velocity shifts of the absorptions compared to optically thin \CSEO\ (3--2) emission are 0.3--2.8~km/s, corresponding to fractions of 3\%\ to 30\% of the free-fall velocities of the clumps. The ammonia infall signature is compared with complementary data of different transitions of HCN, HNC, CS, and HCO$^+$, which are often used to probe infall via their blue-skewed line profiles. The best agreement with the ammonia results is found for the HCO$^+$ (4--3) transitions, but the latter is still strongly blended with emission from associated outflows. This outflow signature is far less prominent in the THz ammonia lines, which confirms it as a powerful probe of infall in molecular clumps. Infall rates in the range from 0.3 to 16~$10^{-3}\,M_\odot/$yr were derived with a tentative correlation with the virial parameters of the clumps. The new observations show that infall on clump scales is ubiquitous through a wide range of evolutionary stages, from $L/M$ covering about ten to several hundreds.

Debris disks are considered to be gas-poor, but recent observations revealed molecular or atomic gas in several 10-40 Myr old systems. We used the APEX and IRAM 30m radiotelescopes to search for CO gas in 20 bright debris disks. In one case, around the 16 Myr old A-type star HD 131835, we discovered a new gas-bearing debris disk, where the CO 3-2 transition was successfully detected. No other individual system exhibited a measurable CO signal. Our Herschel Space Observatory far-infrared images of HD 131835 marginally resolved the disk both at 70 and 100$\mu$m, with a characteristic radius of ~170 au. While in stellar properties HD 131835 resembles $\beta$ Pic, its dust disk properties are similar to those of the most massive young debris disks. With the detection of gas in HD 131835 the number of known debris disks with CO content has increased to four, all of them encircling young ($\leq$40 Myr) A-type stars. Based on statistics within 125 pc, we suggest that the presence of detectable amount of gas in the most massive debris disks around young A-type stars is a common phenomenon. Our current data cannot conclude on the origin of gas in HD 131835. If the gas is secondary, arising from the disruption of planetesimals, then HD 131835 is a comparably young and in terms of its disk more massive analogue of the $\beta$ Pic system. However, it is also possible that this system similarly to HD 21997 possesses a hybrid disk, where the gas material is predominantly primordial, while the dust grains are mostly derived from planetesimals.

OH is a key molecule in H2O chemistry, a valuable tool for probing physical conditions, and an important contributor to the cooling of shock regions. OH participates in the re-distribution of energy from the protostar towards the surrounding ISM. Our aim is to assess the origin of the OH emission from the Cepheus A massive star-forming region and to constrain the physical conditions prevailing in the emitting gas. We thus want to probe the processes at work during the formation of massive stars. We present spectrally resolved observations of OH towards the outflows of Cepheus A with the GREAT spectrometer onboard the SOFIA telescope. Three triplets were observed at 1834.7 GHz, 1837.8 GHz, and 2514.3 GHz (163.4, 163.1, and 119.2 microns), at angular resolutions of 16.3", 16.3", and 11.9", respectively. We present the CO (16-15) spectrum at the same position. We compared the integrated intensities in the redshifted wings to shock models. The two triplets near 163 microns are detected in emission with blending hyperfine structure unresolved. Their profiles and that of CO can be fitted by a combination of 2 or 3 Gaussians. The observed 119.2 microns triplet is seen in absorption, since its blending hyperfine structure is unresolved, but with three line-of-sight components and a blueshifted emission wing consistent with that of the other lines. The OH line wings are similar to those of CO, suggesting that they emanate from the same shocked structure. Under this common origin assumption, the observations fall within the model predictions and within the range of use of our model only if we consider that four shock structures are caught in our beam. Our comparisons suggest that the observations might be consistently fitted by a J-type model with nH > 1e5 cm-3, v > 20 km/s, and with a filling factor of ~1. Such a high density is generally found in shocks associated to high-mass protostars.

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

Context: Infrared dark clouds are the coldest and densest portions of giant molecular clouds. The most massive ones represent some of the most likely birthplaces for the next generation of massive stars in the Milky Way. Because a strong mid-IR background is needed to make them appear in absorption, they are usually assumed to be nearby. Aims: We use THz absorption spectroscopy to solve the distance ambiguity associated with kinematic distances for the IR-dark clouds in the TOP100 ATLASGAL sample, a flux-limited selection of massive clumps in different evolutionary phases of star formation. Methods: The para-H2O ground state transition at 1113.343 GHz, observed with Herschel/HIFI, was used to investigate the occurrence of foreground absorption along the line of sight directly towards infrared-dark clouds. Additional consistency checks were performed using MALT90 and HiGAL archival data and targeted Mopra and APEX spectroscopic observations. Results: We report the first discovery of five IRDCs in the TOP100 lying conclusively at the far kinematic distance, showing that the mere presence of low-contrast mid-IR absorption is not sufficient to unequivocally resolve the near/far ambiguity in favour of the former. All IRDCs are massive and actively forming high-mass stars; four of them also show infall signatures. Conclusions: We give a first estimate of the fraction of dark sources at the far distance (~11% in the TOP100) and describe their appearance and properties. The assumption that all dark clouds lie at the near distance may lead, in some cases, to underestimating masses, sizes, and luminosities, possibly causing clouds to be missed that will form very massive stars and clusters.

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

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

The formation of high mass stars and clusters occurs in giant molecular clouds. Objects in evolved stages of massive star formation such as protostars, hot molecular cores, and ultracompact HII regions have been studied in more detail than earlier, colder objects. With this in mind, the APEX Telescope Large Area Survey of the whole inner Galactic plane at 870 micron (ATLASGAL) has been carried out to provide a global view of cold dust and star formation at submillimetre wavelengths. To derive kinematic distances to a large sample of ATLASGAL clumps we divided them into groups of sources, which are located close together, mostly within a radius of 2 pc, and have velocities in a similar range with a median velocity dispersion of ~ 1 km/s. Using NH3, N2H+ and CS velocities we calculate near and far kinematic distances to 296 groups of ATLASGAL sources in the first quadrant and 393 groups in the fourth quadrant. We analyse HI self-absorption and HI absorption to resolve the kinematic distance ambiguity. We obtain a scale height of ~ 28+/-2 pc and displacement below the Galactic midplane of ~ -7+/-1 pc. Within distances from 2 to 18 kpc ATLASGAL clumps have a broad range of gas masses with a median of 1050 solar masses and a wide distribution of radii with a median of 0.4 pc. Their distribution in galactocentric radii is correlated with spiral arms. Using a statistically significant ATLASGAL sample we derive a power-law exponent of -2.2+/-0.1 of the clump mass function. This is consistent with the slope derived for clusters and with that of the stellar initial mass function. Examining the power-law index for different galactocentric distances and various source samples shows that it is independent of environment and evolutionary phase. Fitting the mass-size relationship by a power law gives a slope of 1.76+/-0.01 for cold sources such as IRDCs and warm clumps associated with HII regions.

We report the first detection of the J = 1 - 0 (102.6 GHz) rotational lines of CF+ (fluoromethylidynium ion) towards CygX-N63, a young and massive protostar of the Cygnus X region. This detection occurred as part of an unbiased spectral survey of this object in the 0.8-3 mm range, performed with the IRAM 30m telescope. The data were analyzed using a local thermodynamical equilibrium model (LTE model) and a population diagram in order to derive the column density. The line velocity (-4 km s-1) and line width (1.6 km s-1) indicate an origin from the collapsing envelope of the protostar. We obtain a CF+ column density of 4.10e11 cm-2. The CF+ ion is thought to be a good tracer for C+ and assuming a ratio of 10e-6 for CF+/C+, we derive a total number of C+ of 1.2x10e53 within the beam. There is no evidence of carbon ionization caused by an exterior source of UV photons suggesting that the protostar itself is the source of ionization. Ionization from the protostellar photosphere is not efficient enough. In contrast, X-ray ionization from the accretion shock(s) and UV ionization from outflow shocks could provide a large enough ionizing power to explain our CF+ detection. Surprisingly, CF+ has been detected towards a cold, massive protostar with no sign of an external photon dissociation region (PDR), which means that the only possibility is the existence of a significant inner source of C+. This is an important result that opens interesting perspectives to study the early development of ionized regions and to approach the issue of the evolution of the inner regions of collapsing envelopes of massive protostars. The existence of high energy radiations early in the evolution of massive protostars also has important implications for chemical evolution of dense collapsing gas and could trigger peculiar chemistry and early formation of a hot core.

Column-density maps of molecular clouds are one of the most important observables in the context of molecular cloud- and star-formation (SF) studies. With the Herschel satellite it is now possible to determine the column density from dust emission. We use observations and simulations to demonstrate how LOS contamination affects the column density probability distribution function (PDF). We apply a first-order approximation (removing a constant level) to the molecular clouds of Auriga, Maddalena, Carina and NGC3603. In perfect agreement with the simulations, we find that the PDFs become broader, the peak shifts to lower column densities, and the power-law tail of the PDF flattens after correction. All PDFs have a lognormal part for low column densities with a peak at Av~2, a deviation point (DP) from the lognormal at Av(DP)~4-5, and a power-law tail for higher column densities. Assuming a density distribution rho~r^-alpha, the slopes of the power-law tails correspond to alpha(PDF)=1.8, 1.75, and 2.5 for Auriga, Carina, and NGC3603 (alpha~1.5-2 is consistent gravitational collapse). We find that low-mass and high-mass SF clouds display differences in the overall column density structure. Massive clouds assemble more gas in smaller cloud volumes than low-mass SF ones. However, for both cloud types, the transition of the PDF from lognormal shape into power-law tail is found at the same column density (at Av~4-5 mag). Low-mass and high-mass SF clouds then have the same low column density distribution, most likely dominated by supersonic turbulence. At higher column densities, collapse and external pressure can form the power-law tail. The relative importance of the two processes can vary between clouds and thus lead to the observed differences in PDF and column density structure.

We report the results of 870-$\mu$m continuum observations, using the Large APEX Bolometer Camera (LABOCA), towards 77 class-II, 6.7-GHz methanol masers identified by the Methanol Multibeam (MMB) survey to map the thermal emission from cool dust towards these objects. These data complement a study of 630 methanol masers associated with compact dense clumps identified from the ATLASGAL survey. Compact dust emission is detected towards 70 sources, which implies a dust-association rate of 99% for the full MMB catalogue. Evaluation of the derived dust and maser properties leads us to conclude that the combined sample represents a single population tracing the same phenomenon. We find median clump masses of a few 10$^3$ M$\odot$ and that all but a handful of sources satisfy the mass-size criterion required for massive star formation. This study provides the strongest evidence of the almost ubiquitous association of methanol masers with massive, star-forming clumps. The fraction of methanol-maser associated clumps is a factor of ~2 lower in the outer Galaxy than the inner Galaxy, possibly a result of the lower metallicity environment of the former. We find no difference in the clump-mass and maser-luminosity distributions of the inner and outer Galaxy. The maser-pumping and clump-formation mechanisms are therefore likely to be relatively invariant to Galactic location. Finally, we use the ratio of maser luminosity and clump mass to investigate the hypothesis that the maser luminosity is a good indicator of the evolutionary stage of the embedded source, however, we find no evidence to support this.

Molecular outflows are a direct consequence of accretion, and therefore they represent one of the best tracers of accretion processes in the still poorly understood early phases of high-mass star formation. Previous studies suggested that the SiO abundance decreases with the evolution of a massive young stellar object probably because of a decay of jet activity, as witnessed in low-mass star-forming regions. We investigate the SiO excitation conditions and its abundance in outflows from a sample of massive young stellar objects through observations of the SiO(8-7) and CO(4-3) lines with the APEX telescope. Through a non-LTE analysis, we find that the excitation conditions of SiO increase with the velocity of the emitting gas. We also compute the SiO abundance through the SiO and CO integrated intensities at high velocities. For the sources in our sample we find no significant variation of the SiO abundance with evolution for a bolometric luminosity-to-mass ratio of between 4 and 50 $L_\odot/M_\odot$. We also find a weak increase of the SiO(8-7) luminosity with the bolometric luminosity-to-mass ratio. We speculate that this might be explained with an increase of density in the gas traced by SiO. We find that the densities constrained by the SiO observations require the use of shock models that include grain-grain processing. For the first time, such models are compared and found to be compatible with SiO observations. A pre-shock density of $10^5\, $cm$^{-3}$ is globally inferred from these comparisons. Shocks with a velocity higher than 25 km s$^{-1}$ are invoked for the objects in our sample where the SiO is observed with a corresponding velocity dispersion. Our comparison of shock models with observations suggests that sputtering of silicon-bearing material (corresponding to less than 10% of the total silicon abundance) from the grain mantles is occurring.

We used PdBI observations of SiO (2-1) to investigate the morphology and profile of the SiO emission within several massive dense clumps (MDCs) in Cygnus-X. We find that most molecular outflows are detected in both SiO and CO, although there are some cases of CO outflows with no SiO counterpart. We find a significant amount of narrow line SiO emission that appears to be unrelated to outflows. The fraction of the total SiO luminosity that is not associated with outflows is highly variable in the different MDCs (from 10% to 90%); this might be a problem when extrapolating outflow properties from SiO luminosities without resolving individual outflows. The extent of the narrow SiO emission varies from rather compact (~ 0.03 pc) to widespread (~0.2 pc), and its kinematics often differs from those found by other high-density tracers such as H13CO+. We find that the least centrally concentrated clumps with the least massive protostellar cores have the most widespread narrow SiO emission. In line with previous evidence of SiO emission associated with low-velocity shocks, we propose an evolutionary picture to explain the existence and distribution of narrow SiO line profiles. In this scenario, the least centrally condensed MDCs are at an early stage where the SiO emission traces shocks from the large-scale collapse of material onto the MDC (e.g. CygX-N40). As the MDC collapses, the SiO emission becomes more confined to the close surroundings of cores, tracing the post-shock material from the infalling MDC against the dense cores (e.g. CygX-N3, N12, and N48). At later stages, when single massive protostars are formed, the SiO luminosity is largely dominated by powerful outflows, and the weaker narrow component shows perhaps the last remnants of the initial collapse (e.g. CygX-N53 and N63).

We have selected the positions of 54 6.7GHz methanol masers from the Methanol Multibeam Survey catalogue, covering a range of longitudes between $20^{\circ}$ and $34^{\circ}$ of the Galactic Plane. These positions were mapped in the J=3-2 transition of both the $\rm{^{13}CO}$ and $\rm{C^{18}O}$ lines. A total of 58 $\rm{^{13}CO}$ emission peaks are found in the vicinity of these maser positions. We search for outflows around all $\rm{^{13}CO}$ peaks, and find evidence for high-velocity gas in all cases, spatially resolving the red and blue outflow lobes in 55 cases. Of these sources, 44 have resolved kinematic distances, and are closely associated with the 6.7GHz masers, a sub-set referred to as Methanol Maser Associated Outflows (MMAOs). We calculate the masses of the clumps associated with each peak using 870 $\rm{\mu m}$ continuum emission from the ATLASGAL survey. A strong correlation is seen between the clump mass and both outflow mass and mechanical force, lending support to models in which accretion is strongly linked to outflow. We find that the scaling law between outflow activity and clump masses observed for low-mass objects, is also followed by the MMAOs in this study, indicating a commonality in the formation processes of low-mass and high-mass stars.

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

In the low-mass regime, it is found that the gas-phase abundances of C-bearing molecules in cold starless cores rapidly decrease with increasing density, as the molecules form mantles on dust grains. We study CO depletion in 102 massive clumps selected from the ATLASGAL 870 micron survey, and investigate its correlation with evolutionary stage and with the physical parameters of the sources. Moreover, we study the gradients in [12C]/[13C] and [18O]/[17O] isotopic ratios across the inner Galaxy, and the virial stability of the clumps. We use low-J emission lines of CO isotopologues and the dust continuum emission to infer the depletion factor fD. RATRAN one-dimensional models were also used to determine fD and to investigate the presence of depletion above a density threshold. The isotopic ratios and optical depth were derived with a Bayesian approach. We find a significant number of clumps with a large fD, up to ~20. Larger values are found for colder clumps, thus for earlier evolutionary phases. For massive clumps in the earliest stages of evolution we estimate the radius of the region where CO depletion is important to be a few tenths of a pc. Clumps are found with total masses derived from dust continuum emission up to ~20 times higher than the virial mass, especially among the less evolved sources. These large values may in part be explained by the presence of depletion: if the CO emission comes mainly from the low-density outer layers, the molecules may be subthermally excited, leading to an overestimate of the dust masses. CO depletion in high-mass clumps seems to behave as in the low-mass regime, with less evolved clumps showing larger values for the depletion than their more evolved counterparts, and increasing for denser sources. The C and O isotopic ratios are consistent with previous determinations, and show a large intrinsic scatter.

The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is the largest and most sensitive systematic survey of the inner Galactic plane in the submillimetre wavelength regime. The observations were carried out with the Large APEX Bolometer Camera (LABOCA), an array of 295 bolometers observing at 870\,$\mu$m (345 GHz). Aim: In this research note we present the compact source catalogue for the 280\degr\ $ <\ell <$ 330\degr\ and 21\degr\ $ <\ell <$ 60\degr\ regions of this survey. Method: The construction of this catalogue was made with the source extraction routine \sex\ using the same input parameters and procedures used to analyse the inner Galaxy region presented in an earlier publication (i.e., 330\degr\ $ <\ell <$ 21\degr). Results: We have identified 3523 compact sources and present a catalogue of their properties. When combined with the regions already published this provides a comprehensive and unbiased database of ~10163 massive, dense clumps located across the inner Galaxy.

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.

Context: Star formation efficiency (SFE) theories are currently based on statistical distributions of turbulent cloud structures and a simple model of star formation from cores. They remain poorly tested, especially at the highest densities. Aims: We investigate the effects of gas density on the SFE through measurements of the core formation efficiency (CFE). With a total mass of $\sim2\times10^4$ M$_\odot$, the W43-MM1 ridge is one of the most convincing candidate precursor of starburst clusters and thus one of the best place to investigate star formation. Methods: We used high-angular resolution maps obtained at 3 mm and 1 mm within W43-MM1 with the IRAM Plateau de Bure Interferometer to reveal a cluster of 11 massive dense cores (MDCs), and, one of the most massive protostellar cores known. An Herschel column density image provided the mass distribution of the cloud gas. We then measured the 'instantaneous' CFE and estimated the SFE and the star formation rate (SFR) within subregions of the W43-MM1 ridge. Results: The high SFE found in the ridge ($\sim$6% enclosed in $\sim$8 pc$^3$) confirms its ability to form a starburst cluster. There is however a clear lack of dense cores in the northern part of the ridge, which may be currently assembling. The CFE and the SFE are observed to increase with volume gas density while the SFR steeply decreases with the virial parameter, $\alpha_{vir}$. Statistical models of the SFR may well describe the outskirts of the W43-MM1 ridge but struggle to reproduce its inner part, which corresponds to measurements at low $\alpha_{vir}$. It may be that ridges do not follow the log-normal density distribution, Larson relations, and stationary conditions forced in the statistical SFR models.

We investigate the physical and chemical processes at work during the formation of a massive protostar based on the observation of water in an outflow from a very young object previously detected in H2 and SiO in the IRAS 17233-3606 region. We estimated the abundance of water to understand its chemistry, and to constrain the mass of the emitting outflow. We present new observations of shocked water obtained with the HIFI receiver onboard Herschel. We detected water at high velocities in a range similar to SiO. We self-consistently fitted these observations along with previous SiO data through a state-of-the-art, one-dimensional, stationary C-shock model. We found that a single model can explain the SiO and H2O emission in the red and blue wings of the spectra. Remarkably, one common area, similar to that found for H2 emission, fits both the SiO and H2O emission regions. This shock model subsequently allowed us to assess the shocked water column density, N(H2O)=1.2x10^{18} cm^{-2}, mass, M(H2O)=12.5 M_earth, and its maximum fractional abundance with respect to the total density, x(H2O)=1.4x10^{-4}. The corresponding water abundance in fractional column density units ranges between 2.5x10^{-5} and 1.2x10^{-5}, in agreement with recent results obtained in outflows from low- and high-mass young stellar objects.

The formation processes and the evolutionary stages of high-mass stars are poorly understood compared to low-mass stars. Large-scale surveys are needed to provide an unbiased census of high column density sites which can potentially host precursors to high-mass stars. Here we use the ATLASGAL survey covering 420 sq. degree of the Galactic plane at 870 $\mu$m; and use the MRE-GLC method to identify the population of embedded sources throughout the inner Galaxy. We identify in total 10861 compact sub-millimeter sources with fluxes above 5 sigma. Completeness tests show that our catalogue is 97% complete above 5 sigma and >99% complete above 7$\sigma$. We correlate this sample with mid-infrared point source catalogues (MSX at 21.3 $\mu$m and WISE at 22 $\mu$m) and determine a lower limit of ~33% that are associated with embedded protostellar objects. We note that the proportion of clumps associated with mid-infrared sources increases with increasing flux density, achieving a rather constant fraction of ~75% of all clumps with fluxes over 5 Jy/beam being associated with star-formation. Examining the source counts as a function of Galactic longitude we are able to identify the most prominent star forming regions in the Galaxy. From the fraction of the likely massive quiescent clumps (~25%) we estimate a formation time-scale of ~7.5+/-2.5 $\times$ 10$^4$yr for the deeply embedded phase before the emergence of luminous YSOs. Such a short duration for the formation of high-mass stars in massive clumps clearly proves that the earliest phases have to be dynamic with supersonic motions.

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

We estimate the turbulent ambipolar diffusion length scale and magnetic field strength in the massive dense cores CygX-N03 and CygX-N53, located in the Cygnus-X star-forming region. The method we use requires comparing the velocity dispersions in the spectral line profiles of the coexistent ion and neutral pair H13CN and H13CO+ (J=1-0) at different length scales. We fit Kolmogorov-type power laws to the lower envelopes of the velocity dispersion spectra of the two species. This allows to calculate the turbulent ambipolar diffusion scale, which in turn determines the plane-of-the-sky magnetic field strength. We find turbulent ambipolar diffusion length scales of 3.8+-0.1 mpc and 21.2+-0.4 mpc, and magnetic field strengths of 0.33 mG and 0.76 mG for CygX-N03 and CygX-N53, respectively. These magnetic field values have uncertainties of a factor of a few. Despite a lower signal-to-noise ratio of the data in CygX-N53 than in CygX-N03, and the caveat that its stronger field might stem in part from projection effects, the difference in field strengths suggests different fragmentation activities of the two cores. Even though the quality of our data, obtained with the IRAM Plateau de Bure Interferometer (PdBI), is somewhat inferior to previous single-dish data, we demonstrate that this method is suited also for observations at high spatial resolution.

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

The 30 Myr old A3-type star HD 21997 is one of the two known debris dust disks having a measurable amount of cold molecular gas. With the goal of understanding the physical state, origin, and evolution of the gas in young debris disks, we obtained CO line observations with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we report on the detection of 12CO and 13CO in the J=2-1 and J=3-2 transitions and C18O in the J=2-1 line. The gas exhibits a Keplerian velocity curve, one of the few direct measurements of Keplerian rotation in young debris disks. The measured CO brightness distribution could be reproduced by a simple star+disk system, whose parameters are r_in < 26 AU, r_out = 138 +/- 20 AU, M_*=1.8 +0.5 -0.2 M_Sun, and i = 32.6 +/- 3.1 degrees. The total CO mass, as calculated from the optically thin C18O line, is about (4-8) x 10^-2 M_Earth, while the CO line ratios suggest a radiation temperature on the order of 6-9 K. Comparing our results with those obtained for the dust component of the HD 21997 disk from the ALMA continuum observations by Mo\'or et al., we conclude that comparable amounts of CO gas and dust are present in the disk. Interestingly, the gas and dust in the HD 21997 system are not co-located, indicating a dust-free inner gas disk within 55 AU of the star. We explore two possible scenarios for the origin of the gas. A secondary origin, which involves gas production from colliding or active planetesimals, would require unreasonably high gas production rates and would not explain why the gas and dust are not co-located. We propose that HD 21997 is a hybrid system where secondary debris dust and primordial gas coexist. HD 21997, whose age exceeds both the model predictions for disk clearing and the ages of the oldest T Tauri-like or transitional gas disks in the literature, may be a key object linking the primordial and the debris phases of disk evolution.

Circumstellar disks around stars older than 10 Myr are expected to be gas-poor. There are, however, two examples of old (30-40 Myr) debris-like disks containing a detectable amount of cold CO gas. Here we present ALMA and Herschel Space Observatory observations of one of these disks, around HD 21997, and study the distribution and origin of the dust and its connection to the gas. Our ALMA continuum images at 886um clearly resolve a broad ring of emission within a diameter of ~4.5 arcsec, adding HD 21997 to the dozen debris disks resolved at (sub)millimeter wavelengths. Modeling the morphology of the ALMA image with a radiative transfer code suggests inner and outer radii of ~55 and ~150 AU, and a dust mass of 0.09 M_Earth. Our data and modeling hints at an extended cold outskirt of the ring. Comparison with the morphology of the CO gas in the disk reveals an inner dust-free hole where gas nevertheless can be detected. Based on dust grain lifetimes, we propose that the dust content of this gaseous disk is of secondary origin produced by planetesimals. Since the gas component is probably primordial, HD 21997 is one of the first known examples of a hybrid circumstellar disk, a so-far little studied late phase of circumstellar disk evolution.

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey aims to characterise the physical and chemical evolution of high-mass star-forming clumps. Exploiting the unique broad frequency range and on-the-fly mapping capabilities of the Australia Telescope National Facility Mopra 22 m single-dish telescope, MALT90 has obtained 3' x 3' maps toward ~2000 dense molecular clumps identified in the ATLASGAL 870 um Galactic plane survey. The clumps were selected to host the early stages of high-mass star formation and to span the complete range in their evolutionary states (from prestellar, to protostellar, and on to HII regions and photodissociation regions). Because MALT90 mapped 16 lines simultaneously with excellent spatial (38") and spectral (0.11 km/s) resolution, the data reveal a wealth of information about the clump's morphologies, chemistry, and kinematics. In this paper we outline the survey strategy, observing mode, data reduction procedure, and highlight some early science results. All MALT90 raw and processed data products are available to the community. With its unprecedented large sample of clumps, MALT90 is the largest survey of its type ever conducted and an excellent resource for identifying interesting candidates for high resolution studies with ALMA.

Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star HD95086. The strong infrared excess of the system indicates that, similarly to HR8799, {\ss} Pic, and Fomalhaut, the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD95086 disk, and to investigate its possible interaction with the newly discovered planet, here we present new optical, infrared and millimeter observations. We detected no CO emission, excluding the possibility of an evolved gaseous primordial disk. Simple blackbody modeling of the spectral energy distribution suggests the presence of two spatially separate dust belts at radial distances of 6 and 64 AU. Our resolved images obtained with the Herschel Space Observatory reveal a characteristic disk size of ~6.0x5.4 arcsec (540x490 AU) and disk inclination of ~25 degree. Assuming the same inclination for the planet candidate's orbit, its re-projected radial distance from the star is 62 AU, very close to the blackbody radius of the outer cold dust ring. The structure of the planetary system at HD95086 resembles the one around HR8799. Both systems harbor a warm inner dust belt and a broad colder outer disk and giant planet(s) between the two dusty regions. Modelling implies that the candidate planet can dynamically excite the motion of planetesimals even out to 270 AU via their secular perturbation if its orbital eccentricity is larger than about 0.4. Our analysis adds a new example to the three known systems where directly imaged planet(s) and debris disks co-exist.

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

We characterize the Millimeter Astronomy Legacy Team 90 GHz (MALT90) Survey and the Mopra telescope at 90 GHz. We combine repeated position-switched observations of the source G300.968+01.145 with a map of the same source in order to estimate the pointing reliability of the position-switched observations and, by extension, the MALT90 survey; we estimate our pointing uncertainty to be 8 arcseconds. We model the two strongest sources of systematic gain variability as functions of elevation and time-of-day and quantify the remaining absolute flux uncertainty. Corrections based on these two variables reduce the scatter in repeated observations from 12-25% down to 10-17%. We find no evidence for intrinsic source variability in G300.968+01.145. For certain applications, the corrections described herein will be integral for improving the absolute flux calibration of MALT90 maps and other observations using the Mopra telescope at 90 GHz.

Using the 870-$\mu$m APEX Telescope Large Area Survey of the Galaxy (ATLASGAL), we have identified 577 submillimetre continuum sources with masers from the methanol multibeam (MMB) survey in the region $280\degr < \ell < 20\degr$; $|\,b\,| < 1.5\degr$. 94\,per\,cent of methanol masers in the region are associated with sub-millimetre dust emission. We estimate masses for ~450 maser-associated sources and find that methanol masers are preferentially associated with massive clumps. These clumps are centrally condensed, with envelope structures that appear to be scale-free, the mean maser position being offset from the peak column density by 0.0 \pm 4". Assuming a Kroupa initial mass function and a star-formation efficiency of ~30\,per\,cent, we find that over two thirds of the clumps are likely to form clusters with masses >20\,\msun. Furthermore, almost all clumps satisfy the empirical mass-size criterion for massive star formation. Bolometric luminosities taken from the literature for ~100 clumps range between ~100 and 10$^6$\,\lsun. This confirms the link between methanol masers and massive young stars for 90\,per\,cent of our sample. The Galactic distribution of sources suggests that the star-formation efficiency is significantly reduced in the Galactic-centre region, compared to the rest of the survey area, where it is broadly constant, and shows a significant drop in the massive star-formation rate density in the outer Galaxy. We find no enhancement in source counts towards the southern Scutum-Centaurus arm tangent at $\ell ~ 315\degr$, which suggests that this arm is not actively forming stars.

Context. The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is the first systematic survey of the inner Galactic plane in the sub-millimetre. The observations were carried out with the Large APEX Bolometer Camera (LABOCA), an array of 295 bolometers observing at 870 microns (345 GHz). Aims. Here we present a first version of the compact source catalogue extracted from this survey. This catalogue provides an unbiased database of dusty clumps in the inner Galaxy. Methods. The construction of this catalogue was made using the source extraction routine SExtractor. We have cross-associated the obtained sources with the IRAS and MSX catalogues, in order to constrain their nature. Results. We have detected 6639 compact sources in the range from 330 < l < 21 degrees and |b| < 1.5 degrees. The catalogue has a 99% completeness for sources with a peak flux above 6 sigma, which corresponds to a flux density of ~0.4 Jy/beam. The parameters extracted for sources with peak fluxes below the 6 sigma completeness threshold should be used with caution. Tests on simulated data find the uncertainty in the flux measurement to be ~12%, however, in more complex regions the flux values can be overestimated by a factor of 2 due to the additional background emission. Using a search radius of 30" we found that 40% of ATLASGAL compact sources are associated with an IRAS or MSX point source, but, ~50% are found to be associated with MSX 21 microns fluxes above the local background level, which is probably a lower limit to the actual number of sources associated with star formation. Conclusions. Although infrared emission is found towards the majority of the clumps detected, this catalogue is still likely to include a significant number of clumps that are devoid of star formation activity and therefore excellent candidates for objects in the coldest, earliest stages of (high-mass) star formation.

(Abridged) We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654--1897 GHz using HIFI on board Herschel, in DR21(OH), a high mass star forming region. Several transitions of CCH and c-C3H2 have also been observed with HIFI and the IRAM 30m telescope. A gas and grain warm-up model was used to identify the primary C3 forming reactions in DR21(OH). We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is \sim6x10^{-10} and \sim3x10^{-9} respectively. For CCH and c-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3, and c-C3H2 abundances are most consistent with a chemical model with n(H2)\sim5x10^{6} cm^-3 post-warm-up dust temperature, T_max =30 K and a time of \sim0.7-3 Myr. Post warm-up gas phase chemistry of CH4 released from the grain at t\sim 0.2 Myr and lasting for 1 Myr can explain the observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. The observed lower C3 abundance in MM1 as compared to MM2 and the envelope could be indicative of destruction of C3 in the more evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived for DR21(OH) in other studies.

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

The hydroxyl radical (OH) is found in various environments within the interstellar medium (ISM) of the Milky Way and external galaxies, mostly either in diffuse interstellar clouds or in the warm, dense environments of newly formed low-mass and high-mass stars, i.e, in the dense shells of compact and ultracompact HII regions (UCHIIRs). Until today, most studies of interstellar OH involved the molecule's radio wavelength hyperfine structure (hfs) transitions. These lines are generally not in LTE and either masing or over-cooling complicates their interpretation. In the past, observations of transitions between different rotational levels of OH, which are at far-infrared wavelengths, have suffered from limited spectral and angular resolution. Since these lines have critical densities many orders of magnitude higher than the radio wavelength ground state hfs lines and are emitted from levels with more than 100 K above the ground state, when observed in emission, they probe very dense and warm material. We probe the warm and dense molecular material surrounding the UCHIIR/OH maser sources W3(OH), G10.62-0.39 and NGC 7538 IRS1 by studying the $^2\Pi_{{1/2}}, J = {3/2} - {1/2}$ rotational transition of OH in emission and, toward the last source also the molecule's $^2\Pi_{3/2}, J = 5/2 - 3/2$ ground-state transition in absorption. We used the Stratospheric Observatory for Infrared Astronomy (SOFIA) to observe these OH lines, which are near 1.84 THz ($163 \mu$m) and 2.51 THz ($119.3 \mu$m). We clearly detect the OH lines, some of which are blended with each other. Employing non-LTE radiative transfer calculations we predict line intensities using models of a low OH abundance envelope versus a compact, high-abundance source corresponding to the origin of the radio OH lines.

Gas-rich primordial disks and tenuous gas-poor debris disks are usually considered as two distinct evolutionary phases of the circumstellar matter. Interestingly, the debris disk around the young main-sequence star 49 Ceti possesses a substantial amount of molecular gas, and possibly represents the missing link between the two phases. Motivated to understand the evolution of the gas component in circumstellar disks via finding more 49 Ceti-like systems, we carried out a CO J=3-2 survey with Atacama Pathfinder EXperiment, targeting 20 infrared-luminous debris disks. These systems fill the gap between primordial and old tenuous debris disks in terms of fractional luminosity. Here we report on the discovery of a second 49 Ceti-like disk around the 30 Myr old A3-type star HD21997, a member of the Columba Association. This system was also detected in the CO(2-1) transition, and the reliable age determination makes it an even clearer example of an old gas-bearing disk than 49 Ceti. While the fractional luminosities of HD21997 and 49 Ceti are not particularly high, these objects seem to harbor the most extended disks within our sample. The double-peaked profiles of HD21997 were reproduced by a Keplerian disk model combined with the LIME radiative transfer code. Based on their similarities, 49 Ceti and HD21997 may be the first representatives of a so far undefined new class of relatively old (>~8 Myr), gaseous dust disks. From our results, neither primordial origin nor steady secondary production from icy planetesimals can unequivocally explain the presence of CO gas in the disk of HD21997.

DR21(OH) is a pc-scale massive, 7000 Msun clump hosting three massive dense cores (MDCs) at an early stage of their evolution. We present a high angular-resolution mosaic, covering 70" by 100", with the IRAM PdBI at 3 mm to trace the dust continuum emission and the N2H+ (J=1-0) and CH3CN (J=5-4) molecular emission. The cold, dense gas traced by the compact emission in N2H+ is associated with the three MDCs and shows several velocity components towards each MDC. These velocity components reveal local shears in the velocity fields which are best interpreted as convergent flows. Moreover, we report the detection of weak extended emission from CH3CN at the position of the N2H+ velocity shears. We propose that this extended CH3CN emission is tracing warm gas associated with the low-velocity shocks expected at the location of convergence of the flows where velocity shears are observed. This is the first detection of low-velocity shocks associated with small (sub-parsec) scale convergent flows which are proposed to be at the origin of the densest structures and of the formation of (high-mass) stars. In addition, we propose that MDCs may be active sites of star-formation for more than a crossing time as they continuously receive material from larger scale flows as suggested by the global picture of dynamical, gravity driven evolution of massive clumps which is favored by the present observations.

Although photometric and spectroscopic surveys with the Spitzer Space Telescope increased remarkably the number of well studied debris disks around A-type and Sun-like stars, detailed analyzes of debris disks around F-type stars remained less frequent. Using the MIPS camera and the IRS spectrograph we searched for debris dust around 82 F-type stars with Spitzer. We found 27 stars that harbor debris disks, nine of which are new discoveries. The dust distribution around two of our stars, HD 50571 and HD 170773, was found to be marginally extended on the 70um MIPS images. Combining the MIPS and IRS measurements with additional infrared and submillimeter data, we achieved excellent spectral coverage for most of our debris systems. We have modeled the excess emission of 22 debris disks using a single temperature dust ring model and of 5 debris systems with two-temperature models. The latter systems may contain two dust rings around the star. In accordance with the expected trends, the fractional luminosity of the disks declines with time, exhibiting a decay rate consistent with the range of model predictions. We found the distribution of radial dust distances as a function of age to be consistent with the predictions of both the self stirred and the planetary stirred disk evolution models. A more comprehensive investigation of the evolution of debris disks around F-type stars, partly based on the presented data set, will be the subject of an upcoming paper.

We study the kinematic properties of dense gas surrounding massive protostars recognized by Bontemps et a. (2010) in a sample of five Massive Dense Cores in Cygnus-X. We investigate whether turbulent support plays a major role in stabilizing the core against fragmentation into Jeans-mass objects or alternatively, the observed kinematics could indicate a high level of dynamics. We present IRAM 30m single-dish (HCO+ and H13CO+) and IRAM PdBI high angular-resolution observations of dense gas tracers (H13CO+ and H13CN) to reveal the kinematics of molecular gas at scales from 0.03 to 0.1 pc. Radiative transfer modeling shows that H13CO+ is depleted within the envelopes of massive protostars and traces the bulk of material surrounding the protostars rather than their inner envelopes. H13CN shows a better correspondence with the peak of the continuum emission, possibly due to abundance anomalies and specific chemistry in the close vicinity of massive protostars. Analyzing the line-widths we show that the observed line-dispersion of H13CO+ at the scale of MDCs is smaller than expected from the quasi-static, turbulent-core model. At large-scales, global organized bulk motions are identified for 3 of the MDCs. At small-scales, several spectral components are identified in all MDCs showing filamentary structures and intrinsic velocity gradients towards the continuum peaks. The dynamics of these flows show diversity among the sample and we link this to the specific fragmentation properties of the MDCs. No clear evidence is found for a turbulence regulated, equilibrium scenario within the sample of MDCs. We propose a picture in which MDCs are not in equilibrium and their dynamics is governed by small-scale converging flows, which may initiate star-formation via their shears.

The formation of massive stars is a highly complex process in which it is not clear whether the star-forming gas is in global gravitational collapse or in an equilibrium state, supported by turbulence. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, the filament containing the well-known sources DR21 and DR21(OH), we expect to find observational signatures that allow to discriminate between the two views. We use molecular line data from our 13CO 1-0, CS 2-1, and N2H+ 1-0 survey of the Cygnus X region obtained with the FCRAO and high-angular resolution observations of CO, CS, HCO+, N2H+, and H2CO, obtained with the IRAM 30m telescope. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of highest column-density, i.e. dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to rotation. Infall signatures in optically thick line profiles of HCO+ and 12CO are observed along and across the whole DR21 filament. From modelling the observed spectra, we obtain a typical infall speed of 0.6 km/s and mass accretion rates of the order of a few 10^-3 Msun/yr for the two main clumps constituting the filament. These massive (4900 and 3300 Msun) clumps are both gravitationally contracting. All observed kinematic features in the DR21 filament can be explained if it is formed by the convergence of flows at large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity at larger scales or from other processes can not be settled with the present study. The observed velocity field and velocity dispersion are consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent flows.

We study the rotational properties of magnetized and self-gravitating molecular cloud cores formed in 2 very high resolution 3D molecular cloud simulations.The simulations have been performed using the code RAMSES at an effective resolution of 4096^3.One simulation represents a mildly magnetically-supercritical cloud and the other a strongly magnetically-supercritical cloud.A noticeable difference between the 2 simulations is the core formation efficiency (CFE) of the high density cores.In the strongly supercritical simulations the CFE is ~33 % per free-fall time of the cloud tff,cl, whereas in the mildly supercritical simulations this value goes down to ~6%/tff,cl. A comparison of the intrinsic specific angular momentum j3D distributions of the cores with the distribitions of j2D derived using synthetic 2D velocity maps of the cores,shows that the synthetic observations tend to overestimate the true value of j by a factor of ~10.The origin of this discrepancy lies in the fact that contrary to the intrinsic determination which sums up the individual gas parcels contributions to j, the determination of j using the observational procedure which is based on a measurement on the global velocity gradient under the hypothesis of uniform rotation smoothes out the complex fluctuations present in the 3D velocity field. Our results provide a natural explanation for the discrepancy by a factor ~10 observed between the intrinsic 3D distributions of j and the corresponding distributions derived in real observations.We suggest that measurements of j which are based on the measurement of the observed global velocity gradients may need to be reduced by a factor of ~10 in order to derive a more accurate estimate of j in the cores.

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

We present Plateau de Bure interferometer observations obtained in continuum at 1.3 and 3.5 mm towards the six most massive and young (IR-quiet) dense cores in Cygnus X. Located at only 1.7 kpc, the Cygnus X region offers the opportunity of reaching small enough scales (of the order of 1700 AU at 1.3 mm) to separate individual collapsing objects. The cores are sub-fragmented with a total of 23 fragments inside 5 cores. Only the most compact core, CygX-N63, could actually be a single massive protostar with an envelope mass as large as 60 Msun. The fragments in the other cores have sizes and separations similar to low-mass pre-stellar and proto-stellar condensations in nearby protoclusters, and are probably of the same nature. A total of 9 out of these 23 protostellar objects are found to be probable precursors of OB stars with envelope masses ranging from 6 to 23 Msun. The level of fragmentation is globally higher than in the turbulence regulated, monolithic collapse scenario, but is not as high as expected in a pure gravo-turbulent scenario where the distribution of mass is dominated by low-mass protostars/stars. Here, the fractions of the total core masses in the high-mass fragments are reaching values as high as 28, 44, and 100 % in CygX-N12, CygX-N53, and CygX-N63, respectively, much higher than what an IMF-like mass distribution would predict. The increase of the fragmentation efficiency as a function of density in the cores is proposed to be due to the increasing importance of self-gravity leading to gravitational collapse at the scale of the dense cores. At the same time, the cores tend to fragment into a few massive protostars within their central regions. We are therefore probably witnessing here the primordial mass segregation of clusters in formation.