Free Access
Issue
A&A
Volume 567, July 2014
Article Number A40
Number of page(s) 31
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/201219692
Published online 07 July 2014

© ESO, 2014

1. Introduction

The high-mass (O and early B-type) stars play a critical role in the mass distribution and chemical evolution of the galaxies. Yet, their forming processes are still not fully understood. The main difficulty in revealing their birth is that, most young high-mass stars are deeply obscured by a dust envelope and therefore, hard to observe at optical wavelengths. Moreover, they are usually at large distances of several kilo parsecs and have a clustered environment.

Masers, such as CH3OH masers (Reid & Menten 1991; Walsh et al. 1998; Beuther et al. 2002), were found to be an ideal tracer for the forming sites of young high-mass stars. In the targeted survey at early decades, the methanol masers were observed exclusively near the Ultra Compact (UC) Hii regions (Menten et al. 1992; Sutton et al. 2001) and totally absent in other types of objects, including the forming regions of low-mass stars, AGB star, and other galaxies. The subsequent observations (Phillips et al. 1998; Walsh et al. 1998) have revealed a greater number of 6.7 GHz masers without a radio continuum counterpart and suggest that they could trace even earlier stages (prior to the UC Hii stage). Based on the accumulated data and newly performed surveys, Pestalozzi et al. (2005) have presented a relatively completed catalogue of the 6.7 GHz methanol masers throughout the Galaxy. After that, more extensive molecular molecular line surveys (Minier et al. 2005; Purcell et al. 2006; Longmore et al. 2006, etc.) were followed as a major procedure to investigate the physical and chemical properties in the methanol-maser sources. As a result, they are detected in many dense molecular tracers such as H13CO+, CH3CN, and sulfur-bearing molecules.

On the other hand, the different methanol maser transitions (classified to be class-I and class-II, Reid & Menten 1991) are found to trace different physical conditions. Besides the widely detected 6.7 GHz masers, which belong to class-II (Reid & Menten 1991), various class-I masers were observed at the interface between the outflow and the ambient gas (Plambeck & Menten 1990; Johnston et al. 1997; Voronkov et al. 2005, etc.). The latter observations with higher resolution showed that the Class-I masers arise from intense outflow and shock emissions (Cyganowski et al. 2008; Wu et al. 2010). These properties are consistent with the theoretical anticipation that the Class-I masers are stimulated by the collisional processes while the class-II transitions are mainly excited by the radiative processes (Cragg et al. 1992).

Using the 13.7-m millimeter telescope at the Qinghai station of the Purple Mountain Observatory, Liu et al. (2010b) observed 98 maser sources (including both type I and II) and studied the relation between their CO-line emissions and infrared luminosities. Since the 13.7-m telescope is appropriate for single-point and mapping observations over a large area, we carried out a more extensive observation towards the Galactic methanol-maser sources, aiming to get a more complete sample of the molecular lines and thereby analyse their physical parameters and statistical properties.

The targeted 6.7 GHz maser sources (referred to as “maser sources” hereafter) for the current observation are taken from Pestalozzi et al. (2005)’s catalogue and mainly in the northern sky (within the Dec. range available for the telescope) in 12CO and 13CO (1 − 0). We estimated the physical parameters from these lines, examined their Galactic distribution and correlations. In addition, we made 13CO\hbox{${\rm ^{13}CO}$} maps for 33 sources in order to study the gas distribution around the maser sources. We report the observation and data reductions in Sect. 2. We show the measurement of the physical parameters in Sect. 3. In Sect. 4, we discuss the statistical properties of the physical parameters and made a comparison with the physical parameters in the infrared dark clouds (IRDCs) that are sampled by Du & Yang (2008). A summary of the major results is presented in Sect. 5.

2. Observation and data reduction

The observations of the J = 1 − 0 lines of the CO isotopologues were performed using the PMO 13.7-m millimeter telescope (Zuo et al. 2004, 2011) from September to December 2009. A superconductor-insulator-superconductor (SIS) receiver with three acoustic-optical spectrometers (AOSs) working at 110.20 and 115.27 GHz covered 13CO\hbox{${\rm ^{13}CO}$} and 12CO\hbox{${\rm ^{12}CO}$}(1 − 0) lines respectively. Every spectrometer had 1024 channels. The bandwidths for 12CO and 13CO lines were 145.330 and 42.762 MHz, corresponding to velocity resolution of ΔV = 0.37 and 0.11 km s-1, respectively. The weather conditions had been good. The system temperature during observations was between 150 and 300 K. The root mean square (rms) noise level per channel was 0.4 and 0.3 K for the 12CO and 13CO bands, respectively. The full-width half-maximum (FWHM) beam width was ~ 57′′ (Zuo et al. 2011). Two objects, W51D or S140 was observed once per hour as flux calibrators (one of them closer to the target source would be adopted). The main-beam brightness temperature of the calibrators (Tmb, calibrator) were estimated from the recently released beam efficiency1, which is ηmb = 0.40 to 0.48 (varying with the frequency and elevation). After correction for ηmb, for S140, the intensity is T(12mbCO)=68.8\hbox{$T{\rm _{mb}(^{12}CO)}=68.8$} K and T(13mbCO)=29.0\hbox{$T{\rm _{mb}(^{13}CO)}=29.0$} K. For W51D, the values are T(12mbCO)=91.7\hbox{$T{\rm _{mb}(^{12}CO)}=91.7$} K and T(13mbCO)=37.5\hbox{$T{\rm _{mb}(^{13}CO)}=37.5$} K. By comparing the spectra of the calibrators at our observing period, we can obtain the correction factor C = Tmb,calibrator/Tant,calibrator and use it to correct the intensity of the target source, which is Tmb,target = C·Tant,target. For S140, its observed T(12antCO)\hbox{$T{\rm _{ant}(^{12}CO)}$} has a daily variation of ± 7% (as due to the elevation change from 45° to 65°), which is adopted as the uncertainty level in the flux calibration. We note that the AOSs also covered the 18CO (1 − 0) line (f = 109.78 GHz). However, spectra often have a low signal-to-noise ratio or fluctuated baseline and are therefore not adopted for the subsequent analysis.

The position-switch mode was used for the observation. The reference positions were examined to have no detectable emission feature within the observed bandwidth. Otherwise, if there is an emission feature, it should be separated from the line profile of the target. For the mapped regions, the mapping steps along the right ascension (RA) and declination (Dec) directions are both 1 arcmin. The integration time was 1 min for each point. The typical angular extent required to map the 13CO emission region is 7 × 7 arcmin2.

Another fraction of the observational data was obtained from the public 13.7-m telescope database2. These data originally came from the observations of the star-forming regions selected from other tracers and Galactic CO line surveys. The 12CO and 13CO (1 − 0) emission within 1.0-arcmin vicinity around the maser sources would be adopted for our study. With the archive data, we have found CO spectra for 160 methanol masers whose Galactic longitudes range from l = 9.6° to 232.6°. All the data were calibrated using the same procedure as that performed for our own observation. The GILDAS software package3 was used for data reduction and plot. The program GraphPad Prism4 was also used for some statistical test and plot. The observed line parameters are shown in Table A.1 in the online material, and the derived physical parameters are shown in Table 1.

3. Results

3.1. Selection for the valid CO line components

All the maser sources in our sample are detected in 13CO (1 − 0) emission with the 13.7-m telescope. It should still be further examined whether the 13CO gas is physically associated with the maser or is merely close in projection. This can be done by comparing the systemic radial velocity between the CO-line profile and the maser emission. In our sample, the velocity difference between the 13CO line and the maser is mainly concentrated in a range from 0 to 5 km s-1, with a small fraction of the sources extending to higher values.

If the maser source and the CO line have a small velocity difference, they should be spatially associated, while larger discrepancies may indicate them to be separated. A threshold velocity difference is needed to distinguish between the two cases. Using the PMO 13.7-m telescope, Guan et al. (2008) measured the 13CO velocity field in a number of molecular clouds and found a spatial variation from 25 km s-1. Taking this result into account, we adopted a selection criteria including three cases:

  • (1)

    if \hbox{$|V_{\rm sys,maser}-V_{\rm sys,^{13}CO}|\leqslant 5$} km s-1, the 13CO emission would be assumed to be associated with the maser source;

  • (2)

    if \hbox{$5<|V_{\rm sys,maser}-V_{\rm sys,^{13}CO}|\leqslant 10$} km s-1, an association is still assumed but labeled as questionable (“q”) in Table 1;

  • (3)

    if | Vsys,maserVsys,13CO | > 10 km s-1 the CO data would be excluded from our sample.

    Table 1

    Physical parameters of the maser sources.

    thumbnail Fig. 1

    Number of sources in each 10-degree Galactic longitude interval. The gray and black-line columns represent the full sample of the methanol masers (Pestalozzi et al. 2005) and the fraction covered by our sample, respectively. Overlaid inset panel: distributions of the maser intensities for the sampled and missed sources in l = 10° to 50° range. For each statistical box, the entire error bar (whisker) represents the overall distribution, the box and its inner segment represent the 25% to 75% range and the median value, respectively.

