# The ATLASGAL survey: The sample of young massive cluster progenitors

T. Csengeri, Sylvain Bontemps, F. Wyrowski, S. T. Megeath, Motte Frédérique, A. Sanna, M. Wienen, K. M. Menten

► To cite this version:

T. Csengeri, Sylvain Bontemps, F. Wyrowski, S. T. Megeath, Motte Frédérique, et al. The ATLASGAL survey: The sample of young massive cluster progenitors. Astronomy & Astrophysics - A&A, 2017, 601, pp.id.A60. 10.1051/0004-6361/201628254. hal-01435306

HAL Id: hal-01435306

https://hal.science/hal-01435306v1

Submitted on 14 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The ATLASGAL survey: The sample of young massive cluster progenitors*

T. Csengeri¹, S. Bontemps², F. Wyrowski¹, S. T. Megeath³ F. Motte⁴,⁵, A. Sanna¹, M. Wienen¹, and K. M. Menten¹

¹ Max Planck Institute for Radioastronomy, Auf dem Hügel 69, 53121 Bonn, Germany
e-mail: ctimea@mpifr-bonn.mpg.de

² OASU/LAB-UMR5804, CNRS, Université Bordeaux 1, 33270 Floirac, France

³ Ritter Astrophysical Research Center, Department of Physics and Astronomy, University of Toledo, 2801 W. Bancroft St., Toledo, OH 43606, USA

⁴ Institut de Planétologie et d'Astrophysique de Grenoble, Univ. Grenoble Alpes - CNRS-INSU, BP 53, 38041 Grenoble Cedex 9, France

⁵ Laboratoire AIM Paris Saclay, CEA-INSU/CNRS-Université Paris Diderot, IRFU/SAp CEA-Saclay, 91191 Gif-sur-Yvette, France

Received 4 February 2016 / Accepted 16 December 2016

## ABSTRACT

**Context.** The progenitors of high-mass stars and clusters are still challenging to recognise. Only unbiased surveys, sensitive to compact regions of high dust column density, can unambiguously reveal such a small population of particularly massive and cold clumps.

**Aims.** Here we use the ATLASGAL survey to identify a sample of candidate progenitors of massive clusters in the inner Galaxy.

**Methods.** We characterise a flux limited sample of compact sources selected from the ATLASGAL survey. Sensitive mid-infrared data at 21–24 µm from the WISE and MIPSGAL surveys were explored to search for embedded objects, and complementary spectroscopic data were used to investigate their stability and their star formation activity.

**Results.** We identify an unbiased sample of infrared-quiet massive clumps in the Galaxy that potentially represent the earliest stages of massive cluster formation. An important fraction of this sample consists of sources that have not been studied in detail before. We first find that clumps hosting more evolved embedded objects and infrared-quiet clumps exhibit similar physical properties in terms of mass and size, suggesting that the sources are not only capable of forming high-mass stars, but likely also follow a single evolutionary track leading to the formation of massive clusters. The majority of the clumps are likely not in virial-equilibrium, suggesting collapse on the clump scale.

**Conclusions.** We identify the precursors of the most massive clusters in the Galaxy within our completeness limit, and argue that these objects are undergoing large-scale collapse. This is in line with the low number of infrared-quiet massive clumps and earlier findings that star formation, in particular for high-mass objects is a fast, dynamic process. We propose a scenario in which massive clumps start to fragment and collapse before their final mass is accumulated indicating that strong self-gravity and global collapse is needed to build up rich clusters and the most massive stars.

**Key words.** surveys – stars: massive – stars: formation

## 1. Introduction

The formation of massive stars and clusters remains an enigma in modern astrophysics. To explain the origin of high-mass stars, two theoretical scenarios are commonly invoked: the turbulent-core model (McKee & Tan 2003) and a more dynamic view with cluster formation (Bonnell et al. 2001; Bonnell & Bate 2006). The observational challenge to these theories is a systematic, Galaxy-wide search for the earliest phases of high-mass stars and massive protoclusters where the initial conditions can best be studied.

With “earliest phases” we refer to the starless/pre-stellar, and the subsequent phase where the submillimetre emission is still dominated by the massive core. This implies that no luminous massive young stellar object (MYSO) has yet emerged, and no ionising emission from embedded UC-H II regions. This phase is difficult to observe, however, mainly because of the

short duration of less than $7.5 \times 10^4$ yr at the scale of ~0.3 pc (Csengeri et al. 2014), and also because of the high complexity of active sites of high-mass star formation. Only a few candidates of massive starless cores have been claimed so far in the literature with masses up to $60 M_\odot$ in a size of ≈0.1 pc (Tan et al. 2013; Cyganowski et al. 2014), although towards the former object outflow activity has recently been discovered (Tan et al. 2016), and the absence of any molecular emission towards the latter object casts doubts on its Galactic origin. Confirmed examples of massive pre-stellar cores, as well as exploration of the even higher mass regime, are of extreme importance in order to constrain star formation models.

Early searches for the birthplaces of high-mass stars focused on IRAS selected infrared-bright sources, where low- and high-luminosity sources were defined (Molinari et al. 1996). High-luminosity sources were associated with UC-H II regions, i.e. the final stage of high-mass star formation, while the lower luminosity sources were associated with MYSOs. UC-H II regions have typical bolometric luminosities of $>10^4 L_\odot$, while (M)YSOs show a broad distribution peaking between $3 \times 10^3-2 \times 10^4 L_\odot$

<footnote>* Full Table 4 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/601/A60</footnote>
(e.g. Molinari et al. 2000; Mottram et al. 2011). In fact, UC-HII regions and lower luminosity sources fill a distinct space in a mass-luminosity diagram, strongly suggesting an evolutionary trend of more luminous objects with similar masses being more evolved than their lower luminosity counterparts (see Molinari et al. 1996; Sridharan et al. 2002; Elia et al. 2013).

Only with the advance of sensitive and large format detector arrays have (sub)millimetre surveys of entire molecular clouds and larger areas of the Galaxy revealed lower luminosity objects associated with very massive gas reservoirs (e.g. Motte et al. 2003, 2007; Mookerjea et al. 2004; López et al. 2011; Russeil et al. 2010; Tackenberg et al. 2012; Traficante et al. 2015; Svoboda et al. 2016). Systematic searches for the massive but colder and lower luminosity sources have begun. In an unbiased survey, Motte et al. (2007) has identified a sample of massive dense cores (MDCs) with masses between 40–200 $M_\odot$, sizes of ~0.1 pc and low luminosities (<10³ $L_\odot$) in the Cygnus-X region. Their high concentration of mass and low bolometric luminosity suggest that they host high-mass protostars, potential precursors to MYSOs. High angular-resolution observations of Bontemps et al. (2010) indeed identified a handful of individual high-mass protostars within MDCs with envelope masses in the range of 8–50 $M_\odot$. Studies of the outflow activity further suggest that these protostars are the high-mass equivalents of Class 0 protostars (Duarte-Cabral et al. 2013). Altogether, this is strong evidence that low luminosity MDCs host objects in the earliest stages of high-mass protostellar evolution. An unbiased, systematic study to identify such objects on Galactic scales is the necessary next step for studying the initial conditions for high-mass star and cluster formation.

In this paper we present the first systematic and unbiased study of the brightest submillimetre sources in the inner Galaxy selected from the 870 µm APEX Telescope LArge Survey of the GALaxy (ATLASGAL, Schuller et al. 2009; Csengeri et al. 2014); we estimate their physical properties and search for the most massive clumps still in an early evolutionary stage among them. This allows us to identify and study the properties of the youngest precursors of massive clusters in the Galaxy which potentially host high-mass protostars. Various samples have been studied from the ATLASGAL survey, mostly based on associations with HII regions and sources from the Red MSX Survey (RMS, Lumsden et al. 2013), and their flux densities at 8 µm (GLIMPSE, Benjamin et al. 2003), and 24 µm (MIPSGAL, Carey et al. 2009). Focusing on the 1st Galactic quadrant, Csengeri et al. (2016a) present the largest sample to date of 430 clumps in different evolutionary stages over a wide mass range. The sample studied in Giannetti et al. (2014) represents a subsample of this source selection limited to the hundred brightest sources in different evolutionary stages in both the 1st and the IVth Galactic quadrant. This strategic selection of ATLASGAL clumps over various evolutionary stages has been targeted by extensive spectroscopic follow-up observations; see details in Csengeri et al. (2016a). Other samples using selective star formation tracers of more evolved objects, such as clumps with HII regions, RMS sources, and associations with masers from the Methanol MultiBeam survey, have been discussed in a series of papers by Urquhart et al. (2013a,b, 2014).

The paper is organised as follows. In Sect. 2 we describe the flux limited sample of the brightest submillimetre sources from the ATLASGAL survey and investigate their statistical properties in Sect. 3. Using complementary 21–24 µm data, we identify a sample of massive clumps which are potentially the youngest massive objects known so far. In Sect. 4 we study their physical properties and discuss their potential to form high-mass stars.

## 2. Brightest submillimetre clumps from the ATLASGAL survey

The ATLASGAL survey (Schuller et al. 2009) imaged the Galactic plane between Galactic longitude $\ell = \pm 60^\circ$ and Galactic latitude $b = \pm 1.5^\circ$ at 870 µm with the LABOCA camera (Siringo et al. 2009) on the APEX Telescope (Güsten et al. 2006). In an extension of the survey, the $-80^\circ < \ell < -60^\circ$ and $-2^\circ < b < +1^\circ$ region was also mapped. Thus, all together the ATLASGAL survey imaged ~420 sq. deg. of the Galactic plane at a 19"/2 spatial resolution. The catalogue of compact sources is 97% complete down to ~260 $M_\odot$ at the distance of the Galactic centre, i.e., 8 kpc, and 1 $M_\odot$ at 500 pc for cold clumps (Csengeri et al. 2014).

We study here a flux limited sample of the ATLASGAL compact sources within the highest sensitivity part of the survey, corresponding to a Galactic longitude range of $|\ell| \le 60^\circ$. The initial sample was selected using a 5 Jy beam-averaged flux density threshold, which translates to ~650 $M_\odot$ (see Csengeri et al. 2014) at an average distance of ATLASGAL sources of $d \approx 4.5$ kpc (Wienen et al. 2015). We estimate the corresponding column density by

$$N(\mathrm{H}_2) = \frac{F_v R}{B_v(T_d) \Omega \kappa_v \mu_{\mathrm{H}_2} m_{\mathrm{H}}}[\mathrm{cm}^{-2}], \quad (1)$$

where $F_v$ is the beam-averaged peak flux density; $\Omega$ is the solid angle of the beam calculated by $\Omega = 1.13 \times \Theta^2$, where $\Theta$ is 19"/2; $\mu_{\mathrm{H}_2}$ is the mean molecular weight per hydrogen molecule and is equal to 2.8 (Kauffmann et al. 2008); $m_{\mathrm{H}}$ is the mass of a hydrogen atom; and $R$ corresponds to the gas-to-dust ratio of 100. Our flux density threshold thus corresponds to a column density ($N_{\mathrm{H}_2}$) of 4.9–40.8 × 10²³ cm⁻² using a centre frequency of 345 GHz at our resolution of 19"/2, and dust temperatures of 10–40 K. For the highest column density sources this translates to an extremely high visual extinction of $A_V > 521$ mag<fnref n="1" /> (see details in Sect. 2.3).

From the ATLASGAL compact source catalogue we find a total of 219 sources above the 5 Jy beam-averaged flux density threshold (Csengeri et al. 2014). In the subsequent steps we first provide distance estimates (Sect. 2.1), and then determine their physical properties (Sects. 2.2 and 2.3). Using mid-infrared data between 21–24 µm from the WISE point source catalogue (Wright et al. 2010) and photometry based on the MIPSGAL survey (Carey et al. 2009) and complemented with classification from the MYSO and UC-HII regions from the RMS (Urquhart et al. 2007) and CORNISH surveys (Hoare et al. 2012; Purcell et al. 2013), we then estimate the evolutionary stage based on the luminosity of the embedded protostars associated with the dust peaks (Sects. 2.4 and 2.5). This leads to the identification of infrared-quiet massive clumps in Sect. 2.6. In Sect. 2.7 we present molecular line data to characterise their associations with shocks from outflow activity to probe ongoing star formation.

### 2.1. Distances and association with massive molecular complexes

The large statistics of the ATLASGAL survey revealed that most of the compact sources, and correspondingly most of the star formation activity is associated with giant molecular complexes (GMCs; e.g. Csengeri et al. 2014). As a first step, therefore, we

