# Ammonia from cold high-mass clumps discovered in the inner Galactic disk by the ATLASGAL survey

M. Wienen, F. Wyrowski, F. Schuller, K. M. Menten, C. M. Walmsley, L.

Bronfman, F. Motte

► To cite this version:

M. Wienen, F. Wyrowski, F. Schuller, K. M. Menten, C. M. Walmsley, et al.. Ammonia from cold high-mass clumps discovered in the inner Galactic disk by the ATLASGAL survey. Astronomy and Astrophysics - A&A, 2012, 544, 10.1051/0004-6361/201118107. insu-03622102

HAL Id: insu-03622102

https://insu.hal.science/insu-03622102

Submitted on 29 Mar 2022

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
# Ammonia from cold high-mass clumps discovered in the inner Galactic disk by the ATLASGAL survey*

M. Wienen<fnref n="1" />,**, F. Wyrowski<fnref n="1" />, F. Schuller<fnref n="1" />,<fnref n="2" />, K. M. Menten<fnref n="1" />, C. M. Walmsley<fnref n="3" />,<fnref n="4" />, L. Bronfman<fnref n="5" />, and F. Motte<fnref n="6" />

Received 16 September 2011 / Accepted 26 May 2012

## ABSTRACT

**Context.** The APEX Telescope Large Area Survey: the GALaxy (ATLASGAL) is an unbiased continuum survey of the inner Galactic disk at 870 µm. It covers ±60° in Galactic longitude and aims to find all massive clumps at various stages of high-mass star formation in the inner Galaxy, particularly the earliest evolutionary phases.

**Aims.** We aim to determine properties such as the gas kinetic temperature and dynamics of new massive cold clumps found by ATLASGAL. Most importantly, we derived their kinematical distances from the measured line velocities.

**Methods.** We observed the ammonia (J, K) = (1, 1) to (3, 3) inversion transitions toward 862 clumps of a flux-limited sample of submm clumps detected by ATLASGAL and extracted ¹³CO (1–0) spectra from the Galactic Ring Survey (GRS). We determined distances for a subsample located at the tangential points (71 sources) and for 277 clumps whose near/far distance ambiguity is resolved.

**Results.** Most ATLASGAL clumps are cold with rotational temperatures from 10–30 K with a median of 17 K. They have a wide range of NH₃ linewidths (1–7 km s⁻<fnref n="1" />) with 1.9 km s⁻<fnref n="1" /> as median, which by far exceeds the thermal linewidth, as well as a broad distribution of high column densities from 10¹<fnref n="4" /> to 10¹<fnref n="6" /> cm⁻<fnref n="2" /> (median of 2 × 10¹<fnref n="5" /> cm⁻<fnref n="2" />) with an NH₃ abundance in the range of 5 to 30 × 10⁻⁸. ATLASGAL sources are massive, ≥100 M☉, and a fraction of clumps with a broad linewidth is in virial equilibrium. We found an enhancement of clumps at Galactocentric radii of 4.5 and 6 kpc. The comparison of the NH₃ lines as high-density probes with the GRS ¹³CO emission as low-density envelope tracer yields broader linewidths for ¹³CO than for NH₃. The small differences in derived clump velocities between NH₃ (representing dense core material) and ¹³CO (representing more diffuse molecular cloud gas) suggests that the cores are essentially at rest relative to the surrounding giant molecular cloud.

**Conclusions.** The high detection rate (87%) confirms ammonia as an excellent probe of the molecular content of the massive, cold clumps revealed by ATLASGAL. A clear trend of increasing rotational temperatures and linewidths with evolutionary stage is seen for source samples ranging from 24 µm dark clumps to clumps with embedded HII regions. The survey provides the largest ammonia sample of high-mass star forming clumps and thus presents an important repository for the characterization of statistical properties of the clumps and the selection of subsamples for detailed, high-resolution follow-up studies.

**Key words.** surveys – submillimeter: general – radio lines: ISM – ISM: molecules – ISM: kinematics and dynamics – stars: formation

## 1. Introduction

High-mass stars are essential for the evolution of galaxies because they energize the interstellar medium and release heavy elements, which are determining Galactic cooling mechanisms. Although the understanding of high-mass star formation is consequently of great importance, only little is known about the early stages of the formation of massive stars in contrast to the well-established evolutionary sequence of isolated low-mass stars (André et al. 2000). This is partly because of the difficult observational conditions that one encounters. High-mass stars form in clusters and are rarer according to the initial mass function (Kroupa et al. 2011), therefore they are usually found at

larger distances than molecular clouds, in which low-mass stars form. Moreover, high-mass stars evolve on shorter time scales and within dynamic and complex environments and are still deeply embedded in their earliest phases.

It is well known that high-mass star formation starts in giant molecular clouds (Dame et al. 2001). In analogy to low-mass star formation, one expects a cold and massive starless clump of ≥500 M☉ consisting of molecular gas for the initial conditions of massive star/cluster formation (Beuther et al. 2007). Young massive (proto) stars evolving in dense cores within the clumps emit ultraviolet radiation that leads to the formation of ultracompact HII regions (UCHIIRs). Many of them were found in radio continuum surveys with the VLA (e.g. Wood & Churchwell 1989b), which revealed their distribution and number. Those regions were very helpful in tracing a more evolved stage of massive star formation in the Galaxy. Often UCHIIRs are associated with so-called hot molecular cores (Cesaroni et al. 1992), which are possible precursors of UCHIIRs, also called high-mass protostellar objects (HMPOs) or massive young stellar

<footnote n="1">¹ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
e-mail: mwienen@mpifr-bonn.mpg.de</footnote>

<footnote n="2">² Alonso de Cordova 3107, Casilla 19001, Santiago 19, Chile</footnote>

<footnote n="3">³ Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy</footnote>

<footnote n="4">⁴ Dublin Institute of Advanced Studies, Fitzwilliam Place 31, Dublin 2, Ireland</footnote>

<footnote n="5">⁵ Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile</footnote>

<footnote n="6">⁶ Laboratoire AIM, CEA/IRFU - CNRS/INSU - Université Paris Diderot, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France</footnote>

<footnote>* Full Tables 1–5 and reduced spectra are 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/544/A146</footnote>

<footnote>** Member of the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne.</footnote>
objects (MYSOs), which have also been discovered recently (Beuther et al. 2002). Still earlier phases are deeply embedded in a massive envelope and are so cold that mostly they cannot be probed by means of mid-infrared emission (e.g. Sridharan et al. 2005; Peretto & Fuller 2010; Motte et al. 2007). Under certain circumstances, these clumps can be seen in absorption against a bright background provided mostly by emission from polycyclic aromatic hydrocarbon molecules (PAHs) as infrared dark clouds (IRDCs, Perault et al. 1996; Egan et al. 1998). Simon et al. (2006b) and Peretto & Fuller (2010) have presented large cata- logues of these IRDCs. But to date, only a small fraction of them has been investigated in detail (Carey et al. 1998, 2000; Pillai et al. 2006; Rathborne et al. 2006). Previous studies revealed that IRDCs are cold (<25 K), dense (>10$^5$ cm$^{-3}$) and have high col- umn densities (~10$^{23}$–10$^{25}$ cm$^{-2}$, Rathborne et al. 2006; Simon et al. 2006a). Their sizes (~2 pc) and masses (~10$^2$–10$^4$ M$_\odot$) are similar to those of cluster-forming molecular clumps. This, together with the cold temperatures and fragmented substruc- ture of IRDCs, suggests that they are protoclusters (Rathborne et al. 2006, 2005). Kinematic distances to IRDC samples in the first and fourth Galactic quadrants have been determined using $^{13}$CO (1–0) and CS (2–1) emission (Jackson et al. 2006, 2008). However, surveys targeted at the early phases of high-mass star formation (e.g. Beuther et al. 2002; Wood & Churchwell 1989a) usually only trace one of those stages, e.g. IRDCs, HMPOs or UCHIRs. In contrast, a first attempt to obtain an unbiased sam- ple using a complete dust continuum imaging at the scale of the molecular complex Cygnus X (Motte et al. 2007) was able to derive some more systematic results on the existence of a cold phase (IR-quiet massive dense cores and Class 0-like massive protostars) for high-mass star formation (Bontemps et al. 2010; Csengeri et al. 2011). Another survey of the molecular complex NGC 6334/NGC 6357 was conducted by Russeil et al. (2010), who identified IR-quiet massive dense cores and estimated high- mass protostellar lifetimes. The IR-quiet massive dense cores are suggested to host massive class 0 protostars as seen by the Herschel Space Telescope (Nguyen Luong et al. 2011a).

Only a few hundred high-mass proto- or young stellar ob- jects have been studied up to now. Hence, these and especially objects in still earlier phases of massive star formation need to be investigated in more detail. To achieve significant progress in that, a large-area survey of the Galactic plane was conducted, which now provides a global view of star formation at submil- limeter wavelengths. The APEX Telescope Large Area Survey of the GALaxy at 870 µm (ATLASGAL) is the first unbiased continuum survey of the whole inner Galactic disk at 870 µm (Schuller et al. 2009). It aims to find all massive clumps that form high-mass stars in the inner Galaxy by using the Large APEX BOlometer CAmera (LABOCA), which is an array with 295 bolometer pixels operated at the APEX telescope (Siringo et al. 2007, 2008) with a field of view of 11' and a beamwidth of 19.2" FWHM at the wavelength of 870 µm.

ATLASGAL, which reaches a Galactic longitude of ±60° and latitude of ±1.5°, is able to detect objects associated with massive star formation at various stages and to compare them (Schuller et al. 2009). Other Galactic plane surveys, conducted over a similar range, are the MSX survey covering 8 to 21 µm (Price et al. 2001), 2MASS around 2 µm (Skrutskie et al. 2006), GLIMPSE from 3 to 8 µm (Benjamin et al. 2003), MIPSGAL at 24 and 70 µm (Carey et al. 2009) and HiGAL from 70 to 500 µm (Molinari et al. 2010). While submillimeter dust con- tinuum surveys are essential for identifying high-mass clumps, they lack information on important parameters, especially the distances to the newly found sources. These are needed to

determine other properties such as their masses and luminosities. But all of this can be addressed by molecular line observations of the high-mass star forming regions. Because their density is at about the critical NH$_3$ density, this molecule is appropriate for determining properties of the clumps without contamination from large-scale cloud structures with lower density. Since the molecular gas of the cores is very dense (~10$^5$ cm$^{-3}$, Beuther et al. 2002) and cold, many molecules such as CS and CO are partly frozen onto dust grains. In contrast to those, the frac- tional gas abundance of ammonia remains constant in prestel- lar cores (Tafalla et al. 2002). Hence, we carried out follow- up observations of northern ATLASGAL sources in the lowest NH$_3$ inversion transitions.

NH$_3$ is known as a reliable temperature probe of interstellar clouds (Tafalla et al. 2004; Walmsley & Ungerechts 1983). The rotational energy levels are given by the total angular momen- tum, *J*, and its projection along the molecular axis, *K*. Radiative transitions between different *K*-ladders are forbidden and the lowest metastable energy levels, for which *J* = *K*, are thus colli- sionally excited. The intensity ratio of their inversion transitions therefore provides the rotational temperature of the gas (Ho & Townes 1983), which can be used to estimate the kinetic tem- perature of the cores. This is needed for proper mass estimates from the submm data. Independently, the linewidth can be used to derive virial masses. Moreover, the inversion transitions are split into distinct hyperfine components and their ratio provides a measure of the optical depth, knowledge of which leads to reli- able column density and rotational temperature determinations.

As this article was ready for submission, we became aware of the recent study by Dunham et al. (2011), who observed NH$_3$ lines of sources in the first Galactic quadrant as well, iden- tified by the Bolocam Galactic Plane Survey (BGPS), which measures the 1.1 mm continuum emission. In contrast to our sample, the observations of Dunham et al. (2011) do not cover the whole northern Galactic longitude range up to 60°, but are conducted only within four different ranges in longitude, and within a smaller Galactic latitude range of ±0.5° compared to our NH$_3$ survey. Our analysis is thus complementary to the re- sults of Dunham et al. (2011).

In Sect. 2, we describe which ATLASGAL sources were se- lected for our NH$_3$ follow-up observations and provide details of the measurements. In Sect. 3, we present how we reduced the data and derived different clump properties from fits to the spectra of NH$_3$ inversion transitions. In Sect. 4, different corre- lations of the NH$_3$ line parameters such as the velocity distri- bution, linewidth, rotational temperature and column density are shown. Then, we compare these with submillimeter dust con- tinuum properties of the clumps such as the H$_2$ column den- sity in Sect. 5 and estimate gas and virial masses. We analyse additional line parameters derived from $^{13}$CO (1–0) emission (Jackson et al. 2006) in Sect. 6 to investigate clump-to-cloud motions etc. In Sect. 7, the determined gas properties are com- pared with those of other high-mass dust-selected star forming regions. A summary and conclusions of this NH$_3$ investigation are given in Sect. 8.

This paper analyses statistics of the NH$_3$ data obtained for the northern ATLASGAL sources. We provide their near and far distance, but do not distinguish between both, although some of the molecular clouds have known distances from a previous study (Roman-Duval et al. 2009). In a second paper (Wienen et al., in prep.) ammonia lines of ATLASGAL sources in the southern hemisphere will be studied and distances to southern molecular clouds will be derived.
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Relative number distribution of NH₃ observations as solid black curve and of the whole northern ATLASGAL sample as dashed red curve vs. the logarithm of the 870 µm peak flux density.</figcaption></figure>

## 2. Observations

### 2.1. Source selection

From a preliminary ATLASGAL point source catalogue, we selected a flux-limited subsample of compact (smaller than 50′′) clumps down to about 0.4 Jy/beam peak flux density, whose positions were selected using the Miriad task sfind (Sault et al. 1995), and observed ~63% of these 1361 sources. From 2008 to 2010 we observed the ammonia emission of a total of 862 dust clumps in a Galactic longitude range from 5° to 60° and latitude within ±1.5°. The 870 µm peak flux density distribution of observed sources is plotted as a solid black curve in Fig. 1 and that of the whole ATLASGAL sample of 1361 sources as a dashed red curve. Their comparison reveals that the 870 µm peak flux distributions of both samples are similar and therefore we expect the NH₃ subsample to be representative of the ATLASGAL sample as a whole. Most NH₃ observations as well as most ATLASGAL sources from the whole sample have a peak flux density of ~0.8 Jy/beam, the distributions above those values can be described by a power law.

We estimate the ATLASGAL sample to be complete above 1 Jy/beam peak flux density, which results in 644 observed clumps. Of those the NH₃ sample contains 406 sources (63%) within this flux limit. Seventy-three clumps (11%) were observed in previous ammonia surveys (Sridharan et al. 2002; Churchwell et al. 1990; Molinari et al. 1996; Pillai et al. 2006).

### 2.2. Observational setup

