Search

We have used data from the Outer Galaxy High-Resolution Survey (OGHReS) to refine the velocities, distances, and physical properties of a large sample of 3584 clumps detected in far infrared/submillimetre emission in the HiGAL survey located in the $\ell = 250^\circ-280^\circ$ region of the Galactic plane. Using $^{12}$CO and $^{13}$CO spectra, we have determined reliable velocities to 3412 clumps (95% of the sample). In comparison to the velocities from the HiGAL catalogue, we find good agreement for 80% of the sample (within 5 km/s). Using the higher resolution and sensitivity of OGHReS has allowed us to correct the velocity for 632 clumps and provide velocities for 687 clumps for which no velocity had been previously allocated. The velocities are used with a rotation curve to refine the distances to the clumps and to calculate the clumps' properties using a distance-dependent gas-to-dust ratio. We have determined reliable physical parameters for 3200 outer Galaxy dense clumps (~90% of the HiGAL sources in the region). We find a trend of decreasing luminosity-to-mass ratio with increasing Galactocentric distance, suggesting the star formation efficiency is lower in the outer Galaxy or that it is resulting in more lower mass stars than in the inner Galaxy. We also find a similar surface density for protostellar clumps located in the inner and outer Galaxy, revealing that the surface density requirements for star formation are the same across the Galactic disc., Comment: 15 pages, 13 figures. Accepted for publication by MNRAS

Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions., Comment: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

The evolutionary sequence for high-mass star formation starts with massive starless clumps that go on to form protostellar, young stellar objects and then compact HII regions. While there are many examples of the three later stages, the very early stages have proved to be elusive. We follow-up a sample of 110 mid-infrared dark clumps selected from the ATLASGAL catalogue with the IRAM telescope in an effort to identify a robust sample of massive starless clumps. We have used the HCO+ (1-0) and HNC (1-0) transitions to identify clumps associated with infall motion and the SiO (2-1) transition to identity outflow candidates. We have found blue asymmetric line profile in 65% of the sample, and have measured the infall velocities and mass infall rates (0.6-$36 \times 10^{-3}$ Msun/yr) for 33 of these clumps. We find a trend for the mass infall rate decreasing with an increase of bolometric luminosity to clump mass i.e. star formation within the clumps evolves. Using the SiO 2-1 line, we have identified good outflow candidates. Combining the infall and outflow tracers reveals that 67% of quiescent clumps are already undergoing gravitational collapse or are associated with star formation; these clumps provide us with our best opportunity to determined the initial conditions and study the earliest stages of massive star formation. Finally, we provide an overview of a systematic high-resolution ALMA study of quiescent clumps selected that allows us to develop a detailed understanding of earliest stages and their subsequent evolution., Comment: 10 pages and 7 figures

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