<footnote n="1">¹ For the conversion we use $A_V = N_{\mathrm{H}_2}/0.94 \times 10^{21} \mathrm{cm}^{-2}$ (Bohlin et al. 1978).</footnote>
<table><caption>**Table 1.** Association of the brightest ATLASGAL sources with GMCs.</caption><thead><tr><th>Associated GMC or molecular cloud</th><th><i>d</i> [kpc]</th><th>Number of sources</th></tr></thead><tbody><tr><td>W51</td><td>5.41<sup>b</sup></td><td>8</td></tr><tr><td>W49</td><td>11.11<sup>c</sup></td><td>3</td></tr><tr><td>IRDC G34.4+0.2</td><td>1.56<sup>p</sup></td><td>3</td></tr><tr><td>W43</td><td>5.49<sup>d</sup></td><td>4</td></tr><tr><td>W31</td><td>4.95<sup>f</sup></td><td>3</td></tr><tr><td>W33</td><td>2.4<sup>e</sup></td><td>3</td></tr><tr><td>M17</td><td>1.98<sup>m</sup></td><td>4</td></tr><tr><td>G5.89</td><td>1.28<sup>r</sup></td><td>2</td></tr><tr><td>Galactic centre</td><td>8.5<sup>j</sup></td><td>9</td></tr><tr><td>NGC 6334</td><td>1.35<sup>i</sup></td><td>9</td></tr><tr><td>RCW 122</td><td>1.8<sup>i,j</sup></td><td>4</td></tr><tr><td>RCW 108-109</td><td>1.8<sup>a,l</sup></td><td>3</td></tr><tr><td>G333.1-0.4</td><td>4.2<sup>a,g</sup></td><td>6</td></tr><tr><td>RCW106</td><td>4.2<sup>a</sup></td><td>5</td></tr><tr><td>G327</td><td>2.5<sup>a</sup></td><td>3</td></tr><tr><td>RCW 95</td><td>2.5<sup>a,g,y</sup></td><td>5</td></tr><tr><td>RCW 92</td><td>3.2<sup>a</sup></td><td>2</td></tr><tr><td>G305</td><td>3.5<sup>q</sup></td><td>2</td></tr><tr><td>G10.47</td><td>8.47<sup>f</sup></td><td>1</td></tr><tr><td>RCW 120</td><td>1.3<sup>a</sup></td><td>1</td></tr><tr><td>RCW 131/NGC 6537</td><td>1.7<sup>h</sup></td><td>1</td></tr><tr><td>G345.5+1.0</td><td>1.8<sup>a,s</sup></td><td>1</td></tr><tr><td>RCW 116b</td><td>1.8<sup>a,g</sup></td><td>1</td></tr><tr><td>RCW 117</td><td>1.8<sup>a,k</sup></td><td>1</td></tr><tr><td>RCW 99</td><td>2.5<sup>a</sup></td><td>1</td></tr><tr><td>G9.62+0.19</td><td>5.15<sup>n</sup></td><td>1</td></tr><tr><td>G351.77</td><td>1.0<sup>o</sup></td><td>1</td></tr><tr><td>SDC335</td><td>3.8<sup>a,‡</sup></td><td>1</td></tr></tbody></table>

**Notes.** (a) This work; (b) Sato et al. (2010); (c) Zhang et al. (2013); (d) Zhang et al. (2014); (e) Immer et al. (2013); (f) Sanna et al. (2014); (g) Moisés et al. (2011); (h) we revised the spectrophotometric distance of 1.81 kpc by Moisés et al. (2011); (i) Wu et al. (2014); (j) Wu et al. (2012); (k) Reid et al. (2014); (l) Dutra et al. (2003); (m) Straw et al. (1987); (n) Xu et al. (2011); (o) Sanna et al. (2009); (p) Leurini et al. (2011); (q) Kurayama et al. (2011); (r) Russeil (2003); (s) Motogi et al. (2011); (t) López et al. (2011); (u) Peretto et al. (2013) adopts 3.25 kpc.

associated the sample with known GMCs (see Table 1). In fact, our selection corresponds to the brightest submillimetre sources in the inner Galactic plane, and it covers well-characterised massive star-forming complexes such as Sgr B2, W51 Main, W49A, W33, and W43. Other examples are found at the edge of expanding bubbles like G348.1825+0.4829, which is the brightest compact source in the vicinity of RCW 120 (Zavagno et al. 2010). We find a few examples of relatively isolated sources, which are not associated with any known GMCs, but smaller molecular clouds, such as G0.5464-0.8521, G333.0164+0.7654, and G43.7950-0.1270, towards which our distance estimates may be less reliable than for the majority of the sample.

We collected distance estimates from the literature using primarily maser parallax measurements (see Reid et al. 2014, for an overview) or spectrophotometric distance estimates (e.g. Russeil 2003; Moisés et al. 2011) for the known GMCs. For the remaining unassociated sources we used the high-density tracer NH₃ to determine kinematic distances from Wienen et al. (2012, 2015).

In short, the rest velocity of the source is compared to the Galactic rotation curve using the parameters derived by Reid et al. (2009) for the Ist and Brand & Blitz (1993) for the IVth Galactic quadrants, respectively. We attempt to solve the distance ambiguity by making use of Hα emission and 21 cm continuum absorption against bright HII regions (see Wienen et al. 2015, for details). We then visually inspected near-IR extinction maps (Schneider et al. 2010) and data from the Galactic Ring Survey (GRS, Jackson et al. 2006) to look for spatial and kinematic associations with nearby molecular clouds also considering the larger scale Galactic structure such as spiral arms, and associations with extinction features as well. We also inspected sensitive mid-infrared data from GLIMPSE and MIPSGAL to look for possible associations with infrared-dark dust lanes appearing in high contrast absorption at mid-infrared wavelengths.

In total we associate 78 (36%) sources with available maser parallax measurements. In the Ist Galactic quadrant, such measurements are available for 59 out of 85 sources (70%), providing a robust estimate of their physical parameters. Sources in the IVth Galactic quadrant have less reliable distance measurements owing to the lack of currently available maser parallax distances, except for G339.8841-1.2568 (Krishnan et al. 2015). We were able to provide distance estimates to all but one source in the sample (see Table 4). For simplicity, in the following analysis we exclude the 9 sources associated with the Galactic centre, leaving a total of 210 sources in our sample. An overview of their Galactic distribution is shown in Fig. 1.

## 2.2. Temperature estimates

Our sample includes clumps in various evolutionary stages such as active star-forming sites with clusters of known (UC-)HII regions (e.g. W51 Main, W49N, G5.89) and typical hot molecular cores (HMCs; e.g. G10.47+0.3, G34.26+0.15) associated with infrared-bright MYSOs. Over such a large star formation activity, the estimated physical parameters from the 870 µm dust emission can be highly sensitive to the used dust temperature. As a simple first-order approximation, a single dust temperature is commonly used for which a $T_d$ of 18 K is a reasonable choice (see e.g. Csengeri et al. 2014) based on the Herschel measured dust temperatures dominated by the interstellar radiation-field (Bernard et al. 2010). However, the brightest submillimetre sources are more likely to host embedded heating source(s), which could imply higher dust temperatures. This is supported by statistical studies; for instance, Battersby et al. (2011) show that tracers of star formation are more likely associated with higher dust temperatures, while König et al. (2017) show an increase in dust temperature with evolutionary stage by fitting the spectral energy distribution of ATLASGAL clumps over a range of evolutionary stages. Hence, a temperature correction is necessary for clumps hosting UC-HII regions and infrared-bright protostars, and we discuss below our choice of a higher dust temperature.

At the high densities probed by the ATLASGAL survey it is reasonable to assume that the gas and dust are thermalised. Therefore, to investigate which dust temperature is appropriate for the warmer clumps, we can also rely on complementary molecular line data available in the literature. The NH₃ measurements by Wienen et al. (2012) give a gas kinetic temperature between 12–35 K for a large sample of ATLASGAL sources within a beam of ~40". Similarly, Dunham et al. (2011) derive a mean gas kinetic temperature of 17 K for a sample of 199 sources from the Bolocam Galactic Plane Survey (BGPS) with a beam of 31". This corresponds well to the average dust
<figure><img src="image_2.png" /><figcaption>Fig. 1. Overview of the Galactic distribution of the sample and the associated molecular clouds. Figure adopted from Csengeri et al. (2016b) showing the ATLASGAL-Planck combined view of the inner galaxy. Magenta crosses show the location of the sources. The size of the symbols scales linearly with the peak flux density.</figcaption></figure>

temperature also used in various studies of ATLASGAL selected sources (e.g. Csengeri et al. 2016a).

A total of 127 sources (corresponding to 58% of the sample) have kinetic temperature estimates from Wienen et al. (2012) and Wienen et al. (2017), and the bulk of these temperature measurements range between 17 and 35 K. Since these measurements do not cover the full sample, they only allow us to estimate a reasonable temperature correction. Furthermore, as these studies mostly rely on the (1, 1) to (2, 2) inversion transitions of NH$_3$, they are more sensitive to the cold gas rather than the hot gas component. Therefore, for the more evolved objects these temperatures should be considered as lower limits on the dust temperature, and thus an upper limit on the mass. When discussing the statistical properties of the sample (Sect. 3), as a second-order approximation, we account for the increasing temperature towards sources with embedded massive stars using a higher T$_d$ of 30 K. This is consistent with the typical values used so far of 15–20 K for cold sources and up to 30–40 K for warm sources (e.g. Motte et al. 2007; Russeil et al. 2010).

## 2.3. Estimating physical parameters: mass, size, density, and surface density

We use the same formula as described in Schuller et al. (2009) and Csengeri et al. (2014) to estimate the gas mass from the optically thin 870 µm dust emission

$$M_{\text{gas}} = \frac{S_{\nu} R d^2}{B_{\nu}(T_d) \kappa_{\nu}} \approx \frac{M_{\text{gas}}}{[M_{\odot}]} = C \times \frac{S_{\nu}}{[\text{Jy}]} \times \left( \frac{d}{[\text{kpc}]} \right)^2, \quad (2)$$

where $S_ν$ is the integrated flux and d corresponds to the distance of the sources. The numerical constant, $C_{18\text{K}} = 6.33$

for $ν = 345$ GHz, $T_d = 18$ K, and $κ_ν = 0.0185 \text{ cm}^2 \text{ g}^{-1}$ from Ossenkopf & Henning (1994) accounting for a gas-to-dust ratio (R) of 100. We note that recent studies suggest a higher gas-to-dust ratio of 150 (Draine 2011), which would directly increase our mass estimates by 50%. Peretto et al. (2013) adopt a 50% lower $κ_ν$, which would also lead to an increase in the estimated mass. Following the analysis in Sect. 2.2, we introduce a temperature corrected mass estimate ($M_{\text{corr}}$) for the more evolved and thus likely warmer objects adopting $T_d = 30$ K, corresponding to $C_{30\text{K}} = 3.08$ from Sect. 3.

An additional factor of uncertainty in the mass determination originates from optically thick free-free emission which may contribute to the 870 µm flux density towards the clumps hosting UC-H II regions. For a uniform density H II region this contribution can be negligible because the thermal dust emission dominates the spectral energy distribution (Hunter et al. 2000). However, for a non-uniform density source with high temperature and density, a more significant contribution from free-free emission is possible (e.g. Beuther et al. 2004; Keto 2002, 2003). This may lead to an overestimate of the mass. Schuller et al. (2009), however, estimate this contribution to be considerably lower than the thermal emission from dust at 345 GHz, not more than 20% for the most extreme cases. In addition, the already formed ionising sources are still compact and therefore are expected to have a minor influence on the mass estimates. To conclude, for sources associated with known embedded UC-H II regions we conservatively consider the derived masses as upper limits (Table 4).

In Csengeri et al. (2014) the ATLASGAL sources have been fitted with a 2D Gaussian, from which we estimate the source size following the method of Bontemps et al. (2010). We define the source radius ($R_{90}$) which contains 90% of the total
mass, assuming a 2D Gaussian profile, by $R_{90} = 1.95 \times FWHM_{deconvolved}/\sqrt{8 \ln 2}$, where $FWHM_{deconvolved}$ corresponds to the deconvolved source size<fnref n="2" />. From the mass and size, we use the following equations to calculate the mean volume density and surface density:

$$ \frac{\bar{n}}{[\text{cm}^{-3}]} = \frac{M_{\text{gas}}}{\frac{4}{3} R_{90}^3 \pi} \qquad (3) $$

$$ \frac{\Sigma}{[\text{g cm}^{-2}]} = \frac{M_{\text{gas}}}{R_{90}^2 \pi} \qquad (4) $$

The estimated physical parameters are listed in Table 4.

## 2.4. Mid-infrared tracers of embedded protostars

Here we investigate the properties of the embedded YSOs closely associated with the dust peaks in order to reveal their star formation activity and, in particular, the youngest clumps within the sample. To classify an embedded source as a high-mass young stellar object, we rely on the IRAS-based colour properties for embedded high-mass stars from Wood & Churchwell (1989) following the method of Motte et al. (2007) and Russeil et al. (2010). Based on stellar models we adopt a luminosity threshold corresponding to objects with $L_{bol} > 10^4 L_\odot$, which potentially host a >15 $M_\odot$ zero age main sequence (ZAMS) star with a spectral type earlier than B0.5 (see also Mottram et al. 2011). In practice, using the expression from Motte et al. (2007), we translate this luminosity to a flux density threshold of ~289 Jy at 21 µm flux for an object at 1 kpc. This limit is a factor of ten higher than was used for Cygnus-X (Motte et al. 2007) and NGC 6334 (Russeil et al. 2010), and we evaluate the implication of the higher luminosity limit in Sect. 2.5. We extrapolate between the 22 µm and 24 µm data assuming blackbody radiation for the embedded protostar following Russeil et al. (2010).

Below this luminosity threshold we find high-mass protostars or their potential precursors, which are still deeply embedded in their mass reservoir, while their energetics is dominated by accretion processes. Although theoretical models are strongly dependent on the accretion rates, they predict Class 0-like high-mass protostars up to $L_{bol} < 10^4 L_\odot$ with an envelope mass of >50 $M_\odot$. Our threshold is also consistent with evolutionary tracks of high-mass stars assuming constant or varying accretion with rates of $10^{-4}-10^{-3} M_\odot \text{yr}^{-1}$ (Duarte-Cabral et al. 2013; Kuiper & Yorke 2013). In fact, the potentially highest mass protostar in Cygnus-X, N63, has an envelope mass of ~44 $M_\odot$ and $L_{bol} \sim 339 L_\odot$ estimated from SED models using Herschel data (Duarte-Cabral et al. 2013).

To extract the mid-infrared flux density of the sample, we first use the most sensitive and highest angular-resolution data from the MIPSGAL survey (Carey et al. 2009). The MIPS detector saturates at 2 Jy for point sources, and confusion due to nearby point sources or bright extended emission may hinder the photometry at the position of the dust peaks. Considering only the point source saturation limit, the MIPSGAL survey reveals B3 type (or earlier) stars only beyond 3.8 kpc. Less luminous but closer objects are likely to appear saturated in the MIPS images. Therefore, in order to cover a continuous sensitivity range it is necessary to complement the MIPSGAL data with data from the WISE satellite. The PSF photometry of the WISE point source catalogue is sensitive up to a ~330 Jy flux density threshold (Cutri et al. 2012).