Observations of the NH₃ (1, 1), (2, 2) and (3, 3) inversion transitions were made from 2008 to 2010 with the Effelsberg 100-m telescope. The frontend was a 1.3 cm cooled HEMT receiver, which is used in a frequency range between 18 GHz and 26 GHz and mounted at the primary focus. In 2008 January, two different spectrometers were used, the 8192 channel autocorrelator (AK 90) and the Fast Fourier Transform Spectrometer (FFTS). In May 2008, the AK90 was replaced by the FFTS, which possesses a greater bandwidth. The autocorrelator contains eight individual modules. A bandwidth of 20 MHz was chosen for each of them, which results in a spectral resolution of about 0.5 km s⁻<fnref n="1" />, while the FFTS consists of two modules with a chosen bandwidth of 500 MHz each and a spectral resolution of ~0.7 km s⁻<fnref n="1" />. Both spectrometers are able to measure two polarizations of the three observed NH₃ lines simultaneously.

The scans observed in 2008 January with the two spectrometers were summed to obtain a better signal-to-noise ratio (S/N). The beamwidth (FWHM) at the frequencies of the (1, 1) to (3, 3) ammonia inversion transitions at about 24 GHz is 40′′. Pointed observations towards the dust emission peaks were conducted in frequency-switching mode with a frequency throw of 7.5 MHz. The median system temperature was about 70 K, with a zenith optical depth that usually varied between 0.02 and 0.09 in the beginning of the year in about January and was more constant, about 0.07 some months later, in about May and October. The total integration time for each source was about 5 min. For some clumps at low declination, which could only be observed at low elevations and consequently have high system temperatures of ~150 K, longer integration times of up to 25 min were chosen.

## 3. Data reduction and analysis

To reduce the NH₃ spectra we used the CLASS software<fnref n="1" />. Because the ammonia lines from each source are detected in two polarizations, all scans that belong to one polarization of an inversion transition were summed and then the spectra of both polarizations were averaged. The frequency switching results in two lines, which were inverted and separated in frequency by 7.5 MHz. This procedure was reversed by a folding, which shifts the two lines by 3.75 MHz and subtracts one of them from the other. The large frequency throw led to fluctuations of the baseline, which needed to be corrected. For this we subtracted a polynomial baseline of the order of 3 to 7, which produces inaccuracies in the subsequent analysis of the spectra. The NH₃(1, 1) hyperfine structure of very many ammonia sources shows different line profiles with various linewidths, which had to be reduced in the same way. Hence, we developed a strategy to set windows, spaces around NH₃ lines, that were excluded from the baseline fitting: the width of the window around each hyperfine component of the (1, 1) line was chosen according to the width of the main line; the order of the baseline depends on the linewidth as well. A polynomial of high order, up to 7, was subtracted for lines with narrow width, <2 km s⁻<fnref n="1" />, while a low order, mostly 3, was chosen for lines with a broad width, greater than 4 km s⁻<fnref n="1" />. To test if this introduced any systematic errors, we compared the linewidths of simulated NH₃ (1, 1) lines with the linewidths resulting from a fit of the same spectra added on observed baselines. For linewidths more narrow than 1.5 km s⁻<fnref n="1" /> the model and the fitted observation agreed, the deviation of the two becomes slightly higher with increasing linewidths. However, the contribution of systematic errors for lines with broad linewidths is only small, ~6% on average. Some examples for reduced and calibrated spectra of observed (1, 1) to (3, 3) inversion transitions are given in Fig. 2.

The hyperfine structure of the NH₃(1, 1) line was fitted with CLASS. Assuming that all hyperfine components have the same excitation temperature, a nonlinear least-squares fit was made to the observed spectra with the CLASS software, which yielded as independent parameters radial velocity, $v_{LSR}$, the linewidth, $\Delta v$, at the full width at half maximum of a Gaussian profile and the optical depth of the main line, $\tau_m$, with their errors.

Because the hyperfine satellite lines of the (2, 2) and (3, 3) transitions are mostly too weak to be detected, their optical depth could not be determined. A single Gaussian was fitted to the main line of the (2, 2) and (3, 3) transitions with CLASS. The three independent fit parameters were the integrated area, A, the radial velocity of the main line, $v_{LSR}$, and the FWHM linewidth,

<footnote n="1">¹ Available at http://www.iram.fr/IRAMFR/GILDAS</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 2. Reduced and calibrated spectra of observed NH₃(1, 1), (2, 2) and (3, 3) inversion transitions, the fit is shown in green for some sources. Towards G10.36–0.15 and G23.27–0.26 two lines are seen at different LSR velocities.</figcaption></figure>

Δv. The temperature of the NH₃ (1, 1) line was derived from the peak intensity of the Gaussian fit to the main line. For the (2, 2) and (3, 3) transitions the temperatures were calculated with

$$T_{\text{MB}} = \frac{A}{\Delta v \times 1.064} \qquad (1)$$

with the error given by the rms value. If lines were not detected, the noise level was used as an upper limit. We list the positions of the sources, the optical depths of the NH₃ (1, 1) lines, τ(1, 1), their LSR velocities, v(1, 1), linewidths, Δv(1, 1), and main beam

brightness temperatures, T<sub>MB</sub>(1, 1), together with their formal fit errors in Table 1. The LSR velocities, linewidths, and main beam brightness temperatures of the (2, 2) and (3, 3) lines together with their formal fit errors are given in Table 2. The other physical parameters such as the rotational temperature (T<sub>rot</sub>), the kinetic temperature (T<sub>kin</sub>) and the ammonia column density (N<sub>NH₃</sub>) were derived using the standard formulation for NH₃ spectra (Ho & Townes 1983; Ungerechts et al. 1986), see Sect. 4.3. They are given with the errors (1σ) calculated from Gaussian error propagation in Table 3.
<table><caption>**Table 1.** NH₃(1, 1) line parameters.</caption><thead><tr><th>Name</th><th>RA<sup>1</sup><br>(J2000)</th><th>Dec<sup>1</sup><br>(J2000)</th><th>τ(1, 1)</th><th>ν(1, 1)<br>(km s<sup>-1</sup>)</th><th>Δν(1, 1)<br>(km s<sup>-1</sup>)</th><th>T<sub>MB</sub>(1, 1)<br>(K)</th></tr></thead><tbody><tr><td>G5.39+0.19</td><td>17 57 12.9150</td><td>-24 12 28.150</td><td>1.71 (±0.29)</td><td>10.58 (±0.03)</td><td>1.42 (±0.06)</td><td>1.7 (±0.11)</td></tr><tr><td>G5.35+0.10</td><td>17 57 29.2160</td><td>-24 16 57.670</td><td>2.11 (±0.31)</td><td>11.21 (±0.03)</td><td>1.48 (±0.07)</td><td>1.73 (±0.13)</td></tr><tr><td>G5.64+0.24</td><td>17 57 34.6840</td><td>-23 58 01.950</td><td>2.2 (±0.13)</td><td>7.94 (±0.01)</td><td>1.61 (±0.03)</td><td>4.22 (±0.12)</td></tr><tr><td>G5.62-0.08</td><td>17 58 44.6080</td><td>-24 08 39.020</td><td>3.14 (±0.14)</td><td>-25.74 (±0.02)</td><td>2.43 (±0.03)</td><td>3.24 (±0.11)</td></tr><tr><td>G6.22-0.05</td><td>17 59 55.2530</td><td>-23 36 10.150</td><td>3.13 (±0.63)</td><td>19 (±0.09)</td><td>2.78 (±0.15)</td><td>0.62 (±0.09)</td></tr><tr><td>G6.19-0.12</td><td>18 00 07.4100</td><td>-23 39 51.830</td><td>1.5 (±0.41)</td><td>10.88 (±0.06)</td><td>1.99 (±0.14)</td><td>0.99 (±0.12)</td></tr><tr><td>G5.89-0.29</td><td>18 00 07.5280</td><td>-24 00 29.480</td><td>2.09 (±0.23)</td><td>10.11 (±0.02)</td><td>1.42 (±0.04)</td><td>1.88 (±0.1)</td></tr><tr><td>G5.89-0.32</td><td>18 00 14.4580</td><td>-24 01 23.810</td><td>0.9 (±0.13)</td><td>10.11 (±0.02)</td><td>1.74 (±0.05)</td><td>2.6 (±0.09)</td></tr><tr><td>G6.25-0.12</td><td>18 00 16.1890</td><td>-23 37 09.950</td><td>1.71 (±0.25)</td><td>14.48 (±0.06)</td><td>3.31 (±0.11)</td><td>1.02 (±0.08)</td></tr><tr><td>G5.83-0.40</td><td>18 00 24.8380</td><td>-24 06 50.090</td><td>0.91 (±0.21)</td><td>7.64 (±0.03)</td><td>1.67 (±0.07)</td><td>1.54 (±0.08)</td></tr><tr><td>G5.89-0.39</td><td>18 00 30.3380</td><td>-24 04 00.209</td><td>0.69 (±0.09)</td><td>9.18 (±0.03)</td><td>3.96 (±0.06)</td><td>4.05 (±0.13)</td></tr><tr><td>G5.90-0.43</td><td>18 00 40.3520</td><td>-24 04 18.700</td><td>1.8 (±0.06)</td><td>6.64 (±0.01)</td><td>2.53 (±0.03)</td><td>3.8 (±0.09)</td></tr><tr><td>G5.83-0.51</td><td>18 00 51.2040</td><td>-24 10 23.290</td><td>1.99 (±0.12)</td><td>16.22 (±0.01)</td><td>1.74 (±0.03)</td><td>2.8 (±0.08)</td></tr><tr><td>G6.19-0.36</td><td>18 01 01.7140</td><td>-23 47 05.130</td><td>2.11 (±0.12)</td><td>-33.81 (±0.02)</td><td>2.2 (±0.04)</td><td>3.72 (±0.12)</td></tr><tr><td>G5.91-0.54</td><td>18 01 08.2520</td><td>-24 07 14.800</td><td>2.09 (±0.22)</td><td>14.79 (±0.02)</td><td>1.37 (±0.05)</td><td>2.32 (±0.11)</td></tr></tbody></table>

**Notes.** Errors are given in parentheses. The full table is available at CDS. (1) Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds.

<table><caption>**Table 2.** NH₃(2, 2) and (3, 3) line parameters.</caption><thead><tr><th>Name</th><th>ν(2, 2)<br>(km s<sup>-1</sup>)</th><th>Δν(2, 2)<br>(km s<sup>-1</sup>)</th><th>T<sub>MB</sub>(2, 2)<br>(K)</th><th>ν(3, 3)<br>(km s<sup>-1</sup>)</th><th>Δν(3, 3)<br>(km s<sup>-1</sup>)</th><th>T<sub>MB</sub>(3, 3)<br>(K)</th></tr></thead><tbody><tr><td>G5.39+0.19</td><td>10.6 (±0.13)</td><td>1.53 (±0.38)</td><td>0.47 (±0.12)</td><td>—</td><td>—</td><td>&lt;0.09</td></tr><tr><td>G5.35+0.10</td><td>10.85 (±0.17)</td><td>1.8 (±0.35)</td><td>0.38 (±0.11)</td><td>—</td><td>—</td><td>&lt;0.08</td></tr><tr><td>G5.64+0.24</td><td>7.95 (±0.02)</td><td>1.89 (±0.06)</td><td>2.88 (±0.12)</td><td>8.29 (±0.09)</td><td>2.81 (±0.2)</td><td>0.82 (±0.09)</td></tr><tr><td>G5.62-0.08</td><td>-25.57 (±0.05)</td><td>2.93 (±0.11)</td><td>1.93 (±0.14)</td><td>-25.47 (±0.05)</td><td>3.54 (±0.13)</td><td>1.85 (±0.11)</td></tr><tr><td>G6.22-0.05</td><td>—</td><td>—</td><td>&lt;0.1</td><td>—</td><td>—</td><td>&lt;0.08</td></tr><tr><td>G6.19-0.12</td><td>10.63 (±0.26)</td><td>2.91 (±0.62)</td><td>0.3 (±0.1)</td><td>—</td><td>—</td><td>&lt;0.1</td></tr><tr><td>G5.89-0.29</td><td>10.13 (±0.06)</td><td>1.65 (±0.16)</td><td>0.85 (±0.1)</td><td>—</td><td>—</td><td>&lt;0.1</td></tr><tr><td>G5.89-0.32</td><td>10.09 (±0.06)</td><td>1.87 (±0.14)</td><td>1.21 (±0.13)</td><td>10.04 (±0.06)</td><td>1.51 (±0.23)</td><td>1.01 (±0.1)</td></tr><tr><td>G6.25-0.12</td><td>13.97 (±0.2)</td><td>3.02 (±0.42)</td><td>0.44 (±0.12)</td><td>14.35 (±0.39)</td><td>4.06 (±1.18)</td><td>0.31 (±0.1)</td></tr><tr><td>G5.83-0.40</td><td>7.77 (±0.11)</td><td>2.48 (±0.3)</td><td>0.67 (±0.11)</td><td>—</td><td>—</td><td>&lt;0.09</td></tr><tr><td>G5.89-0.39</td><td>9.28 (±0.04)</td><td>4.58 (±0.09)</td><td>3.37 (±0.15)</td><td>—</td><td>—</td><td>—</td></tr><tr><td>G5.90-0.43</td><td>6.58 (±0.04)</td><td>3.29 (±0.08)</td><td>2.92 (±0.14)</td><td>6.6 (±0.05)</td><td>3.87 (±0.12)</td><td>1.76 (±0.09)</td></tr><tr><td>G5.83-0.51</td><td>16.16 (±0.06)</td><td>2.35 (±0.14)</td><td>1.13 (±0.1)</td><td>—</td><td>—</td><td>&lt;0.09</td></tr><tr><td>G6.19-0.36</td><td>-33.84 (±0.04)</td><td>2.83 (±0.1)</td><td>2.71 (±0.14)</td><td>-33.98 (±0.04)</td><td>3.41 (±0.1)</td><td>1.82 (±0.09)</td></tr><tr><td>G5.91-0.54</td><td>14.9 (±0.05)</td><td>1.59 (±0.13)</td><td>1.32 (±0.13)</td><td>14.43 (±0.33)</td><td>5.19 (±0.79)</td><td>0.29 (±0.1)</td></tr></tbody></table>

**Notes.** Errors are given in parentheses. The full table is available at CDS.

### 3.1. Line parameters derived from ammonia

Of the 862 observed sources, 752 were detected in NH₃ (1, 1) (87%), 710 clumps (82%) in the (2, 2) line and 415 sources (48%) in NH₃ (3, 3) with a S/N > 3. The hyperfine structure of most (1, 1) lines is clearly detected, while that of the (2, 2) and (3, 3) lines is mostly too weak to be visible. The sample contains sources of different intensities; examples of them are shown in Fig. 2. There are many strong sources such as G5.62–0.08 with an S/N of 30. An example for a weak line is the (3, 3) line in G10.99–0.08. Its S/N is 3, its fit is shown together with that of some other sources in Fig. 2. The comparison of peak intensities using Tables 1 and 2 yields an average peak intensity of the (2, 2) line of 53% of the average (1, 1) peak line intensity and an average (3, 3) main beam brightness temperature of 40% of the average (1, 1) peak intensity.