Based on this selection criteria, 26 spectra lines have been excluded. We note that this criteria may not be fully adequate or complete. First, even if the CO lines and the maser source have similar velocities, the projection effect still exists. Second, another fraction of sources with a large velocity difference may be mistakenly excluded. After the velocity-based selection, 26 13CO components have been excluded, which take up a relatively small fraction of the total number (203 components). Moreover, methanol masers are exclusively inhabited in the young high-mass star-forming sites which are abundant in molecular species. In this case, if a CO component coincides with it in both position and velocity, it should be less likely for them to be separated. A more accurate selection criteria is still to be considered. This is beyond the scope of the current work.

Taking the multiple components into account, we obtained 177 validated 13CO gas components in 160 methanol maser sources in our sample. In the same Galactic longitude range there are 259 masers altogether in Pestalozzi et al. (2005)’s catalogue. In some regions, there are multiple maser spots located in a close vicinity (within 1 arcmin) with similar velocities. If these nearby maser spots are regarded as one single source, the total number of the maser sources would be 205. Then our spectra have covered 79 per cent of these sources. The coordinates and systemic velocity of all the sources and components are presented in Table 1. In Fig. 1, we present the source distribution in each Galactic longitude interval for the maser sources and the fraction included in our sample. As a major insufficiency, 38 maser sources in the Galactic range from l = 10° to 50° (26 per cent of all) are not covered in our observation. In the catalogue of Pestalozzi et al. (2005), the only parameter that can characterize the 6.7 GHz maser properties is their emission intensities. In the overlaid panel in Fig. 1, we show the intensity distributions of the sampled and missed maser sources in this Galactic range. The two groups exhibit similar ranges and average levels for the maser intensities. The “missed” group tends to have slightly lower intensities than the “sampled” group. This indicates that our sample has a slightly better coverage for the strong masers. This is unlikely to significantly affect the sampling unless those missed ones have particularly special properties among the entire sample. In general, we suggest that our data should have an unbiased sampling for the 6.7 GHz maser sources in the corresponding Galactic longitude range.

3.2. Calculation of the physical parameters

We derived the excitation temperature from the 12CO (1 − 0) lines by assuming a high optical depth and a local thermal equilibrium (LTE), using Tex=k[ln(1+/kTMB(12CO)/ηbf+Jν(Tbg))]-1\begin{equation} T_{\rm ex}=\frac{h\nu}{k}[\ln(1+\frac{h\nu/k}{T_{\rm MB}({\rm ^{12}CO})/\eta_{\rm bf}+J_{\nu}(T_{\rm bg})})]^{-1} \end{equation}(1)where TMB(12CO) is the main-beam brightness temperature of the 12CO (1 − 0), Tbg = 2.73 K is the brightness temperature of the cosmic microwave background. Jν(T) = (/k)/ [exp(/kT) − 1], and ηbf is the beaming filling factor. The variable k is the Boltzmann constant, and h is the Plank constant. As shown in the mapped regions (Fig. A.1 in online material) in most cases, the 13CO emissions extend to several arc minutes (Fig. A.1). In this case, ηbf ≃ 1 should be a reasonable assumption.

The 13CO optical depth is derived using (Garden et al. 1991) TMB(12CO)TMB(13CO)1exp(τ12)1exp(τ13),\begin{equation} \frac{T_{\rm MB}({\rm ^{12}CO)}}{T_{\rm MB}({\rm ^{13}CO})}\approx\frac{1-\exp({\tau_{12}})}{1-\exp({\tau_{13}})}, \end{equation}(2)wherein we assumed an optical depth ratio between 12CO\hbox{${\rm ^{12}CO}$} and 13CO\hbox{${\rm ^{13}CO}$} equal to their average Galactic abundance ratio, which is τ12/τ13 = [13CO] / [12CO] = 50. The 13CO\hbox{${\rm ^{13}CO}$} column density is calculated from N(13CO)=3k8π3Bμ2exp[hBJ(J+1)/kTex]J+1×Tex+hB/3k1exp(/kTex)τvdv,\begin{eqnarray} N({\rm ^{13}CO})=\frac{3k}{8\pi^3B\mu^2} \frac{\exp[hBJ(J+1)/kT_{\rm ex}]}{J+1}\nonumber \\ \times\frac{T_{\rm ex}+hB/3k}{1-\exp(-h\nu/kT_{\rm ex})}\int \tau_v \, {\rm d}v, \end{eqnarray}(3)where B = 54.89 GHz is the rotation constant of 13CO\hbox{${\rm ^{13}CO}$}, μ = 0.11079 Debye is the permanent dipole moment and J = 0 is the rotational quantum number of the lower state. The H2 column density is calculated from N(H2) = [H2] / [13CO] × N(13CO). In the calculation, we assumed a 13CO\hbox{${\rm ^{13}CO}$} abundance of [13CO] / [H2] = 1.0 × 10-6, which is also adopted by the Du & Yang (2008) for the Galactic IRDC survey using the 13.7-m telescope. The error in N(H2) is estimated from the errors in T(12mbCO)\hbox{$T{\rm _{mb}(^{12}CO)}$} and T(13mbCO)\hbox{$T{\rm _{mb}(^{13}CO)}$} (both ± 7%). The calculated Tex and N(H2) are presented in Table 1.

3.3. Determination of the distances

The distances to the sources are important for revealing their Galactic distribution and examining whether the measurement of the physical parameters is biased due the distance effect. For each source in our sample, we have searched in literature for the distance measurement. As a result, a few individual sources have a direct distance measurement from interstellar extinction, trigonometric parallax, or the luminosity of the exciting stars in the nearby (UC) Hii regions. For a major fraction of the maser sources, however, we have to rely on the Galactic rotation curve to derive the kinematic distances Dkin. The parameter Dkin is calculated from the Galactic rotation curve of Brand & Blitz (1993) and using the newly measured solar distance R0 = 8.3 ± 0.23 kpc, and the velocity Θ0 = 239 ± 7 km s-1 (Brunthaler et al. 2011). For all the valid 13CO components (including the “questionable” ones), the velocity difference between the maser and 13CO line and between the different 13CO line components are all assumed to reflect the variation in the same region. Consequently, an average radial velocity from the methanol masers and all the 13CO components is adopted to calculate Dkin. In calculation, most sources in the first Galactic quadrant are subjected to the distance ambiguity, where there are two possible solutions for a given radial velocity (commonly known as near and far distances). There are several properties which can assist to resolve the ambiguity, including:

  • (1)

    If the CO emission is associated with a (UC) Hii region, which is observed in 21-cm HI absorption, and if the HI absorption has a systemic velocity larger than the maser source, the maser sources should be located at the far distance (Pandian et al. 2008; Anderson & Bania 2009; Urquhart et al. 2011).

  • (2)

    If the CO emission is associated with the 21-cm HI self absorption in both position and velocity, it should be at a near distance. In contrast, if the CO emission is associated with the HI emission, far distance would be adopted (specified in Roman-Duval et al. 2009, Figs. 1 to 3 therein).

  • (3)

    If the maser source overlaps with an infrared absorption feature, it is then assumed to be at the near distance, considering that there is a small possibility to see the CO gas, which is not associated with the IR absorption but located further behind.

  • (4)

    If the three properties above are unavailable, the one distance value closer to the nearby spiral arm of the two would be adopted. This strategy was used in the previous studies, including Schwartz et al. (1988) and Faúndez et al. (2004). For example, G60.56-0.17 has two distance values of Dkin = 0.3 and 8.0 kpc. Its near distance is even shorter than the Orion molecular cloud which is known to be the nearest high-mass, star-forming site (414 pc, Menten et al. 2007). Instead, if we assume the far distance (8.0 kpc), it would be well located on the Perseus Arm.

