In recent years, a plethora of high spectral resolution observations of sub-mm and FIR transitions of methylidene (CH), have demonstrated this radical to be a valuable proxy for H2, that can be used for characterising molecular gas within the interstellar medium (ISM) on a Galactic scale, including the CO-dark component. Here we report the discovery of the 13CH isotopologue in the ISM using the upGREAT receiver on board SOFIA. We have detected the three hyperfine structure components of the 2THz frequency transition from its ground-state toward four high-mass star-forming regions and determine 13CH column densities. The ubiquity of molecules containing carbon in the ISM has turned the determination of the ratio between the abundances of carbon's two stable isotopes, 12C/13C, into a cornerstone for Galactic chemical evolution studies. Whilst displaying a rising gradient with Galactocentric distance, this ratio, when measured using observations of different molecules (CO, H2CO, and others) shows systematic variations depending on the tracer used. These observed inconsistencies may arise from optical depth effects, chemical fractionation or isotope-selective photo-dissociation. Formed from C+ either via UV-driven or turbulence-driven chemistry, CH reflects the fractionation of C+, and does not show any significant fractionation effects unlike other molecules previously used to determine the 12C/13C isotopic ratio which make it an ideal tracer for the 12C/13C ratio throughout the Galaxy. Therefore, by comparing the derived column densities of 13CH with previously obtained SOFIA data of the corresponding transitions of the main isotopologue 12CH, we derive 12C/13C isotopic ratios toward Sgr B2(M), G34.26+0.15, W49(N) and W51E. Adding our values derived from 12/13CH to previous calculations of the Galactic isotopic gradient we derive a revised value of 12C/13C = 5.85(0.50)R_GC + 15.03(3.40).

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.

Aims: Our aim is to identify the dominant molecular cooling lines and characteristic emission features in the 1.3 mm window of distinct regions in the northern part of the Orion A molecular cloud. By defining and analysing template regions, we also intend to help with the interpretation of observations from more distant sources which cannot be easily spatially resolved. Methods: We analyse an imaging line survey covering the area of OMC-1 to OMC-3 from 200.2 to 281.8 GHz obtained with the PI230 receiver at the APEX telescope. Masks are used to define regions with distinct properties (e.g. column density or temperature ranges) from which we obtain averaged spectra. Lines of 29 molecular species (55 isotopologues) are fitted for each region to obtain the respective total intensity. Results: We find that strong sources like Orion KL have a clear impact on the emission on larger scales. Although not spatially extended, their line emission contributes substantially to spectra averaged over large regions. Conversely, the emission signatures of dense, cold regions like OMC-2 and OMC-3 (e.g. enhanced N$_{2}$H$^{+}$ emission and low HCN/HNC ratio) seem to be difficult to pick up on larger scales, where they are eclipsed by signatures of stronger sources. In all regions, HCO$^{+}$ appears to contribute between 3% and 6% to the total intensity, the most stable value for all bright species. N$_{2}$H$^{+}$ shows the strongest correlation with column density, but not with typical high-density tracers like HCN, HCO$^{+}$, H$_{2}$CO, or HNC. Common line ratios associated with UV illumination, CN/HNC and CN/HCO$^{+}$, show ambiguous results on larger scales, suggesting that the identification of UV illuminated material may be more challenging. The HCN/HNC ratio may be related to temperature over varying scales.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

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 (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

We have constructed the largest sample of dust-associated class II 6.7 GHz methanol masers yet obtained. New measurements from the the Methanol MultiBeam (MMB) Survey were combined with the 870 $\mu$m APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) and the 850 $\mu$m JCMT Plane Survey (JPS). Together with two previous studies we have now identified the host clumps for 958 methanol masers across the Galactic Plane, covering approximately 99% of the MMB catalogue and increasing the known sample of dust-associated masers by over 30%. We investigate correlations between the physical properties of the clumps and masers using distances and luminosities drawn from the literature. Clumps hosting methanol masers are significantly more compact and have higher volume densities than the general population of clumps. We determine a minimum volume density threshold of $n$(H$_2$) $\geq 10^4$ cm$^{-3}$ for the efficient formation of intermediate- and high-mass stars. We find 6.7 GHz methanol masers are associated with a distinct part of the evolutionary process ($L_{\rm bol}$/$M_{\rm fwhm}$ ratios of between 10$^{0.6}$ and 10$^{2.2}$) and have well defined turning on and termination points. We estimate the lower limit for the mass of embedded objects to be $\geq$ 6 M$_{\odot}$ and the statistical lifetime of the methanol maser stage to be $\sim$ 3.3$\times$10$^{4}$ yrs. This suggests that methanol masers are indeed reliable tracers of high mass star formation, and indicates that the evolutionary period traced by this marker is relatively rapid.

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.

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

OMC-2/3 is one of the nearest embedded cluster-forming region that includes intermediate-mass protostars at early stages of evolution. A previous CO (3--2) mapping survey towards this region revealed outflow activity related to sources at different evolutionary phases. The present work aims to study the warm gas in the high-velocity emission from several outflows found in CO (3--2) emission by previous observations, determine their physical conditions, and make comparison with previous results in low-mass star-forming regions. We use the CHAMP+ heterodyne array on the APEX telescope to map the CO (6--5) and CO (7--6) emission in the OMC-2 FIR 6 and OMC-3 MMS 1-6 regions, and to observe $^{13}$CO (6--5) at selected positions. We analyze these data together with previous CO (3--2) observations. In addition, we mapped the SiO (5--4) emission in OMC-2 FIR 6. The CO (6--5) emission was detected in most of the outflow lobes in the mapped regions, while the CO (7--6) was found mostly from the OMC-3 outflows. In the OMC-3 MMS 5 outflow, a previously undetected extremely high velocity gas was found in CO (6--5). This extremely high velocity emission arises from the regions close to the central object MMS 5. Radiative transfer models revealed that the high-velocity gas from MMS 5 outflow consists of gas with $n_{\rm H_2}=$10$^4$--10$^{5}$ cm$^{-3}$ and $T>$ 200 K, similar to what is observed in young Class 0 low-mass protostars. For the other outflows, values of $n_{\rm H_2}>$10$^{4}$ cm$^{-3}$ were found. The physical conditions and kinematic properties of the young intermediate-mass outflows presented here are similar to those found in outflows from Class 0 low-mass objects. Due to their excitation requirements, mid$-J$ CO lines are good tracers of extremely high velocity gas in young outflows likely related to jets.

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.

In this study we analyze 70 radio continuum sources associated with dust clumps and considered to be candidates for the earliest stages of high-mass star formation. The detection of these sources was reported by Rosero et al. (2016), who found most of them to show weak (${\scriptstyle <}$1 mJy) and compact (${\scriptstyle <}\,$0.6$^{\prime \prime}$) radio emission. Herein, we used the observed parameters of these sources to investigate the origin of the radio continuum emission. We found that at least $\sim 30\%$ of these radio detections are most likely ionized jets associated with high-mass protostars, but for the most compact sources we cannot discard the scenario that they represent pressure-confined HII regions. This result is highly relevant for recent theoretical models based on core accretion that predict the first stages of ionization from high-mass stars to be in the form of jets. Additionally, we found that properties such as the radio luminosity as a function of the bolometric luminosity of ionized jets from low and high-mass stars are extremely well-correlated. Our data improve upon previous studies by providing further evidence of a common origin for jets independently of luminosity.

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.

We present Submillimeter Array (SMA) observations in the CO J=3-2, SiO J=5-4 and 8-7, and SO 9_8-8_7 lines, as well as Atacama Pathfinder EXperiment (APEX) observations in the CO J=6-5 line, of an extremely high-velocity and jet-like outflow in high-mass star-forming region HH 80--81. The outflow is known to contain two prominent molecular bullets, namely B1 and B2, discovered from our previous SMA CO J=2-1 observations. While B1 is detected in all the CO, SiO, and SO lines, B2 is only detected in CO lines. The CO 3-2/2-1 line ratio in B1 is clearly greater than that in B2. We perform a large velocity gradient analysis of the CO lines and derive a temperature of 70--210 K for B1 and 20--50 K for B2. Taking into account the differences in the velocity, distance from the central source, excitation conditions, and chemistry between the two bullets, we suggest that the bullets are better explained by direct ejections from the innermost vicinity of the central high-mass protostar, and that we are more likely observing the molecular component of a primary wind rather than entrained or swept-up material from the ambient cloud. These findings further support our previous suggestions that the molecular bullets indicate an episodic, disk-mediated accretion in the high-mass star formation process.

(Abridged) We use HIFI maps of the 987 GHz H2O 2(02)-1(11) emission to measure the sizes and shapes of 19 high-mass protostellar envelopes. To identify infall, we use HIFI spectra of the optically thin C18O 9-8 and H2O-18 1(11)-0(00) lines. The high-J C18O line traces the warm central material and redshifted H2O-18 1(11)-0(00) absorption indicates material falling onto the warm core. We probe small-scale chemical differentiation by comparing H2O 752 and 987 GHz spectra with those of H2O-18. Our measured radii of the central part of the H2O 2(02)-1(11) emission are 30-40% larger than the predictions from spherical envelope models, and axis ratios are <2, which we consider good agreement. For 11 of the 19 sources, we find a significant redshift of the H2O-18 1(11)-0(00) line relative to C18O 9-8. The inferred infall velocities are 0.6-3.2 km/s, and estimated mass inflow rates range from 7e-5 to 2e-2 M0/yr, with the highest mass inflow rates occurring toward the sources with the highest masses, and possibly the youngest ages. The other sources show either expanding motions or H2O-18 lines in emission. The H2O-18 1(11)-0(00) line profiles are remarkably similar to the "differences" between the H2O 2(02)-1(11) and 2(11)-2(02) profiles, suggesting that the H2O-18 line and the H2O 2(02)-1(11) absorption originate just inside the radius where water evaporates from grains, typically 1000-5000 au from the center. In some sources, the H2O-18 line is detectable in the outflow, where no C18O emission is seen. Together, the H2O-18 absorption and C18O emission profiles show that the water abundance around high-mass protostars has at least three levels: low in the cool outer envelope, high within the 100 K radius, and very high in the outflowing gas. Thus, despite the small regions, the combination of lines presented here reveals systematic inflows and chemical information about the outflows.

Hydrocarbons are ubiquitous in the interstellar medium, but their formation is still not well understood, depending on the physical environment they are found in. M8 is host to one of the brightest HII regions and PDRs in our galaxy. Using the APEX, and the IRAM 30 m telescopes, we performed a line survey toward Herschel 36 (Her 36), which is the main ionizing stellar system in M8, and an imaging survey within 1.3 $\times$ 1.3 pc around Her 36 of various transitions of C$_{2}$H and c-C$_{3}$H$_{2}$. We used both LTE and non-LTE methods to determine the physical conditions of the emitting gas along with the column densities and abundances of the observed species, which we compared with (updated) gas phase photochemical PDR models. In order to examine the role of PAHs in the formation of small hydrocarbons and to investigate their association with M8, we compared archival GLIMPSE 8 $\mu$m and the SPIRE 250 $\mu$m continuum images with the C$_{2}$H emission maps. We observed a total of three rotational transitions of C$_{2}$H with their hyperfine structure components and four rotational transitions of c-C$_{3}$H$_{2}$ with ortho and para symmetries toward M8. Fragmentation of PAHs seems less likely to contribute to the formation of small hydrocarbons as the 8 $\mu$m emission does not follow the distribution of C$_{2}$H emission, which is more associated with the molecular cloud. From the quantitative analysis, we obtained abundances of $\sim$ 10$^{-8}$ and 10$^{-9}$ for C$_{2}$H and c-C$_{3}$H$_{2}$ respectively, and volume densities of the hydrocarbon emitting gas in the range $n(\rm H_2)$ $\sim$ 5 $\times$ 10$^{4}$--5 $\times$ 10$^{6}$ cm$^{-3}$. The observed column densities of C$_{2}$H and c-C$_3$H$_{2}$ are reproduced reasonably well by our PDR models. This supports the idea that in high-UV flux PDRs, gas phase chemistry is sufficient to explain hydrocarbon abundances.

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.

Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive, the abundance of these species is likely linked to the evolutionary stage of the source. Aims: We investigate how the abundances of o-H$_2$D$^+$ and N$_2$D$^+$ vary with evolution in high-mass clumps. Methods: We observed with APEX the ground-state transitions of o-H$_2$D$^+$ near 372 GHz, and N$_2$D$^+$(3-2) near 231 GHz for three massive clumps in different evolutionary stages. The sources were selected within the G351.77-0.51 complex to minimise the variation of initial chemical conditions, and to remove distance effects. We modelled their dust continuum emission to estimate their physical properties, and also modelled their spectra under the assumption of local thermodynamic equilibrium to calculate beam-averaged abundances. Results: We find an anticorrelation between the abundance of o-H$_2$D$^+$ and that of N$_2$D$^+$, with the former decreasing and the latter increasing with evolution. With the new observations we are also able to provide a qualitative upper limit to the age of the youngest clump of about 10$^5$ yr, comparable to its current free-fall time. Conclusions: We can explain the evolution of the two tracers with simple considerations on the chemical formation paths, depletion of heavy elements, and evaporation from the grains. We therefore propose that the joint observation and the relative abundance of o-H$_2$D$^+$ and N$_2$D$^+$ can act as an efficient tracer of the evolutionary stages of the star-formation process.

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.

Stretching a length reaching 10 pc projected in the plane of sky, the radio jet associated with Herbig-Haro objects 80 and 81 (HH 80-81) is known as the largest and best collimated protostellar jet in our Galaxy. The nature of the molecular outflow associated with this extraordinary jet remains an unsolved question which is of great interests to our understanding of the relationship between jets and outflows in high-mass star formation. Here we present Atacama Pathfinder EXperiment CO(6-5) and (7-6), James Clerk Maxwell Telescope CO(3-2), Caltech Submillimeter Observatory CO(2-1), and Submillimeter Array CO and $^{13}$CO(2-1) mapping observations of the outflow. We report on the detection of a two-component outflow consisting of a collimated component along the jet path and a wide-angle component with an opening angle of about $30^{\circ}$. The gas velocity structure suggests that each of the two components traces part of a primary wind. From LVG calculations of the CO lines, the outflowing gas has a temperature around 88 K, indicating that the gas is being heated by shocks. Based on the CO(6-5) data, the outflow mass is estimated to be a few $M_{\odot}$, which is dominated by the wide-angle component. A comparison between the HH 80-81 outflow and other well shaped massive outflows suggests that the opening angle of massive outflows continues to increase over time. Therefore, the mass loss process in the formation of early-B stars seems to be similar to that in low-mass star formation, except that a jet component would disappear as the central source evolves to an ultracompact HII region.

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.

The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR). THOR is a Karl G. Jansky Very Large Array large program of atomic, molecular and ionized gas in the range 15{\deg}$\leq$l$\leq$67{\deg} and |b|$\leq$1{\deg}. It is the highest-resolution unbiased OH absorption survey to date towards this region. We combine the derived optical depths with literature 13CO(1-0) and HI observations to determine the OH abundance. We detect absorption in the 1665 and 1667 MHz transitions for continuum sources stronger than $F_{\rm cont}\geq$0.1 Jy/beam. OH absorption is found against $\sim$15% of these continuum sources with increasing fractions for stronger sources. Most of the absorption is associated with Galactic HII regions. We find OH and 13CO gas to have similar kinematic properties. The OH abundance decreases with increasing hydrogen column density. The OH abundance with respect to the total hydrogen nuclei column density (atomic and molecular phase) is in agreement with a constant abundance for $A_V$ < 10-20. Towards the lowest column densities, we find sources that exhibit OH absorption but no 13CO emission, indicating that OH is a well suited tracer of the low column density molecular gas. We present spatially resolved OH absorption towards W43. The unbiased nature of the THOR survey opens a new window onto the gas properties of the ISM. The characterization of the OH abundance over a large range of hydrogen gas column densities contributes to the understanding of OH as a molecular gas tracer and provides a starting point for future investigations.

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.

We present Atacama Pathfinder EXperiment (APEX) observations toward the massive star-forming region G240.31+0.07 in the CO J = 3--2, 6--5, and 7--6 lines. We detect a parsec-sized, bipolar, and high velocity outflow in all the lines, which allow us, in combination with the existing CO J = 2--1 data, to perform a multi-line analysis of physical conditions of the outflowing gas. The CO 7--6/6--5, 6--5/3--2, and 6--5/2--1 ratios are found to be nearly constant over a velocity range of $\sim$5--25 km s$^{-1}$ for both blueshifted and redshifted lobes. We carry out rotation diagram and large velocity gradient (LVG) calculations of the four lines, and find that the outflow is approximately isothermal with a gas temperature of $\sim$50 K, and that the the CO column density clearly decreases with the outflow velocity. If the CO abundance and the velocity gradient do not vary much, the decreasing CO column density indicates a decline in the outflow gas density with velocity. By comparing with theoretical models of outflow driving mechanisms, our observations and calculations suggest that the massive outflow in G240.31+0.07 is being driven by a wide-angle wind and further support a disk mediated accretion at play for the formation of the central high-mass star.

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.

Messier 8 is one of the brightest HII regions in the sky. Using SOFIA, APEX 12 m and IRAM 30 m telescopes, we have performed a comprehensive imaging survey of the emission from the [CII], [CI] and CO isotopologues within 1.3 $\times$ 1.3 pc around the dominant Herschel 36 (Her 36) stellar system. To further explore the morphology of the region we compare archival infrared, optical and radio images of the nebula with our newly obtained [CII] and CO data. We perform a quantitative analysis, using both LTE and non-LTE methods to determine the abundances of some of the observed species as well as kinetic temperatures and volume densities. Bright CO, [CII] and [CI] emission has been found towards the HII region and the photodissociation region (PDR) in M8. Our analysis places the bulk of the molecular material in the background of the nebulosity illuminated by the bright stellar systems Her 36 and 9 Sagitarii. Since the emission from all observed atomic and molecular tracers peaks at or close to the position of Her 36, we conclude that the star is still physically close to its natal dense cloud core and heats it. A veil of warm gas moves away from Her 36 toward the Sun and its associated dust contributes to the foreground extinction in the region. One of the most prominent star forming regions in M8, the Hourglass Nebula, is particularly bright due to cracks in this veil close to Her 36. By using radiative transfer modeling of different transitions of CO isotopologues, we obtain H$_{2}$ densities ranging from $\sim 10^{4} - 10^{6}$ cm$^{-3}$ and kinetic temperatures of 100 $-$ 150 K in the bright PDR caused by Her 36.

Context. HCN is a major constituent of the circumstellar envelopes of carbon-rich evolved stars, and rotational lines from within its vibrationally excited states probe parts of these regions closest to the stellar surface. A number of such lines are known to show maser action. Historically, in one of them, the 177 GHz $J=2\rightarrow1$ line in the $l$-doubled bending mode has been found to show relatively strong maser action, with results only published for a single object, the archetypical high-mass loss asymptotic giant branch (AGB) star IRC+10216. Aims. To examine how common 177 GHz HCN maser emission is, we conducted an exploratory survey for this line toward a select sample of carbon-rich asymptotic giant branch stars that are observable from the southern hemisphere. Methods. We used the Atacama Pathfinder Experiment 12 meter submillimeter Telescope (APEX) equipped with a new receiver to simultaneously observe three $J=2\rightarrow1$ HCN rotational transitions, the $(0,1^{{1}_{\rm c}},0)$ and $(0,1^{{1}_{\rm d}},0)$ $l$-doublet components, and the line from the (0,0,0) ground state. Results. The $(0,1^{{1}_{\rm c}},0)$ maser line is detected toward 11 of 13 observed sources, which all show emission in the (0,0,0) transition. In most of the sources, the peak intensity of the $(0,1^{{1}_{\rm c}},0)$ line rivals that of the (0,0,0) line; in two sources, it is even stronger. Except for the object with the highest mass-loss rate, IRC+10216, the $(0,1^{{1}_{\rm c}},0)$ line covers a smaller velocity range than the (0,0,0) line. Conclusions. Maser emission in the 177 GHz $J=2\rightarrow1$ $(0,1^{{1}_{\rm c}},0)$ line of HCN appears to be common in carbon-rich AGB stars. (Abbreviated)

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.

(Abridged) Context: The assumption of a gas-to-dust mass ratio (\gamma) is a common approach to estimate the basic properties of molecular clouds, such as total mass and column density of molecular hydrogen, from (sub)mm continuum observations of the dust. In the Milky Way a single value is used at all galactocentric radii, independently of the observed metallicity gradients. Both models and extragalactic observations suggest that this quantity increases for decreasing metallicity Z, typical of the outer regions in disks, where fewer heavy elements are available to form dust grains. Aims: We aim to investigate the variation of the gas-to-dust ratio as a function of galactocentric radius and metallicity, to allow a more accurate characterisation of the quantity of molecular gas across the galactic disk, as derived from observations of the dust. Methods: Observations of the optically thin C\$^{18}\$O (2-1) transition were obtained with the APEX telescope for a sample of 23 massive and dense star-forming regions in the far outer Galaxy (galactocentric distance greater than 14 kpc). From the modelling of this line and of the spectral energy distribution of the selected clumps we computed the gas-to-dust ratio and compared it to that of well-studied sources from the ATLASGAL TOP100 sample in the inner galactic disk. Results: The gradient in gas-to-dust ratio is found to be 0.087 dex/kpc (or equivalently \gamma\ proportional to Z\$^{-1.4}\$). The dust-to-metal ratio, decreases with galactocentric radius, which is the most common situation also for external late-type galaxies. This suggests that grain growth dominates over destruction. The predicted gas-to-dust ratio is in excellent agreement with the estimates in Magellanic clouds, for the appropriate value of Z.

