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 inner flat region is possibly caused by fragmentation of the large scale cores into smaller sub-sources, whereas the steeper outer power law distributions indicate finite sizes of the cores.Separating the sources into sub-samples suggests that in the earliest stages prior to the onset of massive star formation the intensity radial distributions are rather flat resembling the structure of intensity peaks in more quiescent molecular clouds. Then in the subsequent collapse and accretion phase the intensity distributions become centrally peaked with steep power law indices. In this evolutionary stage the sources show also the broadest C 34 S linewidth. During the following phase, when ultracompact Hii regions evolve, the intensity power law radial distributions flatten out again. This is probably caused by the ignited massive stars in the center which disrupt the surrounding cores.The mean inner power law intensity index m i (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.

Context. Thanks to its excellent 5100 m high site in Chajnantor, the Atacama Pathfinder Experiment (APEX) systematically explores the southern sky at submillimeter wavelengths, in both continuum and spectral line emission. Studying continuum emission from interstellar dust is essential to locating the highest density regions in the interstellar medium, and deriving their masses, column densities, density structures, and large-scale morphologies. In particular, the early stages of (massive) star formation remain poorly understood, mainly because only small samples of high-mass proto-stellar or young stellar objects have been studied in detail so far. Aims. 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. Only a systematic survey of the Galactic Plane can provide the statistical basis for unbiased studies. A well characterized sample of Galactic star-forming sites will deliver an evolutionary sequence and a mass function of high-mass, star-forming clumps. This systematic survey at submillimeter wavelengths also represents a preparatory work for Herschel and ALMA. Methods. 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 μm, with a beam size of 19. 2. Taking advantage of its large field of view (11. 4) and excellent sensitivity, we started an unbiased survey of the entire Galactic Plane accessible to APEX, with a typical noise level of 50−70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). Results. As a first step, we covered ∼95 deg 2 of the Galactic Plane. These data reveal ∼6000 compact sources brighter than 0.25 Jy, or 63 sources per square degree, 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 protoclusters. Other compact sources harbor hot cores, compact H ii regions, or young embedded clusters, thus tracing more evolved stages after massive stars have formed. 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 . In this first introductory paper, we show preliminary results from these ongoing observations, and discuss the mid-and long-term perspectives of the survey.

Abstract. With the aim of understanding the role of massive outflows in high-mass star formation, we mapped in the 12 CO J = 2 − 1 transition 26 high-mass star-forming regions at very early stages of their evolution. At a spatial resolution of 11 bipolar molecular outflows were found in 21 of them. The other five sources show confusing morphology but strong line wings. This high detection rate of bipolar structure proves that outflows common in low-mass sources are also ubiquitous phenomena in the formation process of massive stars. The flows are large, very massive and energetic, and the data indicate stronger collimation than previously thought. The dynamical timescales of the flows correspond well to the free-fall timescales of the associated cores. Comparing with correlations known for low-mass flows, we find continuity up to the high-mass regime suggesting similar flow-formation scenarios for all masses and luminosities. Accretion rate estimates in the 10 4 L range are around 10 −4 M yr −1 , higher than required for low-mass star formation, but consistent with high-mass star formation scenarios. Additionally, we find a tight correlation between the outflow mass and the core mass over many orders of magnitude. The strong correlation between those two quantities implies that the product of the accretion efficiency facc =Ṁacc/(Mcore/t ff ) and fr (the ratio between jet mass loss rate and accretion rate), which equals the ratio between jet and core mass (faccfr = Mjet/Mcore), is roughly constant for all core masses. This again indicates that the flow-formation processes are similar over a large range of masses. Additionally, we estimate median fr and facc values of approximately 0.2 and 0.01, respectively, which is consistent with current jet-entrainment models. To summarize, the analysis of the bipolar outflow data strongly supports theories which explain massive star formation by scaled up, but otherwise similar physical processes -mainly accretion -to their low-mass counterparts.