Moreover, a few sources like G10.36–0.15 and G23.27–0.26 (see Fig. 2) exhibit two different velocity components. Those are probably two NH₃ clumps, that lie on the same line of sight to the observer with different velocities, e.g. 11 km s⁻¹ and 44 km s⁻¹

for G10.36–0.15, at different distances. The two main components of the inversion transitions can be separated, while the satellite lines partly interfere with each other. These sources were fitted with a model with two velocity components.

### 4. Results and analysis

We first plot different correlations of line parameters, which are derived from the NH₃ spectra observed with the Effelsberg 100 m telescope (FWHM beamwidth = 40′′) (Figs. 3–13). An ammonia (3, 3) maser line together with the thermal emission is displayed in Fig. 14. Then, ammonia results are compared with the 870 µm flux integrated over the Effelsberg beam, extracted from submillimeter dust continuum maps of the ATLASGAL project (Schuller et al. 2009) (Figs. 15 and 16), and with line parameters from the ¹³CO (1–0) emission detected by the Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey (GRS Jackson et al. 2006). Figure 17 compares NH₃ (1, 1) and ¹³CO (1–0) lines. Figure 18 shows our analysis of ¹³CO and NH₃ line-centre velocities, both molecular lines for
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** LSR velocities of detected sources are plotted against the Galactic longitude with CO (1–0) emission (Dame et al. 2001), which is shown in the background. The straight line denotes the 5 kpc molecular ring (Simon et al. 2006b).</figcaption></figure>

<table><caption>**Table 3.** Parameters derived from the NH₃(1, 1) to (3, 3) inversion transitions.</caption><thead><tr><th>Name</th><th>T<sub>rot</sub><br>(K)</th><th>T<sub>kin</sub><br>(K)</th><th>log<sub>10</sub>(N<sub>NH<sub>3</sub></sub>)<br>(cm<sup>-2</sup>)</th></tr></thead><tbody><tr><td>G5.39+0.19</td><td>13.12 (±1.27)</td><td>14.39 (±1.67)</td><td>15.12 (±0.08)</td></tr><tr><td>G5.35+0.10</td><td>11.66 (±1.15)</td><td>12.52 (±1.43)</td><td>15.23 (±0.07)</td></tr><tr><td>G5.64+0.24</td><td>19.53 (±0.82)</td><td>23.94 (±1.39)</td><td>15.34 (±0.03)</td></tr><tr><td>G5.62-0.08</td><td>16.08 (±0.81)</td><td>18.5 (±1.19)</td><td>15.63 (±0.02)</td></tr><tr><td>G6.22-0.05</td><td>—</td><td>—</td><td>—</td></tr><tr><td>G6.19-0.12</td><td>13.92 (±1.99)</td><td>15.45 (±2.7)</td><td>15.21 (±0.12)</td></tr><tr><td>G5.89-0.29</td><td>15.45 (±1.02)</td><td>17.58 (±1.47)</td><td>15.21 (±0.05)</td></tr><tr><td>G5.89-0.32</td><td>18.16 (±1.08)</td><td>21.67 (±1.73)</td><td>14.97 (±0.07)</td></tr><tr><td>G6.25-0.12</td><td>15.79 (±2.13)</td><td>18.07 (±3.09)</td><td>15.5 (±0.07)</td></tr><tr><td>G5.83-0.40</td><td>17.51 (±1.53)</td><td>20.66 (±2.38)</td><td>14.95 (±0.1)</td></tr><tr><td>G5.89-0.39</td><td>27.33 (±1.33)</td><td>39.65 (±3.21)</td><td>15.36 (±0.06)</td></tr><tr><td>G5.90-0.43</td><td>22.69 (±1.15)</td><td>29.63 (±2.22)</td><td>15.5 (±0.03)</td></tr><tr><td>G5.83-0.51</td><td>14.82 (±0.64)</td><td>16.68 (±0.89)</td><td>15.28 (±0.03)</td></tr><tr><td>G6.19-0.36</td><td>20.85 (±1.17)</td><td>26.22 (±2.08)</td><td>15.48 (±0.03)</td></tr><tr><td>G5.91-0.54</td><td>17.44 (±1.22)</td><td>20.55 (±1.89)</td><td>15.22 (±0.05)</td></tr></tbody></table>

**Notes.** Errors are given in parentheses. The full table is available at CDS.

extreme sources of this plot are displayed in Fig. 19. NH₃ (1, 1) and ¹³CO (1–0) linewidths are compared in Fig. 20.

## 4.1. Ammonia velocities

The distribution of NH₃ velocities of the ATLASGAL clumps with Galactic longitude is shown in Fig. 3. For most clumps, the velocities of the ammonia (1, 1) line lie between −10 km s⁻¹ and 120 km s⁻¹, only few exhibit extreme velocities of between −20 km s⁻¹ and −30 km s⁻¹ and between 120 km s⁻¹ to 150 km s⁻¹. They are compared with CO (1–0) emission, observed with the CfA 1.2 m telescope (Dame et al. 2001), which is shown in the background. Most clumps are strongly correlated with CO emission, which probes the larger giant molecular clouds, but some with extreme velocities seem to be related to only weak CO emission. These probably consist of compact clumps, which are detected with the Effelsberg beamwidth of 40”, but are not visible in the larger CfA telescope beam of 8.4’. The straight line indicates the 5 kpc molecular ring that represents the most massive concentration of molecular gas and

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Number distribution of northern ATLASGAL sources with Galactocentric distance.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Correlation plot of the NH₃ (1, 1) and (2, 2) linewidths. The straight green dashed line corresponds to equal widths, the red solid line includes the correction for the hyperfine structure of the (2, 2) line $\Delta v_{\text{Gauss}}(2, 2) = 1.15 \cdot \Delta v_{\text{HFS}}(2, 2) + 0.13$. The (1, 1) and (2, 2) linewidths are mostly equally distributed around the red solid line, as shown by the contour plot in the top panel. For this upper plot we counted the number of sources in each (1, 1) and (2, 2) linewidth bin of 0.1 km s⁻¹. The contours give 10 to 90% in steps of 10% of the peak source number per bin, these levels are used for all contour plots in this article.</figcaption></figure>

star formation activity in the Milky Way (Simon et al. 2006b). Its high CO intensity results from integration over many molecular clouds, which are identified by ¹³CO observations of the Galactic
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** NH₃ (1, 1) main beam brightness temperature plotted against the (1, 1) linewidth. The *upper panel* displays the contour plot with a binning of 1 km s⁻¹ for the (1, 1) linewidth and of 1 K for the (1, 1) temperature, while the *lower panel* shows the scatter plot.</figcaption></figure>

<caption>**Fig. 7.** NH₃ (1, 1) optical depth compared with the (1, 1) linewidth. Contour lines of both parameters are plotted in the *upper panel* with the range of (1, 1) linewidths divided into bins of 1 km s⁻¹ and that of (1, 1) optical depths into bins of 1.</caption>

Ring Survey (GRS, Simon et al. 2006b). For the ATLASGAL sample the NH₃ velocities are crucial to obtain kinematical distances to the clumps, which were calculated using the revised rotation parameters of the Milky Way presented by Reid et al. (2009). We have derived near and far distances (cf. Table 4) and we will distinguish between them by comparing ATLASGAL with extinction maps and absorption in HI lines (Wienen et al., in prep.). Sources associated with infrared extinction peaks are most likely at the near distance. Hence, the addition of the ammonia measurements to the ATLASGAL data reveals the three-dimensional distribution of massive star forming clouds in the first quadrant of our Galaxy. Some clumps are located within GRS molecular clouds, whose kinematic distances have already been determined by using absorption features in the HI 21 cm line against continuum emission (Roman-Duval et al. 2009). These known distances are marked by a star in Table 4. In addition, sources located near the tangential points have similar near and far distances, which are marked by a star in Table 4.

The distribution of Galactocentric distances can be studied independently of near/far kinematic distance ambiguity. Figure 4 shows an enhancement of the number of ATLASGAL sources at Galactocentric radii of 4.5 and 6 kpc, where the Scutum arm and the Sagittarius arm are located. This agrees with previous results. Nguyen Luong et al. (2011b) used ATLASGAL and complementary data to show that the region around Galactic longitude

of 30°, where the W43 Molecular Complex is located, has a very high-mass concentration of sources and star formation activity at a Galactocentric radius of 4.5 kpc. The two peaks in Fig. 4 are also detected in Bronfman et al. (2000), based on a CS(2–1) survey of IRAS point-like sources with FIR colours of UCHII regions (see Fig. 3 and Table 2). Moreover, the same features of Galactic structure were discovered already by the first H II region radio recombination line survey (Reifenstein et al. 1970) and more recently in Anderson & Bania (2009) in a sample of H II regions over the extent of the GRS (see Fig. 9).

## 4.2. Linewidth

The NH₃ (1, 1) inversion transitions exhibit widths from 0.7 km s⁻¹ to 6.5 km s⁻¹ (cf. Fig. 5, *upper panel*), while those of the (2, 2) lines range up to 7.5 km s⁻¹. Linewidths narrower than 0.7 km s⁻¹ cannot be resolved with the backend spectral resolution. One of the clumps with broad linewidths is G10.47+0.03 with Δv(1, 1) = 6.4 km s⁻¹ and Δv(2, 2) = 6.7 km s⁻¹, which is a bright, well-known hot core (Cesaroni et al. 1992). Because the hyperfine structure of the (2, 2) line of most clumps could not be detected, we determined the linewidth by a Gaussian fit. The NH₃(1, 1) and (2, 2) linewidths are compared in the *lower panel* of Fig. 5, the straight dashed
<figure><img src="image_8.png" /><figcaption>**Fig. 8.** Dependence of the rotational temperature between the (1, 1) and (2, 2) inversion transition on the NH₃(1, 1) linewidth is shown as contour plot in the top panel and as a scatter plot in the lower panel. The binning of the (1, 1) linewidth in the contour plot is 1 km s⁻¹ and that of the rotational temperature is 2 K.</figcaption></figure>

<caption>**Fig. 9.** Correlation plot of the logarithm of the column density and the (1, 1) linewidth. The contour plot is illustrated in the upper panel. For the logarithm of the column density bins of 0.4 cm⁻² and for the (1, 1) linewidth bins of 1 km s⁻¹ are chosen.</caption>

line in green shows equal values. In addition, we created a contour plot by dividing the (1, 1) and (2, 2) linewidth ranges into bins of 0.1 km s⁻¹ and counted the number of sources in each bin (upper panel of Fig. 5). The widths of the (2, 2) transition, fitted by a Gaussian, are on average broader than those of the (1, 1) line. The fitting of the (2, 2) transition by a Gaussian does not take optical depths effects and the underlying hyperfine structure splitting into account. For the (2, 2) line the Gaussian fits can therefore result in broader linewidths than for the (1, 1) transition, for which we used a hyperfine structure fit. To investigate this difference, we used a subsample of southern ATLASGAL sources with higher S/N data (Wienen et al., in prep.), for which we were able to detect the hyperfine structure of the (2, 2) transition, and compared their linewidths derived from Gaussian and hyperfine structure fits. This resulted in the correlation $\Delta v_{\text{Gauss}}(2, 2) = 1.15 \cdot \Delta v_{\text{HFS}}(2, 2) + 0.13$, which is a better fit to the data as shown by the solid red line in Fig. 5. The contribution to the line broadening from the optical depth of the line is 15% with a (2, 2) optical depth of 1.24, an additional broadening results from the magnetic hyperfine structure of the (2, 2) line, which is only of about 0.1 km s⁻¹ and therefore cannot be resolved. The contour plot in the top panel of Fig. 5 shows that the contour lines are distributed equally around the red line. However, the lower panel reveals

that some sources still lie above the red line, which indicates broader (2, 2) than (1, 1) linewidths. Because we assumed that the same volume of gas emits the (1, 1) and (2, 2) inversion lines, the same beam filling factor was used for all transitions in calculating the temperature. The difference in the (1, 1) and (2, 2) linewidths now shows that some measured NH₃ lines do not exactly trace the same gas and using equal beam filling factors for those is therefore only an approximation.

The widths of the NH₃ inversion lines consist of a thermal and a non-thermal contribution. The thermal linewidth is given by

$$ \Delta v_{\text{therm}} = \sqrt{\frac{8 \ln(2) k T_{\text{kin}}}{m_{\text{NH}_3}}} \quad (2) $$

with the kinetic temperature $T_{\text{kin}}$, the Boltzmann constant k and the mass of the ammonia molecule $m_{\text{NH}_3}$. Assuming a mean $T_{\text{kin}}$ of 20 K, the thermal linewidth is about 0.22 km s⁻¹. On average measured NH₃(1, 1) linewidths are ~2 km s⁻¹, i.e. they are dominated by non-thermal contributions.

In comparison to Jijina et al. (1999), who measured NH₃ lines of low-mass cores and derived a mean linewidth of about 0.74 km s⁻¹, the linewidths of the ATLASGAL sources are broader on average. The beamwidth of the Effelsberg 100 m
<figure><img src="image_9.png" /><figcaption>**Fig. 10.** Logarithm of the column density compared with the kinetic temperature. There is no correlation between both. Contour lines of the NH₃ parameters are shown in the upper plot. Bins of 0.4 cm⁻² are used for the logarithm of the column density and bins of 2 K for the kinetic temperature.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 11.** Histogram of beam filling factor computed according to Eq. (7). Most sources have a beam filling factor of 0.1.</figcaption></figure>

telescope is 40" at frequencies of the NH₃ lines, corresponding to a linear scale of about 0.8 pc at a typical near distance of our sources of ~4 kpc. Hence, several cores might fill the Effelsberg beam (see discussion of beam filling factors in Sect. 4.5). Assuming that each core possesses a narrow

linewidth, they may add up to the observed linewidths due to a velocity dispersion of the cores within the telescope beam. High-resolution measurements using interferometers are important for investigating individual cores, because they might reveal possible non-thermal contributions to the linewidth of each of them.

To characterize the ammonia sample, Figs. 6 and 7 illustrate the main beam brightness temperature and the optical depth of the NH₃ (1, 1) lines plotted against the (1, 1) linewidth. The binning used for the (1, 1) linewidth is 1 km s⁻¹, the (1, 1) peak intensity values are divided into bins of 1 K and the (1, 1) optical depth into bins of 1 as well. The (1, 1) peak intensity ranges between 0.2 and 6 K with a peak at 1 K, as displayed by the contour plot. Values of the (1, 1) optical depth lie between 0.5 and 5.5 with an average error of 14% and a peak at 1.8. There are no clear trends seen in the (1, 1) main beam brightness temperature and the optical depth.

### 4.3. Rotational temperature

