The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.

Knowledge about the distribution of CO emission in the Milky Way is essential to understand the impact of Galactic environment on the formation and evolution of structures in the interstellar medium. However, currently our insight about the fraction of CO in spiral arm and interarm regions is still limited by large uncertainties in assumed rotation curve models or distance determination techniques. In this work we use a Bayesian approach to obtain the current best assessment of the distribution of 13CO from the Galactic Ring Survey. We performed two different distance estimates that either included or excluded a model for Galactic features. We also include a prior for the solution of the kinematic distance ambiguity that was determined from a compilation of literature distances and an assumed size-linewidth relationship. We find that the fraction of 13CO emission associated with spiral arm features varies from 76% to 84% between the two distance runs. The vertical distribution of the gas is concentrated around the Galactic midplane showing FWHM values of ~75 pc. We do not find any significant difference between gas emission properties associated with spiral arm and interarm features. In particular the distribution of velocity dispersion values of gas emission in spurs and spiral arms is very similar. We detect a trend of higher velocity dispersion values with increasing heliocentric distance, which we attribute to beam averaging effects caused by differences in spatial resolution. We argue that the true distribution of the gas emission is likely more similar to a combination of our two distance results, and highlight the importance of using complementary distance estimations to safeguard against the pitfalls of any single approach. We conclude that the methodology presented in this work is a promising way to determine distances to gas emission features in Galactic plane surveys.