The circumstellar ammonia (NH$_3$) chemistry in evolved stars is poorly understood. Previous observations and modelling showed that NH$_3$ abundance in oxygen-rich stars is several orders of magnitude above that predicted by equilibrium chemistry. In this article, we characterise the spatial distribution and excitation of NH$_3$ in the O-rich circumstellar envelopes (CSEs) of four diverse targets: IK Tau, VY CMa, OH 231.8+4.2, and IRC +10420 with multi-wavelength observations. We observed the 1.3-cm inversion line emission with the Very Large Array (VLA) and submillimetre rotational line emission with the Heterodyne Instrument for the Far-Infrared (HIFI) aboard Herschel from all four targets. For IK Tau and VY CMa, we observed the rovibrational absorption lines in the $\nu_2$ band near 10.5 $\mu$m with the Texas Echelon Cross Echelle Spectrograph (TEXES) at the NASA Infrared Telescope Facility (IRTF). We also attempted to search for the rotational transition within the $v_2=1$ state near 2 mm with the IRAM 30m Telescope towards IK Tau. Non-LTE radiative transfer modelling, including radiative pumping to the vibrational state, was carried out to derive the radial distribution of NH$_3$ in these CSEs. Our modelling shows that the NH$_3$ abundance relative to molecular hydrogen is generally of the order of $10^{-7}$, which is a few times lower than previous estimates that were made without considering radiative pumping and is at least 10 times higher than that in the C-rich CSE of IRC +10216. Incidentally, we also derived a new period of IK Tau from its $V$-band light curve. NH$_3$ is again detected in very high abundance in O-rich CSEs. Its emission mainly arises from localised spatial-kinematic structures that are probably denser than the ambient gas. Circumstellar shocks in the accelerated wind may contribute to the production of NH$_3$. (Abridged abstract)

We mapped the kinetic temperature structure of the Orion molecular cloud 1 with para-H2CO(303-202, 322-221, and 321-220) using the APEX 12m telescope. This is compared with the temperatures derived from the ratio of the NH3(2,2)/(1,1) inversion lines and the dust emission. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured averaged line ratios of para-H2CO 322-221/303-202 and 321-220/303-202. The gas kinetic temperatures derived from the para-H2CO line ratios are warm, ranging from 30 to >200 K with an average of 62 K at a spatial density of 10$^5$ cm$^{-3}$. These temperatures are higher than those obtained from NH3(2,2)/(1,1) and CH3CCH(6-5) in the OMC-1 region. The gas kinetic temperatures derived from para-H2CO agree with those obtained from warm dust components measured in the mid infrared (MIR), which indicates that the para-H2CO(3-2) ratios trace dense and warm gas. The cold dust components measured in the far infrared (FIR) are consistent with those measured with NH3(2,2)/(1,1) and the CH3CCH(6-5) line series. With dust at MIR wavelengths and para-H2CO(3-2) on one side and dust at FIR wavelengths, NH3(2,2)/(1,1), and CH3CCH(6-5) on the other, dust and gas temperatures appear to be equivalent in the dense gas of the OMC-1 region, but provide a bimodal distribution, one more directly related to star formation than the other. The non-thermal velocity dispersions of para-H2CO are positively correlated with the gas kinetic temperatures in regions of strong non-thermal motion (Mach number >2.5) of the OMC-1, implying that the higher temperature traced by para-H2CO is related to turbulence on a 0.06 pc scale. Combining the temperature measurements with para-H2CO and NH3(2,2)/(1,1) line ratios, we find direct evidence for the dense gas along the northern part of the OMC-1 10 km s$^{-1}$ filament heated by radiation from the central Orion nebula.

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.

(abridged) Our aim is to determine the radial abundance profile of SiO and HCN throughout the stellar outflow of R Dor, an oxygen-rich AGB star with a low mass-loss rate. We have analysed molecular transitions of CO, SiO, and HCN measured with the APEX telescope and all three instruments on the Herschel Space Observatory, together with literature data. Photometric data and the infrared spectrum measured by ISO-SWS were used to constrain the dust component of the outflow. Using both continuum and line radiative transfer methods, a physical envelope model of both gas and dust was established. We have performed an analysis of the SiO and HCN molecular transitions in order to calculate their abundances. We have obtained an envelope model that describes the dust and the gas in the outflow, and determined the abundance of SiO and HCN throughout the region of the outflow probed by our molecular data. For SiO, we find that the initial abundance lies between $5.5 \times 10^{-5}$ and $6.0 \times 10^{-5}$ w.r.t. H$_2$. The abundance profile is constant up to $60\ \pm 10\ R_*$, after which it declines following a Gaussian profile with an $e$-folding radius of $3.5 \pm 0.5 \times 10^{13}$ cm. For HCN, we find an initial abundance of $5.0 \times 10^{-7}$ w.r.t. H$_2$. The Gaussian profile that describes the decline starts at the stellar surface and has an $e$-folding radius $r_e$ of $1.85 \pm 0.05 \times 10^{15}$ cm. We cannot to unambiguously identify the mechanism by which SiO is destroyed at $60\ \pm 10\ R_*$. The initial abundances found are larger than previously determined (except for one previous study on SiO), which might be due to the inclusion of higher-$J$ transitions. The difference in abundance for SiO and HCN compared to high mass-loss rate Mira star IK Tau might be due to different pulsation characteristics of the central star and/or a difference in dust condensation physics.

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.

The formation of massive stars is still not well understood. Accumulating a large amount of mass infalling within a single entity in spite of radiation pressure is possible if, among several other conditions, enough thermal energy is released. Despite numerous water line observations, with the Herschel Space Observatory, in most of the sources observations were not able to trace the emission from the hot core around the newly forming protostellar object. We want to probe the physical conditions and water abundance in the inner layers of the host protostellar object NGC7538-IRS1 using a highly excited H2O line. Water maser models predict that several THz water masers should be detectable in these objects. We present SOFIA observations of the o-H2O 8(2,7)-7(3,4) line at 1296.41106 GHz and a 6(1,6)-5(2,3) 22 GHz e-MERLIN map of the region (first-ever 22 GHz images made after the e-MERLIN upgrade). In order to be able to constrain the nature of the emission - thermal or maser - we use near-simultaneous observations of the 22 GHz water maser performed with the Effelsberg radiotelescope and e-MERLIN. A thermal water model using the RATRAN radiative transfer code is presented based on HIFI pointed observations. Molecular water abundances are derived for the hot core. The H2O 8(2,7)- 7(3,4) line is detected toward NGC7538-IRS1 with one feature at the source velocity (-57.7 km/s) and another one at -48.4 km/s. We propose that the emission at the source velocity is consistent with thermal excitation and is excited in the innermost part of the IRS1a massive protostellar object's closest circumstellar environment. The other emission is very likely the first detection of a water THz maser line, pumped by shocks due to IRS1b outflow, in a star-forming region. Assuming thermal excitation of the THz line, the water abundance in NGC7538-IRS1's hot core is estimated to be 5.2x10^{-5} with respect to H2.

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