The distance determination and the references for all the sources are specified in Table 1. The major uncertainty in Dkin and Rgc is from the uncertainty in the fitted Galactic rotation curve (Brand & Blitz 1993). From Fig. 4b in Brand & Blitz (1993), we can see that the rotation velocity has a deviation from the fitted curve with a scale of ± 20 km s-1 for Rgc< 8.5 kpc. Based on this feature, we estimate the error by varying the Galactic rotation curve with this value of velocity. For Rgc> 8.5 kpc, the rotation curve begins to have large uncertainties (>100 km s-1). To give a more reasonable estimate, we assume that the distance error of a source should be less than 1/2 times the interval between two nearby spiral arms along its Galactic longitude, as measured from Fig. 2.

thumbnail Fig. 2

Distribution of our observed sources on the Milky Way disk. The source sizes are proportional to their 13CO\hbox{${\rm ^{13}CO}$} line intensity. The thick faint blue lines represent the distribution of the spiral arms (Russeil 2003) in our observed field. The background image is the artist’s view of the Galactic morphology (R. Hurt: NASA/JPL-Caltech/SSC, with a more detailed introduction in Urquhart et al. 2011).

thumbnail Fig. 3

Distribution of the physical parameters with the heliocentric distance. Left panels: all the data points. Right panels: statistical box plots for the data distribution in each 2-kpc interval. The dots represent the outliers, which are apparently isolated from the major fraction of the data points. For each statistical box, the entire whisker represents the overall distribution. If there are isolated data points (based on eye selection), the whisker then represents the remaining data range. The box and its inner segment represent the 25% to 75% range and the median value, respectively. We selected interval pairs, which may typically reflect the variation trend. The interval pairs are indicated with gray segments above the whiskers. Labeled on the segment are the p values from the t test and M-W test respectively, and the test result (“yes” or “no” in the brackets) for whether the two intervals have significantly different distributions.

3.4. Distance-related errors in the physical parameters

The large variation in our source distance may potentially influence our sampling and the measurement of the physical parameters. Towards longer distances, the beam area covers a larger spatial region, leading to two main effects: (1) the emission from a region smaller than the beam size would be weakened due to the beam dilution; (2) the different gas components with the beam area would all be averaged together. The second effect may cause ΔV(13CO) to be overestimated at larger distances. To investigate the significance of the “distance effects”, we plot the variation of the physical parameters including Tex, N(H2) and ΔV(13CO), as a function of the Heliocentric distance in Fig. 3. We performed a linear least-square fit for each parameter. As a result, the variation of the physical parameters as a function of the distance is insignificant with a correlation factor of <0.2.

As another examination, which can somewhat eliminate the errors in the distance measurement, we divided the data points into each 2-kpc interval from short to long distances and compare the statistical distributions in different intervals. The results are shown in the right panels in Fig. 3. For each parameter, the data points in different intervals appear to be similar with each other. As quantitative examination, we carry out both the Student’s t-test and the Mann-Whitney (M-W) U test. Both tests provide an estimate of whether one data set tends to have larger or smaller values than another one. The difference is that the t-test assumes a Gaussian distribution for the data set, while the M-W test assumes a random distribution. In Fig. 3, we show the test results (p-values for the t-test and the M-W test) between some intervals, which are typical to illustrate the variation trend. The major results for the three parameters are as following:

  • (1)

    For the Tex distribution, the 1012 kpc interval shows significantly smaller Tex compared to the other intervals. For example, it is different from the distribution in the 02 kpc interval with p = 0.004 and 0.003 for the two tests. This suggests that the beam dilution may have weakened the distant sources. All the other intervals show similar Tex distributions in both M-W and t tests (p> 0.05) with none of them having significantly higher or lower values than the other intervals. This indicates that the beam dilution, if exists, may have only slightly influenced the Tex measurement towards the most distant sources.

  • (2)

    The N(H2) (H2 column density) distribution shows a similar result with that for Tex, that is, the 1012 kpc interval tends to have smaller N(H2) values than the other intervals (with p = 0.001 and 0.005). This may also be explained as the beam weakening to the distant sources.

  • (3)

    For the ΔV(13CO) distribution, the 46 and 68 kpc intervals show higher ΔV(13CO) than the other intervals, as indicated by the both tests. The other intervals have similar ΔV(13CO) distributions with each other. From a theoretical point of view, if ΔV(13CO) is subjected to the distance effect, the observed ΔV(13CO) should roughly increase linearly at longer distances, because the beam area would statistically cover more gas components at larger distances. The observed trend is different. It shows a decreased ΔV(13CO) in 812 kpc compared to the 46 and 68 kpc intervals. From Fig. 2, we actually found that a large fraction of the sources in this distance range are located around the 5-kpc molecular ring and the edge of the central Galactic bar. These regions may actually have increased gas turbulence.

We note that the resolved distances are still suggestive more than definitive. If the ambiguity still exists, the distance effect can be evaluated at two extreme cases, which assumes that all the sources have near and far distances, respectively. The first case shows a very small difference compared to Fig. 3 mainly because the near distances are already adopted in most cases. In the second case, which totally adopts the far distances, the resultant parameter distributions are shown in Fig. A.5 online material. The Tex distribution appears to have a decline from D = 0 to 10 kpc. The trend is more clearly shown on the statistical bars (right panel). However, the values rise again from D = 10 to 18 Kpc, so that Tex tends to have the smallest values in 810 kpc interval. A similar trend is also seen in the N(H2) distribution (middle-left and -right panels). The figure shows that the 810 kpc interval has the lowest N(H2) on average, as well as the largest distribution range. This feature is inconsistent with the expected beam dilution at larger distances which is expected to cause a linear decrease in Tex and N(H2). Instead, it may reflect a situation that many sources are mistakenly put into the 812 kpc interval. In the ΔV(13CO) distribution (bottom panels), there is a slight increase for the average level with the distance similar as shown in Fig. 3. In the mean time, 812 kpc intervals also show an abnormally crowded population compared to the other intervals, suggesting that the distances may be unreasonably determined.

As a summary for the analysis above, Fig. 3 shows that the beam dilution might slightly weaken the intensities for the sources at large distances, and the beam convolution might somewhat blend the different velocity components, but these effects are still weak. It is unlikely that the physical parameters have a noticeable systemic bias as due to the distance effect.

3.5. Emission features in the mapped regions

3.5.1. Identification of the molecular cores

The 13CO (1 − 0) maps of all the observed regions are presented in Fig. A.1 in the online material. From the 33 mapped regions, we have identified 51 molecular cores altogether which are associated with or close to the maser sites within several arc minutes. Here, we define a 13CO core, as the gas condensation enclosed by the 50-per cent contour in the 13CO (1 − 0) integrated intensity map. In most regions the molecular cores can be well identified from the 13CO contours. There are a few sources with extended emission morphologies or overlapped multiple sources. In those regions we use ellipses to encircle the plausible 13CO cores. The ellipse is drawn to overlap with the part of the 50% contour that is not contaminated by the nearby cores or diffuse emissions. The sources with contamination include G14.33-0.64, G29.86-0.05, G30.82-0.05, G35.20-0.74, G35.20-1.74, G80.85+0.43, and G136.84+1.12.

We measured the effective radius of the core using πRcore2=Acore\hbox{$\pi R_{\rm core}^2=A_{\rm core}$}, where Acore is the core area encircled by the 50 % intensity contour level in units of pc2. We note that the core radius is subjected to a large uncertainty mainly because the spatial resolution largely differs with the distance. The values are therefore listed in Table A.3 (online material) only for reference and not adopted for further quantitative analysis. The Tex and N(H2) of the molecular cores are estimated from the spectra at the position closest to the core center, and the values are also presented in Table A.2.

3.5.2. Infrared and radio continuum emissions

To investigate the infrared (IR) emission of our molecular cores, we looked for the images from the database of GLIMPSE, MSX, and IRAS sky surveys5. In Fig. A.1, the 13CO (1 − 0) contours are overlaid on the infrared image (gray scale). The IRAC 8 μm images from the GLIMPSE survey would be first considered if they are available. If a region has no GLIMPSE image, the MSX 8.28 μm data would be used instead. There is only one source (G213.70-12.6), which has neither MSX nor IRAC data. For this source, we adopted the IRAS 12 μm image to show the IR emissions. The IR emissions of our sources exhibit various morphologies. In some regions such as G10.62-0.38 and G12.79-0.19, the IR emissions clearly reveal multiple IR sources. In other regions, such as G15.03-0.68 and G30.82-0.05, the IR emissions are very extended and do not reveal any point-like sources. A major fraction of the IR-faint cores are actually not associated with the masers but the “additional” objects in the same region. Among the 33 cores, which coincide with the masers, there are only one source (G173.71) with no evident IR emissions. This suggests that the stellar heating is already significant and often generates extended IR emissions in most maser sources.