The MIPSGAL survey at 24 µm has a $3\sigma$ point source sensitivity of 2 mJy at a 6'' angular resolution, while the WISE survey provides a $5\sigma$ point source sensitivity of 6 mJy at a 12'' resolution at 22 µm in unconfused regions of the Galaxy. Owing to extended emission from nebulosity and polycyclic aromatic hydrocarbons (PAHs), these sensitivity levels may be higher towards confused regions. The factor of 4 difference in beam areas may also result in different flux densities determined for the same source from the MIPSGAL and WISE data. In Appendix A we discuss the cross-calibration between the MIPS and WISE photometry.

Although ultimately the sensitivity threshold is only a factor of few better for the MIPSGAL survey, its angular-resolution (6''), which is more than three times higher, is a key benefit in our study. This angular resolution corresponds to a physical scale of ~0.14 pc at a 4.5 kpc distance, which is the size scale of massive dense cores embedded in massive clumps (Motte et al. 2007; Zhang et al. 2009; Csengeri et al. 2011b). This angular resolution allows us to resolve and distinguish between infrared-bright and -quiet massive dense cores for objects located at distances closer than ~4.5 kpc.

All sources were therefore first visually inspected using the publicly available calibrated data tiles from the MIPSGAL survey, and where no point source was found within 6'' of the ATLASGAL peak position, we extracted an upper limit on the flux density (for more details see Appendix A.1). Aperture photometry was performed on point sources where no saturated pixels are present (Appendix A.2). Where the MIPSGAL data was not useable, we used the flux density from the WISE point source catalogue Wright et al. (2010, Appendix A.3)<fnref n="3" />.

In Fig. 2 we show two examples of massive clumps with no (or weak) embedded mid-infrared sources, the upper panel showing the WISE 22 µm image of G340.9698-1.0212 and the lower panel showing the MIPSGAL 24 µm image towards G34.43+00.24 (Garay et al. 2004; Rathborne et al. 2005); the rest of the sample is shown in Appendix B, Fig. B.1. This particular example shows the importance of angular resolution because this region hosts an already formed UC-H II region in the vicinity of G34.4+0.2-MM1 and MM2, while the brightest dust peak itself is clearly shifted from the position of the UC-H II region (Rathborne et al. 2005). The other embedded dust peaks in this clump are active sites of high-mass star formation in their earliest phase (Rathborne et al. 2008).

In Fig. 3, we show the resulting photometry measurements versus clump mass using $T_d = 18$ K (see Sect. 2.3). The photometry is scaled to a common wavelength of 22 µm and to a distance of 1 kpc. Using a homogeneous dust temperature facilitates the comparison with the larger sample of sources discussed in Sect. 2.5. While the brightest mid-infrared sources correspond to active high-mass star-forming sites in the inner Galaxy, we reveal a population of massive clumps which are clearly fainter at mid-infrared wavelengths (see also Sect. 2.6).

## 2.5. Tracers of embedded high-mass stars

According to stellar evolution models, stars of types B3-B0.5 may also have luminosities between $L_{bol} = 10^3-10^4 L_\odot$, hence they fall below our threshold introduced in Sect. 2.4. We assess their potential contamination of the sample of the youngest massive clumps by using sensitive radio-surveys such as CORNISH (Hoare et al. 2012; Purcell et al. 2013) and

<footnote n="2">² We use $19''/2$ for the resolution of the ATLASGAL maps.</footnote>

<footnote n="3">³ The WISE 22 µm flux density is calculated from the W4 band magnitudes by using the formula $F_{22 \mu m} = 8.2839 \times 10^{\text{magw4}} / -2.5$.</footnote>
<figure><img src="image_3.png" /></figure>

**Fig. 2.** *Top:* example of an infrared-quiet massive clump with a single dust peak. The *left panel* shows the 22–24 µm image from WISE/MIPSGAL with contours from ATLASGAL. The white arrows mark the position of infrared-quiet ATLASGAL clumps (Sect. 2.6). The *right panel* shows the 870 µm dust emission from ATLASGAL with the same contours as in the *left panel*. *Bottom:* example of a massive clump towards IRDC 34.43+00.24 (Garay et al. 2004) with a positional shift between the dust and the 24 µm point source in the MIPSGAL image, which also coincides with a UC-H II region from CORNISH marked by a blue cross.

Urquhart et al. (2007) to search for ionising emission from stars in our sample of types later than B0.5. Since there is no homogeneous high angular-resolution radio dataset available covering the entire sample towards the IVth Galactic quadrant, we estimate a contamination level by performing the same analysis as above for a larger sample of CORNISH sources.

In Fig. 4 we show the same plot of distance scaled photometry at 22 µm versus mass as in Fig. 3 (see Sect. 2.4), but for a larger sample of ATLASGAL sources hosting embedded UC-H II regions from Urquhart et al. (2013b), as well as MYSOs from the RMS survey (Lumsden et al. 2013). This sample has been studied in multiwavelength follow-up observations (for a summary, see Lumsden et al. 2013), and the sources have been classified according whether they host (M)YSOs, H II regions, or a mixture of both. We extracted the 22 µm fluxes for these samples from the WISE catalogue considering matches between the ATLASGAL dust peaks and mid-infrared point sources with a small angular separation (<6′′). We considered only good quality flag photometry and a high signal-to-noise ratio of >10 in the 22 µm band. More importantly, for the RMS and CORNISH sources we only show those with the most reliable distance estimates based on maser parallax measurements (Reid et al. 2014).

We find that massive clumps ($M_{\text{clump}} > 650\,\text{M}_\odot$) hosting embedded UC-H II regions occupy a well-defined region, and mostly lie above our threshold of $S_{22\,\mu\text{m}} > 289\,\text{Jy}$ at 1 kpc. Out of the 78 sources of the shown ATLASGAL-CORNISH matches, 5 lie below this threshold, corresponding to 6.4% of the entire sample. Such weak, compact radio continuum emission has been found towards low bolometric luminosity ($L_{\text{bol}} \le 5 \times 10^4\,\text{L}_\odot$)

<figure><img src="image_4.png" /></figure>

**Fig. 3.** Distance scaled 22 µm flux density based on Spitzer/MIPS photometry and the WISE point source catalogue versus mass estimated from the ATLASGAL flux density. The mid-infrared flux densities have been scaled to a common wavelength. Symbols with arrows pointing left are sources with embedded UC-H II regions and upper limits for their mass estimate. Symbols with arrows pointing down correspond to upper limits for the 22 µm photometry. The dashed line shows the 289 Jy flux density threshold discussed in Sect. 2.6, and the dotted line corresponds to 650 $M_\odot$.

<figure><img src="image_5.png" /></figure>

**Fig. 4.** Mass versus 22 µm flux density of the sources covered by the CORNISH survey normalised to a distance of 1 kpc. Black filled circles show our sample with a left arrow indicating upper limits for the masses.

MYSOs (see e.g. Guzmán et al. 2010; Lumsden et al. 2013; Avison et al. 2015). They are likely associated with shock-ionisation from an early outflow activity (e.g. Hoare & Franco 2007).

Our flux limited selection of ATLASGAL sources with embedded UC-H II regions are also shown as a comparison, and find that 3 sources out of the 31 fall below our mid-infrared threshold, which is 9.7%. We inspected these sources visually, and find that one source lies in a confused region where either the mid-infrared photometry may be affected by extended emission, or the radio continuum emission may be associated with the extended H II region rather than the dust peak. The rest could correspond to the above discussed genuine low bolometric luminosity radio sources associated with MYSOs. To summarise, the
lowest luminosity massive clumps likely lack embedded UC-H II regions and bright MYSOs, and their contamination below our luminosity threshold is estimated to be less than ten per cent.

## 2.6. Classification of the sample

Using the information described in the previous sections, we classify our sample primarily based on their mid-infrared flux density and indicate the presence of UC-H II regions if radio continuum data is available. We do not rely on radio continuum data, however, because sensitive unbiased surveys are not yet available in the IVth quadrant, and thus the sample of UC-H II regions associated with our clumps is incomplete. We introduced the threshold of 289 Jy at 22 µm at 1 kpc as discussed in Sect. 2.4, and as empirically justified in Sect. 2.5. Sources above and below this threshold are termed “infrared bright” and “infrared quiet”, respectively. This classification is consistent, but simplified compared to previous works (e.g. Csengeri et al. 2014; Giannetti et al. 2014; Csengeri et al. 2016a). In addition, it quantitatively selects the young clumps based on the luminosity of the embedded protostar(s). Contamination of more evolved objects is estimated to be less than 10% (Sect. 2.5). Nevertheless, in the following we indicate the subsample of massive clumps hosting genuine UC-H II regions for the infrared-bright clumps.

From the sample of 210 sources we were able to classify 209 sources: 139 (66%) are infrared-bright including 51 (24%) unambiguously associated with UC-H II regions, and 80 (34%) are infrared quiet. Considering associations with infrared sources only from the WISE catalogue, we find 56 sources (26%) with no counterpart; the remaining 74% are associated with a mid-infrared source, which is consistent with earlier findings from Csengeri et al. (2014). Some objects are located in regions with complex extended emission at mid-infrared wavelengths, which may affect the photometry, hence their nature being quiet or bright may remain somewhat uncertain. As shown before, this definition of mid-infrared-quiet clumps picks up massive clumps hosting high-mass protostars or their precursors; where no ancillary radio data is available, the contamination from 8–15 $M_\odot$ main sequence stars is minor. We will investigate the nature of these sources further in Sect. 4. We show the images of infrared-quiet sources in Appendix B.

## 2.7. SiO (8–7) emission towards infrared-quiet massive clumps

To further investigate the properties of the infrared-quiet clumps (see also Sect. 4.1), we used spectroscopic observations of a typical shock tracer, SiO, as an indirect probe of ongoing accretion processes. Of the 80 infrared-quiet clumps, we selected 46 objects to follow up with the FLASH$^+$ receiver (Klein et al. 2014) at the APEX telescope (Güsten et al. 2006). The source selection was driven by the distance distribution of the sample to favour close-by objects ($d \le 4.5$ kpc), and focused on the newly discovered objects (see Sect. 4.2).

Using the APEX/FLASH$^+$ receiver we covered 8 GHz bandwidth simultaneously centred at 347.15 GHz in USB setup, from which we exploit here the SiO (8–7) transition at 347.33058 GHz. The observations were carried out in 2013 and 2014 August, and we obtained 90" × 90" OTF maps. The reference position was offset by 10' in right ascension and declination. Pointing and focus were regularly checked during the observations.

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** SiO (8–7) line observed towards the infrared-quiet massive clump G340.9698-1.0212. The dashed line shows the rest velocity of the source from the $C^{17}O$ (3–2) observations at 337.06110 GHz from the same dataset.</figcaption></figure>

The data have been reduced using the GILDAS software following standard procedures. The spectra at the centre, i.e. the dust peak position, have been extracted from the maps. Baselines of order zero have been subtracted and the data have been scaled to the $T_{mb}$ temperature using 73% efficiency$^4$. The average noise level is 48 mK in 1.3 km s$^{-1}$ resolution.

A total of 27 out of the 46 objects (59%) show emission of this line above $4\sigma$ as measured on the line area with Gaussian fitting. We show an example in Fig. 5, while all the detections are shown in Appendix C, Fig. C.1. The observed detection rate is similar to that of Klaassen et al. (2012) despite the lower sensitivity, and is somewhat lower than reported by Leurini et al. (2014), who had, on the other hand, a factor of ~2 higher sensitivity for the same transition. These studies focused on more evolved, mostly infrared-bright clumps, many of them hosting UC-H II regions. The fact that the high detection rates are broadly similar towards the less luminous thus likely younger clumps is consistent with the findings of Csengeri et al. (2016a) using lower-$J$ SiO lines.

# 3. Statistical properties of the brightest submillimetre sources

In the following we investigate the distribution and the statistics of the estimated physical properties of the entire sample. As defined in Sect. 2.6 we distinguish between infrared-quiet massive clumps and those classified as infrared-bright, and also show the subsample hosting MYSOs and/or UC-H II regions. Their average physical properties are summarised in Table 2.

## 3.1. Mass-distance distribution and completeness

In Fig. 6 we show the distance distribution of the sample. As discussed in Sect. 2.1, a substantial fraction of the sources are part of active star-forming complexes. Their typical distance is around 3–4 kpc, close to the average distance of 4.5 kpc for ATLASGAL clumps (Wienen et al. 2012). Overall 60% of the sample is located within 4.5 kpc; a total of 86 sources are located at greater distances, of which 29 are infrared quiet. The

<footnote>$^4$ http://www.apex-telescope.org/telescope/efficiency/</footnote>
<figure><img src="image_7.png" /></figure>

**Fig. 6.** Temperature corrected mass versus distance distribution. Dotted lines show the mass sensitivity corresponding to the flux density threshold of 5 Jy using $T_d = 30$ K for infrared-bright (red), and $T_d = 18$ K for infrared-quiet clumps (green). Coloured symbols are shown in the figure legend; filled circles correspond to sources where our selection is complete, while filled triangles show sources where the sample is incomplete above 650 $M_\odot$. The dashed horizontal line shows the 650 $M_\odot$ limit for massive clumps potentially capable of forming high-mass stars from Csengeri et al. (2014). Known (proto)clusters from the literature are labelled with magenta stars, as are the estimated initial clump mass at half-radius for the galactic (super) star clusters like the Orion Nebula Cluster (ONC, e.g. Da Rio et al. 2012) and Arches (Habibi et al. 2013), following Tan et al. (2014). The approximate distance limit for detecting sources corresponding to cores, clumps, and cloud fragments is shown as a dotted line (Bergin & Tafalla 2007).

more distant sources are therefore dominated by infrared-bright clumps. We also show a threshold of 650 $M_\odot$, introduced in Csengeri et al. (2014), which is an extrapolation of the massive dense core definition by Motte et al. (2007) to the ~0.3 pc scale of ATLASGAL clumps. Empirically this selection has proved to select cores hosting high-mass protostars by Bontemps et al. (2010).