The rotational temperature between the (1, 1) and (2, 2) inversion transition can be determined via the relation (Ho & Townes 1983)

$$T_{\text{rot}} = \frac{-41.5}{\ln\left(\frac{-0.282}{\tau_{\text{m}}(1, 1)} \ln\left(1 - \frac{T_{\text{MB}}(2, 2)}{T_{\text{MB}}(1, 1)}(1 - \exp(-\tau_{\text{m}}(1, 1)))\right)\right)} \quad (3)$$
<figure><img src="image_11.png" /><figcaption>**Fig. 13.** Ratio of observed to calculated (3,3) main beam brightness temperatures compared to the ratio of the (3,3) to (1,1) linewidths as a contour plot in the *top panel* and as a scatter plot in the *middle panel*. The correlation plot of the rotational temperature and observed to calculated (3,3) main beam brightness temperature ratio is displayed in the *lowest panel*. Red points show sources that are not detected in the NH₃ (3,3) line. Green dashed curves indicate very low hot core beam filling factors between 0.008% and 0.15%.</figcaption></figure>

with the optical depth of the (1, 1) main line, $\tau_m(1, 1)$, and the main beam brightness temperatures, $T_{MB}$, of the (1, 1) and (2, 2) inversion transitions. The NH₃ lines of most sources are well fitted, although the maximum intensity of some lines is slightly underestimated by the fit. The rotational temperature lies

<caption>**Fig. 14.** NH₃ (1, 1) to (3, 3) lines of two ATLASGAL clumps, G25.82-0.18 at the *top* and G30.72-0.08 at the *bottom*: the (3, 3) line profile is different from those of the (1, 1) and (2, 2) transitions, indicating a narrow maser line in addition to the thermal emission.</caption>

between 10 K and 28 K with an average error of 9% (Table 3) and a peak at ~16 K; it is plotted against the (1, 1) linewidth in Fig. 8. We chose only sources with reasonably well-determined rotational temperatures (errors smaller than ~50%) for the correlation plot. There is a trend of increasing rotational temperature with broader width of the (1, 1) lines. This, together with the broadening of the (2, 2) lines, which probe higher temperature regions compared to the (1, 1) lines, shows an increase of turbulence with temperature and is an indication for star formation feedback such as stellar winds or outflows. In contrast, some clumps exhibit similar (2, 2) and (1, 1) linewidths and low rotational temperatures. Those sources could be prestellar clumps without any star formation yet and therefore low temperatures. There are also some clumps with broad (1, 1) lines and low rotational temperatures that have large virial masses. The virial mass is a measure of the kinetic energy of a cloud (cf. Sect. 5.2). Some examples of the most extreme sources in Fig. 8 are, e.g. the clump with an extremely narrow (1, 1) linewidth and/or low rotational temperature, G34.37−0.66 with $\Delta v(1, 1) = 0.78$ km s⁻¹ and $T_{rot} = 10$ K, a known IRDC (Peretto & Fuller 2009). Among the clumps with broad linewidth are well-known sources, such as the UCHIR/hot core G10.47+0.03 (Cesaroni et al. 1992; see Sect. 4.2) with $\Delta v(1, 1) = 6.38$ km s⁻¹ and $T_{rot} = 21$ K. Another
<figure><img src="image_12.png" /><figcaption>**Fig. 15.** Correlation plot of the ammonia column density and H₂ column density. Straight lines indicate a high ammonia abundance $\chi_{NH_3} > 3 \times 10^{-7}$, an average value of ~1.2 × 10⁻⁷ and a low NH₃ abundance of 5 × 10⁻⁸. As the contour plot illustration the same correlation is shown in the *upper panel*, the range of NH₃ column densities is divided into bins of 10¹⁵ cm⁻² and that of H₂ column densities into bins of 10²² cm⁻².</figcaption></figure>

source with a high $T_{rot}$ of 26 K and $\Delta v(1, 1)$ of 5.06 km s⁻¹ is G43.80–0.13, well-known for harbouring OH (Braz & Epchtein 1983) and H₂O maser emission (Lekht 2000). G12.21–0.10 has already been observed with the VLA at 5 GHz by Becker et al. (1994), who identified it as a possible UCHII region. Hence, this clump is probably in a more evolved phase of high-mass star formation and therefore has a high $\Delta v(1, 1)$ of 5.73 km s⁻¹ as well and $T_{rot}$ of 22 K, it is also associated with an IRAS point source (Zoonematkermani et al. 1990).

### 4.4. Source-averaged NH₃ column density and kinetic temperature

To calculate the source-averaged ammonia column density, the optical depth and linewidth of the (1, 1) inversion transition are needed as well as the rotational temperature, which is derived from Eq. (3). Because the optical depth and the rotational temperature depend only on line ratios, the resulting column density is a source-averaged quantity. The values of the rotational temperature, kinetic temperature, and the logarithm of the NH₃ column density are given in Table 3. The total column density was derived from the column density of the (1, 1) level ($N(1, 1)$) assuming that the energy levels are populated via a Boltzmann

<figure><img src="image_13.png" /><figcaption>**Fig. 16.** Dynamical masses plotted against gas masses for our subsamples, which are associated with GRS clouds (black points) or located at tangential points (red points). The solid black line denotes equal gas and dynamical masses, while the dashed green fit shows $M_{gas} = 2M_{dyn}$. Note that the gas masses are uncertain by a factor 2 due to uncertain dust emissivity.</figcaption></figure>

distribution, defined by the rotational temperature. Another assumption is that only the four lowest metastable levels are populated because of the low rotational temperatures of the cloud clump, which is ranging from 10 to 25 K, and H₂ densities of 10⁵ cm⁻³. Hence, the total column density is given by Rohlf & Wilson (2004)

$$N_{\text{tot}} \approx N(1, 1) \left( \frac{1}{3} \exp \left( \frac{23.1}{T_{\text{rot}}(1, 2)} \right) + 1 + \frac{5}{3} \exp \left( -\frac{41.2}{T_{\text{rot}}(1, 2)} \right) + \frac{14}{3} \exp \left( -\frac{99.4}{T_{\text{rot}}(1, 2)} \right) \right) \quad (4)$$

with

$$N(1, 1) = 2.07 \times 10^3 \frac{g_1 \nu^2 T_{\text{rot}}}{g_u A_{\text{ul}}} \tau \Delta v, \quad (5)$$

where $g_1$ and $g_u$ are the statistical weights of the lower and upper inversion level, $A_{ul}$ the Einstein A coefficient, $\nu$ the inversion frequency in GHz and $\Delta v$ the linewidth in km s⁻¹. Weak sources, whose errors of the column density are greater than 50%, were excluded from correlation plots (cf. Fig. 9).

The kinetic temperature is obtained from the rotational temperature (Tafalla et al. 2004):

$$T_{\text{kin}} = \frac{T_{\text{rot}}(1, 2)}{1 - \frac{T_{\text{rot}}(1, 2)}{42 \text{ K}} \ln \left( 1 + 1.1 \exp \left( -\frac{16 \text{ K}}{T_{\text{rot}}(1, 2)} \right) \right)}, \quad (6)$$

where the energy difference between the NH₃ (1, 1) and (2, 2) inversion levels is 42 K. Tafalla et al. (2004) ran different Monte Carlo models to compare their derived rotational temperatures with the kinetic temperatures used in the models. Equation (6) is derived from their fit of the kinetic and rotational temperatures. These authors predicted kinetic temperatures in the range of 5 K to 20 K, in which most of our sources fall with a peak of $T_{kin}$ at 19 K, to better than 5%. The logarithm of the averaged NH₃ column density is compared with the (1, 1) linewidth in Fig. 9 and with the kinetic temperature in Fig. 10. The column density ranges between $4 \times 10^{14}$ cm⁻² and $10^{16}$ cm⁻² with an average error of 17% (cf. Table 3) and a peak at
<figure><img src="image_14.png" /><figcaption>Fig. 17. NH₃(1, 1) inversion transition in red overlaid on the ¹³CO(1–0) line in black of a few sources. Some clumps exhibit similar widths of the main line of the NH₃ (1, 1) inversion transitions and ¹³CO lines, while most ¹³CO linewidths are broader than those of the NH₃ (1, 1) transitions.</figcaption></figure>

2 × 10¹⁵ cm⁻². The source with the highest column density of 10¹⁶ cm⁻² is G10.47+0.03 (Cesaroni et al. 1992). The range of column densities of ATLASGAL sources compares well with values found for UCHIIRs by Wood & Churchwell (1989b) and HMPOs, that are considered to be in a pre-UCHII region phase (Molinari et al. 2002; Beuther et al. 2002). It is somewhat above that of low-mass, dense cores from a NH₃ database (Jijina et al. 1999), whose column densities lie between 10¹⁴ cm⁻² and 3 × 10¹⁵ cm⁻².

## 4.5. Beam filling factor

The beam filling factor, η, which gives the fraction of the beam filled by the observed source, is derived from equations of radiative transfer

$$ \eta = \frac{T_{\text{MB}}(1, 1)}{(T_{\text{ex}} - T_{\text{bg}}) \cdot (1 - \exp(-\tau(1, 1)))}, \quad (7) $$

where T<sub>MB</sub>(1, 1) is the main beam brightness temperature of the (1, 1) inversion transition, τ(1, 1) the optical depth of the (1, 1) main line and T<sub>bg</sub> = 2.73 K. Because the H₂ density of the ATLASGAL clumps is ~10⁵ cm⁻³ (Beuther et al. 2002; Motte et al. 2003), we assumed local thermodynamic equilibrium (LTE) and accordingly used the kinetic temperature as excitation temperature in Eq. (7). We found filling factors mostly in the range between 0.02 and 0.6 with one higher value of ~0.9 (cf. Table 4). Their distribution is shown as a histogram in Fig. 11. Equation (7) can lead to underestimation of the filling factor in cases of sub-thermal excitation. To investigate the deviation of the excitation temperature from the kinetic temperature, we used the non LTE molecular radiative transfer program RADEX (van der Tak et al. 2007). For a background temperature of 2.73 K, an average kinetic temperature of our sample of 19 K,

an H₂ density of 10⁵ cm⁻³, a column density of 10¹⁵ cm⁻² and a linewidth of 1 km s⁻¹ we obtained for para-NH₃ an excitation temperature between the lower and upper (1, 1) inversion level of ~18 K. Hence, the excitation and kinetic temperature agree within 94%, which does not hint at sub-thermal excitation. Since typical ATLASGAL source sizes, obtained from Gaussian fits of the 870 µm dust continuum flux, are ~40'' (cf. Sect. 5.2), low beam filling factors indicate clumpiness on smaller scales, as is also evident from interferometer NH₃ observations (e.g. Devine et al. 2011; Ragan et al. 2011; Olmi et al. 2010).

## 4.6. NH₃ (3, 3) line parameters

In Fig. 12 we plot the widths of the ammonia (3, 3) line, Δν(3, 3), against the (1, 1) line, Δν(1, 1). The straight line shows equal widths. A binning of 0.7 km s⁻¹ is chosen for Δν(1, 1) and of 1 km s⁻¹ for Δν(3, 3) for the contour plot. While the (1, 1) linewidths lie between 1 and 6 km s⁻¹, the widths of the (3, 3) line range from 1 to 10 km s⁻¹ (cf. Table 2) and, with a peak at 3.5 km s⁻¹, are broader than those of the (1, 1) inversion transition for most clumps. Because the excitation of the (3, 3) inversion transition requires a high temperature, the (3, 3) emitting gas is likely heated by embedded already formed or forming stars that generate the turbulence observed in the increased linewidth.

Furthermore, we investigated whether the observed NH₃ (3, 3) temperatures are consistent with the temperatures found from the (1, 1) and (2, 2) lines or whether they probe an additional, embedded hot core component of the clumps. We calculated the expected envelope (3, 3) main beam brightness temperature from equations of radiative transfer under LTE conditions using the results from Sects. 4.3 and 4.4 and compared it to the measured T<sub>MB</sub>(3, 3).
<figure><img src="image_15.png" /><figcaption>**Fig. 18.** Absolute value of the difference in line-centre velocities of $^{13}$CO (1–0) and NH$_3$ lines detected in the ATLASGAL sources plotted against the ammonia linewidth in the *top panel*. The differences of most clumps are similar to the typical NH$_3$ linewidth, which is indicated by the dashed red line, and smaller than the width of the $^{13}$CO transition labelled by the solid green line, the rms of the absolute value of the difference in $^{13}$CO and NH$_3$ (1, 1) velocities is shown by the dashed-dotted line. The *bottom panel* shows the relative number distribution of the ATLASGAL clumps with the velocity differences of both molecules. The $^{13}$CO and NH$_3$ velocities of most sources differ by more than the sound speed at 0.23 km s$^{-1}$.</figcaption></figure>

Sources towards which the NH$_3$ (3,3) line is not detected have an average (1,1) linewidth of 1.8 km s$^{-1}$, which is slightly below the mean value of the whole ATLASGAL sample of 2 km s$^{-1}$. Their average rotational temperature is 15.3 K, which is also lower than the mean of 16.5 K.

Assuming LTE (and thus $T_{\text{kin}} = T_{\text{ex}}$), the (3,3) main beam brightness temperature is

$$T_{\text{MB}}(3, 3)_{\text{calc}} = \eta \cdot T_{\text{kin}}(1 - \exp(-\tau(3, 3))) \quad (8)$$

assuming the same beam filling factor as for the (1,1) line (see Eq. (7)) and the optical depth $\tau(3, 3)$ of the NH$_3$(3,3) line given by

$$\tau(3, 3) = \alpha \times \left(1 - \exp\left(-\frac{h\nu(3, 3)}{kT_{\text{kin}}}\right)\right), \quad (9)$$

where $\nu(3, 3)$ is the frequency of the (3,3) inversion transition (23.87 GHz), $k$ the Boltzmann constant and $T_{\text{kin}}$ the measured kinetic temperatures. $\alpha$ is defined as

$$\alpha = \frac{c^3}{8\pi\nu(3, 3)^3} \frac{2\sqrt{\ln2}}{\sqrt{\pi}} \frac{1}{\Delta\nu(1, 1)} N(3, 3)A_{\text{ul}}(3, 3).$$

<figure><img src="image_16.png" /><figcaption>**Fig. 19.** NH$_3$(1,1) and $^{13}$CO(1–0) lines of G30.68–0.03 are plotted in the *upper panel*, the *lower* one shows these transitions of the clump G37.76–0.22. The Gaussian fit of the spectra is shown.</figcaption></figure>

Here, $c$ denotes the speed of light, $\Delta\nu(1,1)$ the observed (1,1) linewidth, $N(3,3)$ the (3,3) level population number and $A_{\text{ul}}(3,3)$ the Einstein A coefficient of the (3,3) transition. From the Boltzmann distribution of energy levels we determined the (3,3) level population number using derived column densities $N_{\text{tot}}$ (cf. Table 3) (Rohlfs & Wilson 2004)

$$N(3, 3) = N_{\text{tot}}/Q, \quad (10)$$