We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA band, with $>$87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C$^{34}$S, H$_2$CS and CH$_{3}$CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source ($<$2000 $\times$ 760 au) at 1.2 cm, associated with mm1, of which $>$96% is from ionized gas. The ATCA NH$_3$(1,1) and (2,2) emission traces a large-scale (r$\sim$0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest.

Complex organic molecules (COMs) are thought to form on icy dust grains in the earliest phase of star formation. The evolution of these COMs from the youngest Class 0/I protostellar phases toward the more evolved Class II phase is still not fully understood. Since planet formation seems to start early, and mature disks are too cold for characteristic COM emission lines, studying the inventory of COMs on solar system scales in the Class 0/I stage is relevant. ALMA Band 3 (3 mm) and Band 6 (1 mm) observations are obtained of seven Class 0 protostars in the Perseus and Serpens star-forming regions. By modeling the inner protostellar region using 'LTE' models, the excitation temperature and column densities are determined for several O-bearing COMs. B1-c, B1-bS, and Serpens S68N show COM emission, i.e, three out of the seven sources. No clear correlation seems to exist between the occurrence of COMs and source luminosity. The abundances of several COMs with respect to CH3OH are remarkably similar for the three COM-rich sources, and to IRAS 16293-2422B and HH 212. For other COMs the abundances differ by up to an order of magnitude, indicating that local source conditions are case determining. B1-c hosts a cold ($T_{ex}\approx60$ K), more extended component of COM emission with a column density of typically a few % of the warm/hot ($T_{ex}\sim 200$ K), central component. A D/H ratio of 1-3 % is derived based on the CH2DOH/CH3OH ratio suggesting a temperature of $\sim$15~K during the formation of methanol. This ratio is consistent with other low-mass protostars. Future mid-infrared facilities such as JWST/MIRI will be essential to directly observe COM ices. Combining this with a larger sample of COM-rich sources with ALMA will allow for directly linking ice and gas-phase abundances in order to constrain the routes that produce and maintain chemical complexity during the star formation process.

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

Context. The [CII] 158 {\mu}m far-infrared (FIR) fine-structure line is one of the most important cooling lines of the star-forming interstellar medium (ISM). High spectral resolution observations have shown complex structures in the line profiles of the [CII] emission. Aims. Our aim is to determine whether the complex profiles observed in [^{12}CII] are due to individual velocity components along the line-of-sight or to self-absorption based on a comparison of the [^{12}CII] and isotopic [^{13}CII] line profiles. Methods. Deep integrations with the SOFIA/upGREAT 7-pixel array receiver in M43, Horsehead~PDR, Monoceros~R2, and M17~SW allow for the detection of optically thin [^{13}CII] emission lines, along with the [^{12}CII] emission lines, with a high signal-to-noise ratio. We first derived the [^{12}CII] optical depth and the [CII] column density from a single component model. However, the complex line profiles observed require a double layer model with an emitting background and an absorbing foreground. A multi-component velocity fit allows us to derive the physical conditions of the [CII] gas: column density and excitation temperature. Results. We find moderate to high [^{12}CII] optical depths in all four sources and self-absorption of [^{12}CII] in Mon R2 and M17 SW. The high column density of the warm background emission corresponds to an equivalent Av of up to 41 mag. The foreground absorption requires substantial column densities of cold and dense [CII] gas, with an equivalent Av ranging up to about 13 mag. Conclusions. The column density of the warm background material requires multiple photon-dominated region (PDR) surfaces stacked along the line of sight and in velocity. The substantial column density of dense and cold foreground [CII] gas detected in absorption cannot be explained with any known scenario and we can only speculate on its origins

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

The continuum emission from 1 to 2 GHz of The HI/OH/Recombination line survey of the inner Milky Way (THOR) at $\lesssim$18" resolution covers $\sim 132$ square degrees of the Galactic plane and detects 10387 sources. Similarly, the first data release of the Global View of Star Formation in the Milky Way (GLOSTAR) surveys covers $\sim 16$ square degrees of the Galactic plane from 4-8 GHz at 18" resolution and detects 1575 sources. However, a large fraction of the unresolved discrete sources detected in these radio continuum surveys of the Galactic plane remain unclassified. Here, we study the Euclidean-normalized differential source counts of unclassified and unresolved sources detected in these surveys and compare them with simulated extragalactic radio source populations as well as previously established source counts. We find that the differential source counts for THOR and GLOSTAR surveys are in excellent agreement with both simulation and previous observations. We also estimate the angular two-point correlation function of unclassified and unresolved sources detected in THOR survey. We find a higher clustering amplitude in comparison with the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey up to the angular separation of $5^{\circ}$. The decrease in angular correlation with increasing flux cut and the excellent agreement of clustering pattern of sources above 1 mJy with high $z$ samples ($z >0.5$) of the FIRST survey indicates that these sources might be high $z$ extragalactic compact objects. The similar pattern of one-point and two-point statistics of unclassified and compact sources with extragalactic surveys and simulations confirms the extragalactic origin of these sources.

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$n\alpha$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $\tau$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys.

The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information contained in the shape and Doppler-shifted frequency of the emission lines contained in these enormous data sets. In this work we present an overview and first results of a Gaussian decomposition of the entire Galactic Ring Survey (GRS) 13CO (1-0) data that consists of about 2.3 million spectra. We performed the decomposition with the fully automated GaussPy+ algorithm and fitted about 4.6 million Gaussian components to the GRS spectra. We discuss the statistics of the fit components and relations between the fitted intensities, velocity centroids, and velocity dispersions. We find that the magnitude of the velocity dispersion values increase toward the inner Galaxy and around the Galactic midplane, which we speculate is partly due to the influence of the Galactic bar and regions with higher non-thermal motions located in the midplane, respectively. We also use our decomposition results to infer global properties of the gas emission and find that the number of fit components used per spectrum is indicative for the amount of structure along the line of sight. We find that the emission lines from regions located on the far side of the Galaxy show increased velocity dispersion values, likely due to beam averaging effects. We demonstrate how this trend has the potential to aid in characterising Galactic structure by disentangling emission that is belonging to the nearby Aquila Rift molecular cloud from emission that is more likely associated with the Perseus and Outer spiral arms.

In the disk-mediated accretion scenario for the formation of the most massive stars, gravitational instabilities in the disk can force it to fragment. We investigate the effects of inclination and spatial resolution on observable kinematics and stability of disks in high-mass star formation. We study a high-resolution 3D radiation-hydrodynamic simulation that leads to the fragmentation of a massive disk. Using RADMC-3D we produce 1.3 mm continuum and CH3CN line cubes at different inclinations. The model is set to different distances and synthetic observations are created for ALMA at ~80 mas resolution and NOEMA at ~0.3''. The synthetic ALMA observations resolve all fragments and their kinematics well. The synthetic NOEMA observations at 800 pc (~300 au resolution) are able to resolve the fragments, while at 2000 pc (~800 au resolution) only a single slightly elongated structure is observed. The position-velocity (PV) plots show the differential rotation of material best in the edge-on views. As the observations become less resolved, the inner high-velocity components of the disk become blended with the envelope and the PV plots resemble rigid-body-like rotation. Protostellar mass estimates from PV plots of poorly resolved observations are therefore overestimated. We fit the emission of CH3CN lines and produce maps of gas temperature with values in the range of 100-300 K. Studying the Toomre stability of the disks in the resolved observations, we find Q values below the critical value for stability against gravitational collapse at the positions of the fragments and the arms connecting the fragments. For the poorly resolved observations we find low Q values in the outskirts of the disk. Therefore we are able to predict that the disk is unstable and fragmenting even in poorly resolved observations. This conclusion is true regardless of knowledge about the inclination of the disk.

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.

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario.

The chemical structure of high-mass star nurseries is important for a general understanding of star formation. Deuteration is a key chemical process in the earliest stages of star formation because its efficiency is sensitive to the environment. Using the IRAM-30 m telescope at 1.3--4.3 mm wavelengths, we have imaged two parsec-scale high-mass protostellar clumps (P1 and S) that show different evolutionary stages but are located in the same giant filamentary {infrared dark cloud} G28.34+0.06. Deep spectral images at subparsec resolution reveal the dust and gas physical structures of both clumps. We find that (1) the low-$J$ lines of $\rm N_2H^+$, HCN, HNC, and $\rm HCO^+$ isotopologues are subthermally excited; and (2) the deuteration of $\rm N_2H^+$ is more efficient than that of $\rm HCO^+$, HCN, and HNC by an order of magnitude. The deuterations of these species are enriched toward the chemically younger clump S compared with P1, indicating that this process favors the colder and denser environment ($\rm T_{kin} \sim14 K$, $\rm N(NH_3) \sim 9\times 10^{15}\,cm^{-2}$). In contrast, single deuteration of $\rm NH_3$ is insensitive to the environmental difference between P1 and S; and (3) single deuteration of $\rm CH_3OH$ ($\rm > 10\%$) is detected toward the location where CO shows a depletion of $\sim10$. This comparative chemical study between P1 and S links the chemical variations to the environmental differences and shows chemical similarities between the early phases of high- and low-mass star-forming regions.

(Abridged) The initial physical conditions of high-mass stars and protoclusters remain poorly characterized. To this end we present the first targeted ALMA 1.3mm continuum and spectral line survey towards high-mass starless clump candidates, selecting a sample of 12 of the most massive candidates ($400-4000\, M_\odot$) within 5 kpc. The joint 12+7m array maps have a high spatial resolution of $\sim 3000\, \mathrm{au}$ ($\sim 0.8^{\prime\prime}$) and have point source mass-completeness down to $\sim 0.3\, M_\odot$ at $6\sigma$ (or $1\sigma$ column density sensitivity of $1.1\times10^{22}\, \mathrm{cm^{-2}}$). We discover previously undetected signposts of low-luminosity star formation from CO (2-1) and SiO (5-4) bipolar outflows and other signatures towards 11 out of 12 clumps, showing that current MIR/FIR Galactic Plane surveys are incomplete to low- and intermediate-mass protostars ($\lesssim 50\, L_\odot$). We compare a subset of the observed cores with a suite of radiative transfer models of starless cores. We find a high-mass starless core candidate with a model-derived mass consistent with $29^{52}_{15}\, M_\odot$ when integrated over size scales of $2\times10^4\, \mathrm{au}$. Unresolved cores are poorly fit by starless core models, supporting the interpretation that they are protostellar even without detection of outflows. Substantial fragmentation is observed towards 10 out of 12 clumps. We extract sources from the maps using a dendrogram to study the characteristic fragmentation length scale. Nearest neighbor separations when corrected for projection are consistent with being equal to the clump average thermal Jeans length. Our findings support a hierarchical fragmentation process, where the highest density regions are not strongly supported against thermal gravitational fragmentation by turbulence or magnetic fields.

Context: OH masers trace diverse physical processes, from the expanding envelopes around evolved stars to star-forming regions or supernovae remnants. Aims: We identify the ground-state OH masers at 18cm wavelength in the area covered by ``The HI/OH/Recombination line survey of the Milky Way (THOR)''. We present a catalogue of all OH maser features and their possible associated environments. Methods: The THOR survey covers longitude and latitude ranges of 14.3<l<66.8 and b<1.25 deg. All OH ground state lines at 1612, 1665, 1667 and 1720MHz have been observed, employing the Very Large Array. The spatial resolution of the data varies between 12.5'' and 19'', the spectral resolution is 1.5km/s, and the rms sensitivity of the data is ~10mJy/beam per channel. Results: We identify 1585 individual maser spots distributed over 807 maser sites. Based on different criteria from spectral profiles to literature comparison, we try to associate the maser sites with astrophysical source types. Approximately 51\% of the sites exhibit the double-horned 1612MHz spectra typically emitted from the expanding shells of evolved stars. The separations of the two main velocity features of the expanding shells typically vary between 22 and 38km/s. In addition to this, at least 20% of the maser sites are associated with star-forming regions. While the largest fraction of 1720MHz maser spots (21 out of 53) is associated with supernova remnants, a significant fraction of the 1720MHz maser spots (17) are also associated with star-forming regions. We present comparisons to the thermal 13CO(1-0) emission as well as to other surveys of class II CH3OH and H2O maser emission. The catalogue attempts to present associations to astrophysical sources where available, and the full catalogue is available in electronic form.

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

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

Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

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

Context. The role of magnetic fields during the formation of high-mass stars is not yet fully understood, and the processes related to the early fragmentation and collapse are largely unexplored today. The high-mass star forming region G9.62+0.19 is a well known source, presenting several cores at different evolutionary stages. Aims. We determine the magnetic field morphology and strength in the high-mass star forming region G9.62+0.19, to investigate its relation to the evolutionary sequence of the cores. Methods. We use Band 7 ALMA observations in full polarisation mode and we analyse the polarised dust emission. We estimate the magnetic field strength via the Davis-Chandrasekhar-Fermi and the Structure Function methods. Results. We resolve several protostellar cores embedded in a bright and dusty filamentary structure. The polarised emission is clearly detected in six regions. Moreover the magnetic field is oriented along the filament and appears perpendicular to the direction of the outflows. We suggest an evolutionary sequence of the magnetic field, and the less evolved hot core exhibits a magnetic field stronger than the more evolved one. We detect linear polarisation from thermal line emission and we tentatively compared linear polarisation vectors from our observations with previous linearly polarised OH masers observations. We also compute the spectral index, the column density and the mass for some of the cores. Conclusions. The high magnetic field strength and the smooth polarised emission indicate that the magnetic field could play an important role for the fragmentation and the collapse process in the star forming region G9.62+019 and that the evolution of the cores can be magnetically regulated. On average, the magnetic field derived by the linear polarised emission from dust, thermal lines and masers is pointing in the same direction and has consistent strength.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs.

The past two decades have seen extensive surveys of the far-infrared to submillimeter continuum emission in the plane of our Galaxy. We line out prospects for the coming decade for corresponding molecular and atomic line surveys which are needed to fully understand the formation of the dense structures that give birth to clusters and stars out of the diffuse interstellar medium. We propose to work towards Galaxy wide surveys in mid-J CO lines to trace shocks from colliding clouds, Galaxy-wide surveys for atomic Carbon lines in order to get a detailed understanding of the relation of atomic and molecular gas in clouds, and to perform extensive surveys of the structure of the dense parts of molecular clouds to understand the importance of filaments/fibers over the full range of Galactic environments and to study how dense cloud cores are formed from the filaments. This work will require a large (50m) Single Dish submillimeter telescope equipped with massively multipixel spectrometer arrays, such as envisaged by the AtLAST project.

In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $\sigma_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $\sigma_{v_{los}, {\rm 1D}}/\sigma_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ.

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

Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separations vary between sub-100au to several 100 and 1000au. The core separations and approximate masses are largely consistent with thermal Jeans fragmentation of a dense gas core. Due to the high continuum optical depth, most spectral lines are seen in absorption. However, a few exceptional emission lines are found that most likely stem from transitions with excitation conditions above1000K. Toward the main continuum source, these emission lines exhibit a velocity gradient across scales of 100-200au aligned with the molecular outflow and perpendicular to the previously inferred disk orientation. While we cannot exclude that these observational features stem from an inner hot accretion disk, the alignment with the outflow rather suggests that it stems from the inner jet and outflow region. The highest-velocity features are found toward the peak position, and no Hubble-like velocity structure can be identified. Therefore, these data are consistent with steady-state turbulent entrainment of the hot molecular gas via Kelvin-Helmholtz instabilities at the interface between the jet and the outflow. Conclusions: Resolving this high-mass star-forming region at sub-50au scales indicates that the hierarchical fragmentation process in the framework of thermal Jeans fragmentation can continue down to the smallest accessible spatial scales. Velocity gradients on these small scales have to be treated cautiously and do not necessarily stem from disks, but may be better explained with outflow emission.

We have carried out the largest and most unbiased search for hypercompact (HC) HII regions. Our method combines four interferometric radio continuum surveys (THOR, CORNISH, MAGPIS and White2005) with far-infrared and sub-mm Galactic Plane surveys to identify embedded HII regions with positive spectral indices. 120 positive spectrum HII regions have been identified from a total sample of 534 positive spectral index radio sources. None of these HII regions, including the known HCHII regions recovered in our search, fulfills the canonical definition of an HCHII region at 5 GHz. We suggest that the current canonical definition of HCHII regions is not accurate and should be revised to include a hierarchical structure of ionized gas that results in an extended morphology at 5 GHz. Correlating our search with known ultracompact (UC) HII region surveys, we find that roughly half of detected UCHII regions have positive spectral indices, instead of more commonly assumed flat and optically thin spectra. This implies a mix of optically thin and thick emission and has important implications for previous analyses which have so far assumed optically thin emission for these objects. Positive spectrum HII regions are statistically more luminous and possess higher Lyman continuum fluxes than HII regions with flat or negative indices. Positive spectrum HII regions are thus more likely to be associated with more luminous and massive stars. No differences are found in clump mass, linear diameter or luminosity-to-mass ratio between positive spectrum and non-positive spectrum HII regions.

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

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.

Radio continuum surveys of the Galactic plane can find and characterize HII regions, supernova remnants (SNRs), planetary nebulae (PNe), and extragalactic sources. A number of surveys at high angular resolution (<25") at different wavelengths exist to study the interstellar medium (ISM), but no comparable high-resolution and high-sensitivity survey exists at long radio wavelengths around 21cm. We observed a large fraction of the Galactic plane in the first quadrant of the Milky Way (l=14.0-67.4deg and |b| < 1.25deg) with the Karl G. Jansky Very Large Array (VLA) in the C-configuration covering six continuum spectral windows. These data provide a detailed view on the compact as well as extended radio emission of our Galaxy and thousands of extragalactic background sources. We used the BLOBCAT software and extracted 10916 sources. After removing spurious source detections caused by the sidelobes of the synthesised beam, we classified 10387 sources as reliable detections. We smoothed the images to a common resolution of 25" and extracted the peak flux density of each source in each spectral window (SPW) to determine the spectral indices $\alpha$ (assuming $I(\nu)\propto\nu^\alpha$). By cross-matching with catalogs of HII regions, SNRs, PNe, and pulsars, we found radio counterparts for 840 HII regions, 52 SNRs, 164 PNe, and 38 pulsars. We found 79 continuum sources that are associated with X-ray sources. We identified 699 ultra-steep spectral sources ($\alpha < -1.3$) that could be high-redshift galaxies. Around 9000 of the sources we extracted are not classified specifically, but based on their spatial and spectral distribution, a large fraction of them is likely to be extragalactic background sources. More than 7750 sources do not have counterparts in the SIMBAD database, and more than 3760 sources do not have counterparts in the NED database.

Anderson et al. (2017) identified seventy six candidate supernova remnants (SNRs) using data from The HI, OH, Recombination line survey of the Milky Way (THOR). The spectral index and polarization properties can help distinguish between SNRs and H II regions, which are often confused. We confirm two SNR candidates using spectral index data and morphology. However, we observe that the fractional linear polarization cannot distinguish between SNRs and H II regions, likely due to contamination by diffuse Galactic synchrotron emission. We also comment on the association of SNR candidates with pulsars through geometric and age considerations.

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

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

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.

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.

An increasing number of hundred-parsec scale, high line-mass filaments have been detected in the Galaxy. Their evolutionary path, including fragmentation towards star formation, is virtually unknown. We characterize the fragmentation within the Nessie filament, covering size-scales between $\sim$ 0.1-100 pc. We also connect the small-scale fragments to the star-forming potential of the cloud. We combine near-infrared data from the VVV survey with mid-infrared GLIMPSE data to derive a high-resolution dust extinction map and apply a wavelet decomposition technique on it to analyze the fragmentation characteristics of the cloud, which are compared with predictions from fragmentation models. We compare the detected objects to those identified in $\sim$ 10 times coarser resolution from ATLASGAL data. We present a high-resolution extinction map of Nessie. We estimate the mean line-mass of Nessie to be $\sim$ 627 M$_\odot$/pc and the distance to be $\sim$ 3.5 kpc. We find that Nessie shows fragmentation at multiple size scales. The nearest-neighbour separations of the fragments at all scales are within a factor of 2 of the Jeans' length at that scale. However, the relationship between the mean densities of the fragments and their separations is significantly shallower than expected for Jeans' fragmentation. The relationship is similar to the one predicted for a filament that exhibits a Larson-like scaling between size-scale and velocity dispersion; such a scaling may result from turbulent support. Based on the number of YSOs in Nessie, we estimate that the star formation rate is $\sim$ 371 M$_\odot$/Myr; similar values result if using the number of dense cores, or the amount of dense gas, as the proxy of star formation. The star formation efficiency is 0.017. These numbers indicate that Nessie's star-forming content is comparable to the Solar neighborhood giant molecular clouds like Orion A.

Context: The importance of magnetic fields at the onset of star formation related to the early fragmentation and collapse processes is largely unexplored today. Aims: We want to understand the magnetic field properties at the earliest evolutionary stages of high-mass star formation. Methods: The Atacama Large Millimeter Array is used at 1.3mm wavelength in full polarization mode to study the polarized emission and by that the magnetic field morphologies and strengths of the high-mass starless region IRDC18310-4. Results: The polarized emission is clearly detected in four sub-cores of the region. In general it shows a smooth distribution, also along elongated cores. Estimating the magnetic field strength via the Davis-Chandrasekhar-Fermi method and following a structure function analysis, we find comparably large magnetic field strengths between ~0.6 and 3.7mG. Comparing the data to spectral line observations, the turbulent-to-magnetic energy ratio is low, indicating that turbulence does not significantly contribute to the stability of the gas clump. A mass-to-flux ratio around the critical value 1.0 - depending on column density - indicates that the region starts to collapse which is consistent with the previous spectral line analysis of the region. Conclusions: While this high-mass region is collapsing and thus at the verge of star formation, the high magnetic field values and the smooth spatial structure indicate that the magnetic field is important for the fragmentation and collapse process. This single case study can only be the starting point for larger sample studies of magnetic fields at the onset of star formation.

We report ALMA observations with resolution $\approx0.5$" at 3 mm of the extended Sgr B2 cloud in the Central Molecular Zone (CMZ). We detect 271 compact sources, most of which are smaller than 5000 AU. By ruling out alternative possibilities, we conclude that these sources consist of a mix of hypercompact HII regions and young stellar objects (YSOs). Most of the newly-detected sources are YSOs with gas envelopes which, based on their luminosities, must contain objects with stellar masses $M_*\gtrsim8$ M$_\odot$. Their spatial distribution spread over a $\sim12\times3$ pc region demonstrates that Sgr B2 is experiencing an extended star formation event, not just an isolated `starburst' within the protocluster regions. Using this new sample, we examine star formation thresholds and surface density relations in Sgr B2. While all of the YSOs reside in regions of high column density ($N(H_2)\gtrsim2\times10^{23}$ cm$^{-2}$), not all regions of high column density contain YSOs. The observed column density threshold for star formation is substantially higher than that in solar vicinity clouds, implying either that high-mass star formation requires a higher column density or that any star formation threshold in the CMZ must be higher than in nearby clouds. The relation between the surface density of gas and stars is incompatible with extrapolations from local clouds, and instead stellar densities in Sgr B2 follow a linear $\Sigma_*-\Sigma_{gas}$ relation, shallower than that observed in local clouds. Together, these points suggest that a higher volume density threshold is required to explain star formation in CMZ clouds.

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

The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.

We present Submillimeter Array (SMA) observations toward the high-mass star-forming region IRAS 18566+0408. Observations at 1.3 mm continuum and in several molecular line transitions were performed in the compact (2."4 angular resolution) and very-extended (~0."4 angular resolution) configurations. The continuum emission from the compact configuration shows a dust core of 150 Msun, while the very-extended configuration reveals a dense (2.6 x 10^7 cm^-3) and compact (~4,000 AU) condensation of 8 Msun. We detect 31 molecular transitions from 14 species including CO isotopologues, SO, CH3OH, OCS, and CH3CN. Using the different k-ladders of the CH3CN line, we derive a rotational temperature at the location of the continuum peak of 240 K. The 12CO(2-1), 13CO(2-1), and SO(6_5-5_4) lines reveal a molecular outflow at PA ~135^o centered at the continuum peak. The extended 12CO(2-1) emission has been recovered with the IRAM 30 m telescope observations. Using the combined data set, we derive an outflow mass of 16.8 Msun. The chemically rich spectrum and the high rotational temperature confirm that IRAS 18566+0408 is harboring a hot molecular core. We find no clear velocity gradient that could suggest the presence of a rotational disk-like structure, even at the high resolution observations obtained with the very-extended configuration.

Context. The importance of the magnetic field in high-mass-star formation is not yet fully clear and there are still many open questions concerning its role in the accretion processes and generation of jets and outflows. In the past few years, masers have been successfully used to probe the magnetic field morphology and strength at scales of a few au around massive protostars, by measuring linear polarisation angles and Zeeman splitting. The massive protostar IRAS 18089-1732 is a well studied high-mass-star forming region, showing a hot core chemistry and a disc-outflow system. Previous SMA observations of polarised dust revealed an ordered magnetic field oriented around the disc of IRAS 18089-1732. Aims. We want to determine the magnetic field in the dense region probed by 6.7 GHz methanol maser observations and compare it with observations in dust continuum polarisation, to investigate how the magnetic field in the compact maser region relates to the large-scale field around massive protostars. Methods. We reduced MERLIN observations at 6.7 GHz of IRAS 18089-1732 and we analysed the polarised emission by methanol masers. Results. Our MERLIN observations show that the magnetic field in the 6.7 GHz methanol maser region is consistent with the magnetic field constrained by the SMA dust polarisation observations. A tentative detection of circularly polarised line emission is also presented. Conclusions. We found that the magnetic field in the maser region has the same orientation as in the disk. Thus the large-scale field component, even at the au scale of the masers, dominates over any small-scale field fluctuations. We obtained, from the circular polarisation tentative detection, a field strength along the line of sight of 5.5 mG which appeared to be consistent with the previous estimates.

There is a considerable deficiency in the number of known supernova remnants (SNRs) in the Galaxy compared to that expected. Searches for extended low-surface brightness radio sources may find new Galactic SNRs, but confusion with the much larger population of HII regions makes identifying such features challenging. SNRs can, however, be separated from HII regions using their significantly lower mid-infrared (MIR) to radio continuum intensity ratios. We use the combination of high-resolution 1-2 GHz continuum data from The HI, OH, Recombination line survey of the Milky Way (THOR) and lower-resolution VLA 1.4 GHz Galactic Plane Survey (VGPS) continuum data, together with MIR data from the Spitzer GLIMPSE, Spitzer MIPSGAL, and WISE surveys to identify SNR candidates. To ensure that the candidates are not being confused with HII regions, we exclude radio continuum sources from the WISE Catalog of Galactic HII Regions, which contains all known and candidate H II regions in the Galaxy. We locate 76 new Galactic SNR candidates in the THOR and VGPS combined survey area of 67.4deg>l>17.5deg, |b|<1.25deg and measure the radio flux density for 52 previously-known SNRs. The candidate SNRs have a similar spatial distribution to the known SNRs, although we note a large number of new candidates near l=30deg, the tangent point of the Scutum spiral arm. The candidates are on average smaller in angle compared to the known regions, 6.4'+/-4.7' versus 11.0'+/-7.8', and have lower integrated flux densities. If the 76 candidates are confirmed as true SNRs, for example using radio polarization measurements or by deriving radio spectral indices, this would more than double the number of known Galactic SNRs in the survey area. This large increase would still, however, leave a discrepancy between the known and expected SNR populations of about a factor of two.

OTS44 is one of only four free-floating planets known to have a disk. We have previously shown that it is the coolest and least massive known free-floating planet ($\sim$12 M$_{\rm Jup}$) with a substantial disk that is actively accreting. We have obtained Band 6 (233 GHz) ALMA continuum data of this very young disk-bearing object. The data shows a clear unresolved detection of the source. We obtained disk-mass estimates via empirical correlations derived for young, higher-mass, central (substellar) objects. The range of values obtained are between 0.07 and 0.63 M$_{\oplus}$ (dust masses). We compare the properties of this unique disk with those recently reported around higher-mass (brown dwarfs) young objects in order to infer constraints on its mechanism of formation. While extreme assumptions on dust temperature yield disk-mass values that could slightly diverge from the general trends found for more massive brown dwarfs, a range of sensible values provide disk masses compatible with a unique scaling relation between $M_{\rm dust}$ and $M_{*}$ through the substellar domain down to planetary masses.

Context: High-mass stars have a high degree of multiplicity and most likely form via disk accretion processes. The detailed physics of the binary and disk formation are still poorly constrained. Methods: Using the VLA in its most extended configuration at ~24GHz toward the prototypical high-mass star-forming region NGC7538IRS1 has allowed us to study the NH3 and thermal CH3OH emission and absorption as well as the cm continuum emission at an unprecedented spatial resolution of 0.06"x0.05", corresponding to a linear resolution of ~150AU at a distance of 2.7kpc. Results: A comparison of these new cm continuum data with previous VLA observations from 23yrs ago reveals no recognizable proper motions. If the emission were caused by a protostellar jet, proper motion signatures should have been easily identified. In combination with the high spectral indices S~nu^{alpha} (alpha between 1 and 2), this allows us to conclude that the continuum emission is from two hypercompact HII regions separated in projection by about 430AU. The NH3 spectral line data reveal a common rotating envelope indicating a bound high-mass binary system. In addition to this, the thermal CH3OH data show two separate velocity gradients across the two hypercompact HII regions. This indicates two disk-like structures within the same rotating circumbinary envelope. Disk and envelope structures are inclined by ~33deg, which can be explained by initially varying angular momentum distributions within the natal, turbulent cloud. Conclusions: Studying high-mass star formation at sub-0.1" resolution allows us to isolate multiple sources as well as to separate circumbinary from disk-like rotating structures.

Aims: We resolve the small-scale structure around the high-mass hot core region G351.77-0.54 to investigate its disk and fragmentation properties. Methods: Using ALMA at 690GHz with baselines exceeding 1.5km, we study the dense gas, dust and outflow emission at an unprecedented spatial resolution of 0.06" (130AU@2.2kpc). Results: Within the inner few 1000AU, G351.77 fragments into at least four cores (brightness temperatures between 58 and 197K). The central structure around the main submm source #1 with a diameter of ~0.5" does not show additional fragmentation. While the CO(6-5) line wing emission shows an outflow lobe in the north-western direction emanating from source #1, the dense gas tracer CH3CN shows a velocity gradient perpendicular to the outflow that is indicative of rotational motions. Absorption profile measurements against the submm source #2 indicate infall rates on the order of 10^{-4} to 10^{-3}M_sun/yr which can be considered as an upper limit of the mean accretion rates. The position-velocity diagrams are consistent with a central rotating disk-like structure embedded in an infalling envelope, but they may also be influenced by the outflow. Using the CH_3CN(37_k-36_k) k-ladder with excitation temperatures up to 1300K, we derive a gas temperature map of source #1 exhibiting temperatures often in excess of 1000K. Brightness temperatures of the submm continuum never exceed 200K. This discrepancy between gas temperatures and submm dust brightness temperatures (in the optically thick limit) indicates that the dust may trace the disk mid-plane whereas the gas could be tracing a hotter gaseous disk surface layer. In addition, we conduct a pixel-by-pixel Toomre gravitational stability analysis of the central rotating structure. The derived high Q values throughout the structure confirm that this central region appears stable against gravitational instability.

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.

We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass $\sim$ 400 M$_\odot$ pc$^{-1}$) in the Orion A molecular cloud. We have observed a 1.6 pc long section of the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of $\sim$3" (1 200 AU). We identify from the region 43 dense cores with masses about a solar mass. 60% of the ALMA cores are protostellar and 40\% are starless. The nearest neighbour separations of the cores do not show a preferred fragmentation scale; the frequency of short separations increases down to 1 200 AU. We apply a two-point correlation analysis on the dense core separations and show that the ALMA cores are significantly grouped at separations below $\sim$17 000 AU and strongly grouped below $\sim$6 000 AU. The protostellar and starless cores are grouped differently: only the starless cores group strongly below $\sim$6 000 AU. In addition, the spatial distribution of the cores indicates periodic grouping of the cores into groups of $\sim$30 000 AU in size, separated by $\sim$50 000 AU. The groups coincide with dust column density peaks detected by Herschel. These results show hierarchical, two-mode fragmentation in which the maternal filament periodically fragments into groups of dense cores. Critically, our results indicate that the fragmentation models for lower line-mass filaments ($\sim$ 16 M$_\odot$ pc$^{-1}$) fail to capture the observed properties of the ISF. We also find that the protostars identified with Spitzer and Herschel in the ISF are grouped at separations below $\sim$17 000 AU. In contrast, young stars with disks do not show significant grouping. This suggests that the grouping of dense cores is partially retained over the protostar lifetime, but not over the lifetime of stars with disks.

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 process of atomic-to-molecular (HI-to-H$_2$) gas conversion is fundamental for molecular-cloud formation and star formation. 21 cm observations of the star-forming region W43 revealed extremely high HI column densities, of 120-180 M$_{\odot}$ pc$^{-2}$, a factor of 10-20 larger than predicted by HI-to-H$_2$ transition theories. We analyze the observed HI with an HI-to-H$_2$ transition theoretical model, and show that the theory-observation discrepancy cannot be explained by the intense radiation in W43, nor by variations of the assumed volume density or H$_2$ formation-rate coefficient. We show that the large observed HI columns are naturally explained by several ($9-22$) HI-to-H$_2$ transition layers, superimposed along the sightlines of W43. We discuss other possible interpretations such as a non-steady-state scenario, and inefficient dust absorption. The case of W43 suggests that HI thresholds reported in extra-galactic observations are probably not associated with a single HI-to-H$_2$ transition, but are rather a result of several transition layers (clouds) along the sightlines, beam-diluted with diffuse inter-cloud gas.

We infer the absolute time dependence of kinematic gas temperature along a proposed orbit of molecular clouds in the Central Molecular Zone (CMZ) of the Galactic Center (GC). Ammonia gas temperature maps are one of the results of the "Survey of Water and Ammonia in the Galactic Center" (SWAG, PI: J. Ott); the dynamical model of molecular clouds in the CMZ was taken from Kruijssen et al. (2015). We find that gas temperatures increase as a function of time in both regimes before and after the cloud passes pericenter on its orbit in the GC potential. This is consistent with the recent proposal that pericenter passage triggers gravitational collapse. Other investigated quantities (line width, column density, opacity) show no strong sign of time dependence but are likely dominated by cloud-to-cloud variations.

Aims: To gain insight into the expected gas dynamics at the interface of the Galactic bar and spiral arms in our own Milky Way galaxy, we examine as an extragalactic counterpart the evidence for multiple distinct velocity components in the cold, dense molecular gas populating a comparable region at the end of the bar in the nearby galaxy NGC3627. Methods: We assemble a high resolution view of molecular gas kinematics traced by CO(2-1) emission and extract line-of-sight velocity profiles from regions of high and low gas velocity dispersion. Results: The high velocity dispersions arise with often double-peaked or multiple line-profiles. We compare the centroids of the different velocity components to expectations based on orbital dynamics in the presence of bar and spiral potential perturbations. A model of the region as the interface of two gas-populated orbits families supporting the bar and the independently rotating spiral arms provides an overall good match to the data. An extent of the bar to the corotation radius of the galaxy is favored. Conclusions: Using NGC3627 as an extragalactic example, we expect situations like this to favor strong star formation events such as observed in our own Milky Way since gas can pile up at the crossings between the orbit families. The relative motions of the material following these orbits is likely even more important for the build up of high density in the region. The surface densities in NGC3627 are also so high that shear at the bar end is unlikely to significantly weaken the star formation activity. We speculate that scenarios in which the bar and spiral rotate at two different pattern speeds may be the most favorable for intense star formation at such interfaces.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

The mass of prestellar cores is an essential ingredient to understand the onset of star formation in the core. The low level of emission from cold dust may keep parts of it hidden from observation. We aim to determine the fraction of core mass in the temperature range < 8 K that can be expected for typical low- and high-mass star formation regions. We calculate the dust temperature within standard spherically symmetric prestellar cores for a grid of density powerlaws in the outer core regions, core masses, and variation in the external multi-component radiation field. The dust is assumed to be composed of amorphous silicate and carbon, and variations of its optical properties are discussed. As measure for the distribution of cores and clumps, we use core mass functions derived for various environments. In view of the high densities in very cold central regions, dust and gas temperatures are assumed to be equal. We find that the fraction of mass with temperatures < 8 K in typical low- and high-mass cores is < 20\%. It is possible to obtain higher fractions of very cold gas by placing intermediate- or high-mass cores in a typical low-mass star formation environment. We show that the mass uncertainty arising from FIR to mm modeling of very cold dust emission is smaller than the mass uncertainty due to the unknown dust opacities.

The dense, cold regions where high-mass stars form are poorly characterised, yet they represent an ideal opportunity to learn more about the initial conditions of high-mass star formation (HMSF), since high-mass starless cores (HMSCs) lack the violent feedback seen at later evolutionary stages. We present continuum maps obtained from the Submillimeter Array (SMA) interferometry at 1.1 mm for four infrared dark clouds (IRDCs, G28.34S, IRDC 18530, IRDC 18306, and IRDC 18308). We also present 1 mm/3 mm line surveys using IRAM 30 m single-dish observations. Our results are: (1) At a spatial resolution of 10^4 AU, the 1.1 mm SMA observations resolve each source into several fragments. The mass of each fragment is on average >10 Msun, which exceeds the predicted thermal Jeans mass of the whole clump by a factor of up to 30, indicating that thermal pressure does not dominate the fragmentation process. Our measured velocity dispersions in the 30 m lines imply that non-thermal motions provides the extra support against gravity in the fragments. (2) Both non-detection of high-J transitions and the hyperfine multiplet fit of N2H+(1-0), C2H(1-0), HCN(1-0), and H13CN(1-0) indicate that our sources are cold and young. However, obvious detection of SiO and the asymmetric line profile of HCO+(1-0) in G28.34S indicate a potential protostellar object and probable infall motion. (3) With a large number of N-bearing species, the existence of carbon rings and molecular ions, and the anti-correlated spatial distributions between N2H+/NH2D and CO, our large-scale high-mass clumps exhibit similar chemical features as small-scale low-mass prestellar objects. This study of a small sample of IRDCs illustrates that thermal Jeans instability alone cannot explain the fragmentation of the clump into cold (~15 K), dense (>10^5 cm-3) cores and that these IRDCs are not completely quiescent.

Large-scale gaseous filaments with length up to the order of 100 pc are on the upper end of the filamentary hierarchy of the Galactic interstellar medium. Their association with respect to the Galactic structure and their role in Galactic star formation are of great interest from both observational and theoretical point of view. Previous "by-eye" searches, combined together, have started to uncover the Galactic distribution of large filaments, yet inherent bias and small sample size limit conclusive statistical results to be drawn. Here, we present (1) a new, automated method to identify large-scale velocity-coherent dense filaments, and (2) the first statistics and the Galactic distribution of these filaments. We use a customized minimum spanning tree algorithm to identify filaments by connecting voxels in the position-position-velocity space, using the Bolocam Galactic Plane Survey spectroscopic catalog. In the range of $7.^{\circ}5 \le l \le 194^{\circ}$, we have identified 54 large-scale filaments and derived mass ($\sim 10^3 - 10^5 \, M_\odot$), length (10-276 pc), linear mass density (54-8625 $M_\odot \, \rm{pc}^{-1}$), aspect ratio, linearity, velocity gradient, temperature, fragmentation, Galactic location and orientation angle. The filaments concentrate along major spiral arms. They are widely distributed across the Galactic disk, with 50% located within $\pm$20 pc from the Galactic mid-plane and 27% run in the center of spiral arms (aka "bones"). An order of 1% of the molecular ISM is confined in large filaments. Massive star formation is more favorable in large filaments compared to elsewhere. This is the first comprehensive catalog of large filaments useful for a quantitative comparison with spiral structures and numerical simulations.

To unambiguously diagnose the detailed physical mechanisms and the evolutionary status of high-mass star-forming regions (HMSFRs), we have performed 0.4" ($\sim1000$ AU) resolution observations towards NGC7538S and IRS1, using the Plateau de Bure Interferometer (PdBI). This paper presents a joint analysis of the 1.37 mm continuum emission and the line intensity of 15 molecular species. Assuming local thermal equilibrium, we derived molecular column densities and molecular abundances for each internal gas substructure which is spatially resolved. These derived quantities are compared with a suite of 1-D gas-grain models. NGC7538S is resolved into at least three dense gas condensations. Despite the comparable continuum intensity of these condensations, their differing molecular line emission is suggestive of an overall chemical evolutionary trend from the northeast to the southeast. Since these condensations are embedded within the same parent gas core, their differing chemical properties are the most likely due to the different warm-up histories, rather than the different dynamic timescales. Despite remaining spatially unresolved, in IRS1 we detect abundant complex organic molecules, indicating that IRS1 is the most chemically evolved hot molecular core presented here. We observe a continuum that is dominated by absorption features with at least three strong emission lines, potentially from $\rm CH_3OH$. The $\rm CH_3OH$ lines which are purely in emission have higher excitation than the ones showing purely absorption. Potential reasons for that difference are discussed. This is the first comprehensive comparison of observations of the high-mass cores NGC7538S and IRS1 and a chemical model. We have found that different chemical evolutionary stages can coexist in the same natal gas core. Our achievement illustrates the strength of chemical analysis for understanding HMSFRs.

We present observations towards a high-mass ($\rm >40\,M_{\odot}$), low luminosity ($\rm <10\,L_{\odot}$) $\rm 70\,\mu$m dark molecular core G 28.34 S-A at 3.4 mm, using the IRAM 30 m telescope and the NOEMA interferometer. We report the detection of $\rm SiO$ $J=\rm 2\rightarrow1$ line emission, which is spatially resolved in this source at a linear resolution of $\sim$0.1 pc, while the 3.4 mm continuum image does not resolve any internal sub-structures. The SiO emission exhibits two W-E oriented lobes centring on the continuum peak. Corresponding to the red-shifted and blue-shifted gas with velocities up to $\rm 40\,km\,s^{-1}$ relative to the quiescent cloud, these lobes clearly indicate the presence of a strong bipolar outflow from this $\rm 70\,\mu$m dark core, a source previously considered as one of the best candidates of "starless" core. Our SiO detection is consistent with ALMA archival data of $\rm SiO$ $J=\rm 5\rightarrow4$, whose high-velocity blue-shifted gas reveals a more compact lobe spatially closer to the dust center. This outflow indicates that the central source may be in an early evolutionary stage of forming a high-mass protostar. We also find that the low-velocity components (in the range of $\rm V_{lsr}$$\rm_{-5}^{+3}\,km\,s^{-1}$) have an extended, NW-SE oriented distribution. Discussing the possible accretion scenarios of the outflow-powering young stellar object, we argue that the molecular line emission and the molecular outflows may provide a better indication of the accretion history when forming young stellar object, than that from a snapshot observations of the present bolometric luminosity. This is particularly significant for the cases of episodic accretion, which may occur during the collapse of the parent molecular core.

Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search of GMFs in the fourth Galactic quadrant and determine their basic physical properties. We identify GMFs based on their dust extinction signatures in near- and mid-infrared and velocity structure probed by ^{13}CO line emission. We use the ^{13}CO line emission and ATLASGAL dust emission data to estimate the total and dense gas masses of the GMFs. We combine our sample with an earlier sample from literature and study the Galactic environment of the GMFs. We identify nine GMFs in the fourth Galactic quadrant; six are located in the Centaurus spiral arm and three in inter-arm regions. Combining this sample with an earlier study using the same identification criteria in the first Galactic quadrant results in 16 GMFs, nine of which are located within spiral arms. The GMFs have sizes of 80-160 pc and ^{13}CO-derived masses between 5-90 x 10^{4} Msun. Their dense gas mass fractions are between 1.5-37%, being higher in the GMFs connected to spiral arms. We also compare the different GMF-identification methods and find that emission and extinction based techniques overlap only partially, highlighting the need to use both to achieve a complete census.

Aims: Previous observations of the young massive star formation region IRAS 19410+2336 have revealed strong outflow activity with several interacting outflows. We aim to get a better understanding of the outflow activity in this region by observing the SiO and H$^{13}$CO$^+$ emission with high angular resolution. SiO is known to trace shocked gas, which is often associated with young energetic outflows. With the H$^{13}$CO$^+$ data, we intend to better understand the distribution of the quiescent gaseous component of the region. Methods: The SiO observations in the J=2-1 v=0 transition and H$^{13}$CO$^+$ J=1-0 observations were performed by the Plateau de Bure Interferometer, combined with IRAM 30 m single-dish observations, in order to get the missing short-spacing information. We complement this new high-resolution observation with earlier CO and H$_2$ data. Results: The SiO observations do not trace the previously in CO and H$_2$ identified outflows well. Although we identify regions of highly increased SiO abundance indicative of shock interaction, there are hardly any bipolar structures in the data. The southern part of the region, which exhibits strong H$_2$ emission, shows almost no SiO. The CO and SiO data show only weak similarities, and the main SiO emission lies between the two dominating dust clumps of the region. Conclusions: Most SiO emission is likely to be a result of high-velocity shocks due to protostellar jets. However, this does not explain all the emission features and additional effects; for example, colliding gas flows at the interface of the two main regions may play an important role in the origin of the emission. The present SiO data show that several different effects can influence SiO emission, which makes the interpretation of SiO data more difficult than often assumed.

The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a high-mass star-forming environment. For the first time, we complemented 1.3 mm Submillimeter Array (SMA) interferometric line survey with IRAM 30 m single-dish observations of the Orion-KL region. Covering a 4 GHz bandwidth in total, this survey contains over 160 emission lines from 20 species (25 isotopologues), including 11 complex organic molecules (COMs). At a spatial resolution of 1200 AU, the continuum substructures are resolved. Extracting the spectra from individual substructures and providing the intensity-integrated distribution map for each species, we studied the small-scale chemical variations in this region. Our main results are: (1) We identify lines from the low-abundance COMs CH3COCH3 and CH3CH2OH, as well as tentatively detect CH3CHO and long carbon-chains C6H and HC7N. (2) We find that while most COMs are segregated by type, peaking either towards the hot core (e.g., N-bearing species) or the compact ridge (e.g., O-bearing species like HCOOCH3 and CH3OCH3), while the distributions of others do not follow this segregated structure (e.g., CH3CH2OH, CH3OH, CH3COCH3). (3) We find a second velocity component of HNCO, SO2, 34SO2, and SO lines, which may be associated with a strong shock event in the low-velocity outflow. (4) Temperatures and molecular abundances show large gradients between central condensations and the outflow regions, illustrating a transition between hot molecular core and shock-chemistry dominated regimes. Our observations of spatially resolved chemical variations in Orion-KL provide the nearest reference source for hot molecular core and outflow chemistry, which will be an important example for interpreting the chemistry of more distant HMSFRs.

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

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

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.

Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its $\sim$6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e., $\sigma_{NT}/c_s\lesssim1$) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the ISM. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson's velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales.

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.