The mass estimation corresponding to the threshold of 5 Jy used for our flux limit is indicated in Fig. 6 as green and red dotted lines for infrared-bright ($T_d = 30$ K) and infrared-quiet clumps ($T_d = 18$ K), respectively. It shows that our selection is complete for infrared-quiet clumps above 650 $M_\odot$ up to a distance of 4.5 kpc, while for infrared-bright clumps above 310 $M_\odot$. Since we expect the more evolved clumps to have lost a part of their mass in the central region through cluster formation and dissipation, we can argue that this sample is roughly complete for all clumps having 650 $M_\odot$ at birth at all their stages of evolution from infrared quiet to infrared bright and possibly also UC-H II regions. A total of 28 infrared-quiet and 55 infrared-bright clumps are found above 650 $M_\odot$ and 310 $M_\odot$ at a distance closer than 4.5 kpc.

The area covered by the sector of our Galaxy from $\ell = -60^\circ$ to $60^\circ$ and for a distance lower than 4.5 kpc is found to represent ~10% of the total area of the galactic disk within the solar circle after removing the inner 3 kpc. In a rough approximation, we can consider this 10% as the fraction of star formation in the Galactic disk which is occurring in our distance limited sample.

## 3.2. Temperature distribution

From the $\text{NH}_3$ data discussed in Sect. 2.2 we determine a median kinetic temperature of 25.5 K for 60% of the sample. Figure 7

<figure><img src="image_8.png" /></figure>

**Fig. 7.** Distribution of kinetic temperature from $\text{NH}_3$ (1, 1) and (2, 2) measurements from Wienen et al. (2012) and (2017). The infrared-bright and infrared-quiet sources are shown in different colours, as explained in the figure legend. We also show massive clumps hosting UC-H II regions as a subsample of the infrared-bright sources.

shows the temperature distribution of these sources. In this sample 28 have a firm association with UC-H II-regions, and the majority of the sources are classified as infrared bright, in total 71. These types of sources have a median kinetic temperature of 29.6 and 27.2 K, with a standard deviation of 8–11 K. Towards infrared-quiet sources we find a median of 23.8 K with a standard deviation of 4.7 K. This suggests a trend of lower temperatures towards younger sources; however, given their large standard deviation, this difference is not statistically significant. Following our argumentation in Sect. 2.4, and considering the similar temperatures, here we do not make a distinction between massive clumps with UC-H II regions and those hosting infrared-bright sources. A more accurate temperature assignment to the group of infrared-quiet massive clumps still remains uncertain since the sources with available $\text{NH}_3$ measurements are dominated by the more evolved sources, and the targeted infrared-quiet clumps may not be representative of the entire subsample.

## 3.3. Distribution of physical properties

In Fig. 8 we show the mass distribution of the sample, where we adopted $T_d = 30$ K for the infrared-bright sources and kept $T_d = 18$ K for the infrared-quiet clumps (see Sect. 2.3). We give an overview of the statistical properties of the physical parameters of all sources and the two source types in Table 2. The global statistical properties of the infrared-bright and -quiet samples are similar, and we find a remarkably flat mass distribution for all source types, the majority of the sources lying between 650 $M_\odot$ and 8000 $M_\odot$. Lower mass objects correspond to nearby massive dense cores with masses between ~86–120 $M_\odot$. Towards sources at greater distances we find masses of ~$10^5 M_\odot$ corresponding to cloud fragments rather than clumps. The mean clump mass is ~4.5 × $10^3 M_\odot$, which is about a factor of two higher compared to the star formation selected larger sample of ATLASGAL clumps studied by Urquhart et al. (2014).

The size distribution of the sample is shown in Fig. 9, and as above, there is no significant difference between the different source types. Adopting the nomenclature of Williams et al. (2000, see also Bergin & Tafalla 2007), we find a small fraction of sources that correspond to cores or entire cloud fragments,
<table><caption>**Table 2.** Statistics of the physical properties of all sources, and the infrared-bright and infrared-quiet samples.</caption><thead><tr><th></th><th colspan="5">All sources</th></tr><tr><th></th><th>Min</th><th>Max</th><th>Mean</th><th>Median</th><th>Stdev</th></tr></thead><tbody><tr><td>Distance [kpc]</td><td>1.00</td><td>15.6</td><td>4.9</td><td>3.9</td><td>3.30</td></tr><tr><td>R<sub>90</sub> [pc]</td><td>0.07</td><td>1.63</td><td>0.40</td><td>0.32</td><td>0.28</td></tr><tr><td>M<sub>18 K</sub> [M<sub>☉</sub>]</td><td>120</td><td>0.9 × 10<sup>5</sup></td><td>4.5 × 10<sup>3</sup></td><td>1.8 × 10<sup>3</sup></td><td>7.9 × 10<sup>3</sup></td></tr><tr><td>M<sub>corr</sub> [M<sub>☉</sub>]</td><td>58</td><td>4.2 × 10<sup>5</sup></td><td>2.8 × 10<sup>3</sup></td><td>1.3 × 10<sup>3</sup></td><td>4.6 × 10<sup>3</sup></td></tr><tr><td>Σ<sub>corr</sub> [g cm<sup>-2</sup>]</td><td>0.27</td><td>5.9</td><td>1.11</td><td>0.75</td><td>1.02</td></tr><tr><td>n&#773;<sub>corr</sub> [×10<sup>5</sup> cm<sup>-3</sup>]</td><td>0.14</td><td>44.15</td><td>3.17</td><td>1.56</td><td>5.25</td></tr><tr><td colspan="6">Infrared-bright massive clumps</td></tr><tr><td>Distance [kpc]</td><td>1.0</td><td>14.20</td><td>5.14</td><td>4.20</td><td>3.43</td></tr><tr><td>R<sub>90</sub></td><td>0.07</td><td>1.41</td><td>0.42</td><td>0.35</td><td>0.27</td></tr><tr><td>Mass<sub>18 K</sub> [M<sub>☉</sub>]</td><td>120</td><td>8.7 × 10<sup>4</sup></td><td>5.2 × 10<sup>3</sup></td><td>2.4 × 10<sup>3</sup></td><td>9.3 × 10<sup>3</sup></td></tr><tr><td>Mass<sub>corr</sub> [M<sub>☉</sub>]</td><td>58</td><td>4.2 × 10<sup>4</sup></td><td>2.7 × 10<sup>3</sup></td><td>1.2 × 10<sup>3</sup></td><td>4.7 × 10<sup>3</sup></td></tr><tr><td>Σ<sub>18 K</sub> [g cm<sup>-2</sup>]</td><td>0.27</td><td>6.00</td><td>0.99</td><td>0.56</td><td>1.14</td></tr><tr><td>n&#773;<sub>18 K</sub> [×10<sup>5</sup> cm<sup>-3</sup>]</td><td>0.14</td><td>44.15</td><td>2.78</td><td>1.03</td><td>5.89</td></tr><tr><td colspan="6">Infrared-quiet massive clumps</td></tr><tr><td>Distance [kpc]</td><td>1.15</td><td>15.60</td><td>4.47</td><td>3.70</td><td>3.05</td></tr><tr><td>R<sub>90</sub></td><td>0.08</td><td>1.63</td><td>0.36</td><td>0.28</td><td>0.28</td></tr><tr><td>Mass<sub>18 K</sub> [M<sub>☉</sub>]</td><td>200</td><td>2.1 × 10<sup>4</sup></td><td>3.1 × 10<sup>3</sup></td><td>1.4 × 10<sup>3</sup></td><td>4.5 × 10<sup>3</sup></td></tr><tr><td>Σ<sub>18 K</sub> [g cm<sup>-2</sup>]</td><td>0.51</td><td>4.05</td><td>1.31</td><td>1.03</td><td>0.74</td></tr><tr><td>n&#773;<sub>18 K</sub> [×10<sup>5</sup> cm<sup>-3</sup>]</td><td>0.20</td><td>18.46</td><td>3.83</td><td>2.39</td><td>3.91</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Temperature corrected mass distribution (using $T_d = 30$ K for the infrared-bright clumps). Colours are the same as in Fig. 7. The clumps hosting genuine UC-H II regions are a subsample of the infrared-bright clumps.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Distribution of the deconvolved physical size ($R_{90}$) of the sample. Colours are the same as in Fig. 7. The dashed lines indicate an approximate distinction between cores, clumps, and cloud fragments (Bergin & Tafalla 2007).</figcaption></figure>

while the majority of the sources correspond to clumps. They have typical sizes of 0.3–0.4 pc, which is comparable to the mean size of ~0.36 pc found by Russeil et al. (2010) in NGC 6334/NGC 6357. Compared to the population of massive dense cores in Cygnus-X, the sizes are ~3 times larger, while the masses are on average ~5 times larger, suggesting their potential to form high-mass stars.

In Fig. 10 we show the mass surface density distribution of the sample. Lower surface densities are found towards the infrared-bright sources, a considerable fraction of which host UC-H II regions. This is likely due to an observational bias, driven both by a larger fraction of such sources at a larger distance, which then correspond to larger structures with lower average surface density, and the lower mass selection threshold for

the infrared-bright sources. The mean surface density is above the theoretical value of $\Sigma_{cl} \ge 1 \text{ g cm}^{-2}$ by Krumholz & McKee (2008), suggesting that a substantial fraction of our sources have the potential to form high-mass stars. Correspondingly, the median volume-averaged density of $3.13 \times 10^5 \text{ cm}^{-3}$ is higher than that found in other samples in the literature (e.g. Russeil et al. 2010), and suggests that we probe here high density, star-forming gas.

## 4. Nature of the brightest submillimetre sources in the inner Galaxy

The brightest submillimetre sources in the inner Galaxy are generally associated with active sites of star formation. Here we
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Mass surface density (Σ) distribution of the sample using $T_d = 30$ K for the infrared-bright clumps. Colours and symbols are the same as in Fig. 7.</figcaption></figure>

first investigate the star formation potential and evolutionary aspect of the entire sample (Sect. 4.1), and argue that the majority of the infrared-quiet massive clumps are potential precursors of massive stars and/or clusters. We then show that our sample of infrared-quiet clumps contains sources previously not recognised in other surveys using infrared emission, or absorption (Sect. 4.2). Finally, we search for signposts of embedded protostars in the infrared-quiet subsample to assess the stage of ongoing star formation (Sect. 4.3). We then discuss the stability of the entire sample (Sect. 4.4), and provide predictions for their small-scale fragmentation (Sect. 4.5).

## 4.1. Precursors of massive clusters

Comparing the physical properties in terms of mass, size, and density, in Table 2 we show that our sample exhibits rather homogeneous properties, and no statistically significant differences are found between the two source types. This includes well-studied Galactic mini-starburst regions, for instance the W43 ridge (e.g. Motte et al. 2003; Bally et al. 2010; Nguyen Luong et al. 2011), W51 Main (e.g. Mufson & Liszt 1979; Zapata et al. 2008; Ginsburg et al. 2015), and W49N (e.g. Peng et al. 2010; Galván-Madrid et al. 2013). While indications of ongoing star formation is prevalent in the sample, are all clumps prone to form mini-starbursts and rich clusters? We compare our results with the empirical relation of Kauffmann & Pillai (2010) in Fig. 11, which suggests that all of our sources have the potential to form high-mass stars. While the comparative statistical analysis of infrared-quiet and -bright sources reveals similar physical characteristics in terms of mass and size, they seem to differ only in their embedded (proto-)stellar population. Low bolometric luminosity sources with similar physical characteristics to their more evolved counterparts have been frequently interpreted as being in an earlier evolutionary stage (e.g. Motte et al. 2007; Molinari et al. 2008). This suggests that infrared-quiet clumps are the precursors of the infrared-bright clumps and UC-H II regions, and the sample can be arranged along a similar evolutionary path.

The large mass reservoir of infrared-bright sources suggests that they likely continue producing high-mass stars. Unless specific conditions are required for the formation of rich clusters, we shall have selected here their true precursors.

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Temperature corrected mass estimates (using $T_d = 30$ K for the infrared-bright clumps) versus deconvolved physical size ($R_{90}$) from Sect. 2.3. Dashed lines show the relation of Kauffmann et al. (2010) scaled to the Ossenkopf & Henning (1994) dust opacity following Dunham et al. (2011). In the upper left corner we show the factor of two uncertainty in the mass estimate. Colours are the same as in Fig. 7.</figcaption></figure>

## 4.2. Massive cluster progenitors in isolation?

Because the clustered nature of high-mass star formation was recognised early on, several studies searched for younger objects in the vicinity of bright, luminous sources (e.g. Garay et al. 2004; Thompson et al. 2005). To look for neighbouring, more evolved objects with UC-H II regions in the close vicinity of the infrared-quiet sample, we used the CS (2–1) survey of IRAS sources (Bronfman et al. 1996), which is complete towards IRAS selected UC-H II regions in the Galactic plane. In total we found 92 associations (42%) within 0.8' offset, which corresponds to 0.26 pc for the nearest source and 2.9 pc for the most distant source in our sample. However, in the infrared-quiet sample only 29 (36%) are found in the vicinity of an IRAS source, and their offsets correspond to a physical separation between 0.3–16.3 pc. Compared with the more sensitive and higher angular resolution Red MSX survey (RMS, Lumsden et al. 2013), we find six infrared-quiet clumps in the vicinity of a bright MYSO with $L_{bol} > 2 \times 10^4 L_{\odot}$.

Surveys sensitive to cold dust are better suited to directly probe the youngest clumps. To evaluate the fraction of previously identified infrared-quiet clumps, we first looked for associations with cloud structures identified as infrared-dark clouds (IRDCs) by correlating the peak dust positions with the catalogues of IRDCs revealed by MSX (Jackson et al. 2008; Simon et al. 2006) and *Spitzer* (Peretto & Fuller 2009). In total we find 20 associations with an MSX Dark Cloud within 120", 15 of which are infrared-quiet sources. For the association with *Spitzer* Dark Clouds we used the Herschel-confirmed clouds from Peretto et al. (2016) with a smaller positional offset of 50" corresponding to their average equivalent size of ~25" + 1σ, where σ is their standard deviation. In total, 65 of our clumps have associations, 31 of which are infrared-quiet clumps. This means that ~61% of the infrared-quiet clumps have not been recognised before in these surveys.