We exploit the increased sensitivity of the recently installed AO SOUL at the LBT to obtain new high-spatial-resolution NIR images of the massive young stellar object IRAS20126+4104 and its outflow. We aim to derive the jet proper motions and kinematics, as well as to study its photometric variability by combining the novel performances of SOUL together with previous NIR images. We used both broad-band ($K_{s}$, $K'$) and narrow-band (Br$\gamma$, H2) observations from a number of NIR cameras (UKIRT/UFTI,SUBARU/CIAO,TNG/NICS,LBT/PISCES,and LBT/LUCI1) to derive maps of the continuum and the H$_2$ emission in the 2.12 $\mu$m line. Three sets of images, obtained with AO systems (CIAO,2003; FLAO,2012; SOUL,2020), allowed us to derive the proper motions of a large number of H$_2$ knots along the jet. Photometry from all images was used to study the jet variability. We derived knot proper motions in the range of 1.7-20.3 mas yr$^{-1}$ (i.e. 13-158 km s$^{-1}$ at 1.64 kpc, avg. outflow tangential velocity $\sim$ 80 km s$^{-1}$). The derived knot dynamical age spans a $\sim$ 200-4000 yr interval. A ring-like H$_2$ feature near the protostar location exhibits peculiar kinematics and may represent the outcome of a wide-angle wind impinging on the outflow cavity. Both H$_2$ geometry and velocities agree with those inferred from proper motions of the H$_2$O masers, located at a smaller distance from the protostar. Although the total H$_2$ line emission from the knots does not exhibit time variations at a $\widetilde{>}$ 0.3 mag level, we have found a clear continuum flux variation (radiation scattered by the dust in the cavity opened by the jet) which is anti-correlated between the blue-shifted and red-shifted lobes and may be periodic (with a period of $\sim$ 12-18 yr). We suggest that the continuum variability might be related to inner-disc oscillations which have also caused the jet precession., Comment: 26 pages, 22 figures, 2 mpeg files, accepted by Astronomy & Astrophysics

We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr$^{-1}$, of which $1.7 \pm 0.6$~M$_{\odot}$~yr$^{-1}$ coming from clumps with reliable heliocentric distance assignment. This value is in general agreement with estimates found in the literature of last decades. The profile of SFR density averaged in Galactocentric rings is found to be qualitatively similar to others previously computed, with a peak corresponding to the Central Molecular Zone and another one around Galactocentric radius $R_\mathrm{gal} \sim 5$~kpc, followed by an exponential decrease as $\log(\Sigma_\mathrm{SFR}/[\mathrm{M}_\odot~\mathrm{yr}^{-1}~\mathrm{kpc}^{-2}])=a\,R_\mathrm{gal}/[\mathrm{kpc}]+b $, with $a=-0.28 \pm 0.01$. In this regard, the fraction of SFR produced within and outside the Solar circle is 84\% and 16\%, respectively; the fraction corresponding to the far outer Galaxy ($R_\mathrm{gal} > 13.5$~kpc) is only 1\%. We also find that, for $R_\mathrm{gal}>3$~kpc, our data follow a power law as a function of density, similarly to the Kennicutt-Schmidt relation. Finally, we compare the distribution of the SFR density across the face-on Galactic plane and those of median parameters, such as temperature, luminosity/mass ratio and bolometric temperature, describing the evolutionary stage of Hi-GAL clumps. We found no clear correlation between the SFR and the clump evolutionary stage., Comment: To be published in ApJ

The Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium (SEDIGISM) survey has produced high (spatial and spectral) resolution $^{13}$CO (2-1) maps of the Milky Way. It has allowed us to investigate the molecular interstellar medium in the inner Galaxy at an unprecedented level of detail and characterise it into molecular clouds. In a previous paper, we have classified the SEDIGISM clouds into four morphologies. However, how the properties of the clouds vary for these four morphologies is not well understood. Here, we use the morphological classification of SEDIGISM clouds to find connections between the cloud morphologies, their integrated properties, and their location on scaling relation diagrams. We observe that ring-like clouds show the most peculiar properties, having, on average, higher masses, sizes, aspect ratios and velocity dispersions compared to other morphologies. We speculate that this is related to the physical mechanisms that regulate their formation and evolution, for example, turbulence from stellar feedback can often results in the creation of bubble-like structures. We also see a trend of morphology with virial parameter whereby ring-like, elongated, clumpy and concentrated clouds have virial parameters in a decreasing order. Our findings provide a foundation for a better understanding of the molecular cloud behaviour based on their measurable properties., Comment: Accepted for publication in A&A. 28 pages (17 of Appendices), 32 figures, 6 tables

We present one of the very first extensive classifications of a large sample of molecular clouds based on their morphology. This is achieved using a recently published catalogue of 10663 clouds obtained from the first data release of the SEDIGISM survey. The clouds are classified into four different morphologies by visual inspection and using an automated algorithm -- J plots. The visual inspection also serves as a test for the J plots algorithm, as this is the first time it has been used on molecular gas. Generally, it has been found that the structure of molecular clouds is highly filamentary and our observations indeed verify that most of our molecular clouds are elongated structures. Based on our visual classification of the 10663 SEDIGISM clouds, 15% are ring-like, 57% are elongated, 15% are concentrated and 10% are clumpy clouds. The remaining clouds do not belong to any of these morphology classes and are termed unclassified. We compare the SEDIGISM molecular clouds with structures identified through other surveys, i.e. ATLASGAL elongated structures and the bubbles from Milky Way Project (MWP). We find that many of the ATLASGAL and MWP structures are velocity coherent. ATLASGAL elongated structures overlap with ~ 21% of the SEDIGISM elongated structures (elongated and clumpy clouds) and MWP bubbles overlap with ~ 25% of the SEDIGISM ring-like clouds. We also analyse the star-formation associated with different cloud morphologies using two different techniques. The first technique examines star formation efficiency (SFE) and the dense gas fraction (DGF), based on SEDIGISM clouds and ATLASGAL clumps data. The second technique uses the high-mass star formation (HMSF) threshold for molecular clouds. The results indicate that clouds with ring-like and clumpy morphologies show a higher degree of star formation., Comment: Accepted for publication in A&A. 40 pages (26 of Appendices), 55 figures, 13 tables. The updated SEDIGISM cloud catalogue, containing cloud morphology, will be available as part of the SEDIGISM database

Massive stars form deeply embedded in their parental clouds, making it challenging to directly observe these stars and their immediate environments. It is known that accretion and ejection processes are intrinsically related, thus observing massive protostellar outflows can provide crucial information about the processes governing massive star formation close to the central engine. We aim to probe the IRAS 18264-1152 (G19.88-0.53) high-mass star-forming complex in the near infrared (NIR) through its molecular hydrogen (H2) jets to analyse the morphology and composition of the line emitting regions and to compare with other outflow tracers. We observed the H2 NIR jets via K-band (1.9-2.5um) observations obtained with the integral field units VLT/SINFONI and VLT/KMOS. SINFONI provides the highest NIR angular resolution achieved so far for the central region (~0.2''). We compared the geometry of the NIR outflows with that of the associated molecular outflow probed by CO (2-1) emission mapped with SMA. We identify nine point sources. Four of these display a rising continuum in the K-band and are BrG emitters, revealing that they are young, potentially jet-driving sources. The spectro-imaging analysis focusses on the H2 jets, for which we derived visual extinction, temperature, column density, area, and mass. The intensity, velocity, and excitation maps based on H2 emission strongly support the existence of a protostellar cluster, with at least two (and up to four) different large-scale outflows. The literature is in agreement with the outflow morphology found here. We derived a stellar density of ~4000 stars pc^-3. Our study reveals the presence of several outflows driven by young sources from a forming cluster of young, massive stars. The derived stellar number density together with the geometry of the outflows suggest that stars can form in a relatively ordered manner in this cluster., Comment: 23 pages, 14 figures. Published in A&A. Updated to correct typos and match proofs

ATLASGAL is a 870-mircon dust survey of 420 square degrees of the inner Galactic plane and has been used to identify ~10 000 dense molecular clumps. Dedicated follow-up observations and complementary surveys are used to characterise the physical properties of these clumps, map their Galactic distribution and investigate the evolutionary sequence for high-mass star formation. The analysis of the ATLASGAL data is ongoing: we present an up-to-date version of the catalogue. We have classified 5007 clumps into four evolutionary stages (quiescent, protostellar, young stellar objects and HII regions) and find similar numbers of clumps in each stage, suggesting a similar lifetime. The luminosity-to-mass (L/M) ratio curve shows a smooth distribution with no significant kinks or discontinuities when compared to the mean values for evolutionary stages indicating that the star-formation process is continuous and that the observational stages do not represent fundamentally different stages or changes in the physical mechanisms involved. We compare the evolutionary sample with other star-formation tracers (methanol and water masers, extended green objects and molecular outflows) and find that the association rates with these increases as a function of evolutionary stage, confirming that our classification is reliable. This also reveals a high association rate between quiescent sources and molecular outflows, revealing that outflows are the earliest indication that star formation has begun and that star formation is already ongoing in many of the clumps that are dark even at 70 micron., Comment: Accepted for publication in MNRAS. Consists of 20 pages, 15 figures, 8 table. The complete tables will be available from CDS and upon request

Filaments are a ubiquitous morphological feature of the molecular interstellar medium and are identified as sites of star formation. In recent years, more than 100 large-scale filaments (with a length $>10$\,pc) have been observed in the inner Milky Way. As they appear linked to Galactic dynamics, studying those structures represents an opportunity to link kiloparsec-scale phenomena to the physics of star formation, which operates on much smaller scales. In this letter, we use newly acquired Outer Galaxy High Resolution Survey (OGHReS) $^{12}$CO(2-1) data to demonstrate that a significant number of large-scale filaments are present in the outer Galaxy as well. The 37 filaments identified appear tightly associated with inter-arm regions. In addition, their masses and linear masses are, on average, one order of magnitude lower than similar-sized molecular filaments located in the inner Galaxy, showing that Milky Way dynamics is able to create very elongated features in spite of the lower gas supply in the Galactic outskirts., Comment: 5 pages (+ appendices, 9 pages), 8 figures, 1 table, accepted for publication in A&A Letters

Innovative developments in data processing, archiving, analysis, and visualization are nowadays unavoidable to deal with the data deluge expected in next-generation facilities for radio astronomy, such as the Square Kilometre Array (SKA) and its precursors. In this context, the integration of source extraction and analysis algorithms into data visualization tools could significantly improve and speed up the cataloguing process of large area surveys, boosting astronomer productivity and shortening publication time. To this aim, we are developing a visual analytic platform (CIRASA) for advanced source finding and classification, integrating state-of-the-art tools, such as the CAESAR source finder, the ViaLactea Visual Analytic (VLVA) and Knowledge Base (VLKB). In this work, we present the project objectives and the platform architecture, focusing on the implemented source finding services., Comment: 16 pages, 6 figures

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources., Comment: 16 pages, 11 figures, Accepted for publication in A&A

(Abridged) Understanding the details of the formation process of massive (i.e. M<8-10M$_\odot$) stars is a long-standing problem in astrophysics. [...] We present a method to derive accurate timescales of the different evolutionary phases of the high-mass star formation process. We model a representative number of massive clumps of the ATLASGAL-TOP100 sample which cover all the evolutionary stages. The models describe an isothermal collapse and the subsequent warm-up phase, for which we follow their chemical evolution. The timescale of each phase is derived by comparing the results of the models with the properties of the sources of the ATLASGAL-TOP100 sample, taking into account the mass and luminosity of the clumps, and the column densities of methyl acetylene (CH$_3$CCH), acetonitrile (CH$_3$CN), formaldehyde (H$_2$CO) and methanol (CH$_3$OH). We find that the chosen molecular tracers are affected by the thermal evolution of the clumps, showing steep ice evaporation gradients from 10$^3$ to 10$^5$ AU during the warm-up phase. We succeed in reproducing the observed column densities of CH$_3$CCH and CH$_3$CN, while H$_2$CO and CH$_3$OH show a poorer agreement with the observed values. The total (massive) star formation time is found to be $\sim5.2\times10^5$ yr, which is defined by the timescales of the individual evolutionary phases of the ATLASGAL-TOP100 sample: $\sim5\times10^4$ yr for 70-$\mu$m weak, $\sim1.2\times10^5$ yr for mid-IR weak, $\sim2.4\times10^5$ yr for mid-IR bright and $\sim1.1\times10^5$ yr for HII-regions phases. Our models, with an appropriate selection of molecular tracers that can act as chemical clocks, allow to get robust estimates of the duration of the individual phases of the high-mass star formation process, with the advantage of being capable to include additional tracers aimed at increasing the accuracy of the estimated timescales., Comment: Published on A&A (19 pages, 9 figures, 7 tables)

We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again., Comment: 29 pages, 14 figures, accepted for publication in A&A

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively., Comment: 25 pages, 16 figures

The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical 1-sigma sensitivity of 0.8-1.0 K at 0.25 km/s velocity resolution, it gives access to a wide range of structures, from individual star-forming clumps to giant molecular clouds and complexes. The coverage includes a good fraction of the first and fourth Galactic quadrants, allowing us to constrain the large scale distribution of cold molecular gas in the inner Galaxy. In this paper we provide an updated overview of the full survey and the data reduction procedures used. We also assess the quality of these data and describe the data products that are being made publicly available as part of this first data release (DR1). We present integrated maps and position-velocity maps of the molecular gas and use these to investigate the correlation between the molecular gas and the large scale structural features of the Milky Way such as the spiral arms, Galactic bar and Galactic centre. We find that approximately 60 per cent of the molecular gas is associated with the spiral arms and these appear as strong intensity peaks in the derived Galactocentric distribution. We also find strong peaks in intensity at specific longitudes that correspond to the Galactic centre and well known star forming complexes, revealing that the 13CO emission is concentrated in a small number of complexes rather than evenly distributed along spiral arms., Comment: MNRAS, in press

We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the global properties of clouds using a science sample, consisting of 6664 well resolved sources and for which the distance estimates are reliable. In particular, we compare the scaling relations retrieved from SEDIGISM to those of other surveys, and we explore the properties of clouds with and without high-mass star formation. Our results suggest that there is no single global property of a cloud that determines its ability to form massive stars, although we find combined trends of increasing mass, size, surface density and velocity dispersion for the sub-sample of clouds with ongoing high-mass star formation. We then isolate the most extreme clouds in the SEDIGISM sample (i.e. clouds in the tails of the distributions) to look at their overall Galactic distribution, in search for hints of environmental effects. We find that, for most properties, the Galactic distribution of the most extreme clouds is only marginally different to that of the global cloud population. The Galactic distribution of the largest clouds, the turbulent clouds and the high-mass star-forming clouds are those that deviate most significantly from the global cloud population. We also find that the least dynamically active clouds (with low velocity dispersion or low virial parameter) are situated further afield, mostly in the least populated areas. However, we suspect that part of these trends may be affected by some observational biases, and thus require further follow up work in order to be confirmed., Comment: 25 pages (+ appendices, 15 pages), 26 figures, MNRAS

By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurrent paper to estimate the dense gas fractions (DGF$_{\rm gmc}=\sum M_{\rm clump}/M_{\rm gmc}$) and the instantaneous star forming efficiencies (i.e., SFE$_{\rm gmc} = \sum L_{\rm clump}/M_{\rm gmc}$). We find that the molecular material associated with ATLASGAL clumps is concentrated in the spiral arms ($\sim$60% found within $\pm$10 km s$^{-1}$ of an arm). We have searched for variations in the values of these physical parameters with respect to their proximity to the spiral arms, but find no evidence for any enhancement that might be attributable to the spiral arms. The combined results from a number of similar studies based on different surveys indicate that, while spiral-arm location plays a role in cloud formation and HI to H$_2$ conversion, the subsequent star formation processes appear to depend more on local environment effects. This leads us to conclude that the enhanced star formation activity seen towards the spiral arms is the result of source crowding rather than the consequence of a any physical process., Comment: 15 pages, 12 figures. Full list of affiliations can be found at the end of the paper

We have used catalogues from several Galactic plane surveys and dedicated observations to investigate the relationship between various maser species and Galactic star forming clumps, as identified by the ATLASGAL survey. The maser transitions of interest are the 6.7 & 12.2 GHz methanol masers, 22.2 GHz water masers, and the masers emitting in the four ground-state hyperfine structure transitions of hydroxyl. We find clump association rates for the water, hydroxyl and methanol masers to be 56, 39 and 82 per cent respectively, within the Galactic longitude range of 60{\deg} > $l$ > -60{\deg}. We investigate the differences in physical parameters between maser associated clumps and the full ATLASGAL sample, and find that clumps coincident with maser emission are more compact with increased densities and luminosities. However, we find the physical conditions within the clumps are similar for the different maser species. A volume density threshold of $n$(H$_{2}$) > 10$^{4.1}$ cm$^{-3}$ for the 6.7 GHz methanol maser found in our previous study is shown to be consistent across for all maser species investigated. We find limits that are required for the production of maser emission to be 500 L$_{\odot}$ and 6 M$_{\odot}$ respectively. The evolutionary phase of maser associated clumps is investigated using the L/M ratio of clumps coincident with maser emission, and these have similar L/M ranges (~10$^{0.2}$ - 10$^{2.7}$ L$_{\odot}$/M$_{\odot}$) regardless of the associated transitions. This implies that the conditions required for the production of maser emission only occur during a relatively narrow period during a star's evolution. Lower limits of the statistical lifetimes for each maser species are derived, ranging from ~0.4 - 2 x 10$^{4}$ yrs and are in good agreement with the "straw man" evolutionary model previously presented., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 17 pages, 14 figures, 5 tables

Methanol (CH3OH) is one of the most abundant interstellar molecules, offering a vast number of transitions to be studied, including many maser lines. While the strongest Galactic CH3OH lines, the so-called class II masers, show no indications for the presence of superluminous counterparts in external galaxies, the less luminous Galactic class I sources appear to be different. Here we report class I 36GHz CH3OH 4(-1) - 3(0) E line emission from the nearby galaxies Maffei2 and IC342, measured with the 100-m telescope at Effelsberg at three different epochs within a time span of about five weeks. The 36GHz methanol line of Maffei2 is the second most luminous among the sources detected with certainty outside the Local Group of galaxies. This is not matched by the moderate infrared luminosity of Maffei2. Higher resolution data are required to check whether this is related to its prominent bar and associated shocks. Upper limits for M82, NGC4388, NGC5728 and Arp220 are also presented. The previously reported detection of 36GHz maser emission in Arp220 is not confirmed. Non-detections are reported from the related class I 44GHz methanol transition towards Maffei2 and IC342, indicating that this line is not stronger than its 36GHz counterpart. In contrast to the previously detected 36GHz CH3OH emission in NGC253 and NGC4945, our 36GHz profiles towards Maffei2 and IC342 are similar to those of previously detected non-masing lines from other molecular species. However, by analogy to our Galactic center region, it may well be possible that the 36GHz methanol lines in Maffei~2 and IC~342 are composed of a large number of faint and narrow maser features that remain spatially unresolved. In view of this, a search for a weak broad 36GHz line component would also be desirable in NGC253 and NGC4945., Comment: Accepted by Astronomy & Astrophysics, 8 pages, 1 table, 3 figures, slight changes w.r.t. version 1 due to language editing

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., Comment: 22 pages, 13 figures, accepted for publication in A&A

An estimate of the degree of CO-depletion ($f_D$) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen onto dust grains. A strong depletion state (i.e. $f_D>10$, as assumed in our models) is highly favoured in the innermost regions of dark clouds, where the temperature is $<20$ K and the number density of molecular hydrogen exceeds a few $\times$10$^{4}$ cm$^{-3}$. In this work, we estimate the size of the depleted region by studying the Infrared Dark Cloud (IRDC) G351.77-0.51. Continuum observations performed with the $Herschel$ $Space$ $Observatory$ and the $LArge$ $APEX$ $BOlometer$ $CAmera$, together with APEX C$^{18}$O and C$^{17}$O J=2$\rightarrow$1 line observations, allowed us to recover the large-scale beam- and line-of-sight-averaged depletion map of the cloud. We built a simple model to investigate the depletion in the inner regions of the clumps in the filament and the filament itself. The model suggests that the depletion radius ranges from 0.02 to 0.15 pc, comparable with the typical filament width (i.e.$\sim$0.1 pc). At these radii, the number density of H$_2$ reaches values between 0.2 and 5.5$\times$10$^{5}$ cm$^{-3}$. These results provide information on the approximate spatial scales on which different chemical processes operate in high-mass star-forming regions and also suggest caution when using CO for kinematical studies in IRDCs., Comment: accepted to MNRAS (14 pages, 10 figures, 1 table)

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., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

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

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

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.

Nitrogen (N) fractionation is used as a tool to search for a link between the chemical history of the Solar System and star-forming regions. A large variation of $^{14}$N/$^{15}$N is observed towards different astrophysical sources, and current chemical models cannot reproduce it. With the advent of high angular resolution radiotelescopes it is now possible to search for N-fractionation at core scales. We present IRAM NOEMA observations of the J=1-0 transition of N$_{2}$H$^{+}$, $^{15}$NNH$^{+}$ and N$^{15}$NNH$^{+}$ towards the high-mass protocluster IRAS 05358+3543. We find $^{14}$N/$^{15}$N ratios that span from $\sim$100 up to $\sim$220 and these values are lower or equal than those observed with single-dish observations towards the same source. Since N-fractionation changes across the studied region, this means that it is regulated by local environmental effects. We find also the possibility, for one of the four cores defined in the protocluster, to have a more abundant $^{15}$NNH$^{+}$ with respect to N$^{15}$NNH$^{+}$. This is another indication that current chemical models may be missing chemical reactions or may not take into account other mechanisms, like photodissociation or grain surface chemistry, that could be important., Comment: 19 pages, 8 figures, 6 tables, 3 appendices Accepted in Monthly Notices of the Royal Astronomical Society Letters

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

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

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

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., Comment: 45 pages, 20 pages for Table of derived parameters

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., Comment: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region., Comment: Astronomy & Astrophysics, in press

It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars., Comment: 12 pages, 9 figures, 2 tables, minor revision included, accepted by Astronomy & Astrophysics

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., Comment: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

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., Comment: Manuscript submitted to A&A, and revised after first referee report (abstract abbreviated)

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

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

Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the $^{13}$CH$_3$OH and CH$_2$DOH distribution in the inner regions ($\sim$100 au) of the HH212 system in Orion B. To this end, we used ALMA Cycle 1 and Cycle 4 observations in Band 7 with angular resolution down to $\sim$0.15$"$. Results: We detected 6 lines of $^{13}$CH$_3$OH and 13 lines of CH$_2$DOH with upper level energies up to 438 K in temperature units. We derived a rotational temperature of (171 $\pm$ 52) K and column densities of 7$\times$10$^{16}$ cm$^{-2}$ ($^{13}$CH$_3$OH) and 1$\times$10$^{17}$ cm$^{-2}$ (CH$_2$DOH), respectively. Consequently, the D/H ratio is (2.4 $\pm$ 0.4)$\times$10$^{-2}$, a value lower by an order of magnitude with respect to what was previously measured using single dish telescopes toward protostars located in Perseus. Our findings are consistent with the higher dust temperatures in Orion B with respect to that derived for the Perseus cloud. The emission is tracing a rotating structure extending up to 45 au from the jet axis and elongated by 90 au along the jet axis. So far, the origin of the observed emission appears to be related with the accretion disk. Only higher spatial resolution measurements however, will be able to disentangle between different possible scenarios: disk wind, disk atmosphere, or accretion shocks., Comment: Astronomy & Astrophysics Letter, in press

The dependence of the proton-to-electron mass ratio, mu, on the local matter density was investigated using methanol emission in the dense dark cloud core L1498. Towards two different positions in L1498, five methanol transitions were detected and an extra line was tentatively detected at a lower confidence level in one of the positions. The observed centroid frequencies were then compared with their rest frame frequencies derived from least-squares fitting to a large data set. Systematic effects, as the underlying methanol hyperfine structure and the Doppler tracking of the telescope, were investigated and their effects were included in the total error budget. The comparison between the observations and the rest frame frequencies constrains potential mu variation at the level of Dmu/mu < 6 x 10^(-8), at a 3 sigma confidence level. For the dark cloud we determine a total CH3OH (A+E) beam averaged column density of 3-4 x 10^(12) cm(-2) (within roughly a factor of two), an E- to A-type methanol column density ratio of N(A-CH3OH)/N(E-CH3OH) = 1.00 +/- 0.15, a density of n(H2) = 3 x 10^5 cm^(-3) (again within a factor of two), and a kinetic temperature of Tkin = 6 +/- 1 K. In a kinetic model including the line intensities observed for the methanol lines, the n(H2) density is higher and the temperature is lower than that derived in previous studies based on different molecular species; the intensity of the 1_0 --> 1_-1 E line strength is not well reproduced., Comment: 10 pages, 6 figures, accepted for publication in MNRAS

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

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar., Comment: 9 pages, 4 figures, 6 tables, Astronomy and astrophysics, in press

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

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

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

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"., Comment: 30 pages. Accepted for publication in A&A. Abstract slightly shortened due to arXiv requirements

In the BHR71 region, two low-mass protostars drive two distinguishable outflows. They constitute an ideal laboratory to investigate the effects of shock chemistry and the mechanisms that led to their formation. We aim to define the morphology of the warm gas component of the BHR 71 outflow and at modelling its shocked component. We present the first far infrared Herschel images of the BHR71 outflow in the CO(14-13), H$_2$O (2$_{21}$-1$_{10}$), H$_2$O (2$_{12}$-1$_{01}$) and [OI] 145 $\mu$m, lines, revealing the presence of several knots of warm, shocked gas associated with fast outflowing gas. In two of these knots we performed a detailed study of the physical conditions by comparing a large set of transitions from several molecules to a grid of shock models. Herschel lines ratios in the outflow knots are quite similar, showing that the excitation conditions of the fast moving gas do not change significantly within the first $\sim$ 0.068 pc of the outflow, apart at the extremity of the southern blue-shifted lobe that is expanding outside the molecular cloud. Rotational diagram, spectral line profile and LVG analysis of the CO lines in knot A show the presence of two gas components: one extended, cold ($T\sim$80 K) and dense ($n$(H$_2$) = 3$\times$10$^5$-4$\times$10$^6$ cm$^{-3}$) and another compact (18 arcsec), warm ($T$ = 1700-2200 K) with slightly lower density ($n$(H$_2$) = (2-6)$\times$10$^4$ cm$^{-3}$). In the two brightest knots (where we performed shock modelling) we found that H$_2$ and CO are well fitted with non-stationary (young) shocks. These models, however, significantly underestimate the observed fluxes of [OI] and OH lines, but are not too far off those of H$_2$O, calling for an additional, possibly dissociative, J-type shock component. Our modelling indirectly suggests that an additional shock component exists, possibly a remnant of the primary jet, Comment: 16 pages, 15 figures

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