IRAS 22134+5834 was observed in the centimeter with (E)VLA, 3~mm with CARMA, 2~mm with PdBI, and 1.3~mm with SMA, to study the continuum emission as well as the molecular lines, that trace different physical conditions of the gas to study the influence of massive YSOs on nearby starless cores, and the possible implications in the clustered star formation process. The multi-wavelength centimeter continuum observations revealed two radio sources within the cluster, VLA1 and VLA2. VLA1 is considered to be an optically thin UCHII region with a size of 0.01~pc and sits at the edge of the near-infrared (NIR) cluster. The flux of ionizing photons of the VLA1 corresponds to a B1 ZAMS star. VLA2 is associated with an infrared point source and has a negative spectral index. We resolved six millimeter continuum cores at 2~mm, MM2 is associated with the UCHII region VLA1, and other dense cores are distributed around the UCH{\sc ii} region. Two high-mass starless clumps (HMSC), HMSC-E (east) and HMSC-W (west), are detected around the NIR cluster with N$_2$H$^+$(1--0) and NH$_3$ emission, and show different physical and chemical properties. Two N$_2$D$^+$ cores are detected on an NH$_3$ filament close to the UCHII region, with a projected separation of $\sim$8000~AU at the assumed distance of 2.6~kpc. The kinematic properties of the molecular line emission confirm the expansion of the UCHII region and that the molecular cloud around the near infrared (NIR) cluster is also expanding. Our multi-wavelength study has revealed different generations of star formation in IRAS 22134+5834. The formed intermediate- to massive stars show strong impact on nearby starless clumps. We propose that while the stellar wind from the UCHII region and the NIR cluster drives the large scale bubble, the starless clumps and HMPOs formed at the edge of the cluster.

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.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

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.

Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature.

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.

We report the results of our observations of the S255IR area with the SMA at 1.3 mm in the very extended configuration and at 0.8 mm in the compact configuration as well as with the IRAM-30m at 0.8 mm. The best achieved angular resolution is about 0.4 arcsec. The dust continuum emission and several tens of molecular spectral lines are observed. The majority of the lines is detected only towards the S255IR-SMA1 clump, which represents a rotating structure (probably disk) around the young massive star. The achieved angular resolution is still insufficient for conclusions about Keplerian or non-Keplerian character of the rotation. The temperature of the molecular gas reaches 130-180 K. The size of the clump is about 500 AU. The clump is strongly fragmented as follows from the low beam filling factor. The mass of the hot gas is significantly lower than the mass of the central star. A strong DCN emission near the center of the hot core most probably indicates a presence of a relatively cold ($\lesssim 80$ K) and rather massive clump there. High velocity emission is observed in the CO line as well as in lines of high density tracers HCN, HCO+, CS and other molecules. The outflow morphology obtained from combination of the SMA and IRAM-30m data is significantly different from that derived from the SMA data alone. The CO emission detected with the SMA traces only one boundary of the outflow. The outflow is most probably driven by jet bow shocks created by episodic ejections from the center. We detected a dense high velocity clump associated apparently with one of the bow shocks. The outflow strongly affects the chemical composition of the surrounding medium.

Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecular cloud that is an ideal candidate for a filament at an early stage of fragmentation. We employ dust extinction mapping in conjunction with near-infrared data from the NEWFIRM instrument, and 870 um dust continuum emission data from the LABOCA instrument, to estimate column densities. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud's morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, 1.6 pc long filament in its Center region. The line mass of the filament is 21-31 Ms pc^-1 and FWHM 0.07 pc. Its radial profile can be fitted with a Plummer profile that has the power-index of 2.6 \pm 11%, flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of 5-10 x 10^4 cm^-3. The fragments in the cloud have a mean separation of 0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the cloud ends, leaving its center yet relatively non-fragmented, possibly because of gravitational focusing in a finite geometry.

We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with HII regions. We find that the N(H2) probability density functions (N-PDFs) of these three classes of objects are clearly different: the N-PDFs of starless clumps are narrowest and close to log-normal in shape, while star-forming clouds and HII regions exhibit a power-law shape over a wide range of column densities and log-normal-like components only at low column densities. We use the N-PDFs to estimate the evolutionary time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the N-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds contain greater relative amounts of dense gas across all evolutionary classes.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

We want to understand the kinematic and thermal properties of young massive gas clumps prior to and at the earliest evolutionary stages of high-mass star formation. Do we find signatures of gravitational collapse? Do we find temperature gradients in the vicinity or absence of infrared emission sources? Do we find coherent velocity structures toward the center of the dense and cold gas clumps? To determine kinematics and gas temperatures, we used ammonia, because it is known to be a good tracer and thermometer of dense gas. We observed the NH$_3$(1,1) and (2,2) lines within seven very young high-mass star-forming regions with the VLA and the Effelsberg 100m telescope. This allows us to study velocity structures, linewidths, and gas temperatures at high spatial resolution of 3-5$"$, corresponding to $\sim$0.05 pc. We find on average cold gas clumps with temperatures in the range between 10 K and 30 K. The observations do not reveal a clear correlation between infrared emission peaks and ammonia temperature peaks. We report an upper limit for the linewidth of $\sim$1.3 km s$^{-1}$, at the spectral resolution limit of our VLA observation. This indicates a relatively low level of turbulence on the scale of the observations. Velocity gradients are present in almost all regions with typical velocity differences of 1 to 2 km s$^{-1}$ and gradients of 5 to 10 km s$^{-1}$ pc$^{-1}$. These velocity gradients are smooth in most cases, but there is one exceptional source (ISOSS23053), for which we find several velocity components with a steep velocity gradient toward the clump centers that is larger than 30 km s$^{-1}$ pc$^{-1}$. This steep velocity gradient is consistent with recent models of cloud collapse. Furthermore, we report a spatial correlation of ammonia and cold dust, but we also find decreasing ammonia emission close to infrared emission sources.

The chemical evolution in high-mass star-forming regions is still poorly constrained. Studying the evolution of deuterated molecules allows to differentiate between subsequent stages of high-mass star formation regions due to the strong temperature dependence of deuterium isotopic fractionation. We observed a sample of 59 sources including 19 infrared dark clouds, 20 high-mass protostellar objects, 11 hot molecular cores and 9 ultra-compact HII regions in the (3-2) transitions of the four deuterated molecules, DCN, DNC, DCO+ and N2D+ as well as their non-deuterated counterpart. The overall detection fraction of DCN, DNC and DCO+ is high and exceeds 50% for most of the stages. N2D+ was only detected in a few infrared dark clouds and high-mass protostellar objects. It can be related to problems in the bandpass at the frequency of the transition and to low abundances in the more evolved, warmer stages. We find median D/H ratios of ~0.02 for DCN, ~0.005 for DNC, ~0.0025 for DCO+ and ~0.02 for N2D+. While the D/H ratios of DNC, DCO+ and N2D+ decrease with time, DCN/HCN peaks at the hot molecular core stage. We only found weak correlations of the D/H ratios for N2D+ with the luminosity of the central source and the FWHM of the line, and no correlation with the H2 column density. In combination with a previously observed set of 14 other molecules (Paper I) we fitted the calculated column densities with an elaborate 1D physico-chemical model with time-dependent D-chemistry including ortho- and para-H2 states. Good overall fits to the observed data have been obtained the model. It is one of the first times that observations and modeling have been combined to derive chemically based best-fit models for the evolution of high-mass star formation including deuteration.

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

The massive infrared dark cloud G0.253+0.016 projected 45pc from the Galactic centre contains ~10^5Msun of dense gas whilst being mostly devoid of observed star-formation tracers. To scrutinise the physical properties, dynamics and structure of this cloud with reference to its star-forming potential, we have carried out a concerted SMA and IRAM 30m study of this cloud in dust continuum, CO isotopologues, shock tracing molecules, as well as H$_2$CO to trace the gas temperature. We detect and characterise the dust cores within G0.253+0.016 at ~1.3 mm and find that the kinetic temperature of the gas is >320K on size-scales of ~0.15 pc. Analysis of the position-velocity diagrams of our observed lines show broad linewidths and strong shock emission in the south of the cloud, indicating that G0.253+0.016 is colliding with another cloud at v(LSR)~70 km/s. We confirm via an analysis of the observed dynamics in the CMZ that it is an elongated structure, orientated with Sgr B2 closer to the Sun than Sgr A*, however our results suggest that the actual geometry may be more complex than an elliptical ring. We find that the column density PDF of G0.253+0.016 is log-normal with no discernible power-law tail, consistent with little star formation, and that its width can be explained in the framework of theory predicting the density structure of clouds created by supersonic, magnetised turbulence. We also present the delta-variance spectrum of this region, and show it is consistent with that expected for clouds with no star formation. Using G0.253+0.016 as a test-bed of the conditions required for star formation in a different physical environment to that of nearby clouds, we also conclude that there is not one column density threshold for star formation, but instead this value is dependant on the local physical conditions. [Abbrv.]

We present new Plateau de Bure Interferometer observations of a region in the filamentary infrared-dark cloud (IRDC) G011.11-0.12 containing young, star-forming cores. In addition to the 3.2mm continuum emission from cold dust, we map this region in the N$_2$H$^+$(1-0) line to trace the core kinematics with an angular resolution of 2" and velocity resolution of 0.2km s$^{-1}$. These data are presented in concert with recent {\em Herschel} results, single-dish N$_2$H$^+$(1-0) data, SABOCA 350$\mu$m continuum data, and maps of the C$^{18}$O (2-1) transition obtained with the IRAM 30m telescope. We recover the star-forming cores at 3.2mm continuum, while in N$_2$H$^+$ they appear at the peaks of extended structures. The mean projected spacing between N$_2$H$^+$ emission peaks is 0.18pc, consistent with simple isothermal Jeans fragmentation. The 0.1pc-sized cores have low virial parameters on the criticality borderline, while on the scale of the whole region, we infer that it is undergoing large-scale collapse. The N$_2$H$^+$ linewidth increases with evolutionary stage, while CO isotopologues show no linewidth variation with core evolution. Centroid velocities of all tracers are in excellent agreement, except in the starless region where two N$_2$H$^+$ velocity components are detected, one of which has no counterpart in C$^{18}$O. We suggest that gas along this line of sight may be falling into the quiescent core, giving rise to the second velocity component, possibly connected to the global collapse of the region.

Context: How do molecular clouds form out of the atomic phase? And what are the relative fractions of carbon in the ionized, atomic and molecular phase? These are questions at the heart of cloud and star formation. Methods: Using multiple observatories from Herschel and SOFIA to APEX and the IRAM 30m telescope, we mapped the ionized, atomic and molecular carbon ([CII]@1900GHz, [CI]@492GHz and C18O(2-1)@220GHz) at high spatial resolution (12"-25") in four young massive infrared dark clouds (IRDCs). Results: The three carbon phases were successfully mapped in all four regions, only in one source the [CII] line remained a non-detection. Both the molecular and atomic phases trace the dense structures well, with [CI] also tracing material at lower column densities. [CII] exhibits diverse morphologies in our sample, from compact to diffuse structures probing the cloud environment. In at least two out of the four regions, we find kinematic signatures strongly indicating that the dense gas filaments have formed out of a dynamically active and turbulent atomic/molecular cloud, potentially from converging gas flows. The atomic-to-molecular carbon gas mass ratios are low between 7% and 12% with the lowest values found toward the most quiescent region. In the three regions where [CII] is detected, its mass is always higher by a factor of a few than that of the atomic carbon. The ionized carbon emission depends as well on the radiation field, however, we also find strong [CII] emission in a region without significant external sources, indicating that other processes, e.g., energetic gas flows can contribute to the [CII] excitation as well.

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.

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.

Throughout the Milky Way, molecular clouds typically appear filamentary, and mounting evidence indicates that this morphology plays an important role in star formation. What is not known is to what extent the dense filaments most closely associated with star formation are connected to the surrounding diffuse clouds up to arbitrarily large scales. How are these cradles of star formation linked to the Milky Way's spiral structure? Using archival Galactic plane survey data, we have used multiple datasets in search of large-scale, velocity-coherent filaments in the Galactic plane. In this paper, we present our methods employed to identify coherent filamentary structures first in extinction and confirmed using Galactic Ring Survey data. We present a sample of seven Giant Molecular Filaments (GMFs) that have lengths of order $\sim$100 pc, total masses of 10$^4$ - 10$^5$ M$_{\odot}$, and exhibit velocity coherence over their full length. The GMFs we study appear to be inter-arm clouds and may be the Milky Way analogues to spurs observed in nearby spiral galaxies. We find that between 2 and 12% of the total mass (above $\sim$10$^{20}$ cm$^{-2}$) is "dense" (above 10$^{22}$ cm$^{-2}$), where filaments near spiral arms in the Galactic midplane tend to have higher dense gas mass fractions than those further from the arms.

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.

Numerical simulations have explored the possibility to form molecular clouds through either a quasi-static, self-gravitating mechanism or the collision of gas streams or lower-density clouds. They also quantitatively predict the distribution of matter at the transition from atomic to molecular gases. We aim to observationally test these models by studying the environment of W43, a molecular cloud complex near the tip of the Galactic long bar. Using Galaxy-wide HI and 12CO surveys we searched for gas flowing toward the W43 molecular cloud complex. We also estimated the HI and H2 mass surface densities to constrain the transition from atomic to molecular gas around and within W43. We found 3 cloud ensembles within the position-velocity diagrams of 12CO and HI gases. They are separated by 20km/s along the line of sight and extend into the 13CO velocity structure of W43. Since their velocity gradients are consistent with free-fall, they could be nearby clouds attracted by, and streaming toward, the W43 10^7Msun potential well. We show that the HI surface density, Sigma_HI=45-85Msun/pc2, does not reach any threshold level but increases when entering the 130pc-wide molecular complex previously defined. This suggests that an equilibrium between H2 formation and photodissociation has not yet been reached. The H2-to-HI ratio measured over the W43 region and its surroundings, R_H2~3.5, is high, indicating that most of the gas is already in molecular form in W43 and in structures several hundreds of parsecs downstream along the Scutum-Centaurus arm. The W43 molecular cloud complex may have formed, and in fact may still be accreting mass from the agglomeration of clouds. Already in the molecular-dominated regime, most of these clouds are streaming from the Scutum-Centaurus arm. This is in clear disagreement with quasi-static and steady-state models of molecular cloud formation.

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.

Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying physical processes and 2) to characterize the evolution of the chemical composition. We observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact HII regions at 1mm and 3mm with the IRAM 30m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit chi^2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10,000 years for the IRDC stage, ~60,000 years for the HMPO stage, ~40,000 years for the HMC stage, and ~10,000 years for the UCHII stage. The total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 10^5 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.

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.

The mid- and far-infrared view on high-mass star formation, in particular with the results from the Herschel space observatory, has shed light on many aspects of massive star formation. However, these continuum studies lack kinematic information. We study the kinematics of the molecular gas in high-mass star-forming regions. We complemented the PACS and SPIRE far-infrared data of 16 high-mass star-forming regions from the Herschel key project EPoS with N2H+ molecular line data from the MOPRA and Nobeyama 45m telescope. Using the full N2H+ hyperfine structure, we produced column density, velocity, and linewidth maps. These were correlated with PACS 70micron images and PACS point sources. In addition, we searched for velocity gradients. For several regions, the data suggest that the linewidth on the scale of clumps is dominated by outflows or unresolved velocity gradients. IRDC18454 and G11.11 show two velocity components along several lines of sight. We find that all regions with a diameter larger than 1pc show either velocity gradients or fragment into independent structures with distinct velocities. The velocity profiles of three regions with a smooth gradient are consistent with gas flows along the filament, suggesting accretion flows onto the densest regions. We show that the kinematics of several regions have a significant and complex velocity structure. For three filaments, we suggest that gas flows toward the more massive clumps are present.

Infrared-dark clouds (IRDCs) harbor the early phases of cluster and high-mass star formation and are comprised of cold (~20 K), dense (n > 10$^4$ cm$^{-3}$) gas. The spectral energy distribution (SED) of IRDCs is dominated by the far-infrared and millimeter wavelength regime, and our initial Herschel study examined IRDCs at the peak of the SED with high angular resolution. Here we present a follow-up study using the SABOCA instrument on APEX which delivers 7.8" angular resolution at 350 micron, matching the resolution we achieved with Herschel/PACS, and allowing us to characterize substructure on ~0.1pc scales. Our sample of 11 nearby IRDCs are a mix of filamentary and clumpy morphologies, and the filamentary clouds show significant hierarchical structure, while the clumpy IRDCs exhibit little hierarchical structure. All IRDCs, regardless of morphology, have about 14% of their total mass in small scale core-like structures which roughly follow a trend of constant volume density over all size scales. Out of the 89 protostellar cores we identified in this sample with Herschel, we recover 40 of the brightest and re-fit their SEDs and find their properties agree fairly well with our previous estimates (<T> ~ 19K). We detect a new population of "cold cores" which have no 70 micron counterpart, but are 100 and 160 micron-bright, with colder temperatures (<T> ~ 16K). This latter population, along with SABOCA-only detections, are predominantly low-mass objects, but their evolutionary diagnostics are consistent with the earliest starless or prestellar phase of cores in IRDCs.