where $Q$ is given by

$$Q = \frac{1}{14} \exp\left(\frac{123.57}{T_{\text{kin}}}\right) + \frac{3}{14} \exp\left(\frac{99.4}{T_{\text{kin}}}\right) + \frac{5}{14} \exp\left(\frac{59.12}{T_{\text{kin}}}\right) + 1.$$

The ratio of observed to calculated (3,3) main beam brightness temperature is compared with the ratio of (3,3) to (1,1) linewidths in Fig. 13. The ranges of both ratios are divided into bins of 0.4 to create the contour plot. We found values between 0.3 and 5 for the linewidth ratio and from 0.3 to 6 for the (3,3) line temperature ratio. The median of $\Delta\nu(3,3)/\Delta\nu(1,1)$ is 1.6 with a dispersion of 0.8 and the median of $T_{\text{MB}}(3,3)_{\text{obs}}/T_{\text{MB}}(3,3)_{\text{calc}}$ is 1.2 with a dispersion of 1.1. Many clumps have higher observed (3,3) line temperatures than is expected from the calculation using the (1,1) and (2,2) line temperatures, which can have various reasons:

1. Non-thermal excitation can overproduce or invert the upper level of the (3,3) doublet. Different examples of these maser lines have already been observed: Mangum & Wootten (1994) detected $^{14}$NH$_3$ (3,3) maser emission towards the
<figure><img src="image_17.png" /><figcaption>**Fig. 20.** Widths of the ¹³CO (1–0) transition compared to those of the NH₃ (1, 1) line as contour plot in the *upper panel* and as scatter plot in the *lower panel*. The straight line corresponds to equal widths. The binning used for the NH₃ (1, 1) linewidth is 0.4 km s⁻¹ and for the ¹³CO linewidth is 1 km s⁻¹.</figcaption></figure>

DR 21(OH) star forming region. Hofner et al. (1994) discovered the NH₃ (5,5) maser line with interferometric observations towards the H II region complex G9.62+0.19. Some ATLASGAL clumps have (3,3) line profiles, which are different from those of the (1,1) and (2,2) lines and indicate maser emission; two examples are shown in Fig. 14: the (3,3) line of G25.82–0.18 shows two narrow components at lower and higher velocities than v(3, 3) in addition to the thermal emission. This clump also harbours CH₃OH and water masers (Szymczak et al. 2005). In the NH₃ (3,3) spectrum of G30.72–0.08 we see an additional peak at a slightly lower velocity than that of the (3,3) thermal emission. This clump is located in the W43 high-mass star forming region and is revealed as a compact cloud fragment detected at 1.3 mm and 350 µm by Motte et al. (2003). However, they did not find an association of this source with an OH, H₂O or CH₃OH maser, but instead a bright compact H II region with a flux of ~1 Jy at 3.5 cm. Recent Herschel observations show that it is centred on a 1 pc radius IRDC seen in absorption at 70 µm and bright at wavelengths longer than 160 µm (Bally et al. 2010). However, because we did not find strong maser lines, but only weak narrow components in the (3,3) spectra of a small subsample of ATLASGAL clumps, it is

unlikely that population inversion of the (3,3) transition explains the trend of increased observed $T_{MB}$(3, 3) values compared to calculations.

2. Since radiative and collisional transitions let the spin orientations remain constant, transitions between ortho- and para-NH₃ are not allowed. There are processes allowing interconversion between both species proposed by Cheung et al. (1969), but they are very slow (~10⁶ yr). The distribution between ortho- and para states is a strong function of temperature (Takano et al. 2002): the abundance ratio will be about 1 if the formation of NH₃ occurs in high-temperature regions with more than 40 K, but it will rise if NH₃ is produced at low temperatures, because the lowest level (0,0) is an ortho state. Our average rotational temperature is 19 K, for which Takano et al. (2002) calculated an ortho-to-para abundance ratio of 2. This ratio would result in an increased population of the (3,3) level and a decrease of observed to calculated (3,3) line temperature ratio.

3. High $T_{MB}$(3, 3)obs values might originate from a hot embedded source, which is revealed by broad (3,3) linewidths leading to a high ratio of $\Delta v$(3,3) to $\Delta v$(1,1). We can exclude any calibration issue, because the NH₃ (1,1) to (3,3) lines are observed simultaneously and a relative calibration error would be the same for all sources, but a wide range of $T_{MB}$(3, 3)obs/$T_{MB}$(3, 3)calc is found. Since only few typical maser line profiles are observed and the ortho-to-para ratio would only lead to variations of about a factor two and does not explain the increase in linewidth, we favour the hot core explanation. To quantify this effect, we computed the resulting rotational and (3,3) temperatures when adding a hot core with filling factor $\eta_{hc}$ and temperature of 100 K. In Fig. 13 the resulting increase in the (3,3) line temperatures is shown for various filling factors. The effect is diminished for higher rotational temperatures, since then an observable (3,3) line is produced already by the large-scale cold clump alone. Sources that are not detected in the NH₃ (3,3) line are plotted as upper limits (corresponding to three times the rms noise in $T_{MB}$) in red.

# 5. Comparison with the 870 µm dust continuum

The ATLASGAL data were reduced and calibrated as described in Schuller et al. (2009). For the sources of our sample, selected as described in Sect. 2.1, we integrated the 870 µm dust continuum flux within an aperture with a diameter of 40", which corresponds to the beam of the Effelsberg telescope at the ammonia inversion frequencies of ~24 GHz. The integrated dust continuum fluxes $S_{870\mu m}^{40''}$ are given in Table 4.

## 5.1. Submillimeter flux and NH₃

The integrated submillimeter dust continuum flux ranges between 0.6 Jy and 47 Jy and in this section it is compared with the ammonia column density. Because the 870 µm continuum emission is optically thin, the submm flux is proportional to the dust column density at a given temperature. Assuming a constant gas-to-dust mass ratio, the dust column density is proportional to the H₂ density ($N_{H_2}$), which is related to the ammonia column density ($N_{NH_3}$) and abundance ($\chi_{NH_3}$) via the relation $N_{H_2} = N_{NH_3}/\chi_{NH_3}$. Hence, assuming a constant NH₃ abundance, a correlation between the integrated 870 µm flux and the ammonia column density is expected, which is confirmed by this study. However, a few clumps exhibit high NH₃ column densities and
<table><caption>**Table 4.** Parameters from the 870 µm dust continuum and NH₃ lines.</caption><thead><tr><th>Name</th><th>η</th><th>N<sub>H2</sub><br>(10<sup>22</sup> cm<sup>-2</sup>)</th><th>χ<br>(10<sup>-7</sup>)</th><th>S<sub>870 µm</sub><sup>40''</sup><br>(Jy)</th><th>Distance (near)<br>(kpc)</th><th>Distance (far)<br>(kpc)</th></tr></thead><tbody><tr><td>G5.39+0.19</td><td>0.18</td><td>1.18</td><td>1.12</td><td>1.21</td><td>2.8</td><td>13.9</td></tr><tr><td>G5.35+0.10</td><td>0.2</td><td>1.79</td><td>0.95</td><td>1.44</td><td>2.9</td><td>13.8</td></tr><tr><td>G5.64+0.24</td><td>0.22</td><td>3.79</td><td>0.58</td><td>8.42</td><td>2.2</td><td>14.5</td></tr><tr><td>G5.62-0.08</td><td>0.21</td><td>2.59</td><td>1.65</td><td>3.96</td><td>—</td><td>—</td></tr><tr><td>G6.22-0.05</td><td>—</td><td>—</td><td>—</td><td>2.74</td><td>3.7</td><td>13.0</td></tr><tr><td>G6.19-0.12</td><td>0.1</td><td>1.32</td><td>1.23</td><td>1.52</td><td>2.7</td><td>14.1</td></tr><tr><td>G5.89-0.29</td><td>0.14</td><td>2.1</td><td>0.77</td><td>2.98</td><td>2.5</td><td>14.2</td></tr><tr><td>G5.89-0.32</td><td>0.23</td><td>2.15</td><td>0.43</td><td>4.15</td><td>2.6</td><td>14.1</td></tr><tr><td>G6.25-0.12</td><td>0.08</td><td>2.16</td><td>1.46</td><td>3.19</td><td>3.1</td><td>13.6</td></tr><tr><td>G5.83-0.40</td><td>0.14</td><td>1.45</td><td>0.61</td><td>2.62</td><td>2.1</td><td>14.6</td></tr><tr><td>G5.89-0.39</td><td>0.22</td><td>9.74</td><td>0.24</td><td>41.64</td><td>2.4</td><td>14.4</td></tr><tr><td>G5.90-0.43</td><td>0.17</td><td>6.88</td><td>0.46</td><td>20.37</td><td>1.9</td><td>14.8</td></tr><tr><td>G5.83-0.51</td><td>0.23</td><td>3.8</td><td>0.5</td><td>4.96</td><td>3.4</td><td>13.3</td></tr><tr><td>G6.19-0.36</td><td>0.18</td><td>3.73</td><td>0.81</td><td>9.38</td><td>—</td><td>—</td></tr><tr><td>G5.91-0.54</td><td>0.15</td><td>1.37</td><td>1.21</td><td>2.45</td><td>3.3</td><td>13.5</td></tr></tbody></table>

**Notes.** Columns are beam filling factor, H₂ column density, NH₃ abundance, 870 µm dust continuum flux integrated over 40'', near kinematic distance, far kinematic distance. * Are sources with known distances (cf. Sect. 5.2.). The full table is available at CDS.

low submm fluxes, which can be explained by an increased ammonia abundance or an overestimation of the temperature. To investigate the two effects we compared the NH₃ column density with the H₂ column density, computed using the NH₃ kinetic temperatures.

Models of dust grains with thick ice mantles give an absorption coefficient, $\kappa$, of 1.85 cm²/g at 870 µm at a gas density $n(H) = 10^6$ cm⁻³ (Ossenkopf & Henning 1994). We estimated the H₂ column density per the relation (Kauffmann et al. 2008)

$$N(\text{H}_2) = \frac{2.03 \times 10^{16} S_{870 \text{ µm}} \lambda^3 \left( \exp\left(\frac{1.44 \times 10^4}{T_d \lambda}\right) - 1 \right) Z_{\odot}}{\kappa(\lambda)\theta^2} Z, \quad (11)$$

where $S_{870 \text{ µm}}^{40''}$ is the 870 µm flux density in Jy integrated within 40'', $\lambda$ the wavelength in µm and $T_d$ the dust temperature with $T_d = T_{\text{kin}}$ under the assumption of equal gas and dust temperatures. $\theta$ is the 100-m telescope beamwidth of 40'' and $Z/Z_{\odot}$ the ratio of metallicity to the solar metallicity, assuming $Z/Z_{\odot} = 1$. We found that most H₂ column densities lie between $2.6 \times 10^{21}$ and $1.8 \times 10^{23}$ cm⁻² with a peak at $1.7 \times 10^{22}$ cm⁻². For each clump we can now use the ammonia column density $N_{\text{NH}_3}$ to derive the NH₃ abundance via $\chi_{\text{NH}_3} = N_{\text{NH}_3}/N_{\text{H}_2}$, which lies in the range from $7 \times 10^{-9}$ to $1.2 \times 10^{-6}$ with an average of $\sim 1.2 \times 10^{-7}$. The H₂ column densities, NH₃ abundances and the integrated dust continuum fluxes $S_{870 \text{ µm}}^{40''}$ are given in Table 4. NH₃ column densities are plotted against H₂ column densities in Fig. 15. They show a better correlation than NH₃ column densities and the 870 µm fluxes, which are affected by variations of the temperature of the sources. Since the kinetic temperature was used in the calculation of the H₂ column density, Fig. 15 shows the remaining variation in abundance, which can give us a better understanding of the chemical history of the clumps. The scatter, given by the ratio of the rms to the mean value, improved from 0.039 for $S_{870 \text{ µm}}^{40''}/N_{\text{NH}_3}$ to 0.0271 for $N_{\text{H}_2}/N_{\text{NH}_3}$. Deviations from the trend in Fig. 15 are due to a higher ammonia abundance, clumps with $\chi_{\text{NH}_3} > 3 \times 10^{-7}$ are below the lower straight line, or a decreased abundance, values lower than $5 \times 10^{-8}$ are above the upper straight line. We note that the abundances are upper limits since they are rather

derived from the ratio of a source and a beam-averaged column density, hence affected by the filling factor of the cores. With the derived beam filling factors (see Fig. 11) the abundances would on average be higher by about a factor ten. With this uncertainty, we find abundances in the range of $5 \times 10^{-9}$ to $3 \times 10^{-7}$. Dunham et al. (2011) found NH₃ abundances between $2 \times 10^{-9}$ and $5 \times 10^{-7}$, which is consistent with our results. Moreover, Pillai et al. (2006) obtained NH₃ abundances from $10^{-8}$ to $10^{-7}$ for a sample of IRDCs, their $\chi_{\text{NH}_3}$ range is narrower compared to our values and similar to our higher abundances. The average of our derived NH₃ abundances agrees with those predicted for low-mass pre-protostellar cores by chemical models (Bergin & Langer 1997). They show that ammonia is more abundant in cold dense cores than other molecules such as CO, which are depleted from the gas phase at high densities.

## 5.2. Virial masses and gas masses

We can estimate the masses of the clumps using parameters derived from the 870 µm dust continuum and from NH₃ lines. The gas mass of a source obtained from the dust is (Kauffmann et al. 2008)

$$M_{\text{gas}} = 1.2 \times 10^{-14} \frac{S_{870 \text{ µm}}^{40''} \lambda^3 d^2 \left( \exp\left(\frac{1.44 \times 10^4}{T_d \lambda}\right) - 1 \right) Z_{\odot}}{\kappa(\lambda)} Z, \quad (12)$$

with the distance to the source, $d$, in pc; for the other parameters see Eq. (11). The masses can be determined only for the clumps for which we have a preliminary distance estimation: 277 clumps are associated with clouds from the GRS survey with known distances (Roman-Duval et al. 2009), and 71 clumps, which are located at the tangential points, have similar near and far distances, marked in Table 4. To determine distances of the clumps studied by Roman-Duval et al. (2009), we followed their choice of near/far distance and used the ammonia velocity together with the rotation curve of the Milky Way given by Reid et al. (2009). The mass can also be calculated using the NH₃ (1, 1) linewidth under the assumption of virial equilibrium (Rohlfss & Wilson 2004)

