Ammonia from cold high-mass clumps discovered in the inner Galactic disk by the ATLASGAL survey
Key Words.:Surveys — Submillimeter — Radio lines: ISM — ISM: molecules — ISM: kinematics and dynamics — Stars: formation
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 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 (,) = (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) with 1.9 km s as median, which by far exceeds the thermal linewidth, as well as a broad distribution of high column densities from to cm (median of cm) with an NH abundance in the range of 5 to . 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.
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 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 catalogues 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 ( cm) and have high column densities ( 10-10 cm, Rathborne et al. 2006; Simon et al. 2006a). Their sizes ( 2 pc) and masses ( 10-10 M) are similar to those of cluster-forming molecular clumps. This, together with the cold temperatures and fragmented substructure 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 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 UCHIIRs. In contrast, a first attempt to obtain an unbiased sample 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 objects 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 submillimeter wavelengths. ATLASGAL (The APEX Telescope Large Area Survey of the Galaxy at 870 m) 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 and latitude of , 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 continuum 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 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 cm Beuther et al. 2002) and cold, many molecules such as CS and CO are partly frozen onto dust grains. In contrast to those, the fractional gas abundance of ammonia remains constant in prestellar cores (Tafalla et al. 2002). Hence, we carried out follow-up observations of northern ATLASGAL sources in the lowest NH inversion transitions.
NH 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 momentum, , and its projection along the molecular axis, . Radiative transitions between different -ladders are forbidden and the lowest metastable energy levels, for which , are thus collisionally 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 temperature 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 reliable 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 lines of sources in the first Galactic quadrant as well, identified 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 compared to our NH survey. Our analysis is thus complementary to the results of Dunham et al. (2011).
In Section 2, we describe which ATLASGAL sources were selected for our NH follow-up observations and provide details of the measurements. In Section 3, we present how we reduced the data and derived different clump properties from fits to the spectra of NH inversion transitions. In Section 4, different correlations of the NH line parameters such as the velocity distribution, linewidth, rotational temperature and column density are shown. Then, we compare these with submillimeter dust continuum properties of the clumps such as the H column density in Section 5 and estimate gas and virial masses. We analyse additional line parameters derived from CO (1-0) emission (Jackson et al. 2006) in Section 6 to investigate clump-to-cloud motions etc. In Section 7, the determined gas properties are compared with those of other high-mass dust-selected star forming regions. A summary and conclusions of this NH investigation are given in Section 8.
This paper analyses statistics of the NH 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.
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 to 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, while the FFTS consists of two modules with a chosen bandwidth of 500 MHz each and a spectral resolution of km s. 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. 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 minutes. For some clumps at low declination, which could only be observed at low elevations and consequently have high system temperatures of K, longer integration times of up to 25 minutes were chosen.
3 Data reduction and analysis
To reduce the NH spectra we used the CLASS software111available at http://www.iram.fr/IRAMFR/GILDAS. 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, while a low order, mostly 3, was chosen for lines with a broad width, greater than 4 km s.
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 the model and the fitted observation agreed, the deviation of the two becomes sligthly 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, , the linewidth, , at the full width at half maximum of a Gaussian profile and the optical depth of the main line, , 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, , the radial velocity of the main line, , and the FWHM linewidth, . 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 temperature were calculated with
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, (1,1), linewidths, (1,1), and main beam brightness temperatures, (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 (), the kinetic temperature () and the ammonia column density () 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.
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 signal-to-noise (S/N) ratio 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.620.08 with an S/N ratio of 30. An example for a weak line is the (3,3) line in G10.990.08. Its signal-to-noise ratio 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.360.15 and G23.270.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.360.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 = ) (Fig. 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). Fig. 17 compares NH (1,1) and CO (1-0) lines. Fig. 18 shows our analysis of CO and NH line-centre velocities, both molecular lines for 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 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 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. Fig. 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 complimentary 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).
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 (1,1) = 6.4 km s and (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 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 , 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
with the kinetic temperature , the Boltzmann constant and the mass of the ammonia molecule . Assuming a mean of 20 K, the thermal linewidth is about 0.22 km s. On average measured NH(1,1) linewidths are 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 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)
with the optical depth of the (1,1) main line, , and the main beam brightness temperatures, , 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 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 %) 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.370.66 with (1,1) = 0.78 km s and = 10 K, a known IRDC (Peretto & Fuller 2009). Among the clumps with broad linewidth are well-known sources, such as the UCHIIR/hot core G10.47+0.03 (Cesaroni et al. 1992) (see Sect. 4.2) with (1,1) = 6.38 km s and = 21 K. Another source with a high of 26 K and (1,1) of 5.06 km s is G43.800.13, well-known for harbouring OH (Braz & Epchtein 1983) and HO maser emission (Lekht 2000). G12.210.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 (1,1) of 5.73 km s as well and 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 equation 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 () assuming that the energy levels are populated via a Boltzmann 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 Rohlfs & Wilson (2004)
where and are the statistical weights of the lower and upper inversion level, the Einstein coefficient, the inversion frequency in GHz and 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):
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 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 cm and cm with an average error of 17% (cf. Table 3) and a peak at cm. The source with the highest column density of 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 high-mass protostellar objects (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 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
where is the the main beam brightness temperature of the (1,1) inversion transition, the optical depth of the (1,1) main line and =2.73 K. Because the H density of the ATLASGAL clumps is 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 (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 10cm, 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 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 (3,3).
Sources towards which the NH (3,3) line is not detected have an average (1,1) linewidth of 1.8 km s, which is slightly below the mean value of the whole ATLASGAL sample of 2 km s. Their average rotational temperature is 15.3 K, which is also lower than the mean of 16.5 K.
Assuming LTE (and thus =), the (3,3) main beam brightness temperature is
assuming the same beam filling factor as for the (1,1) line (see Eq. 7) and the optical depth of the NH(3,3) line given by
where is the frequency of the (3,3) inversion transition (23.87 GHz), the Boltzmann constant and the measured kinetic temperatures. is defined as
Here, denotes the speed of light, (1,1) the observed (1,1) linewidth, the (3,3) level population number and the Einstein 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 (cf. Table 3) (Rohlfs & Wilson 2004)
where Q is given by
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 (3,3)/(1,1) is 1.6 with a dispersion of 0.8 and the median of / 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:
Non-thermal excitation can overpopulate or invert the upper level of the (3,3) doublet. Different examples of these maser lines have already been observed: Mangum & Wootten (1994) detected NH (3,3) maser emission towards the 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.820.18 shows two narrow components at lower and higher velocities than in addition to the thermal emission. This clump also harbours CHOH and water masers (Szymczak et al. 2005). In the NH (3,3) spectrum of G30.720.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, HO or CHOH 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 values compared to calculations.
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 ( 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.
High values might originate from a hot embedded source, which is revealed by broad (3,3) linewidths leading to a high ratio of (3,3) to (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 / 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 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) 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 GHz. The integrated dust continuum fluxes S 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 (), which is related to the ammonia column density () and abundance () via the relation = /. 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 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, , of 1.85 cm/g at 870 m at a gas density n(H) = 10 cm (Ossenkopf & Henning 1994). We estimated the H column density per the relation (Kauffmann et al. 2008)
where is the 870 m flux density in Jy integrated within 40, the wavelength in m and the dust temperature with = under the assumption of equal gas and dust temperatures. is the 100-m telescope beamwidth of 40 and Z/Z the ratio of metallicity to the solar metallicity, assuming Z/Z = 1. We found that most H column densities lie between and cm with a peak at cm. For each clump we can now use the ammonia column density to derive the NH abundance via =/, which lies in the range from to with an average of . The H column densities, NH abundances and the integrated dust continuum fluxes 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 / to 0.0271 for /. Deviations from the trend in Fig. 15 are due to a higher ammonia abundance, clumps with are below the lower straight line, or a decreased abundance, values lower than 5 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 to . Dunham et al. (2011) found NH abundances between and , which is consistent with our results. Moreover, Pillai et al. (2006) obtained NH abundances from 10 to 10 for a sample of IRDCs, their 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)
with the distance to the source, , in pc; for the other parameters see equation 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 (Rohlfs & Wilson 2004)
with the linewidth , corrected for the resolution of the spectrometer . The radius in pc is obtained from Gaussian fits of the 870 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 M and virial masses from 20 to 5.5 M. 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 ==2. 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)
which is 1 for clumps supported against gravitational collapse. When we use = from equation 13, sources with narrow linewidths have virial parameters on average much smaller (mean ) than broad linewidth sources with a mean, , of 0.45. With higher density tracers such as HCO and CS 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 1 for clumps with broad linewidth. Therefore, these sources are in virial equilibrium. In contrast, sources with small (1,1) still have a low virial parameter of 0.42.
Hill et al. (2010), who observed NH 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 NH (1,1) linewidth of 2.9 km s in contrast to our mean (1,1) of 2 km s, which results in larger virial masses than we found.
Dunham et al. (2010) obtained an average virial parameter of 1 from their NH 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 and therefore probed much smaller masses compared to the ATLASGAL clumps, whose gas masses range up to 1.6 M.
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 (McKee et al. 1993), where is the hydrogen column density in units of 10 cm. For our values of the hydrogen column density between and cm, we obtain magnetic fields between 5 and 450 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 () 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. (2012b, 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 NH and CO (1-0) emission
While ammonia is a high-density tracer ( 10 cm; Ungerechts et al. 1986), the CO molecule can be used to probe low-density regions ( cm; Ungerechts et al. 1997), which are surrounding the dense cores. In this section, the ammonia observations are compared with the CO J = 1 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 km s, comparable or better angular resolution (46) and sampling (22) than previous molecular line surveys of the inner Galaxy. CO is used instead of the commonly observed CO J = 1 0 emission, because 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 is covered by the whole survey with Galactic longitudes of l = 18 - 55.7, limited data are available for 14, and latitudes of . The velocity coverage ranges from -5 to 135 km s up to a Galactic longitude of 40 and from -5 to 85 km s for l . We extracted 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 CO LSR velocity of sources detected by GRS and ATLASGAL, their linewidths and integrated temperatures T, each with errors. The results are given in Table 5, NH and CO line parameters are compared in correlation plots.
6.1 Comparison of ammonia and CO lines
We extracted the CO (1-0) line from GRS intensity images at the position of the ammonia observations, fitted a Gaussian to the CO line profile and related the CO line with the velocity closest to the NH velocity to the observed ATLASGAL source. Within the ammonia sample we found 517 ATLASGAL clumps associated with CO emission. The spectra of CO and NH (1,1) lines are overlaid, some examples of strong sources are shown in Fig. 17. The CO linewidths are similar to the widths of the NH (1,1) inversion transitions for a few sources such as G16.910.08, G21.390.25, G28.82+0.36, G29.600.62, while most clumps possess broader CO than NH (1,1) linewidths, e.g. G15.43+0.19, G18.820.47, G18.880.51, G22.37+0.38, G18.830.48, and G28.920.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.910.08, G21.390.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 NH as high-density tracer and CO as well as CO 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 NGC1333 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 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.680.03 with a deviation in both velocities of -2.4 km s has a lower CO velocity, while G37.760.22 exhibits a lower NH velocity with a difference to CO of 2.24 km s. Fig. 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 NH and CO or CO in cores observed by Walsh et al. (2004, 2007) are also similar to their NH 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 NH and CO linewidth of km s and 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 NH and CO velocities similar to the NH 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 compared to the difference in NH and CO 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 CO and ammonia linewidths
The width of the CO (1-0) transition is plotted against that of the NH (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 is chosen for the NH (1,1) linewidth and of 1 km/s for the CO linewidth. The CO line, exhibiting widths from 1 km s to 10 km s with a peak at 3.3 km s, is mostly broader than the ammonia line, which may indicate that the lines arise from different parts of the clumps. It is expected that 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. CO traces densities of cm, but is frozen out at densities of 10 cm, where NH is still existing in the gas phase. Hence, ammonia probes the massive, denser inner parts of a clump, while CO is also sensitive to the outer less dense regions of the cloud. Turbulence within large structures can cause broader CO than NH linewidths found in small-scale regions of the clumps.
Another explanation of the difference in linewidths could be that multiple CO emitting clumps may lie on the same line of sight to the observer, since CO traces low densities and its abundance of 10 (Du & Yang 2008) is higher compared to the ammonia molecule with to (cf. Sect. 5.1). Fig. 17 also shows that some ATLASGAL source lines of sight have more than one CO component. Those different clumps may also have narrow linewidths, which would add up to the observed increased CO widths due to velocity dispersion within the beam.
7 Temperature and linewidth distributions in different evolutionary stages
To investigate trends of NH 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 (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 with different evolutionary stages of massive star formation (Sridharan et al. 2005; Pillai et al. 2006). Fig. 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 sligthly 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 (1,1) to (3,3) lines of IRDCs (G18.090.30, G10.670.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), CHOH 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 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 CHOH samples, might be directly related to the maser (CHOH 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 (CHOH maser large as dashed-dotted blue curve). The range of small and large search radii is given in Table 6. The CHOH samples found within both search radii have a similar range of (1,1) and , although they have broader linewidths and higher 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 and (1,1) by averaging over the large search radius. Fig. 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 CHOH 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 rotational temperatures of the IRDCs with those of the H II regions envelope from Hoare et al. (2005) reveals that those are significantly different as well. Cumulative distribution plots (Fig. 22) show the NH (1,1) linewidth and rotational temperature distribution for the different subsamples. They illustrate a trend of broader linewidths and higher temperatures for sources with CHOH 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 CHOH sample and the H II regions envelope sample are drawn from the same parent population.
The (1,1) to (3,3) ammonia inversion transitions of 898 ATLASGAL sources, located between 5 and in Galactic longitude and 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 parameters and compared them to the submm dust continuum and CO emission.
The NH (1,1) line’s LSR velocity of most clumps lies between -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 kinematical distances, that were calculated for 750 unknown sources using the revised rotation parameters of the Milky Way presented 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 temperature. Most ATLASGAL clumps have rotational temperatures between 10 K and 25 K, some exhibit up to 30 K, while their average 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 derived 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 column density leads to ammonia abundances relative to H mostly from 510 to 310.
For ATLASGAL clumps located at tangential points and associated 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 parameter 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 evolutionary stages, whose NH line parameters were compared. Clumps with low temperatures, narrow linewidths, and high column densities may be IRDCs, which are assumed to be the cold precursors to hot cores and are consequently in an early evolutionary phase of high-mass star formation (Rathborne et al. 2006). Moreover, few ATLASGAL sources show higher temperatures, 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 massive star formation and the determined velocities are important especially for the so far unknown ATLASGAL sources, we extended 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. LB 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 NGC7027, 1741-038, G10.62, and NRAO530 every hour to correct for pointing offsets. NGC7027 was also used to calibrate the temperature scale. The focal point was determined by measuring NGC7027 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 brightness temperature during the calibration procedure. Because many sources could only be observed at low elevations between 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 NGC7027, 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 NGC7027 as relative calibration.
The sources were calibrated by the total calibration factor , which is given by a constant factor belonging to the primary flux calibration divided by a gain elevation curve for each observation day according to
The term y1 with the elevation describes the distortion of the telescope caused by its weight, when it is inclined, as well as thermal conditions. y2 with the zenith opacity and the factor airmass, , takes the weakening of the radiation penetrating the atmosphere into account.
During the primary flux calibration, the factor was determined using continuum scans of NGC7027, which possesses a flux density Jy at 1.3 cm (Ott et al. 1994). The main beam brightness temperature of this calibrator for a given beamwidth and wavelength is given by (Rohlfs & Wilson 2004)
With a beamwidth of and wavelength of cm, of NGC7027 is 7.5 K. We used only one calibration factor for the NH (1,1) to (3,3) inversion transitions, because the variation of the noise diode temperature222see http://www.mpifr-bonn.mpg.de/div/effelsberg/calibration/1.3cmpf.html is between 4 and 6% and of the flux of NGC7027 (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
with the system temperature for an atmospheric temperature , the zenith opacity and the temperature , which contains the temperature of the receiver and the ground. G15.660.50 was observed every day with the Effelsberg telescope. Hence, the variation of its NH(1,1) line could be used to estimate the rms uncertainty of the calibration, which was determined to be , which is small considering the low elevations of the sources up to 25. But this does not influence important NH 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.
|RAUnits of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds.||DecUnits of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds.||(1,1)||(1,1)||(1,1)||T(1,1)|
|Name||(J2000)||(J2000)||(km s)||(km s)||(K)|
|G5.39+0.19||17 57 12.9150||-24 12 28.150||1.71 0.29)||10.58 0.03)||1.42 0.06)||1.7 0.11)|
|G5.35+0.10||17 57 29.2160||-24 16 57.670||2.11 0.31)||11.21 0.03)||1.48 0.07)||1.73 0.13)|
|G5.64+0.24||17 57 34.6840||-23 58 01.950||2.2 0.13)||7.94 0.01)||1.61 0.03)||4.22 0.12)|
|G5.62-0.08||17 58 44.6080||-24 08 39.020||3.14 0.14)||-25.74 0.02)||2.43 0.03)||3.24 0.11)|
|G6.22-0.05||17 59 55.2530||-23 36 10.150||3.13 0.63)||19 0.09)||2.78 0.15)||0.62 0.09)|
|G6.19-0.12||18 00 07.4100||-23 39 51.830||1.5 0.41)||10.88 0.06)||1.99 0.14)||0.99 0.12)|
|G5.89-0.29||18 00 07.5280||-24 00 29.480||2.09 0.23)||10.11 0.02)||1.42 0.04)||1.88 0.1)|
|G5.89-0.32||18 00 14.4580||-24 01 23.810||0.9 0.13)||10.11 0.02)||1.74 0.05)||2.6 0.09)|
|G6.25-0.12||18 00 16.1890||-23 37 09.950||1.71 0.25)||14.48 0.06)||3.31 0.11)||1.02 0.08)|
|G5.83-0.40||18 00 24.8380||-24 06 50.090||0.91 0.21)||7.64 0.03)||1.67 0.07)||1.54 0.08)|
|G5.89-0.39||18 00 30.3380||-24 04 00.209||0.69 0.09)||9.18 0.03)||3.96 0.06)||4.05 0.13)|
|G5.90-0.43||18 00 40.3520||-24 04 18.700||1.8 0.06)||6.64 0.01)||2.53 0.03)||3.8 0.09)|
|G5.83-0.51||18 00 51.2040||-24 10 23.290||1.99 0.12)||16.22 0.01)||1.74 0.03)||2.8 0.08)|
|G6.19-0.36||18 01 01.7140||-23 47 05.130||2.11 0.12)||-33.81 0.02)||2.2 0.04)||3.72 0.12)|
|G5.91-0.54||18 01 08.2520||-24 07 14.800||2.09 0.22)||14.79 0.02)||1.37 0.05)||2.32 0.11)|
|Name||(km s)||(km s)||(K)||(km s)||(km s)||(K)|
|G5.39+0.19||10.6 0.13)||1.53 0.38)||0.47 0.12)||-||-||0.09|
|G5.35+0.10||10.85 0.17)||1.8 0.35)||0.38 0.11)||-||-||0.08|
|G5.64+0.24||7.95 0.02)||1.89 0.06)||2.88 0.12)||8.29 0.09)||2.81 0.2)||0.82 0.09)|
|G5.62-0.08||-25.57 0.05)||2.93 0.11)||1.93 0.14)||-25.47 0.05)||3.54 0.13)||1.85 0.11)|
|G6.19-0.12||10.63 0.26)||2.91 0.62)||0.3 0.1)||-||-||0.1|
|G5.89-0.29||10.13 0.06)||1.65 0.16)||0.85 0.1)||-||-||0.1|
|G5.89-0.32||10.09 0.06)||1.87 0.14)||1.21 0.13)||10.04 0.06)||1.51 0.23)||1.01 0.1)|
|G6.25-0.12||13.97 0.2)||3.02 0.42)||0.44 0.12)||14.35 0.39)||4.06 1.18)||0.31 0.1)|
|G5.83-0.40||7.77 0.11)||2.48 0.3)||0.67 0.11)||-||-||0.09|
|G5.89-0.39||9.28 0.04)||4.58 0.09)||3.37 0.15)||-||-||-|
|G5.90-0.43||6.58 0.04)||3.29 0.08)||2.92 0.14)||6.6 0.05)||3.87 0.12)||1.76 0.09)|
|G5.83-0.51||16.16 0.06)||2.35 0.14)||1.13 0.1)||-||-||0.09|
|G6.19-0.36||-33.84 0.04)||2.83 0.1)||2.71 0.14)||-33.98 0.04)||3.41 0.1)||1.82 0.09)|
|G5.91-0.54||14.9 0.05)||1.59 0.13)||1.32 0.13)||14.43 0.33)||5.19 0.79)||0.29 0.1)|
|G5.39+0.19||13.12 1.27)||14.39 1.67)||15.12 0.08)|
|G5.35+0.10||11.66 1.15)||12.52 1.43)||15.23 0.07)|
|G5.64+0.24||19.53 0.82)||23.94 1.39)||15.34 0.03)|
|G5.62-0.08||16.08 0.81)||18.5 1.19)||15.63 0.02)|
|G6.19-0.12||13.92 1.99)||15.45 2.7)||15.21 0.12)|
|G5.89-0.29||15.45 1.02)||17.58 1.47)||15.21 0.05)|
|G5.89-0.32||18.16 1.08)||21.67 1.73)||14.97 0.07)|
|G6.25-0.12||15.79 2.13)||18.07 3.09)||15.5 0.07)|
|G5.83-0.40||17.51 1.53)||20.66 2.38)||14.95 0.1)|
|G5.89-0.39||27.33 1.33)||39.65 3.21)||15.36 0.06)|
|G5.90-0.43||22.69 1.15)||29.63 2.22)||15.5 0.03)|
|G5.83-0.51||14.82 0.64)||16.68 0.89)||15.28 0.03)|
|G6.19-0.36||20.85 1.17)||26.22 2.08)||15.48 0.03)|
|G5.91-0.54||17.44 1.22)||20.55 1.89)||15.22 0.05)|
|N||S||Distance (near)||Distance (far)|
|Name||(km s)||(km s)||(K)|
|G14.62+0.33||26.73 0.03)||2.56 0.09)||3.26 0.17)|
|G14.63+0.31||26.03 0.04)||2.34 0.1)||3.65 0.16)|
|G14.92+0.07||25.93 0.05)||3.35 0.12)||3.76 0.26)|
|G15.01+0.01||25.38 0.13)||4.35 0.42)||2.65 0.15)|
|G15.43+0.19||47.85 0.03)||3.13 0.06)||3.31 0.17)|
|G14.99-0.12||49.46 0.06)||3.15 0.19)||1.67 0.18)|
|G15.46-0.00||66.92 0.06)||3.3 0.13)||1.66 0.12)|
|G14.89-0.40||61.64 0.13)||2.65 0.32)||2.74 0.66)|
|G15.20-0.44||20.86 0.03)||2.38 0.08)||3.69 0.13)|
|G15.53-0.41||39.27 0.07)||4.36 0.18)||1.97 0.19)|
|G15.56-0.46||15.24 0.03)||1.78 0.06)||3.92 0.28)|
|G16.36-0.07||47.25 0.05)||2.72 0.13)||5.75 0.19)|
|G16.58-0.05||59.41 0.02)||3.27 0.06)||6.53 0.17)|
|G16.83+0.08||62.43 0.03)||2.48 0.07)||2.35 0.15)|
|G16.58-0.08||40.85 0.17)||3.74 0.57)||2.32 0.13)|
|Sample||Search radius ()||Number||T (K)||(T)||(1,1) (km s)||((1,1))|
|CHOH maser small||5 - 25||61||19.23||1.61||2.77||0.1|
|CHOH maser large||40 - 47||137||18.97||1.65||2.65||0.1|
|HII regions embedded object||5||57||18.98||0.37||2.79||0.18|
|HII regions envelope||40.3||194||18.63||0.2||2.63||0.08|
- Anderson & Bania (2009) Anderson, L. D. & Bania, T. M. 2009, ApJ, 690, 706
- André et al. (2000) André, P., Ward-Thompson, D., & Barsony, M. 2000, Protostars and Planets IV, 59
- Bally et al. (2010) Bally, J., Anderson, L. D., Battersby, C., et al. 2010, A&A, 518, L90+
- Becker et al. (1994) Becker, R. H., White, R. L., Helfand, D. J., & Zoonematkermani, S. 1994, ApJS, 91, 347
- Benjamin et al. (2003) Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953
- Bergin & Langer (1997) Bergin, E. A. & Langer, W. D. 1997, ApJ, 486, 316
- Bertoldi & McKee (1992) Bertoldi, F. & McKee, C. F. 1992, ApJ, 395, 140
- Beuther et al. (2007) Beuther, H., Churchwell, E. B., McKee, C. F., & Tan, J. C. 2007, Protostars and Planets V, 165
- Beuther et al. (2002) Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945
- Bonnell et al. (1997) Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 1997, MNRAS, 285, 201
- Bonnell et al. (2001) Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 2001, MNRAS, 323, 785
- Bontemps et al. (2010) Bontemps, S., Motte, F., Csengeri, T., & Schneider, N. 2010, A&A, 524, A18+
- Braz & Epchtein (1983) Braz, M. A. & Epchtein, N. 1983, A&AS, 54, 167
- Bronfman et al. (2000) Bronfman, L., Casassus, S.,