We aim at characterizing the large-scale distribution of H2O in G327.3-0.6, a massive star-forming region made of individual objects in different evolutionary phases. We investigate variations of H2O abundance as function of evolution. We present Herschel continuum maps at 89 and 179 $\mu$m of the whole region and an APEX map at 350 {\mu}m of the IRDC. New spectral HIFI maps toward the IRDC region covering low-energy H2O lines at 987 and 1113 GHz are also presented and combined with HIFI pointed observations of the G327 hot core. We infer the physical properties of the gas through optical depth analysis and radiative transfer modeling. The continuum emission at 89 and 179 {\mu}m follows the thermal continuum emission at longer wavelengths, with a peak at the position of the hot core, a secondary peak in the Hii region, and an arch-like layer of hot gas west of the Hii region. The same morphology is observed in the 1113 GHz line, in absorption toward all dust condensations. Optical depths of ~80 and 15 are estimated and correspond to column densities of 10^15 and 2 10^14 cm-2, for the hot core and IRDC position. These values indicate an H2O to H2 ratio of 3 10^-8 toward the hot core; the abundance of H2O does not change along the IRDC with values of some 10^-8. Infall (over ~ 20") is detected toward the hot core position with a rate of 1-1.3 10^-2 M_sun /yr, high enough to overcome the radiation pressure due to the stellar luminosity. The source structure of the hot core region is complex, with a cold outer gas envelope in expansion, situated between the outflow and the observer, extending over 0.32 pc. The outflow is seen face-on and centered away from the hot core. The distribution of H2O along the IRDC is roughly constant with an abundance peak in the more evolved object. These water abundances are in agreement with previous studies in other massive objects and chemical models.

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

Recent submillimeter and far-infrared wavelength observations of absorption in the rotational ground-state lines of various simple molecules against distant Galactic continuum sources have opened the possibility of studying the chemistry of diffuse molecular clouds throughout the Milky Way. In order to calculate abundances, the column densities of molecular and atomic hydrogen, HI, must be known. We aim at determining the atomic hydrogen column densities for diffuse clouds located on the sight lines toward a sample of prominent high-mass star-forming regions that were intensely studied with the HIFI instrument onboard Herschel. Based on Jansky Very Large Array data, we employ the 21 cm HI absorption-line technique to construct profiles of the HI opacity versus radial velocity toward our target sources. These profiles are combined with lower resolution archival data of extended HI emission to calculate the HI column densities of the individual clouds along the sight lines. We employ Bayesian inference to estimate the uncertainties of the derived quantities. Our study delivers reliable estimates of the atomic hydrogen column density for a large number of diffuse molecular clouds at various Galactocentric distances. Together with column densities of molecular hydrogen derived from its surrogates observed with HIFI, the measurements can be used to characterize the clouds and investigate the dependence of their chemistry on the molecular fraction, for example.

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.

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

The dynamical processes that control star formation in molecular clouds are not well understood, and in particular, it is unclear if rotation plays a major role in cloud evolution. We investigate the importance of rotation in cloud evolution by studying the kinematic structure of a spiral-shaped Galactic molecular cloud G052.24+00.74. The cloud belongs to a large filament, and is stretching over ~ 100 pc above the Galactic disk midplane. The spiral-shaped morphology of the cloud suggests that the cloud is rotating. We have analysed the kinematic structure of the cloud, and study the fragmentation and star formation. We find that the cloud exhibits a regular velocity pattern along west-east direction - a velocity shift of ~ 10 km/s at a scale of ~ 30 pc. The kinematic structure of the cloud can be reasonably explained by a model that assumes rotational support. Similarly to our Galaxy, the cloud rotates with a prograde motion. We use the formalism of Toomre (1964) to study the cloud's stability, and find that it is unstable and should fragment. The separation of clumps can be consistently reproduced assuming gravitational instability, suggesting that fragmentation is determined by the interplay between rotation and gravity. Star formation occurs in massive, gravitational bound clumps.

We present a high sensitivity radio continuum survey at 6 and 1.3$\,$cm using the Karl G. Jansky Very Large Array towards a sample of 58 high-mass star forming regions. Our sample was chosen from dust clumps within infrared dark clouds with and without IR sources (CMC-IRs, CMCs, respectively), and hot molecular cores (HMCs), with no previous, or relatively weak radio continuum detection at the $1\,$mJy level. Due to the improvement in the continuum sensitivity of the VLA, this survey achieved map rms levels of $\sim$ 3-10 $\mu$Jy beam$^{-1}$ at sub-arcsecond angular resolution. We extracted 70 centimeter continuum sources associated with 1.2$\,$mm dust clumps. Most sources are weak, compact, and are prime candidates for high-mass protostars. Detection rates of radio sources associated with the mm dust clumps for CMCs, CMC-IRs and HMCs are 6$\%$, 53$\%$ and 100$\%$, respectively. This result is consistent with increasing high-mass star formation activity from CMCs to HMCs. The radio sources located within HMCs and CMC-IRs occur close to the dust clump centers with a median offset from it of 12,000$\,$AU and 4,000$\,$AU, respectively. We calculated 5 - 25$\,$GHz spectral indices using power law fits and obtain a median value of 0.5 (i.e., flux increasing with frequency), suggestive of thermal emission from ionized jets. In this paper we describe the sample, observations, and detections. The analysis and discussion will be presented in Paper II.

We carried out a millimeter wavelength line survey between 79 and 356 GHz with the IRAM-30m telescope to investigate the physical and chemical properties of the molecular envelope of the oxygen-rich AGB star IK Tau. We analysed the molecular lines detected using the population diagram technique to derive rotational temperatures and column densities. Additionally, we conducted a radiative transfer analysis of the SO$_2$ lines detected. For the first time in this source we detected rotational lines in the ground vibrational state of HCO$^+$, NS, NO, and H$_2$CO, as well as several isotopologues of molecules previously identified. We also detected several rotational lines in vibrationally excited states of SiS and SiO isotopologues, and rotational lines of H$_2$O ($\nu_{\mathrm{2}}$=2). We have also increased the number of rotational lines detected of molecules that were previously identified, enabling a detailed study of the molecular abundances and excitation temperatures. IK Tau displays a rich chemistry for an oxygen-rich circumstellar envelope. We highlight the detection of NS and H$_2$CO with fractional abundances of f(NS)$\sim$10$^{-8}$ and f(H$_2$CO)$\sim$[10$^{-7}$--10$^{-8}$]. Most of the molecules display rotational temperatures between 15 and 40 K. NaCl and SiS isotopologues display rotational temperatures higher than the average (~65 K). In the case of SO$_2$ a warm component with $T_{\mathrm{rot}}$$\sim$290 K is also detected. This SO$_2$ warm component is probably arising from the inner regions of the envelope (at $\sim$8$R_{*}$) where SO$_2$ has a fractional abundance of f(SO$_2$)$\sim$10$^{-6}$. This result should be considered for future investigation of the main formation channels of this, and other, parent species in the inner winds of O-rich AGB stars, which at present are not well reproduced by current chemistry models.

Massive stars have a major influence on their environment yet their formation is difficult to study. W43 is a highly luminous galactic massive star forming region at a distance of 5.5 kpc and the MM1 part hosts a very massive dense core (1000 M$_{\odot}$ within 0.05 pc). We present new Herschel HIFI maps of the W43 MM1 region covering the main low-energy water lines at 557, 987, and 1113 GHz, their H$_2^{18}$O counterparts, and other lines such as $^{13}$CO(10-9) and C$^{18}$O(9-8) which trace warm gas. These water lines are, with the exception of line wings, observed in absorption. Herschel SPIRE and JCMT 450 $\mu$m data have been used to make a model of the continuum emission at these wavelengths. Analysis of the maps, and in particular the optical depth maps of each line and feature, shows that a velocity gradient, possibly due to rotation, is present in both the envelope and the protostellar core. Velocities increase in both components from SW to NE, following the general source orientation. While the H$_2$O lines trace essentially the cool envelope, we show that the envelope cannot account for the H$_2^{18}$O absorption, which traces motions close to the protostar. The core has rapid infall, 2.9 kms, as manifested by the H$_2^{18}$O absorption features which are systematically red-shifted with respect to the $^{13}$CO(10-9) emission line which also traces the inner material with the same angular resolution. Some H$_2^{18}$O absorption is detected outside the central core and thus outside the regions expected to be above 100 K - we attribute this to warm gas associated with the other massive dense cores in W43 MM1. Using the maps to identify absorption from cool gas on large scales, we subtract this component to model spectra for the inner envelope. Modeling the new spectra results in a lower water abundance, decreased from $8 10^{-8}$ to $8 10^{-9}$ , with no change in infall rate.

We present the first detection of a jet in the far-IR [O I] lines from an intermediate mass protostar. We have carried out a Herschel/PACS spectral mapping study in the [O I] lines of OMC-2 FIR 3 and FIR 4, two of the most luminous protostars in Orion outside of the Orion Nebula. The spatial morphology of the fine structure line emission reveals the presence of an extended photodissociation region (PDR) and a narrow, but intense jet connecting the two protostars. The jet seen in [O I] emission is spatially aligned with the Spitzer/IRAC 4.5 micron jet and the CO (6-5) molecular outflow centered on FIR 3. The mass loss rate derived from the total [O I] 63 micron line luminosity of the jet is 7.7 x 10^-6 M_sun/yr, more than an order of magnitude higher than that measured for typical low mass class 0 protostars. The implied accretion luminosity is significantly higher than the observed bolometric luminosity of FIR 4, indicating that the [O I] jet is unlikely to be associated with FIR 4. We argue that the peak line emission seen toward FIR 4 originates in the terminal shock produced by the jet driven by FIR 3. The higher mass-loss rate that we find for FIR 3 is consistent with the idea that intermediate mass protostars drive more powerful jets than their low-mass counterparts. Our results also call into question the nature of FIR 4.

We present CARMA CO (J=1-0) observations and Herschel PACS spectroscopy, characterizing the outflow properties toward extremely young and deeply embedded protostars in the Orion molecular clouds. The sample comprises a subset of the Orion protostars known as the PACS Bright Red Sources (PBRS) (Stutz et al. 2013). We observed 14 PBRS with CARMA and 8 of these 14 with Herschel, acquiring full spectral scans from 55 micron to 200 micron. Outflows are detected in CO (J=1-0) from 8 of 14 PBRS, with two additional tentative detections; outflows are also detected from the outbursting protostar HOPS 223 (V2775 Ori) and the Class I protostar HOPS 68. The outflows have a range of morphologies, some are spatially compact, <10000 AU in extent, while others extend beyond the primary beam. The outflow velocities and morphologies are consistent with being dominated by intermediate inclination angles (80 deg > i > 20 deg). This confirms the interpretation of the very red 24 micron to 70 micron colors of the PBRS as a signpost of high envelope densities, with only one (possibly two) cases of the red colors resulting from edge-on inclinations. We detect high-J (J_up > 13) CO lines and/or H_2O lines from 5 of 8 PBRS and only for those with detected CO outflows. The far-infrared CO rotation temperatures of the detected PBRS are marginally colder (~230 K) than those observed for most protostars (~300 K), and only one of these 5 PBRS has detected [OI] 63 micron emission. The high envelope densities could be obscuring some [OI] emission and cause a ~20 K reduction to the CO rotation temperatures.

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

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $\leq0.1\,M_{\odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{\rm gas}=(0.1-1.5)\times 10^{-4}$. HD 100546 itself is a $\lambda\,$Bo\"{o}tis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{\rm gas}=(0.2-5.0)\times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{\rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems.

The formation mechanism of brown dwarfs (BDs) is one of the long-standing problems in star formation because the typical Jeans mass in molecular clouds is too large to form these substellar objects. To answer this question, it is crucial to study a BD at the embedded phase. IRAS 16253-2429 is classified as a very low luminosity object (VeLLO) with internal luminosity 0.1 Lsun. VeLLOs are believed to be very low-mass protostars or even proto-BDs. We observed the jet/outflow driven by IRAS 16253-2429 in CO (2-1), (6-5), and (7-6) using the IRAM 30 m and APEX telescopes and the SMA in order to study its dynamical features and physical properties. Our SMA map reveals two protostellar jets, indicating the existence of a proto-binary system as implied by the precessing jet detected in H2 emission. We detect a wiggling pattern in the position-velocity diagrams along the jet axes, which is likely due to the binary orbital motion. Based on this, we derive the current mass of the binary as ~0.032 Msun. Given the low envelope mass, IRAS 16253-2429 will form a binary that probably consist of one or two BDs. Furthermore, we found that the outflow force as well as the mass accretion rate are very low based on the multi-transition CO observations, which suggests that the final masses of the binary components are at the stellar/substellar boundary. Since IRAS 16253 is located in an isolated environment, we suggest that BDs can form through fragmentation and collapse like low-mass stars.

Aims. The aim of the study is to constrain the physics of gamma-ray bursts (GRBs) by analysing the multi-wavelength afterglow data set of GRB 121024A that covers the full range from radio to X-rays. Methods. Using multi-epoch broad-band observations of the GRB 121024A afterglow, we measured the three characteristic break frequencies of the synchrotron spectrum. We used six epochs of combined XRT and GROND data to constrain the temporal slopes, the dust extinction, the X-ray absorption, and the spectral slope with high accuracy. Two more epochs of combined data from XRT, GROND, APEX, CARMA, and EVLA were used to set constraints on the break frequencies and therefore on the micro-physical and dynamical parameters. Results. The XRT and GROND light curves show a simultaneous and achromatic break at around 49 ks. As a result, the crossing of the synchrotron cooling break is no suitable explanation for the break in the light curve. The multi wavelength data allow us to test two plausible scenarios explaining the break: a jet break, and the end of energy injection. The jet-break scenario requires a hard electron spectrum, a very low cooling break frequency, and a non-spreading jet. The energy injection avoids these problems, but requires $\epsilon_e > 1 (k = 2)$, spherical outflow, and $\epsilon_B < 10^{-9}$. Conclusions. In light of the extreme microphysical parameters required by the energy-injection model, we favour a jet-break scenario where $\nu_m < \nu_{sa}$ to explain the observations. This scenario gives physically meaningful microphysical parameters, and it also naturally explains the reported detection of linear and circular polarisation.

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

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

Gravity plays important roles in the evolution of molecular clouds. We present an acceleration mapping method to estimate the acceleration induced by gravitational interactions in molecular clouds based on observational data. We find that the geometry of a region has a significant impact on the behavior of gravity. In the Pipe nebula which can be approximated as a gas filament, we find that gravitational acceleration can effectively compress the end of this filament, which may have triggered star formation. We identify this as the "gravitational focusing" effect proposed by Burkert & Hartman (2004). In the sheet-like IC348-B3 region, gravity can lead to collapse at its edge, while in the centrally condensed NGC1333 cluster-forming region gravity can drive accretion towards the center. In general, gravitational acceleration tends to be enhanced in the localized regions around the ends of the filaments and the edges of sheet-like structures. Neglecting magnetic fields, these "gravitational focusing" and "edge collapse" effects can promote the formation of dense gas in a timescale that is much shorter than the global dynamical time. Since the interstellar medium is in general structured, these edge effects should be prevalent.

The prototypical Mira variable, $o$ Cet (Mira), has been observed as a Science Verification target in the 2014 ALMA Long Baseline Campaign with a longest baseline of 15 km. ALMA clearly resolves the images of the continuum and molecular line emission/absorption at an angular resolution of ~30 mas at 220 GHz. We image the data of the $^{28}$SiO v=0, 2 $J$=5-4 and H$_2$O $\nu_2$=1 $J(K_a,K_c)$=5(5,0)-6(4,3) transitions and extract spectra from various lines-of-sight towards Mira's extended atmosphere. In the course of imaging, we encountered ambiguities in the resulting images and spectra that appear to be related to the performance of the CLEAN algorithm. We resolve Mira's millimetre continuum emission and our data are consistent with a radio photosphere with a brightness temperature of 2611$\pm$51 K, in agreement with recent results obtained with the VLA. We do not confirm the existence of a compact region (<5 mas) of enhanced brightness. We derive the gas density, kinetic temperature, molecular abundance and outflow/infall velocities in Mira's extended atmosphere by modelling the SiO and H$_2$O lines. We find that SiO-bearing gas starts to deplete beyond 4$R_\star$ and at a kinetic temperature of $\lesssim$600 K. The inner dust shells are probably composed of grain types other than pure silicates. During this observation, Mira's atmosphere generally exhibited infall motion, with a shock front of velocity $\lesssim$12 km/s outside the radio photosphere. The structures predicted by the hydrodynamical model CODEX can reproduce the observed spectra in astonishing detail; while some other models fail when confronted with the new data. Combined with radiative transfer modelling, ALMA successfully demonstrates the ability to reveal the physical conditions of the extended atmospheres and inner winds of AGB stars in unprecedented detail. (Abbreviated abstract)

Context. The chemistry of the diffuse interstellar medium rests upon three pillars: exothermic ion-neutral reactions (" cold chemistry "), endothermic neutral-neutral reactions with significant activation barriers (" warm chemistry "), and reactions on the surfaces of dust grains. While warm chemistry becomes important in the shocks associated with turbulent dissipation regions, the main path for the formation of interstellar OH and H2O is that of cold chemistry. Aims. The aim of this study is to observationally confirm the association of atomic oxygen with both atomic and molecular gas phases, and to understand the measured abundances of OH and OH + as a function of the available reservoir of H2. Methods. We obtained absorption spectra of the ground states of OH, OH+ and OI with high-velocity resolution, with GREAT on-board SOFIA, and with the THz receiver at the APEX. We analyzed them along with ancillary spectra of HF and CH from HIFI. To deconvolve them from the hyperfine structure and to separate the blend that is due to various velocity components on the sightline, we fit model spectra consisting of an appropriate number of Gaussian profiles using a method combining simulated annealing with downhill simplex minimization. Together with HF and/or CH as a surrogate for H2, and HI $\lambda$21 cm data, the molecular hydrogen fraction f^N\_H2 = N(H 2)/(N(H) + 2N(H 2)) can be determined. We then investigated abundance ratios as a function of f^N\_H2. Results. The column density of OI is correlated at a high significance with the amount of available molecular and atomic hydrogen, with an atomic oxygen abundance of $3 \times 10 ^{-4}$ relative to H nuclei. While the velocities of the absorption features of OH and OH+ are loosely correlated and reflect the spiral arm crossings on the sightline, upon closer inspection they display an anticorrespondence. The arm-to-interarm density contrast is found to be higher in OH than in OH+. While both species can coexist, with a higher abundance in OH than in OH+, the latter is found less frequently in absence of OH than the other way around, which is a direct consequence of the rapid destruction of OH+ by dissociative recombination when not enough H2 is available. This conjecture has been substantiated by a comparison between the OH/OH+ ratio with f^N\_H2, showing a clear correlation. The hydrogen abstraction reaction chain OH+ (H2,H) H2O+ (H2,H)H3O+ is confirmed as the pathway for the production of OH and H 2 O. Our estimate of the branching ratio of the dissociative recombination of H3O+ to OH and H2O is confined within the interval of 84 to 91%, which matches laboratory measurements (74 to 83%). -- A correlation between the linewidths and column densities of OH+ features is found to be significant with a false-alarm probability below 5%. Such a correlation is predicted by models of interstellar MHD turbulence. For OH the same correlation is found to be insignificant because there are more narrow absorption features. Conclusions. While it is difficult to assess the contributions of warm neutral-neutral chemistry to the observed abundances, it seems fair to conclude that the predictions of cold ion-neutral chemistry match the abundance patterns we observed.

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.

We present Herschel/HIFI observations (WISH KP) of 14 water lines in a small sample of galactic massive protostellar objects: NGC6334I(N), DR21(OH), IRAS16272-4837, and IRAS05358+3543. We analyze the gas dynamics from the line profiles. Through modeling of the observations using RATRAN, we estimate outflow, infall, turbulent velocities, molecular abundances, and investigate any correlation with the evolutionary status of each source. The molecular line profiles exhibit a broad component coming from the shocks along the cavity walls associated with the protostars, and an infalling (or expansion for IRAS05358+3543) and passively heated envelope component, with highly supersonic turbulence likely increasing with the distance from the center. Accretion rates between 6.3 10^{-5} and 5.6 10^{-4} \msun yr^{-1} are derived from the infall observed in three of our sources. The outer water abundance is estimated to be at the typical value of a few 10^{-8} while the inner abundance varies from 1.7 10^{-6} to 1.4 10^{-4} with respect to H2 depending on the source. We confirm that regions of massive star formation are highly turbulent and that the turbulence likely increases in the envelope with the distance to the star. The inner abundances are lower than the expected 10^{-4} perhaps because our observed lines do not probe deep enough into the inner envelope, or because photodissociation through protostellar UV photons is more efficient than expected. We show that the higher the infall/expansion velocity in the protostellar envelope, the higher is the inner abundance, maybe indicating that larger infall/expansion velocities generate shocks that will sputter water from the ice mantles of dust grains in the inner region. High-velocity water must be formed in the gas-phase from shocked material.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

OB type stars have strong ionizing radiation, and drive energetic winds. The ultraviolet (UV) radiation from ionizing stars may heat dust and ionize gas to sweep up an expanding bubble shell. This shell may be the result of feedback leading to a new generation of stars. N131 is an infrared dust bubble residing in a molecular filament. We study the formation and fragmentation of this bubble with multi-wavelength dust and gas observations. Towards the bubble N131, we analyzed archival multi-wavelength observations including 3.6, 4.5, 5.8, 8.0, 24, 70, 160, 250, 350, 500 $\mu$m, 1.1 mm, and 21 cm. In addition, we performed new observations of CO (2-1), CO (1-0), and $^{13}$CO (1-0) with the IRAM 30-m telescope. Multi-wavelength dust and gas observations reveal a ringlike shell with compact fragments, two filamentary structures, and a secondary bubble N131-A. The bubble N131 is a rare object with a large hole at 24 $\mu$m and 21 cm in the direction of its center. The dust and gas clumps are compact and might have been compressed at the inner edge of the ringlike shell, while they are extended and might be pre-existing at the outer edge. The column density, excitation temperature, and velocity show a potentially hierarchical distribution from the inner to outer edge of the ringlike shell. We also detected the front and back sides of the secondary bubble N131-A in the direction of its center. The derived Lyman-continuum ionizing photon flux within N131-A is equivalent to an O9.5 star. Based on the above, we suggest that the bubble N131 might be triggered by the strong stellar winds from a group of massive stars inside the bubble. We propose a scenario in which the bubble N131 forms from the disruption of a gas filament by expansion of HII region, strong stellar winds, and fragments under self-gravity.

We have used the Herschel-HIFI instrument to observe both nuclear spin symmetries of amidogen (NH2) towards the high-mass star-forming regions W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4) and G34.3+0.1. The aim is to investigate the ratio of nuclear spin types, the ortho-to-para ratio (OPR), of NH2. The excited NH2 transitions are used to construct radiative transfer models of the hot cores and surrounding envelopes in order to investigate the excitation and possible emission of the ground state rotational transitions of ortho-NH2 N_(K_a,K_c} J=1_(1,1) 3/2 - 0_(0,0) 1/2 and para-NH2 2_(1,2) 5/2 - 1_(0,1) 3/2$ used in the OPR calculations. Our best estimate of the average OPR in the envelopes lie above the high temperature limit of three for W49N, specifically 3.5 with formal errors of \pm0.1, but for W31C, W51, and G34.3+0.1 we find lower values of 2.5\pm0.1, 2.7\pm0.1, and 2.3\pm0.1, respectively. Such low values are strictly forbidden in thermodynamical equilibrium since the OPR is expected to increase above three at low temperatures. In the translucent interstellar gas towards W31C, where the excitation effects are low, we find similar values between 2.2\pm0.2 and 2.9\pm0.2. In contrast, we find an OPR of 3.4\pm0.1 in the dense and cold filament connected to W51, and also two lower limits of >4.2 and >5.0 in two other translucent gas components towards W31C and W49N. At low temperatures (T \lesssim 50 K) the OPR of H2 is <10^-1, far lower than the terrestrial laboratory normal value of three. In such a "para-enriched H2" gas, our astrochemical models can reproduce the variations of the observed OPR, both below and above the thermodynamical equilibrium value, by considering nuclear-spin gas-phase chemistry. The models suggest that values below three arise in regions with temperatures >20-25 K, depending on time, and values above three at lower temperatures.

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.

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

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.

[Abridged] Our aim is to explore the gas dynamics and the accretion process in the early phase of high-mass star formation. The inward motion of molecular gas in the massive star forming region G34.26+0.15 is investigated by using high-resolution profiles of seven transitions of ammonia at THz frequencies observed with Herschel-HIFI. The shapes and intensities of these lines are interpreted in terms of radiative transfer models of a spherical, collapsing molecular envelope. An accelerated Lambda Iteration (ALI) method is used to compute the models. The seven ammonia lines show mixed absorption and emission with inverse P-Cygni-type profiles that suggest infall onto the central source. A trend toward absorption at increasingly higher velocities for higher excitation transitions is clearly seen in the line profiles. The $J = 3\leftarrow2$ lines show only very weak emission, so these absorption profiles can be used directly to analyze the inward motion of the gas. This is the first time a multitransitional study of spectrally resolved rotational ammonia lines has been used for this purpose. Broad emission is, in addition, mixed with the absorption in the $1_0-0_0$ ortho-NH$_3$ line, possibly tracing a molecular outflow from the star forming region. The best-fitting ALI model reproduces the continuum fluxes and line profiles, but slightly underpredicts the emission and absorption depth in the ground-state ortho line $1_0-0_0$. The derived ortho-to-para ratio is approximately 0.5 throughout the infalling cloud core similar to recent findings for translucent clouds in sight lines toward W31C and W49N. We find evidence of two gas components moving inwards toward the central region with constant velocities: 2.7 and 5.3 km$\,$s$^{-1}$, relative to the source systemic velocity. The inferred mass accretion rates derived are sufficient to overcome the expected radiation pressure from G34.26+0.15.

We report on observations of the hydroxyl radical (OH) within The H{\sc I}, OH Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l=29.2 to 31.5$^\circ$ and b=-1.0 to +1.0$^\circ$ and includes the high-mass star forming region W43. We identify 103 maser sites, including 72 with 1612\,MHz masers, 42 showing masers in either of the main line transitions at 1665 and 1667\,MHz and four showing 1720\,MHz masers. Most maser sites with either main-line or 1720\,MHz emission are associated with star formation, whereas most of the 1612\,MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We also find diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the well known star forming complex W43 Main.

The fine-structure line of [OI] at 63micron is an important diagnostic tool in different fields of astrophysics. However, our knowledge of this line relies on observations with low spectral resolution, and the real contribution of each component (PDR, jet) in complex environment of star-forming regions (SFRs) is poorly understood. We investigate the contribution of jet and PDR emission, and of absorption to the [OI]63micron line towards the ultra-compact H{\sc ii} region G5.89--0.39 and study its far-IR line luminosity in different velocity regimes through [OI], [CII], CO, OH, and H2O. We mapped G5.89--0.39 in [OI] and in CO(16--15) with the GREAT receiver onboard SOFIA. We observed the central position of the source in the OH^2\Pi_{3/2}, J=5/2\toJ=3/2 and ^2\Pi_{1/2}, J=3/2\toJ=1/2 lines. These data were complemented with APEX CO(6-5) and CO(7-6) and HIFI maps and single-pointing observations in [CII], H2O, and HF. The [OI] spectra in G5.89--0.39 are severely contaminated by absorptions from the envelope and from different clouds along the line of sight. Emission is detected only at HV, clearly associated with the compact north-south outflows traced by extremely HV low-J CO. The mass-loss rate and energetics of derived from [OI] agree well with estimates from CO, suggesting that the molecular outflows in G5.89--0.39 are driven by the jet system seen in [OI]. The far-IR line luminosity of G5.89--0.39 is dominated by [OI] at HV; the second coolant in this velocity regime is CO, while [CII], OH and H2O are minor contributors to the cooling in the outflow. Our study shows the importance of spectroscopically resolved data of [OI]63micron for using this line as diagnostic of SFRs. While this was not possible until now, the GREAT receiver onboard SOFIA has recently opened the possibility of detailed studies of this line to investigate its potential for probing different environments.

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.

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.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

Orion KL has served as a benchmark for spectral line searches throughout the (sub)millimeter regime. The main goal is to systematically study spectral characteristics of Orion KL in the 1.3 cm band. We carried out a spectral line survey (17.9 GHz to 26.2 GHz) with the Effelsberg-100 m telescope towards Orion KL. We find 261 spectral lines, yielding an average line density of about 32 spectral features per GHz above 3$\sigma$. The identified lines include 164 radio recombination lines (RRLs) and 97 molecular lines. A total of 23 molecular transitions from species known to exist in Orion KL are detected for the first time in the interstellar medium. Non-metastable 15NH3 transitions are detected in Orion KL for the first time. Based on the velocity information of detected lines and the ALMA images, the spatial origins of molecular emission are constrained and discussed. A narrow feature is found in SO2 ($8_{1,7}-7_{2,6}$), possibly suggesting the presence of a maser line. Column densities and fractional abundances relative to H2 are estimated for 12 molecules with LTE methods. Rotational diagrams of non-metastable 14NH3 transitions with J=K+1 to J=K+4 yield different results; metastable 15NH3 is found to have a higher excitation temperature than non-metastable 15NH3, indicating that they may trace different regions. Elemental and isotopic abundance ratios are estimated: 12C/13C=63+-17, 14N/15N=100+-51, D/H=0.0083+-0.0045. The dispersion of the He/H ratios derived from H$\alpha$/He$\alpha$ pairs to H$\delta$/He$\delta$ pairs is very small, which is consistent with theoretical predictions that the departure coefficients bn factors for hydrogen and helium are nearly identical. Based on a non-LTE code neglecting excitation by the infrared radiation field and a likelihood analysis, we find that the denser regions have lower kinetic temperature, which favors an external heating of the Hot Core.

OH231.8+4.2 is a well studied O-rich CSE around an intermediate-mass evolved star that displays bipolar molecular outflows accelerated up to 400 km/s. OH231 also presents an exceptional molecular richness probably due to shock-induced chemical processes. We report the first detection in this source of HNCO, HNCS, HC3N and NO, with the IRAM-30m telescope in a sensitive mm-wavelength survey towards this target. HNCO and HNCS are also first detections in CSEs. The observed line profiles show that the emission arises in the massive central component of the envelope and at the base of the fast lobes. The NO profiles are broader than those of HNCO, HNCS, and HC3N, and most importantly, broader than the line profiles of 13CO. This indicates that the NO abundance is enhanced in the fast lobes relative to the slow, central parts. From LTE and non-LTE excitation analysis, we estimate beam-average rotational temperatures of 15-30 K (and, maybe, up to 55 K for HC3N) and fractional abundances of X(HNCO)=[0.8-1]E-7, X(HNCS)=[0.9-1]E-8, X(HC3N)=[5-7]E-9 and X(NO)=[1-2]E-6. NO is, therefore, amongst the most abundant N-bearing species in OH231. We have performed thermodynamical chemical equilibrium and chemical kinetics models to investigate the formation of these N-bearing species in OH231. The model underestimates the observed abundances for HNCO, HNCS, and HC3N by several orders of magnitude, which indicates that these molecules can hardly be products of standard UV-photon and/or cosmic-ray induced chemistry in OH231, and that other processes (e.g. shocks) play a major role in their formation. For NO, the model abundance is compatible with the observed average value. The new detections presented in this work corroborate the particularly rich chemistry of OH231, which is likely profoundly influenced by shock-induced processes, as proposed in earlier works.

We present the G-virial method (available at http://gxli.github.io/G-virial/) which aims to quantify (1) the importance of gravity in molecular clouds in the position-position-velocity (PPV) space, and (2) properties of the gas condensations in molecular clouds. Different from previous approaches that calculate the virial parameter for different regions, our new method takes gravitational interactions between all the voxels in 3D PPV data cubes into account, and generates maps of the importance of gravity. This map can be combined with the original data cube to derive relations such as the mass-radius relation. Our method is important for several reasons. First, it offers the the ability to quantify the centrally condensed structures in the 3D PPV data cubes, and enables us to compare them in an uniform framework. Second, it allows us to understand the importance of gravity at different locations in the data cube, and provides a global picture of gravity in clouds. Third, it offers a robust approach to decomposing the data into different regions which are gravitationally coherent. To demonstrate the application of our method we identified regions from the Perseus and Ophiuchus molecular clouds, and analyzed their properties. We found an increase in the importance of gravity towards the centers of the individual molecular condensations. We also quantified the properties of the regions in terms of mass-radius and mass-velocity relations. Through evaluating the virial parameters based on the G-virial, we found that all our regions are almost gravitationally bound. Cluster-forming regions appear are more centrally condensed.

OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar envelopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO+ and H13CO+ molecular ions and the first detection of SO+, N2H+, and (tentatively) H3O+ in this source. SO+ and H3O+ have not been detected before in CSEs around evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 30 m radio telescope and with Herschel/HIFI. Except for H3O+, all the molecular ions detected in this work display emission lines with broad profiles (FWHM 50-90 km/s), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The narrow profile (FWHM 14 km/s) and high critical densities (>1e6cm-3 ) of the H3O+ transitions observed are consistent with this ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. From rotational diagram analysis, we deduce excitation temperatures of Tex 10-20 K for all ions except for H3O+, which is most consistent with Tex 100 K. Although uncertain, the higher excitation temperature suspected for H3O+ is similar to that recently found for H2O and a few other molecules, which selectively trace a previously unidentified, warm nebular component.The column densities of the molecular ions reported here are in the range Ntot [1-8]x1e13 cm-2, leading to beam-averaged fractional abundances relative to H2 of X(HCO+) 1e-8, X(H13CO+) 2e-9, X(SO+) 4e-9, X(N2H+) 2e-9, and X(H3O+) 7e-9 cm-2. We have performed chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated by cosmic-ray and UV-photon ionization. (abridged).

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 present a survey toward 16 Low Luminosity Objects (LLOs with an internal luminosity, Lint, lower than 0.2 Lsun) with N2H+ (1-0), N2H+ (3-2), N2D+ (3-2), HCO+ (3-2) and HCN (3-2) using the Arizona Radio Observatory Kitt Peak 12m Telescope and Submillimeter Telescope. Our goal is to probe the nature of these faint protostars which are believed to be either very low mass or extremely young protostars. We find that the N2D+/N2H+ column density ratios of LLOs are similar to those of typical starless cores and Class 0 objects. The N2D+/N2H+ column density ratios are relatively high (> 0.05) for LLOs with kinetic temperatures less than 10 K in our sample. The distribution of N2H+ (1-0) line widths spreads between that of starless cores and young Class 0 objects. If we use the line width as a dynamic evolutionary indicator, LLOs are likely young Class 0 protostellar sources. We further use the optically thick tracers, HCO+ (3-2) and HCN (3-2), to probe the infall signatures of our targets. We derive the asymmetry parameters from both lines and estimate the infall velocities by fitting the HCO+ (3-2) spectra with two-layer models. As a result, we identify eight infall candidates based on the infall velocities and seven candidates have infall signatures supported by asymmetry parameters from at least one of HCO+ (3-2) and HCN (3-2).

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

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

Water is a sensitive tracer of physical conditions in star-forming regions because of its large abundance variations between hot and cold regions. We use spectrally resolved observations of rotational lines of H$_2$O and its isotopologs with Herschel/HIFI to constrain the physical conditions of the water emitting region toward the high-mass protostar AFGL2591. We use analytical estimates and rotation diagrams to estimate T$_{ex}$ and column densities of H$_2$O of the envelope, the outflow, and a foreground cloud. Furthermore, we use the non-LTE radiative transfer code to estimate the temperature and volume density of the H$_2$O emitting gas. Assuming LTE, we estimate an T$_{ex}$ of 42 K and a column density of 2$\times$10$^{14}$ cm$^{-2}$ for the envelope and 45 K and 4$\times$10$^{13}$ cm$^{-2}$ for the outflow, in beams of 4" and 30", respectively. Non-LTE models indicate a kinetic temperature of 60-230 K and a volume density of 7$\times$10$^6$-10$^8$ cm$^{-3}$ for the envelope, and a kinetic temperature of 70-90 K and a gas density of 10$^7$-10$^8$ cm$^{-3}$ for the outflow. The o/p ratio of the foreground absorption is 1.9$\pm$0.4, suggesting a low temperature. In contrast, the o/p ratio seen in absorption by the outflow is 3.5$\pm$1.0, as expected for warm gas. The water abundance in the envelope is 10$^{-9}$, similar to the low values found for other high- and low-mass protostars, suggesting that this abundance is constant during the embedded phase of high-mass star formation. The water abundance in the outflow is 10$^{-10}$, which is 10$\times$ lower than in the envelope and in the outflows of high- and low-mass protostars. Since beam size effects can only increase this estimate by a factor of 2, we suggest that the water in the outflow is affected by dissociating UV radiation due to the low extinction in the outflow lobe.

IRC +10216 is the prototypical carbon star exhibiting an extended molecular circumstellar envelope. Its spectral properties are therefore the template for an entire class of objects. The main goal is to systematically study the $\lambda$ $\sim$1.3 cm spectral line characteristics of IRC +10216. We carried out a spectral line survey with the Effelsberg-100 m telescope toward IRC +10216. It covers the frequency range between 17.8 GHz and 26.3 GHz (K-band). In the circumstellar shell of IRC +10216, we find 78 spectral lines, among which 12 remain unidentified. The identified lines are assigned to 18 different molecules and radicals. A total of 23 lines from species known to exist in this envelope are detected for the first time outside the Solar System and there are additional 20 lines first detected in IRC +10216. The potential orgin of "U" lines is also discussed. Assuming local thermodynamic equilibrium (LTE), we then determine rotational temperatures and column densities of 17 detected molecules. Molecular abundances relative to H$_{2}$ are also estimated. A non-LTE analysis of NH$_{3}$ shows that the bulk of its emission arises from the inner envelope with a kinetic temperature of 70$\pm$20 K. Evidence for NH$_{3}$ emitting gas with higher kinetic temperature is also obtained, and potential abundance differences between various $^{13}$C-bearing isotopologues of HC$_{5}$N are evaluated. Overall, the isotopic $^{12}$C/$^{13}$C ratio is estimated to be 49$\pm$9. Finally, a comparison of detected molecules in the $\lambda$ $\sim$1.3 cm range with the dark cloud TMC-1 indicates that silicate-bearing molecules are more predominant in IRC +10216.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

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.

We report 1.3 cm and 6 cm continuum observations toward the massive proto-stellar candidate G11.11-0.12P1 using the Karl G. Jansky Very Large Array (VLA). We detect a string of four unresolved radio continuum sources coincident with the mid-IR source in G11P1. The continuum sources have positive spectral indices consistent with a thermal (free-free) ionized jet. The most likely origin of the ionized gas are shocks due to the interaction of a stellar wind with the surrounding high-density material. We also present NIR United Kingdom Infrared Telescope (UKIRT) archival data which show an extended structure detected only at K-band (2.2 $\mu$m), which is oriented perpendicular to the jet, and that may be scattered light from a circumstellar disk around the massive protostar. Our observations plus the UKIRT archival data thus provide new evidence that a disk/jet system is present in the massive protostellar candidate located in the G11.11-0.12P1 core.

The formation of hypercompact (HC) HII regions is an important stage in massive star formation. Spectral line and continuum observations can explore its dynamic conditions. We present high angular resolution observations carried out with the Submillimeter Array (SMA) and the Very Large Array (VLA) toward the HC HII region G35.58-0.03. With the 1.3 mm SMA and 1.3 cm VLA, we detected a total of about 25 transitions of 8 different species and their isotopologues (CO, CH$_3$CN, SO$_2$, CH$_3$CCH, OCS, CS, H, and NH$_{3}$). G35.58-0.03 consists of an HC HII core with electron temperature $T_e^*\geqq5500$ K, emission measure EM $\approx 1.9\times10^{9}$ pc cm$^{-6}$, local volume electron density $n_e=3.3\times10^{5}$ cm$^{-3}$, and a same width of radio recombination line FWHM $\approx$ 43.2 km s$^{-1}$ for both H30$\alpha$ and H38$\beta$ at its intrinsic core size $\sim$3714 AU. The H30$\alpha$ line shows evidence of an ionized outflow driving a molecular outflow. Based on the derived Lyman continuum flux, there should be an early-type star equivalent to O6.5 located inside the HII region. From the continuum spectral energy distribution from 3.6 cm, 2.0 cm, 1.3 cm, 1.3 mm and 0.85 mm to 0.45 mm, we distinguished the free-free emission (25% $\sim$ 55%) from the warm dust component (75% $\sim$ 45%) at 1.3 mm. The molecular envelope shows evidence of infall and outflow with an infall rate 0.05 M$_{ \odot}$ yr$^{-1}$ and a mass loss rate $5.2\times10^{-3}$ ${\rm M_{ \odot} yr^{-1}}$. The derived momentum ($\sim0.05$ M$_{\odot}$ km s$^{-1}$) is consistent between the infalling and outflowing gas per year. It is suggested that the infall is predominant and the envelope mass of dense core is increasing rapidly, but the accretion in the inner part might already be halted.

Context. The earliest phases of high-mass star formation are poorly understood. Aims. Our goal is to determine the physical conditions and kinematic structure of massive star-forming cloud clumps. Methods. We analyze H$_2$O 557 GHz line profiles observed with HIFI toward four positions in two infrared-dark cloud clumps . By comparison with ground-based C$^{17}$O, N$_2$H$^+$, CH$_3$OH and NH$_3$ line observations, we constrain the volume density and kinetic temperature of the gas and estimate the column density and abundance of H$_2$O and N$_2$H$^+$. Results. The observed water lines are complex with emission and absorption components. The absorption is red shifted and consistent with a cold envelope, while the emission is interpreted as resulting from protostellar outflows. The gas density in the clumps is $\sim$ 10$^7$ cm$^{-3}$. The o-H2O outflow column density is 0.3 to 3.0 10$^{14}$ cm$^{-2}$, the o-H2O absorption column density is between 1.5 10$^{14}$ and 2.6 10$^{15}$ cm$^-2$ with cold o-H2O abundances between 1.5 10$^{-9}$ and 3.1 10$^{-8}$. Conclusions. All clumps have high gas densities ($\sim$ 10$^7$ cm$^{-3}$) and display infalling gas. Three of the four clumps have outflows. The clumps form an evolutionary sequence as probed by H$_2$O N$_2$H$^+$, NH$_3$ and CH$_3$OH. We find that G28-MM is the most evolved, followed by G11-MM then G28-NH3. The least evolved clump is G11-NH3 which shows no sign-posts of star-formation. G11-NH3 is a high-mass pre-stellar core.

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.

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.

We present the results of Herschel HOBYS photometric mapping combined with BIMA observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra compact (UC) HII region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow, and finally the NH_2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UCHII region and discuss the Aquila supershell expansion as a mayor influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large scale structures through kinematic information.

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.

We present Submillimeter Array (SMA) $\lambda =$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

We present an analysis of the dust continuum emission at 870 um in order to investigate the mass distribution of clumps within infrared dark clouds (IRDCs). We map six IRDCs with the Large APEX BOlometer CAmera (LABOCA) at APEX, reaching an rms noise level of 28-44 mJy/beam. The dust continuum emission coming from these IRDCs was decomposed by using two automated algorithms, Gaussclumps and Clumpfind. Moreover, we carried out single-pointing observations of the N_2H^+ (3-2) line toward selected positions to obtain kinematic information. The mapped IRDCs are located in the range of kinematic distances of 2.7-3.2 kpc. We identify 510 and 352 sources with Gaussclumps and Clumpfind, respectively, and estimate masses and other physical properties assuming a uniform dust temperature. The mass ranges are 6-2692 Msun (Gaussclumps) and 7-4254 Msun (Clumpfind) and the ranges in effective radius are around 0.10-0.74 pc (Gaussclumps) and 0.16-0.99 pc (Clumpfind). The mass distribution, independent of the decomposition method used, is fitted by a power law, dN/dM propto M^alpha, with an index of -1.60 +/- 0.06, consistent with the CO mass distribution and other high-mass star-forming regions.

(Abridged) We study the response of the gas to energetic processes associated with high-mass star formation and compare it with studies on low- and intermediate-mass young stellar objects (YSOs) using the same methods. The far-IR line emission and absorption of CO, H$_2$O, OH, and [OI] reveals the excitation and the relative contribution of different species to the gas cooling budget. Herschel-PACS spectra covering 55-190 um are analyzed for ten high-mass star forming regions of various luminosities and evolutionary stages at spatial scales of ~10^4 AU. Radiative transfer models are used to determine the contribution of the envelope to the far-IR CO emission. The close environments of high-mass YSOs show strong far-IR emission from molecules, atoms, and ions. Water is detected in all 10 objects even up to high excitation lines. CO lines from J=14-13 up to typically 29-28 show a single temperature component, Trot~300 K. Typical H$_2$O temperatures are Trot~250 K, while OH has Trot~80 K. Far-IR line cooling is dominated by CO (~75 %) and to a smaller extent by OI (~20 %), which increases for the most evolved sources. H$_2$O is less important as a coolant for high-mass sources because many lines are in absorption. Emission from the envelope is responsible for ~45-85 % of the total CO luminosity in high-mass sources compared with only ~10 % for low-mass YSOs. The highest-J lines originate most likely from shocks, based on the strong correlation of CO and H$_2$O with physical parameters of the sources from low- to high-masses. Excitation of warm CO is very similar for all mass regimes, whereas H$_2$O temperatures are ~100 K higher for high-mass sources than the low-mass YSOs. Molecular cooling is ~4 times more important than cooling by [OI]. The total far-IR line luminosity is about 10$^{-3}$ and 10$^{-5}$ times lower than the dust luminosity for the low- and high-mass YSOs.

Star formation occurs in molecular gas. In previous studies, the structure of the molecular gas has been studied in terms of molecular clouds, but has been overlooked beyond the cloud scale. We present an observational study of the molecular gas at 49.5 degree <l<52.5 degree and -5.0 km/s <v_lsr <17.4 km/s. The molecular gas is found in the form of a huge (>= 500 pc) filamentary gas wisp. This has a large physical extent and a velocity dispersion of ~5 km/s. The eastern part of the filamentary gas wisp is located ~130 pc above the Galactic disk (which corresponds to 1.5-4 e-folding scale-heights), and the total mass of the gas wisp is >= 1 X 10^5 M_sun. It is composed of two molecular clouds and an expanding bubble. The velocity structure of the gas wisp can be explained as a smooth quiescent component disturbed by the expansion of a bubble. That the length of the gas wisp exceeds by much the thickness of the molecular disk of the Milky Way is consistent with the cloud-formation scenario in which the gas is cold prior to the formation of molecular clouds. Star formation in the filamentary gas wisp occurs at the edge of a bubble (G52L nebula), which is consistent with some models of triggered star formation.

Stars are born within dense clumps of giant molecular clouds, constituting young stellar agglomerates known as embedded clusters, which only evolve into bound open clusters under special conditions. We statistically study all embedded clusters (ECs) and open clusters (OCs) known so far in the inner Galaxy, investigating particularly their interaction with the surrounding molecular environment and the differences in their evolution. We first compiled a merged list of 3904 clusters from optical and infrared clusters catalogs in the literature, including 75 new (mostly embedded) clusters discovered by us in the GLIMPSE survey. From this list, 695 clusters are within the Galactic range |l| < 60 deg and |b| < 1.5 deg covered by the ATLASGAL survey, which was used to search for correlations with submm dust continuum emission tracing dense molecular gas. We defined an evolutionary sequence of five morphological types: deeply embedded cluster (EC1), partially embedded cluster (EC2), emerging open cluster (OC0), OC still associated with a submm clump in the vicinity (OC1), and OC without correlation with ATLASGAL emission (OC2). Together with this process, we performed a thorough literature survey of these 695 clusters, compiling a considerable number of physical and observational properties in a catalog that is publicly available. We found that an OC defined observationally as OC0, OC1, or OC2 and confirmed as a real cluster is equivalent to the physical concept of OC (a bound exposed cluster) for ages in excess of ~16 Myr. Some observed OCs younger than this limit can actually be unbound associations. We found that our OC and EC samples are roughly complete up to ~1 kpc and ~1.8 kpc from the Sun, respectively, beyond which the completeness decays exponentially. Using available age estimates for a few ECs, we derived an upper limit of 3 Myr for the duration of the embedded phase... (Abridged)

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.

We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas.

Protostellar outflow is a prominent process that accompanies the formation of stars. It is generally agreed that wide-angled protostellar outflows come from the interaction between the wind from a forming star and the ambient gas. However, it is still unclear how the interaction takes place. In this work, we theoretically investigate the possibility that the outflow results from interaction between the wind and the ambient gas in the form of turbulent entrainment. In contrast to the previous models, turbulent motion of the ambient gas around the protostar is taken into account. In our model, the ram-pressure of the wind balances the turbulent ram-pressure of the ambient gas, and the outflow consists of the ambient gas entrained by the wind. The calculated outflow from our modelling exhibits a conical shape. The total mass of the outflow is determined by the turbulent velocity of the envelope as well as the outflow age, and the velocity of the outflow is several times higher than the velocity dispersion of the ambient gas. The outflow opening angle increases with the strength of the wind and decreases with the increasing ambient gas turbulence. The outflow exhibits a broad line width at every position. We propose that the turbulent entrainment process, which happens ubiquitously in nature, plays a universal role in shaping protostellar outflows.

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

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.

In order to constrain the broad-band spectral energy distribution of the afterglow of GRB 100621A, dedicated observations were performed in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at the 2.2m MPG/ESO telescope, in the sub-millimeter band with the large bolometer array LABOCA at APEX, and at radio frequencies with ATCA. Utilizing also Swift X-ray observations, we attempt an interpretation of the observational data within the fireball scenario. The afterglow of GRB 100621A shows a very complex temporal as well as spectral evolution. We identify three different emission components, the most spectacular one causing a sudden intensity jump about one hour after the prompt emission. The spectrum of this component is much steeper than the canonical afterglow. We interpret this component using the prescription of Vlasis et al. (2011) for a two-shell collision after the first shell has been decelerated by the circumburst medium. We use the fireball scenario to derive constraints on the microphysical parameters of the first shell. Long-term energy injection into a narrow jet seems to provide an adequate description. Another noteworthy result is the large ($A_V$ = 3.6 mag) line-of-sight host extinction of the afterglow in an otherwise extremely blue host galaxy.

To understand the origin of water line emission and absorption during high-mass star formation, we decompose high-resolution Herschel-HIFI line spectra toward 19 high-mass star-forming regions into three distinct physical components. Protostellar envelopes are usually seen as narrow absorptions or emissions in the H2O 1113 and 1669 GHz ground-state lines, the H2O 987 GHz excited-state line, and the H2O-18 1102 GHz ground-state line. Broader features due to outflows are usually seen in absorption in the H2O 1113 and 1669 GHz lines, in 987 GHz emission, and not seen in H2O-18, indicating a low column density and a high excitation temperature. The H2O 1113 and 1669 GHz spectra show narrow absorptions by foreground clouds along the line of sight, which have a low column density and a low excitation temperature, although their H2O ortho/para ratios are close to 3. The intensities of the H2O 1113 and 1669 GHz lines do not show significant trends with luminosity, mass, or age. In contrast, the 987 GHz line flux increases with luminosity and the H2O-18 line flux decreases with mass. Furthermore, appearance of the envelope in absorption in the 987 GHz and H2O-18 lines seems to be a sign of an early evolutionary stage. We conclude that the ground state transitions of H2O trace the outer parts of the envelopes, so that the effects of star formation are mostly noticeable in the outflow wings. These lines are heavily affected by absorption, so that line ratios of H2O involving the ground states must be treated with caution. The average H2O abundance in high-mass protostellar envelopes does not change much with time. The 987 GHz line appears to be a good tracer of the mean weighted dust temperature of the source, which may explain why it is readily seen in distant galaxies.

The goal of this paper is to determine the respective importance of solid state vs. gas phase reactions for the formation of dimethyl ether. This is done by a detailed analysis of the excitation properties of the ground state and the torsionally excited states, v11=1 and v15=1, toward the high-mass star-forming region G327.3-0.6. With the Atacama Pathfinder EXperiment 12 m submillimeter telescope, we performed a spectral line survey. The observed spectrum is modeled assuming local thermal equilibrium. CH3OCH3 has been detected in the ground state, and in the torsionally excited states v11=1 and v15=1, for which lines have been detected here for the first time. The emission is modeled with an isothermal source structure as well as with a non-uniform spherical structure. For non-uniform source models one abundance jump for dimethyl ether is sufficient to fit the emission, but two components are needed for the isothermal models. This suggests that dimethyl ether is present in an extended region of the envelope and traces a non-uniform density and temperature structure. Both types of models furthermore suggest that most dimethyl ether is present in gas that is warmer than 100 K, but a smaller fraction of 5%-28% is present at temperatures between 70 and 100 K. The dimethyl ether present in this cooler gas is likely formed in the solid state, while gas phase formation probably is dominant above 100 K. Finally, the v11=1 and v15=1 torsionally excited states are easily excited under the density and temperature conditions in G327.3-0.6 and will thus very likely be detectable in other hot cores as well.

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.

(Abridged) Far-infrared Herschel-PACS spectra of 18 low-mass protostars of various luminosities and evolutionary stages are studied. We quantify their far-infrared line emission and the contribution of different atomic and molecular species to the gas cooling budget during protostellar evolution. We also determine the spatial extent of the emission and investigate the underlying excitation conditions. Most of the protostars in our sample show strong atomic and molecular far-infrared emission. Water is detected in 17 objects, including 5 Class I sources. The high-excitation H2O line at 63.3 micron is detected in 7 sources. CO transitions from J=14-13 up to 49-48 are found and show two distinct temperature components on Boltzmann diagrams with rotational temperatures of ~350 K and ~700 K. H2O has typical excitation temperatures of ~150 K. Emission from both Class 0 and I sources is usually spatially extended along the outflow direction but with a pattern depending on the species and the transition. The H2O line fluxes correlate strongly with those of the high-J CO lines, as well as with the bolometric luminosity and envelope mass. They correlate less strongly with OH and not with [OI] fluxes. The PACS data probe at least two physical components. The H2O and CO emission likely arises in non-dissociative (irradiated) shocks along the outflow walls with a range of pre-shock densities. Some OH is also associated with this component, likely resulting from H2O photodissociation. UV-heated gas contributes only a minor fraction to the CO emission observed by PACS, based on the strong correlation between the shock-dominated CO 24-23 line and the CO 14-13 line. [OI] and some of the OH emission probe dissociative shocks in the inner envelope. The total far-infrared cooling is dominated by H2O and CO, with [OI] increasing for Class I sources.

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

We present the results of combined NH3(1,1) and (2,2) line emission observed with the Very Large Array and the Effelsberg 100m telescope of the Infrared Dark Cloud G14.225-0.506. The NH3 emission reveals a network of filaments constituting two hub-filament systems. Hubs are associated with gas of rotational temperature Trot \sim 25 K, non-thermal velocity dispersion ~1.1 km/s, and exhibit signs of star formation, while filaments appear to be more quiescent (Trot \sim 11 K, non-thermal velocity dispersion ~0.6 km/s). Filaments are parallel in projection and distributed mainly along two directions, at PA \sim 10 deg and 60 deg, and appear to be coherent in velocity. The averaged projected separation between adjacent filaments is between 0.5 pc and 1pc, and the mean width of filaments is 0.12 pc. Cores within filaments are separated by ~0.33 pc, which is consistent with the predicted fragmentation of an isothermal gas cylinder due to the 'sausage'-type instability. The network of parallel filaments observed in G14.225-0.506 is consistent with the gravitational instability of a thin gas layer threaded by magnetic fields. Overall, our data suggest that magnetic fields might play an important role in the alignment of filaments, and polarization measurements in the entire cloud would lend further support to this scenario.

Most studies of high-mass star formation focus on massive luminous clumps, but the physical properties of their larger scale environment are poorly known. In this work, we aim at characterising the effects of clustered star formation and feedback of massive stars on the surrounding medium by studying the distribution of warm gas through mid-J 12CO and 13CO data. We present APEX 12CO(6-5), (7-6), 13CO(6-5), (8-7) and HIFI 13CO(10-9) maps of the star forming region G327.36-0.6. We infer the physical properties of the emitting gas on large scales through a LTE analysis, while we apply a more sophisticated LVG approach on selected positions. Maps of all lines are dominated in intensity by the PDR around the Hii region G327.3-0.5. Mid-J 12CO emission is detected over the whole extent of the maps with excitation temperatures ranging from 20K up to 80K in the gas around the Hii region, and H2 column densities from few 10^21 cm-2 in the inter-clump gas to 3 10^22 cm-2 towards the hot core G327.3-0.6. The warm gas (traced by 12 and 13CO(6-5) emission) is only a small percentage (10%) of the total gas in the IRDC, while it reaches values up to 35% of the total gas in the ring surrounding the Hii region. The 12CO ladders are qualitatively compatible with PDR models for high density gas, but the much weaker than predicted 13CO emission suggests that it comes from a large number of clumps along the line of sight. All lines are detected in the inter-clump gas when averaged over a large region with an equivalent radius of 50" (~0.8pc), implying that the mid-J 12CO and 13CO inter-clump emission is due to high density components with low filling factor. Finally, the detection of the 13CO(10-9) line allows to disentangle the effects of gas temperature and gas density on the CO emission, which are degenerate in the APEX observations alone.

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.

The onset of massive star formation is not well understood because of observational and theoretical difficulties. To find the dense and cold clumps where massive star formation can take place, we compiled a sample of high infrared extinction clouds, which were observed previously by us in the 1.2 mm continuum emission and ammonia. We try to understand the star-formation stages of the clumps in these high extinction clouds by studying the infall and outflow properties, the presence of a young stellar object (YSO), and the level of the CO depletion through a molecular line survey with the IRAM 30m and APEX 12m telescopes. Moreover, we want to know if the cloud morphology, quantified through the column density contrast between the clump and the clouds, has an impact on the star formation occurring inside it. We find that the HCO+(1-0) line is the most sensitive for detecting infalling motions. SiO, an outflow tracer, was mostly detected toward sources with infall, indicating that infall is accompanied by collimated outflows. The presence of YSOs within a clump depends mostly on its column density; no signs of YSOs were found below 4E22 cm-2. Star formation is on the verge of beginning in clouds that have a low column density contrast; infall is not yet present in the majority of the clumps. The first signs of ongoing star formation are broadly observed in clouds where the column density contrast between the clump and the cloud is higher than two; most clumps show infall and outflow. Finally, the most evolved clumps are in clouds that have a column density contrast higher than three; almost all clumps have a YSO, and in many clumps, the infall has already halted. Hence, the cloud morphology, based on the column density contrast between the cloud and the clumps, seems to have a direct connection with the evolutionary stage of the objects forming inside.

High spatial resolution low-J 12CO observations have shown that the wide-angle outflow seen in the Orion BN/KL region correlates with the famous H2 fingers. Recently, high-resolution large-scale mappings of mid- and higher-J CO emissions have been reported toward the Orion molecular cloud 1 core region using the APEX telescope. Therefore, it is of interest to investigate this outflow in the higher-J 12CO emission, which is likely excited by shocks. The observations were carried out using the dual-color heterodyne array CHAMP+ on the APEX telescope. The images of the Orion BN/KL region were obtained in the 12CO J=6-5 and J=7-6 transitions with angular resolutions of 8.6 and 7.4 arcsec, respectively. The results show a good agreement between our higher-J 12CO emission and SMA low-J 12CO data, which indicates that this wide-angle outflow in Orion BN/KL is likely the result of an explosive event that is related to the runaway objects from a dynamically decayed multiple system. From our observations, we estimate that the kinetic energy of this explosive outflow is about 1-2x10^47 erg. In addition, a scenario has been proposed where part of the outflow is decelerated and absorbed in the cloud to explain the lack of CO bullets in the southern part of BN/KL, which in turn induces the methanol masers seen in this region.

The APEX Telescope Large Area Survey: The Galaxy (ATLASGAL) is an unbiased continuum survey of the inner Galactic disk at 870 \mu m. It covers +/- 60 deg in Galactic longitude and aims to find all massive clumps at various stages of high-mass star formation in the inner Galaxy, particularly the earliest evolutionary phases. We aim to determine properties such as the gas kinetic temperature and dynamics of new massive cold clumps found by ATLASGAL. Most importantly, we derived their kinematical distances from the measured line velocities. We observed the ammonia (J,K) = (1,1) to (3,3) inversion transitions toward 862 clumps of a flux-limited sample of submm clumps detected by ATLASGAL and extracted 13CO (1-0) spectra from the Galactic Ring Survey (GRS). We determined distances for a subsample located at the tangential points (71 sources) and for 277 clumps whose near/far distance ambiguity is resolved. Most ATLASGAL clumps are cold with rotational temperatures from 10-30 K. They have a wide range of NH3 linewidths, which by far exceeds the thermal linewidth, as well as a broad distribution of high column densities with an NH3 abundance in the range of 5 to 30 * 10^{-8}. We found an enhancement of clumps at Galactocentric radii of 4.5 and 6 kpc. The high detection rate (87%) confirms ammonia as an excellent probe of the molecular content of the massive, cold clumps revealed by ATLASGAL. A clear trend of increasing rotational temperatures and linewidths with evolutionary stage is seen for source samples ranging from 24 \mu m dark clumps to clumps with embedded HII regions. The survey provides the largest ammonia sample of high-mass star forming clumps and thus presents an important repository for the characterization of statistical properties of the clumps and the selection of subsamples for detailed, high-resolution follow-up studies.

Methanol masers at 6.7 GHz are known to be tracers of high-mass star formation in our Galaxy. In this paper, we study the large scale physical conditions in the star forming clumps/cores associated with 6.7 GHz methanol masers using observations of the (1,1), (2,2) and (3,3) inversion transitions of ammonia with the Effelsberg telescope. The gas kinetic temperature is found to be higher than in infrared dark clouds, highlighting the relatively evolved nature of the maser sources. Other than a weak correlation between maser luminosity and the ammonia line width, we do not find any differences between low and high luminosity methanol masers.

We present Herschel/HIFI observations of fourteen water lines in W43-MM1, a massive protostellar object in the luminous star cluster-forming region W43. We analyze the gas dynamics from the line profiles using Herschel-HIFI observations (WISH-KP) of fourteen far-IR water lines (H2O, H217O, H218O), CS(11-10), and C18O(9-8) lines, and using our modeling of the continuum spectral energy distribution. As for lower mass protostellar objects, the molecular line profiles are a mix of emission and absorption, and can be decomposed into 'medium', and 'broad' velocity components. The broad component is the outflow associated with protostars of all masses. Our modeling shows that the remainder of the water profiles can be well fitted by an infalling and passively heated envelope, with highly supersonic turbulence varying from 2.2 km/s in the inner region to 3.5 km/s in the outer envelope. Also, W43-MM1 has a high accretion rate, between 4.0 x 10^{-4} and 4.0 x 10^{-2} \msun /yr, derived from the fast (0.4-2.9 km/s) infall observed. We estimate a lower mass limit of gaseous water of 0.11 \msun and total water luminosity of 1.5 \lsun (in the 14 lines presented here). The central hot core is detected with a water abundance of 1.4 x 10^{-4} while the water abundance for the outer envelope is 8 x10^{-8}. The latter value is higher than in other sources, most likely related to the high turbulence and the micro-shocks created by its dissipation. Examining water lines of various energies, we find that the turbulent velocity increases with the distance to the center. While not in clear disagreement with the competitive accretion scenario, this behavior is predicted by the turbulent core model. Moreover, the estimated accretion rate is high enough to overcome the expected radiation pressure.

(Abridged) Water is a key tracer of dynamics and chemistry in low-mass protostars, but spectrally resolved observations have so far been limited in sensitivity and angular resolution. In this first systematic survey of spectrally resolved water emission in low-mass protostellar objects, H2O was observed in the ground-state transition at 557 GHz with HIFI on Herschel in 29 embedded Class 0 and I protostars. Complementary far-IR and sub-mm continuum data (including PACS data from our program) are used to constrain the spectral energy distribution of each source. H2O intensities are compared to inferred envelope and outflow properties and CO 3-2 emission. H2O emission is detected in all objects except one. The line profiles are complex and consist of several kinematic components. The profiles are typically dominated by a broad Gaussian emission feature, indicating that the bulk of the water emission arises in outflows, not the quiescent envelope. Several sources show multiple shock components in either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles in 7 sources (6 Class 0, 1 Class I) indicative of infalling envelopes, and regular P-Cygni profiles in 4 sources (3 Class I, 1 Class 0) indicative of expanding envelopes. "Bullets" moving at >50 km/s are seen in 4 Class 0 sources; 3 of these are new detections. In the outflow, the H2O/CO abundance ratio as a function of velocity is nearly the same for all sources, increasing from 10^-3 at <5 km/s to >10^-1 at >10 km/s. The H2O abundance in the outer envelope is low, ~10^-10. The different H2O profile components show a clear evolutionary trend: in the Class 0 sources, emission is dominated by outflow components originating inside an infalling envelope. When the infall diminishes during the Class I phase, the outflow weakens and H2O emission disappears.

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.

Cloud contraction and infall are the fundamental processes of star formation. While "blue-skewed" line profiles observed in high- mass star forming regions are commonly taken as evidence of infall by an ever increasing number of studies, their interpretation offers many pitfalls. Detecting infall via redshifted absorption in front of continuum sources is a much more direct and reliable method but so far mostly restricted toward absorption in the centimeter toward strong HII regions. Here we present a novel approach by probing absorption of rotational ammonia transitions in front of the strong dust emission of massive star-forming regions. A carefully selected sample of three regions with different stages of evolution is selected to study infall through the evolution of massive star-forming clumps. Redshifted absorption is detected toward all three sources and infall rates between 3-10x10-3 Msol yr-1 are derived.

We present and analyze two spectrally resolved high-J CO lines towards the molecular outflow Cep E, driven by an intermediate-mass class 0 protostar. Using the GREAT receiver on board SOFIA, we observed the CO (12--11) and (13--12) transitions (E_u ~ 430 and 500 K, respectively) towards one position in the blue lobe of this outflow, that had been known to display high-velocity molecular emission. We detect the outflow emission in both transitions, up to extremely high velocities (~ 100 km/s with respect to the systemic velocity). We divide the line profiles into three velocity ranges that each have interesting spectral features: standard, intermediate, and extremely high-velocity. One distinct bullet is detected in each of the last two. A large velocity gradient analysis for these three velocity ranges provides constraints on the kinetic temperature and volume density of the emitting gas, >~ 100 K and > ~ 10^4 cm^-3, respectively. These results are in agreement with previous ISO observations and are comparable with results obtained by Herschel for similar objects. We conclude that high-J CO lines are a good tracer of molecular bullets in protostellar outflows. Our analysis suggests that different physical conditions are at work in the intermediate velocity range compared with the standard and extremely high-velocity gas at the observed position.

In this study, we present an unbiased sample of the earliest stages of massive star formation across 20 square-degree of the sky. Within the region 10deg < l < 20deg and |b| < 1deg, we search the ATLASGAL survey at 870 micron for dense gas condensations. These clumps are carefully examined for indications of ongoing star formation using YSOs from the GLIMPSE source catalog as well as sources in the 24 micron MIPSGAL images, to search for starless clumps. We calculate the column densities as well as the kinematic distances and masses for sources where the v_lsr is known from spectroscopic observations. Within the given region, we identify 210 starless clumps with peak column densities > 1 x 10e23 cm^(-2). In particular, we identify potential starless clumps on the other side of the Galaxy. The sizes of the clumps range between 0.1 pc and 3 pc with masses between a few tens of solar masses up to several ten thousands of solar masses. Most of them may form massive stars, but in the 20 square-degree we only find 14 regions massive enough to form stars more massive than 20 solar masses and 3 regions with the potential to form stars more massive than 40 40 solar masses. The slope of the high-mass tail of the clump mass function for clumps on the near side of the Galaxy is 2.2 and, therefore, Salpeter-like. We estimate the lifetime of the most massive starless clumps to be 60000 yr. The sample offers a uniform selection of starless clumps. In the large area surveyed, we only find a few potential precursors of stars in the excess of 40 solar masses. It appears that the lifetime of these clumps is somewhat shorter than their free-fall times, although both values agree within the errors. In addition, these are ideal objects for detailed studies and follow-up observations.

The ATLASGAL 870mum continuum survey conducted with the APEX telescope is the first survey covering the whole inner Galactic plane (60deg>l>-60deg & b<1.5deg) in submm continuum emission tracing the cold dust of dense and young star-forming regions. Here, we present the overall distribution of sources within our Galactic disk. The submm continuum emission is confined to a narrow range around the galactic plane, but shifted on average by ~0.07deg below the plane. Source number counts show strong enhancements toward the Galactic center, the spiral arms and toward prominent star-forming regions. Comparing the distribution of ATLASGAL dust continuum emission to that of young intermediate- to high-mass young stellar objects (YSOs) derived from Spitzer data, we find similarities as well as differences. In particular, the distribution of submm dust continuum emission is significantly more confined to the plane than the YSO distribution (FWHM of 0.7 and 1.1deg, corresponding to mean physical scale heights of approximately 46 and 80pc, respectively). While this difference may partly be caused by the large extinction from the dense submm cores, gradual dispersal of stellar distributions after their birth could also contribute to this effect. Compared to other tracers of Galactic structure, the ATLASGAL data are strongly confined to a narrow latitude strip around the Galactic plane.

A high density portion of the Orion Molecular Cloud 1 (OMC-1) contains the prominent, warm Kleinmann-Low (KL) nebula that is internally powered by an energetic event plus a farther region in which intermediate to high mass stars are forming. Its outside is affected by ultraviolet radiation from the neighboring Orion Nebula Cluster and forms the archetypical photon-dominated region (PDR) with the prominent bar feature. Its nearness makes the OMC-1 core region a touchstone for research on the dense molecular interstellar medium and PDRs. Using the Atacama Pathfinder Experiment telescope (APEX), we have imaged the line emission from the multiple transitions of several carbon monoxide (CO) isotopologues over the OMC-1 core region. Our observations employed the 2x7 pixel submillimeter CHAMP+ array to produce maps (~ 300 arcsec x 350 arcsec) of 12CO, 13CO, and C18O from mid-J transitions (J=6-5 to 8-7). We also obtained the 13CO and C18O J=3-2 images toward this region. The 12CO line emission shows a well-defined structure which is shaped and excited by a variety of phenomena, including the energetic photons from hot, massive stars in the nearby Orion Nebula's central Trapezium cluster, active high- and intermediate-mass star formation, and a past energetic event that excites the KL nebula. Radiative transfer modeling of the various isotopologic CO lines implies typical H2 densities in the OMC-1 core region of ~10^4-10^6 cm^-3 and generally elevated temperatures (~ 50-250 K). We estimate a warm gas mass in the OMC-1 core region of 86-285 solar masses.

Atacama Pathfinder Experiment (APEX) mapping observations in CO (9-8) and (4-3) toward a high-mass star-forming region, NGC 6334 I, are presented. The CO (9-8) map has a 6.4 arcsec resolution, revealing a ~0.5 pc, jet-like, and bipolar outflow. This is the first map of a molecular outflow in a THz line. The CO (9-8) and (4-3) lines arising from the outflow lobes both show extremely high-velocity line wings, and their ratios indicate a gas temperature greater than 100 K and a density higher than 10^4 cm^{-3}. The spatial-velocity structure of the CO (9-8) data is typical of a bow-shock-driven flow, which is consistent with the association between the bipolar outflow and the infrared bow-shaped tips. In short, the observations unveil a highly-excited and collimated component in a bipolar outflow that is powered by a high-mass protostar, and provide insights into the driving mechanism of the outflow. Meanwhile, the observations demonstrate that high-quality mapping observations can be performed with the new THz receiver on APEX.

We describe a pilot survey conducted with the Mopra 22-m radio telescope in preparation for the Millimeter Astronomy Legacy Team Survey at 90 GHz (MALT90). We identified 182 candidate dense molecular clumps using six different selection criteria and mapped each source simultaneously in 16 different lines near 90 GHz. We present a summary of the data and describe how the results of the pilot survey shaped the design of the larger MALT90 survey. We motivate our selection of target sources for the main survey based on the pilot detection rates and demonstrate the value of mapping in multiple lines simultaneously at high spectral resolution.

Infrared dark clouds are massive, dense clouds seen in extinction against the IR Galactic background. Many of these objects appear to be on the verge of star and star cluster formation. Our aim is to understand the physical properties of IRDCs in very early evolutionary phases. We selected the filamentary IRDC G351.77 - 0.51, which is remarkably IR quiet at 8{\mu}m. As a first step, we observed mm dust continuum emission and rotational lines of moderate and dense gas tracers to characterise different condensations along the IRDC and study the velocity field of the filament. Our initial study confirms coherent velocity distribution along the infrared dark cloud ruling out any coincidental projection effects. Excellent correlation between MIR extinction, mm continuum emission and gas distribution is found. Large-scale turbulence and line profiles throughout the filament is indicative of a shock in this cloud. Excellent correlation between line width, and MIR brightness indicates turbulence driven by local star formation.

We report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.

We present a new high-resolution study of pre-protocluster regions in tracers exclusively probing the coldest and dense gas (NH_2D). The data are used to constrain the chemical, thermal, kinematic, and physical conditions (i.e., densities) in G29.96e and G35.20w. NH_3, NH_2D, and continuum emission were mapped using the VLA, and PdBI. In particular, NH_2D is a unique tracer of cold, precluster gas at high densities, while NH_3 traces both the cold and warm gas of modest-to-high densities. In G29.96e, Spitzer images reveal two massive filaments, one of them in extinction (infrared dark cloud). We observe very low line widths in NH_3 (FWHM <1km/s). These multi-wavelength, high-resolution observations of high-mass pre-protocluster regions show that the target regions are characterized by (i) turbulent Jeans fragmentation of massive clumps into cores (from a Jeans analysis); (ii) cores and clumps that are "over-bound/subvirial", i.e. turbulence is too weak to support them against collapse, meaning that (iii) some models of monolithic cloud collapse are quantitatively inconsistent with data; (iv) accretion from the core onto a massive star, which can (for observed core sizes and velocities) be sustained by accretion of envelope material onto the core, suggesting that (similar to competitive accretion scenarios) the mass reservoir for star formation is not necessarily limited to the natal core; (v) high deuteration ratios ([NH_2D/NH_3]>6%), which make the above discoveries possible; (vi) and the destruction of NH_2D toward embedded stars. [abridged]

Context: High-mass star formation has been a very active field over the last decade, however, most studies targeted regions of luminosities between 10^4 and 10^5 L_sun. Methods: We selected the W31 star-forming complex with a total luminosity of ~6x10^6 L_sun for a multi-wavelength spectral line and continuum study covering wavelengths from the near- and mid-infrared via (sub)mm wavelength observations to radio data in the cm regime. Results: While the overall structure of the multi-wavelength continuum data resembles each other well, there are several intriguing differences. The 24mum emission stemming largely from small dust grains follows tightly the spatial structure of the cm emission tracing the ionized free-free emission. Hence warm dust resides in regions that are spatially associated with the ionized hot gas (~10^4 K) of the HII regions. Furthermore, we find several evolutionary stages within the same complexes, ranging from infrared-observable clusters, via deeply embedded regions associated with active star formation traced by 24\,$\mu$m and cm emission, to at least one high-mass gas clump devoid of any such signature. The 13CO(2-1) and C18O(2-1) spectral line observations reveal a large kinematic breadth in the entire region with a total velocity range of approximately 90 km/s. While the average virial mass ratio for W31 is close to unity, the line width analysis indicates large-scale evolutionary differences between the southern and northern sub-regions (G10.2-0.3 and G10.3-0.1) of the whole W31 complex. The clump mass function - tracing cluster scales and not scales of individual stars - derived from the 875mum continuum data has a slope of 1.5+-0.3, consistent with previous cloud mass functions.

Massive stars form deeply embedded in dense molecular gas, which they stir and heat up and ionize. During an early phase, the ionization is confined to hypercompact HII regions, and the stellar radiation is entirely absorbed by dust, giving rise to a hot molecular core. To investigate the innermost structure of such high-mass star-forming regions, we observed vibrationally excited HCN (via the direct $\ell$-type transition of v2=1, $\Delta J$=0, J=13, which lies 1400 K above ground) toward the massive hot molecular cores G10.47+0.03, SgrB2-N, and SgrB2-M with the Very Large Array (VLA) at 7 mm, reaching a resolution of about 1000 AU (0.1"). We detect the line both in emission and in absorption against HII regions. The latter allows to derive lower limits on the column densities of hot HCN, which are several times $10^{19}$ cm$^{-2}$. We see indication of expansion motions in G10.47+0.03 and detect velocity components in SgrB2-M at 50, 60, and 70 km/s relative to the Local Standard of Rest. The emission originates in regions of less than 0.1 pc diameter around the hypercompact HII regions G10.47+0.03 B1 and SgrB2-N K2, and reaches brightness temperatures of more than 200 K. Using the three-dimensional radiative transfer code RADMC-3D, we model the sources as dense dust cores heated by stars in the HII regions, and derive masses of hot (>300 K) molecular gas of more than 100 solar masses (for an HCN fractional abundance of 10$^{-5}$), challenging current simulations of massive star formation. Heating only by the stars in the HII regions is sufficient to produce such large quantities of hot molecular gas, provided that dust is optically thick to its own radiation, leading to high temperatures through diffusion of radiation.

Recent studies suggest that infrared dark clouds (IRDCs) have the potential of harboring the earliest stages of massive star formation and indeed evidence for this is found toward distinct regions within them. We present a study with the Plateau de Bure Interferometer of a core in the archetypal filamentary IRDC G11.11-0.12 at few arcsecond resolution to determine its physical and chemical structure. The data consist of continuum and line observations covering the C34S 2-1 line and the methanol 2_k-1_k v_t=0 lines at 3mm and the methanol 5_k-4_k v_t =0 lines at 1mm. Our observations show extended emission in the continuum at 1 and 3 mm. The methanol 2_k-1_k v_t=0 emission presents three maxima extending over 1 pc scale (when merged with single-dish short-spacing observations); one of the maxima is spatially coincident with the continuum emission. The fitting results show enhanced methanol fractional abundance (~3x10^-8) at the central peak with respect to the other two peaks, where it decreases by about an order of magnitude (~4-6x10^-9). Evidence of extended 4.5 microns emission, "wings" in the CH3OH 2_k-1_k spectra, and CH3OH abundance enhancement point to the presence of an outflow in the East-West direction. In addition, we find a gradient of ~4 km/s in the same direction, which we interpret as being produced by an outflow(s)-cloud interaction.

We used a large database extracted from Galaxy-wide surveys of H {\scriptsize I}, $^{13}$CO 1-0, $8 \micron$ and $870 \micron$ continuum to trace diffuse atomic gas, low- to medium-density molecular gas, high-density molecular gas, and star formation activity which we complemented by dedicated $^{12}$CO 2--1, 3--2 observations of the W43 region. From the detailed 3D (space-space-velocity) analysis of the molecular and atomic cloud tracers through the region and despite its wide velocity range (\emph{FWHM}$\sim22.3 \kms$ around $\vlsr\sim 95.9 \kms$), we identified W43 as a large (equivalent diameter $\sim 140$ pc) and coherent complex of molecular clouds which is surrounded by an atomic gas envelope (equivalent diameter $\sim 290$ pc). We measured the total mass of this newly-identified molecular complex ($M_{{\tiny total}}\sim 7.1 \times 10^6 \msun$), the mass contained in dense $870 \micron$ clumps ($<5$ pc dense cloud structures, $M_{{\tiny clumps}}\sim 8.4 \times 10^5 \msun$) and conclude that W43 is particularly massive and concentrated. The distance we assume for the W43 complex is 6 kpc from the Sun, which may place it at the meeting point of the Scutum-Centaurus (or Scutum-Crux) Galactic arm and the Bar, a dynamically complex region where high-velocity streams could easily collide. The star formation rate of W43 is suggested not to be steady but it is increasing from $\, \sim 0.01 \msun\, {yr}^{-1}$ (measured from its $8\,\mu$m luminosity) to $\sim 0.1 \msun\, {yr}^{-1}$ (measured from its molecular content). From the global properties of W43, we claim it is an extreme molecular complex in the Milky Way and it could even be forming starburst clusters in the near future.

Hydrides are key ingredients of interstellar chemistry since they are the initial products of chemical networks that lead to the formation of more complex molecules. The fundamental rotational transitions of light hydrides fall into the submillimeter wavelength range. Using the APEX telescope, we observed the long sought hydrides SH+ and OH+ in absorption against the strong continuum source Sagittarius B2(M). Both, absorption from Galactic center gas as well as from diffuse clouds in intervening spiral arms over a large velocity range is observed. The detected absorption of a continuous velocity range on the line of sight shows these hydrides to be an abundant component of diffuse clouds. In addition, we used the strongest submillimeter dust continuum sources in the inner Galaxy to serve as background candles for a systematic census of these hydrides in diffuse clouds and massive star forming regions of our Galaxy and initial results of this survey are presented.

This paper reviews the first results of observations of H2O line emission with Herschel-HIFI towards high-mass star-forming regions, obtained within the WISH guaranteed time program. The data reveal three kinds of gas-phase H2O: `cloud water' in cold tenuous foreground clouds, `envelope water' in dense protostellar envelopes, and `outflow water' in protostellar outflows. The low H2O abundance (1e-10 -- 1e-9) in foreground clouds and protostellar envelopes is due to rapid photodissociation and freeze-out on dust grains, respectively. The outflows show higher H2O abundances (1e-7 -- 1e-6) due to grain mantle evaporation and (probably) neutral-neutral reactions.

`Water In Star-forming regions with Herschel' (WISH) is a key program on the Herschel Space Observatory designed to probe the physical and chemical structure of young stellar objects using water and related molecules and to follow the water abundance from collapsing clouds to planet-forming disks. About 80 sources are targeted covering a wide range of luminosities and evolutionary stages, from cold pre-stellar cores to warm protostellar envelopes and outflows to disks around young stars. Both the HIFI and PACS instruments are used to observe a variety of lines of H2O, H218O and chemically related species. An overview of the scientific motivation and observational strategy of the program is given together with the modeling approach and analysis tools that have been developed. Initial science results are presented. These include a lack of water in cold gas at abundances that are lower than most predictions, strong water emission from shocks in protostellar environments, the importance of UV radiation in heating the gas along outflow walls across the full range of luminosities, and surprisingly widespread detection of the chemically related hydrides OH+ and H2O+ in outflows and foreground gas. Quantitative estimates of the energy budget indicate that H2O is generally not the dominant coolant in the warm dense gas associated with protostars. Very deep limits on the cold gaseous water reservoir in the outer regions of protoplanetary disks are obtained which have profound implications for our understanding of grain growth and mixing in disks.

W49A is a giant molecular cloud which harbors some of the most luminous embedded clusters in the Galaxy. However, the explanation for this starburst-like phenomenon is still under debate. Methods. We investigated large-scale Spitzer mid-infrared images together with a Galatic Ring Survey 13CO J = 1-0 image, complemented with higher resolution (~ 11 arcsec) 13CO J = 2-1 and C18O J = 2-1 images over a ~ 15 x 13 pc^2 field obtained with the IRAM 30m telescope. Two expanding shells have been identified in the mid-infrared images, and confirmed in the position-velocity diagrams made from the 13CO J = 2-1 and C18O J = 2-1 data. The mass of the averaged expanding shell, which has an inner radius of ~ 3.3 pc and a thickness of ~ 0.41 pc, is about 1.9 x 10^4 M*. The total kinetic energy of the expanding shells is estimated to be ~ 10^49 erg which is probably provided by a few massive stars, whose radiation pressure and/or strong stellar winds drive the shells. The expanding shells are likely to have a common origin close to the two ultracompact Hii regions (source O and source N), and their expansion speed is estimated to be ~ 5 km/s, resulting in an age of ~ 3-7 x 10^5 years. In addition, on larger (~ 35 x 50 pc^2) scales, remnants of two gas ejections have been identified in the 13CO J = 1 - 0 data. Both ejections seem to have the same center as the expanding shells with a total energy of a few times 10^50 erg. The main driving mechanism for the gas ejections is unclear, but likely related to the mechanism which triggers the starburst in W49A.

Based on color-selected IRAS point sources, we have started to conduct a survey of 47 high-mass star-forming regions in the southern hemisphere in 870um dust continuum and molecular line emission in several frequency ranges between 290 GHz and 806 GHz. This paper describes the pilot study of the three sources IRAS12326-6245, IRAS16060-5146, and IRAS16065-5158. To characterize the physical and chemical properties of southern massive star-forming regions, the three high-luminosity southern hemisphere hot cores were observed with APEX in five frequency setups aimed at groups of lines from the following molecules: CH3OH, H2CO, and CH3CN. Using the LTE approximation, temperatures, source sizes, and column densities were determined through modeling of synthetic spectra with the XCLASS program. Dust continuum observations were done with the Large APEX BOlometer CAmera (LABOCA) at 870um and the 3mm continuum was imaged with the ATCA. Based on the detection of high-excitation CH3CN lines and lines from complex organic species, the three sources are classified as line rich, hot core type sources. For all three, the modeling indicates that the line emission emerges from a combination of an extended, cooler envelope, and a hot compact component. All three sources show an overabundance of oxygen-bearing species compared to nitrogen-bearing species. Based on the results obtained in the three sources, which served as templates for the survey, the most promising (and feasible) frequency setups for the remaining 44 sources were decided upon.

We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

We have used the Atacama Pathfinder Experiment 12 m telescope (APEX) to carry out an absorption study of submillimeter wavelength rotational ground-state lines of H35Cl, H37Cl, 13CH+, and, for the first time, of the SH+ radical (sulfoniumylidene or sulfanylium). We detected the quartet of ground-state hyperfine structure lines of SH+ near 683 GHz with the CHAMP+ array receiver against the strong continuum source Sagittarius B2, which is located close to the center of our Galaxy. In addition to absorption from various kinematic components of Galactic center gas, we also see absorption at the radial velocities belonging to intervening spiral arms. This demonstrates that SH+ is a ubiquitous component of the diffuse interstellar medium. We do not find clear evidence for other SH+ lines we searched for, which is partially due to blending with lines from other molecules. In addition to SH+, we observed absorption from H35Cl, H37Cl, and 13CH+. The observed submillimeter absorption is compared in detail with absorption in 3 mm transitions of H13CO+ and c-C3H2 and the CO J = 1-0 and 3-2 transitions.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

We attempt to determine the nature of the bubbles observed by Spitzer in the Galactic plane, mainly to establish if possible their association with massive stars. We take advantage of the very simple morphology of these objects to search for star formation triggered by HII regions, and to estimate the importance of this mode of star formation. We consider a sample of 102 bubbles detected by Spitzer-GLIMPSE, and catalogued by Churchwell et al.(2006). We use mid-infrared and radio-continuum public data to discuss their nature. We use the ATLASGAL survey at 870 micron to search for dense neutral material collected on their borders. Results: We find that 86% of the bubbles contain ionized gas detected by means of its radio-continuum emission at 20-cm. Thus, most of the bubbles observed at 8.0 micron enclose HII regions ionized by O-B2 stars. Ninety-eight percent of the bubbles exhibit 24 micron emission in their central regions. The ionized regions at the center of the 8.0 micron bubbles seem to be devoid of PAHs but contain hot dust. Among the 65 regions for which the angular resolution of the observations is high enough to resolve the spatial distribution of cold dust at 870 micron, we find that 40% are surrounded by cold dust, and that another 28% contain interacting condensations. The former are good candidates for the collect and collapse process, as they display an accumulation of dense material at their borders. The latter are good candidates for the compression of pre-existing condensations by the ionized gas. Eighteen bubbles exhibit associated ultracompact HII regions and/or methanol masers in the direction of dust condensations adjacent to their ionization fronts. Our results suggest that more than a quarter of the bubbles may have triggered the formation of massive objects.

Due to the recent dramatic technological advances, infrared interferometry can now be applied to new classes of objects, resulting in exciting new science prospects, for instance, in the area of high-mass star formation. Although extensively studied at various wavelengths, the process through which massive stars form is still only poorly understood. For instance, it has been proposed that massive stars might form like low-mass stars by mass accretion through a circumstellar disk/envelope, or otherwise by coalescence in a dense stellar cluster. After discussing the technological challenges which result from the special properties of these objects, we present first near-infrared interferometric observations, which we obtained on the massive YSO IRAS 13481-6124 using VLTI/AMBER infrared long-baseline interferometry and NTT speckle interferometry. From our extensive data set, we reconstruct a model-independent aperture synthesis image which shows an elongated structure with a size of 13x19 AU, consistent with a disk seen under an inclination of 45 degree. The measured wavelength-dependent visibilities and closure phases allow us to derive the radial disk temperature gradient and to detect a dust-free region inside of 9.5 AU from the star, revealing qualitative and quantitative similarities with the disks observed in low-mass star formation. In complementary mid-infrared Spitzer and sub-millimeter APEX imaging observations we detect two bow shocks and a molecular out ow which are oriented perpendicular to the disk plane and indicate the presence of a bipolar outflow emanating from the inner regions of the system.

Circumstellar disks are an essential ingredient of the formation of low-mass stars. It is unclear, however, whether the accretion-disk paradigm can also account for the formation of stars more massive than about 10 solar masses, in which strong radiation pressure might halt mass infall. Massive stars may form by stellar merging, although more recent theoretical investigations suggest that the radiative-pressure limit may be overcome by considering more complex, nonspherical infall geometries. Clear observational evidence, such as the detection of compact dusty disks around massive young stellar objects, is needed to identify unambiguously the formation mode of the most massive stars. Here we report near-infrared interferometric observations that spatially resolve the astronomical unit-scale distribution of hot material around a high-mass (approx. 20 solar masses) young stellar object. The image shows an elongated structure with a size of about 13 x 19 astronomical units, consistent with a disk seen at an inclination angle of 45 degree. Using geometric and detailed physical models, we found a radial temperature gradient in the disk, with a dust-free region less than 9.5 astronomical units from the star, qualitatively and quantitatively similar to the disks observed in low-mass star formation. Perpendicular to the disk plane we observed a molecular outflow and two bow shocks, indicating that a bipolar outflow emanates from the inner regions of the system.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predictions from conventional chemical models by many orders of magnitude. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen-rich stars using the NH3 ortho J_K = 1_0 - 0_0 ground-state line. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 J_K = 1_0 - 0_0 transition near 572.5 GHz, simultaneously with the ortho-H2O J_Ka,Kc = 1_1,0 -1_0,1 transition, toward VY CMa, OH 26.5+0.6, IRC+10420, and IK Tau. We conducted non-LTE radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the latter two stars, Very Large Array imaging of NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those obtained for the ground state H2O line. The NH3 abundances relative to H2 are very high and range from 2 x 10-7 to 3 x 10-6 for the objects we have studied. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield satisfactory fits to the observed line profiles, it leads to abundance estimates that confirm the very high values found in earlier studies. New ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

We derive the dense core structure and the water abundance in four massive star-forming regions which may help understand the earliest stages of massive star formation. We present Herschel-HIFI observations of the para-H2O 1_11-0_00 and 2_02-1_11 and the para-H2-18O 1_11-0_00 transitions. The envelope contribution to the line profiles is separated from contributions by outflows and foreground clouds. The envelope contribution is modelled using Monte-Carlo radiative transfer codes for dust and molecular lines (MC3D and RATRAN), with the water abundance and the turbulent velocity width as free parameters. While the outflows are mostly seen in emission in high-J lines, envelopes are seen in absorption in ground-state lines, which are almost saturated. The derived water abundances range from 5E-10 to 4E-8 in the outer envelopes. We detect cold clouds surrounding the protostar envelope, thanks to the very high quality of the Herschel-HIFI data and the unique ability of water to probe them. Several foreground clouds are also detected along the line of sight. The low H2O abundances in massive dense cores are in accordance with the expectation that high densities and low temperatures lead to freeze-out of water on dust grains. The spread in abundance values is not clearly linked to physical properties of the sources.

The Heterodyne Instrument for the Far Infrared (HIFI) onboard the Herschel Space Observatory allows the first observations of light diatomic molecules at high spectral resolution and in multiple transitions. Here, we report deep integrations using HIFI in different lines of hydrides towards the high-mass star forming region AFGL 2591. Detected are CH, CH+, NH, OH+, H2O+, while NH+ and SH+ have not been detected. All molecules except for CH and CH+ are seen in absorption with low excitation temperatures and at velocities different from the systemic velocity of the protostellar envelope. Surprisingly, the CH(JF,P = 3/2_2,- - 1/2_1,+) and CH+(J = 1 - 0, J = 2 - 1) lines are detected in emission at the systemic velocity. We can assign the absorption features to a foreground cloud and an outflow lobe, while the CH and CH+ emission stems from the envelope. The observed abundance and excitation of CH and CH+ can be explained in the scenario of FUV irradiated outflow walls, where a cavity etched out by the outflow allows protostellar FUV photons to irradiate and heat the envelope at larger distances driving the chemical reactions that produce these molecules.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

Water in Star-forming regions with Herschel (WISH) is a Herschel Key Program investigating the water chemistry in young stellar objects (YSOs) during protostellar evolution. Hydroxyl (OH) is one of the reactants in the chemical network most closely linked to the formation and destruction of H2O. High-temperature chemistry connects OH and H2O through the OH + H2 <-> H2O + H reactions. Formation of H2O from OH is efficient in the high-temperature regime found in shocks and the innermost part of protostellar envelopes. Moreover, in the presence of UV photons, OH can be produced from the photo-dissociation of H2O. High-resolution spectroscopy of the OH 163.12 micron triplet towards HH 46 and NGC 1333 IRAS 2A was carried out with the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel. The low- and intermediate-mass YSOs HH 46, TMR 1, IRAS 15398-3359, DK Cha, NGC 7129 FIRS 2, and NGC 1333 IRAS 2A were observed with the Photodetector Array Camera and Spectrometer (PACS) in four transitions of OH and two [OI] lines. The OH transitions at 79, 84, 119, and 163 micron and [OI] emission at 63 and 145 micron were detected with PACS towards the class I low-mass YSOs as well as the intermediate-mass and class I Herbig Ae sources. No OH emission was detected from the class 0 YSO NGC 1333 IRAS 2A, though the 119 micron was detected in absorption. With HIFI, the 163.12 micron was not detected from HH 46 and only tentatively detected from NGC 1333 IRAS 2A. The combination of the PACS and HIFI results for HH 46 constrains the line width (FWHM > 11 km/s) and indicates that the OH emission likely originates from shocked gas. This scenario is supported by trends of the OH flux increasing with the [OI] flux and the bolometric luminosity. Similar OH line ratios for most sources suggest that OH has comparable excitation temperatures despite the different physical properties of the sources.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

We are carrying out a physical and chemical study of the protostellar envelopes in a representative sample of IM Class 0 protostars. In our first paper (Crimier et al. 2010), we determined the physical structure (density-temperature radial profiles) of the protostellar envelopes. Here, we study the CO depletion and N2H+ deuteration. We observed the millimeter lines of C18O, C17O, N2H+ and N2D+ toward the protostars using the IRAM 30m telescope. Based on these observations, we derived the C18O, N2H+ and N2D+ radial abundance profiles across their envelopes using a radiative transfer code. In addition, we modeled the chemistry of the protostellar envelopes. All the C18O 1-0 maps are well fit assuming that the C18O abundance decreases inwards within the protostellar envelope until the gas and dust reach the CO evaporation temperature, 20-25K, where the CO is released back to the gas phase. The N2H+ deuterium fractionation in Class 0 IMs is [N2D+]/[N2H+]=0.005-0.014, two orders of magnitude higher than the elemental [D/H] value in the interstellar medium, but a factor of 10 lower than in pre-stellar clumps. Chemical models account for the C18O and N2H+ observations if we assume the CO abundance is 2 times lower than the canonical value in the inner envelope. This could be the consequence of the CO being converted into CH3OH on the grain surfaces prior to the evaporation and/or the photodissociation of CO by the stellar UV radiation. The deuterium fractionation is not fitted by chemical models. This discrepancy is very likely caused by the simplicity of our model that assumes spherical geometry and neglects important phenomena like the effect of bipolar outflows and UV radiation from the star. More important, the deuterium fractionation is dependent on the ortho-to-para H2 ratio, which is not likely to reach the steady-state value in the dynamical time scales of these protostars.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

The Atacama Pathfinder Experiment (APEX) 12m telescope was used to observe the N=1-0, J=0-1 ground state transitions of OH+ at 909.1588 GHz with the CHAMP+ heterodyne array receiver. Two blended hyperfine structure transitions were detected in absorption against the strong continuum source Sagittarius B2(M) and in several pixels offset by 18". Both, absorption from Galactic center gas as well as absorption from diffuse clouds in intervening spiral arms in a velocity range from -116 to 38.5 km/s is observed. The total OH+ column density of absorbing gas is 2.4 \times 10^15 cm-2. A column density local to Sgr B2(M) of 2.6 \times 10^14 cm-2 is found. On the intervening line-of-sight the column density per unit velocity interval are in the range from 1 to 40 \times 10^12 cm-2/(km/s). OH+ is found to be on average more abundant than other hydrides such as SH+ and CH+. Abundance ratios of OH and atomic oxygen to OH+ are found to be in the range of 10^1-2 and 10^3-4, respectively. The detected absorption of a continuous velocity range on the line-of-sight shows OH+ to be an abundant component of diffuse clouds.

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

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

Water is a key molecule in the star formation process, but its spatial distribution in star-forming regions is not well known. We study the distribution of dust continuum and H2O and 13CO line emission in DR21, a luminous star-forming region with a powerful outflow and a compact HII region. Herschel-HIFI spectra near 1100 GHz show narrow 13CO 10-9 emission and H2O 1(11)-0(00) absorption from the dense core and broad emission from the outflow in both lines. The H2O line also shows absorption by a foreground cloud known from ground-based observations of low-J CO lines. The dust continuum emission is extended over 36" FWHM, while the 13CO and H2O lines are confined to ~24" or less. The foreground absorption appears to peak further North than the other components. Radiative transfer models indicate very low abundances of ~2e-10 for H2O and ~8e-7 for 13CO in the dense core, and higher H2O abundances of ~4e-9 in the foreground cloud and ~7e-7 in the outflow. The high H2O abundance in the warm outflow is probably due to the evaporation of water-rich icy grain mantles, while the H2O abundance is kept down by freeze-out in the dense core and by photodissociation in the foreground cloud.

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

The earliest phases of massive star formation are found in cold and dense infrared dark clouds (IRDCs). Since the detection method of IRDCs is very sensitive to the local properties of the background emission, we present here an alternative method to search for high column density in the Galactic plane by using infrared extinction maps. We find clouds between 1 and 5 kpc, of which many were missed by previous surveys. By studying the physical conditions of a subsample of these clouds, we aim at a better understanding of the initial conditions of massive star formation. We made extinction maps of the Galactic plane based on the 3.6-4.5 microns color excess between the two shortest wavelength Spitzer IRAC bands, reaching to visual extinctions of ~100 mag and column densities of 9x10^22 cm^-2. From this we compiled a new sample of cold and compact high extinction clouds. We used the MAMBO array at the IRAM 30m telescope to study the morphology, masses and densities of the clouds and the dense clumps within them. The latter were followed up by pointed ammonia observations with the 100m Effelsberg telescope, to determine rotational temperatures and kinematic distances. Extinction maps of the Galactic plane trace large scale structures such as the spiral arms. The 1.2 mm emission maps reveal that the high extinction clouds contain extended cold dust emission, from filamentary structures to still diffuse clouds. Most of the clouds are dark in 24 microns, but several show already signs of star formation via maser emission or bright infrared sources, suggesting that the high extinction clouds contain a variety of evolutionary stages. The observations suggest an evolutionary scheme from dark, cold and diffuse clouds, to clouds with a stronger 1.2 mm peak and to finally clouds with many strong 1.2 mm peaks, which are also warmer, more turbulent and already have some star formation signposts.

Context: Intermediate mass protostars provide a bridge between low- and high-mass protostars. Furthermore, they are an important component of the UV interstellar radiation field. Despite their relevance, little is known about their formation process. Aims: We present a systematic study of the physical structure of five intermediate mass, candidate Class 0 protostars. Our two goals are to shed light on the first phase of intermediate mass star formation and to compare these protostars with low- and high-mass sources. Methods: We derived the dust and gas temperature and density profiles of the sample. We analysed all existing continuum data on each source and modelled the resulting SED with the 1D radiative transfer code DUSTY. The gas temperature was then predicted by means of a modified version of the code CHT96. Results: We found that the density profiles of five out of six studied intermediate mass envelopes are consistent with the predictions of the "inside-out" collapse theory.We compared several physical parameters, like the power law index of the density profile, the size, the mass, the average density, the density at 1000 AU and the density at 10 K of the envelopes of low-, intermediate, and high-mass protostars. When considering these various physical parameters, the transition between the three groups appears smooth, suggesting that the formation processes and triggers do not substantially differ.

The aim of this paper is to study the molecular composition in the circumstellar envelope around the oxygen-rich star IK Tau. We observed IK Tau in several (sub)millimeter bands using the APEX telescope during three observing periods. To determine the spatial distribution of the $\mathrm{^{12}CO(3-2)}$ emission, mapping observations were performed. To constrain the physical conditions in the circumstellar envelope, multiple rotational CO emission lines were modeled using a non local thermodynamic equilibrium radiative transfer code. The rotational temperatures and the abundances of the other molecules were obtained assuming local thermodynamic equilibrium. An oxygen-rich Asymptotic Giant Branch star has been surveyed in the submillimeter wavelength range. Thirty four transitions of twelve molecular species, including maser lines, were detected. The kinetic temperature of the envelope was determined and the molecular abundance fractions of the molecules were estimated. The deduced molecular abundances were compared with observations and modeling from the literature and agree within a factor of 10, except for SO$_2$, which is found to be almost a factor 100 stronger than predicted by chemical models. From this study, we found that IK Tau is a good laboratory to study the conditions in circumstellar envelopes around oxygen-rich stars with (sub)millimeter-wavelength molecular lines. We could also expect from this study that the molecules in the circumstellar envelope can be explained more faithful by non-LTE analysis with lower and higher transition lines than by simple LTE analysis with only lower transition lines. In particular, the observed CO line profiles could be well reproduced by a simple expanding envelope model with a power law structure.

Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope around IK~Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also the line shapes, are used as diagnostic tool to study the envelope structure. Results: The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the envelope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance determination; the modeling of high-excitation SO$_2$ lines is cumbersome, possibly related to line misidentifications or problems with the collisional rates. The SiO isotopic ratios ($^{29}$SiO/$^{28}$SiO and $^{30}$SiO/$^{28}$SiO) point toward an enhancement in $^{28}$SiO compared to results of classical stellar evolution codes. Predictions for H$_2$O lines in the spectral range of the Herschel/HIFI mission are performed. [abbreviated]

Molecular outflows from high-mass young stellar objects provide an excellent way to study the star formation process, and investigate if they are scaled-up versions of their low-mass counterparts. We selected the nearby massive star forming region IRAS 17233-3606 in order to study the kinematics and physics along the molecular outflow(s) originating from this source. We observed IRAS 17233-3606 in CO, a typical tracer of gas associated with molecular outflow, with the Submillimeter Array in the (2-1) transition, and with the APEX telescope in the higher excitation (6-5) line. Additional infrared H2 observations were performed with the UKIRT telescope. The CO data were analysed using a LVG approach. Our data resolve the previously detected molecular outflow in at least three different components, one of them with a high collimation factor ~4, and characterised by emission at extremely high velocities (|v-v_{LSR}|>120 km s^{-1}). The estimate of the kinematical outflow parameters are typical of massive YSOs, and in agreement with the measured bolometric luminosity of the source. The kinematic ages of the flows are in the range 10^2-10^3 yr, and therefore point to young objects that still did not reach the main sequence.

The origin and heating mechanisms of warm (50<T<200 K) molecular gas in low-mass young stellar objects (YSOs) are strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by shocks and by UV associated with the outflows or accretion have been proposed. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, bipolar outflows and the quiescent cloud. High-J CO maps (12CO J=6--5 and 7--6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of the 13CO J=6--5 transition and [C I] 3P_2-3P_1 line were taken. Strong quiescent narrow-line 12CO 6--5 and 7--6 emission is seen toward all protostars. In the case of HH~46 and Ced 110 IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. Shock-heated warm gas is only detected for Class 0 flows and the more massive Class I sources such as HH~46. Outflow temperatures, estimated from the CO 6--5 and 3--2 line wings, are ~100 K, close to model predictions, with the exception of the L~1551 IRS 5 and IRAS 12496-7650, for which temperatures <50 K are found. APEX-CHAMP+ is uniquely suited to directly probe a protostar's feedback on its accreting envelope gas in terms of heating, photodissociation, and outflow dispersal by mapping 1'x1' regions in high-J CO and [C I] lines.

Aims. Expanding HII regions and propagating shocks are common in the environment of young high-mass star-forming complexes. They can compress a pre-existing molecular cloud and trigger the formation of dense cores. We investigate whether these phenomena can explain the formation of high-mass protostars within an infrared dark cloud located at the position of G327.3-0.6 in the Galactic plane, in between two large infrared bubbles and two HII regions. Methods: The region of G327.3-0.6 was imaged at 450 ? m with the CEA P-ArT\'eMiS bolometer array on the Atacama Pathfinder EXperiment telescope in Chile. APEX/LABOCA and APEX-2A, and Spitzer/IRAC and MIPS archives data were used in this study. Results: Ten massive cores were detected in the P-ArT\'eMiS image, embedded within the infrared dark cloud seen in absorption at both 8 and 24 ?m. Their luminosities and masses indicate that they form high-mass stars. The kinematical study of the region suggests that the infrared bubbles expand toward the infrared dark cloud. Conclusions: Under the influence of expanding bubbles, star formation occurs in the infrared dark areas at the border of HII regions and infrared bubbles.

AIMS Our aim is to characterize the size, mass, density and temperature profiles of the protostellar envelope of HH~46 IRS 1 and its surrounding cloud material as well as the effect the outflow has on its environment.METHODS The CHAMP+ and LABOCA arrays on the APEX telescope, combined with lower frequency line receivers, are used to obtain a large continuum map and smaller heterodyne maps in various isotopologues of CO and HCO+. The high-J lines of CO (6--5 and 7--6) and its isotopologues together with [C I] 2--1, observed with CHAMP+, are used to probe the warm molecular gas in the inner few hundred AU and in the outflowing gas. The data are interpreted with continuum and line radiative transfer models. RESULTS Broad outflow wings are seen in CO low- and high-J lines at several positions, constraining the gas temperatures to a constant value of ~100 K along the red outflow axis and to ~60 K for the blue outflow. The derived outflow mass is of order 0.4--0.8 M_sol, significantly higher than previously found. The bulk of the strong high-J CO line emission has a surprisingly narrow width, however, even at outflow positions. These lines cannot be fit by a passively heated model of the HH 46 IRS envelope. We propose that it originates from photon heating of the outflow cavity walls by ultraviolet photons originating in outflow shocks and the accretion disk boundary layers. At the position of the bow shock itself, the UV photons are energetic enough to dissociate CO. The envelope mass of ~5 M_sol is strongly concentrated towards HH 46 IRS with a density power law of -1.8.

(Abridged) Studying continuum emission from interstellar dust is essential to locating and characterizing the highest density regions in the interstellar medium. In particular, the early stages of massive star formation remain poorly understood. Our goal is to produce a large-scale, systematic database of massive pre- and proto-stellar clumps in the Galaxy, to understand how and under what conditions star formation takes place. A well characterized sample of star-forming sites will deliver an evolutionary sequence and a mass function of high-mass, star-forming clumps. This systematic survey at submm wavelengths also represents a preparatory work for Herschel and ALMA. The APEX telescope is ideally located to observe the inner Milky Way. The Large APEX Bolometer Camera (LABOCA) is a 295-element bolometer array observing at 870 microns, with a beam size of 19". Taking advantage of its large field of view (11.4') and excellent sensitivity, we started an unbiased survey of the Galactic Plane, with a noise level of 50-70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). As a first step, we covered 95 sq. deg. These data reveal 6000 compact sources brighter than 0.25 Jy, as well as extended structures, many of them filamentary. About two thirds of the compact sources have no bright infrared counterpart, and some of them are likely to correspond to the precursors of (high-mass) proto-stars or proto-clusters. Other compact sources harbor hot cores, compact HII regions or young embedded clusters. Assuming a typical distance of 5 kpc, most sources are clumps smaller than 1 pc with masses from a few 10 to a few 100 M_sun. In this introductory paper, we show preliminary results from these ongoing observations, and discuss the perspectives of the survey.

NGC 2024, a sites of massive star formation, have complex internal structures caused by cal heating by young stars, outflows, and stellar winds. These complex cloud structures lead to intricate emission line shapes. The goal of this paper is to show that the complex line shapes of 12 CO lines in NGC 2024 can be explained consistently with a model, whose temperature and velocity structure are based on the well-established scenario of a PDR and the Blister model. We present velocity-resolved spectra of seven CO lines ranging from J=3 to J=13, and we combined these data with CO high-frequency data from the ISO satellite. We find that the bulk of the molecular cloud associated with NGC 2024 consists of warm (75 K) and dense (9e5 cm-3) gas. An additional hot (~ 300 K) component, located at the interface of the HII region and the molecular cloud, is needed to explain the emission of the high-J CO lines. Deep absorption notches indicate that very cold material (20 K) exists in front of the warm material, too. A temperature and column density structure consistent with those predicted by PDR models, combined with the velocity structure of a Blister model, appropriately describes the observed emission line profiles of this massive star forming region. This case study of NGC 2024 shows that, with physical insights into these complex regions and careful modeling, multi-line observations of CO can be used to derive detailed physical conditions in massive star forming regions.

We present molecular line observations of the massive star forming region IRAS 17233-3606 aimed at studying the molecular core associated with the source. The observations were made using the Atacama Pathfinder Experiment telescope in the CO (3-2) and HCO^+ (4-3) transitions, and in the CH_3OH (6_K-5_K), (7_K-6_K) and CH_3CN (16_K-15_K) bands. For the CO(3-2) and HCO^+ (4-3) transitions, we obtained maps with a size of 70''\times 70''. The typical angular resolution of the data is ~18''. Our observations reveal an exceptionally rich molecular spectrum, a signpost of hot core activity. Comparisons with two other prominent southern hot cores were made through observations in the same frequency setups. We also detected a bipolar outflow in CO (3-2) and HCO^+ (4-3) lines. Modelling reveals a hot core of size ~3'' and a temperature of 150 K in the IRAS17233-3606 region. The parameters of the molecular outflow are derived through the analysis of the CO (3-2) emission, and are typical of outflows driven by high-mass young stellar objects.

In this contribution, our knowledge of the initial conditions under which massive star formation takes place is reviewed. Massive stars are born in massive clumps of giant molecular clouds (GMCs), hence first the properties of GMCs are summarized. As a potentially early stage of molecular clouds, infrared dark clouds have been discovered a decade ago as dark patches in mid-infrared (MIR) images of the Galactic plane and many studies of the physical conditions within them have been conducted recently. Without the guidance of MIR absorption, large scale, unbiased cold dust surveys can be used as well to identify massive cold clumps. In the absence of indicators of ongoing massive star formation, like compact HII regions and bright IR sources, these clumps are the most promising objects for the study of the initial conditions of massive star formation. Current observational approaches to find IR quiet clumps and their physical and chemical properties are summarized.

We present APEX observations of HCN (9-8) and (4-3) lines toward a sample of hot cores. The spectral shapes of the main transitions are asymmetric and self-absorbed, as expected for high optical depth in a possibly infalling envelope. For spherical symmetry, the large column densities of these sources would mean that the central region is obscured by dust above a certain frequency. However, we detected the vibrationally excited satellite lines (v_2=1; J=9-8) at 797 GHz, which originate from the inner regions. This indicates that high-frequency ALMA observations of hot core centers will be feasible.

We present observations of the \nh3 (J,K) = (1,1) and (2,2) inversion transitions toward the infrared dark cloud G28.34+0.06, using the Very Large Array. Strong NH3 emission is found to coincide well with the infrared absorption feature in this cloud. The northern region of G28.34+0.06 is dominated by a compact clump (P2) with a high rotation temperature (29 K), large line width (4.3 km s$^{-1}$), and is associated with strong water maser (240 Jy) and a 24 $\mu$m point source with far IR luminosity of $10^3$ \lsun. We infer that P2 has embedded massive protostars although it lies in the 8 $\mu$m absorption region. The southern region has filamentary structures. The rotation temperature in the southern region decreases with the increase of the integrated NH3 intensity, which indicates an absence of strong internal heating in these clumps. In addition, the compact core P1 in the south has small line width (1.2 km s$^{-1}$) surrounded by extended emission with larger line width (1.8 km s$^{-1}$), which suggests a dissipation of turbulence in the dense part of the cloud. Thus, we suggest that P1 is at a much earlier evolutionary stage than P2, possibly at a stage that begins to form a cluster with massive stars.

Since most high- and intermediate-mass protostars are at great distance and form in clusters, high linear resolution observations are needed to investigate their physical properties. To study the gas in the innermost region around the protostars in the proto-cluster IRAS 05358+3543, we observed the source in several transitions of methanol and other molecular species with the Plateau de Bure Interferometer and the Submillimeter Array, reaching a linear resolution of 1100 AU. We determine the kinetic temperature of the gas around the protostars through an LVG and LTE analysis of their molecular emission; the column densities of CH3OH, CH3CN and SO2 are also derived. Constrains on the density of the gas are estimated for two of the protostellar cores. We find that the dust condensations are in various evolutionary stages. The powerhouse of the cluster, mm1a, harbours a hot core with T~220 (75<T<330) K. A double-peaked profile is detected in several transitions toward mm1a, and we found a velocity gradient along a linear structure which could be perpendicular to one of the outflows from the vicinity of mm1a. Since the size of the double-peaked emission is less than 1100 AU, we suggest that mm1a might host a massive circumstellar disk. The other sources are in earlier stages of star formation. The least active source, mm3, could be a starless massive core, since it is cold (T<20 K), with a large reservoir of accreting material (M ~ 19 M_\odot), but no molecular emission peaks on it.

NGC6334I and I(N) have been observed with the Swedish-ESO Submillimetre Telescope, SEST, at wavelengths of 3, 2, and 1.3 mm. Especially NGC6334 I shows rich emission from many different molecules, comparable in line density to prototypical hot cores such as Orion-KL and SgrB2(N). In addition, a 4' by 4' region enfolding NGC 6334 I and I(N) has been mapped at a wavelength of 3 mm (75 to 116 GHz) with the Mopra telescope.

A systematic study of a large sample of sources, covering a wide range in galactocentric distances, masses and luminosities, is a fast and efficient way of obtaining a good overview of the different stages of high-mass star formation. With these goals in mind, we have started a survey of 40 color selected IRAS sources south of -20 degrees declination with the APEX telescope on Chajnantor, Chile. Our first APEX results already demonstrate that the selection criteria were successful, since some of the sources are very rich in molecular lines.

We investigate the mm-morphology of IC1396N with unprecedented spatial resolution to analyze its dust and molecular gas properties, and draw comparisons with objects of similar mass. We have carried out sensitive observations in the most extended configurations of the IRAM Plateau de Bure interferometer, to map the thermal dust emission at 3.3 and 1.3mm, and the emission from the $J$=13$_k\to12_k$ hyperfine transitions of methyl cyanide (CH$_3$CN). We unveil the existence of a sub-cluster of hot cores in IC1396N, distributed in a direction perpendicular to the emanating outflow. The cores are embedded in a common envelope of extended and diffuse dust emission. We find striking differences in the dust properties of the cores ($\beta\simeq$ 0) and the surrounding envelope ($\beta\simeq$ 1), very likely testifying to differences in the formation and processing of dust material. The CH$_3$CN emission peaks towards the most massive hot core and is marginally extended in the outflow direction.

We have observed several emission lines of two Nitrogen-bearing (C2H5CN and C2H3CN) and two Oxygen-bearing (CH3OCH3 and HCOOCH3) molecules towards a sample of well-known hot molecular cores (HMCs) in order to check whether the chemical differentiation seen in the Orion-HMC and W3(H_2O) between O- and N-bearing molecules is a general property of HMCs. With the IRAM-30m telescope we have observed 12 HMCs in 21 bands, centered at frequencies from 86250 to 258280 MHz. The rotational temperatures obtained range from ~100 to ~150 K in these HMCs. Single Gaussian fits performed to unblended lines show a marginal difference in the line peak velocities of the C2H5CN and CH3OCH3 lines, indicating a possible spatial separation between the region traced by the two molecules. On the other hand, neither the linewidths nor the rotational temperatures and column densities confirm such a result. By comparing the abundance ratio of the pair C2H5CN/C2H3CN with the predictions of theoretical models, we derive that the age of our cores ranges between 3.7 and 5.9x10^{4} yrs. The abundances of C2H5CN and C2H3CN are strongly correlated, as expected from theory which predicts that C2H3CN is formed through gas phase reactions involving C2H5CN. A correlation is also found between the abundances of C2H3CN and CH3OCH3, and C2H5CN and CH3OCH3. In all tracers the fractional abundances increase with the H_2 column density while they are not correlated with the gas temperature.

We performed interferometric observations of a high-mass protostellar candidate IRAS 18566+0408 in the \nh3 (J,K)=(1,1), (2,2) and (3,3) inversion transitions, the SiO J=2-1 and HCN J=1-0 lines, and the 43 and 87 GHz continuum emission using the VLA and OVRO. The 87 GHz continuum emission reveals two continuum peaks MM-1 and MM-2 along a molecular ridge. The dominant peak MM-1 coincides with a compact emission feature at 43 GHz, and arises mostly from the dust emission. For dust emissivity index $\beta$ of 1.3, the masses in the dust peaks amount to 70 \msun for MM-1, and 27 \msun for MM-2. Assuming internal heating, the central luminosities of MM-1 and MM-2 are $6 \times 10^4$ and $8 \times 10^3$ \lsun, respectively. The SiO emission reveals a well collimated outflow emanating from MM-1. The jet-like outflow is also detected in \nh3 at velocities similar to the SiO emission. The outflow, with a mass of 27 \msun, causes significant heating in the gas to temperatures of 70 K, much higher than the temperature of $\lsim 15$ K in the extended core. Compact ($< 3''$) and narrow line ($<1.5$ \kms-1) \nh3 (3,3) emission features are found associated with the outflow. They likely arise from weak population inversion in \nh3 similar to the maser emission. Toward MM-1, there is a compact \nh3 structure with a linewidth that increases from 5.5 \kms-1 FHWM measured at 3$''$ resolution to 8.7 \kms-1\ measured at 1$''$ resolution. This linewidth is much larger than the FWHM of $<$ 2 \kms-1 in the entire core, and does not appear to originate from the outflow. This large linewidth may arise from rotation/infall, or relative motions of unresolved protostellar cores.

The transition between the low density groups of T Tauri stars and the high density clusters around massive stars occurs in the intermediate-mass (IM) range (M$_*$$\sim$2--8 M$_\odot$). High spatial resolution studies of IM young stellar objects (YSO) can provide important clues to understand the clustering in massive star forming regions. Aims: Our aim is to search for clustering in IM Class 0 protostars. The high spatial resolution and sensitivity provided by the new A configuration of the Plateau de Bure Interferometer (PdBI) allow us to study the clustering in these nearby objects. Methods: We have imaged three IM Class 0 protostars (Serpens-FIRS 1, IC 1396 N, CB 3) in the continuum at 3.3 and 1.3mm using the PdBI. The sources have been selected with different luminosity to investigate the dependence of the clustering process on the luminosity of the source. Results: Only one millimeter (mm) source is detected towards the low luminosity source Serpens--FIRS 1. Towards CB 3 and IC1396 N, we detect two compact sources separated by $\sim$0.05 pc. The 1.3mm image of IC 1396 N, which provides the highest spatial resolution, reveal that one of these cores is splitted in, at least, three individual sources.

The young massive star-forming region IRAS05358+3543 was observed at high-spatial resolution in the continuum emission at 3.1 and 1.2mm with the Plateau de Bure Interferometer, and at 875 and 438mum with the Submillimeter Array. We resolve at least four continuum sub-sources that are likely of protostellar nature. Two of them are potentially part of a proto-binary system with a projected separation of 1700AU. Additional (sub)mm continuum peaks are not necessarily harboring protostars but maybe caused by the multiple molecular outflows. The spectral energy distributions (SEDs) of the sub-sources show several features. The main power house mm1, which is associated with CH3OH maser emission, a hypercompact HII region and a mid-infrared source, exhibits a typical SED with a free-free emission component at cm and long mm wavelengths and a cold dust component in the (sub)mm part of the spectrum (spectral index between 1.2mm and 438mum alpha~3.6). The free-free emission corresponds to a Lyman continuum flux of an embedded 13Msun B1 star. The coldest source of the region, mm3, has alpha~3.7 between 1.2mm and 875mum, but has lower than expected fluxes in the shorter wavelength 438mum band. This turnover of the Planck-function sets an upper limit on the dust temperature of mm3 of approximately 20K. The uv-data analysis of the density structure of individual sub-cores reveals distributions with power-law indices between 1.5 and 2. This resembles the density distributions of the larger-scale cluster-forming clump as well as those from typical low-mass cores.

UltraCompact HII regions are signposts of high mass star formation. Since high-mass star formation occurs in clusters, one expects to find even earlier phases of massive star formation in the vicinity of UltraCompact HII regions. Here, we study the amount of deuteration and depletion toward pre/protocluster clumps found in a wide-field (10 X 10 arcmin) census of clouds in 32 massive star-forming regions that are known to harbour UCHII regions. We find that 65% of the observed sources have strong NH2D emission and more than 50% of the sources exhibit a high degree of deuteration, (0.1 < NH2D/NH3 < 0.7), 0.7 being the highest observed deuteration of NH3 reported to date. Our search for NHD2 in two sources did not result in a detection. The enhancement in deuteration coincides with moderate CO depletion onto dust grains. There is no evidence of a correlation between the two processes, though an underlying correlation cannot be ruled out as the depletion factor is very likely to be only a lower limit. In summary, we find CO depletion and high deuteration towards cold cores in massive star forming regions. Therefore, these are good candidates for sources at the early phases of massive star formation. While our sensitive upper limits on NHD2 do not prove the predictions of the gas-phase and grain chemistry models wrong, an enhancement of ~10^4 over the cosmic D/H ratio from NH2D warrants explanation.

Unbiased molecular line surveys are a powerful tool for analyzing the physical and chemical parameters of astronomical objects and are the only means for obtaining a complete view of the molecular inventory for a given source. The present work stands for the first such investigation of a photon-dominated region. The first results of an ongoing millimeter-wave survey obtained towards the Orion Bar are reported. The APEX telescope in combination with the APEX-2A facility receiver was employed in this investigation. We derived the physical parameters of the gas through LVG analyses of the methanol and formaldehyde data. Information on the sulfur and deuterium chemistry of photon-dominated regions is obtained from detections of several sulfur-bearing molecules and DCN.

We observed the high velocity outflow originating from NGC6334I in several CO transitions with the APEX telescope, with the goal of deriving the physical parameters of the gas. Using an LVG analysis, we studied line ratios between the CO(3-2), CO(4-3), and CO(7-6) data as a function of the density and of the kinetic temperature of the gas. An upper limit on the CO column density is derived by comparison with 13CO data. We constrained the temperature to be higher than 50 K and the H_2 density to values higher than n~10^4 cm^{-3} towards the peak position in the red lobe, while T>15 K and n>10^3 cm^{-3} are derived towards the peak position in the blue lobe. The physical parameters of the outflow, its mass and its energetics, have been computed using the temperatures derived from this analysis. We conclude that high kinetic temperatures are present in the outflow and traced by high excitation CO lines. Observations of high-J CO lines are thus needed to infer reliable values of the kinetic temperatures and of the other physical parameters in outflows.

The THz atmospheric windows centered at roughly 1.3 and 1.5~THz, contain numerous spectral lines of astronomical importance, including three high-J CO lines, the N+ line at 205 microns, and the ground transition of para-H2D+. The CO lines are tracers of hot (several 100K), dense gas; N+ is a cooling line of diffuse, ionized gas; the H2D+ line is a non-depleting tracer of cold (~20K), dense gas. As the THz lines benefit the study of diverse phenomena (from high-mass star-forming regions to the WIM to cold prestellar cores), we have built the CO N+ Deuterium Observations Receiver (CONDOR) to further explore the THz windows by ground-based observations. CONDOR was designed to be used at the Atacama Pathfinder EXperiment (APEX) and Stratospheric Observatory For Infrared Astronomy (SOFIA). CONDOR was installed at the APEX telescope and test observations were made to characterize the instrument. The combination of CONDOR on APEX successfully detected THz radiation from astronomical sources. CONDOR operated with typical Trec=1600K and spectral Allan variance times of 30s. CONDOR's first light observations of CO 13-12 emission from the hot core Orion FIR4 (= OMC1 South) revealed a narrow line with T(MB) = 210K and delta(V)=5.4km/s. A search for N+ emission from the ionization front of the Orion Bar resulted in a non-detection. The successful deployment of CONDOR at APEX demonstrates the potential for making observations at THz frequencies from ground-based facilities.

We present a submm study of the massive hot core G327.3-0.6 that constrains its physical parameters and environment. The APEX telescope was used to image CO and N2H+ emission, to observe lines from other molecules toward a hot and a cold molecular core, and to measure the continuum flux density of the hot core. In the C18O J=3-2 line, two clumps were found, one associated with the HII region G327.3-0.5 and the other associated with the hot core. An additional cold clump is found 30 arcsec (0.4 pc) northeast of the hot core in bright N2H+ emission. From the the continuum data, we calculate a mass of 420 Msol and a size of 0.1 pc for the hot core. A new, more accurate position of the hot core is reported, which allows the association of the core with a bright mid-infrared source. The luminosity of the hot core is estimated to be between 5 and 15 10^4 Lsol. This study revealed several different evolutionary stages of massive star formation in the G327.3-0.6 region.

A study of 12 ultracompact HII regions was conducted to probe the physical conditions and kinematics in the inner envelopes of the molecular clumps harboring them. The APEX telescope was used to observe the sources in the CO (4-3) and 13CO (8-7) lines. Line intensities were modeled with the RATRAN radiative transfer code using power laws for the density and temperature to describe the physical structure of the clumps. All sources were detected in both lines. The optically thick CO (4-3) line shows predominantly blue skewed profiles reminiscent of infall. Line intensities can be reproduced well using the physical structure of the clumps taken from the literature. The optically thick line profiles show that CO is a sensitive tracer of ongoing infall in the outer envelopes of clumps harboring ultracompact HII regions and hot molecular cores.

High-mass star-forming regions are known to have a rich molecular spectrum from many species. Some of the very highly excited lines are emitted from very hot and dense gas close to the central object(s). The physics and chemistry of the inner cores of two high mass star forming regions, NGC6334(I) and G327.3-0.6, shall be characterized. Submillimeter line surveys with the APEX telescope provide spectra which sample many molecular lines at high excitation stages. Partial spectral surveys were obtained, the lines were identified, physical parameters were determined through fitting of the spectra. Both sources show similar spectra that are comparable to that of the only other high mass star forming region ever surveyed in this frequency range}, Orion-KL, but with an even higher line density. Evidence for very compact, very hot sources is found.

Context : Highly-collimated outflows are believed to be the earliest stage in outflow evolution, so their study is essential for understanding the processes driving outflows. The BHR71 Bok globule is known to harbour such a highly-collimated outflow, which is powered by a protostar belonging to a protobinary system. Aims : We aimed at investigating the interaction of collimated outflows with the ambient molecular cloud by using molecular tracers. Methods : We mapped the BHR71 highly-collimated outflow in CO(3-2) with the APEX telescope, and observed several bright points of the outflow in the molecular transitions CO(4-3), 13CO(3-2), C18O(3-2), and CH3OH(7-6). We use an LVG code to characterise the temperature enhancements in these regions. Results : In our CO(3-2) map, the second outflow driven by IRS2, which is the second source of the binary system, is completely revealed and shown to be bipolar. We also measure temperature enhancements in the lobes. The CO and methanol LVG modelling points to temperatures between 30 and 50K in the IRS1 outflow, while the IRS2 outflow seems to be warmer (up to 300K).

APECS is the distributed control system of the new Atacama Pathfinder EXperiment (APEX) telescope located on the Llano de Chajnantor at an altitude of 5107 m in the Atacama desert in northern Chile. APECS is based on Atacama Large Millimeter Array (ALMA) software and employs a modern, object-oriented design using the Common Object Request Broker Architecture (CORBA) as the middleware. New generic device interfaces simplify adding instruments to the control system. The Python based observer command scripting language allows using many existing software libraries and facilitates creating more complex observing modes. A new self-descriptive raw data format (Multi-Beam FITS or MBFITS) has been defined to store the multi-beam, multi-frequency data. APECS provides an online pipeline for initial calibration, observer feedback and a quick-look display. APECS is being used for regular science observations in local and remote mode since August 2005.

Infrared Dark Clouds appear to be the long sought population of cold and dense aggregations with the potential of harbouring the earliest stages of massive star formation. Up to now there has been no systematic study on the temperature distribution, velocity fields, chemical and physical state toward this new cloud population. Knowing these properties is crucial for understanding the presence, absence and the very potential of star formation. The present paper aims at addressing these questions. We analyse temperature structures and velocity fields and gain information on their chemical evolution. The gas emission is remarkably coextensive with the extinction seen at infrared wavelengths and with the submillimeter dust emission. Our results show that IRDCs are on average cold (T < 20 K) and have variations among the different cores. IRDC cores are in virial equilibrium, are massive (M > 100 M_sun), highly turbulent (1 -- 3 km/s) and exhibit significant velocity structure (variations around 1 -- 2 km/s over the cloud). We find an increasing trend in temperature from IRDCs with high ammonia column density to high mass protostellar objects and hot core/Ultracompact Hii regions stages of early warm high-mass star formation while linewidths of IRDCs are smaller. On the basis of this sample, we infer that while active star formation is not yet pervasive in most IRDCs, local condensations might collapse in the future or have already begun forming stars.

We report detection of moderate to high-mass star formation in an infrared dark cloud (G11.11-0.12) where we discovered class II methanol and water maser emissions at 6.7 GHz and 22.2 GHz, respectively. We also observed the object in ammonia inversion transitions. Strong emission from the (3,3) line indicates a hot (~60 K) compact component associated with the maser emission. The line width of the hot component (4 km/s), as well as the methanol maser detection, are indicative of high mass star formation. To further constrain the physical parameters of the source, we derived the spectral energy distribution (SED) of the dust continuum by analysing data from the 2MASS survey, HIRAS, MSX, the Spitzer Space Telescope, and interferometric 3mm observations. The SED was modelled in a radiative transfer program: a) the stellar luminosity equals 1200 L_sun corresponding to a ZAMS star of 8 M_sun; b) the bulk of the envelope has a temperature of 19 K; c) the mass of the remnant protostellar cloud in an area 8x10^17 cm or 15 arcsec across amounts to 500M_sun, if assuming standard dust of the diffuse medium, and to about 60 M_sun, should the grains be fluffy and have ice mantles; d) the corresponding visual extinction towards the star is a few hundred magnitudes. The near IR data can be explained by scattering from tenuous material above a hypothetical disk. The class II methanol maser lines are spread out in velocity over 11 km/s. To explain the kinematics of the masing spots, we propose that they are located in a Kepler disk at a distance of about 250 AU. The dust temperatures there are around 150 K, high enough to evaporate methanol--containing ice mantles.

We report the identification of a sample of potential High-Mass Starless Cores (HMSCs). The cores were discovered by comparing images of the fields containing candidate High-Mass Protostellar Objects (HMPOs) at 1.2mm and mid-infrared (8.3um; MIR) wavelengths. While the HMPOs are detected at both wavelengths, several cores emitting at 1.2mm in the same fields show absorption or no emission at the MIR wavelength. We argue that the absorption is caused by cold dust. The estimated masses of a few 10^2Msun - 10^3 Msun and the lack of IR emission suggests that they may be massive cold cores in a pre-stellar phase, which could presumably form massive stars eventually. Ammonia (1,1) and (2,2) observations of the cores indicate smaller velocity dispersions and lower rotation temperatures compared to HMPOs and UCHII regions suggesting a quiescent pre-stellar stage. We propose that these newly discovered cores are good candidates for the HMSC stage in high-mass star-formation. This sample of cores will allow us to study the high-mass star and cluster formation processes at the earliest evolutionary stages.

An extensive search for the complex molecule trans-ethyl methyl ether towards several hot core regions has been performed. Using the IRAM 30m telescope and the SEST 15m we looked at several frequencies where trans-ethyl methyl ether has strong transitions, as well as lines which are particularly sensitive to the physical conditions in which the molecule can be found. We included G34.26, NGC6334(I), Orion KL, and W51e2 which have previously been proven to have a rich chemistry of complex molecules. Our observations cannot confirm the tentative Orion KL detection made by Charnley et al. (2001) within their stated column density limits, but we confirm the existence of the trans-ethyl methyl ether towards W51e2 with a column density of 2x10^14 cm-2. The dimethyl ether/methanol ratio of 0.6 as well as the newly found ethyl methyl ether/ethanol ratio of 0.13 indicate relative high abundances of ethers toward W51e2. Furthermore, the observation of ethyl methyl ether also confirms the importance of ethanol as a grain mantle constituent. We present new upper limits of around 8x10^13 cm-2 for the column densities of the molecule toward Orion KL, G34.26, NGC6334(I) and estimate the column density towards SgrB2(N) to be of the same order. The W51e2 observations are discussed in more detail.

We summarize the properties of Infrared Dark Clouds, massive, dense, and cool aggregations of interstellar gas and dust that are found througout the Galaxy in projection against the strong mid-infrared background. We describe their distribution and give an overview of their physical properties and chemistry. These objects appear to be the progenitors of high-mass stars and star clusters, but seem to be largely devoid of star formation, which however is taking place in localized spots.

Molecular line CN, CS and mm continuum observations of two intermediate- to high-mass star-forming regions - IRAS20293+3952 and IRAS19410+2336 - obtained with the Plateau de Bure Interferometer at high spatial resolution reveal interesting characteristics of the gas and dust emission. In spite of the expectation that the CN and CS morphology might closely follow the dense gas traced by the dust continuum, both molecules avoid the most central cores. Comparing the relative line strengths of various CN hyperfine components, this appears not to be an opacity effect but to be due to chemical and physical effects. The CN data also indicate enhanced emission toward the different molecular outflows in the region. Regarding CS, avoiding the central cores can be due to high optical depth, but the data also show that the CS emission is nearly always associated with the outflows of the region. Therefore, neither CS nor CN appear well suited for dense gas and disk studies in these two sources, and we recommend the use of different molecules for future massive disk studies. An analysis of the 1 and 3mm continuum fluxes toward IRAS20293+3952 reveals that the dust opacity index beta is lower than the canonical value of 2. Tentatively, we identify a decreasing gradient of beta from the edge of the core to the core center. This could be due to increasing optical depth toward the core center and/or grain growth within the densest cores and potential central disks. We detect 3mm continuum emission toward the collimated outflow emanating from IRAS20293+3952. The spectral index of alpha ~ 0.8 in this region is consistent with standard models for collimated ionized winds.

We present observations made with the Berkeley-Illinois-Maryland Association millimeter array of the H2S 2(2,0)-2(1,1) and C18O 2-1 transitions toward a sample of four hot molecular cores associated with ultracompact HII regions: G9.62+0.19, G10.47+0.03, G29.96-0.02 and G31.41+0.31. The angular resolution varies from 1.5 to 2.4 arcsec, corresponding to scales of ~0.06 pc at the distance of these sources. High-velocity wings characteristic of molecular outflows are detected toward all four sources in the H2S line. In two cases (G29.96 and G31.41) red- and blueshifted lobes are clearly defined and spatially separate, indicating that the flows are collimated. We also confirm the previous detection of the outflow in G9.62F. Although the gas-phase H2S abundance is not well constrained, assuming a value of 10^-7 yields lower limits to total outflow masses of ~8 Msun, values which are consistent with the driving sources being massive protostars. Linear velocity gradients are detected in both C18O and H2S across G9.62, G29.96 and, to a lesser extent, G31.41. These gradients are observed to be at a different position angle to the outflow in G9.62F and G29.96, suggestive of a rotation signature in these two hot cores. Our observations show that these hot cores contain embedded massive protostellar objects which are driving bipolar outflows. Furthermore, the lack of strong centimeter-wave emission toward the outflow centers in G29.96 and G31.41 indicates that the outflow phase begins prior to the formation of a detectable ultracompact HII region.

We present subarcsecond observations at 2.7 and 1.4 mm of a sample of massive young stellar objects made with the BIMA millimetre array. For most sources the continuum emission on the smallest scales at 2.7 mm is dominated by free-free emission from the stellar wind or jet. Strong emission at 1.4 mm shows the presence of significant dust associated with Cep A and GL 490 but our resolution is not sufficient to resolve any structure. The 2.7-mm emission from GL 490 is resolved but it is not clear whether we are seeing a single circumstellar disk or a secondary companion, although near-infrared data support the disk hypothesis. Estimates of the dust mass yield values of ~1--4 Msun within radii of 150 to 1000 AU.

Millimeter and mid-infrared observations have been made of the dense clumps of dust and gas and of young stellar objects (YSOs) associated with the bright, compact submillimeter source G79.3+0.3 P1 in the relatively nearby MSX infrared-dark cloud G79.3+0.3. The Gemini mid-infrared observations reported here indicate the presence of three YSOs within the cloud. BIMA 3 mm continuum observations show that the brightest of the YSOs is likely to be a Herbig Ae/Be star. High-angular-resolution molecular-line observations suggest that a wind from this star may be triggering collapse in the adjacent molecular cloud. The submillimeter source G79.3+0.3 P1 itself does not contain infrared sources and may represent an earlier stage of star formation.

We present 1.4-mm BIMA observations with subarcsecond resolution of the bright dust and molecular line emission from hot cores associated with a sample of four ultracompact HII regions: G9.62+0.19, G10.47-0.03, G29.96-0.02, G31.41+0.31. Density power laws can reproduce the observed continuum emission but break down on scales smaller than 2000 AU. A total of 38 transitions from 18 species are detected, with G10.47 and G31.41 showing the greatest number of lines. In particular, these sources display emission from two collisionally-excited transitions of methanol lying more than 950 K above the ground state. Outflows traced by H2S emission provide evidence for embedded exciting sources and the observed morphology of molecular lines is consistent with internal heating of the cores.

We report the detection of direct l-type transitions towards the proto-planetary nebula CRL 618 during a study of vibrationally excited carbon chains. The J=8,10,11,12,13,14 DJ=0 transitions of HCN in its first excited bending mode v2=1 were detected in absorption against the continuum of the central HII region making use of the Effelsberg 100m telescope and the Very Large Array. Additionally, the J=9 direct l-type transition was detected in emission presumably indicating a weak maser. All lines are blueshifted with respect to the systemic velocity of CRL 618 indicating that the lines originate from a hot, expanding circumstellar envelope. The HCN column density along the line of sight in front of the continuum is 2E18 cm^-2.

We used the Effelsberg 100m and IRAM 30m telescopes to observe vibrationally excited cyanoacetylene (HCCCN) in several rotational transitions toward the proto-planetary nebula CRL618. Lines from 9 different vibrationally excited states with energies ranging up to 1600 K above ground were detected. The lines show P Cygni profiles indicating that the HCCCN emission originates from an expanding and accelerating molecular envelope. The HCCCN rotational temperature varies with velocity, peaks at 520 K, 3 km/s blue-shifted from the systemic velocity and decreases with higher blueshift of the gas. The column density of the absorbing HCCCN is 3-6 x 1E17 cm^2. We modeled spectra based on spherical models of the expanding envelope which provide an excellent fit to the observations, and discuss the implications of the models. Additionally, lines from 13C substituted cyanoacetylene were observed. They can be used to constrain the 12C/13C ratio in this source to 10+-2.

We present a comparison of ClassII CH3OH and H2O masers at high spatial resolution in a sample of 29 massive star-forming regions. Absolute positions of both maser types are compared with mm dust continuum, cm continuum and mid-infrared sources. Spatial correlations of the different tracers and kinematic signatures are discussed.

Well established correlations of low-mass flows show a continuity in the high-mass regime, and mean accretion rates for L>1000 Lsun are dM/dt~10^(-4) Msun/yr. We also present a tight correlation between outflow masses and core masses, which implies that the accretion rate is to first order a linear function of the core mass. All those results support the idea that similar physical processes drive outflows of all masses.

We present a detailed 1.2 mm continuum and CS spectral line study of a large sample of 69 massive star forming regions in very early stages of evolution, most of them prior to building up an ultracompact HII region. The continuum data show a zoo of different morphologies and give detailed information on the spatial distributions, the masses, column densities and average densities of the whole sample. Fitting the radial intensity profiles shows that three parameters are needed to describe the spatial distribution of the sources: constant emission from the center out to a few arcsec radius followed by a first power law intensity distribution which steepens further outside into a second power law distribution. The mean inner power law intensity index mi (I~r^(-mi)) is 1.2 corresponding to density indices p (n~r^(-p)) of 1.6. In total the density distribution of our massive star formations sites seem to be not too different from their low-mass counterparts, but we show that setting tight constrains on the density indices is very difficult and subject to many possible errors. The local densities we derive from CS calculations are higher (up to one order of magnitude) than the mean densities we find via the mm-continuum. Such inhomogeneous density distribution reflects most likely the ubiquitous phenomenon of clumping and fragmentation in molecular clouds. Linewidth-mass relations show a departure from virial equilibrium in the stages of strongly collapsing cores.

We describe a systematic program aimed at identifying and characterizing candidate high-mass proto-stellar objects (HMPOs). Our candidate sample consists of 69 objects selected by criteria based on those established by Ramesh & Sridharan (1997) using far-infrared, radio-continuum and molecular line data. Infrared-Astronomical-Satellite (IRAS) and Midcourse-Space-Experiment (MSX) data were used to study the larger scale environments of the candidate sources and to determine their total luminosities and dust temperatures. To derive the physical and chemical properties of our target regions, we observed continuum and spectral line radiation at millimeter and radio wavelengths. We imaged the free-free and dust continuum emission at wavelengths of 3.6 cm and 1.2 mm, respectively, searched for H2O and CH3OH maser emission and observed the CO 2-1 and several NH3 lines toward all sources in our sample. Other molecular tracers were observed in a subsample. The presented results indicate that a substantial fraction of our sample harbors HMPOs in a pre-UCHII region phase, the earliest known stage in the high-mass star formation process.

The VLA and the BIMA array were used to obtain high resolution (10-20") observations of C+, traced by the C91a recombination line at 8.6 GHz, and the dense molecular gas, traced by HCN and HCO+(1-0), of the photon dominated region (PDR) associated with the reflection nebula NGC 2023. Using the VLA, continuum emission is detected at 8.6 GHz from a faint HII region associated with HD 37903. The C91a emission originates from a 0.4 pc long filament, extending from the east to the south of the exciting star HD 37903. Within the filament three C91a clumps can be distinguished, each associated with filamentary vibrationally excited H2 emission in the direction toward HD 37903. The HCO+ emission has a clumpy appearance superimposed on a more extended component. C91a is, in general, closer to the exciting star than HCO+ emission as expected from PDR models. The morphologies of HCO+ and HCN are quite similar. Based on the C91a linewidth towards one of the clumps a limit of 170 K on the kinetic temperature in the ionized carbon layer can be derived. This value is consistent with PDR models with H2 densities of about 10^5/cm^3. However, this result suggests surprisingly low limits on the turbulence in the PDR. We detected a compact 3 mm continuum source in the PDR, which appears to be a cold "core" of density 10^7/cm^3, 0.03 parsec diameter, and 6 solar masses. We conclude that it may have formed within the PDR. In an appendix, observations of the C91a recombination line toward five additional PDRs using the Effelsberg 100m telescope are described.

We used the IRAM Interferometer to obtain sub-arcsecond resolution observations of the high-mass star-forming region W3(OH) and its surroundings at a frequency of 220 GHz. With the improved angular resolution, we distinguish 3 peaks in the thermal dust continuum emission originating from the hot core region about 6 arcsec (0.06 pc) east of W3(OH). The dust emission peaks are coincident with known radio continuum sources, one of which is of non-thermal nature. The latter source is also at the center of expansion of a powerful bipolar outflow observed in water maser emission. We determine the hot core mass to be 15 solar masses based on the integrated dust continuum emission. Simultaneously many molecular lines are detected allowing the analysis of the temperature structure and the distribution of complex organic molecules in the hot core. From HNCO lines, spanning a wide range of excitation, two 200 K temperature peaks are found coincident with dust continuum emission peaks suggesting embedded heating sources within them.

We have carried out VLA D-array observations of the C91alpha carbon recombination line as well as Effelsberg 100-m observations of the C65alpha line in a 5 arcmin square region centered between the Bar and the Trapezium stars in the Orion Nebula with spatial resolutions of 10 arcsec and 40 arcsec, respectively. The results show the ionized carbon in the PDR associated with the Orion Bar to be in a thin, clumpy layer sandwiched between the ionization front and the molecular gas. From the observed line widths we get an upper limit on the temperature in the C+ layer of 1500 K and from the line intensity a hydrogen density between 5 10^4 and 2.5 10^5 cm-3 for a homogeneous medium. The observed carbon level population is not consistent with predictions of hydrogenic recombination theory but could be explained by dielectronic recombination. The layer of ionized carbon seen in C91alpha is found to be essentially coincident with emission in the v=1-0 S(1) line of vibrationally excited molecular hydrogen. This is surprising in the light of current PDR models and some possible explanations of the discrepancy are discussed.

Interferometric observations of the CO(1-0) line in the elliptical galaxy NGC759 shows that it contains 2 10^9 Mo of molecular gas, confined to a circumnuclear ring with an average radius of only 650pc. The maximum gas surface density is 750 Mo/pc^2, which is comparable to the critical gas surface density for large scale gravitational instabilities. The low CO(2-1)/CO(1-0) ratio of 0.4 is consistent with a two-component molecular gas, consisting of a cold and dense phase, containing most of the mass, and a warm and diffuse component dominating the CO emission. We use a mass model for the underlying galaxy which is applicable to spherical galaxies with an r^1/4 luminosity profile. This leads to more modest estimates of the molecular gas fraction and gas surface density then would have been derived using simpler models, suggesting that many of the spectacular molecular gas properties of ultraluminous FIR galaxies, which could be described by similar mass models, may have to be revised. We discuss the possibility that NGC759 is in a late stage of a merging between two gas-rich disk galaxies.