Embedded clusters like W3 Main are complex and dynamically evolving systems that represent an important phase of the star formation process. We aim at the characterization of the entire stellar content of W3 Main in a statistical sense to identify possible differences in evolutionary phase of the stellar populations and find clues about the formation mechanism of this massive embedded cluster. Methods. Deep JHKs imaging is used to derive the disk fraction, Ks-band luminosity functions and mass functions for several subregions in W3 Main. A two dimensional completeness analysis using artificial star experiments is applied as a crucial ingredient to assess realistic completeness limits for our photometry. We find an overall disk fraction of 7.7 $\pm$ 2.3%, radially varying from 9.4 $\pm$ 3.0 % in the central 1 pc to 5.6 $\pm$ 2.2 % in the outer parts of W3 Main. The mass functions derived for three subregions are consistent with a Kroupa and Chabrier mass function. The mass function of IRSN3 is complete down to 0.14 Msun and shows a break at M $\sim$ 0.5 Msun. We interpret the higher disk fraction in the center as evidence for a younger age of the cluster center. We find that the evolutionary sequence observed in the low-mass stellar population is consistent with the observed age spread among the massive stars. An analysis of the mass function variations does not show evidence for mass segregation. W3 Main is currently still actively forming stars, showing that the ionizing feedback of OB stars is confined to small areas ($\sim$ 0.5 pc). The FUV feedback might be influencing large regions of the cluster as suggested by the low overall disk fraction.

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.

Aims: Revealing the fragmentation, infall, and outflow processes in the immediate environment around massive young stellar objects is crucial for understanding the formation of the most massive stars. Methods: With this goal in mind we present the so far highest spatial-resolution thermal submm line and continuum observations toward the young high-mass protostar NGC7538 IRS1. Using the Plateau de Bure Interferometer in its most extended configuration at 843mum wavelength, we achieved a spatial resolution of 0.2"x0.17", corresponding to ~500AU at a distance of 2.7\,kpc. Results: For the first time, we have observed the fragmentation of the dense inner core of this region with at least three subsources within the inner 3000 AU. The outflow exhibits blue- and red-shifted emission on both sides of the central source indicating that the current orientation has to be close to the line-of-sight, which differs from other recent models. We observe rotational signatures in northeast-southwest direction; however, even on scales of 500 AU, we do not identify any Keplerian rotation signatures. This implies that during the early evolutionary stages any stable Keplerian inner disk has to be very small (<=500 AU). The high-energy line HCN(4-3)v2=1 (E_u/k=1050K) is detected over an extent of approximately 3000 AU. In addition to this, the detection of red-shifted absorption from this line toward the central dust continuum peak position allows us to estimate infall rates of ~1.8x10^(-3)Msun/yr on the smallest spatial scales. Although all that gas will not necessarily be accreted onto the central protostar, nevertheless, such inner core infall rates are among the best proxies of the actual accretion rates one can derive during the early embedded star formation phase. These data are consistent with collapse simulations and the observed high multiplicity of massive stars.

We present near-infrared spectroscopy and 1 mm line and continuum observations of a recently identified site of high mass star formation likely to be located in the Central Molecular Zone near Sgr C. Located on the outskirts of the massive evolved HII region associated with Sgr C, the area is characterized by an Extended Green Object measuring ~10" in size (0.4 pc), whose observational characteristics suggest the presence of an embedded massive protostar driving an outflow. Our data confirm that early-stage star formation is taking place on the periphery of the Sgr C HII region, with detections of two protostellar cores and several knots of H2 and Brackett gamma emission alongside a previously detected compact radio source. We calculate the cores' joint mass to be ~10^3 Msun, with column densities of 1-2 x 10^24 cm-2. We show the host molecular cloud to hold ~10^5 Msun of gas and dust with temperatures and column densities favourable for massive star formation to occur, however, there is no evidence of star formation outside of the EGO, indicating that the cloud is predominantly quiescent. Given its mass, density, and temperature, the cloud is comparable to other remarkable non-star-forming clouds such as G0.253 in the Eastern CMZ.

Observations of atomic or molecular lines can provide important information about the physical state of star forming regions. In order to investigate the line profiles from dynamical collapsing massive star forming regions (MSFRs), we model the emission from hydrodynamic simulations of a collapsing cloud in the absence of outflows. By performing radiative transfer calculations, we compute the optically thick HCO+ and optically thin N2H+ line profiles from two collapsing regions at different epochs. Due to large-scale collapse, the MSFRs have large velocity gradients, reaching up to 20 km/s/pc across the central core. The optically thin lines typically contain multiple velocity components resulting from the superposition of numerous density peaks along the line-of-sight. The optically thick lines are only marginally shifted to the blue side of the optically thin line profiles, and frequently do not have a central depression in their profiles due to self-absorption. As the regions evolve the lines become brighter and the optically thick lines become broader. The lower order HCO+ (1-0) transitions are better indicators of collapse than the higher order (4-3) transitions. We also investigate how the beam sizes affect profile shapes. Smaller beams lead to brighter and narrower lines that are more skewed to the blue in HCO+ relative to the true core velocity, but show multiple components in N2H+. High resolution observations (e.g. with ALMA) can test these predictions and provide insights into the nature of MSFRs.

Aims: We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. Methods: Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100mum wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45m single-dish spectral line and continuum observations. Results: The massive gas reservoir fragments at spatial scales of ~18000AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation properties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. Conclusions: This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic collapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings.

Despite increasing research in massive star formation, little is known about its earliest stages. Infrared Dark Clouds (IRDCs) are cold, dense and massive enough to harbour the sites of future high-mass star formation. But up to now, mainly small samples have been observed and analysed. To understand the physical conditions during the early stages of high-mass star formation, it is necessary to learn more about the physical conditions and stability in relatively unevolved IRDCs. Thus, for characterising IRDCs studies of large samples are needed. We investigate a complete sample of 218 northern hemisphere high-contrast IRDCs using the ammonia (1,1)- and (2,2)-inversion transitions. We detected ammonia (1,1)-inversion transition lines in 109 of our IRDC candidates. Using the data we were able to study the physical conditions within the star-forming regions statistically. We compared them with the conditions in more evolved regions which have been observed in the same fashion as our sample sources. Our results show that IRDCs have, on average, rotation temperatures of 15 K, are turbulent (with line width FWHMs around 2 km s$^{-1}$), have ammonia column densities on the order of $10^{14}$ cm$^{-2}$ and molecular hydrogen column densities on the order of $10^{22}$ cm$^{-2}$. Their virial masses are between 100 and a few 1000 M$_\odot$. The comparison of bulk kinetic and potential energies indicate that the sources are close to virial equilibrium. IRDCs are on average cooler and less turbulent than a comparison sample of high-mass protostellar objects, and have lower ammonia column densities. Virial parameters indicate that the majority of IRDCs are currently stable, but are expected to collapse in the future.

The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

We present continuum observations of the infrared dark cloud (IRDC) G48.66-0.22 (G48) obtained with Herschel, Spitzer, and APEX, in addition to several molecular line observations. The Herschel maps are used to derive temperature and column density maps of G48 using a model based on a modified blackbody. We find that G48 has a relatively simple structure and is relatively isolated; thus this IRDC provides an excellent target to study the collapse and fragmentation of a filamentary structure in the absence of complicating factors such as strong external feedback. The derived temperature structure of G48 is clearly non-isothermal from cloud to core scale. The column density peaks are spatially coincident with the lowest temperatures (~ 17.5 K) in G48. A total cloud mass of ~390Msun is derived from the column density maps. By comparing the luminosity-to-mass ratio of 13 point sources detected in the Herschel/PACS bands to evolutionary models, we find that two cores are likely to evolve into high-mass stars (M>8 Msun). The derived mean projected separation of point sources is smaller than in other IRDCs but in good agreement with theoretical predications for cylindrical collapse. We detect several molecular species such as CO, HCO+, HCN, HNC and N2H+. CO is depleted by a factor of ~3.5 compared to the expected interstellar abundance, from which we conclude that CO freezes out in the central region. Furthermore, the molecular clumps, associated with the sub-millimeter peaks in G48, appear to be gravitationally unbound or just pressure confined. The analysis of critical line masses in G48 show that the entire filament is collapsing, overcoming any internal support.

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

We review our recent studies demonstrating that the radiation pressure problem in the formation of massive stars can be circumvented via an anisotropy of the thermal radiation field. Such an anisotropy naturally establishes with the formation of a circumstellar disk. The required angular momentum transport within the disk can be provided by developing gravitational torques. Radiative Rayleigh-Taylor instabilities in the cavity regions - as previously suggested in the literature - are not required and are shown to be not occurring in the context of massive star formation.

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.

(Abriged) In the framework of the Herschel GTKP "The earliest phases of star formation", we have imaged B68 between 100 and 500 um. Ancillary (sub)mm data, spectral line maps of the 12/13CO(2-1) transitions as well as a NIR extinction map were added to the analysis. We employed a ray-tracing algorithm to derive the 2D mid-plane dust temperature and volume density distribution without suffering from LoS averaging effects of simple SED fitting procedures. Additional 3D radiative transfer calculations were employed to investigate the connection between the external irradiation and the peculiar crescent shaped morphology found in the FIR maps. For the first time, we spatially resolve the dust temperature and density distribution of B68. We find T_dust dropping from 16.7 K at the edge to 8.2 K in the centre, which is about 4 K lower than the result of the simple SED fitting approach. N_H peaks at 4.3x10^22 cm^-2 and n_H at 3.4x10^5 cm^-3 in the centre. B68 has a mass of 3.1 M_sun of material with A_K > 0.2 mag for an assumed distance of 150 pc. We detect a compact source in the southeastern trunk, which is also seen in extinction and CO. We find the radial density distribution from the edge of the inner plateau outward to be n_H ~ r^-3.5. Such a steep profile can arise from either or both of the following: external irradiation with a significant UV contribution or the fragmentation of filamentary structures. Our 3D radiative transfer model of an externally irradiated core by an anisotropic ISRF reproduces the crescent morphology. Our CO observations show that B68 is part of a chain of globules in both space and velocity, which may indicate that it was once part of a filament which dispersed. We also resolve a new compact source in the SE trunk and find that it is slightly shifted in centroid velocity from B68, lending qualitative support to core collision scenarios.

(Abridged) We present an overview of the sample of high-mass star and cluster forming regions observed as part of the Earliest Phases of Star Formation (EPoS) Herschel Guaranteed Time Key Program. A sample of 45 infrared-dark clouds (IRDCs) were mapped at PACS 70, 100, and 160 micron and SPIRE 250, 350, and 500 micron. In this paper, we characterize a population of cores which appear in the PACS bands and place them into context with their host cloud and investigate their evolutionary stage. We construct spectral energy distributions (SEDs) of 496 cores which appear in all PACS bands, 34% of which lack counterparts at 24 micron. From single-temperature modified blackbody fits of the SEDs, we derive the temperature, luminosity, and mass of each core. These properties predominantly reflect the conditions in the cold, outer regions. Taking into account optical depth effects and performing simple radiative transfer models, we explore the origin of emission at PACS wavelengths. The core population has a median temperature of 20K and has masses and luminosities that span four to five orders of magnitude. Cores with a counterpart at 24 micron are warmer and bluer on average than cores without a 24 micron counterpart. We conclude that cores bright at 24 micron are on average more advanced in their evolution, where a central protostar(s) have heated the outer bulk of the core, than 24 micron-dark cores. The 24 micron emission itself can arise in instances where our line of sight aligns with an exposed part of the warm inner core. About 10% of the total cloud mass is found in a given cloud's core population. We uncover over 300 further candidate cores which are dark until 100 micron. These are candidate starless objects, and further observations will help us determine the nature of these very cold cores.

We present a study of outflow, infall, and rotation in a ~10^5 Lsun (solar luminosity) star-forming region, IRAS 18360-0537, with Submillimeter Array (SMA) and IRAM 30m observations. The 1.3 mm continuum map shows a 0.5 pc dust ridge, of which the central compact part has a mass of ~80 Msun (solar mass) and harbors two condensations, MM1 and MM2. The CO (2--1) and SiO (5--4) maps reveal a biconical outflow centered at MM1, which is a hot molecular core (HMC) with a gas temperature of 320+/-50 K and a mass of ~13 Msun. The outflow has a gas mass of 54 Msun and a dynamical timescale of 8,000 yr. The kinematics of the HMC is probed by high-excitation CH3OH and CH3CN lines, which are detected at sub-arcsecond resolution and unveil a velocity gradient perpendicular to the outflow axis, suggesting a disk-like rotation of the HMC. An infalling envelope around the HMC is evidenced by CN lines exhibiting a profound inverse P-Cygni profile, and the estimated mass infall rate, 1.5x10^{-3} Msun/yr, is well comparable to that inferred from the mass outflow rate. A more detailed investigation of the kinematics of the dense gas around the HMC is obtained from the 13CO and C18O (2--1) lines; the position-velocity diagrams of the two lines are consistent with the model of a free-falling and Keplerian-like rotating envelope. The observations suggest that the protostar of a current mass ~10 Msun embedded within MM1 will develop into an O star via disk accretion and envelope infall.

The Core Mass Functions (CMFs) of low-mass star-forming regions are found to resemble the shape of the Initial Mass Function (IMF). A similar result is observed for the dust clumps in high-mass star forming regions, although at spatial scales of clusters that do not resolve the substructure found in them. The region IRAS 19410+2336 is one exception, having been observed at spatial scales on the order of ~2500AU, resolving the clump substructure into individual cores. We mapped that region with the PdBI in the 1.4mm and 3mm continuum and several transitions of H2CO and CH3CN. The H2CO transitions were also observed with the IRAM 30m Telescope. We detected 26 continuum sources at 1.4mm with a spatial resolution down to ~2200 AU, distributed in two protoclusters. With the lines emission we derived the temperature structure of the region, ranging from 35 to 90K. With them we calculated the core masses of the detected sources, ranging from ~0.7 to ~8 M_sun. These masses were strongly (~90%) affected by the interferometer spatial filtering. Considering only the detected dense cores we derived a CMF with a power-law index b=-2.3+-0.2. We resolve the Jeans length of the protoclusters by one order of magnitude, and only find little velocity dispersion between the different subsources. Since we cannot unambiguously differentiate protostellar and prestellar cores, the derived CMF is not prestellar. Also, because of the large missing flux, we cannot establish a firm link between the CMF and the IMF. This implies that future high-mass CMF studies will need to complement the interferometer continuum data with the short spacing data, a task suitable for ALMA. We note that the method of extracting temperatures using H2CO lines becomes less applicable when reaching the dense core scales of the interferometric observations because most of the H2CO appears to originate in the envelope structure.

We present near-infrared JHKs imaging as well as K-band multi-object spectroscopy of the massive stellar content of W3 Main using LUCI at the LBT. We confirm 13 OB stars by their absorption line spectra in W3 Main and spectral types between O5V and B4V have been found. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find an age spread of 2-3 Myr.

We observed three high-mass star-forming regions in the W3 high-mass star formation complex with the Submillimeter Array and IRAM 30 m telescope. These regions, i.e. W3 SMS1 (W3 IRS5), SMS2 (W3 IRS4) and SMS3, are in different evolutionary stages and are located within the same large-scale environment, which allows us to study rotation and outflows as well as chemical properties in an evolutionary sense. While we find multiple mm continuum sources toward all regions, these three sub-regions exhibit different dynamical and chemical properties, which indicates that they are in different evolutionary stages. Even within each subregion, massive cores of different ages are found, e.g. in SMS2, sub-sources from the most evolved UCHII region to potential starless cores exist within 30 000 AU of each other. Outflows and rotational structures are found in SMS1 and SMS2. Evidence for interactions between the molecular cloud and the HII regions is found in the 13CO channel maps, which may indicate triggered star formation.

Context: The nature of embedded accretion disks around forming high-mass stars is one of the missing puzzle pieces for a general understanding of the formation of the most massive and luminous stars. Methods: Using the Plateau de Bure Interferometer at 1.36mm wavelengths in its most extended configuration we probe the dust and gas emission at ~0.3",corresponding to linear resolution elements of ~800AU. Results: NGC7538IRS1 remains a single compact and massive gas core with extraordinarily high column densities, corresponding to visual extinctions on the order of 10^5mag, and average densities within the central 2000AU of ~2.1x10^9cm^-3 that have not been measured before. We identify a velocity gradient across in northeast-southwest direction that is consistent with the mid-infrared emission, but we do not find a gradient that corresponds to the proposed CH3OH maser disk. The spectral line data toward NGC7538IRS1 reveal strong blue- and red-shifted absorption toward the mm continuum peak position. The red-shifted absorption allows us to estimate high infall rates on the order of 10^-2 Msun/yr. Although we cannot prove that the gas will be accreted in the end, the data are consistent with ongoing star formation activity in a scaled-up low-mass star formation scenario. Compared to that, NGC7538S fragments in a hierarchical fashion into several sub-sources. While the kinematics of the main mm peak are dominated by the accompanying jet, we find rotational signatures from a secondary peak. Furthermore, strong spectral line differences exist between the sub-sources which is indicative of different evolutionary stages within the same large-scale gas clump.

We have started a campaign to identify massive star clusters inside bright molecular bubbles towards the Galactic Center. The CN15/16/17 molecular complex is the first example of our study. The region is characterized by the presence of two young clusters, DB10 and DB11, visible in the NIR, an ultra-compact HII region identified in the radio, several young stellar objects visible in the MIR, a bright diffuse nebulosity at 8\mu m coming from PAHs and sub-mm continuum emission revealing the presence of cold dust. Given its position on the sky (l=0.58, b=-0.85) and its kinematic distance of ~7.5 kpc, the region was thought to be a very massive site of star formation in proximity of the CMZ. The cluster DB11 was estimated to be as massive as 10^4 M_sun. However the region's properties were known only through photometry and its kinematic distance was very uncertain given its location at the tangential point. We aimed at better characterizing the region and assess whether it could be a site of massive star formation located close to the Galactic Center. We have obtained NTT/SofI JHKs photometry and long slit K band spectroscopy of the brightest members. We have additionally collected data in the radio, sub-mm and mid infrared, resulting in a quite different picture of the region. We have confirmed the presence of massive early B type stars and have derived a spectro-photometric distance of ~1.2 kpc, much smaller than the kinematic distance. Adopting this distance we obtain clusters masses of M(DB10) ~ 170 M_sun and M(DB11) ~ 275 M_sun. This is consistent with the absence of any O star, confirmed by the excitation/ionization status of the nebula. No HeI diffuse emission is detected in our spectroscopic observations at 2.113\mu m, which would be expected if the region was hosting more massive stars. Radio continuum measurements are also consistent with the region hosting at most early B stars.

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.

Context: The earliest stages of high-mass star formation are still poorly characterized. It is also unknown whether the initial conditions vary with environment. Aims: We want to investigate the youngest massive gas clumps in the environment of extremely active star formation. Methods: We selected the IRDC18454 complex associated with the W43 Galactic mini-starburst, and observed it in the continuum emission between 70mu and 1.2mm, and in spectral line emission of N2H+ and 13CO with single-dish instruments and an interferometer. Results: The multi-wavelength continuum study allows us to identify clumps that are infrared dark even at 70mu and hence the best candidates to be genuine high-mass starless gas clumps. The spectral energy distributions reveal elevated temperatures and luminosities compared to more quiescent environments. Furthermore, we identify a temperature gradient from the W43 mini-starburst toward the starless clumps. We discuss whether the radiation impact of the nearby mini-starburst changes the fragmentation properties of the gas clumps and by that maybe favors more high-mass star formation in such an environment. The spectral line data reveal two different velocity components of the gas at 100 and 50km/s. While chance projection is a possibility to explain these components, the projected associations of the emission sources as well as the prominent location at the Galactic bar -- spiral arm interface also allow the possibility that these two components may be spatially associated and even interacting. Conclusions: High-mass starless gas clumps can exist in the close environment of very active star formation without being destroyed. The impact of the active star formation sites may even allow for more high-mass stars to form in these 2nd generation gas clumps.

(Context) Many physical parameters change with time in star forming regions. Here we attempt to correlate changes in infall and outflow motions in high mass star forming regions with evolutionary stage using JCMT observations. (Aims) From a sample of 45 high mass star forming regions in three phases of evolution, we investigate the presence of established infall and outflow tracers to determine whether there are any trends attributable to the age of the source. (Methods) We obtained JCMT observations of HCO+/H13CO+ J=4-3 to trace large scale infall, and SiO J=8-7 to trace recent outflow activity. We compare the infall and outflow detections to the evolutionary stage of the host source (high mass protostellar objects, hypercompact HII regions and ultracompact HII regions). We also note that the integrated intensity of SiO varies with the full width at half maximum of the H13CO+. (Results) We find a surprising lack of SiO detections in the middle stage (Hypercompact HII regions), which may be due to an observational bias. When SiO is detected, we find that the integrated intensity of the line increases with evolutionary stage. We also note that all of the sources with infall signatures onto Ultracompact HII regions have corresponding outflow signatures as well.