Thus, we find that the infrared-quiet sample is found in relative isolation compared to the bright clumps which already undergo active star formation. Less than half of the ATLASGAL selected massive clumps have IRDC counterparts, and comparison with the star formation selected ATLASGAL clumps of
Urquhart et al. (2014) shows that 25 of the infrared-quiet sources are newly identified. In summary, we find that ~23% of our infrared-quiet sample has not been recognised by these previous studies.

### 4.3. Signposts of already formed embedded high-mass protostars

To tackle the stage of ongoing star formation, here we compare our sample with tracers of high-mass star formation. We found a total of 23 clumps coinciding with extended green objects (EGOs); (Cyganowski et al. 2008; and Chen et al. 2013) most of them are infrared-quiet clumps, 2 sources are infrared-bright, and 3 are UC-H II regions. Since EGOs are believed to trace powerful outflows mostly from (presumably massive) Class 0 YSOs, it confirms that the infrared-quiet clumps host the youngest protostellar phases. However, only 25% of the infrared-quiet clumps are associated with EGOs. This may indicate that the EGO phase is shorter than the infrared-quiet clump phase, and therefore does not represent well the entire early evolutionary phase of massive clumps. The few EGOs associated with infrared-bright clumps may indicate that young massive protostars can still form at the time of the infrared-bright phase, and even perhaps at the UC-H II phase if star formation is a continuous process in clumps.

On the other hand, the lack of EGOs does not necessarily imply a lack of outflow activity, since high extinction, even at 4.5 µm, may prevent their detection. In our SiO (8–7) spectroscopic study, ~60% of infrared-quiet clumps show emission, which is a much larger fraction than that of EGOs. Using the lower energy SiO (2–1) line, Csengeri et al. (2016a) find that an even higher fraction of the infrared-quiet massive clumps exhibit emission in this shocked gas tracer. This is consistent with a picture that infall and accretion processes may already be occurring in at least some of the infrared-quiet clumps, indirectly pointing at an underlying population of deeply embedded protostars (see the example of NGC 6334 I(N) by Megeath & Tieftrunk 1999).

Another tracer of the nature of embedded protostars is Class II methanol maser emission. These masers are excited by radiative pumping requiring the intense mid-infrared emission from the warm, dense material of deeply embedded hot-cores or UC-H II regions (Menten 1991; Sobolev et al. 1997). We therefore cross-correlated our sample with the Methanol MultiBeam survey (Caswell et al. 2010, 2011; Green et al. 2010; Breen et al. 2015), and found 138 (66%) associations within 20", of which 57 are infrared-quiet massive clumps corresponding to ~71% of the subsample. Methanol masers thus seem to trace a larger fraction of the infrared-quiet phase, although less selectively since they are also found towards the more evolved stages of clump evolution. Altogether we find that a considerable fraction of the infrared-quiet sample show signs of embedded (high-mass) protostars and ongoing star formation process.

### 4.4. Stability estimates and timescales

Using high-density tracers we estimate here the virial mass and compare it to the mass estimates based on the dust emission. As a first step we used the NH₃ (1, 1) line width measurements of Wienen et al. (2012) and Wienen et al. (2017) where available, and then complemented these data with measurements of the CS (2–1) line from Bronfman et al. (1996), and finally of the H¹³CO⁺ (1–0) line from recent surveys using the Mostra (Wyrowski et al., in prep.) and the IRAM 30 m telescopes (Csengeri et al. 2016a). We have line width estimates for a total

<figure><img src="image_13.png" /><figcaption>Fig. 12. Temperature corrected clump mass estimates compared to the virial mass. The dashed line indicates where the virial mass equals the dust mass. Colours are the same as in Fig. 7, filled circles correspond to line width estimates from Wienen et al. (2012, 2015), filled upward triangle to that of Csengeri et al. (2016a), filled downward triangle to that of Wyrowski et al. (in prep.), filled diamonds to that of Bronfman et al. (1996). Sources with d > 5 kpc are plotted in smaller symbols.</figcaption></figure>

of 173 sources, 80% of the whole sample. While this is a smaller sample than studied in Wienen et al. (2012) and Urquhart et al. (2014), it is more focused on the brightest ATLASGAL clumps.

We estimate the virial mass assuming a density profile of $n(r) \sim r^{-2}$, and using $M_{Vir} = 3 \times \frac{R_{90} \times \sigma_{tot}^2}{G}$ (Bertoldi & McKee 1992), where $R_{90}$ is the radius of a sphere, $\sigma_{tot}$ is the velocity dispersion of the gas, and G is the gravitational constant,

$$ \frac{M_{\text{Vir}}}{[M_{\odot}]} = 697 \frac{R_{90}}{[\text{pc}]} \frac{\sigma_{\text{tot}}^2}{[\text{km s}^{-1}]} \quad (5) $$

As shown in Fig. 12, we find surprisingly small $M_{Vir}/M$ ratios towards the majority of the sources, on average 0.47 with a median of 0.29, strongly suggesting that they are gravitationally unstable. While the used molecular tracers selectively probe the high-density gas, uncertainties due to the different tracers and beam sizes could affect our virial mass estimates. However, despite the various datasets used, we see no systematics in the estimated virial masses. As Fig. 12 shows, 155 sources (89%) have a smaller virial mass than their estimated gas mass. We find the same trend for the more distant sources of the sample, which are beyond our completeness limit. Uncertainties in the mass estimates due to the dust opacity and gas-to-dust ratio would lower the $M_{Vir}/M$ ratio by a factor of two. Even considering the case of a uniform $n(r) \sim r^0$ density profile, we would still find a large fraction (67%) of the clumps having $M_{Vir}/M < 1$, which could thus still be considered as gravitationally unstable. Magnetic fields could, however, slow down the rate of collapse. For a large number of clouds Kauffmann et al. (2013) provides estimates of the virial parameter defined by $M_{Vir}/M$, and find similarly low values. They suggest that magnetic fields as strong as 1 mG would be required to slow down the collapse, otherwise they undergo rapid collapse.

The size-scale of clumps, 0.2–1 pc, is large compared to the scale of MDCs (<0.15 pc) and cores (<0.05 pc) for which collapse signatures are frequently observed. Our results therefore suggest that these massive cluster progenitors are already collapsing on a global scale, i.e. a clump scale. This is particularly intriguing because it suggests that gravity could be the
<table><caption>**Table 3.** Comparison of the properties of the infrared-quiet sources with the literature.</caption><thead><tr><th></th><th>ATLASGAL infrared-quiet<sup>†</sup><br>clumps</th><th>Herschel-IRDC<sup>a</sup><br>clumps</th><th>IRDC<br>clumps</th><th>CygX<br>clumps</th></tr></thead><tbody><tr><td>Distance [kpc]</td><td>1.15–15.60 (4.47)</td><td>0.12–16.45 (4.18)</td><td>1.8–7.1 (3.7)<sup>b</sup></td><td>1.4<sup>d</sup></td></tr><tr><td>Size [pc]</td><td>0.08–1.63 (0.36)</td><td>0.02–2.6 (0.6)</td><td>1–3<sup>b</sup></td><td>0.56<sup>e</sup></td></tr><tr><td><i>M</i> [M<sub>☉</sub>]<sup>‡</sup></td><td>200–21 000 (1400)</td><td>0.11–16 000 (700)</td><td>120–16 000 (950)<sup>b</sup></td><td>680<sup>e</sup></td></tr><tr><td><i>n̄</i><sub>corr</sub> [×10<sup>5</sup> cm<sup>−3</sup>]</td><td>0.20–18.46 (2.39)</td><td>0–0.58 (0.01)</td><td>0.001–0.1<sup>b</sup></td><td>0.12<sup>e</sup></td></tr></tbody></table>

**Notes.** (†) Range of minimum to maximum values, the median is given between parentheses. (‡) For the adopted dust temperature see their original work. (a) Adopted from the values listed in Table 2 from Traficante et al. (2015). (b) Adopted from the values listed in Table 3 from Rathborne et al. (2006). (c) Adopted from the values listed in Table 1 from Rathborne et al. (2006). (d) Distance from Rygl et al. (2012). (e) Scaled to 1.4 kpc from Motte et al. (2007).

dominant force in the mass assembly process over size-scales of entire clumps, which are larger than previously considered (e.g. McKee & Ostriker 2007).

We also calculate the free-fall time of these clumps by $\tau_{ff} = \sqrt{3\pi/32Gn}$, where G is the gravitational constant, and n is the volume density. We find on average $5.64 \times 10^4$ yr, which is somewhat shorter than the lifetime estimate of $7.5 \times 10^4$ yr by Csengeri et al. (2014). Such relative collapse timescales (i.e. a few times the free-fall timescale) are typically observed towards low-mass cores (see Ward-Thompson et al. 2007, and references therein), and have recently been suggested for the formation of high-mass cores by Duarte-Cabral et al. (2013).

Collapse timescales of a few times $10^4$ yr suggests that all massive clumps, precursors of massive clusters, are transient entities and should disappear on short timescales, typically less than $10^5$ yr (see also Zinnecker & Yorke 2007; Motte et al. 2007). Embedded MYSOs and UC-H II evolve, however, on longer timescales, of the order of $10^5–10^6$ yr according to stellar evolution models with accretion rates of $M_{acc} \sim 10^{-3} M_{\odot} \text{yr}^{-1}$ (e.g. Churchwell 2002; Hoare et al. 2007; Kuiper & Yorke 2013). This is at least up to an order of magnitude longer than estimated for the collapse of their parent clump. This means that the clump material should have already collapsed leaving a considerably less dense and less massive clump towards the more evolved objects. This apparent contradiction to the observations can only be explained if there is mass replenishment onto the clump on large scales fuelling a continuous star formation, or if we selectively pick here the youngest of the galactic UC-H II regions still surrounded by their natal clump. Case studies of massive cluster progenitors show evidence both for global collapse and for flows of dense gas (e.g. Schneider et al. 2010; Csengeri et al. 2011a,b; Peretto et al. 2013; Tackenberg et al. 2014). This points to a scenario where massive clumps undergo star formation, while gas replenishment arrives on larger scales, and maintains the observed similar mass range for these clumps.

### 4.5. Initial stages of massive proto-clusters: prediction for fragmentation properties

We overview the physical properties of our infrared-quiet sample in comparison with typical IRDCs (Rathborne et al. 2006), Herschel clumps embedded in IRDCs (Traficante et al. 2015), and clumps in the Cygnus-X high-mass star-forming region (Motte et al. 2007) in Table 3. While IRDCs are typically more extended structures, they have lower masses and exhibit smaller volume densities. Herschel selected clumps span a much broader range of physical conditions and cover similarly high-mass clumps as our sample. They are found to be, however, more

extended and are therefore on average lower density structures. The ATLASGAL counterparts of IRDCs selected by Kainulainen & Tan (2013) correspond to sources below the flux density threshold used for this study. The ATLASGAL sources used here are the brightest and therefore the most massive among the compact millimetre sources in the inner Galaxy. Their typical size of ~0.4 pc corresponds to protoclusters. We can therefore discuss our sample in terms of the precursors of the richest clusters in the Galaxy.

In Fig. 13 we show the temperature corrected mass versus the surface density distribution of the entire sample following Tan et al. (2013), with different colours corresponding to the infrared-bright and infrared-quiet subsamples. It is clear that the sample discussed here lies, on average, close to the theoretical limit for high-mass star formation around $\Sigma \approx 1 \text{ g cm}^{-2}$ (Krumholz & McKee 2008) supporting the scenario that it is capable of forming high-mass stars and clusters. Taking the two extreme profiles of constant density $n(r) \sim r^0$ and power law with $n(r) \sim r^{-2}$, we mark with stars the predicted fragmentation on scales smaller by a factor of ten, approaching the size-scales of a few thousand AU corresponding to individual collapsing objects (e.g. Bontemps et al. 2010). The uniform density profile, where the surface density stays constant, predicts fragments below solar mass, which would correspond to clumps forming a rich cluster of low-mass stars. The case of the steeper density profile predicts more massive structures on smaller scales, reaching a mass range of several tens of solar masses. Envelope masses of this order are expected to form high-mass stars.

Some of the well-studied sources of our selection have been subject to high angular resolution follow-ups. The extreme clumps of the W43 complex, MM1 to MM3 (Motte et al. 2003; Sridharan et al. 2014) hosts one of the most extreme sites of active high-mass star formation in terms of mass and youth (Louvet et al. 2014). Massive envelopes have been identified towards other clumps as well, for example G35.03+0.35 (G35.0250+0.3501) (Beltrán et al. 2014), G35.20-0.74 (G35.1976-0.7427) (Sánchez-Monge et al. 2013), and SDC335-MM1(G335.5857-0.2906) (Peretto et al. 2013). High angular resolution observations with ALMA towards a large fraction of the infrared-quiet sources of the presented sample indeed reveal a significant population of massive dense cores (Csengeri et al. 2017).

### 5. Summary and conclusions