Among the 30 13CO cores with noticeable infrared emissions, G8.68-0.37 and G15.03-0.68 have IR emissions at the edge of the 13CO cores, and have much weaker IR emissions toward the core center. The observing center (zero offset) of each region is located at the mazer point which is frequently but not always coincident with the IR emission peak. In the regions of G29.86-0.05, G35.20-0.74, G045.44+0.07 and so on, the maser site does not have intense IR emission. Since the IR emission is necessary for the pumping the Class-II methanol masers (Reid & Menten 1991), there may still be stellar emissions in those cores, which are probably highly obscured by the dust.

To inspect the possible radio continuum emission in the molecular cores, we referred to the MAGPIS6 6 cm survey and other literature. The MAGPIS covers a Galactic longitude range of 350°<l< 49.5° with a detection limit of 2.3 mJy. There is another unbiased 6 cm survey performed by Gregory & Condon (1991), which covers the sky field between 0°< Dec < 75° with a detection limit of ~5 mJy. The combination of MAGPIS and Gregory’s observation, have well covered our mapped regions. In addition, we also sought for other radio continuum observations for our fields. If a radio continuum source is located within the 50% contour of a 13CO core, it is then regarded to be associated with the core. The pointing accuracy of all the free-free emission surveys is better than 10 arcsec. The radio continuum fluxes for the 13CO cores are listed in Table A.2.

In our mapped region there are 17 cores which are likely associated with radio continuum emissions. In these regions, the radio emissions all coincide well with the IR emission, indicating that it should originate from the ionized gas near the young massive stars. The radio continuum intensities are presented in Table A.2. If a core does not have radio emission as shown in MAGPIS, Gregory & Condon (1991), or other published observational studies, it would be regarded to be radio-faint. We note that the detection of the radio continuum highly depends on the observational sensitivity. For instance, G188.79+1.02, G188.95+0.89, G189.03+0.76, and G232.62+0.99 are below the detection limit of Gregory & Condon (1991). However, they actually have faint emission around 1 mJy (see Table A.3). However, due to their weakness, we assume these cores to be radio-faint. The locations of radio continuum sources are labeled with open squares in Fig. A.1. With the radio-bright sources excluded, there are 18 masers sources, or 54% that are left, which are likely pre-UC Hii due to weak or undetected radio emissions.

3.6. Infall and outflow signatures

As an indication of high-velocity gas, some cores show plausible non-Gaussian line wings in 12CO (1 − 0) spectra. To further identify the potential outflow motion originating in the 13CO\hbox{${\rm ^{13}CO}$} cores, we made p-v plot across each core along both the RA and Dec directions. If the contours on the p-v plot shows a significant swell toward the blue- and/or redshifted velocities, it is then suggested to be an outflow candidate. Using this selection rule, high-velocity wings are found in six of the 51 cloud cores: G29.86-0.06(1), G35.20-0.74(1), G45.44+0.07, G59.78+0.06(1), G94.58-1.79, and G173.49+2.42. The PV diagrams and spectra at the core center are shown together in Fig. A.2. The outflow in G35.20-0.74 has been previously detected in CO (2 − 1) (Wu et al. 2005). On the other hand, we did not observe evident line wings in some sources with outflows detected, including G8.68-0.37 (Longmore et al. 2011), G10.62-0.38 (Liu et al. 2010a), G10.47+0.03 (Gibb et al. 2004), and G31.41+0.31 (Olmi & Cesaroni 1999). In these regions, the outflow may have much smaller spatial scales, as in G10.62-0.38 where the outflow shows up as compact gas condensations at 10-arcsec scales (Liu et al. 2010a). These outflows should have been diluted or contaminated with the gas turbulence on our observational scale due to the dilution by the large beam size. In Fig. A.3, we present four sources which have relatively well-shaped outflow lobes. In the other two outflow candidates, the outflow lobes are largely contaminated with the surrounding gas and cannot be clearly identified. Their maps are thus not shown in the figure.

Since class-I methanol masers show a broad association with the outflow and shock emissions (Chen et al. 2009, and the references therein), we also examined the association between the class-I masers and the outflows. As a result, we found that two outflow candidates, G45.44 and G59.78, are indeed associated with the class-I masers. The non-detection of outflow in other class-I masers may also be tentatively attributed to the beam dilution and gas contamination. In the recent extensive investigation for class-I masers (Val’tts & Larionov 2007), 25 per cent of the class-I sources are found to have outflows. This detection rate is much higher than that shown in our sample. In our sample, the low outflow detection towards class-I sources (75% of all) should be due to the small-scale outflows being less well traced by the CO (1 − 0) emission or diluted by our beam.

Plausible infall signature, which is a 12CO profile with blueshifted absorption and a single peak 13CO line (Zhou et al. 1993), are detected in four cores. We estimated their infall rate from the 12CO line profiles. We first use Eq. (9) in Myers et al. (1996) to derive the typical infall velocity from the 12CO line profile and then use Eq. (3) in Klaassen & Wilson (2007) to estimate the infall rate. In calculation, the typical infall radius (Rgm) is adopted as the core radius (Rcore) measured from the 13CO emission region. The derived infall parameters are listed in Table A.3. We note that the low resolution and insufficient sampling may cause great uncertainty in Rcore and, hence, the infall rate. The infall rates are presented only as tentative results and comparison to subsequent studies.

Klaassen & Wilson (2007) observed 23 massive star-forming regions in HCO+ (4 − 3) and detected infall motion in nine objects in the sample. The derived infall rates are on average around 10-4M year-1, which are much smaller than the tentative measurement here. In G8.68-0.37 where the infall signature is detected in both CO and HCO+, we estimated an infall velocity of 0.7 km s-1, which appears consistent with their value, and our different infall rate is mainly due to the large core radius estimated from the low-resolution 13CO map. Except for G34.26, all the infall candidates in Klaassen & Wilson (2007), including G8.68, G10.47, G15.03, G109.86, G133.94, and G192.6 are also in our sample. Besides G8.68-0.37, our CO spectra do not show evidence of infall signatures in any of the other sources. The object G35.20-0.74 was detected with a blue profile in CO (2 − 1) line (Wu et al. 2005), which is not evident in the current CO (1 − 0) line. The possible reason is that the two molecular species are tracing different gas components. The CO, which mainly traces the low-density gas, may have a flatter temperature profile, so that self-absorption is not significant. The difference between the CO and HCO+ may imply an inside-out mode for the core collapse in these regions, which can be tested by future high-resolution studies.

3.7. Comments on individual objects