Context: When massive stars exert a radiation pressure onto their environment that is higher than their gravitational attraction, they launch a radiation-pressure-driven outflow. It has been claimed that a radiative Rayleigh-Taylor instability should lead to the collapse of the outflow cavity and foster the growth of massive stars. Aims: We investigate the stability of idealized radiation-pressure-dominated cavities, focusing on its dependence on the radiation transport approach for the stellar radiation feedback. Methods: We compare two different methods for stellar radiation feedback: gray flux-limited diffusion (FLD) and ray-tracing (RT). We also derive simple analytical models to support our findings. Results: Only the FLD cases lead to prominent instability in the cavity shell. The RT cases do not show such instability. The gray FLD method underestimates the opacity at the location of the cavity shell and leads to extended epochs of marginal Eddington equilibrium in the cavity shell, making them prone to the radiative Rayleigh-Taylor instability. In the RT cases, the radiation pressure exceeds gravity by 1-2 orders of magnitude. The radiative Rayleigh-Taylor instability is then consequently suppressed. Conclusions: Treating the stellar irradiation in the gray FLD approximation underestimates the radiative forces acting on the cavity shell. This can lead artificially to situations that are affected by the radiative Rayleigh-Taylor instability. The proper treatment of direct stellar irradiation by massive stars is crucial for the stability of radiation-pressure-dominated cavities.

Context: Infrared Dark Clouds (IRDCs) harbor progenitors of high-mass stars. Little is known of the parental molecular clouds of the IRDCs. Aims: We demonstrate the feasibility of the near-infrared (NIR) dust extinction mapping in tracing the parental molecular clouds of IRDCs at the distances of D = 2.5 - 8 kpc. Methods: We derive NIR extinction maps for 10 prominent IRDC complexes using a color-excess mapping technique and NIR data from the UKIDSS/Galactic Plane Survey. We compare the resulting maps to the 13CO emission line data, to the 8 \mu m dust opacity data, and to the millimeter dust emission data. We derive distances for the clouds by comparing the observed NIR source densities to the Besancon stellar distribution model and compare them to the kinematic distance estimates. Results: The NIR extinction maps provide a view to the IRDC complexes over the dynamical range of Av = 2 - 40 mag, in spatial resolution of 30". The NIR extinction data correlate well with the 13CO data and probe a similar gas component, but also extend to higher column densities. The NIR data reveal a wealth of extended structures surrounding the dense gas traced by the 8 \mu m shadowing features and sub-mm dust emission, showing that the clouds contain typically > 10 times more mass than traced by those tracers. The IRDC complexes of our sample contain relatively high amount of high-column density material, and their cumulative column density distributions resemble active nearby star-forming clouds like Orion rather than less active clouds like California. Conclusions: NIR dust extinction data provide a new powerful tool to probe the mass distribution of the parental molecular clouds of IRDCs up to the distances of D = 8 kpc. This encourages for deeper NIR observations of IRDCs, because the sensitivity and resolution of the data can be directly enhanced with dedicated observations.

We present near-infrared multi-object spectroscopy and JHKs imaging of the massive stellar content of the Galactic star-forming region W3 Main, obtained with LUCI at the Large Binocular Telescope. We confirm 15 OB stars in W3 Main and derive spectral types between O5V and B4V from their absorption line spectra. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. The color-color diagram of the detected sources allows a detailed investigation of the slope of the near-infrared extinction law towards W3 Main. Analysis of the Hertzsprung Russell diagram suggests that the Nishiyama extinction law fits the stellar population of W3 Main best (E(J-H)/E(H-Ks) = 1.76 and R_(Ks) = 1.44). From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive star (IRS2) is already evolved, indications for high-mass pre-main-sequence evolution is found for another star (IRS N1), deeply embedded in an ultra compact HII region, in line with the different evolutionary phases observed in the corresponding HII regions. We derive a stellar mass of W3 Main of (4 +- 1) 10^3 Msun, by extrapolating from the number of OB stars using a Kroupa IMF and correcting for our spectroscopic incompleteness. We have detected the photospheres of OB stars from the more evolved diffuse HII region to the much younger UCHII regions, suggesting that these stars have finished their formation and cleared away their circumstellar disks very fast. Only in the hyper-compact HII region (IRS5), the early type stars seem to be still surrounded by circumstellar material.

Observations indicate that outflows from massive young stars are more collimated during their early evolution compared to later stages. Our paper investigates various physical processes that impacts the outflow dynamics, i.e. its acceleration and collimation. We perform axisymmetric MHD simulations particularly considering the radiation pressure exerted by the star and the disk. We have modified the PLUTO code to include radiative forces in the line-driving approximation. We launch the outflow from the innermost disk region (r < 50 AU) by magneto-centrifugal acceleration. In order to disentangle MHD effects from radiative forces, we start the simulation in pure MHD, and later switch on the radiation force. We perform a parameter study considering different stellar masses (thus luminosity), magnetic flux, and line-force strength. For our reference simulation - assuming a 30 Msun star, we find substantial de-collimation of 35 % due to radiation forces. The opening angle increases from 20 deg to 32 deg for stellar masses from 20 Msun to 60 Msun. A small change in the line-force parameter 'alpha' from 0.60 to 0.55 changes the opening angle by ~ 8 deg. We find that it is mainly the stellar radiation which affects the jet dynamics. Unless the disk extends very close to the star, its pressure is too small to have much impact. Essentially, our parameter runs with different stellar mass can be understood as a proxy for the time evolution of the star-outflow system. Thus, we have shown that when the stellar mass (thus luminosity) increases (with age), the outflows become less collimated.

We made a detailed observational analysis of a well known hot molecular core lying in the high-mass star-forming region G31.41+0.31. This core is believed to contain deeply embedded massive stars and presents a velocity gradient that has been interpreted either as rotation or as expansion, depending on the authors. Our aim was to shed light on this question and possibly prepare the ground for higher resolution ALMA observations which could directly detect circumstellar disks around the embedded massive stars. Observations at sub-arcsecond resolution were performed with the Submillimeter Array in methyl cyanide, a typical hot molecular core tracer, and 12CO and 13CO, well known outflow tracers. We also obtained sensitive continuum maps at 1.3 mm. Our findings confirm the existence of a sharp velocity gradient across the core, but cannot confirm the existence of a bipolar outflow perpendicular to it. The improved angular resolution and sampling of the uv plane allow us to attain higher quality channel maps of the CH3CN lines with respect to previous studies and thus significantly improve our knowledge of the structure and kinematics of the hot molecular core. While no conclusive argument can rule out any of the two interpretations (rotation or expansion) proposed to explain the velocity gradient observed in the core, in our opinion the observational evidence collected so far indicates the rotating toroid as the most likely scenario. The outflow hypothesis appears less plausible, because the dynamical time scale is too short compared to that needed to form species such as CH3CN, and the mass loss and momentum rates estimated from our measurements appear too high.

Context: The large Coalsack dark cloud is one of the rare molecular clouds without clear signs of star formation. Aims: We investigate the dynamical properties of the gas within the Coalsack. Methods: The two highest extinction regions were mapped with the APEX telescope in 13CO(2-1) comprising a region of ~1 square degree. Results: In addition to the well-known, nearby gas component around -4km/s, we identified additional molecular gas components - in particular a second extended molecular cloud at a velocity of ~-30km/s and an estimated distance of ~3.1kpc - that dominate the column density and visual extinction distributions in the northeastern part of the Coalsack. Although comprising ~2600Msun, the mass of this distant cloud is distributed over an extent of ~73pc, much larger than typical high-mass infrared dark clouds. Its filamentary structure is consistent with a compressible gaseous self-gravitating cylinder, and its low mass per length indicates that it may be stable against gravitational collapse. We find barely any mid-infrared emission in archival MSX data, which is indicative of almost no star-formation activity in the near and far cloud complexes. The nearby clouds have narrow, almost thermal velocity dispersions with median values between 0.2 and 0.4km/s, which is also consistent with low star-formation activity. Only Tapia's Globule 2 exhibits a velocity dispersion increase toward the extinction peak and peak-velocity gradients over the core, which is indicative of a state of elevated dynamical properties. Conclusions: The Coalsack is not one single coherent structure, but consists of several cloud complexes nearby as well as at several kpc distance. All studied clouds appear as starless low-turbulence regions that may not even collapse in the future. Only one globule exhibits more dynamical signatures and is a good candidate for present/future star formation.

We present near-infrared JHKs imaging as well as K-band multi-object spectroscopy of the massive stellar content of W3 Main using LUCI at the LBT. We confirm 13 OB stars by their absorption line spectra in W3 Main and spectral types between O5V and B4V have been found. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive star (IRS2) is already evolved, indications for high-mass pre--main-sequence evolution is found for another star (IRS N1), deeply embedded in an ultra compact HII region, in line with the different evolutionary phases observed in the corresponding HII regions. We have detected the photospheres of OB stars from the more evolved diffuse HII region to the much younger UCHII regions, suggesting that the OB stars have finished their formation and cleared away their possible circumstellar disks very fast. Only in the hyper-compact HII region (IRS5), the early type stars are still surrounded by circumstellar material.

Context. This study is part of a large project to study the physics of accretion and molecular outflows towards a selected sample of high-mass star-forming regions that show evidence of infall and rotation from previous studies. Aims. We wish to make a thorough study at high-angular resolution of the structure and kinematics of the HMCs and corresponding molecular outflows in the high-mass star-forming region G24.78+0.08. Methods. We carried out SMA and IRAM PdBI observations at 1.3 and 1.4 mm, respectively, of dust and of typical high-density and molecular outflow tracers with resolutions of <1". Complementary IRAM 30-m 12CO and 13CO observations were carried out to recover the short spacing information of the molecular outflows. Results. The millimeter continuum emission towards cores G24 A1 and A2 has been resolved into 3 and 2 cores, respectively, and named A1, A1b, A1c, A2, and A2b. All these cores are aligned in a southeast-northwest direction coincident with that of the molecular outflows detected in the region, which suggests a preferential direction for star formation in this region. The masses of the cores range from 7 to 22 Msun, and the rotational temperatures from 128 to 180 K. The high-density tracers have revealed the existence of 2 velocity components towards A1, one of them peaks close to the position of the millimeter continuum peak and of the HC HII region, and is associated with the velocity gradient seen in CH3CN towards this core, while the other one peaks southwest of core A1 and is not associated with any millimeter continuum emission peak. The position-velocity plots along outflow A and the 13CO averaged blueshifted and redshifted emission indicate that this outflow is driven by core A2.

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.

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A_V^tail = 2-4 mag, or equivalently, at \Sigma^tail = 40-80 Ms pc^{-2}. The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n = 10^3 cm^{-3}. This is 5-10 times larger than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A_V^tail may be linked to a transition between this population and the surrounding medium. The star formation rates in the clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A_V^tail, dropping abruptly below that threshold. These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A_V^tail.

The most massive stars can form via standard disk accretion - despite of the radiation pressure generated - due to the fact that the massive accretion disk yields a strong anisotropy in the radiation field, releasing most of the radiation pressure perpendicular to the disk accretion flow. Here, we analyze the self-gravity of the forming circumstellar disk as the potential major driver of the angular momentum transport in such massive disks responsible for the high accretion rates needed for the formation of massive stars. For this purpose, we perform self-gravity radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. The formation and evolution of the resulting circumstellar disk is investigated in 1.) axially symmetric simulations using an alpha-shear-viscosity prescription and 2.) a three-dimensional simulation, in which the angular momentum transport is provided self-consistently by developing gravitational torques in the self-gravitating accretion disk. The simulation series of different strength of the alpha-viscosity shows that the accretion history of the forming star is mostly independent of the alpha-viscosity-parameter. The accretion history of the three-dimensional run driven by self-gravity is more time-dependent than the viscous disk evolution in axial symmetry. The mean accretion rate, i.e. the stellar mass growth, is nearly identical to the alpha-viscosity models. We conclude that the development of gravitational torques in self-gravitating disks around forming massive stars provides a self-consistent mechanism to efficiently transport the angular momentum to outer disk radii. Also the formation of the most massive stars can therefore be understood in the standard accretion disk scenario.

We investigate the radiation pressure feedback in the formation of massive stars in 1, 2, and 3D radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. In contrast to previous research, we consider frequency dependent stellar radiation feedback, resolve the dust sublimation front in the vicinity of the forming star down to 1.27 AU, compute the evolution for several 10^5 yrs covering the whole accretion phase of the forming star, and perform a comprehensive survey of the parameter space. The most fundamental result is that the formation of a massive accretion disk in slowly rotating cores preserves a high anisotropy in the radiation field. The thermal radiation escapes through the optically thin atmosphere, effectively diminishing the radiation pressure feedback onto the accretion flow. Gravitational torques in the self-gravitating disk drive a sufficiently high accretion rate to overcome the residual radiation pressure. Simultaneously, the radiation pressure launches an outflow in the bipolar direction, which grows in angle with time and releases a substantial fraction of the initial core mass from the star-disk system. Summarized, for an initial core mass of 60, 120, 240, and 480 Msol these mechanisms allow the star to grow up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