$$M_{\text{vir}} = 250\Delta v(1, 1)_{\text{corr}}^2 \times R, \quad (13)$$
with the linewidth $\Delta v(1, 1)_{corr}$, corrected for the resolution of the spectrometer $\Delta v(1, 1)_{corr} = \sqrt{\Delta v(1, 1)^2 - (0.7 \text{ km s}^{-1})^2}$. The radius $R$ in pc is obtained from Gaussian fits of the 870 $\mu$m continuum flux using the Miriad task sfind (Sault et al. 1995) at the distance of the source. Equation (13) is derived from the virial theorem by neglecting the magnetic energy and the virial mass is equal to the dynamical mass. The gas and dynamical masses are compared in Fig. 16, the black points in the scatter plot indicate sources associated with GRS clouds and the red points clumps at the tangential points. Gas masses range from 60 to $1.6 \times 10^4 M_{\odot}$ and virial masses from 20 to $5.5 \times 10^3 M_{\odot}$. The straight solid black line shows equal gas and dynamical masses. Most data lie below that relation, indicating that for most sources their masses exceed their dynamical masses. Since magnetic fields can contribute to prevent the collapse of molecular clouds (Schulz 2005), the magnetic energy has to be taken into account in the virial theorem. We did not measure the magnetic field and assumed equipartition between magnetic and kinetic energy. Then, virialisation is obtained for $M_{\text{gas}} = M_{\text{vir}} = 2M_{\text{dyn}}$. This relation is shown by the dashed green line, but most sources are still located below it. We calculated the virial parameter (Bertoldi & McKee 1992)

$$ \alpha = \frac{M_{\text{vir}}}{M_{\text{gas}}} \qquad (14) $$

which is ~1 for clumps supported against gravitational collapse. When we use $M_{\text{vir}} = M_{\text{dyn}}$ from Eq. (13), sources with narrow linewidths have virial parameters on average much smaller (mean $\alpha = 0.21$) than broad linewidth sources with a mean, $\alpha$, of 0.45. With higher density tracers such as $H^{13}\text{CO}^+$ and $C^{34}\text{S}$ on average broader linewidths are found, leading to higher virial mass estimates. This will be discussed in Wienen et al. (in prep.). Including again the magnetic energy into the virial theorem gives a factor 2 resulting in $\alpha \sim 1$ for clumps with broad linewidth. Therefore, these sources are in virial equilibrium. In contrast, sources with small $\Delta v(1, 1)$ still have a low virial parameter of 0.42.

Hill et al. (2010), who observed $\text{NH}_3$ lines of high-mass star forming regions with the Parkes telescope, found that their sample is mostly in virial equilibrium. However, they included many sources that have associated infrared sources, UCHII regions, and high-mass protostars, but only a few cores in an early evolutionary phase without any hint at ongoing massive star formation. Hill et al. (2010) found a higher average $\text{NH}_3$ (1, 1) linewidth of $2.9 \text{ km s}^{-1}$ in contrast to our mean $\Delta v(1, 1)$ of $2 \text{ km s}^{-1}$, which results in larger virial masses than we found.

Dunham et al. (2010) obtained an average virial parameter of ~1 from their $\text{NH}_3$ observations of high-mass star forming sources, which form a similar sample as ours, only at a near distance of about 2 kpc compared to a typical near distance of our sources of ~4 kpc. However, Dunham et al. (2010) calculated gas and virial masses between 20 and about 1000 $M_{\odot}$ and therefore probed much smaller masses compared to the ATLASGAL clumps, whose gas masses range up to $1.6 \times 10^4 M_{\odot}$.

Dunham et al. (2011) found virial parameters about a factor 2 larger than our study. While the linewidths and dust masses in both studies are comparable, Dunham et al. (2011) used for their virial mass computation an “effective radius” that on average is about twice as large as the half maximum radii that we used from ATLASGAL, which accounts for this difference.

We estimated the magnetic field strength, which is necessary to prevent gravitational collapse, using $B = 2.5 \times N(H_2) \mu G$ (McKee et al. 1993), where $N(H_2)$ is the hydrogen column density in units of $10^{21} \text{ cm}^{-2}$. For our values of the hydrogen column

density between $2.3 \times 10^{21}$ and $1.8 \times 10^{23} \text{ cm}^{-2}$, we obtain magnetic fields between ~5 and 450 $\mu$G.

Estimates of the gas mass and gas column density as calculated in Eq. (11) can vary strongly. The large uncertainty of Eq. (12) is discussed in Beuther et al. (2002) and Motte et al. (2007), who gave a factor 2 of uncertainty in the mass estimates due to uncertain dust emissivity. In addition, Martin et al. (2012) discussed evidence for changes of the dust opacity in different environments and also found observationally a factor 2 of variation but excluding high-density regions. For denser regions, direct observations of dust opacity changes are difficult. Ormel et al. (2011) presented new calculations of dust opacities in dense regions as a function of grain growth and found for moderate time scales ($<10^7$) a factor 2 increase in the dust opacity with time. The effects of the systematic errors of the linewidths that enter Figs. 15 and 16 increase with linewidths, but are much smaller than the uncertainties in the dust opacity. As discussed already in this section, the main uncertainty here is how representative the measured ammonia linewidths are for the dense material of the dust clumps. As we will discuss in more detail in Wienen et al. (in prep.), observations of higher density probes towards ATLASGAL sources show on average broader linewidths, which in turn lead to higher virial mass estimates.

## 6. Correlation of $\text{NH}_3$ and $^{13}\text{CO}$ (1–0) emission

While ammonia is a high-density tracer (~$10^4 \text{ cm}^{-3}$; Ungerechts et al. 1986), the $^{13}\text{CO}$ molecule can be used to probe low-density regions (~$10^3 \text{ cm}^{-3}$; Ungerechts et al. 1997), which are surrounding the dense cores. In this section, the ammonia observations are compared with the $^{13}\text{CO}$ $J = 1 \rightarrow 0$ line emission measure at 110.2 GHz in the course of the Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey (GRS) (Jackson et al. 2006). It has a greater sensitivity of ~0.13 K, a higher spectral resolution of $0.2 \text{ km s}^{-1}$, comparable or better angular resolution (46") and sampling (22") than previous molecular line surveys of the inner Galaxy. $^{13}\text{CO}$ is used instead of the commonly observed $^{12}\text{CO}$ $J = 1 \rightarrow 0$ emission, because $^{13}\text{CO}$ is much less abundant, and consequently possesses an optically thinner transition with narrower linewidths and its analysis can determine the properties of the clumps more reliably. An area of 75.4 deg$^2$ is covered by the whole survey with Galactic longitudes of $l = 18^\circ - 55.7^\circ$, limited data are available for $14^\circ < l < 18^\circ$, and latitudes of $|b| < 1^\circ$. The velocity coverage ranges from -5 to 135 km s$^{-1}$ up to a Galactic longitude of 40$^\circ$ and from -5 to 85 km s$^{-1}$ for $l > 40^\circ$. We extracted $^{13}\text{CO}$ (1–0) lines from GRS intensity images, which are integrated over the velocity (Jackson et al. 2006). The CLASS software was used to derive the $^{13}\text{CO}$ LSR velocity $v_{^{13}\text{CO}}$ of sources detected by GRS and ATLASGAL, their linewidths $\Delta v_{^{13}\text{CO}}$ and integrated temperatures $T_{\text{MB},^{13}\text{CO}}$, each with errors. The results are given in Table 5, $\text{NH}_3$ and $^{13}\text{CO}$ line parameters are compared in correlation plots.

### 6.1. Comparison of ammonia and $^{13}\text{CO}$ lines

We extracted the $^{13}\text{CO}$ (1–0) line from GRS intensity images at the position of the ammonia observations, fitted a Gaussian to the $^{13}\text{CO}$ line profile and related the $^{13}\text{CO}$ line with the velocity closest to the $\text{NH}_3$ velocity to the observed ATLASGAL source. Within the ammonia sample we found 517 ATLASGAL clumps associated with $^{13}\text{CO}$ emission. The spectra of $^{13}\text{CO}$ and $\text{NH}_3$ (1, 1) lines are overlaid, some examples of strong sources are shown in Fig. 17. The $^{13}\text{CO}$ linewidths are similar to the widths
<table><caption>**Table 5.** Line parameters derived from ¹³CO (1–0) transition.</caption><thead><tr><th>Name</th><th>v<sub>13CO</sub><br>(km s<sup>-1</sup>)</th><th>Δv<sub>13CO</sub><br>(km s<sup>-1</sup>)</th><th>T<sub>MB, 13CO</sub><br>(K)</th></tr></thead><tbody><tr><td>G14.62+0.33</td><td>26.73 (±0.03)</td><td>2.56 (±0.09)</td><td>3.26 (±0.17)</td></tr><tr><td>G14.63+0.31</td><td>26.03 (±0.04)</td><td>2.34 (±0.1)</td><td>3.65 (±0.16)</td></tr><tr><td>G14.92+0.07</td><td>25.93 (±0.05)</td><td>3.35 (±0.12)</td><td>3.76 (±0.26)</td></tr><tr><td>G15.01+0.01</td><td>25.38 (±0.13)</td><td>4.35 (±0.42)</td><td>2.65 (±0.15)</td></tr><tr><td>G15.43+0.19</td><td>47.85 (±0.03)</td><td>3.13 (±0.06)</td><td>3.31 (±0.17)</td></tr><tr><td>G14.99-0.12</td><td>49.46 (±0.06)</td><td>3.15 (±0.19)</td><td>1.67 (±0.18)</td></tr><tr><td>G15.46-0.00</td><td>66.92 (±0.06)</td><td>3.3 (±0.13)</td><td>1.66 (±0.12)</td></tr><tr><td>G14.89-0.40</td><td>61.64 (±0.13)</td><td>2.65 (±0.32)</td><td>2.74 (±0.66)</td></tr><tr><td>G15.20-0.44</td><td>20.86 (±0.03)</td><td>2.38 (±0.08)</td><td>3.69 (±0.13)</td></tr><tr><td>G15.53-0.41</td><td>39.27 (±0.07)</td><td>4.36 (±0.18)</td><td>1.97 (±0.19)</td></tr><tr><td>G15.56-0.46</td><td>15.24 (±0.03)</td><td>1.78 (±0.06)</td><td>3.92 (±0.28)</td></tr><tr><td>G16.36-0.07</td><td>47.25 (±0.05)</td><td>2.72 (±0.13)</td><td>5.75 (±0.19)</td></tr><tr><td>G16.58-0.05</td><td>59.41 (±0.02)</td><td>3.27 (±0.06)</td><td>6.53 (±0.17)</td></tr><tr><td>G16.83+0.08</td><td>62.43 (±0.03)</td><td>2.48 (±0.07)</td><td>2.35 (±0.15)</td></tr><tr><td>G16.58-0.08</td><td>40.85 (±0.17)</td><td>3.74 (±0.57)</td><td>2.32 (±0.13)</td></tr></tbody></table>

Notes. Errors are given in parentheses. The full table is available at CDS.

of the NH₃ (1, 1) inversion transitions for a few sources such as G16.91–0.08, G21.39–0.25, G28.82+0.36, G29.60–0.62, while most clumps possess broader ¹³CO than NH₃ (1, 1) linewidths, e.g. G15.43+0.19, G18.82–0.47, G18.88–0.51, G22.37+0.38, G18.83–0.48, and G28.92–0.23. While most ATLASGAL clump lines of sight show only one NH₃ (1, 1) velocity component, the ¹³CO (1–0) often shows multiple components (see e.g. G16.91–0.08, G21.39–0.25, G28.82+0.36, G15.43+0.19, and G22.37+0.38 in Fig. 17). Concerning those clumps, we assumed the ¹³CO line with the velocity closest to the NH₃ velocity to be associated with the ATLASGAL source, which has been demonstrated to be a reasonable choice by Kirk et al. (2007) for low-mass dense cores.

## 6.2. Clump-to-cloud motion

While it is not yet understood how high-mass stars form, two main theoretical models are currently under consideration: (1) in the turbulent core model massive star forming clouds are supported against gravitational collapse mainly by supersonic turbulence (McKee & Tan 2003; Mac Low & Klessen 2004); (2) according to recent simulations of star formation in clusters (e.g. Bonnell et al. 1997, 2001), the accretion rate and in turn the mass of the forming star are strongly affected by the density of the molecular cloud through which it is moving. In addition, material is supposed to diffuse into the potential well of the cloud at the centre of a cluster, where high-mass stars are preferentially forming. This theory of competitive accretion also agrees with observations, which show that the most massive stars of a cluster are located at its centre (Bontemps et al. 2010). According to Walsh et al. (2004), this moving accretion is only possible if the protostellar core is moving at appreciable velocity with respect to its surrounding envelope, which will be associated with a significant difference in the line-centre velocity of high- and low-density tracers. Previous studies of low-mass star forming regions have analysed the motions of dense cores with respect to their surrounding envelope, using N₂H⁺ as high-density tracer and ¹³CO as well as C¹⁸O as probes of less dense gas (e.g. Walsh et al. 2004, 2007; Kirk et al. 2007). While Walsh et al. (2004) investigated mostly single isolated cores, clustered regions of dense cores such as NGC 1333 and the Perseus molecular cloud were also studied by Walsh et al. (2007) and Kirk et al. (2007).

However, they did not find significant motions of the cores relative to their lower density surroundings. In this section we will compare their results to conditions prevailing in high-mass star forming regions.

Although in the more distant high-mass star forming regions we sample clumps on larger scales than the aforementioned studies, it is still interesting to investigate motions of clouds relative to their host clouds. It is especially useful that the NH₃ and ¹³CO lines, good tracers for the former and the latter, respectively, were measured with similar beamwidths, 40" and 46", respectively. The absolute value of the difference in ¹³CO and NH₃ (1, 1) line-centre velocities is plotted against the linewidth of the ammonia (1, 1) inversion transition in the top panel of Fig. 18. If ammonia clumps showed significant motions with respect to their envelopes, the differences in line-centre velocities should be comparable to the width of the ¹³CO lines (Walsh et al. 2004). Deviations in the range of NH₃ linewidths would be associated with random motions within the ammonia clumps. The dashed red line in the top panel of Fig. 18 at 2 km s⁻¹ indicates the average widths of NH₃ (1, 1) lines and the solid green line at ~4.5 km s⁻¹ illustrates the average ¹³CO (1–0) linewidths. The rms of the difference in both velocities is 0.4 km s⁻¹, which is shown by the dashed dotted black line. Because the distribution of almost all clumps is narrower than the average NH₃ linewidth and also narrower than the average width of the ¹³CO (1–0) transition, random motions of the ammonia clumps are expected and hint at turbulence, which could prevail in the sources. Hence, these low relative velocities between the ammonia clump and ¹³CO cloud put strong constraints on the role of turbulence in the clump formation. There are only two sources whose absolute value of NH₃ (1, 1) and ¹³CO velocity difference is slightly higher than the NH₃ linewidth: G30.68–0.03 with a deviation in both velocities of ~2.4 km s⁻¹ has a lower ¹³CO velocity, while G37.76–0.22 exhibits a lower NH₃ velocity with a difference to ¹³CO of 2.24 km s⁻¹. Figure 19 illustrates that the ¹³CO lines of both clumps might show either self-absorption or each ¹³CO spectrum could be emitted by two clumps, that lie on the same line of sight to the observer, leading to higher uncertainties in the ¹³CO velocities. Concerning low-mass star formation, the line-centre velocity differences of N₂H⁺ and ¹³CO or C¹⁸O in cores observed by Walsh et al. (2004, 2007) are also similar to their N₂H⁺ linewidths, which is consistent with the analysis of the ATLASGAL sources. This result indicates that the velocity deviations of both molecular probes in low-mass sources are much lower than in massive ones, because low-mass cores have a much more narrow average N₂H⁺ and ¹³CO linewidth of ~0.38 km s⁻¹ and ~1.18 km s⁻¹ (Walsh et al. 2007) compared to the broader ammonia and ¹³CO lines of high-mass cores. This also agrees with and is additionally supported by the analysis of Kirk et al. (2007), who measured differences in N₂H⁺ and C¹⁸O velocities similar to the N₂H⁺ thermal linewidths of 0.065 km s⁻¹ for most dense cores in the Perseus molecular cloud assuming a temperature of 15 K.