We investigate the properties of a flux limited sample of massive clumps from the ATLASGAL survey. Based on mid-infrared photometry, we divide the sample into infrared-quiet and infrared-bright clumps. Infrared-bright clumps peak on
<table><caption>**Table 4.** Flux limited sample of ATLASGAL sources with distance estimate and the corresponding dust properties.</caption><thead><tr><th>ATLASGAL<sup>a</sup><br>source name</th><th>Gal. long.<br>[°]</th><th>Gal. lat.<br>[°]</th><th>F<sub>r</sub><br>[Jy/beam]</th><th>S<sub>r</sub><br>[Jy]</th><th>Θ<sub>A</sub>×Θ<sub>B</sub><br>["×"]</th><th>PA<br>[°]</th><th>d<sup>p</sup><br>[kpc]</th><th>d<sub>min</sub><br>[kpc]</th><th>d<sub>max</sub><br>[kpc]</th><th>Ref.</th><th>M<sub>con</sub><br>[M<sub>6</sub>]</th><th>R<sub>90</sub><br>[pc]</th><th>Σ<br>[g·cm<sup>-2</sup>]</th><th>Π<br>[x10<sup>4</sup>cm<sup>-3</sup>]</th><th>IRq</th><th>GMC/SIRDC</th><th>EGO</th></tr></thead><tbody><tr><td>G300/9688+1.1456</td><td>300.9688</td><td>1.1456</td><td>6.5</td><td>16.2</td><td>32.1×28.4</td><td>125</td><td>4.3</td><td>1.6</td><td>7.0</td><td>(1)</td><td>920</td><td>0.4</td><td>0.4</td><td>5.3</td><td>n</td><td>–</td><td>–</td></tr><tr><td>G301/1365-0.2256*</td><td>301.1365</td><td>-0.2256</td><td>21.4</td><td>34.5</td><td>26.3×22.6</td><td>102</td><td>4.2</td><td>3.5</td><td>10.4</td><td>(1)</td><td>1877</td><td>0.3</td><td>1.8</td><td>43.6</td><td>n</td><td>301.147-0.228</td><td>–</td></tr><tr><td>G305/2083+0.2063*</td><td>305.2083</td><td>0.2063</td><td>13.5</td><td>26.1</td><td>29.7×24.0</td><td>15</td><td>3.5</td><td>1.5</td><td>6.6</td><td>(17)</td><td>985</td><td>0.3</td><td>0.9</td><td>21.0</td><td>n</td><td>G305/SIRDC 305.209+0.205</td><td>–</td></tr><tr><td>G305/2703-0.0090</td><td>305.2703</td><td>-0.0090</td><td>6.7</td><td>18.8</td><td>37.3×27.9</td><td>126</td><td>3.5</td><td>1.3</td><td>4.9</td><td>(17)</td><td>710</td><td>0.4</td><td>0.3</td><td>5.6</td><td>n</td><td>G305/SIRDC 305.270+0.018</td><td>–</td></tr><tr><td>G305/3570+0.2033</td><td>305.3570</td><td>0.2033</td><td>7.6</td><td>20.5</td><td>35.0×28.4</td><td>66</td><td>3.5</td><td>1.5</td><td>6.6</td><td>(17)</td><td>774</td><td>0.4</td><td>0.4</td><td>6.7</td><td>n</td><td>G305</td><td>–</td></tr><tr><td>G305/3617+0.1514</td><td>305.3617</td><td>0.1514</td><td>5.4</td><td>9.6</td><td>29.1×22.3</td><td>80</td><td>3.5</td><td>1.5</td><td>6.6</td><td>(361)</td><td>0.2</td><td>0.2</td><td>0.4</td><td>10.5</td><td>n</td><td>G305</td><td>–</td></tr><tr><td>G305/3666+0.2130</td><td>305.3666</td><td>0.2130</td><td>5.7</td><td>10.2</td><td>26.5×24.8</td><td>80</td><td>3.5</td><td>1.3</td><td>6.6</td><td>(17)</td><td>792</td><td>0.2</td><td>0.9</td><td>22.0</td><td>y</td><td>G305/SIRDC 305.364+0.215</td><td>–</td></tr><tr><td>G305/7994+0.2445</td><td>305.7994</td><td>-0.2445</td><td>11.7</td><td>16.1</td><td>23.6×21.6</td><td>92</td><td>3.5</td><td>1.5</td><td>12.0</td><td>(17)</td><td>1251</td><td>0.2</td><td>2.8</td><td>101.6</td><td>y</td><td>G305</td><td>–</td></tr><tr><td>G307/5610-0.5873</td><td>307.5610</td><td>-0.5873</td><td>5.8</td><td>11.9</td><td>29.7×25.6</td><td>105</td><td>3.5</td><td>1.3</td><td>6.6</td><td>(17)</td><td>448</td><td>0.3</td><td>0.4</td><td>7.9</td><td>n</td><td>G305</td><td>–</td></tr><tr><td>G309/9212+0.4779</td><td>309.9212</td><td>0.4779</td><td>6.8</td><td>14.1</td><td>29.7×25.7</td><td>89</td><td>4.3</td><td>2.3</td><td>5.4</td><td>(1)</td><td>802</td><td>0.4</td><td>0.4</td><td>7.5</td><td>n</td><td>G305</td><td>–</td></tr><tr><td>G311/6261+0.2901</td><td>311.6261</td><td>0.2901</td><td>7.0</td><td>13.1</td><td>28.9×23.8</td><td>96</td><td>5.5</td><td>3.5</td><td>8.0</td><td>(1)</td><td>1220</td><td>0.4</td><td>0.5</td><td>7.5</td><td>n</td><td>Ws69</td><td>–</td></tr><tr><td>G314/2198+0.2707</td><td>314.2198</td><td>0.2707</td><td>6.7</td><td>15.9</td><td>35.9×24.6</td><td>106</td><td>7.8</td><td>5.8</td><td>9.3</td><td>(1)</td><td>2987</td><td>0.4</td><td>0.4</td><td>3.1</td><td>n</td><td>Ws88</td><td>–</td></tr><tr><td>G316/7985+0.0554</td><td>316.7985</td><td>-0.0554</td><td>6.7</td><td>21.0</td><td>47.6×24.1</td><td>49</td><td>2.5</td><td>1.5</td><td>6.6</td><td>(1)</td><td>404</td><td>0.3</td><td>0.3</td><td>7.0</td><td>n</td><td>G317</td><td>–</td></tr><tr><td>G316/8131-0.0589</td><td>316.8131</td><td>-0.0589</td><td>5.6</td><td>17.1</td><td>39.8×28.1</td><td>13</td><td>2.5</td><td>1.1</td><td>4.0</td><td>(1)</td><td>329</td><td>0.3</td><td>0.3</td><td>6.0</td><td>n</td><td>G317/SIRDC 316.821+0.062</td><td>–</td></tr><tr><td>G318/0494+0.0861</td><td>318.0494</td><td>0.0861</td><td>7.0</td><td>12.9</td><td>26.7×25.5</td><td>-5</td><td>9.6</td><td>4.0</td><td>9.6</td><td>(1)</td><td>3658</td><td>0.7</td><td>0.5</td><td>4.4</td><td>n</td><td>SIRDC 318.038+0.086</td><td>–</td></tr><tr><td>G318/9471-0.1969</td><td>318.9471</td><td>-0.1969</td><td>7.0</td><td>12.6</td><td>27.1×24.7</td><td>55</td><td>10.7</td><td>4.0</td><td>10.7</td><td>(1)</td><td>4454</td><td>0.8</td><td>0.5</td><td>4.1</td><td>n</td><td>SIRDC 318.957+0.199</td><td>–</td></tr><tr><td>G320/2325-0.2844</td><td>320.2325</td><td>-0.2844</td><td>6.3</td><td>12.1</td><td>28.9×24.6</td><td>2</td><td>3.9</td><td>2.7</td><td>5.4</td><td>(1)</td><td>1162</td><td>0.3</td><td>0.9</td><td>18.1</td><td>y</td><td>Ws120/SIRDC 320.243+0.293</td><td>–</td></tr><tr><td>G321/9348+0.0066</td><td>321.9348</td><td>-0.0066</td><td>6.3</td><td>11.8</td><td>27.1×25.4</td><td>92</td><td>1.8</td><td>0.6</td><td>3.0</td><td>(1)</td><td>242</td><td>0.1</td><td>0.9</td><td>42.5</td><td>y</td><td>SIRDC 321.936+0.003</td><td>G321.94+0.01</td></tr><tr><td>G322/1581+0.6354</td><td>322.1581</td><td>0.6354</td><td>15.1</td><td>30.3</td><td>29.9×24.7</td><td>143</td><td>3.2</td><td>2.1</td><td>4.3</td><td>(1)</td><td>936</td><td>0.3</td><td>1.0</td><td>24.0</td><td>n</td><td>RCW92/SIRDC 322.149+0.639</td><td>–</td></tr><tr><td>G322/1632+0.6221</td><td>322.1632</td><td>0.6221</td><td>9.1</td><td>31.1</td><td>32.8×25.9</td><td>128</td><td>3.2</td><td>2.1</td><td>4.3</td><td>(1)</td><td>664</td><td>0.5</td><td>0.5</td><td>11.2</td><td>n</td><td>RCW92/SIRDC 322.162+0.615</td><td>–</td></tr><tr><td>G322/9314+1.3933</td><td>322.9314</td><td>1.3933</td><td>15.8</td><td>24.2</td><td>25.5×22.1</td><td>118</td><td>2.0</td><td>1.0</td><td>3.0</td><td>(1)</td><td>613</td><td>0.1</td><td>3.1</td><td>164.5</td><td>y</td><td>–</td><td>–</td></tr><tr><td>G323/4589+0.0800</td><td>323.4589</td><td>-0.0800</td><td>5.4</td><td>13.3</td><td>32.6×27.8</td><td>145</td><td>3.8</td><td>2.7</td><td>4.9</td><td>(1)</td><td>593</td><td>0.4</td><td>0.3</td><td>5.1</td><td>n</td><td>SIRDC 323.471+0.081</td><td>–</td></tr><tr><td>G323/7407-0.2635</td><td>323.7407</td><td>-0.2635</td><td>9.6</td><td>15.8</td><td>27.6×21.9</td><td>124</td><td>2.8</td><td>1.7</td><td>3.9</td><td>(1)</td><td>783</td><td>0.2</td><td>1.6</td><td>57.5</td><td>y</td><td>SIRDC 323.738+0.270</td><td>–</td></tr><tr><td>G324/2006+0.1208</td><td>324.2006</td><td>0.1208</td><td>9.6</td><td>17.1</td><td>28.7×23.0</td><td>78</td><td>8.7</td><td>6.9</td><td>9.9</td><td>(1)</td><td>4000</td><td>0.6</td><td>0.7</td><td>7.1</td><td>n</td><td>–</td><td>–</td></tr><tr><td>G326/4745+0.7027</td><td>326.4745</td><td>0.7027</td><td>12.8</td><td>24.4</td><td>28.4×24.8</td><td>80</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>964</td><td>0.2</td><td>1.8</td><td>58.6</td><td>y</td><td>RCW95/SIRDC 326.476+0.706</td><td>–</td></tr><tr><td>G326/6411+0.6127</td><td>326.6411</td><td>0.6127</td><td>8.0</td><td>18.0</td><td>29.9×27.6</td><td>44</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>710</td><td>0.2</td><td>1.0</td><td>27.2</td><td>y</td><td>RCW95</td><td>–</td></tr><tr><td>G326/6577+0.5941*</td><td>326.6577</td><td>0.5941</td><td>6.1</td><td>26.2</td><td>44.5×35.8</td><td>36</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>503</td><td>0.4</td><td>0.3</td><td>4.4</td><td>n</td><td>RCW95</td><td>–</td></tr><tr><td>G326/6706+0.5539</td><td>326.6706</td><td>0.5539</td><td>5.9</td><td>10.0</td><td>27.4×22.8</td><td>130</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>396</td><td>0.2</td><td>1.0</td><td>36.2</td><td>y</td><td>RCW95/SIRDC 326.678+0.549</td><td>–</td></tr><tr><td>G326/7236+0.6148</td><td>326.7236</td><td>0.6148</td><td>7.6</td><td>15.2</td><td>29.3×25.1</td><td>105</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>292</td><td>0.2</td><td>0.5</td><td>15.6</td><td>n</td><td>RCW95/SIRDC 326.726+0.598</td><td>–</td></tr><tr><td>G327/2921+0.5792*</td><td>327.2921</td><td>-0.5792</td><td>49.2</td><td>71.2</td><td>24.8×21.5</td><td>0</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>2815</td><td>0.1</td><td>10.6</td><td>498.4</td><td>y</td><td>G327</td><td>–</td></tr><tr><td>G327/3005-0.5509</td><td>327.3005</td><td>-0.5509</td><td>11.3</td><td>32.3</td><td>37.6×28.1</td><td>63</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>621</td><td>0.3</td><td>0.6</td><td>12.9</td><td>n</td><td>G327/SIRDC 327.306+0.566</td><td>–</td></tr><tr><td>G327/309+0.5775</td><td>327.309</td><td>-0.5775</td><td>7.1</td><td>21.3</td><td>50.0×22.0</td><td>147</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>843</td><td>0.3</td><td>0.7</td><td>15.9</td><td>y</td><td>SIRDC 327.397+0.439</td><td>G327.40+0.44</td></tr><tr><td>G327/4024+0.4436</td><td>327.4024</td><td>0.4436</td><td>10.4</td><td>18.6</td><td>26.7×24.7</td><td>102</td><td>5.0</td><td>3.5</td><td>6.0</td><td>(1)</td><td>2943</td><td>0.4</td><td>1.6</td><td>27.7</td><td>y</td><td>MSDCX G328.25+0.05.1</td><td>–</td></tr><tr><td>G328/2553-0.5481</td><td>328.2553</td><td>-0.5481</td><td>8.0</td><td>19.9</td><td>34.3×26.7</td><td>143</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>788</td><td>0.2</td><td>0.9</td><td>23.0</td><td>y</td><td>MSDCX G328.25+0.05.1</td><td>–</td></tr><tr><td>G328/2551-0.5321</td><td>328.2551</td><td>-0.5321</td><td>8.3</td><td>15.0</td><td>28.3×23.4</td><td>-31</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>591</td><td>0.2</td><td>1.3</td><td>43.9</td><td>y</td><td>MSDCX G328.25+0.05.1</td><td>–</td></tr><tr><td>G328/3073+0.4301</td><td>328.3073</td><td>0.4301</td><td>9.3</td><td>16.7</td><td>28.6×23.2</td><td>79</td><td>10.0</td><td>8.2</td><td>10.8</td><td>(1)</td><td>5140</td><td>0.7</td><td>0.7</td><td>5.9</td><td>n</td><td>–</td><td>–</td></tr><tr><td>G328/5670+0.5344</td><td>328.5670</td><td>-0.5344</td><td>11.3</td><td>24.4</td><td>33.6×23.7</td><td>-21</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>469</td><td>0.2</td><td>0.7</td><td>19.9</td><td>n</td><td>RCW99</td><td>–</td></tr><tr><td>G328/5763+0.5303</td><td>328.5763</td><td>-0.5303</td><td>5.4</td><td>10.1</td><td>28.6×24.1</td><td>23</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>194</td><td>0.2</td><td>0.4</td><td>12.6</td><td>n</td><td>RCW99</td><td>–</td></tr><tr><td>G328/8087+0.6328*</td><td>328.8087</td><td>0.6328</td><td>23.2</td><td>40.6</td><td>26.8×24.1</td><td>138</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>1607</td><td>0.3</td><td>3.6</td><td>130.1</td><td>n</td><td>SIRDC 328.798+0.625</td><td>–</td></tr><tr><td>G329/0303+0.2022*</td><td>329.0303</td><td>-0.2022</td><td>8.8</td><td>29.4</td><td>51.0×24.1</td><td>101</td><td>2.5</td><td>2.0</td><td>4.2</td><td>(1)</td><td>1162</td><td>0.3</td><td>0.8</td><td>17.2</td><td>y</td><td>SIRDC 329.028+0.202</td><td>–</td></tr><tr><td>G329/1835-0.3147</td><td>329.1835</td><td>-0.3147</td><td>6.4</td><td>11.3</td><td>28.5×22.7</td><td>86</td><td>4.2</td><td>3.5</td><td>5.0</td><td>(1)</td><td>1257</td><td>0.3</td><td>1.0</td><td>21.3</td><td>y</td><td>RCW99</td><td>–</td></tr><tr><td>G329/3382+0.1471</td><td>329.3382</td><td>0.1471</td><td>8.1</td><td>17.8</td><td>30.5×26.6</td><td>15</td><td>5.7</td><td>4.7</td><td>8.8</td><td>(1)</td><td>1775</td><td>0.5</td><td>0.5</td><td>6.0</td><td>n</td><td>–</td><td>–</td></tr><tr><td>G329/4058+0.4589</td><td>329.4058</td><td>-0.4589</td><td>6.4</td><td>14.8</td><td>29.2×29.2</td><td>97</td><td>10.5</td><td>9.6</td><td>11.4</td><td>(1)</td><td>5029</td><td>1.0</td><td>0.4</td><td>2.4</td><td>n</td><td>–</td><td>–</td></tr><tr><td>G330/2948+0.3945</td><td>330.2948</td><td>-0.3945</td><td>5.8</td><td>10.9</td><td>26.6×25.9</td><td>29</td><td>10.2</td><td>9.3</td><td>11.1</td><td>(1)</td><td>3494</td><td>0.8</td><td>0.4</td><td>3.4</td><td>n</td><td>–</td><td>–</td></tr><tr><td>G330/8788-0.3681*</td><td>330.8788</td><td>-0.3681</td><td>16.9</td><td>28.1</td><td>26.1×23.4</td><td>82</td><td>4.2</td><td>3.5</td><td>5.0</td><td>(1)</td><td>3132</td><td>0.3</td><td>2.8</td><td>65.5</td><td>y</td><td>–</td><td>–</td></tr></tbody></table>