We describe a systematic program aimed at identifying and characterizing candidate high-mass protostellar 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 H 2 O and CH 3 OH maser emission and observed the CO J = 2 → 1 and several NH 3 lines toward all sources in our sample. Other molecular tracers were observed in a subsample.While dust continuum emission was detected in all sources, most of them show only weak or no emission at 3.6 cm. Where detected, the cm emission is frequently found to be offset from the mm emission, indicating that the free-free and dust emissions arise from different subsources possibly belonging to the same (proto)cluster. A comparison of the luminosities derived from the cm emission with bolometric luminosities calculated from the IRAS far-infrared fluxes shows that the cm emission very likely traces the most massive source, whereas the whole cluster contributes to the far-infrared luminosity. Estimates of the accretion luminosity indicate that a significant fraction of the bolometric luminosity is still due to accretion processes. The earliest stages of HMPO evolution we seek to identify are represented by dust cores without radio emission.Line wings due to outflow activity are nearly omnipresent in the CO observations, and the molecular line data indicate the presence of hot cores for several sources, where the abundances of various molecular species are elevated due to evaporation of icy grain mantles. Kinetic gas temperatures of 40 sources are derived from NH 3 (1,1) and (2,2) data, and we compare the results with the dust temperatures obtained from the IRAS data.Comparing the amount of dust, and hence the gas, associated with the HMPOs and with ultracompact Hii regions (UCHiis) we find that the two types of sources are clearly separated in mass-luminosity diagrams: for the same dust masses the UCHii regions have higher bolometric luminosities than HMPOs. We suggest that this is an evolutionary trend with the HMPOs being younger and reprocessing less (stellar) radiation in the IR than the more evolved UCHiis regions.These 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.

Aims. Our current knowledge of high-mass star formation is mainly based on follow-up studies of bright sources found by IRAS, and is thus biased against its earliest phases, inconspicuous at infrared wavelengths. We therefore started searching, in an unbiased way and in the closest high-mass star-forming complexes, for the high-mass analogs of low-mass pre-stellar cores and class 0 protostars. Methods. We have made an extensive 1.2 mm continuum mosaicing study of the Cygnus X molecular cloud complex using the MAMBO cameras at the IRAM 30 m telescope. The ∼ 3 • 2 imaged areas cover all the high-column density (A V ≥ 15 mag) clouds of this nearby (∼ 1.7 kpc) cloud complex actively forming OB stars. We then compared our millimeter maps with mid-infrared images, and have made SiO(2-1) follow-up observations of the best candidate progenitors of high-mass stars. Results. Our complete study of Cygnus X with ∼ 0.09 pc resolution provides, for the first time, an unbiased census of massive young stellar objects. We discover 129 massive dense cores (FWHM size ∼ 0.1 pc, M 1.2 mm = 4 − 950 M ⊙ , volume-averaged density ∼ 10 5 cm −3 ), among which ∼ 42 are probable precursors of high-mass stars. A large fraction of the Cygnus X dense cores (2/3 of the sample) remain undetected by the MSX satellite, regardless of the mass range considered. Among the most massive (> 40 M ⊙ ) cores, infrared-quiet objects are driving powerful outflows traced by SiO emission. Our study qualifies 17 cores as good candidates for hosting massive infrared-quiet protostars, while up to 25 cores potentially host high-luminosity infrared protostars. We fail to discover in the high-mass analogs of pre-stellar dense cores (∼ 0.1 pc, > 10 4 cm −3 ) in Cygnus X, but find several massive starless clumps (∼ 0.8 pc, 7 × 10 3 cm −3 ) that might be gravitationally bound. Conclusions. Since our sample is derived from a single molecular complex and covers every embedded phase of high-mass star formation, it gives the first statistical estimates of their lifetime. In contrast to what is found for low-mass class 0 and class I phases, the infrared-quiet protostellar phase of high-mass stars may last as long as their better-known high-luminosity infrared phase. The statistical lifetimes of high-mass protostars and pre-stellar cores (∼ 3 × 10 4 yr and < 10 3 yr) in Cygnus X are one and two order(s) of magnitude smaller, respectively, than what is found in nearby, low-mass star-forming regions. We therefore propose that high-mass pre-stellar and protostellar cores are in a highly dynamic state, as expected in a molecular cloud where turbulent processes dominate.