The distribution of differences in line-centre velocities of ¹³CO and NH₃ in the bottom panel of Fig. 18 illustrates that both, higher and lower ¹³CO velocities than ammonia velocities are present. The differences can be compared to the sound speed of the ambient medium, which is 0.23 km s⁻¹ assuming gas that consists of 90% H₂ and 10% He with a mean molecular weight of 2.33 at a temperature of 15 K. Most clumps show differences larger than the sound speed, only 32% of them have subsonic motions. A high percentage of sources shows a deviation in clump-to-cloud velocities between 1 km s⁻¹ and ~0.7 km s⁻¹, some range to 2.2 km s⁻¹ and ~2.9 km s⁻¹. This result can be
<table><caption>**Table 6.** Properties of ATLASGAL subsamples belonging to different evolutionary stages of high-mass star formation.</caption><thead><tr><th>Sample</th><th>Search radius (")</th><th>Number</th><th>T<sub>rot</sub> (K)</th><th>σ(T<sub>rot</sub>)</th><th>Δν(1, 1) (km s<sup>-1</sup>)</th><th>σ(Δν(1, 1))</th></tr></thead><tbody><tr><td>24 µm dark</td><td>40</td><td>98</td><td>15.25</td><td>0.21</td><td>1.85</td><td>0.09</td></tr><tr><td>IRDCs</td><td>32.27</td><td>264</td><td>16.31</td><td>0.19</td><td>1.96</td><td>0.04</td></tr><tr><td>IRDCs24 µm</td><td>32.27</td><td>117</td><td>16.89</td><td>0.27</td><td>2.08</td><td>0.06</td></tr><tr><td>YSOs</td><td>40.05</td><td>57</td><td>17.31</td><td>0.35</td><td>2.14</td><td>0.08</td></tr><tr><td>CH<sub>3</sub>OH maser small</td><td>5–25</td><td>61</td><td>19.23</td><td>1.61</td><td>2.77</td><td>0.1</td></tr><tr><td>CH<sub>3</sub>OH maser large</td><td>40–47</td><td>137</td><td>18.97</td><td>1.65</td><td>2.65</td><td>0.1</td></tr><tr><td>HII regions embedded object</td><td>5</td><td>57</td><td>18.98</td><td>0.37</td><td>2.79</td><td>0.18</td></tr><tr><td>HII regions envelope</td><td>40.3</td><td>194</td><td>18.63</td><td>0.2</td><td>2.63</td><td>0.08</td></tr></tbody></table>

compared to the difference in N<sub>2</sub>H<sup>+</sup> and C<sup>18</sup>O line-centre velocities of low-mass, starless and protostellar cores observed in the Perseus molecular cloud by Kirk et al. (2007). Although the motion of protostellar cores with respect to the envelope could be influenced by turbulence generated by heating or outflows, most of them as well as the majority of the starless cores show relative velocities of both molecules slower than the sound speed, while the remaining sources have only slightly higher differences. This agrees with the analysis of Walsh et al. (2004), who measured a core-to-envelope velocity exceeding the sound speed in only 3% of their sample, which included mostly isolated low-mass cores. Hence, those comparisons also illustrate that the ATLASGAL clumps as a high-mass star forming region sample show much more distinct velocity differences between low- and high-density tracers than low-mass cores.

### 6.3. <sup>13</sup>CO and ammonia linewidths

The width of the <sup>13</sup>CO (1–0) transition is plotted against that of the NH<sub>3</sub> (1, 1) inversion line in the lower panel of Fig. 20, the straight line represents equal widths. For the contour plot in the upper panel of Fig. 20 a binning of 0.4 km s<sup>-1</sup> is chosen for the NH<sub>3</sub> (1, 1) linewidth and of 1 km s<sup>-1</sup> for the <sup>13</sup>CO linewidth. The <sup>13</sup>CO line, exhibiting widths from 1 km s<sup>-1</sup> to 10 km s<sup>-1</sup> with a peak at 3.3 km s<sup>-1</sup>, is mostly broader than the ammonia line, which may indicate that the lines arise from different parts of the clumps. It is expected that <sup>13</sup>CO traces a larger scale than ammonia, although the beamwidth used for the GRS of 46" is only slightly broader than that of the Effelsberg telescope. <sup>13</sup>CO traces densities of ~10<sup>3</sup> cm<sup>-3</sup>, but is frozen out at densities of 10<sup>5</sup> cm<sup>-3</sup>, where NH<sub>3</sub> is still existing in the gas phase. Hence, ammonia probes the massive, denser inner parts of a clump, while <sup>13</sup>CO is also sensitive to the outer less dense regions of the cloud. Turbulence within large structures can cause broader <sup>13</sup>CO than NH<sub>3</sub> linewidths found in small-scale regions of the clumps.

Another explanation of the difference in linewidths could be that multiple <sup>13</sup>CO emitting clumps may lie on the same line of sight to the observer, since <sup>13</sup>CO traces low densities and its abundance of 10<sup>-6</sup> (Du & Yang 2008) is higher compared to the ammonia molecule with 7 × 10<sup>-9</sup> to 1.2 × 10<sup>-6</sup> (cf. Sect. 5.1). Figure 17 also shows that some ATLASGAL source lines of sight have more than one <sup>13</sup>CO component. Those different clumps may also have narrow linewidths, which would add up to the observed increased <sup>13</sup>CO widths due to velocity dispersion within the beam.

## 7. Temperature and linewidth distributions in different evolutionary stages

To investigate trends of NH<sub>3</sub> line parameters with evolutionary phases, we divided most clumps into subsamples by searching

for their associations with other samples that cover the main evolutionary stages of high-mass star formation. We used search radii dependent on the type of sample. Search radius, number of associated sources, mean and error of the NH<sub>3</sub> (1, 1) linewidth, and the rotational temperature of the resulting ATLASGAL subsamples are given in Table 6. While a few sources do not match any sample, some clumps show characteristics of more than one phase of high-mass star formation and are therefore associated with different subsamples. Similar studies using smaller samples have already analysed the association of sources traced by NH<sub>3</sub> with different evolutionary stages of massive star formation (Sridharan et al. 2005; Pillai et al. 2006). Figure 21 shows the number distribution with the (1, 1) linewidth and the rotational temperature of the different ATLASGAL subsamples for clarity as curves and not histograms, although the data sampling is not as good as suggested.

The lowest panel compares the whole ATLASGAL sample as a solid curve in green with a subsample of ATLASGAL clumps not detected at 24 µm, denoted as “24 µm dark” as a dotted curve in black, IRDCs as a dashed red curve from Peretto & Fuller (2009) and a subsample of those IRDCs emitting at 24 µm as a dashed-dotted curve in blue. As search radius for associations of IRDCs with ATLASGAL sources we used the equivalent radius given by Peretto & Fuller (2009), which agrees with the radius of a disk, whose area is the same as that of the IRDC. An average of the search radii used is given in Table 6. The 24 µm dark ATLASGAL clumps are the coldest sources with the narrowest linewidths, suggesting that they are in a very early phase without any hint of star formation. The distribution of the IRDCs and the ATLASGAL sample is similar in linewidth and rotational temperature. Although the IRDCs that emit at 24 µm show signs of star formation activity, they have only slightly broader linewidths than the average of ATLASGAL sources and similar rotational temperatures. The 24 µm emitting star embedded in those IRDCs is a local heating component of little influence over the volume of the whole clump and thus has no significant influence on the source properties. Examples of the NH<sub>3</sub> (1, 1) to (3, 3) lines of IRDCs (G18.09–0.30, G10.67–0.22) are displayed in Fig. 2.

We plot the distribution of young stellar objects (YSOs as dotted black curve) from the RMS survey (Hoare et al. 2005), CH<sub>3</sub>OH masers from Pandian et al. (2011), Szymczak et al. (2002), and Caswell et al. (2010) (dashed and dashed-dotted curves) and the ATLASGAL sample as a solid curve in green in the middle panel. The YSOs are MSX sources (Price et al. 2001), which emit at infrared wavelengths and are therefore in a more evolved phase than the IRDCs. We searched for them within the ATLASGAL beam in addition to the MSX positional uncertainty. The embedded infrared sources have broader (1, 1) linewidths than the mean value of ATLASGAL clumps. This hints at star formation activity in the YSOs, which results
<figure><img src="image_18.png" /></figure>

**Fig. 21.** Relative number distribution of ATLASGAL subsamples in different evolutionary phases of high-mass star formation showing the NH₃ (1, 1) linewidth in the upper part and the rotational temperature in the lower part. The whole ATLASGAL sample is shown as solid green curve in each panel, while the other subsamples are indicated by different patterns and colours: *lower panel:* dotted black curve: ATLASGAL clumps not detected at 24 µm, dashed red: IRDCs from Peretto & Fuller (2009), dashed-dotted blue: an IRDC subsample emitting at 24 µm; *Middle panel:* dotted black curve: young stellar objects (YSOs) from the rms survey (Hoare et al. 2005), dashed red: CH₃OH masers from Pandian et al. (2011), Szymczak et al. (2002), and Caswell et al. (2010) found within a small search radius (see Sect. 7), dashed-dotted blue: same CH₃OH maser samples found in a large radius (see Sect. 7); *upper panel:* UCHII regions from Hoare et al. (2005), dotted black curve: associations with the envelope in a large search radius (see Sect. 7), dashed red curve: associations with embedded object in a small search radius (see Sect. 7).

in higher rotational temperatures than found for the ATLASGAL sample on average (Fig. 21) and at turbulence, e.g. due to molecular outflows. Concerning the methanol masers, we used two search radii for each of the three samples: sources found within a small radius, consisting of the pointing error of the ATLASGAL survey of 5" and those of the telescopes used to observe the CH₃OH samples, might be directly related to the maser (CH₃OH maser small as a dashed red curve). ATLASGAL clumps located in the same complex as the methanol maser source can be revealed within the large radius that includes the Effelsberg beamwidth in addition to the pointing errors (CH₃OH maser large as dashed-dotted blue curve). The range of small and large

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

**Fig. 22.** Cumulative distribution functions illustrating the NH₃ (1, 1) linewidth (*upper panel*) and rotational temperature (*lower panel*) for the three largest subsamples in Table 6. The distribution of IRDCs from Peretto & Fuller (2009) is shown as a solid red line, the CH₃OH sample found within the large search radius is plotted as a dashed blue curve and the H II regions envelope as a dotted black curve.

search radii is given in Table 6. The CH₃OH samples found within both search radii have a similar range of Δv(1, 1) and Trot, although they have broader linewidths and higher Trot than the ATLASGAL sample on average. An ammonia spectrum of a YSO (G27.02+0.20) is shown in Fig. 2.

The upper plot compares the ATLASGAL sample (the solid green curve) with UCHII regions from Hoare et al. (2005). We again used a small and large search radius as described above to distinguish between the embedded object as a dashed red curve, and the envelope of the H II region (the dotted black curve). Because they already contain a protostar heating up the molecular cloud, they are at a more evolved stage of high-mass star formation than the IRDCs and YSOs and are therefore the warmest sources with the broadest linewidths. Since the embedded object traces the central part of the molecular cloud, it covers slightly higher temperatures and broader linewidths than the envelope. In contrast, the envelope includes cold regions as well, which results in a decrease of Trot and Δv(1, 1) by averaging over the large search radius. Figure 2 illustrates an example of NH₃ lines of an embedded UCHIIR (G12.89+0.49).

To investigate if the differences in physical properties given in Table 6 are statistically significant, we performed a Kolmogorov-Smirnov (KS) test using the three largest subsamples from different evolutionary phases. The comparison of the NH₃ (1, 1) linewidth and rotational temperature distributions of the IRDCs from Peretto & Fuller (2009) and the CH₃OH sample found within the large search radius using the KS test showed that both are measurably different (probability >99.9%).
In addition, the KS test used for the (1,1) linewidths and rota-
tional temperatures of the IRDCs with those of the H II regions
envelope from Hoare et al. (2005) reveals that those are signif-
icantly different as well. Cumulative distribution plots (Fig. 22)
show the NH₃ (1,1) linewidth and rotational temperature dis-
tribution for the different subsamples. They illustrate a trend
of broader linewidths and higher temperatures for sources with
CH₃OH masers (dashed curve in blue) and with H II regions
(dotted black line) compared to IRDCs (solid red curve). The
distributions of methanol masers and H II regions are similar.
Our KS test results for the NH₃ (1,1) linewidths and rotational
temperatures do not contradict that the CH₃OH sample and the
H II regions envelope sample are drawn from the same parent
population.

8. Conclusions

The (1,1) to (3,3) ammonia inversion transitions of 898 ATLASGAL sources, located between 5° and 60° in Galactic longitude and ±1.5° in Galactic latitude, were observed with the Effelsberg telescope. They build the largest available ammonia sample of high-mass star forming clumps. Our results are summarised in this section.

We analysed the correlation of different ammonia line pa-
rameters and compared them to the submm dust continuum and
¹³CO emission.

The NH₃ (1,1) line's LSR velocity of most clumps lies be-
tween −10 km s⁻¹ and 120 km s⁻¹, most clumps show a strong
correlation with CO emission (Dame et al. 2001), probing the
larger giant molecular clouds. The velocities yield near/far kine-
matical distances, that were calculated for 750 unknown sources
using the revised rotation parameters of the Milky Way pre-
sented by Reid et al. (2009). Galactocentric distances show an
enhancement of ATLASGAL sources at 4.5 and 6 kpc, tracing
the Scutum arm and Sagittarius arm, which agrees with previous
results for other tracers of massive star formation (Anderson &
Bania 2009).

Most sources exhibit NH₃ (1,1) linewidths between
0.7 km s⁻¹ and 3.5 km s⁻¹, some range up to 7 km s⁻¹, which
is much broader than the thermal NH₃ linewidth and common
linewidths observed in cores of low-mass (Jijina et al. 1999).
The rotational temperature between the (1,1) and (2,2) inversion
levels was derived and was used to calculate the kinetic temper-
ature. Most ATLASGAL clumps have rotational temperatures
between 10 K and 25 K, some exhibit up to 30 K, while their av-
erage kinetic temperatures range from 12 K to 35 K, only few up
to 40 K. The increasing rotational temperature and the broader
widths of the (3,3) line, which probes higher temperature, hint at
star formation activity inside the clumps. Indeed, many clumps
show a stronger than expected intensity of the NH₃ (3,3) line as
estimated for a homogeneous, cold clump with a temperature de-
rived from the (1,1) and (2,2) lines. This results probably from
an embedded hot core component with a low beam filling factor
and high temperature.

We used the kinetic temperature and submm dust continuum
flux to determine the H₂ density. Its correlation with the NH₃ col-
umn density leads to ammonia abundances relative to H₂ mostly
from 5 × 10⁻⁸ to 3 × 10⁻⁷.

For ATLASGAL clumps located at tangential points and as-
sociated with clouds from the GRS survey, preliminary distance
estimates were given and gas masses were calculated that lie in
the range from 60 to 10⁴ M☉. Comparing them with the virial
masses yields a virial parameter of ~1 for sources with broad
linewidths, which are supported against gravitational collapse.

In contrast, clumps with narrow linewidths have a virial param-
eter of only ~0.5.

The comparison of the NH₃ lines as high-density probe with the GRS ¹³CO emission as low-density envelope tracer yields broader ¹³CO than ammonia linewidths. This can result from turbulence within large structures probed by ¹³CO. We also analysed the relative motion of the clump and its cloud using NH₃, which represents the dense core material and ¹³CO to trace the more diffuse molecular cloud gas. Since our comparison yields small differences between NH₃ and ¹³CO velocities, less than the average ¹³CO linewidth, random motions of the ammonia clumps are expected and hint at turbulence. In contrast to low-mass cores, which show mainly subsonic core-to-envelope motions (Walsh et al. 2004, 2007; Kirk et al. 2007), the relative motions of our sample mostly exceed the sound speed.

Recent NH₃ (1,1) to (3,3) line observations of northern BGPS sources were carried out by Dunham et al. (2011) at the Robert F. Byrd Green Bank Telescope (GBT) with a lower FWHM of 31'' compared to the Effelsberg beamwidth. Our results agree well with those of Dunham et al. (2011), yielding similar ranges of kinetic temperature, NH₃ column density, H₂ column density, thus NH₃ abundance, and gas masses.

Association of the ATLASGAL sources with other high-
mass star forming samples led to subsamples in various evo-
lutionary stages, whose NH₃ line parameters were compared.
Clumps with low temperatures, narrow linewidths, and high col-
umn densities may be IRDCs, which are assumed to be the cold
precursors to hot cores and are consequently in an early evo-
lutionary phase of high-mass star formation (Rathborne et al.
2006). Moreover, few ATLASGAL sources show higher temper-
ature, greater linewidths and lower column densities than the
IRDCs. Those properties hint at objects in a more evolved stage
of massive star formation, such as YSOs that contain a luminous
infrared source or clumps with embedded UCHII regions.

Because these NH₃ observations and analyses can also build
a valuable database for the selection of different phases of mas-
sive star formation and the determined velocities are important
especially for the so far unknown ATLASGAL sources, we ex-
tended the NH₃ survey to the south using the Parkes telescope.
The comparison of first and fourth quadrant properties and the
distinction between near and far distances will be published in a
forthcoming paper (Wienen et al., in prep.).

Acknowledgements. We thank the anonymous referee for very useful comments which improved the paper. L.B. acknowledges support from CONICYT projects FONDAP 15010003 and Basal PFB-06. M. Wienen was supported for this research through a stipend from the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne.

**Appendix A: Pointing and calibration**

In this section, we describe the pointing and calibration of the
NH₃ data obtained with the Effelsberg telescope.

Continuum drift scans were taken on nearby pointing sources
such as NGC 7027, 1741-038, G10.62, and NRAO530 every
hour to correct for pointing offsets. NGC 7027 was also used to
calibrate the temperature scale. The focal point was determined
by measuring NGC 7027 at the beginning of each observational
period, before and after sunrise and sunset.

The data taken with the Effelsberg telescope are given in
noise tube units and must be converted into main beam bright-
ness temperature during the calibration procedure. Because
many sources could only be observed at low elevations be-
tween 15° and 25° because of their southern position, several
flux calibrators at different elevations had to be observed. These are the continuum sources, whose flux density was measured during the pointing. In addition, the atmospheric opacity has to be included in the calibration procedure owing to the low elevations of the sources. The primary flux calibration was performed using NGC 7027, since its flux density is precisely known. The other pointing sources can be used as secondary calibrators, which were observed over the entire elevation range and had to be adapted to NGC 7027 as relative calibration.

The sources were calibrated by the total calibration factor $T_{cal,tot}$, which is given by a constant factor $T_{cal}$ belonging to the primary flux calibration divided by a gain elevation curve for each observation day according to

$$T_{cal,tot} = \frac{T_{cal}}{(a + b \cdot El + c \cdot El^2) \cdot \exp(-\tau \cdot amass)} \quad (A.1)$$

The term $y1 = a + b \cdot El + c \cdot El^2$ with the elevation $El$ describes the distortion of the telescope caused by its weight, when it is inclined, as well as thermal conditions. $y2 = e^{-\tau \cdot amass}$ with the zenith opacity $\tau$ and the factor $airmass$, $amass = 1/\sin(El)$, takes the weakening of the radiation penetrating the atmosphere into account.

During the primary flux calibration, the factor $T_{cal}$ was determined using continuum scans of NGC 7027, which possesses a flux density $S_{\lambda} \approx 5.6$ Jy at 1.3 cm (Ott et al. 1994). The main beam brightness temperature $T_{MB}$ of this calibrator for a given beamwidth $\theta$ and wavelength $\lambda$ is given by (Rohlfs & Wilson 2004)

$$(T_{MB}/K) = \frac{(S_{\lambda}/Jy)(\lambda/cm)^2}{(\theta/\text{arcmin})^2 2.65} \quad (A.2)$$

With a beamwidth of $\theta = 40''$ and wavelength of $\lambda = 1.25$ cm, $T_{MB}$ of NGC 7027 is 7.5 K. We used only one calibration factor for the NH$_3$ (1, 1) to (3, 3) inversion transitions, because the variation of the noise diode temperature<fnref n="2" /> is between 4 and 6% and of the flux of NGC 7027 (Ott et al. 1994) is only ~0.1%, both are therefore smaller than the average calibration uncertainty of 10%.

The opacity is given by the water vapour radiometer, although it could not be used in May 2008, because it did not work during this month. Hence, the opacity of the sources observed in that time was obtained by plotting their system temperatures against elevation and fitting them. The system temperature is related to the elevation according to

$$T_{sys} = T_{0} + T_{atm}(1 - \exp(-\tau \cdot amass)), \quad (A.3)$$

with the system temperature $T_{sys}$ for an atmospheric temperature $T_{atm}$, the zenith opacity $\tau$ and the temperature $T_{0}$, which contains the temperature of the receiver and the ground. G15.66–0.50 was observed every day with the Effelsberg telescope. Hence, the $T_{MB}$ variation of its NH$_3$ (1, 1) line could be used to estimate the rms uncertainty of the calibration, which was determined to be ~10%, which is small considering the low elevations of the sources up to 25°. But this does not influence important NH$_3$ line parameters such as the optical depth or the rotational and kinetic temperature of the clumps, because their determination does not require calibration. The optical depth was obtained by the ratio of the hyperfine components of one inversion transition. The rotational and thus kinetic temperature were calculated by the intensity ratio of the ammonia (1, 1) and (2, 2) lines, which can be observed simultaneously. However, the data needed to be calibrated to estimate e.g. the source-averaged column density.

## References

Anderson, L. D., & Bania, T. M. 2009, ApJ, 690, 706

André, P., Ward-Thompson, D., & Barsony, M. 2000, Protostars and Planets IV, 59

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

Becker, R. H., White, R. L., Helfand, D. J., & Zoonematkermani, S. 1994, ApJS, 91, 347

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

Bergin, E. A., & Langer, W. D. 1997, ApJ, 486, 316

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

Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945

Beuther, H., Churchwell, E. B., McKee, C. F., & Tan, J. C. 2007, Protostars and Planets V, 165

Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 1997, MNRAS, 285, 201

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

Braz, M. A., & Epchtein, N. 1983, A&AS, 54, 167

Bronfman, L., Casassus, S., May, J., & Nyman, L.-Å. 2000, A&A, 358, 521

Carey, S. J., Clark, F. O., Egan, M. P., et al. 1998, ApJ, 508, 721

Carey, S. J., Feldman, P. A., Redman, R. O., et al. 2000, ApJ, 543, L157

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

Cesaroni, R., Walmsley, C. M., & Churchwell, E. 1992, A&A, 256, 618

Cheung, A. C., Rank, D. M., Townes, C. H., Knowles, S. H., & Sullivan, III, W. T. 1969, ApJ, 157, L13

Churchwell, E., Walmsley, C. M., & Cesaroni, R. 1990, A&AS, 83, 119

Csengeri, T., Bontemps, S., Schneider, N., & Motte, F. 2011, in IAU Symp. 270, eds. J. Alves, B. G. Elmegreen, J. M. Girart, & V. Trimble, 53

Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792

Devine, K. E., Chandler, C. J., Brogan, C., et al. 2011, ApJ, 733, 44

Du, F., & Yang, J. 2008, ApJ, 686, 384

Dunham, M. K., Rosolowsky, E., Evans, II, N. J., et al. 2010, ApJ, 717, 1157

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

Egan, M. P., Shipman, R. F., Price, S. D., et al. 1998, ApJ, 494, L199

Hill, T., Longmore, S. N., Pinte, C., et al. 2010, MNRAS, 402, 2682

Ho, P. T. P., & Townes, C. H. 1983, ARA&A, 21, 239

Hoare, M. G., Lumsden, S. L., Oudmaijer, R. D., et al. 2005, in Massive Star Birth: A Crossroads of Astrophysics, eds. R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley, IAU Symp., 227, 370

Hofner, P., Kurtz, S., Churchwell, E., Walmsley, C. M., & Cesaroni, R. 1994, ApJ, 429, L85

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

Jijina, J., Myers, P. C., & Adams, F. C. 1999, ApJS, 125, 161

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

Kirk, H., Johnstone, D., & Tafalla, M. 2007, ApJ, 668, 1042

Kroupa, P., Weidner, C., Pflamm-Altenburg, J., et al. 2011 [arXiv:1112.3340]

Lekht, E. E. 2000, A&AS, 141, 185

Mac Low, M.-M., & Klessen, R. S. 2004, Rev. Mod. Phys., 76, 125

Mangum, J. G., & Wootten, A. 1994, ApJ, 428, L33

Martin, P. G., Roy, A., Bontemps, S., et al. 2012, ApJ, 751, 28

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

Mckee, C. F., Zweibel, E. G., Goodman, A. A., & Heiles, C. 1993, in Protostars and Planets III, eds. E. H. Levy, & J. I. Lunine, 327

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

Molinari, S., Testi, L., Rodríguez, L. F., & Zhang, Q. 2002, ApJ, 570, 758

Molinari, S., Swinyard, B., Bally, J., et al. 2010, PASP, 122, 314

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

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

Nguyen Luong, Q., Motte, F., Hennemann, M., et al. 2011a, A&A, 535, A76

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

Olmi, L., Araya, E. D., Chapin, E. L., et al. 2010, ApJ, 715, 1132

Ormel, C. W., Min, M., Tielens, A. G. G. M., Dominik, C., & Paszun, D. 2011, A&A, 532, A43

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

Ott, M., Witzel, A., Quirrenbach, A., et al. 1994, A&A, 284, 331

Pandian, J. D., Momjian, E., Xu, Y., Menten, K. M., & Goldsmith, P. F. 2011, ApJ, 730, 55

Perault, M., Omont, A., Simon, G., et al. 1996, A&A, 315, L165

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

Peretto, N., & Fuller, G. A. 2010, ApJ, 723, 555

Pillai, T., Wyrowski, F., Carey, S. J., & Menten, K. M. 2006, A&A, 450, 569

<footnote n="2">² See http://www.mpifr-bonn.mpg.de/div/effelsberg/calibration/1.3cmpf.html</footnote>
Price, S. D., Egan, M. P., Carey, S. J., Mizuno, D. R., & Kuchar, T. A. 2001, AJ, 121, 2819

Ragan, S. E., Bergin, E. A., & Wilner, D. 2011, ApJ, 736, 163

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

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

Reifenstein, E. C., Wilson, T. L., Burke, B. F., Mezger, P. G., & Altenhoff, W. J. 1970, A&A, 4, 357

Rohlfs, K., & Wilson, T. L. 2004, Tools of radio astronomy, 4th rev. and enl., eds. K. Rohlfs, & T. L. Wilson (Berlin: Springer)

Roman-Duval, J., Jackson, J. M., Heyer, M., et al. 2009, ApJ, 699, 1153

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

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data Analysis Software and Systems IV, eds. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, ASP Conf. Ser., 77, 433

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

Schulz, N. S. 2005, From Dust To Stars Studies of the Formation and Early Evolution of Stars

Simon, R., Jackson, J. M., Rathborne, J. M., & Chambers, E. T. 2006a, ApJ, 639, 227

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

Siringo, G., Weiss, A., Kreysa, E., et al. 2007, The Messenger, 129, 2

Siringo, G., Kreysa, E., Kovacs, A., et al. 2008, in SPIE Conf. Ser., 7020

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

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

Sridharan, T. K., Beuther, H., Saito, M., Wyrowski, F., & Schilke, P. 2005, ApJ, 634, L57

Szymczak, M., Kus, A. J., Hrynek, G., Kěpa, A., & Pazderski, E. 2002, A&A, 392, 277

Szymczak, M., Pillai, T., & Menten, K. M. 2005, A&A, 434, 613

Tafalla, M., Myers, P. C., Caselli, P., Walmsley, C. M., & Comito, C. 2002, ApJ, 569, 815

Tafalla, M., Myers, P. C., Caselli, P., & Walmsley, C. M. 2004, A&A, 416, 191

Takano, S., Nakai, N., & Kawaguchi, K. 2002, PASJ, 54, 195

Ungerechts, H., Winnewisser, G., & Walmsley, C. M. 1986, A&A, 157, 207

Ungerechts, H., Bergin, E. A., Goldsmith, P. F., et al. 1997, ApJ, 482, 245

van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627

Walmsley, C. M., & Ungerechts, H. 1983, A&A, 122, 164

Walsh, A. J., Myers, P. C., & Burton, M. G. 2004, ApJ, 614, 194

Walsh, A. J., Myers, P. C., Di Francesco, J., et al. 2007, ApJ, 655, 958

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

Wood, D. O. S., & Churchwell, E. 1989b, ApJS, 69, 831

Zoonematkermani, S., Helfand, D. J., Becker, R. H., White, R. L., & Perley, R. A. 1990, ApJS, 74, 181