Notes. The full table is available only in electronic form at the CDS. Column 1 gives the corresponding ATLASGAL source name from Csengeri et al. (2014) including the Galactic coordinates, while Cols. 2 and 3 list the Galactic coordinates. Columns 4, 5 give the peak flux density, and the integrated flux density assuming a Gaussian flux distribution from Csengeri et al. (2014), respectively. Columns 6–7 list the major and minor axis of the fitted Gaussian FWHM to the source, and the position angle. Column 8 gives the used distance estimate. Columns 9, 10 give the minimum and the maximum distances. Column 11 gives the reference for the distance. Columns 12–15 give the physical parameters of the sources (see text for details). Column 16 indicates whether the source is infrared quiet or not. Column 17 lists the assigned GMC complex or Spitzer Dark Clouds, and Col. 18 lists the EGO associations (see text for details.) (a) We revised the spectrophotometric distance of 1.81 kpc by Moisés et al. (2011). (b) Sources also included in the ATLASGAL top100 selection Giannetti et al. (2014).

**References.** (1) This work; (2) Sato et al. (2010); (3) Zhang et al. (2013); (4) Zhang et al. (2014); (5) Immer et al. (2013); (6) Sanna et al. (2014); (7) Moisés et al. (2011); (8) Wu et al. (2014); (9) Wu et al. (2012); (10) Reid et al. (2014); (11) Dutra et al. (2003); (12) Straw et al. (1987); (13) Xu et al. (2011); (14) Sanna et al. (2009); (15) Leurini et al. (2011); (16) Kurayama et al. (2011); (17) Russeil (2003); (18) Reid et al. (2009); (19) Brunthaler et al. (2009); (20) Urquhart et al. (2014b); (21) Krishnan et al. (2015); (22) Motogi et al. (2011).
<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Surface-density, Σ, versus temperature corrected mass estimates (using $T_d = 30$ K for the infrared-bright clumps). The coloured dashed lines in different shades show constant radius (green), $n_H$ number density (red), and escape speed (blue; see Tan et al. 2014). Colours and symbols are the same as in Fig. 8. We mark two massive cores with $M_{\text{MDC}} = 60 M_{\odot}$ (C1-S; Tan et al. 2013) and $55 M_{\odot}$ (CygX-N63; Bontemps et al. 2010). For comparison IRDC clumps (Kainulainen & Tan 2013) and cores are shown (Butler & Tan 2012). The dotted line marks the $1 \text{ g cm}^{-3}$ value of Krumholz & McKee (2008).</figcaption></figure>

embedded objects with $L_{\text{bol}} > 10^4 L_{\odot}$ corresponding to MYSOs or high-mass stars with UC-H II regions, while infrared-quiet clumps lack such bright embedded objects. From the initial sample of 210 sources above a 5 Jy beam-averaged peak flux density, we identify 80 infrared-quiet massive clumps, of which 50 are located within $d < 4.5$ kpc. About one-fourth of them have not been recognised by other surveys.

We estimate their physical properties based on the dust emission, and find that statistically both samples exhibit similar properties in terms of mass and size. Only their embedded (proto)stellar content is different, the infrared-quiet sample lacks luminous embedded objects. This suggests that infrared-quiet clumps are the precursors of the bright clumps, and thus they likely follow the same evolutionary path. A substantial fraction of infrared-quiet clumps exhibit signatures of ongoing star formation activity such as emission from shocked gas associated with material ejection. In particular, selective star formation tracers, such as EGO activity and Class II methanol maser emission, suggest that at least some infrared-quiet massive clumps host high-mass protostars. This makes them excellent candidates to be the potential precursors of the most massive stars forming in our Galaxy.

The overall low ratio of virial mass to mass of the sample suggests a collapse on the clump scale, which is consistent with a fast, dynamic star formation process. In this scenario massive clumps start to fragment and collapse before their final mass is accumulated indicating that strong self-gravity and global collapse is needed to build up rich clusters and the most massive stars.

High angular-resolution observations of infrared-quiet massive clumps are the next step to reveal their protostellar content and uncover the initial conditions for high-mass star formation.

**Acknowledgements.** We thank the referee for the careful reading of the manuscript. This work was partially funded by the ERC Advanced Investigator Grant GLOSTAR (247078). T.C.s. acknowledges support from the Deutsche Forschungsgemeinschaft, DFG via the SPP (priority programme) 1573 "Physics of the ISM". T.C.s. thanks R. Gutermuth for making available the PhotVis tool. This paper is based on data acquired with the Atacama Pathfinder

Experiment (APEX). APEX is a collaboration between the Max Planck Institute for Radioastronomy, the European Southern Observatory, and the Onsala Space Observatory. This research made use of data products from the Midcourse Space Experiment. Processing of the data was funded by the Ballistic Missile Defense Organization with additional support from NASA Office of Space Science. This research has also made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This research made use of Montage, funded by the National Aeronautics and Space Administration's Earth Science Technology Office, Computation Technologies Project, under Cooperative Agreement Number NCC5-626 between NASA and the California Institute of Technology. Montage is maintained by the NASA/IPAC Infrared Science Archive.

## References

Avison, A., Peretto, N., Fuller, G. A., et al. 2015, A&A, 577, A30

Bally, J., Anderson, L. D., Battersby, C., et al. 2010, A&A, 518, L90

Battersby, C., Bally, J., Ginsburg, A., et al. 2011, A&A, 535, A128

Beltrán, M. T., Sánchez-Monge, Á., Cesaroni, R., et al. 2014, A&A, 571, A52

Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953

Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339

Bernard, J.-P., Paradis, D., Marshall, D. J., et al. 2010, A&A, 518, L88

Bertoldi, F., & McKee, C. F. 1992, ApJ, 395, 140

Beuther, H., Zhang, Q., Greenhill, L. J., et al. 2004, ApJ, 616, L31

Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132

Bonnell, I. A., & Bate, M. R. 2006, MNRAS, 370, 488

Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 2001, MNRAS, 323, 785

Bontemps, S., Motte, F., Csengeri, T., & Schneider, N. 2010, A&A, 524, A18

Brand, J., & Blitz, L. 1993, A&A, 275, 67

Breen, S. L., Fuller, G. A., Caswell, J. L., et al. 2015, MNRAS, 450, 4109

Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 81

Butler, M. J., & Tan, J. C. 2012, ApJ, 754, 5

Carey, S. J., Noriega-Crespo, A., Mizuno, D. R., et al. 2009, PASP, 121, 76

Caswell, J. L., Fuller, G. A., Green, J. A., et al. 2010, MNRAS, 404, 1029

Caswell, J. L., Fuller, G. A., Green, J. A., et al. 2011, MNRAS, 417, 1964

Chen, X., Gan, C.-G., Ellingsen, S. P., et al. 2013, ApJS, 206, 9

Churchwell, E. 2002, ARA&A, 40, 27

Csengeri, T., Bontemps, S., Schneider, N., Motte, F., & Dib, S. 2011a, A&A, 527, A135

Csengeri, T., Bontemps, S., Schneider, N., et al. 2011b, ApJ, 740, L5

Csengeri, T., Urquhart, J. S., Schuller, F., et al. 2014, A&A, 565, A75

Csengeri, T., Leurini, S., Wyrowski, F., et al. 2016a, A&A, 586, A149

Csengeri, T., Weiss, A., Wyrowski, F., et al. 2016b, A&A, 585, A104

Csengeri, T., Bontemps, S., Wyrowski, F., et al. 2017, A&A, 600, L10

Cutri, R. M., Wright, E. L., Conrow, T., et al. 2012, Explanatory Supplement to the WISE All-Sky Data Release Products, Tech. rep.

Cyganowski, C. J., Whitney, B. A., Holden, E., et al. 2008, AJ, 136, 2391

Cyganowski, C. J., Brogan, C. L., Hunter, T. R., et al. 2014, ApJ, 796, L2

Da Rio, N., Robberto, M., Hillenbrand, L. A., Henning, T., & Stassun, K. G. 2012, ApJ, 748, 14

Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium (Princeton University Press)

Duarte-Cabral, A., Bontemps, S., Motte, F., et al. 2013, A&A, 558, A125

Dunham, M. K., Rosolowsky, E., Evans, II, N. J., Cyganowski, C., & Urquhart, J. S. 2011, ApJ, 741, 110

Dutra, C. M., Bica, E., Soares, J., & Barbuy, B. 2003, A&A, 400, 533

Elia, D., Molinari, S., Fukui, Y., et al. 2013, ApJ, 772, 45

Engelbracht, C. W., Blaylock, M., Su, K. Y. L., et al. 2007, PASP, 119, 994

Galván-Madrid, R., Liu, H. B., Zhang, Z.-Y., et al. 2013, ApJ, 779, 121

Garay, G., Faúndez, S., Mardones, D., et al. 2004, ApJ, 610, 313

Giannetti, A., Wyrowski, F., Brand, J., et al. 2014, A&A, 570, A65

Ginsburg, A., Bally, J., Battersby, C., et al. 2015, A&A, 573, A106

Green, J. A., Caswell, J. L., Fuller, G. A., et al. 2010, MNRAS, 409, 913

Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13

Gutermuth, R. A., Myers, P. C., Megeath, S. T., et al. 2008, ApJ, 674, 336

Guzmán, A. E., Garay, G., & Brooks, K. J. 2010, ApJ, 725, 734

Habibi, M., Stolte, A., Brandner, W., Hußmann, B., & Motohara, K. 2013, A&A, 556, A26

Hoare, M. G., & Franco, J. 2007, Astrophys. Space Sci. Proc., 1, 61

Hoare, M. G., Kurtz, S. E., Lizano, S., Keto, E., & Hofner, P. 2007, Protostars and Planets V, 181
Hoare, M. G., Purcell, C. R., Churchwell, E. B., et al. 2012, PASP, 124, 939

Hunter, T. R., Churchwell, E., Watson, C., et al. 2000, AJ, 119, 2711

Immer, K., Reid, M. J., Menten, K. M., Brunthaler, A., & Dame, T. M. 2013, A&A, 553, A117

Jackson, J. M., Rathborne, J. M., Shah, R. Y., et al. 2006, ApJS, 163, 145

Jackson, J. M., Finn, S. C., Rathborne, J. M., Chambers, E. T., & Simon, R. 2008, ApJ, 680, 349

Kainulainen, J., & Tan, J. C. 2013, A&A, 549, A53

Kauffmann, J., & Pillai, T. 2010, ApJ, 723, L7

Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, II, N. J., & Lee, C. W. 2008, A&A, 487, 993

Kauffmann, J., Pillai, T., Shetty, R., Myers, P. C., & Goodman, A. A. 2010, ApJ, 712, 1137

Kauffmann, J., Pillai, T., & Goldsmith, P. F. 2013, ApJ, 779, 185

Keto, E. 2002, ApJ, 580, 980