The properties of the sources noteworthy are described as follows:

  • 1.

    G8.68-0.36. This region contains a single compact 13CO core. The UC Hii region G8.68-0.36 (Wood & Churchwell 1989) is associated with the strong IR emission in the IRAC 8 μm bands. The UC Hii region is located west to the 13CO core center. Recent high-resolution observation revealed a strong bipolar outflow originated from the core (Longmore et al. 2011).

  • 2.

    G10.47+0.03. This region has a single compact 13CO core. The core is elongated from the southeast to northwest. The densest part of the molecular core is associated with an extended 8 μm emission, which shows a clustered environment at the center of the molecular core. The high-resolution study (Pascucci et al. 2004) reveals a group of young high-mass stellar objects at different evolutionary stages.

  • 3.

    G10.62-0.38. This field contains a roughly spherical 13CO core. The IRAC 8 μm image shows several stellar sources concentrated at the center and diffuse emission stretching to the outer part, and a cluster is being formed at the center. Ionized accretion flow was detected in the UC Hii region (Keto 2002; Sollins & Ho 2005). High-resolution study revealed a complicated outflow system arising from the stellar objects at the center (Liu et al. 2010a).

  • 4.

    G12.79-0.19. This region shows a compact 13CO core, which is slightly elongated from northeast to southwest. The 12CO (1 − 0) line at the center shows a blue profile, which indicates the infall motion (Fig. 3). It also has a broad emission feature extended to the redshifted side of the line profile (Vlsr = 40 to 50 km s-1). However, this emission feature also show up in the 13CO (1 − 0), indicating a considerable optical depth. It is therefore more likely to be a different velocity component rather than merely the outflow. In the IRAC 8.0 μm image, we can see three bright stellar sources at the core center with one of them coincident with the radio continuum source. In their vicinity, there are several extended filamental structures which extend further away from the center. They can be hot gas and dust expanding from the center, or the ambient gas compressed and swept by the stellar emission.

  • 5.

    G14.33-0.64. This region shows extended CO emission. The 13CO emission is resolved into four major dense cores. Core-1 and core-2 are not fully separated from each other, but their IR sources are well separated. Thus, we expect them to be different molecular cores. Core-3 and 4 are free of IR emissions. The MAGPIS 6 cm survey shows no detectable radio continuum emission in this (entire) region.

  • 6.

    G29.86-0.05. In this region, the bulk of the CO emission is located to the north of the center (i.e., the CH3OH maser location). Therefore the methanol maser may arise from a less extended core, which is currently unresolved. The emission region is extended and resolved into four major dense cores. Core-1 is associated with a bright IRAC point source. It also drives an outflow, as shown by the broad line wings (Fig. A.1). The remaining three cores are much fainter in the CO (1 − 0) and IR emissions.

  • 7.

    G30.82-0.05. This region exhibits an extended CO emission feature. We have identified three 13CO cores corresponding the strongest emission peaks. We examined the C18O emission and confirmed the location of the three strongest cores. The cores all reside at the edge of the IR emission region. The methanol maser is associated with core-1. Core-2 and core-3 have no methanol masers but are detected in the radio continuum emission.

  • 8.

    G31.41+0.31. This source is a frequently studied high-mass star-forming region. In our observation, it shows a quite extended and irregular emission region with a scale of ~6′ (9 pc). It may therefore contain many unresolved sub-structures. Nevertheless, the 13CO emission peak coincides well with the methanol maser and the hot molecular core observed at higher resolutions (Girart et al. 2009, and reference therein.). This core has a CH3CN rotational temperature of ~164 K (Beltrán et al. 2005), which is much higher than the 12CO excitation temperature (17.5 K). The difference should be mainly because the CO (1 − 0) intrinsically traces colder gas components than the higher-excited lines and the high-density molecular tracers. This should be universal from the highest Tex, the comparison between the sources at different evolutionary stages suggests that the stellar heating can still noticeably affect the CO temperature (see Sect. 4.2).

  • 9.

    G35.20-1.74. In the MSX 8 μm image, this region shows a prominent contrast between the bright IR emission in the east and the broad dark region in the west. The 13CO (1 − 0) line reveals two dense cores near the IR emission. There is also a large area of extended CO emission in the west, which is coincident with the IR-dark region, wherein we do not identify any well-shaped cores. There is a strong radio continuum emission (Gregory & Condon 1991), as marked with the red square. Since the radio continuum sources are located at a position, where the two cores are overlapped, we assume that the both cores are associated with the radio emission. Core-1 exhibits a self-absorption in 12CO (1 − 0), which may indicate an infall motion.

  • 10.

    G080.85+0.43. There are four dense cores in an alignment from northwest to southeast. Two of them are associated with the MSX point sources. Core-1 is detected in radio-continuum emission. The other cores are radio-quiet.

  • 11.

    G183.34+0.59. This region shows a relatively compact core which is associated with a weak MSX point source. A strong MSX source is located on the western edge of the core, which also shows bright radio continuum emission. The 13CO emission shows a sharp decline toward the IR source, indicating that the molecular gas is being destroyed in the UC Hii region.

thumbnail Fig. 4

Distribution of the physical parameters with the Galactocentric distance. The first panel shows the data points for the 12CO\hbox{${\rm ^{12}CO}$} excitation temperature. The other panels show the statistics of the parameters in each 2-kpc interval. For each statistical box, the entire whisker represents the overall distribution. If there are isolated data points (based on eye selection), the whisker then represents the remaining data range. The box and its inner segment represent the 25% to 75% range and the median value, respectively. We selected interval pairs which may typically reflect the variation trend. The interval pairs are indicated with gray segments above the whiskers. Labeled on the segment are the p values from the t test and M-W test respectively, and the test result (“yes” or “no” in the brackets) for whether the two intervals have significantly different distributions.

thumbnail Fig. 5

Variation of the 12CO\hbox{${\rm ^{12}CO}$} excitation temperature and the 13CO\hbox{${\rm ^{13}CO}$} line width as a function the H2 column density. The dashed lines represent the best-fit linear correlation.

4. Discussion

4.1. Variation with the Galactocentric distance