Many questions remain regarding the properties of disks around massive prototstars. Here we present the observations of a high mass protostellar object including an elongated dust continuum structure perpendicular to the outflow. Submillimeter Array 230 GHz line and continuum observations of the high mass protostellar object IRAS 18151-1208 along with single dish IRAM 30m observations afford us high spatial resolution (0.8") as well as recovery of the extended emission that gets filtered out by the interferometer. The observations of 12CO confirm the outflow direction to be in the southeast-northwest direction, and the 1.3 mm continuum exhibits an elongation in the direction perpendicular to the outflow. We model the physical parameters of the elongated structure by simultaneously fitting the observed spectral energy distribution (SED) and the brightness profile along the major axis using the 3D Radiative Transfer code MC3D. Assuming a density profile similar to that of a low mass disk, we can also reproduce the observations of this high mass protostellar object. This is achieved by using the same density distribution and flaring parameters as were used in the low mass case, and scaling up the size parameters that successfully modeled the circumstellar disk of several T Tauri stars. We also calculate that a region within the inner 30 AU of such a high mass disk is stable under the Toomre criterion. While we do not rule out other scenarios, we show here that the observations in the high mass regime are consistent with a scaled up version of a low mass disk. Implications on high mass star formation are discussed.

Massive stars play an important role in shaping the structure of galaxies. Infrared dark clouds (IRDCs), with their low temperatures and high densities, have been identified as the potential birthplaces of massive stars. In order to understand the formation processes of massive stars the physical and chemical conditions in infrared dark clouds have to be characterized. The goal of this paper is to investigate the chemical composition of a sample of southern infrared dark clouds. One important aspect of the observations is to check, if the molecular abuncances in IRDCs are similar to the low-mass pre-stellar cores, or whether they show signatures of more evolved evolutionary stages. We performed observations toward 15 IRDCs in the frequency range between 86 and 93 GHz using the 22-m Mopra radio telescope. We detect HNC, HCO$^+$ and HNC emission in all clouds and N$_2$H$^+$ in all IRDCs except one. In some clouds we detect SiO emission. Complicated shapes of the HCO$^+$ emission line profile are found in all IRDCs. Both signatures indicates the presence of infall and outflow motions and beginning of star formation activity, at least in some parts of the IRDCs. Where possible, we calculate molecular abundances and make a comparison with previously obtained values for low-mass pre-stellar cores and high-mass protostellar objects (HMPOs). We show a tendency for IRDCs to have molecular abundances similar to low-mass pre-stellar cores rather than to HMPOs abundances on the scale of our single-dish observations.

To understand evolutionary and environmental effects during the formation of high-mass stars, we observed three regions of massive star formation at different evolutionary stages that reside in the same natal molecular cloud. Methods. The three regions S255IR, S255N and S255S were observed at 1.3 mm with the Submillimeter Array (SMA) and followup short spacing information was obtained with the IRAM 30m telescope. Near infrared (NIR) H + K-band spectra and continuum observations were taken for S255IR with VLT-SINFONI to study the different stellar populations in this region. The combination of millimeter (mm) and near infrared data allow us to characterize different stellar populations within the young forming cluster in detail. While we find multiple mm continuum sources toward all regions, their outflow, disk and chemical properties vary considerably. The most evolved source S255IR exhibits a collimated bipolar outflow visible in CO and H2 emission, the outflows from the youngest region S255S are still small and rather confined in the regions of the mm continuum peaks. Also the chemistry toward S255IR is most evolved exhibiting strong emission from complex molecules, while much fewer molecular lines are detected in S255N, and in S255S we detect only CO isotopologues and SO lines. Also, rotational structures are found toward S255N and S255IR. Furthermore, a comparison of the NIR SINFONI and mm data from S255IR clearly reveal two different (proto) stellar populations with an estimated age difference of approximately 1 Myr. A multi-wavelength spectroscopy and mapping study reveals different evolutionary phases of the star formation regions. We propose the triggered outside-in collapse star formation scenario for the bigger picture and the fragmentation scenario for S255IR.

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase of several 10^5 yr for the forming massive star and to perform a comprehensive scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. Our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field, yielding to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msol the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase (up to the end of the accretion disk epoch) for the forming massive star and to perform a broad scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. While confirming the upper mass limit of spherically symmetric accretion, our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field. This yields to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations, far beyond the upper mass limit of spherical accretion. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msun the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msun respectively.

We present multi-wavelengths observations and a radiative transfer model of a newly discovered massive circumstellar disk of gas and dust which is one of the largest disks known today. Seen almost edge-on, the disk is resolved in high-resolution near-infrared (NIR) images and appears as a dark lane of high opacity intersecting a bipolar reflection nebula. Based on molecular line observations we estimate the distance to the object to be 3.5 kpc. This leads to a size for the dark lane of ~10500 AU but due to shadowing effects the true disk size could be smaller. In Spitzer/IRAC 3.6 micron images the elongated shape of the bipolar reflection nebula is still preserved and the bulk of the flux seems to come from disk regions that can be detected due to the slight inclination of the disk. At longer IRAC wavelengths, the flux is mainly coming from the central regions penetrating directly through the dust lane. Interferometric observations of the dust continuum emission at millimeter wavelengths with the SMA confirm this finding as the peak of the unresolved mm-emission coincides perfectly with the peak of the Spitzer/IRAC 5.8 micron flux and the center of the dark lane seen in the NIR images. Simultaneously acquired CO data reveal a molecular outflow along the northern part of the reflection nebula which seems to be the outflow cavity. An elongated gaseous disk component is also detected and shows signs of rotation. The emission is perpendicular to the molecular outflow and thus parallel to but even more extended than the dark lane in the NIR images. Based on the dust continuum and the CO observations we estimate a disk mass of up to a few solar masses depending on the underlying assumptions. Whether the disk-like structure is an actual accretion disk or rather a larger-scale flattened envelope or pseudodisk is difficult to discriminate with the current dataset (abridged).

Aims: Our aim is to understand the evolutionary sequence of high-mass star formation from the earliest evolutionary stage of high-mass starless cores, via high-mass cores with embedded low- to intermediate-mass objects, to finally high-mass protostellar objects. Methods: Herschel far-infrared PACS and SPIRE observations are combined with existing data at longer and shorter wavelengths to characterize the spectral and physical evolution of massive star-forming regions. Results: The new Herschel images spectacularly show the evolution of the youngest and cold high-mass star-forming regions from mid-infrared shadows on the Wien-side of the spectral energy distribution (SED), via structures almost lost in the background emission around 100mum, to strong emission sources at the Rayleigh-Jeans tail. Fits of the SEDs for four exemplary regions covering evolutionary stages from high-mass starless cores to high-mass protostellar objects reveal that the youngest regions can be fitted by single-component black-bodies with temperatures on the order of 17K. More evolved regions show mid-infrared excess emission from an additional warmer component, which however barely contributes to the total luminosities for the youngest regions. Exceptionally low values of the ratio between bolometric and submm luminosity additionally support the youth of the infrared-dark sources. Conclusions: The Herschel observations reveal the spectral and physical properties of young high-mass star-forming regions in detail. The data clearly outline the evolutionary sequence in the images and SEDs. Future work on larger samples as well as incorporating full radiative transfer calculations will characterize the physical nature at the onset of massive star formation in even more depth.

The intermediate-mass star-forming core UYSO 1 has previously been found to exhibit intriguing features. While deeply embedded and previously only identified by means of its (sub-)millimeter emission, it drives two powerful, dynamically young, molecular outflows. Although the process of star formation has obviously started, the chemical composition is still pristine. We present Herschel PACS and SPIRE continuum data of this presumably very young region. The now complete coverage of the spectral energy peak allows us to precisely constrain the elevated temperature of 26 - 28 K for the main bulge of gas associated with UYSO1, which is located at the interface between the hot HII region Sh 2-297 and the cold dark nebula LDN 1657A. Furthermore, the data identify cooler compact far-infrared sources of just a few solar masses, hidden in this neighbouring dark cloud.

Infrared-dark clouds (IRDCs) are the precursors to massive stars and stellar clusters. G011.11-0.12 is a well-studied filamentary IRDC, though, to date, the absence of far-infrared data with sufficient spatial resolution has limited the understanding of the structure and star-formation activity. We use Herschel to study the embedded population of young pre- and protostellar cores in this IRDC. We examine the cloud structure, which appears in absorption at short wavelength and in emission at longer wavelength. We derive the properties of the massive cores from the spectral energy distributions of bright far-infrared point sources detected with the PACS instrument aboard Herschel. We report on the detection and characterization of pre- and protostellar cores in a massive filamentary infrared-dark cloud G011.11-0.12 using PACS. We characterize 18 cores directly associated with the filament, two of which have masses over 50 Msun, making them the best candidates to become massive stars in G011.11-0.12. These cores are likely at various stages of protostar formation, showing elevated temperature (<T> ~ 22 K) with respect to the ambient gas reservoir. The core masses (<M> ~ 24 Msun) are small compared to that in the cold filament. The mean core separation is 0.9 pc, well in excess of the Jeans length in the filament. We confirm that star formation in IRDCs is underway and diverse, and IRDCs have the capability of forming massive stars and clusters.

We present VLT/SINFONI integral field spectroscopy of RCW 34 along with Spitzer/IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class II objects is found. At the northern edge of this bubble, an HII region is located, ionized by 3 OB stars. Intermediate mass stars (2 - 3 Msun) are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. North of the HII region, a photon-dominated region is present, marking the edge of a dense molecular cloud traced by H2 emission. Several class 0/I objects are associated with this cloud, indicating that star formation is still taking place. The distance to RCW 34 is revised to 2.5 +- 0.2 kpc and an age estimate of 2 - 1 Myrs is derived from the properties of the PMS stars inside the HII region. The most likely scenario for the formation of the three regions is that star formation propagates from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, after that, the HII region is formed from a dense core at the edge of the molecular cloud, resulting in the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud which in turn induces the current star-formation in the molecular cloud. (b) An external triggering is responsible for the star-formation propagating from South to North. [abridged]

A 5'x5' region encompassing NGC6334I and I(N) has been mapped at a wavelength of 3mm (from 83.5 to 115.5GHz) with the Mopra telescope at an angular resolution between 33 arcsec and 36 arcsec. This investigation has made use of the recently installed 3mm MMIC receiver and the Mopra Spectrometer (MOPS) with broadband capabilities permitting total coverage of the entire frequency range with just five different observations. In total, the spatial distribution of nineteen different molecules, ions and radicals, along with additional selected isotopologues have been studied. Whilst most species trace the sites of star formation, CH_3CN appears to be most closely associated with NGC6334I and I(N). Both CN and C_2H appear to be widespread, tracing gas that is not associated with active star formation. Both N_2H^+ and HC_3N closely resemble dust continuum emission, showing they are reliable tracers of dense material, as well as the youngest stages of high mass star formation. Hot (E_u/k>100K) thermal CH_3OH emission is preferentially found towards NGC6334I, contrasting with I(N), where only cold (E_u/k<22K) thermal CH_3OH emission is found.

We used the CRIRES spectrograph on the VLT to study the ro-vibrational 12CO/13CO, the Pfund beta and H2 emission between 4.59 and 4.72mu wavelengths toward the BN object, the disk candidate source n, and a proposed dust density enhancement IRC3. We detected CO absorption and emission features toward all three targets. Toward the BN object, the data partly confirm the results obtained more than 25 years ago by Scoville et al., however, we also identify several new features. While the blue-shifted absorption is likely due to outflowing gas, toward the BN object we detect CO in emission extending in diameter to ~3300AU. Although at the observational spectral resolution limit, the 13CO line width of that feature increases with energy levels, consistent with a disk origin. If one attributes the extended CO emission also to a disk origin, its extent is consistent with other massive disk candidates in the literature. For source n, we also find the blue-shifted CO absorption likely from an outflow. However, it also exhibits a narrower range of redshifted CO absorption and adjacent weak CO emission, consistent with infalling motions. We do not spatially resolve the emission for source n. For both sources we conduct a Boltzmann analysis of the 13CO absorption features and find temperatures between 100 and 160K, and H2 column densities of the order a few times 10^23cm^-2. The observational signatures from IRC3 are very different with only weak absorption against a much weaker continuum source. However, the CO emission is extended and shows wedge-like position velocity signatures consistent with jet-entrainment of molecular gas, potentially associated with the Orion-KL outflow system. We also present and discuss the Pfund beta and H2 emission in the region.

Aims: We derive the probability density functions (PDFs) of column density for a complete sample of prominent molecular cloud complexes closer than 200 pc. Methods: We derive near-infrared dust extinction maps for 23 molecular cloud complexes, using the "nicest" colour excess mapping technique and data from the 2MASS archive. The extinction maps are then used to examine the column density PDFs in the clouds. Results: The column density PDFs of most molecular clouds are well-fitted by log-normal functions at low column densities (0.5 mag < A_v < 3-5 mag). However, at higher column densities prominent, power-law-like wings are common. In particular, we identify a trend among the PDFs: active star-forming clouds always have prominent non-log-normal wings. In contrast, clouds without active star formation resemble log-normals over the whole observed column density range, or show only low excess of higher column densities. This trend is also reflected in the cumulative PDFs, showing that the fraction of high column density material is significantly larger in star-forming clouds. These observations are in agreement with an evolutionary trend where turbulent motions are the main cloud-shaping mechanism for quiescent clouds, but the density enhancements induced by them quickly become dominated by gravity (and other mechanisms) which is strongly reflected by the shape of the column density PDFs. The dominant role of the turbulence is restricted to the very early stages of molecular cloud evolution, comparable to the onset of active star formation in the clouds.

We observed with the VLA, PdBI, and SMA the centimeter and millimeter continuum, N2H+(1-0), and CO(2-1) emission associated with a dusty cloud harboring a nascent cluster with intermediate-mass protostars. At centimeter wavelengths we found a strong source, tracing a UCHII region, at the eastern edge of the dusty cloud, with a shell-like structure, and with the near-infrared counterpart falling in the center of the shell. This is presumably the most massive source of the forming cluster. About 15'' to the west of the UCHII region and well embedded in the dusty cloud, we detected a strong millimeter source, MM1, associated with centimeter and near-infrared emission. MM1 seems to be driving a prominent high-velocity CO bipolar outflow, and is embedded in a ridge of dense gas traced by N2H+. We estimated that MM1 is an intermediate-mass source in the Class 0/I phase. About 15'' to the south of MM1, and still more deeply embedded in the dusty cloud, we detected a compact millimeter source, MM2, with neither centimeter nor near-infrared emission, but with water maser emission. MM2 is associated with a clump of N2H+, whose kinematics reveal a clear velocity gradient and additionally we found signposts of infall motions. MM2, being deeply embedded within the dusty cloud, with an associated water maser but no hints of CO outflow emission, is an intriguing object, presumably of intermediate mass. In conclusion, the UCHII region is found at the border of a dusty cloud which is currently undergoing active star formation. Two intermediate-mass protostars in the dusty cloud seem to have formed after the UCHII region and have different properties related to the outflow phenomenon.

We present a high-resolution, multi-wavelength study of the massive protostellar cluster NGC 6334I(N) that combines new spectral line data from the Submillimeter Array (SMA) and VLA with a reanalysis of archival VLA continuum data, 2MASS and Spitzer images. As shown previously, the brightest 1.3 mm source SMA1 contains substructure at subarcsecond resolution, and we report the first detection of SMA1b at 3.6 cm along with a new spatial component at 7 mm (SMA1d). We find SMA1 (aggregate of sources a, b, c, and d) and SMA4 to be comprised of free-free and dust components, while SMA6 shows only dust emission. Our 1.5" resolution 1.3 mm molecular line images reveal substantial hot-core line emission toward SMA1 and to a lesser degree SMA2. We find CH3OH rotation temperatures of 165\pm 9 K and 145\pm 12 K for SMA1 and SMA2, respectively. We estimate a diameter of 1400 AU for the SMA1 hot core emission, encompassing both SMA1b and SMA1d, and speculate that these sources comprise a >800 AU separation binary that may explain the previously-suggested precession of the outflow emanating from the SMA1 region. The LSR velocities of SMA1, SMA2, and SMA4 all differ by 1-2 km/s. Outflow activity from SMA1, SMA2, SMA4, and SMA6 is observed in several molecules including SiO(5--4) and IRAC 4.5 micron emission; 24 micron emission from SMA4 is also detected. Eleven water maser groups are detected, eight of which coincide with SMA1, SMA2, SMA4, and SMA6. We also detect a total of 83 Class I CH3OH 44GHz maser spots which likely result from the combined activity of many outflows. Our observations paint the portrait of multiple young hot cores in a protocluster prior to the stage where its members become visible in the near-infrared.

To better understand the physical properties of accretion disks in high-mass star formation, we present a study of a 12 high-mass accretion disk candidates observed at high spatial resolution with the Australia Telescope Compact Array (ATCA) in the NH3 (4,4) and (5,5) lines. Almost all sources were detected in NH3, directly associated with CH3OH Class II maser emission. From the remaining eleven sources, six show clear signatures of rotation and/or infall motions. These signatures vary from velocity gradients perpendicular to the outflows, to infall signatures in absorption against ultracompact HII regions, to more spherical infall signatures in emission. Although our spatial resolution is ~1000AU, we do not find clear Keplerian signatures in any of the sources. Furthermore, we also do not find flattened structures. In contrast to this, in several of the sources with rotational signatures, the spatial structure is approximately spherical with sizes exceeding 10^4 AU, showing considerable clumpy sub-structure at even smaller scales. This implies that on average typical Keplerian accretion disks -- if they exist as expected -- should be confined to regions usually smaller than 1000AU. It is likely that these disks are fed by the larger-scale rotating envelope structure we observe here. Furthermore, we do detect 1.25cm continuum emission in most fields of view.

To characterize the initial conditions for intermediate- to high-mass star formation, we observed two Infrared Dark Clouds (IRDCs) that remain absorption features up to 70mum wavelength, with the PdBI in the 3.23mm dust continuum as well as the N2H+(1--0) and 13CS(2-1) line emission. While IRDC19175-4 is clearly detected in the 3.23mm continuum, the second source in the field, IRDC19175-5, is only barely observable above the 3sigma continuum detection threshold. However, the N2H+(1-0) observations reveal 17 separate sub-sources in the vicinity of the two IRDCs. Most of them exhibit low levels of turbulence (dv \leq 1km/s), indicating that the fragmentation process in these cores may be dominated by the interplay of thermal pressure and gravity, but not so much by turbulence. Combining the small line widths with the non-detection up to 70mum and the absence of other signs of star formation activity, most of these 17 cores with masses between sub-solar to ~10M_sun are likely still in a starless phase. Furthermore, we find a large CS depletion factor of the order 100. Although the strongest line and continuum peak is close to virial equilibrium, its slightly broader line width compared to the other cores is consistent with it being in a contraction phase potentially at the verge of star formation. The relative peak velocities between neighboring cores are usually below 1km/s, and we do not identify streaming motions along the filamentary structures. Average densities are between 10^5 and 10^6cm^{-3} (one to two orders of magnitude larger than for example in the Pipe nebula) implying relatively small Jeans-lengths that are consistent with the observed core separations of the order 5000AU. The quest for high-mass starless cores prior to any star formation activity remains open.

We examine an Infrared Dark Cloud (IRDC) at high spatial resolution as a means to study rotation, outflow, and infall at the onset of massive star formation. Submillimeter Array observations combined with IRAM 30 meter data in 12CO(2--1) reveal the outflow orientation in the IRDC 18223-3 region, and PdBI 3 mm observations confirm this orientation in other molecular species. The implication of the outflow's presence is that an accretion disk is feeding it, so using high density tracers such as C18O, N2H+, and CH3OH, we looked for indications of a velocity gradient perpendicular to the outflow direction. Surprisingly, this gradient turns out to be most apparent in CH3OH. The large size (28,000 AU) of the flattened rotating object detected indicates that this velocity gradient cannot be due solely to a disk, but rather from inward spiraling gas within which a Keplerian disk likely exists. From the outflow parameters, we derive properties of the source such as an outflow dynamical age of ~37,000 years, outflow mass of ~13 M_sun, and outflow energy of ~1.7 x 10^46 erg. While the outflow mass and energy are clearly consistent with a high-mass star forming region, the outflow dynamical age indicates a slightly more evolved evolutionary stage than previous spectral energy distribution (SED) modeling indicates. The calculated outflow properties reveal that this is truly a massive star in the making. We also present a model of the observed methanol velocity gradient. The rotational signatures can be modeled via rotationally infalling gas. These data present evidence for one of the youngest known outflow/infall/disk systems in massive star formation. A tentative evolutionary picture for massive disks is discussed.

We investigate the structure of accretion disks around massive protostar applying steady state models of thin disks. The thin disk equations are solved with proper opacities for dust and gas taking into account the huge temperature variation along the disk. We explore a wide parameter range concerning stellar mass, accretion rate, and viscosity parameter \alpha . The most essential finding is a very high temperature of the inner disk. For e.g. a 10 M_sun protostar and an accretion rate of 10^-4 M_sun/yr, the disk midplane temperature may reach almost 10^5 K. The disk luminosity in this case is about 10^4 L_sun and, thus, potentially higher than that of a massive protostar. We motivate our disk model with similarly hot disks around compact stars. We calculate a dust sublimation radius by turbulent disk self-heating of more than 10AU, a radius, which is 3 times larger than caused by stellar irradiation. We discuss implications of this result on the flashlight effect and the consequences for the radiation pressure of the central star. In difference to disks around low mass protostars our models suggest rather high values for the disk turbulence parameter \alpha close to unity. However, disk stability to fragmentation due to thermal effects and gravitational instability would require a lower \alpha value. For \alpha = 0.1 we find stable disks out to 80AU. Essentially, our model allows to compare the outer disk to some of the observed massive protostellar disk sources, and from that, extrapolate on the disk structure close to the star which is yet impossible to observe.

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

We combine multifrequency observations from the millimeter to near infrared wavelengths that demonstrate the spatial distributions of H2, CO, and NH3 emission, which are all manifestations of various shocks driven by outflows of deeply embedded sources in NGC6334I. In addition to the well-known northeast-southwest outflow we detect at least one more outflow in the region by combining observations from APEX, ATCA, SMA, Spitzer and VLT/ISAAC. Potential driving sources will be discussed. NGC6334I exhibits several signs of active star formation and will be a major target for future observatories such as Herschel and ALMA.

Massive star formation exhibits an extremely rich chemistry. However, not much evolutionary details are known yet, especially at high spatial resolution. Therefore, we synthesize previously published Submillimeter Array high-spatial-resolution spectral line observations toward four regions of high-mass star formation that are in various evolutionary stages with a range of luminosities. Estimating column densities and comparing the spatially resolved molecular emission allows us to characterize the chemical evolution in more detail. Furthermore, we model the chemical evolution of massive warm molecular cores to be directly compared with the data. The four regions reveal many different characteristics. While some of them, e.g., the detection rate of CH3OH, can be explained by variations of the average gas temperatures, other features are attributed to chemical effects. For example, C34S is observed mainly at the core-edges and not toward their centers because of temperature-selective desorption and successive gas-phase chemistry reactions. Most nitrogen-bearing molecules are only found toward the hot molecular cores and not the earlier evolutionary stages, indicating that the formation and excitation of such complex nitrogen-bearing molecules needs significant heating and time to be fully developed. Furthermore, we discuss the observational difficulties to study massive accretion disks in the young deeply embedded phase of massive star formation. The general potential and limitations of such kind of dataset are discussed, and future directions are outlined. The analysis and modeling of this source sample reveals many interesting features toward a chemical evolutionary sequence. However, it is only an early step, and many observational and theoretical challenges in that field lie ahead.

Although most young massive stars appear to be part of multiple systems, it is poorly understood how this multiplicity influences the formation of massive stars. The high-mass star-forming region W3 IRS5 is a prime example of a young massive cluster where the cluster center is resolved into multiple subsources at cm and infrared wavelengths, a potential proto-Trapezium system. The region W3 IRS5 was mapped with the PdBI at 1.4mm and 3.4mm in the AB configurations, observing shock-tracing SiO and SO_2 emission. In the continuum we detect five sources, one of them for the first time, while counterparts were detected in the NIR, MIR or at radio wavelengths for the remaining four sources. Three of the detected sources are within the inner 2100AU, where the protostellar number density exceeds 10^6 protostars pc^-3 assuming spherical symmetry. Lower limits for the circumstellar masses of the detected sources were calculated, although they were strongly affected by the spatial filtering of the interferometer. However, the projected separations of the sources ranging between ~750 and ~4700AU indicate a multiple, Trapezium-like system. We detected five molecular outflows in SiO, two of them nearly in the line of sight direction, which allowed us to see the collapsing protostars in the NIR through the cavities carved by the outflows. The SO_2 velocity structure indicates a rotating, bound system, and we find tentative signatures of converging flows as predicted by the gravoturbulent star formation and converging flow theories. The obtained data strongly indicate that the clustered environment has a major influence on the formation of high-mass stars; however, our data do not clearly allow us to distinguish whether the ongoing star-forming process follows a monolithic collapse or a competitive accretion mechanism.

Supermassive black holes in giant elliptical galaxies are remarkably faint given their expected accretion rates. This motivates models of radiatively inefficient accretion, due to either ion-electron thermal decoupling, generation of outflows that inhibit accretion, or settling of gas to a gravitationally unstable disk that forms stars in preference to feeding the black hole. The latter model predicts the presence of cold molecular gas in a thin disk around the black hole. Here we report Submillimeter Array observations of the nucleus of the giant elliptical galaxy M87 that probe 230 GHz continuum and CO(J=2--1) line emission. Continuum emission is detected from the nucleus and several knots in the jet, including one that has been undergoing flaring behavior. We estimate a conservative upper limit on the mass of molecular gas within ~100pc and +-400km/s line of sight velocity of the central black hole of ~8x10^6Msun, which includes an allowance for possible systematic errors associated with subtraction of the continuum. Ignoring such errors, we have a 3 sigma sensitivity to about 3x10^6Msun. In fact, the continuum-subtracted spectrum shows weak emission features extending up to 4 sigma above the RMS dispersion of the line-free channels. These may be artifacts of the continuum subtraction process. Alternatively, if they are interpreted as CO emission, then the implied molecular gas mass is ~5x10^6Msun spread out over a velocity range of 700km/s. These constraints on molecular gas mass are close to the predictions of the model of self-gravitating, star-forming accretion disks fed by Bondi accretion (Tan & Blackman 2005).

We present sensitive Very Large Array observations with an angular resolution of a few arcseconds of the $J= 1 - 0$ line of SiO in the $v$=1 and 2 vibrationally excited states toward a sample of 60 Galactic regions in which stars of high or intermediate mass are currently forming and/or have recently formed. We report the detection of SiO maser emission in \textit{both} vibrationally excited transitions toward only three very luminous regions: Orion-KL, W51N and Sgr B2(M). Toward all three, SiO maser emission had previously been reported, in Orion-KL in both lines, in W51N only in the $v=2$ line and in Sgr B2(M) only in the $v=1$ line. Our work confirms that SiO maser emission in star-forming regions is a rare phenomenon, indeed, that requires special, probably extreme, physical and chemical conditions not commonly found. In addition to this SiO maser survey, we also present images of the simultaneously observed 7 mm continuum emission from a subset of our sample of star-forming regions where such emission was detected. This is in most cases likely to be free-free emission from compact- and ultracompact-HII regions.

We present Spitzer Space Telescope IRAC and MIPS observations toward a sample of nine high-mass star forming regions at a distance of around 2 kpc. Based on IRAC and MIPS 24 $\mu$m photometric results and 2MASS JHKs data, we carry out a census of young stellar objects (YSOs) in a 5' by 5' field toward each region. Toward seven out of the nine regions, we detect parsec sized clusters with around 20 YSOs surrounded by a more extended and sparse distribution of young stars and protostars. For the other two regions, IRAS 20126+4104 and IRAS 22172+5549, the former has the lowest number of YSOs in the sample and shows no obvious cluster, and the latter appears to be part of a larger, potentially more evolved cluster. The deep IRAC imaging reveals at least twelve outflows in eight out of the nine regions, with nine outflows prominent in the 4.5 $\mu$m band most probably attributed to shocked H$_2$ emission, two outflows dominated by scattered light in the 3.6 and 4.5 $\mu$m bands, and one outflow standing out from its hydrocarbon emission in the 8.0 $\mu$m band. In comparison with previous ground-based observations, our IRAC observations reveal new outflow structures in five regions. The dramatically different morphologies of detected outflows can be tentatively interpreted in terms of possible evolution of massive outflows. The driving sources of these outflows are deeply embedded in dense dusty cores revealed by previous millimeter interferometric observations. We detect infrared counterparts of these dusty cores in the IRAC or MIPS 24 $\mu$m bands. Reflection nebulae dominated by the emission from UV heated hydrocarbons in the 8 $\mu$m band can be found in most regions and they may imply the presence of young B stars.

The enigmatic outflows of the Orion-KL region have raised discussions about their potential driving sources for several decades. Here, we present C18O(2-1) observations combined from the Submillimeter Array and the IRAM30m telescope. The molecular gas is associated on large scales with the famous northwest-southeast high-velocity outflow whereas the high-velocity gas on small spatial scales traces back to the recently identified submm source SMA1. Therefore, we infer that SMA1 may host the driving source of this outflow. Based on the previously published thermal and maser SiO data, source I is the prime candidate to drive the northeast-southwest low-velocity outflow. The source SMA1 is peculiar because it is only detected in several submm wavelength bands but neither in the infrared nor cm regime. We discuss that it may be a very young intermediate- to high-mass protostar. The estimated outflow masses are high whereas the dynamical time-scale of the outflow is short of the order 10^3yrs.

Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy.

This paper presents a summary of a panel discussion to more directly confront theoretical models and observational constraints in massive star formation research. The panel was moderated by Hans Zinnecker, and panel members included Ian Bonnell, Chris McKee, Francesco Palla, Malcolm Walmsley, and Harold Yorke. Additional ample discussion with the audience is recorded.

APEX single-dish observations at sub-millimeter wavelengths toward a sample of massive star-forming regions reveal that C2H is almost omni-present toward all covered evolutionary stages from Infrared Dark Clouds via High-Mass Protostellar Objects to Ultracompact HII regions. High-resolution data from the Submillimeter Array toward one hot-core like High-Mass Protostellar Object show a shell-like distribution of C2H with a radius of ~9000AU around the central submm peak position. These observed features are well reproduced by a 1D cloud model with power-law density and temperature distributions and a gas-grain chemical network. The reactive C2H radical (ethynyl) is abundant from the onset of massive star formation, but later it is rapidly transformed to other molecules in the core center. In the outer cloud regions the abundance of C2H remains high due to constant replenishment of elemental carbon from CO being dissociated by the interstellar UV photons. We suggest that C2H may be a molecule well suited to study the initial conditions of massive star formation.

Aims: Investigation of the dense gas, the outflows and the continuum emission from the massive twin cores NGC6334I and I(N) at high spatial resolution. Methods: We imaged the region with the Australia Telescope Compact Array (ATCA) at 3.4mm wavelength in continuum as well as CH3CN(5_K-4_K) and HCN(1-0) spectral line emission. Results: While the continuum emission in NGC6334I mainly traces the UCHII region, toward NGC6334I(N) we detect line emission from four of the previously identified dust continuum condensations that are of protostellar or pre-stellar nature. The CH3CN(5_K-4_K) lines are detected in all K-components up to energies of 128K above ground toward two protostellar condensations in both regions. We find line-width increasing with increasing K for all sources, which indicates a higher degree of internal motions closer to the central protostars. Toward the main mm and CH3CN source in NGC6334I we identify a velocity gradient approximately perpendicular to the large-scale molecular outflow. This may be interpreted as a signature of an accretion disk, although other scenarios, e.g., an unresolved double source, could produce a similar signature as well. No comparable signature is found toward any of the other sources. HCN does not trace the dense gas well but it is dominated by the molecular outflows. While the outflow in NGC6334I exhibits a normal Hubble-law like velocity structure, the data indicate a precessing outflow close to the plane of the sky for NGC6334I(N). Furthermore, we observe a wide (~15.4km/s) HCN absorption line, much broader than the previously observed CH3OH and NH3 absorption lines. Several explanations for the difference are discussed.

We present a spatially resolved 894 micron map of the circumstellar disk of the Butterfly star in Taurus (IRAS 04302+2247), obtained with the Submillimeter Array (SMA). The predicted and observed radial brightness profile agree well in the outer disk region, but differ in the inner region with an outer radius of ~80-120 AU. In particular, we find a local minimum of the radial brightness distribution at the center, which can be explained by an increasing density / optical depth combined with the decreasing vertical extent of the disk towards the center. Our finding indicates that young circumstellar disks can be optically thick at wavelengths as long as 894 micron. While earlier modeling lead to general conclusions about the global disk structure and, most importantly, evidence for grain growth in the disk (Wolf, Padgett, & Stapelfeldt 2003), the presented SMA observations provide more detailed constraints for the disk structure and dust grain properties in the inner, potentially planet-forming region (inside ~80-120 AU) vs. the outer disk region (~120-300 AU).

This article presents recent work to constrain the physical and chemical properties in high-mass star formation based largely on interferometric high-spatial-resolution continuum and spectral line studies at (sub)mm wavelengths. After outlining the concepts, potential observational tests, a proposed evolutionary sequence and different possible definitions for massive protostars, four particular topics are highlighted: (a) What are the physical conditions at the onset of massive star formation? (b) What are the characteristics of potential massive accretion disks and what do they tell us about massive star formation in general? (c) How do massive clumps fragment, and what does it imply to high-mass star formation? (d) What do we learn from imaging spectral line surveys with respect to the chemistry itself as well as for utilizing molecules as tools for astrophysical investigations?

Characterizing rotation, infall and accretion disks around high-mass protostars is an important topic in massive star formation research. With the Australia Telescope Compact Array and the Very Large Array we studied a massive disk candidate at high angular resolution in ammonia (NH3(4,4) & (5,5)) tracing the warm disk but not the envelope. The observations resolved at ~0.4'' resolution (corresponding to ~1400AU) a velocity gradient indicative of rotation perpendicular to the molecular outflow. Assuming a Keplerian accretion disk, the estimated protostar-disk mass would be high, similar to the protostellar mass. Furthermore, the position-velocity diagram exhibits additional deviation from a Keplerian rotation profile which may be caused by infalling gas and/or a self-gravitating disk. Moreover, a large fraction of the rotating gas is at temperatures >100K, markedly different to typical low-mass accretion disks. In addition, we resolve a central double-lobe cm continuum structure perpendicular to the rotation. We identify this with an ionized, optically thick jet.

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.

We present arcsecond-scale mid-ir photometry (in the 10.5 micron N band and at 24.8 microns), and low resolution spectra in the N band (R~100) of a candidate high mass protostellar object (HMPO) in IRAS 18151-1208 and of two HMPO candidates in IRAS 20343+4129, IRS 1 and IRS 3. In addition we present high resolution mid-ir spectra (R~80000) of the two HMPO candidates in IRAS 20343+4129. These data are fitted with simple models to estimate the masses of gas and dust associated with the mid-ir emitting clumps, the column densities of overlying absorbing dust and gas, the luminosities of the HMPO candidates, and the likely spectral type of the HMPO candidate for which [Ne II] 12.8 micron emission was detected (IRAS 20343+4129 IRS 3). We suggest that IRAS 18151-1208 is a pre-ultracompact HII region HMPO, IRAS 20343+4129 IRS 1 is an embedded young stellar object with the luminosity of a B3 star, and IRAS 20343+4129 IRS 3 is a B2 ZAMS star that has formed an ultracompact HII region and disrupted its natal envelope.

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.

In order to elucidate the nature of the brightest infrared sources associated with IRAS 20343+4129, IRS1 and IRS3, we observed with the Submillimeter Array (SMA) the 1.3 mm continuum and CO(2-1) emission of the region. Faint millimeter dust continuum emission was detected toward IRS1, and we derived an associated gas mass of ~0.8 Msun. The IRS1 spectral energy distribution agrees with IRS1 being an intermediate-mass Class I source of about 1000 Lsun, whose circumstellar material is producing the observed large infrared excess. We have discovered a high-velocity CO bipolar outflow in the east-west direction, which is clearly associated with IRS1, and the outflow parameters are similar to those of intermediate-mass young stellar objects. Associated with the blue large scale CO outflow lobe, detected with single-dish observations, we only found two elongated low-velocity structures on either side of IRS3. The large-scale outflow lobe is almost completely resolved out by the SMA. Our detected low-velocity CO structures are coincident with elongated H2 emission features. The strongest millimeter continuum condensations in the region are found on either side of IRS3, where the infrared emission is extremely weak, and the CO and H2 elongated structures follow the border of the millimeter continuum emission that is facing IRS3. All these results suggest that the dust is associated with the walls of an expanding cavity driven by IRS3, estimated to be a B2 star. Within and beyond the expanding cavity, the millimeter continuum sources can be sites of future low-mass star formation.

While most sources above 10^5Lsun have already formed an Ultracompact HII region (UCHII), this is not necessarily the case for sources of lower luminosity. Characterizing sources in the transition phase, i.e., very luminous objects without any detectable free-free emission, is important for a general understanding of massive star formation. Therefore, we observed the luminous High-Mass Protostellar Object IRAS23151+5912 with the Submillimeter Array at 875mum in the submm continuum and spectral line emission at sub-arcsecond resolution. The 875mum submm continuum emission has been resolved into at least two condensations. The previously believed driving source of one of the outflows, the infrared source IRS1, is ~0.9'' offset from the main submm peak. Over the entire 4GHz bandwidth we detect an intermediate dense spectral line forest with 27 lines from 8 different species, isotopologues or vibrationally-torsionally excited states. Temperature estimates based on the CH3OH line series result in values of T(Peak1)~150+-50K and T(Peak2)~80~30K for the two submm peak positions, respectively. The SiO(8-7) red- and blue-shifted line maps indicate the presence of two molecular outflows. In contrast, the vibrationally-torsionally excited CH3OH line exhibits a velocity gradient approximately perpendicular to one of the outflows. With a size of approximately 5000AU and no Keplerian rotation signature, this structure does not resemble a genuine accretion disk but rather a larger-scale rotating toroid that may harbor a more common accretion disk at its so far unresolved center.

Infrared Dark Clouds are expected to harbor sources in different, very young evolutionary stages. To better characterize these differences, we observed a sample of 43 massive Infrared Dark Clouds, originally selected as candidate high-mass starless cores, with the IRAM 30m telescope covering spectral line tracers of low-density gas, high-density gas, molecular outflows/jets and temperature effects. The SiO(2-1) observations reveal detections toward 18 sources. Assuming that SiO is exclusively produced by sputtering from dust grains, this implies that at least in 40% of this sample star formation is on-going. A broad range of SiO line-widths is observed (between 2.2 and 65km/s), and we discuss potential origins for this velocity spread. While the low-density tracers 12CO(2-1) and 13CO(1-0) are detected in several velocity components, the high-density tracer H13CO+(1--0) generally shows only a single velocity component and is hence well suited for kinematic distance estimates of IRDCs. Furthermore, the H13CO+ line-width is on average 1.5 times larger than that of previously observed NH3(1,1). This is indicative of more motion at the denser core centers, either due to turbulence or beginning star formation activity. In addition, we detect CH3CN toward only six sources whereas CH3OH is observed toward approximately 40% of the sample. Estimates of the CH3CN and CH3OH abundances are low with average values of 1.2x10^{-10} and 4.3x10^{-10}, respectively. These results are consistent with chemical models at the earliest evolutionary stages of high-mass star formation. Furthermore, the CH3OH abundances compare well to recently reported values for low-mass starless cores.

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.

Aiming at a better understand of the physical and chemical processes in the hot molecular core stage of high-mass star formation, we observed the prototypical hot core G29.96-0.02 in the 862mu band with the Submillimeter Array (SMA) at sub-arcsecond spatial resolution. The observations resolved the hot molecular core into six submm continuum sources with the finest spatial resolution of 0.36''x0.25'' (~1800AU) achieved so far. Four of them located within 7800(AU)^2 comprise a proto-Trapezium system with estimated protostellar densities of 1.4x0^5 protostars/pc^3. The plethora of ~80 spectral lines allows us to study the molecular outflow(s), the core kinematics, the temperature structure of the region as well as chemical effects. The derived hot core temperatures are of the order 300K. We find interesting chemical spatial differentiations, e.g., C34S is deficient toward the hot core and is enhanced at the UCHII/hot core interface, which may be explained by temperature sensitive desorption from grains and following gas phase chemistry. The SiO(8-7) emission outlines likely two molecular outflows emanating from this hot core region. Emission from most other molecules peaks centrally on the hot core and is not dominated by any individual submm peak. Potential reasons for that are discussed. A few spectral lines that are associated with the main submm continuum source, show a velocity gradient perpendicular to the large-scale outflow. Since this velocity structure comprises three of the central protostellar sources, this is not a Keplerian disk. While the data are consistent with a gas core that may rotate and/or collapse, we cannot exclude the outflow(s) and/or nearby expanding UCHII region as possible alternative causes of this velocity pattern.

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.

Aims: The massive twin cores NGC6334I and I(N) are in different evolutionary stages and hence ideal targets to study evolutionary variations within the same larger-scale environment. Here, we study the warm, compact gas components. Methods: We imaged the two regions with the Australia Telescope Compact Array (ATCA) at high angular resolution in the NH3(3,3) to (6,6) inversion lines. Results: Compact emission is detected toward both regions in all observed inversion lines with energy levels up to 407K above ground. This is particularly surprising for NGC6334I(N) since it lacks bright infrared emission and is considered a massive cold core at an early evolutionary stage. High optical depth and multiply-peaked line profiles complicate rotation temperature estimates, and we can only conclude that gas components with temperatures >100K are present in both regions. Toward NGC6334I, we confirm previous reports of NH3(3,3) maser emission toward the outflow bow-shocks. Furthermore, we report the first detection of an NH3(6,6) maser toward the central region of NGC6334I. This maser is centered on the second millimeter (mm) peak and elongated along the outflow axis, indicating that this mm continuum core harbors the driving source of the molecular outflow. Toward the main mm peak in NGC6334I(N), we detect a double-horn line profile in the NH3(6,6) transition. The current data do not allow us to differentiate whether this double-horn profile is produced by multiple gas components along the line of sight, or whether it may trace a potential underlying massive accretion disk. The data to Figures 3 to 7 are also available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/.

We aim at studying the cluster environment surrounding the UCHII region in IRAS 20293+3952, a region in the first stages of formation of a cluster around a high-mass star. BIMA and VLA were used to observe the 3 mm continuum, N2H+(1-0), NH3(1,1) and (2,2), and CH3OH(2-1) emission of the surroundings of the UCHII region. We studied the kinematics of the region and computed the rotational temperature and column density maps by fitting the hyperfine structure of N2H+ and NH3. The dense gas traced by N2H+ and NH3 shows two different clouds, a main cloud to the east of the UCHII region, of 0.5 pc and 250 Msun, and a western cloud, of 0.15 pc and 30 Msun. The dust emission reveals two strong components in the northern side of the main cloud, associated with Young Stellar Objects (YSOs) driving molecular outflows, and two fainter components in the southern side, with no signs of star forming activity. Regarding the CH3OH, we found strong emission in a fork-like structure associated with outflow B, as well as emission associated with outflow A. The rotational temperature is higher in the northern side of the main cloud, around 22 K, where there are most of the YSOs, than in the southern side, around 16 K. There is strong chemical differentiation in the region, since we determined low values of the NH3/N2H+ ratio, around 50, associated with YSOs in the north of the main cloud, and high values, up to 300, associated with cores with no detected YSOs, in the south of the main cloud. Such a chemical differentiation is likely due to abundance/depletion effects. Finally, interaction between the different sources in the region is important, mainly through their UV radiation and outflows (abridged).

HH 212 is a nearby (460 pc) protostellar jet discovered in H$_2$ powered by a Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. It is highly collimated and symmetric with matched pairs of bow shocks on either side of the source. We have mapped it in 850 $\mu$m continuum, SiO ($J=8-7$), CO ($J=3-2$), SO ($N_J=8_9-7_8$), HCO$^+$ ($J = 4-3$), and H$^{13}$CO$^+$ ($J = 4-3$) emission simultaneously at $\sim$ 1$''$ resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum around the source, mainly arising from an inner envelope (i.e., the inner part of a previously seen flattened envelope) and a possible disk. The inner envelope is also seen with rotation in CO, HCO$^+$, and probably SO. Like H$_2$ emission, CO and SiO emission are seen along the jet axis but extending closer to the source, tracing the bow shocks with a broad range of velocities and the continuous structures in between. SO emission is seen only around the source, forming a jetlike structure extending along the jet axis from the source, likely tracing the jet near the launching region. The jet is episodic and bending. It may also be slightly precessing as the jetlike SO structure shows a slight S-shaped symmetry about the source. A hint of jet rotation is also seen across the jet axis. Internal outflow shells are seen in CO and HCO$^+$, associated with the bow shocks in the inner part of the jet. The bases of the HCO$^+$ shells are seen with a hint of rotation similar to that seen in the inner envelope, probably consisted mainly of the material extended from the inner envelope and even the possible disk.

At the onset of high-mass star formation, accreting protostars are deeply embedded in massive cores made of gas and dust. Their spectral energy distribution is still dominated by the cold dust and rises steeply from near-to far-infrared wavelengths. The young massive star-forming region IRDC18223-3 is a prototypical Infrared-Dark-Cloud with a compact mm continuum core that shows no protostellar emission below 8mum. However, based on outflow tracers, early star formation activity was previously inferred for this region. Here, we present recent Spitzer observations from the MIPSGAL survey that identify the central protostellar object for the first time at 24 and 70mum. Combining the mid- to far-infrared data with previous mm continuum observations and the upper limits below 8mum, one can infer physical properties of the central source. At least two components with constant gas mass M and dust temperature T are necessary: one cold component (~15K and ~576M_sun) that contains most of the mass and luminosity, and one warmer component (>=51K and >=0.01M_sun) to explain the 24mum data. The integrated luminosity of ~177L_sun can be used to constrain additional parameters of the embedded protostar from the turbulent core accretion model for massive star formation. The data of IRDC18223-3 are consistent with a massive gas core harboring a low-mass protostellar seed of still less than half a solar mass with high accretion rates of the order 10^-4M_sun/yr. In the framework of this model, the embedded protostar is destined to become a massive star at the end of its formation processes.

We present studies of a massive protocluster AFGL5142 in the J=2-1 transition of the CO isotopologues, SO, CH_3OH and CH_3CN lines, as well as continuum at 225 GHz and 8.4 GHz. The 225 GHz continuum emission reveals three prominent peaks MM-1, MM-2 and MM-3. MM-1 and MM-2 are associated with strong CH_3CN emission with temperatures of 90 \pm 20 and 250 \pm 40 K, respectively, while both MM-1 and MM-3 are associated with faint continuum emission at 8.4 GHz. Additional dust continuum peaks MM-4 and MM-5 appear to be associated with \h2O masers. With many continuum sources at cm and mm wavelengths, and those already identified in the infrared, this region is forming a cluster of stars. The CO and SO emission reveals at least three molecular outflows originating from the center of the dust core. The outflows are well collimated, with terminal velocities up to 50 \kms-1 from the cloud velocity. Outflow A coincides with the SiO jet identified previously by Hunter et al. (1999). Since jet-like outflows and disk-mediated accretion process are physically connected, the well collimated outflows indicate that even in this cluster environment, accretion is responsible for the formation of individual stars in the cluster.

This contribution gives an overview about various applications of (sub)mm interferometry in star formation research. The topics covered are molecular outflows, accretion disks, fragmentation and chemical properties of low- and high-mass star-forming regions. A short outlook on the capabilities of ALMA is given as well.

Young massive star-forming regions are known to produce hot molecular gas cores (HMCs) with a rich chemistry. While this chemistry is interesting in itself, it also allows to investigate important physical parameters. I will present recent results obtained with high-angular-resolution interferometers disentangling the small-scale structure and complexity of various molecular gas components. Early attempts to develop a chemical evolutionary sequence are discussed. Furthermore, I will outline the difficulty to isolate the right molecular lines capable to unambiguously trace potential massive accretion disks.

We present high angular resolution observations toward two massive star forming regions IRAS 18264-1152 and IRAS 23151+5912 with the Plateau de Bure Interferometer (PdBI) in the SiO (J=2-1) and H$^{13}$CO$^{+}$ (J=1-0) lines and at 1.3 mm and 3.4 mm continuum, and with the Very Large Array (VLA) in the NH$_3$ (J,K)=(1,1), (2,2) lines. The NH$_3$ (1,1) and (2,2) emission is detected toward IRAS 18264-1152 only. For IRAS 18264-1152, the SiO observations reveal at least two quasi-perpendicular outflows with high collimation factors, and the most dominant feature is a redshifted jet-like outflow with very high velocities up to about ${\Delta}v=60$ kms$^{-1}$ with respect to the systemic velocity. The very-high-velocity component (${\Delta}v=22 - 60$ kms$^{-1}$) of this outflow is spatially offset from its high-velocity (${\Delta}v=3 - 21$ kms$^{-1}$) component. The SiO line profiles and position-velocity characteristics of these two components suggest that this outflow can be driven by an underlying precessing jet. For IRAS 23151+5912, the bipolar but mainly blueshifted SiO outflow coincides with the inner parts of the single-dish CO outflow. In particular the quasi-parabolic shape of the blueshifted outflow coincides with the near-infrared nebulosity and is consistent with entrainment of the gas by an underlying wide-angle wind. The analysis of the molecular outflow data of the two luminous sources further support high-mass stars forming via a disk-mediated accretion process similar to low-mass stars.

Using the Submillimeter Array (SMA), we have imaged the 1.3 millimeter continuum emission at the center of the massive star-forming regions NGC 6334 I and I(N). In both regions, the SMA observations resolve the emission into multiple millimeter sources, with most of the sources clustered into areas only 10,000 AU in diameter. Toward NGC 6334 I, we find four compact sources: the two brightest (I-SMA1 and I-SMA2) are associated with previously-known ammonia cores; I-SMA3 is coincident with the peak of the compact HII region (NGC 6334 F); and I-SMA4 is a newly-discovered object. While I-SMA3 exhibits a mixture of free-free and dust emission, the rest of the objects are dust cores. Toward NGC 6334 I(N), seven compact dust cores are found, one of which is associated with a faint centimeter source. With the exception of I-SMA3, none of the millimeter sources have infrared counterparts in Spitzer Space Telescope 3-8 micron images. Using a simple physical model for the dust continuum emission, the estimated mass of interstellar material toward each of these compact objects is in the range of 3 to 66 solar masses. The total mass in the compact objects appears to be similar in I and I(N). The small size of these groups of sources suggest that these objects are proto-Trapezia forming in the centers of clusters of low to intermediate-mass stars.

Submillimeter Array 1.3mm line and continuum observations toward the young massive star-forming region IRAS18182-1433 are presented. The data are complemented with short-spacing CO(2-1) observations and SiO(1-0) data from the VLA. Multiple massive outflows emanate from the mm continuum peak. The CO(2-1) data reveal a quadrupolar outflow system consisting of two outflows inclined by \~90 degrees. One outflow exhibits a cone-like red-shifted morphology with a jet-like blue-shifted counterpart where a blue counter-cone can only be tentatively identified. The SiO(1-0) data suggest the presence of a third outflow. Analyzing the 12CO/13CO line ratios indicates decreasing CO line opacities with increasing velocities. The other seven detected molecular species - also high-density tracers like CH3CN, CH3OH, HCOOCH3 - are all ~1-2'' offset from the mm continuum peak, but spatially associated with a strong molecular outflow peak and a cm emission feature indicative of a thermal jet. This spatial displacement between the molecular lines and the mm continuum emission could be either due to an unresolved sub-source at the position of the cm feature, or the outflow/jet itself alters the chemistry of the core enhancing the molecular abundances toward that region. A temperature estimate based on the CH3CN(12_k-11_k) lines suggests temperatures of the order 150K. A velocity analysis of the high-density tracing molecules reveals that at the given spatial resolution none of them shows any coherent velocity structure which would be consistent with a rotating disk. We discuss this lack of rotation signatures and attribute it to intrinsic difficulties to observationally isolate massive accretion disks from the surrounding dense gas envelopes and the molecular outflows.

We review the known properties of molecular outflows from low- and high-mass young stars. General trends among outflows are identified, and the most recent studies on the morphology, kinematics, energetics, and evolution of molecular outflows are discussed, focusing on results from high-resolution millimeter observations. We review the existing four broad classes of outflow models and compare numerical simulations with the observational data. A single class of models cannot explain the range of morphological and kinematic properties that are observed, and we propose a possible solution. The impact of outflows on their cloud is examined, and we review how outflows can disrupt their surrounding environment, through the clearing of gas and the injection of momentum and energy onto the gas at distances from their powering sources from about 0.01 to a few pc. We also discuss the effects of shock-induced chemical processes on the ambient medium, and how these processes may act as a chemical clock to date outflows. Lastly, future outflow research with existing and planned millimeter and submillimeter instruments is presented.

Massive stars have a profound influence on the Universe, but their formation remains poorly understood. We review the current status of observational and theoretical research in this field, describing the various stages of an evolutionary sequence that begins with cold, massive gas cores and ends with the dispersal and ionization of gas by the newly-formed star. The physical processes in massive star formation are described and related to their observational manifestations. Feedback processes and the relation of massive stars to star cluster formation are also discussed. We identify key observational and theoretical questions that future studies should address.

We observed the HH 211 jet in the submillimeter continuum and the CO(3-2) and SiO(8-7) transitions with the Submillimeter Array. The continuum source detected at the center of the outflow shows an elongated morphology, perpendicular to the direction of the outflow axis. The high-velocity emission of both molecules shows a knotty and highly collimated structure. The SiO(8-7) emission at the base of the outflow, close to the driving source, spans a wide range of velocities, from -20 up to 40 km s^{-1}. This suggests that a wide-angle wind may be the driving mechanism of the HH 211 outflow. For distances greater than 5" (1500 AU) from the driving source, emission from both transitions follows a Hubble-law behavior, with SiO(8-7) reaching higher velocities than CO(3-2), and being located upstream of the CO(3-2) knots. This indicates that the SiO(8-7) emission is likely tracing entrained gas very close to the primary jet, while the CO(3-2) is tracing less dense entrained gas. From the SiO(5-4) data of Hirano et al. we find that the SiO(8-7)/SiO(5-4) brightness temperature ratio along the jet decreases for knots far from the driving source. This is consistent with the density decreasing along the jet, from (3-10)x10^6 cm^{-3} at 500 AU to (0.8-4)x10^6 cm^{-3} at 5000 AU from the driving source.

We present bispectrum speckle interferometry of the massive protostellar object IRAS 23151+5912 in the near-infrared K' band. The reconstructed image shows the diffuse nebulosity north-east of two point-like sources in unprecedented detail. The comparison of our near-infrared image with mm continuum and CO molecular line maps shows that the brighter of the two point sources lies near the center of the mm peak, indicating that it is a high-mass protostar. The nebulosity coincides with the blue-shifted molecular outflow component. The most prominent feature in the nebulosity is a bow-shock-like arc. We assume that this feature is associated with a precessing jet which has created an inward-pointed cone in the swept-up material. We present numerical jet simulations that reproduce this and several other features observed in our speckle image of the nebulosity. Our data also reveal a linear structure connecting the central point source to the extended diffuse nebulosity. This feature may represent the innermost part of a jet that drives the strong molecular outflow (PA ~80 degr) from IRAS 23151+5912. With the aid of radiative transfer calculations, we demonstrate that, in general, the observed inner structures of the circumstellar material surrounding high-mass stars are strongly influenced by the orientation and symmetry of the bipolar cavity.

HH 212 is a highly collimated jet discovered in H2 powered by a young Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. We have mapped around it in 1.33 mm continuum, 12CO ($J=2-1$), 13CO ($J=2-1$), C18O ($J=2-1$), and SO ($J_K = 6_5-5_4$) emission at $\sim$ \arcs{2.5} resolution with the Submillimeter Array. A dust core is seen in the continuum around the source. A flattened envelope is seen in C18O around the source in the equator perpendicular to the jet axis, with its inner part seen in 13CO. The structure and kinematics of the envelope can be roughly reproduced by a simple edge-on disk model with both infall and rotation. In this model, the density of the disk is assumed to have a power-law index of $p=-1.5$ or -2, as found in other low-mass envelopes. The envelope seems dynamically infalling toward the source with slow rotation because the kinematics is found to be roughly consistent with a free fall toward the source plus a rotation of a constant specific angular momentum. A 12CO outflow is seen surrounding the H2 jet, with a narrow waist around the source. Jetlike structures are also seen in 12CO near the source aligned with the H2 jet at high velocities. The morphological relationship between the H2 jet and the 12CO outflow, and the kinematics of the 12CO outflow along the jet axis are both consistent with those seen in a jet-driven bow shock model. SO emission is seen around the source and the H2 knotty shocks in the south, tracing shocked emission around them.

We report the first detection of an extragalactic millimeter wavelength H2O maser at 183 GHz towards NGC 3079 using the Submillimeter Array (SMA), and a tentative submillimeter wave detection of the 439 GHz maser towards the same source using the James Clerk Maxwell Telescope (JCMT). These H2O transitions are known to exhibit maser emission in star-forming regions and evolved stars. NGC 3079 is a well-studied nuclear H2O maser source at 22 GHz with a time-variable peak flux density in the range 3 -- 12 Jy. The 183 GHz H2O maser emission, with peak flux density $\sim$0.5 Jy (7$\sigma$ detection), also originates from the nuclear region of NGC 3079 and is spatially coincident with the dust continuum peak at 193 GHz (53 mJy integrated). Peak emission at both 183 and 439 GHz occurs in the same range of velocity as that covered by the 22 GHz spectrum. We estimate the gas to dust ratio of the nucleus of NGC 3079 to be $\approx$150, comparable to the Galactic value of 160. Discovery of maser emission in an active galactic nucleus beyond the long-known 22 GHz transition opens the possibility of future position-resolved radiative transfer modeling of accretion disks and outflows $<1$ pc from massive black holes.

We present 7 mm, 1.3 cm and 3.6 cm continuum observations made with the Very Large Array toward a sample of ten luminous IRAS sources that are believed to be regions of massive star formation. We detect compact 7 mm emission in four of these objects: IRAS 18089-1732(1), IRAS 18182-1433, IRAS 18264-1152 and IRAS 18308-0841 and for the first time find that these IRAS sources are associated with double or triple radio sources separated by a few arcseconds. We discuss the characteristics of these sources based mostly on their spectral indices and find that their nature is diverse. Some features indicate that the 7 mm emission is dominated by dust from disks or envelopes. Toward other components the 7 mm emission appears to be dominated by free-free radiation, both from ionized outflows or from optically thick H II regions. Furthermore, there is evidence of synchrotron contamination in some of these sources. Finally, we found that the sources associated with ionized outflows, or thermal jets are correlated with CH$_3$OH masers. The precise determination of the nature of these objects requires additional multifrequency observations at high angular resolution. The 3.6 cm continuum observations also revealed seven UCHII regions in the vicinity of the sources IRAS 18089-1732(1) and two more in the source IRAS 18182-1433. We show that the small photoionized nebulae of these UCHII regions are produced by early B-type stars.

Combining mid-infrared data from the SPITZER Space Telescope with cold gas and dust emission observations from the Plateau de Bure Interferometer, we characterize the Infrared Dark Cloud IRDC18223-3 at high spatial resolution. The millimeter continuum data reveal a massive `184Msun gas core with a projected size of ~28000AU that has no associated protostellar mid-infrared counterpart. However, the detection of 4.5mum emission at the edge of the core indicates early outflow activity, which is supported by broad CO and CS spectral line-wing emission. Moreover, systematically increasing N2H+(1-0) line-width toward the mm core center can be interpreted as additional evidence for early star formation. Furthermore, the N2H+(1-0) line emission reveals a less massive secondary core which could be in an evolutionary stage prior to any star formation activity.

Submillimeter Array observations of Orion-KL at 1'' resolution in the 440mu/690GHz band reveal new insights about the continuum and line emission of the region. The 440mu continuum flux density measurement from source I allows us to differentiate among the various proposed physical models: Source I can be well modeled by a ``normal'' protostellar SED consisting of a proton-electron free-free emission component at low frequencies and a strong dust component in the submillimeter bands. Furthermore, we find that the protostellar object SMA1 is clearly distinct from the hot core. The non-detection of SMA1 at cm and infrared wavelengths suggests that it may be one of the youngest sources in the entire Orion-KL region. The molecular line maps show emission mainly from the sources I, SMA1 and the hot core peak position. An analysis of the CH$3CN(37_K-36_K) K-ladder (K=0...3) indicates a warm gas component of the order 600+-200K. In addition, we detect a large fraction (~58%) of unidentified lines and discuss the difficulties of line identifications at these frequencies.

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.

We present the first submm (865mum) imaging spectral line survey at one arcsecond resolution conducted with the Submillimeter Array toward Orion-KL. Within the 2x2 GHz bandpasses (lower and upper sidebands, 337.2-339.2GHz and 347.2-349.2GHz), we find about 145 spectral lines from 13 species, 6 isotopologues, and 5 vibrational excited states. Most nitrogen-bearing molecules are strong toward the hot core, whereas the oxygen-bearing molecules peak toward the south-west in the so-called compact ridge. Imaging of spectral lines is shown to be an additional tool to improve the identifications of molecular lines. Arcsecond spatial resolution allows us to distinguish the molecular line emission of the sources I and n from that of the hot core. The only molecular species detected strongly toward source I is SiO, delineating mainly the collimated north-east south-west low-velocity outflow. The two positions close to source I, which have previously been reported to show maser emission in the v=0 28SiO(1-0) and (2-1) lines, show no detectable maser emission in the v=0 28SiO(8-7) line at our spatial resolution. SiO is weak toward source n, and thus source n may not currently be driving a molecular outflow. CH$_3$OH is the molecule with the highest number of identified lines (46) in this spectral window. This ``line forest'' allows us to estimate temperatures in the region, and we find temperatures between 50 and 350K, with the peak temperatures occurring toward the hot core. The detection of strong vibrational excited line emission from the submm continuum peak SMA1 supports the interpretation that the source SMA1 is likely of protostellar nature.

Investigating in more detail the previously suggested massive disk scenario for the High-Mass Protostellar Object IRAS18089-1732, we observed the source in the 860mum band with the Submillimeter Array in various spectral lines and the submm continuum emission at (sub-)arcsecond spatial resolution. Fifty spectral lines from eighteen different species spanning upper level energy states between 17 and 747K were detected. One of the assumed best tracers for massive disks, CH3CN, is optically thick and does not allow a further disk investigation. However, the complex molecule HCOOCH3 appears optically thin and exhibits a velocity shift across the central core perpendicular to the emanating outflow. Assuming equilibrium between centrifugal and gravitational force, the estimated mass for this rotating structure is 16/(sin2(i))Msun (with i the unknown inclination angle), of the same order as the gas mass derived from the continuum emission. A temperature estimate based on the CH3CN(19--18) K-ladder results in ~350K, thus a hot core has already formed in this region. An analysis of the submm continuum emission reveals even at this high spatial resolution only one isolated massive dust core without any detectable companions down to mass limits between 0.2 and 3Msun (depending on the assumed temperature). Potential implications for the massive cluster formation are discussed.

Molecular line observations of NH3 (J,K)=(1,1), (2,2) and CH3OH at 24.93GHz taken with the Australian Telescope Compact Array (ATCA) toward the massive twin cores NGC6334 I & I(N) reveal significant variations in the line emission between the two massive cores. The UCHII region/hot core NGC6334 I exhibits strong thermal NH3 and CH3OH emission adjacent to the UCHII region and coincident with two mm continuum peaks observed by Hunter et al. (in prep.). In contrast, we find neither compact NH3 nor thermal CH3OH line emission toward NGC6334 I(N). There, the NH3 emission is distributed over a broad region (>1') without a clear peak, and we find Class I CH3OH maser emission with peak brightness temperatures up to 7000K. The maser emission peaks appear to be spatially associated with the interfaces between the molecular outflows and the ambient dense gas. Peak NH3(1,1) line brightness temperatures >= 70K in both regions indicate gas temperatures of the same order. NH3 emission is also detected toward the outflow in NGC6334 I resulting in an estimated rotational temperature of Trot~19K. Furthermore, we observe CH3OH and NH3 absorption toward the UCHII region, the velocity structure is consistent with expanding molecular gas around the UCHII region. Thermal and kinematic effects possibly imposed from the UCHII region on the molecular core are also discussed.

Since this contributions was meant to cover two subjects which are both in the field of massive star formation but which in its details can be discussed separately, this paper is divided in two sections. First, we present characteristics of precursors of UCH{\sc ii} regions and their likely evolutionary properties. The second section discusses massive molecular outflows, their implications for high-mass star formation, and a possible evolutionary sequence for massive outflows.

At the earliest evolutionary stages, massive star-forming regions are deeply embedded within their natal cores and not observable at optical and near-infrared wavelengths. Interferometric high-spatial resolution mm dust continuum observations of one very young high-mass star-forming region disentangle its cluster-like nature already at the very beginning of the star formation process. The derived protocluster mass function is consistent with the stellar IMF. Hence, fragmentation of the initial massive cores may determine the IMF and the masses of the final stars. This implies that stars of all masses can form via accretion, and coalescence of protostars appears not to be necessary.

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 the first 865 mu continuum image with sub-arcsecond resolution obtained with the Submillimeter Array. These data resolve the Orion-KL region into the hot core, the nearby radio source I, the sub-mm counterpart to the infrared source n (radio source L), and new sub-mm continuum sources. The radio to submillimeter emission from source I may be modeled as either the result of proton-electron free-free emission that is optically thick to ~100 GHz plus dust emission that accounts for the majority of the submillimeter flux, or H- free-free emission that gives rise to a power-law spectrum with power-law index of ~1.6. The latter model would indicate similar physical conditions as found in the inner circumstellar environment of Mira variable stars. Future sub-arcsecond observations at shorter sub-mm wavelengths should easily discriminate between these two possibilities. The sub-mm continuum emission toward source n can be interpreted in the framework of emission from an accretion disk.

Submillimeter Array (SMA) observations of the high-mass star-forming region IRAS 18089-1732 in the 1 mm and 850 $\mu$m band with 1 GHz bandwidth reveal a wealth of information. We present the observations of 34 lines from 16 different molecular species. Most molecular line maps show significant contributions from the outflow, and only few molecules are confined to the inner core. We present and discuss the molecular line observations and outline the unique capabilities of the SMA for future imaging line surveys at high spatial resolution.

Studying outflows from young massive star-forming clusters allows one to deduce physical processes that lead to the formation of the most massive stars. I will review the current state of high-spatial-resolution interferometric (sub-)mm studies of massive molecular outflows and their implications for high-mass star formation. A possible evolutionary scenario for massive outflows with highly collimated structures at the very beginning and more wind-like outflows at later stages will be outlined.

We present initial results of an ongoing search for interferometric calibrators at submillimeter (sub-mm) wavelengths with the Submillimeter Array (SMA). Powerful radio galaxies are commonly used as calibrators at centimeter and millimeter wavelengths, but many are not strong enough to serve as calibrators at sub-mm wavelengths because of their rapidly declining flux densities toward shorter wavelengths. The inability to find a calibrator close to the target source may limit or even prevent us from imaging many interesting sources at sub-mm wavelengths. Here, we investigate whether high-mass protostellar objects and ultracompact HII regions can serve as useful calibrators for the SMA. The dust emission associated with these objects makes them among the brightest sub-mm sources in the sky. Our observations at 0.85 mm (345 GHz) with an angular resolution of ~3" reveal that although a large fraction of the dust emission originates from an extended ``halo'' component, a compact unresolved component often remains that when sufficiently strong may serve as a useful calibrator. These observations also provide a first glimpse at the small-scale distribution of dust around ultracompact HII regions and high-mass protostellar objects at sub-mm wavelengths. We discuss the origin of the core-halo structure seen in many sources, and conclude with suggestions for future searches for calibrators with the SMA.

We have mapped the ultracompact HII region, G5.89-0.39, and its molecular surroundings with the Submillimeter Array at 2".8 x 1".8 angular resolution in 1.3 mm continuum, SiO(5-4), and eight other molecular lines. We have resolved for the first time the highly energetic molecular outflow in this region. At this resolution, the outflow is definitely bipolar and appears to originate in a 1.3 mm continuum source. The continuum source peaks in the center of the HII region. The axis of the outflow lines up with a recently discovered O5V star.

We present high-spatial resolution Plateau de Bure Interferometer CO(2-1) and SiO(2-1) observations of one intermediate-mass and one high-mass star-forming region. The intermediate-mass region IRAS20293+3952 exhibits four molecular outflows, one being as collimated as the highly collimated jet-like outflows observed in low-mass star formation sources. Furthermore, comparing the data with additional infrared H2 and cm observations we see indications that the nearby ultracompact HII region triggers a shock wave interacting with the outflow. The high-mass region IRAS19217+1651 exhibits a bipolar outflow as well and the region is dominated by the central driving source. Adding two more sources from the literature, we compare position-velocity diagrams of the intermediate- to high-mass sources with previous studies in the low-mass regime. We find similar kinematic signatures, some sources can be explained by jet-driven outflows whereas other are better constrained by wind-driven models. The data also allow to estimate accretion rates varying from a few times 10^{-5}Msun/yr for the intermediate-mass sources to a few times 10^{-4}Msun/yr for the high-mass source, consistent with models explaining star formation of all masses via accretion processes.

SMA observations of the massive star-forming region IRAS 18089-1732 in the 1mm and 850mu band reveal outflow and disk signatures in different molecular lines. The SiO(5--4) data show a collimated outflow in the northern direction. In contrast, the HCOOCH3(20--19) line, which traces high-density gas, is confined to the very center of the region and shows a velocity gradient across the core. The HCOOCH3 velocity gradient is not exactly perpendicular to the outflow axis but between an assumed disk plane and the outflow axis. We interpret these HCOOCH3 features as originating from a rotating disk that is influenced by the outflow and infall. Based on the (sub-)mm continuum emission, the mass of the central core is estimated to be around 38M_sun. The dynamical mass derived from the HCOOCH3 data is 22Msun, of about the same order as the core mass. Thus, the mass of the protostar/disk/envelope system is dominated by its disk and envelope. The two frequency continuum data of the core indicate a low dust opacity index beta ~ 1.2 in the outer part, decreasing to beta ~ 0.5 on shorter spatial scales.

Studies of evolved massive stars indicate that they form in a clustered mode. During the earliest evolutionary stages, these regions are embedded within their natal cores. Here, we show high-spatial-resolution interferometric dust continuum observations disentangling the cluster-like structure of a young massive star-forming region. The derived protocluster mass distribution is consistent with the stellar initial mass function. Thus, fragmentation of the initial massive cores may determine the initial mass function and the masses of the final stars. This implies that stars of all masses can form via accretion processes, and coalescence of intermediate-mass protostars appears not to be necessary.

(Sub-)mm interferometry is the most favorable technique to investigate the earliest stages of massive star formation. I will outline general applications in that field and discuss results of different sub-topics (hot core chemistry and massive molecular outflows). Furthermore, recent data obtained with the Submillimeter Array will be shown to present the unique capabilities of this new instrument. Finally, I will give a short outlook on the main physical topics of massive star formation to be tackled with (sub-)mm interferometry within the next decade.

PdBI high-spatial resolution CO observations combined with near-infrared H2 data disentangle at least 7 (maybe even 9) molecular outflows in the massive star-forming region IRAS19410+2336. Position-velocity diagrams of the outflows reveal Hubble-like relationships similar to outflows driven by low-mass objects. Estimated accretion rates are of the order 10^-4 Msun/yr, sufficiently high to overcome the radiation pressure and form massive stars via disk-mediated accretion processes. The single-dish large-scale mm continuum cores fragment into several compact condensations at the higher spatial resolution of the PdBI which is expected due to the clustering in massive star formation. While single-dish data give a simplified picture of the source, sufficiently high spatial resolution resolves the structures into outflows resembling those of low-mass star-forming cores. We interpret this as further support for the hypothesis that massive stars do form via disk-accretion processes similar to low-mass stars.

We report the detection of hard X-ray emission (>2 keV) from a number of point sources associated with the very young massive star-forming region IRAS 19410+2336. The X-ray emission is detected from several sources located around the central and most deeply embedded mm continuum source, which remains undetected in the X-ray regime. All X-ray sources have K-band counterparts, and those likely belonging to the evolving massive cluster show near-infrared colors in the 2MASS data indicative of pre-main-sequence stages. The X-ray luminosities around 10^{31} erg/s are at the upper end of luminosities known for low-mass pre-main-sequence sources, and mass estimates based on the infrared data indicate that at least some of the X-ray detected sources are intermediate-mass objects. Therefore, we conclude that the X-ray emission is due to intermediate-mass pre-main-sequence Herbig Ae/Be stars or their precursors. The emission process is possibly due to magnetic star-disk interaction as proposed for their low-mass counterparts.

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.

We present a high-angular-resolution molecular line and millimeter continuum study of the massive star formation site IRAS 05358+3543. The most remarkable feature is a highly collimated (collimation factor ~10) and massive (>10 M_sun) bipolar outflow of 1 pc length, which is part of a quadrupolar outflow system. The three observed molecular outflows forming the IRAS 05358+3543 outflow system resemble, in structure and collimation, those typical of low-mass star-forming regions. They might therefore, just like low-mass outflows, be explained by shock entrainment models of jets. We estimate a mass accretion rate of 10^{-4) M_sun/yr, sufficient to overcome the radiative pressure of the central object and to build up a massive star, lending further support to the hypothesis that massive star formation occurs similarly to low-mass star formation, only with higher accretion rates and energetics.

We present a mapping study of massive molecular outflows in 26 high-mass star-forming regions at 11'' spatial resolution. Bipolar morpholgy is found in 80% of the sources and the collimation is higher than previously thought. Additionally, we find that well known low-mass correlations continue up to the high-mass regime, and accretion rates are around 10^(-4) Msun/yr rising as high as 10^(-3) Msun/yr. A tight correlation between the outflow and the core mass is established, implying that the accretion rate is a rather linear function of the initial core mass. All the data suggest that massive outflows are qualitatively similar and just quantitatively enhanced compared to their low mass counterparts.

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.

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