Keto, E. 2003, ApJ, 599, 1196

Klaassen, P. D., Testi, L., & Beuther, H. 2012, A&A, 538, A140

Klein, T., Ciechanowicz, M., Leinz, C., et al. 2014, IEEE Transactions on Terahertz Science and Technology, 4, 588

König, C., Urquhart, J. S., Csengeri, T., et al. 2017, A&A, 599, A139

Krishnan, V., Ellingsen, S. P., Reid, M. J., et al. 2015, ApJ, 805, 129

Krumholz, M. R., & McKee, C. F. 2008, Nature, 451, 1082

Kuiper, R., & Yorke, H. W. 2013, ApJ, 772, 61

Kurayama, T., Nakagawa, A., Sawada-Satoh, S., et al. 2011, PASJ, 63, 513

Leurini, S., Codella, C., Zapata, L., et al. 2011, A&A, 530, A12

Leurini, S., Codella, C., López-Sepulcre, A., et al. 2014, A&A, 570, A49

López, C., Bronfman, L., Nyman, L.-Å., May, J., & Garay, G. 2011, A&A, 534, A131

Louvet, F., Motte, F., Hennebelle, P., et al. 2014, A&A, 570, A15

Lumsden, S. L., Hoare, M. G., Urquhart, J. S., et al. 2013, ApJS, 208, 11

McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565

McKee, C. F., & Tan, J. C. 2003, ApJ, 585, 850

Megeath, S. T., & Tieftrunk, A. R. 1999, ApJ, 526, L113

Megeath, S. T., Guttermuth, R., Muzerolle, J., et al. 2012, AJ, 144, 192

Menten, K. M. 1991, ApJ, 380, L75

Moisés, A. P., Damineli, A., Figuerêdo, E., et al. 2011, MNRAS, 411, 705

Molinari, S., Brand, J., Cesaroni, R., & Palla, F. 1996, A&A, 308, 573

Molinari, S., Brand, J., Cesaroni, R., & Palla, F. 2000, A&A, 355, 617

Molinari, S., Pezzuto, S., Cesaroni, R., et al. 2008, A&A, 481, 345

Mookerjea, B., Kramer, C., Nielbock, M., & Nyman, L.-Å. 2004, A&A, 426, 119

Motogi, K., Sorai, K., Habe, A., et al. 2011, PASJ, 63, 31

Motte, F., Schilke, P., & Lis, D. C. 2003, ApJ, 582, 277

Motte, F., Bontemps, S., Schilke, P., et al. 2007, A&A, 476, 1243

Mottram, J. C., Hoare, M. G., Davies, B., et al. 2011, ApJ, 730, L33

Mufson, S. L., & Liszt, H. S. 1979, ApJ, 232, 451

Nguyen Luong, Q., Motte, F., Schuller, F., et al. 2011, A&A, 529, A41

Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943

Peng, T.-C., Wyrowski, F., van der Tak, F. F. S., Menten, K. M., & Walmsley, C. M. 2010, A&A, 520, A84

Peretto, N., & Fuller, G. A. 2009, A&A, 505, 405

Peretto, N., Fuller, G. A., Duarte-Cabral, A., et al. 2013, A&A, 555, A112

Peretto, N., Lenfestey, C., Fuller, G. A., et al. 2016, A&A, 590, A72

Purcell, C. R., Hoare, M. G., Cotton, W. D., et al. 2013, ApJS, 205, 1

Rathborne, J. M., Jackson, J. M., Chambers, E. T., et al. 2005, ApJ, 630, L181

Rathborne, J. M., Jackson, J. M., & Simon, R. 2006, ApJ, 641, 389

Rathborne, J. M., Jackson, J. M., Zhang, Q., & Simon, R. 2008, ApJ, 689, 1141

Reid, M. J., Menten, K. M., Zheng, X. W., et al. 2009, ApJ, 700, 137

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130

Russeil, D. 2003, A&A, 397, 133

Russeil, D., Zavagno, A., Motte, F., et al. 2010, A&A, 515, A55

Rygl, K. L. J., Brunthaler, A., Sanna, A., et al. 2012, A&A, 539, A79

Sánchez-Monge, Á., Cesaroni, R., Beltrán, M. T., et al. 2013, A&A, 552, L10

Sanna, A., Reid, M. J., Moscadelli, L., et al. 2009, ApJ, 706, 464

Sanna, A., Reid, M. J., Menten, K. M., et al. 2014, ApJ, 781, 108

Sato, M., Reid, M. J., Brunthaler, A., & Menten, K. M. 2010, ApJ, 720, 1055

Schneider, N., Csengeri, T., Bontemps, S., et al. 2010, A&A, 520, A49

Schuller, F., Menten, K. M., Contreras, Y., et al. 2009, A&A, 504, 415

Simon, R., Rathborne, J. M., Shah, R. Y., Jackson, J. M., & Chambers, E. T. 2006, ApJ, 653, 1325

Siringo, G., Kreysa, E., Kovács, A., et al. 2009, A&A, 497, 945

Sobolev, A. M., Cragg, D. M., & Godfrey, P. D. 1997, MNRAS, 288, L39

Sridharan, T. K., Beuther, H., Schilke, P., Menten, K. M., & Wyrowski, F. 2002, ApJ, 566, 931

Sridharan, T. K., Rao, R., Qiu, K., et al. 2014, ApJ, 783, L31

Straw, S., Hyland, A. R., Jones, T. J., et al. 1987, ApJ, 314, 283

Svoboda, B. E., Shirley, Y. L., Battersby, C., et al. 2016, ApJ, 822, 59

Tackenberg, J., Beuther, H., Henning, T., et al. 2014, A&A, 565, A101

Tackenberg, J., Beuther, H., Henning, T., et al. 2012, A&A, 540, A113

Tan, J. C., Kong, S., Butler, M. J., Caselli, P., & Fontani, F. 2013, ApJ, 779, 96

Tan, J. C., Beltrán, M. T., Caselli, P., et al. 2014, Protostars and Planets VI, 149

Tan, J. C., Kong, S., Zhang, Y., et al. 2016, ApJ, 821, L3

Thompson, M. A., Gibb, A. G., Hatchell, J. H., Wyrowski, F., & Pestalozzi, M. R. 2005, in Protostars and Planets V Posters, 8109

Traficante, A., Fuller, G. A., Peretto, N., Pineda, J. E., & Molinari, S. 2015, MNRAS, 451, 3089

Urquhart, J. S., Busfield, A. L., Hoare, M. G., et al. 2007, A&A, 461, 11

Urquhart, J. S., Moore, T. J. T., Schuller, F., et al. 2013a, MNRAS, 431, 1752

Urquhart, J. S., Thompson, M. A., Moore, T. J. T., et al. 2013b, MNRAS, 435, 400

Urquhart, J. S., Moore, T. J. T., Csengeri, T., et al. 2014, MNRAS, 443, 1555

Ward-Thompson, D., André, P., Crutcher, R., et al. 2007, Protostars and Planets V, 33

Wienen, M., Wyrowski, F., Schuller, F., et al. 2012, A&A, 544, A146

Wienen, M., Wyrowski, F., Menten, K. M., et al. 2015, A&A, 579, A91

Wienen, M., Wyrowski, F., & Menten, K. M., et al. 2017, A&A, submitted

Williams, J. P., Blitz, L., & McKee, C. F. 2000, Protostars and Planets IV, 97

Wood, D. O. S., & Churchwell, E. 1989, ApJ, 340, 265

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Wu, Y. W., Xu, Y., Menten, K. M., Zheng, X. W., & Reid, M. J. 2012, IAU Symp. 287, eds. R. S. Booth, W. H. T. Vlemmings, & E. M. L. Humphreys, 425

Wu, Y. W., Sato, M., Reid, M. J., et al. 2014, A&A, 566, A17

Xu, Y., Moscadelli, L., Reid, M. J., et al. 2011, ApJ, 733, 25

Zapata, L. A., Palau, A., Ho, P. T. P., et al. 2008, A&A, 479, L25

Zavagno, A., Russeil, D., Motte, F., et al. 2010, A&A, 518, L81

Zhang, B., Reid, M. J., Menten, K. M., et al. 2013, ApJ, 775, 79

Zhang, B., Moscadelli, L., Sato, M., et al. 2014, ApJ, 781, 89

Zhang, Q., Wang, Y., Pillai, T., & Rathborne, J. 2009, ApJ, 696, 268

Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481
# Appendix A: Mid-infrared photometry

## A.1. MIPSGAL upper limits

Using the publicly available calibrated MIPSGAL tiles, we searched for embedded compact objects spatially coinciding with the ATLASGAL dust peaks. From the 219 brightest submm sources we find that 27 do not show any clear point source within 5" from the dust peaks. After a visual inspection, we extracted an upper limit for these sources by taking the median of the pixel values within a radius of 3", i.e. half of the FWHM of the PSF at 24 µm. We then multiply this value by the beam area assuming a Gaussian beam with a FWHM of 6" to obtain an upper limit on the source flux. We estimate the error by calculating the REDMSQ<fnref n="5" /> value in the same area, in the footprint of the PSF, as was defined in Megeath et al. (2012). We extract upper limits between 0.041–1.368 Jy/beam with an average of 0.4 Jy/beam.

<figure><img src="image_15.png" /><figcaption>Fig. A.1. Cross-calibration between flux densities from the WISE point source catalogue and the MIPSGAL survey.</figcaption></figure>

## A.2. MIPSGAL aperture photometry at 24 µm

We performed aperture photometry towards the dust peaks that are clearly associated with a 24 µm point source and without saturated pixels. In greater detail, we extracted the flux density of point sources from the MIPSGAL data using the PhotVis v1.10 tool (Gutermuth et al. 2008) positionally coinciding with the ATLASGAL source. We adopted a zero point magnitude of 12.448, which was calculated assuming a pixel size of 1"/25, an aperture correction of 2.05, which is used for apertures of 7" with a sky annulus between 20" to 32" (Engelbracht et al. 2007)<fnref n="6" />. We chose these radii because the background emission varies substantially and a larger sky annulus helps to avoid contamination by nearby bright sources. However, our background values still show variations and so may still be contaminated by extended emission. Therefore, we caution that our values are subject to uncertainties.

We then converted the measured magnitudes to flux density using $F_v = c_{corr} \times F_{v0} \times 10^{(-\text{mag}/2.5)}$, where $F_{v0}$ corresponds to the zero magnitude flux level, which is 7.17 Jy and mag is the magnitude given by PhtoVis. The targeted sample consists of deeply embedded objects, where $\beta \sim 1–2$ corresponds to rising SEDs towards longer wavelengths. Therefore, we applied a colour correction ($c_{corr}$) of 0.981 for objects with $\beta \sim 1$ (Engelbracht et al. 2007).

## A.3. Comparing MIPSGAL and WISE flux densities

We measure 24 µm flux densities (and upper limits) between 0.041 and 6.88 Jy for these sources, with a mean of 1.84 and a median of 0.81 Jy. For 75 sources in the sample we also have

photometry from the WISE point source catalogue (Wright et al. 2010). Since there is an overlap in the wavelength ranges of the WISE band at 22 µm and the MIPS detector at 24 µm, we attempted to cross-calibrate the photometry of the two measurements. We found that numerous sources appearing in the WISE catalogue and not saturated with MIPS are flagged as likely spurious detections. This is likely due to source confusion and extended nebulosity surrounding the majority of our weak sources, which results in bad quality flags for these sources in the WISE point source catalogue. Therefore, we only have 21 sources with reliable flux measurement from WISE, of which only 19 have smaller than 5" angular offset from the ATLASGAL peak position.

Our attempt to cross-calibrate the two instruments is shown in Fig. A.1, where we also correct for a spectral index of 1.5 that has been found between MSX and MIPSGAL according to the MIPSGAL data delivery document<fnref n="7" />. With some scatter, we find similar flux density measurements below ~4 Jy with the two instruments. Since the MIPS detector starts to saturate at 2 Jy, the flattening of the MIPSGAL fluxes between ~2–4 Jy is likely due to the saturation or non-linearity of the detector. In addition, since WISE has a beam that is two times larger (the MIPS PSF is 6" at 24 µm, while WISE has a 12" angular resolution), it may include emission from the surrounding diffuse emission, not only the protostar. Since saturation and non-linearity of the detector may significantly influence the determined values, in our analysis for sources above 2 Jy we rely on the WISE point source catalogue values.

<footnote n="5">⁵ REDMSQ = $\sqrt{\text{median}((S[i, j] - \text{median}(S[i, j]))^2)}$.</footnote>

<footnote n="6">⁶ We used the following formula to calculate the zero point magnitude, $\text{mag}_0 = 2.5 \times \log_{10}(\frac{F_{v0}}{10^6 \times \Omega \times \text{app}_{corr}})$, where $\Omega$ corresponds to the pixel size in steradian and $F_{v0}$ is 7.17 Jy.</footnote>

<footnote n="7">⁷ http://irsa.ipac.caltech.edu/data/SPITZER/MIPSGAL/images/mipsgalDeliveryGuide_v3_29aug08.pdf</footnote>
Appendix B: Infrared-quiet massive clumps

<figure><img src="image_16.png" /><figcaption>**Fig. B.1.** *Left:* mid-infrared emission from MIPSGAL at 24 µm or WISE at 22 µm in colour scale with the 870 µm contours from ATLASGAL. The arrows mark the position of infrared-quiet clumps. *Right:* 870 µm emission with contours from ATLASGAL. The blue crosses mark the positions of UC-H II regions from CORNISH. The figure label shows the corresponding ATLASGAL source name.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
Appendix C: SiO (8–7) lines towards infrared-quiet massive clumps with APEX

<figure><img src="image_23.png" /><figcaption>Fig. C.1. SiO (8–7) detections with the APEX/FLASH$^+$ receiver of a selection of infrared-quiet massive clumps. The red dashed line shows the fitted Gaussians, and the grey histogram shows the line area. The blue dashed line indicates the $v_{lsr}$ of the sources.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>Fig. C.1. continued.</figcaption></figure>