In Fig. 4 we plot the variation of ΔV(13CO), N(H2), and T(12exCO)\hbox{$T{\rm _{ex}(^{12}CO)}$} as a function of the Galactocentric distance Rgc for our maser sources. We note that Rgc always has a single solution and is unaffected by the distance ambiguity. As a major result, the line width shows a noticeable decrease with Rgc. At Rgc = 4 kpc, the average line width is ΔV = 8 km s-1, and it drops to ΔV = 2.5 km s-1 at Rgc = 9 kpc. The two statistical tests show that the different intervals vary with each other with p< 0.05 in most cases. The ΔVRgc relation may suggest that the gas becomes more turbulent in the maser sources closer to the Galactic center. The inner Galaxy on average has a higher star-and-gas density than the outer region. In a denser environment, the molecular cores can become more turbulent when affected by two aspects: a more frequent dynamical interaction with nearby molecular cores or clouds and the energy input from the stellar wind and radiation. Besides, the line broadening due to the optical depth should also be considered. Assuming that the optical depth is a Gaussian-shaped function of the radial velocity (Myers et al. 1996), τ(V) = τ0exp [ − (VVsys)2/ 2σ2], where τ0 is the optical depth at the systemic velocity Vsys and σ is the velocity dispersion along the line of sight and can be measured from the intrinsic line width using σ=ΔV/8ln2\hbox{$\sigma=\Delta V/\!\sqrt{8\ln2}$}. The observed line profile would be Tmb(V) = Jν(Tex) [1 − exp( − τ(V)]. From these equations, the ratio between the observed line width (ΔVobs) and the intrinsic line width (broadening factor) is derived to be ΔVobsΔV=[lnτ0/ln[2/(1+eτ0)]ln2]1/2·\begin{equation} \frac{\Delta V_{\rm obs}}{\Delta V}=\left[\frac{\ln{\tau_0/\ln[2/(1+{\rm e}^{-\tau_0})] }}{\ln 2}\right]^{1/2}\cdot \end{equation}(4)For all our 13CO components, the broadening factor in 13CO (1 − 0) turns out to be less than 1.5. Thus its contribution to ΔV13 and the ΔV13Rgc correlation would not be significant, and we suggest that the cores at smaller Rgc may indeed have larger velocity dispersions. A similar trend between ΔV and Rgc was also suggested for the CO gas in the IRDCs (Du & Yang 2008).

The other two parameters, N(H2) and T(12exCO)\hbox{$T{\rm _{ex}(^{12}CO)}$} do not show significant Galactic variation. For the N(H2) distribution, the 24 and 1012 kpc intervals have very large probabilities to have similar distributions. The Tex distribution shows a slight increase toward large Rgc with the difference reaching the significant level (p> 0.05) between 02 and 1012 kpc intervals. However, its variation scale is still much weaker than the trend in ΔV.

thumbnail Fig. 6

Statistical comparison for the excitation temperatures between the different subsamples of the mapped cores. Left panel: the radio-bright and radio-faint cores. Right panel: the IR-bright and IR-faint cores. For each data set, the whisker represents the overall distribution; the box and inner segment represent the 25% to 75% range and the median value, respectively.

4.2. Correlation between the physical parameters

The relation between the different physical parameters are shown in Fig. 5. Figure 5a shows that the 12CO excitation temperature significantly increases with the H2 column density (with a slope of 3.3 ± 0.3 and a correlation factor of R = 0.4). Krumholz et al. (2010) performed a hydrodynamic simulation about the molecular core collapse and star formation therein. The simulation shows that a higher surface density can significantly facilitate the formation of high-mass stars (M> 10 M), thus lead to an increased temperature in the natal gas-and-dust core. In their densest case of Σ = 10 g cm-2 or N(H2) = 2 × 1024 cm-2, the final temperature in the surrounding gas exceeds 1000 K as heated by the newly formed massive stars. In contrast, only a group of low-mass stars can be formed in the low-density case of Σ = 0.1 g cm-2, and the final temperature is ~100 K. The simulation result shows a reasonable consistency with the observed TexN(H2) correlation. Namely, the formation rate of the massive stars and/or heating efficiency would increase in higher-N(H2) regions.

The above explanation have some apparent uncertainties. The major question is whether the internal stellar heating can noticeably change T(12exCO)\hbox{$T{\rm _{ex}(^{12}CO)}$}. Due to the large optical depth in 12CO, its surface temperature might be insensitive to the internal stellar heating and instead, more affected by the external emissions. Based on our observational data, there is one way to examine this problem. Our mapped regions have provided a number of 13CO cores that are not associated with radio continuum and/or IR emissions. These cores can naturally serve as a control group for us to evaluate the significance of the internal heating.

We first compared the temperatures in the radio-bright and radio-faint cores (Fcm< 6 mJy). Their distributions are shown in Fig. 6a. The figure shows that the radio-bright group are apparently inclined to high temperatures. The M-W and t tests both indicate that radio-bright cores have significantly higher Tex than the radio-faint ones. The average temperatures of the two groups are 26.4 K and 20.5 K, respectively. The second comparison is based on their IR emissions, as shown in Fig. 6b. The IR-bright and IR-faint groups have average temperatures of 24 K and 16 K with the contrast being slightly larger. This is within our expectation, since the IR-faint cores are actually a sub-sample within the radio-faint group. The comparison suggests that the cores without IR and radio continuum emissions indeed tend to have the lower T(12exCO)\hbox{$T{\rm _{ex}(^{12}CO)}$}.

In the case that T(12exCO)\hbox{$T{\rm _{ex}(^{12}CO)}$} is sensitive to the internal stellar heating, the observed TexN(H2) correlation would favor the explanation that the stellar emission is actually more intense and/or more efficient in the higher-N(H2) cores. Nevertheless, our current data is still insufficient to confirm the physical connection between N(H2) and star-forming activities. For example, due to the low resolution, we cannot determine whether the IR and radio sources are embedded in or separated from the densest gas in a given core. It calls for denser molecular tracers to more accurately determine the temperature and N(H2) distributions and to better reveal their correlation with the star-forming processes.

The relation between ΔV(13CO) and N(H2) is investigated and shown in Fig. 9. The value of ΔV(13CO) appears to have no significant correlation with the H2 column density. A linear fit provides a slope of 0.28 and the correlation factor of R = 0.09 which are both much smaller than the corresponding values in the TexN(H2) relation. This may be explained as a result from the Larson relation (Larson 1981). Larson (1981) showed that, ΔV is sensitive only to the spatial size, which reflects the nature of the hierarchical turbulence, and in the mean time, the number density is anti-correlated with the spatial size. These two factors would cause the column density to be almost independent on the spatial size, which is consistent with the observed weak correlation between ΔV(13CO) and N(H2).

4.3. Comparison between the maser sources and the infrared-dark cores

Besides the CH3OH masers, the IRDCs also serve as a potential reservoir for high-mass cores at the earliest evolutionary stages (e.g. Rathborne et al. 2006). Du & Yang (2008) have performed a single-point and mapping survey of 61 IRDCs, selected from the prominent absorption feature at the MSX 8-micron band that uses the PMO 13.7 m telescope in the three isotopic CO (1 − 0) lines. Our sample only has two sources (G35.19 and G173.71) in common with Du & Yang (2008). Even in these two regions, our 13CO cores have an offset larger than 5 arcmins from the observing center of Du & Yang (2008). A large fraction of the maser sources are also associated with IR absorption feature, as shown in Peretto & Fuller (2009; see also our Table 1), but do not appear in their IRDC sample probably because the absorption is much less prominent than those shown on the MSX 8 μm band.

In Fig. 7, we plot the distributions of the physical parameters in the maser and IRDC samples. The two samples turn out to have largely different N(H2) distributions. The IRDCs have a major fraction with N(H2) > 4 × 1022 cm-2 and have no distribution below 1022 cm-2. The methanol masers, in contrast, are mainly concentrated in a small N(H2) range, with the distribution largely decreasing at N(H2) > 1022 cm-2. The IRDC sample is selected from the intense absorption features at MSX 8 μm band; thus, they are expected to have large N(H2). Moreover, the difference between the two samples implies that the source selection based on the IR absorption feature could be biased toward high-N(H2) regions, which misses out a considerable amount of high-mass star-forming cores inhabited in the clouds with moderate N(H2).

thumbnail Fig. 7

Comparison of the physical parameters in our maser sources and the IRDCs (Du & Yang 2008). The left panels show the parameter distribution in each interval. The right panels show the overall statistics. For each data set, the whiskers represent the overall distribution; and the box and its inner segment represent the 25% to 75% range and the median value, respectively.

As shown in the middle and bottom panels in Fig. 7, the two samples have similar Tex(12CO) and ΔV(13CO) distributions much more similar, except that the maser sample has larger dispersion in Tex(12CO), which extends from 5 to 40 K. In comparison, the IRDCs show a concentration between 10 and 30 K. The possible explanation is that the methanol masers sources are generally at later evolutionary stages, thus having more various physical conditions. The stellar emission can either increase the temperature or dissipate the molecular gas and drive chemical reactions to consume CO. The latter aspect would cause a decrease in both the optical depth and beam filling factor, and lead to the smaller values in the observed Tex (5 to 10 K), as compared to the IRDCs.

5. Summary

We performed an observational study with the 12CO and 13CO (1 − 0) lines of 160 methanol masers and made mapping observation of a fraction of them (33 sources). We estimated the physical parameters and examined their Galactic variation and statistical properties. Our main results include four aspects:

  • (1)

    We examined physical parameters as a function of the Galactocentric distance (Rgc). As a result, the 13CO line widths show evidence of a decrease from inner to outer Galaxy. This may reflect a feature that star-forming regions toward the inner Galaxy are more turbulent. The other two parameters, Tex(12CO) and N(H2), do not show a prominent variation trend throughout the Galaxy.

  • (2)

    We examined the correlation between the different physical parameters. As a result, the 12CO\hbox{${\rm ^{12}CO}$} excitation temperatures show a noticeable increase with the H2 column densities. This provides consistency with the numerical simulation that suggests an increasing high-mass star formation and heating efficiency of higher-N(H2) regions. This explanation is possible because the IR- and radio-bright cores show higher 12CO\hbox{${\rm ^{12}CO}$} excitation temperatures than those with an absence of heating sources.

  • (3)

    We compared the physical parameters of the molecular cores in our maser sample to those of the infrared dark clouds. The maser sources tend to have noticeably lower NH2 compared to the IRDCs. As an indication, the intense IR absorption features might miss out on some high-mass star-forming cores in low-N(H2) regions. The other parameters show less significant distinctions between the two samples.

  • (4)

    In our mapped regions, we have identified 51 cores from 13CO (1 − 0) emission. Among those cores we found six outflow and four infall candidates. There are 17 cores associated with significant radio continuum emission and 42 cores with significant IR emissions. This indicates that the stellar emissions are already noticeably heating the surrounding gas in the methanol-maser sources, while the ionization is still at a less developed stage.


4

GraphPad Prism version 5.00 for Windows, GraphPad Software, La Jolla California USA, http://www.graphpad.com

Acknowledgments

We are grateful to the staff working in the Qinghai station of the Purple Mountain Observatory 13.7 m telescope for their valuable support to our observations. We would like to thank the anonymous referee for the careful inspection of the manuscript and constructive comments. We also thank the staff who maintain the MIRIAD MacPort package. This work is supported by the China Ministry of Science and Technology under State Key Development Program for Basic Research (973 program) No. 2012CB821802, the National Natural Science Foundation of China No. 11373038, No. 11373045, Strategic Priority Research Program “The Emergence of Cosmological Structures” of the Chinese Academy of Sciences, grant No. XDB09010302, and the Young Researcher Grant of National Astronomical Observatories, Chinese Academy of Sciences.

References

  1. Anderson, L. D., & Bania, T. M. 2009, ApJ, 690, 706 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arzoumanian, Z., Fruchter, A. S., & Taylor, J. H. 1994, ApJ, 426, L85 [NASA ADS] [CrossRef] [Google Scholar]
  3. Beltrán, M. T., Cesaroni, R., Neri, R., et al. 2005, A&A, 435, 901 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945 [NASA ADS] [CrossRef] [Google Scholar]
  5. Blitz, L., Fich, M., & Stark, A. A. 1982, ApJS, 49, 183 [NASA ADS] [CrossRef] [Google Scholar]
  6. Brand, J., & Blitz, L. 1993, A&A, 275, 67 [NASA ADS] [Google Scholar]
  7. Brunthaler, A., Reid, M. J., Menten, K. M., et al. 2011, Astron. Nachr., 332, 461 [NASA ADS] [CrossRef] [Google Scholar]
  8. Chen, X., Ellingsen, S. P., & Shen, Z.-Q. 2009, MNRAS, 396, 1603 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cragg, D. M., Johns, K. P., Godfrey, P. D., & Brown, R. D. 1992, MNRAS, 259, 203 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cyganowski, C. J., Whitney, B. A., Holden, E., et al. 2008, AJ, 136, 2391 [NASA ADS] [CrossRef] [Google Scholar]
  11. Du, F., & Yang, J. 2008, ApJ, 686, 384 [NASA ADS] [CrossRef] [Google Scholar]
  12. Faúndez, S., Bronfman, L., Garay, G., et al. 2004, A&A, 426, 97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Garden, R. P., Hayashi, M., Hasegawa, T., Gatley, I., & Kaifu, N. 1991, ApJ, 374, 540 [NASA ADS] [CrossRef] [Google Scholar]
  14. Gibb, A. G., Wyrowski, F., & Mundy, L. G. 2004, ApJ, 616, 301 [NASA ADS] [CrossRef] [Google Scholar]
  15. Girart, J. M., Beltrán, M. T., Zhang, Q., Rao, R., & Estalella, R. 2009, Science, 324, 1408 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  16. Gregory, P. C., & Condon, J. J. 1991, ApJS, 75, 1011 [NASA ADS] [CrossRef] [Google Scholar]
  17. Guan, X., Wu, Y., & Ju, B. 2008, MNRAS, 391, 869 [NASA ADS] [CrossRef] [Google Scholar]
  18. Gwinn, C. R., Moran, J. M., & Reid, M. J. 1992, ApJ, 393, 149 [NASA ADS] [CrossRef] [Google Scholar]
  19. He, J. H., Takahashi, S., & Chen, X. 2012, ApJS, 202, 1 [NASA ADS] [CrossRef] [Google Scholar]
  20. Helfand, D. J., Becker, R. H., White, R. L., Fallon, A., & Tuttle, S. 2006, AJ, 131, 2525 [NASA ADS] [CrossRef] [Google Scholar]
  21. Johnston, K. J., Gaume, R. A., Wilson, T. L., Nguyen, H. A., & Nedoluha, G. E. 1997, ApJ, 490, 758 [NASA ADS] [CrossRef] [Google Scholar]
  22. Keto, E. 2002, ApJ, 568, 754 [NASA ADS] [CrossRef] [Google Scholar]
  23. Klaassen, P. D., & Wilson, C. D. 2007, ApJ, 663, 1092 [NASA ADS] [CrossRef] [Google Scholar]
  24. Krumholz, M. R., Cunningham, A. J., Klein, R. I., & McKee, C. F. 2010, ApJ, 713, 1120 [NASA ADS] [CrossRef] [Google Scholar]
  25. Kuchar, T. A., & Bania, T. M. 1994, ApJ, 436, 117 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kurtz, S., Churchwell, E., & Wood, D. O. S. 1994, ApJS, 91, 659 [NASA ADS] [CrossRef] [Google Scholar]
  27. Larson, R. B. 1981, MNRAS, 194, 809 [NASA ADS] [CrossRef] [Google Scholar]
  28. Liu, H. B., Ho, P. T. P., & Zhang, Q. 2010a, ApJ, 725, 2190 [NASA ADS] [CrossRef] [Google Scholar]
  29. Liu, T., Wu, Y.-F., & Wang, K. 2010b, Res. Astron. Astrophys., 10, 67 [NASA ADS] [CrossRef] [Google Scholar]
  30. Longmore, S. N., Burton, M. G., Minier, V., & Walsh, A. J. 2006, MNRAS, 369, 1196 [NASA ADS] [CrossRef] [Google Scholar]
  31. Longmore, S. N., Pillai, T., Keto, E., Zhang, Q., & Qiu, K. 2011, ApJ, 726, 97 [NASA ADS] [CrossRef] [Google Scholar]
  32. Matthews, H., Kirk, H., Johnstone, D., et al. 2009, AJ, 138, 1380 [NASA ADS] [CrossRef] [Google Scholar]
  33. Menten, K. M., Reid, M. J., Pratap, P., Moran, J. M., & Wilson, T. L. 1992, ApJ, 401, L39 [NASA ADS] [CrossRef] [Google Scholar]
  34. Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Minier, V., Burton, M. G., Hill, T., et al. 2005, A&A, 429, 945 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Myers, P. C., Mardones, D., Tafalla, M., Williams, J. P., & Wilner, D. J. 1996, ApJ, 465, L133 [NASA ADS] [CrossRef] [Google Scholar]
  37. Olmi, L., & Cesaroni, R. 1999, A&A, 352, 266 [NASA ADS] [Google Scholar]
  38. Pandian, J. D., Momjian, E., & Goldsmith, P. F. 2008, A&A, 486, 191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Pascucci, I., Apai, D., Henning, T., Stecklum, B., & Brandl, B. 2004, A&A, 426, 523 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Peretto, N., & Fuller, G. A. 2009, A&A, 505, 405 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Pestalozzi, M. R., Minier, V., & Booth, R. S. 2005, A&A, 432, 737 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Phillips, C. J., Norris, R. P., Ellingsen, S. P., & McCulloch, P. M. 1998, MNRAS, 300, 1131 [NASA ADS] [CrossRef] [Google Scholar]
  43. Plambeck, R. L., & Menten, K. M. 1990, ApJ, 364, 555 [NASA ADS] [CrossRef] [Google Scholar]
  44. Purcell, C. R., Balasubramanyam, R., Burton, M. G., et al. 2006, MNRAS, 367, 553 [NASA ADS] [CrossRef] [Google Scholar]
  45. Ragan, S., Henning, T., Krause, O., et al. 2012, A&A, 547, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Rathborne, J. M., Jackson, J. M., & Simon, R. 2006, ApJ, 641, 389 [NASA ADS] [CrossRef] [Google Scholar]
  47. Reid, M. J., & Menten, K. 1991, in Atoms, Ions and Molecules: New Results in Spectral Line Astrophysics, eds. A. D. Haschick, & P. T. P. Ho, ASP Conf. Ser., 16, 375 [Google Scholar]
  48. Rodriguez, L. F., Torrelles, J. M., & Moran, J. M. 1981, AJ, 86, 1245 [Google Scholar]
  49. Roman-Duval, J., Jackson, J. M., Heyer, M., et al. 2009, ApJ, 699, 1153 [Google Scholar]
  50. Russeil, D. 2003, A&A, 397, 133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Russeil, D., Adami, C., & Georgelin, Y. M. 2007, A&A, 470, 161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Sato, M., Hirota, T., Reid, M. J., et al. 2010, PASJ, 62, 287 [NASA ADS] [Google Scholar]
  53. Schwartz, P. R., Gee, G., & Huang, Y.-L. 1988, ApJ, 327, 350 [NASA ADS] [CrossRef] [Google Scholar]
  54. Simon, R., Rathborne, J. M., Shah, R. Y., Jackson, J. M., & Chambers, E. T. 2006, ApJ, 653, 1325 [NASA ADS] [CrossRef] [Google Scholar]
  55. Sollins, P. K., & Ho, P. T. P. 2005, ApJ, 630, 987 [NASA ADS] [CrossRef] [Google Scholar]
  56. Sutton, E. C., Sobolev, A. M., Ellingsen, S. P., et al. 2001, ApJ, 554, 173 [NASA ADS] [CrossRef] [Google Scholar]
  57. Urquhart, J. S., Moore, T. J. T., Hoare, M. G., et al. 2011, MNRAS, 410, 1237 [NASA ADS] [CrossRef] [Google Scholar]
  58. Val’tts, I. E., & Larionov, G. M. 2007, Astron. Rep., 51, 519 [NASA ADS] [CrossRef] [Google Scholar]
  59. Voronkov, M. A., Sobolev, A. M., Ellingsen, S. P., & Ostrovskii, A. B. 2005, MNRAS, 362, 995 [NASA ADS] [CrossRef] [Google Scholar]
  60. Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 1998, MNRAS, 301, 640 [NASA ADS] [CrossRef] [Google Scholar]
  61. White, R. L., Becker, R. H., & Helfand, D. J. 2005, AJ, 130, 586 [NASA ADS] [CrossRef] [Google Scholar]
  62. Wood, D. O. S., & Churchwell, E. 1989, ApJ, 340, 265 [NASA ADS] [CrossRef] [Google Scholar]
  63. Wu, Y., Zhang, Q., Chen, H., et al. 2005, AJ, 129, 330 [NASA ADS] [CrossRef] [Google Scholar]
  64. Wu, Y. W., Xu, Y., Pandian, J. D., et al. 2010, ApJ, 720, 392 [NASA ADS] [CrossRef] [Google Scholar]
  65. Xu, Y., Moscadelli, L., Reid, M. J., et al. 2011, ApJ, 733, 25 [NASA ADS] [CrossRef] [Google Scholar]
  66. Yang, J., Umemoto, T., Iwata, T., & Fukui, Y. 1991, ApJ, 373, 137 [NASA ADS] [CrossRef] [Google Scholar]
  67. Zavagno, A., Deharveng, L., Nadeau, D., & Caplan, J. 2002, A&A, 394, 225 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Zhou, S., Evans, II, N. J., Koempe, C., & Walmsley, C. M. 1993, ApJ, 404, 232 [NASA ADS] [CrossRef] [Google Scholar]
  69. Zuo, Y.-X., Yang, J., Shi, S.-C., et al. 2004, Chin. J. Astron. Astrophys., 4, 390 [Google Scholar]
  70. Zuo, Y.-X., Li, Y., Sun, J., & Yang, J. 2011, Chin. Astron. Astrophys., 35, 439 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A

Table A.1

12CO and 13CO line properties.

Table A.2

Physical parameters and radio continuum emission of the identified cores.

thumbnail Fig. A.1

13CO (1 − 0) maps (contours) overlaid on the far-infrared images (gray). The thick contour in each map denote the 50 per cent level of the maximum intensity. For each map, the interpolation is performed upon a grid with 100 × 100 pixels. For all the regions, it is guaranteed that the pixel size is far smaller than the mapping step. And in each panel, the black cross indicates the central position (where the maser locates), the square indicates the radio continuum source (plausible Ultra- or Hyper-Compact Hii regions, if they exist), the dashed ellipses indicate the identified cores if there are overlapped multiple cores.

thumbnail Fig. A.1

continued.

thumbnail Fig. A.1

continued.

thumbnail Fig. A.1

continued.

thumbnail Fig. A.2

CO spectra (12CO\hbox{${\rm ^{12}CO}$} in black and 13CO\hbox{${\rm ^{13}CO}$} in red) at the mapping point closest to the core center and position–velocity diagrams of the outflow candidates. The dashed lines indicate the velocity ranges adopted as the line wings for the four sources (G59.78 only has evident red wings). The swelled line wings that may indicate the outflow are also encircled with the dashed ellipses.

Table A.3

Infall velocities and mass infall rates.

thumbnail Fig. A.3

Outflow emissions. The solid blue and dashed red contours indicate the blue- and red-shifted emissions, respectively. The dashed line in each panel indicates the speculated outflow orientation. The velocity ranges for the integration are shown in each panel. The black cross indicates the central position (maser source).

thumbnail Fig. A.4

CO spectra of the infall candidates. The spectra of another infall candidate.

thumbnail Fig. A.5

Same as Fig. 1 but with the far distance adopted for all the sources with two Dkin solutions. The statistical box is the same as described in Fig. 3.

All Tables

Table 1

Physical parameters of the maser sources.

Table A.1

12CO and 13CO line properties.

Table A.2

Physical parameters and radio continuum emission of the identified cores.

Table A.3

Infall velocities and mass infall rates.

All Figures

thumbnail Fig. 1

Number of sources in each 10-degree Galactic longitude interval. The gray and black-line columns represent the full sample of the methanol masers (Pestalozzi et al. 2005) and the fraction covered by our sample, respectively. Overlaid inset panel: distributions of the maser intensities for the sampled and missed sources in l = 10° to 50° range. For each statistical box, the entire error bar (whisker) represents the overall distribution, the box and its inner segment represent the 25% to 75% range and the median value, respectively.

In the text
thumbnail Fig. 2

Distribution of our observed sources on the Milky Way disk. The source sizes are proportional to their 13CO\hbox{${\rm ^{13}CO}$} line intensity. The thick faint blue lines represent the distribution of the spiral arms (Russeil 2003) in our observed field. The background image is the artist’s view of the Galactic morphology (R. Hurt: NASA/JPL-Caltech/SSC, with a more detailed introduction in Urquhart et al. 2011).

In the text
thumbnail Fig. 3

Distribution of the physical parameters with the heliocentric distance. Left panels: all the data points. Right panels: statistical box plots for the data distribution in each 2-kpc interval. The dots represent the outliers, which are apparently isolated from the major fraction of the data points. For each statistical box, the entire whisker represents the overall distribution. If there are isolated data points (based on eye selection), the whisker then represents the remaining data range. The box and its inner segment represent the 25% to 75% range and the median value, respectively. We selected interval pairs, which may typically reflect the variation trend. The interval pairs are indicated with gray segments above the whiskers. Labeled on the segment are the p values from the t test and M-W test respectively, and the test result (“yes” or “no” in the brackets) for whether the two intervals have significantly different distributions.

In the text
thumbnail Fig. 4

Distribution of the physical parameters with the Galactocentric distance. The first panel shows the data points for the 12CO\hbox{${\rm ^{12}CO}$} excitation temperature. The other panels show the statistics of the parameters in each 2-kpc interval. For each statistical box, the entire whisker represents the overall distribution. If there are isolated data points (based on eye selection), the whisker then represents the remaining data range. The box and its inner segment represent the 25% to 75% range and the median value, respectively. We selected interval pairs which may typically reflect the variation trend. The interval pairs are indicated with gray segments above the whiskers. Labeled on the segment are the p values from the t test and M-W test respectively, and the test result (“yes” or “no” in the brackets) for whether the two intervals have significantly different distributions.

In the text
thumbnail Fig. 5

Variation of the 12CO\hbox{${\rm ^{12}CO}$} excitation temperature and the 13CO\hbox{${\rm ^{13}CO}$} line width as a function the H2 column density. The dashed lines represent the best-fit linear correlation.

In the text
thumbnail Fig. 6

Statistical comparison for the excitation temperatures between the different subsamples of the mapped cores. Left panel: the radio-bright and radio-faint cores. Right panel: the IR-bright and IR-faint cores. For each data set, the whisker represents the overall distribution; the box and inner segment represent the 25% to 75% range and the median value, respectively.

In the text
thumbnail Fig. 7

Comparison of the physical parameters in our maser sources and the IRDCs (Du & Yang 2008). The left panels show the parameter distribution in each interval. The right panels show the overall statistics. For each data set, the whiskers represent the overall distribution; and the box and its inner segment represent the 25% to 75% range and the median value, respectively.

In the text
thumbnail Fig. A.1

13CO (1 − 0) maps (contours) overlaid on the far-infrared images (gray). The thick contour in each map denote the 50 per cent level of the maximum intensity. For each map, the interpolation is performed upon a grid with 100 × 100 pixels. For all the regions, it is guaranteed that the pixel size is far smaller than the mapping step. And in each panel, the black cross indicates the central position (where the maser locates), the square indicates the radio continuum source (plausible Ultra- or Hyper-Compact Hii regions, if they exist), the dashed ellipses indicate the identified cores if there are overlapped multiple cores.

In the text
thumbnail Fig. A.1

continued.

In the text
thumbnail Fig. A.1

continued.

In the text
thumbnail Fig. A.1

continued.

In the text
thumbnail Fig. A.2

CO spectra (12CO\hbox{${\rm ^{12}CO}$} in black and 13CO\hbox{${\rm ^{13}CO}$} in red) at the mapping point closest to the core center and position–velocity diagrams of the outflow candidates. The dashed lines indicate the velocity ranges adopted as the line wings for the four sources (G59.78 only has evident red wings). The swelled line wings that may indicate the outflow are also encircled with the dashed ellipses.

In the text
thumbnail Fig. A.3

Outflow emissions. The solid blue and dashed red contours indicate the blue- and red-shifted emissions, respectively. The dashed line in each panel indicates the speculated outflow orientation. The velocity ranges for the integration are shown in each panel. The black cross indicates the central position (maser source).

In the text
thumbnail Fig. A.4

CO spectra of the infall candidates. The spectra of another infall candidate.

In the text
thumbnail Fig. A.5

Same as Fig. 1 but with the far distance adopted for all the sources with two Dkin solutions. The statistical box is the same as described in Fig. 3.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.