This site uses cookies. By continuing to use this site you agree to our use of cookies. To find out more, see our Privacy and Cookies policy.

AN OUTER ARM IN THE SECOND GALACTIC QUADRANT: STRUCTURE

, , , , , , and

Published 2016 May 13 © 2016. The American Astronomical Society. All rights reserved.
, , Citation Xinyu Du et al 2016 ApJS 224 7 DOI 10.3847/0067-0049/224/1/7

0067-0049/224/1/7

ABSTRACT

The lack of arm tracers, especially remote tracers, is one of the most difficult problems preventing us from studying the structure of the Milky Way. Fortunately, with its high-sensitivity CO survey, the Milky Way Imaging Scroll Painting (MWISP) project offers such an opportunity. Since completing about one-third of its mission, an area of l = [100, 150]°, b = [−3, 5]° has nearly been covered. The Outer arm of the Milky Way first clearly revealed its shape in the second galactic quadrant in the form of molecular gas—this is the first time that the Outer arm has been reported in such a large-scale mapping of molecular gas. Using the 115 GHz 12CO(1–0) data of MWISP at the LSR velocity ≃[−100, −60] km s−1 and in the area mentioned above, we have detected 481 molecular clouds in total, and among them 332 (about 69%) are newly detected and 457 probably belong to the Outer arm. The total mass of the detected Outer arm clouds is ∼3.1 × 106 M. Assuming that the spiral arm is a logarithmic spiral, the pitch angle is fitted as ∼13fdg1. Besides combining both the CO data from MWISP and the 21 cm H i data from the Canadian Galactic Plane Survey (CGPS), the gas distribution, warp, and thickness of the Outer arm are also studied.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

It is well-known that our galaxy has a spiral structure, but (Oort et al. 1958; Bok 1959) a detailed understanding of the Milky Way is still lacking. Since Georgelin & Georgelin (1976) presented a large-scale spiral pattern with H ii regions, nearly 100 Milky Way models have been proposed (Steiman-Cameron 2010). However, most models share the same large structure in spite of differing details (e.g., Russeil et al. 2007; Vallée 2008; Hou & Han 2014).

The fact that there exists a spiral arm beyond the Perseus arm has been generally recognized for a long time. Since Moffat & Vogt (1975) noticed that several young star clusters were located far beyond the position of the Perseus arm, more work has been done to confirm the existence of this external arm. Henderson et al. (1982) first systematically studied the distribution of H i in the outer galaxy beyond the solar circle, which is a vast frontier of our galaxy. Heyer et al. (1998) described the molecular image of the outer galaxy with the Five College Radio Astronomy Observatory (FCRAO) CO survey results, providing evidence for the presence of molecular clouds in this outer area. Negueruela & Marco (2003) traced the external arm using OB stars and concluded that the Cam OB3 association lies on it. On the basis of a maser source (associated with high-mass star-forming region, HMSFR) distances, Reid et al. (2014) first delineated several spiral arms, including the Outer arm, using the trigonometric parallax method. In addition, detailed parallax results and analyses of the Outer arm are also presented by Hachisuka et al. (2015).

However, accurate maser source positions cannot outline the distributions of interstellar gases. The FCRAO CO survey only traced the high-mass molecular clouds because of its poor sensitivity. Fortunately, the Milky Way Imaging Scroll Painting (MWISP) project4 provides a chance to study the molecular gases with no-bias high sensitive 12CO(1 – 0), 13CO(1 − 0), and C18O(1 − 0) observations. Due to such high-sensitivity observations, a new spiral arm (hereafter the New arm) beyond the Outer arm has been discovered by Sun et al. (2015). Combining the 12CO(1 − 0) observations of MWISP and atomic hydrogen data from the Canadian Galactic Plane Survey (CGPS; Taylor et al. 2003), we presented the results of the Outer arm: the arm located between the Perseus arm and the New arm. The results derived from 13CO(1 − 0) and C18O(1 − 0) data and further analyses will be published in future works.

Please note that the Outer arm does not have a particular name yet. Some authors have labeled it the "Cygnus arm," or the "Perseus +I arm," or the "Norma—Cygnus arm" (Vallée 2008). Otherwise, the name "Outer arm" is also widely used (e.g., Dame et al. 2001; Reid et al. 2014). Considering the fact that the name "Cygnus" has been described as the location near the Sun in some early papers, here we choose the name "Outer arm" to avoid confusion.

In Section 2 we introduce our CO observation conditions and archival data of atomic hydrogen. In Section 3 we describe how we pick out the Outer arm clouds and briefly study their properties. In Section 4 we study the properties of the Outer arm, including the pitch angle, the gas distribution, the thickness, and the warp. A10 summary is given in Section 5.

2. OBSERVATIONS AND ARCHIVAL DATA

2.1. CO Observations

The 12CO (1 − 0), 13CO(1 − 0), and C18O(1 − 0) lines were observed simultaneously using the Purple Mountain Observatory Delingha 13.7 m telescope from 2011 September to 2015 March as one of the scientific demonstration regions for the MWISP project, which is the first no-bias high-sensitivity CO survey with such a large-scale aiming at l = [−10, 250]°, b = [−5, 5]°. With the on-the-fly (OTF) observing mode, MWISP now has completed about one-third of its plan, and an area of l = [100, 150]°, b = [−3, 5]° has mostly been covered—the total area was 288 square degrees. A superconductor-insulator-superconductor (SIS) superconducting receiver with a nine-beam array was used as the front end (Shan et al. 2012). A Fast Fourier Transform (FFT) spectrometer with a total bandwidth of 1000 MHz and 16,384 channels was used as the back end. For the 115 GHz 12CO (1 − 0) observations, the main beam width was about 52'', the main beam efficiency (ηMB) was 0.46, and the typical rms noise level was ∼0.5 K, corresponding to a channel width of 0.16 km s−1. (∼0.2 K per 0.8 km s−1 channel, which is 3 times better than the FCRAO OGS sensitivity.) All the data were corrected by ${T}_{\mathrm{MB}}={T}_{{\rm{A}}}^{*}/{\eta }_{\mathrm{MB}}$. The data were sampled every 30''. All the data were reduced using the GILDAS/CLASS package.

2.2. Archival Data of Atomic Hydrogen

The 21 cm line data were retrieved from the CGPS. We downloaded data of l = [63, 155]°, b = [−3, 5]° from the Canadian Astronomy Data Centre.5 The velocity coverage of the data is in the range of −153 to 40 km s−1, with a channel separation of 0.82 km s−1. The survey has a spatial resolution of $58^{\prime\prime} $, which is comparable to our CO observations.

3. ANALYSIS OF CLOUDS

We have detected 481 clouds in total; 332 (about 69%) clouds are newly detected, 457 clouds are identified in the Outer arm, and 24 are identified in the New arm. Among the Outer arm clouds, 7 are reported by Brand & Wouterloot (1994; hereafter BW94), 75 are reported by Heyer et al. (2001; hereafter HCS01 clouds), 125 are reported by Brunt et al. (2003; hereafter BKP03). And all 24 of the New arm clouds are newly detected. (These 24 New arm clouds do not overlap with the clouds detected by Sun et al. 2015.) The parameters of all 481 clouds are summarized in Table 1.

Table 1.  Parameters of Molecular Clouds Derived by CO(1–0)

Name VLSR Tpeak ΔV ICO Area d R Size Mass Z Scale Note
  (km s−1) (K) (km s−1) (K km s−1) (arcmin2) (kpc) (kpc) (pc) (103 M) (kpc)  
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)
MWISP G101.400+03.967 −86.0 2.3 1.6 3.9 15.4 8.2 12.8 10.4 1.3 0.6  
MWISP G101.508+04.092 −90.5 3.2 1.1 3.6 3.0 8.6 13.1 4.4 0.2 0.6  
MWISP G101.633+02.867 −82.6 3.6 1.6 6.2 10.0 7.8 12.5 7.9 1.2 0.4  
MWISP G101.683+02.917 −81.7 4.4 1.3 5.9 21.7 7.7 12.4 11.6 2.5 0.4  
MWISP G101.700+03.833 −89.9 7.9 1.4 11.5 25.8 8.6 13.1 14.2 7.4 0.6  
MWISP G101.708+02.867 −84.3 4.9 1.2 6.3 18.4 8.0 12.6 11.1 2.5 0.4  
MWISP G101.767+02.808 −82.6 16.5 2.9 51.0 43.3 7.8 12.5 16.7 45.8 0.4  
MWISP G101.800+02.958 −81.4 2.5 0.7 2.0 4.8 7.7 12.4 5.2 0.2 0.4  
MWISP G101.808+03.817 −88.0 2.2 1.6 3.9 18.5 8.3 12.9 11.5 1.7 0.6  
MWISP G101.850+02.908 −81.7 7.7 2.0 16.1 40.4 7.7 12.4 15.9 13.1 0.4  

Note. Column (1): source named by the MWISP project and the Tpeak position in Galactic Coordinates. Columns (2)−(5): results of a Gaussian fit to the spectra. Column (6): solid angle defined by the 3σ limits. Columns (7)−(8): heliocentric distance and galactocentric radius; both are derived by galactic parameters obtained by Reid et al. (2014). Column (9): equivalent diameters of the molecular clouds corrected by the beam size of the telescope. Column (10): cloud mass calculated by X = 1.8 × 1020 cm−2 (K km s−1)−1 (Dame et al. 2001). Column (11): scale height (=Dsin(b)). Column (12): BW94: sources detected by Brand & Wouterloot (1994); HCS01: sources detected by Heyer et al. (2001); BKP03: sources detected by Brunt et al. (2003); new Arm: clouds that are supposed to belong to the new spiral arm detected by Sun et al. (2015).

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

3.1. Cloud Identification

Because of the differential rotation of the Milky Way, most LSR velocities that are consistent with circular Galactic rotation are negative in the second quadrant. Starting from VLSR ∼ 0 km s−1, increasingly negative velocities successively trace the Local arm, the Perseus arm, the Outer arm, and the New arm. In order to find all of the Outer arm clouds, we need to know the Outer arm LSR velocity range at every galactic longitude. As mentioned above, the Outer arm LSR velocity is located between the Perseus arm and the New arm, which provides us with a clue for picking out the clouds. First, using the H i data from the CGPS, we plotted the longitude-velocity map of H i, integrated over all latitudes (from −3° to 5°). Second, we projected the spatial (x, y) curves of the Outer arm, the Perseus arm (both fitted by Reid et al. 2014, and hereafter the Reid Outer spiral and the Reid Perseus spiral), and the New arm (fitted by Sun et al. 2015, and hereafter the Sun New spiral) on that H i longitude-velocity map. (In other words, we converted the dashed curves shown in Figure 3 into the ones shown in Figure 5. And the converting method is presented in detail in Section 4.3.) And then we marked all the HCS01 clouds on that map. Combining the H i map, the positions of the longitude-velocity curves, and HCS01 clouds, we can estimate the velocity range of the Outer arm as a function of galactic longitude. Third, we compiled an automatic procedure to list all the positions with emissions greater than 2.5σ of the data cube at that velocity range (where σ is the typical rms noise and =0.5 K). Last, we checked both the list and data cube and identified the clouds via the naked eye. It is necessary to emphasize that we did not separate the isolated clouds or cloud complexes into small pieces of molecular clumps. In other words, some of the clouds have multiple spectral or spatial peaks.

Finally, we have detected 481 clouds in total. But not all of them belong to the Outer arm. Of these clouds, 24 may be located in the New arm, since the velocity gap between the 24 clouds and the other ones is relatively large in the zoomed-in longitude-velocity map.

However, it is necessary to point out that the Outer arm clouds that we we found do not absolutely belong to that arm. From $l\simeq 100^\circ $ to 120°, there exists an arm-blending region. The LSR velocities of H i gas of the Perseus arm and the Outer arm are mixed together (see Figures 5 and 6). Also, the LSR velocity gap of CO between the two arms is not obvious. Since our cloud identification criterion is based on the LSR velocity, it is difficult to pick out Outer arm clouds in that region. We may omit some clouds with an LSR velocity $\gtrsim -70$ km s−1, or falsely pick the clouds that may be located at the inter-arm area.

3.2. Cloud Parameters

Heliocentric distance is one of the most important values for studying both the cloud properties and the spiral structure. Generally there are three methods for measuring the distance: trigonometric parallax, photometry, and the kinematic method. And the accuracy of the trigonometric parallax and photometry methods is better than that for the kinematic method. (see Xu et al. 2006) However, there is only one source with a parallax distance and a luminosity distance in the Outer arm region: the one associated with MWISP G135.267+02.800—its parallax distance is 6.0 kpc (Reid et al. 2014) and its luminosity distance is also 6.0 kpc (Hou & Han 2014). So the kinematic method is the only choice for us.

We chose to use the kinematic distances on the basis of the galactic parameters of Model A5 of Reid et al. (2014) (Θ0 = 240 km s−1, R0 = 8.34 kpc, $\tfrac{d{\rm{\Theta }}}{{dR}}=-0.2$ km s−1 kpc−1, U = 10.7 km s−1, V = 15.6 km s−1, W = 8.9 km s−1, ${\bar{U}}_{s}=2.9$ km s−1, ${\bar{V}}_{s}=-1.6$ km s−1, hereafter Reid model) and the FORTRAN source code provided by Reid et al. (2009). Note that all of the cloud distances we finally used are kinematic distances, including MWISP G135.267+02.800.

The Reid model is much better than the IAU model (namely Θ0 = 220 km s−1, R0 = 8.5 kpc). Take the cloud MWISP G135.267+02.800 as an example. Its parallax and luminosity distances both are 6.0 kpc and the Reid kinematic distance (using Reid model) is 6.7 kpc, whereas the IAU kinematic distance (using IAU model) is 8.6 kpc. A detailed comparison can be seen in Section 4 of Reid et al. (2009)

However, one should also keep in mind that there still exist biases from the assumption of circular motion inherent in this kinematic distance model. In the second galactic quadrant the kinematic distance may be a little larger than the real distance, just as seen in the example of MWISP G135.267+02.800 or the left panel of Figure 3. A detailed discussion about the biases is presented in Section 4.2.

Knowing the distance d and latitude b, we calculated the scale height Z by $Z=d\mathrm{sin}(b)$. And using b, Z, d and longitude l, we calculated the galactocentric radius R with ${R}^{2}={d}^{2}{\mathrm{cos}}^{2}(b)+{R}_{0}^{2}-2{R}_{0}d\mathrm{cos}(b)\mathrm{cos}(l)+{Z}^{2}$.

The cloud solid angle A is defined by the 3σ limits. The cloud diameter D is obtained after the beam deconvolution: $D=d\sqrt{\tfrac{4}{\pi }A-{\theta }_{\mathrm{MB}}^{2}}$ (Ladd et al. 1994), where θMB is the main beam width. Adopting the CO-to-${{\rm{H}}}_{2}\;{\rm{X}}$ factor $1.8\times {10}^{20}\;{\mathrm{cm}}^{-2}{({\rm{K}}\cdot \mathrm{km}{{\rm{s}}}^{-1})}^{-1}$ (Dame et al. 2001), cloud mass is calculated from $M=2\;\mu \;{m}_{{\rm{H}}}\;X\;\pi {\left(\tfrac{D}{2}\right)}^{2}\int {T}_{{\rm{B}}}d\;V$, where μ = 1.36 (Hildebrand 1983) is the mean atomic weight per H atom in the ISM, mH is the H atomic mass, and $\int {T}_{{\rm{B}}}d\;V$ is the integrated intensity of the ${T}_{\mathrm{peak}}$ spectrum. However, the mass is probably underestimated since the X factor adopted is measured in the solar neighborhood, and a recent study has suggested an increase of X in the outer Galaxy. (Abdo et al. 2010).

3.3. Comparison with HCS01

Heyer et al. (2001) have used the CO data of the FCRAO Outer Galaxy Survey to identify molecular clouds. To make a comparison, we plotted the distributions of line width, size, and mass for the HCS01 Outer arm clouds and ours. We have two reasons for not comparing the clouds with BW94 and BKP03: (i) the number of BW94 clouds is too few and (ii) the BKP03 clouds are generated by the same data but from a different method than the HCS01 clouds, and their method tends to find small-size molecular clumps, which is not consistent with our cloud identification criterion. Figure 1 shows the results. One may notice that there are 102 HCS01 clouds in Figure 1, but only 75 clouds are labeled in Table 1. The reason is the different cloud identification methods: HCS01 includes both small, isolated clouds and clumps within larger cloud complexes, but our catalog only includes isolated clouds and cloud complexes. So the clouds are not exactly matched one-to-one.

Figure 1.

Figure 1. Distribution of the Outer arm cloud line width (a), size (b), and mass (c). The red and gray shades indicate the results from MWISP and FCRAO, respectively.

Standard image High-resolution image

Since the cloud heliocentric distances of HCS01 are derived from the IAU model, the sizes and masses are correspondingly biased. In order make our comparison using the same criterion, we revised the distances of the Reid model and the sizes and masses are correspondingly revised.

However, the comparison result is against expectations. The signal-to-noise ratio of our data is higher than that of HCS01. Consequently, HCS01 should have contained more luminous (and therefore more massive) clouds. But the distribution shows that the fractions of small sizes and low masses of HCS01 are higher than those of ours. Two main reasons may have caused this. (i) HCS01 includes smaller clumps within larger cloud complexes, but we did not detach clumps from cloud complexes. (ii) Because of our higher signal-to-noise ratio, we can detect a larger angle area for the same cloud, which results in larger sizes and masses.

4. PROPERTIES OF OUTER ARM

4.1. Pitch Angle

Spiral arms of external galaxies are usually approximated by logarithmic spirals. This is a simple but reasonable assumption and is also applicable to our Milky Way galaxy. That the pitch angle of one spiral is constant is one of the properties of a logarithmic spiral. Usually, the spiral arms of one galaxy are roughly fitted by logarithmic spirals with the same pitch (e.g., Vallée 2008), or more accurately, different spiral arms are fitted by different pitches (e.g., Russeil 2003). And sometimes even one spiral arm is fitted by varying the pitch angle (e.g., the Sagittarius arm in Taylor & Cordes 1993). In fact, Honig & Reid (2015) recently showed that in some external galaxies, the pitch angle in the same spiral arm varies, and the variance is as large as the that among different arms within the same galaxy. So a more accurate assumption may be that the pitch angle of one segment of the arm within one galaxy is constant.

Now we assume that the segment form l = 100° to l = 150° of the Outer arm spiral is a logarithmic spiral, and using the same equation as Reid et al. (2009, 2014), we fitted the spiral of those 457 Outer arm clouds using

Equation (1)

where ψ is the spiral pitch angle, β is the galactocentric azimuth (namely the Source-GC-Sun angle), and Rref and βref are the reference radius and reference azimuth, respectively. The fitting algorithm is a minimizing chi-square error statistic, and the chi-square error statistic is computed as

Equation (2)

where $k=-\mathrm{tan}(\psi )$, $b={\beta }_{\mathrm{ref}}\mathrm{tan}(\psi )+\mathrm{ln}({R}_{\mathrm{ref}})$, xi and ${y}_{i}$ are β and $\mathrm{ln}(R)$ of each cloud, respectively, and Wi is the weight. The weight is defined as $W={\mathrm{log}}_{10}\left(\tfrac{M}{{M}_{\odot }}\right)$.

Figure 2 shows the fitting result. The pitch angle is $\sim 13\fdg 1$, the Rref is ∼13.6 kpc, the βref is ∼26fdg9, and the χ2 is ∼4.1. Our fitting pitch angle is close to the result of Reid et al. (2014; pitch = 13fdg8), and is consistent with the results of Vallée (2015), who summarized large numbers of recent studies about Milky Way pitch angle and yielded a mean global value of 13fdg1.

Figure 2.

Figure 2. Fitting result of the Outer arm pitch angle. R is the Galactocentric radius in units of kpc. β is Galactocentric longitude. Circles indicate molecular clouds of the Outer arm, and their sizes indicate different masses. The fitting line is weighted by mass.

Standard image High-resolution image

4.2. Plan View

Since Georgelin & Georgelin (1976) presented the famous Milky Way plane view, numerous models have been suggested (see Figure 2 in Steiman-Cameron 2010). Our knowledge about the Milky Way structure has advanced since then a little, so at least the existence of Outer arm is now undoubted, though it was absent in the Georgelin model. In contrast, the lack of arm tracers, especially the remote tracers, still prevent us from knowing more about our galaxy. Thanks to the high-sensitivity CO survey of MWISP, so many molecular clouds of the Outer arm have been detected, and this is the first time that we have detected such large numbers of remote molecular clouds, which contributes a lot to studying the Milky Way structure.

The left panel of Figure 3 shows the plane view of the Milky Way. The red thick curve indicates the Outer arm spiral fitted by us (hereafter our Outer spiral). Since the Outer arm pitch angle fitted by us is similar to the one fitted by Reid et al. (2014), and their parallax distance of the Outer arm is more precise, we moved our Outer spiral to be parallel to the position of the Reid Outer spiral. Accordingly, the Outer arm clouds are parallelly moved by the same distance. The detailed moving process is as follows. First, using the Rref and βref fitted by Reid et al. (2014), with the pitch angle fitted by us, we plotted a new spiral curve as the parallel moved curve. Second, we calculated the distance between the new curve and the old curve at every cloud galactic longitude. Third, the new cloud heliocentric distance was calculated as the "old cloud heliocentric distance minus the distance obtained in the second step." This moving process decreased the cloud heliocentric distances but did not change their longitudes. The right panel of Figure 3 is the result after the move.

Figure 3.

Figure 3. (Left panel) Locations of molecular clouds and HMSFRs in the plane view of the Milky Way. The red and blue circles respectively mark the molecular clouds of the Outer arm and New arm, which are summarized in Table 1. The blue triangles mark the molecular clouds of the New arm that were detected by Sun et al. (2015). Different sizes of circles and triangles indicate different masses. The red and green squares respectively mark the HMSFRs (Reid et al. 2014) of the Outer and Perseus arms. Distance error bars of HMSFRs are indicated. The color depths of the red circles indicate their Z heights. The red diamonds indicate the locations of the Outer arm HMSFRs calculated by kinematic method. The green, red, and blue dashed curves respectively indicate the Reid Perseus spiral, the Reid Outer spiral, and the Sun New spiral. The red thick curve indicates our Outer spiral. The gray dotted curves indicate galactocentric radii =10, 12, 14, 16 kpc. The gray dotted lines indicate galactic longitudes =60°, 75°, 90°, 105°, 120°, 135°, 150°. The gray crosses and the words beside them indicate the locations and the corresponding LSR velocities. (Right panel) Translational result of the Outer arm clouds and our Outer spiral. For a better comparison, we have not moved or changed any symbols except the red circles and the red thick curve, and the sizes of the red circles remain the same although their masses are changed because of the revised distances.

Standard image High-resolution image

We mentioned these distance biases in Section 3.2. One may wonder why the kinematic distances obtained from the model that Reid et al. (2014) provided are greater than the distances that they measured. In fact, the biases are mainly caused by errors. Limited parallax data and the parallax measuring errors lead to larger errors of the galactic model. The Reid model narrows the gap between the kinematic distance and the parallax distance but cannot perfectly match them. (And actually, the cloud distances that we moved are roughly below the parallax distance errors. Just see the error bars of the Outer arm HMSFRs in Figure 3.) Besides, the Reid model is a global model and it is fitted by maser sources that distribute in almost three galactic quadrants (see Figure 1 of Reid et al. 2014). But different regions of the Milky Way may have their own peculiar motions. So in different regions the kinematic distances calculated from their model may be a little biased. As a result of these reasons, in the second galactic quadrant (or more exactly in the Outer arm region of the second galactic quadrant) the kinematic distances are a little larger than the parallax distances, just as we see in the left panel of Figure 3.

Figure 4 shows the whole view of the Milky Way. We extended our translational (namely the parallel moved) spiral to the inner galaxy. The Reid Outer spiral is also plotted. Additionally, for comparison we plotted two recent fitting results for this arm, of which the arm tracers or fitting methods are different. The cyan dashed curve is one of the fitting results from Hou & Han (2014; "arm-5" in the third column of Table 4 in their paper). Their arm tracers are H ii regions, and their fitting model is the polynomial-logarithmic spiral arm model. The magenta pecked curve is the fitting result from Bobylev & Bajkova (2014). Their arm tracers are 3 HMSFRs and 12 very young star clusters, and their fitting model is the the logarithmic spiral. The fact that our translational spiral is located closer to the other three spirals may suggest that the translational locations of the clouds are relatively better. However, this suggestion is not adequate enough to make us revise the distances and other parameters such as mass in Table 1. The original parameters derived in Section 3.2 are retained.

Figure 4.

Figure 4. An artist's conception of the Milky Way (R. Hurt: NASA/JPL-Caltech/SSC). The red circles indicate the parallel moved Outer arm clouds with mass >103 M. The white squares indicate the Outer arm HMSFRs and the distance error bars are also plotted. The red thick curve indicates the extended translational position of our Outer spiral, and the red thick dashed curve indicates the original position of our Outer spiral. The white dotted curve indicates the Reid Outer spiral. The cyan dashed curve and magenta pecked curve indicate the Outer arm fitting results of Hou & Han (2014) and Bobylev & Bajkova (2014), respectively.

Standard image High-resolution image

4.3. l–V, V–b and l–b Map

Figure 5 shows the l–V map of H i emission. All the clouds (including the New arm clouds) and HMSFRs are marked on it. Also, the projections of the Reid Outer spiral, the Reid Perseus spiral, the Sun New spiral, and our Outer spiral are plotted. In order to project those spiral curves on the l–V map, or in other words, to convert the dashed curves shown in Figure 3 into the ones shown in Figure 5, we compiled a C program on the basis of the FORTRAN code provided by Reid et al. (2009). Their program can calculate the revised heliocentric kinematic distances when inputting the cloud positions and the LSR velocities, just as mentioned in Section 3.2. Our program is the inverse—the inputs are cloud position (galactic longitude and latitude) and galactocentric radius, and the output is the LSR velocity. Knowing the expression of the arm curves (namely through Equation (1)), we can obtain an array of galactocentric radii at every galactic longitude. And then ,using the Reid model, the cooresponding LSR velocities can be calculated by our C program. Using the arrays of longitudes and velocities, we then can plot the l–V curves of the arms.

Figure 5.

Figure 5. Longitude–velocity diagram of H i from the CGPS, integrated over all latitudes (from −3° to 5°). The red and blue circles respectively mark the molecular clouds of the Outer arm and the New arm, which are summarized in Table 1. The blue triangles mark the molecular clouds of the New arm detected by Sun et al. (2015). Different sizes of circles and triangles indicate different masses. The red and green squares respectively mark the HMSFRs of the Outer and Perseus arms. The green, red, and blue dashed curves respectively indicate the projections of the Reid Perseus spiral, the Reid Outer spiral, and the Sun New spiral. The red thick curve indicates the projection of the Outer arm spiral fitted by us. The pink ribbon indicates the integrated velocity window of Figure 7.

Standard image High-resolution image

Figure 6 shows the V–b montage of H i emission. The high-mass Outer arm clouds (>103 M), the New arm clouds, and the HMSFRs are marked on it.

Figure 6.

Figure 6. Velocity–latitude montage of H i from the CGPS, integrated every 5° of galactic longitude. The Outer arm is embraced in the two dotted red lines in each figure except the figures of the arm-blending region. The meanings of other symbols are consistent with Figure 5. Note that the low-mass clouds (<103 M) are not plotted, and the symbol size does NOT indicate mass.

Standard image High-resolution image

As mentioned in Section 3.1, from longitude l ≃ 100° to 120°, there exists an arm-blending region. This phenomenon can be clearly seen from those two figures and may mainly be caused for the following three reasons. (i) Streaming motions near spiral arms have long been predicted by density wave theory and have been observed in other galaxies. (e.g., Figure 4 in Visser 1980, Figure 5 in Aalto et al. 1999) This can probably happen at this region in our galaxy. (ii) Spiral shock may lead to the condition that one LSR velocity could share two different distances. (e.g., Foster & MacWilliams 2006) (iii) The expansion motion of the H i super-bubble near l = 123°, b = –6° associated with the Perseus arm may also lead to velocity mixing. (Sato et al. 2007) Whatever reason mainly causes the mixing LSR velocities, this at least suggests that in this region the gas motion is very peculiar.

Figure 7 shows the velocity-integrated intensity of H i emission. All the clouds and HMSFRs of the Outer arm are marked on it. Since the velocity ranges of the Outer arm are different at different longitudes, we need to define an lV function as the integrated velocity window. We have adopted two polynomial fitting curves as the outline of integrated range instead of the spiral lV projection. This is because the Outer arm velocity ranges are irregular and the polynomial fitting curve is more appropriate for defining the various ranges. The inset of Figure 7 shows the mass of molecular gas as a function of galactic longitude. The bin is 1 galactic longitude degree. The total mass of all the Outer arm clouds is ∼3.1 × 106 M. More details about gas distributions are presented in Section 4.4.

Figure 7.

Figure 7. Velocity-integrated intensity of H i from the CGPS corresponding to the Outer arm. The integrated velocity window is marked by the pink ribbon in Figure 5. The meanings of the symbols are consistent with Figure 5. The inset shows the mass of molecular gas as a function of galactic longitude: the solid line and the dashed line correspond to different units.

Standard image High-resolution image

Those three figures present a 3D view of the arm. The H i and H2 gases are roughly matched. The warp is obvious.

4.4. Gas Distribution

The gas distribution of the Milky Way is a frequently discussed topic. (e.g., Gordon & Burton 1976; Sodroski et al. 1987; Nakanishi & Sofue 2003 and their serial papers; Duarte-Cabral et al. 2015). Whereas, most works often discuss its global distribution (e.g., Burton et al. 1975), only a few have focused on the individual arm (e.g., Grabelsky et al. 1987). Now we present our results of the Outer arm gas distribution.

Figure 8 shows the Outer arm H i surface density distribution (or column density distribution) along the galactic longitude. Assuming that the 21 cm line is optically thin, the H i surface density is calculated from ${{\rm{\Sigma }}}_{{\rm{H}}{\rm{I}}}=1.82\times {10}^{18}{m}_{{\rm{H}}}\int {T}_{{\rm{B}}}(d)\;V$, and the column density is calculated by the conversion factor of 1 M pc−2 = 1.25 × 1020 cm−2, where mH is the H atom mass, and $\int {T}_{{\rm{B}}}d\;V$ is the integrated intensity. The velocity-integrated range is the pink ribbon shown in Figure 5. Clearly, the H i surface density abruptly descends at $l\simeq 100^\circ $, and ascends a little at $l\simeq 120^\circ $. The arm-blending region may account for the sharp drop in l = [100, 120]°, but cannot explain the low surface density region of l = [120, 155]°. This may suggest that in one spiral arm the gas quantity can largely change. One possible reason may be that: some of the apparent variation in H i surface density could arise from colder, optically thick H i gas, which is known to be widespread in the outer Galaxy, perhaps overlapping with some of the molecular clouds. (see Knee & Brunt 2001, Strasser et al. 2007, and Gibson 2010).

Figure 8.

Figure 8. Distribution of H i surface density (or column density) of the Outer arm along the galactic longitude. The lines with different colors indicate different galactic latitudes, and the black line indicates the mean value.

Standard image High-resolution image

Using the H i surface density and kinematic distance, the H imass is obtained from the following process. First, since we know the H i surface density at every pixel of the H i integrated map (namely Figure 7), we then can calculate the mean surface density at every square degree. (This step is somewhat like smoothing the integrated map into a one square degree per pixel map.) Second, using the Reid model we calculate the kinematic distance of every square degree, and the mean LSR velocity of the integrated window (namely the pink ribbon of Figure 5) at every galactic longitude is used as the LSR velocity of its corresponding square degree. Third, using the kinematic distance and the surface density we can calculate the H i mass in every square degree by ${M}_{{\rm{H}}{\rm{I}}}={{\rm{\Sigma }}}_{{\rm{H}}{\rm{I}}}{d}_{{\rm{H}}\;{\rm{I}}}^{2}$, where ΣH i and dH i respectively indicate the surface density and kinematic distance at the corresponding square degree. Figure 9 shows the H i and H2 mass distribution along the galactic longitude. (Note that H2 is different from molecular gas; the latter includes helium.) The bins of H i and H2 are both 1 galactic longitude degree. The total mass of H i in l = [63, 155]°, b = [−3, 5]° of the Outer arm is ∼1.4 × 108 M. And the total mass of H2 in l = [100, 150]°, b = [−3, 5]° is ∼2.3 × 106 M. Since the observation is not fully covered in b, the mass must be higher. The mean mass ratio of H2 to H i in l = [100, 150]°, b = [−3, 5]° is about 0.1. Figure 10 shows the H i and H2 mass distributions along the Z scale in the region of l = [100, 150]°, b = [−3, 5]°. About 50% of the gas mass is included in Z = [0.2, 0.4] kpc.

Figure 9.

Figure 9. Distribution of H i and H2 masses of the Outer arm along the galactic longitude.

Standard image High-resolution image
Figure 10.

Figure 10. Distribution of the H i and H2 masses along the scale height. This result only shows the region l = [100, 150]°, b = [−3, 5]°. The blue words beside the asterisks are the mass ratios of H2 to H i.

Standard image High-resolution image

Interestingly, there seems to be a trend in the Z scale distribution (Figure 10): as the total gas mass becomes large, the H2 to H i ratio becomes higher. For example, in Figure 10, at Z = 0.3 kpc, the total gas mass is the largest, and accordingly the H2 to H i ratio is the highest. On the other hand, at Z = 0.1 or 0.6 kpc, the gas mass descends and so does the ratio. Also, this trend is visible in the longitude distribution (Figure 9). One possible explanation may be that although the H2 distribution is not as diffuse as H i, once it exists, it contributes a lot to the total mass.

It is necessary to emphasize that we are reporting only the mass of H2 traced by detected CO emission. Significant additional H2 which CO cannot detect (or say the "CO-dark" H2, e.g., Grenier et al. 2005) may also be present. Most of the sight lines plotted in Figure 8 exceed the minimum column density for self-shielding H2 (a few times 1020 cm−2; see Snow & McCall 2006 and Sheffer et al. 2008). This is for individual clouds rather than integrated sight lines, but it still seems likely that H2 could exist in many sightlines where the total H i column is high enough.

4.5. Warp and Thickness

Warp is a general phenomenon among disk galaxies. (e.g., Sancisi 1976; Sanchez-Saavedra et al. 1990; Reshetnikov & Combes 1998; García-Ruiz et al. 2002). Additionally, it is one of the most fascinating and important subjects in the study of the Milky Way structure (e.g., Burke 1957; Reed 1996; Momany et al. 2006; López-Corredoira et al. 2010). Figure 3 in the review of Kalberla & Kerp (2009) shows a distinct three-dimensional (3D) picture of the warped Milky Way plane.

In our paper, one can easily notice the Milk Way warp from Figures 68 and 10. We can roughly obtain a first impression from Figure 7: the Outer arm starts to bend upward at $l\simeq 80^\circ $, and at $l\simeq 105^\circ $ it reaches its peak and then begins to fall. The whole arm from l = 63° to l = 155° just looks like an arch. In addition, almost all segments of this arch are beyond the b = 0° plane. But since the distance is different at every longitude, the l–b map cannot fully represent the scale height distribution. Now a more detailed analysis is presented. First, we plotted H i mass distributions along the Z scale every 2 galactic longitudes, and H2 mass distributions every 4 galactic longitudes. Then we fitted every distribution with a Gaussian curve. We define the centric position as the Z scale and the FWHM is the thickness. The final results are shown in Figure 11. The mean Z scales of the H i layer and the H2 layer are about 0.31 kpc and 0.25 kpc, respectively. The mean thicknesses of the H i and H2 layers are about 550 and 80 pc, respectively. The Z scales of H i and H2 are close (≃0.3 kpc), but the H i thickness is much larger than that of of H2 (about 7 times thicker, which is not obvious in Figure 7).

Figure 11.

Figure 11. Distribution of the Outer arm thickness (upper panel) and scale height (lower panel) along the galactic longitude.

Standard image High-resolution image

The increasing trend of spiral arm thickness with galactocentric radius has been widely observed and accepted. (e.g., Wouterloot et al. 1990; Kalberla & Kerp 2009). Here we detected a similar trend. Assuming a mean distance of 2 kpc,the H i thickness of the Perseus arm is about 200 pc, corresponding to 5°. And as mentioned above, the Outer arm thickness is about 550 pc. Meanwhile, the thickness of the New arm is 400–600 pc (Sun et al. 2015). Obviously the Outer arm is much thicker than the Perseus arm, but is nearly as thick as the New arm. Maybe this provides a little evidence to a trend newly discovered by Honig & Reid (2015) that in the outermost parts of the galaxies some arms become narrow.

5. SUMMARY

Combining 115 GHz 12CO(1–0) data (l = [100, 150]°, b = [−3, 5]°) of MWISP and 21 cm H i data (l = [65, 115]°, b = [−3, 5]°) of CGPS, we present the properties of the Outer arm in the second galactic quadrant of the Milky Way.

(1) Using CO(1–0), We have detected 481 molecular clouds in total, of which 332 (about 69%) are newly detected, 457 clouds are identified in the Outer arm, and 24 are identified in the New arm. The parameters of all 481 clouds are summarized in Table 1.

(2) Assuming that the spiral arm is logarithmic spiral, the pitch angle of the Outer arm is fitted by minimizing the chi-square error statistic algorithm, and the result is ∼13fdg1.

(3) The total masses of molecular gas and H2 from l = 100° to 150° in the Outer arm are about 3.1 × 106 M and 2.3 × 106 M, respectively. And the total mass of H i gas from l = 63° to 155° is about 1.4 × 108 M. Since the observation is not fully covered in b, the mass must be higher. The mean mass ratio of H2 to H i in l = [100, 150]°, b = [−3, 5]° is about 0.1.

(4) The warp of the Outer arm is obvious. The mean Outer arm thicknesses of H i and H2 are about 550 and 80 pc, respectively, while the scale height of both gases is ≃0.3 kpc.

We are grateful to all the members of the Milky Way Scroll Painting survey group, especially the staff of the Qinghai Radio Observing Station at Delingha for technical support. Additionally, we gratefully acknowledge the anonymous referee for the helpful and valuable comments. This work is supported by the National Natural Science Foundation of China (grant numbers 11133008 and 11233007), the Strategic Priority Research Program of the Chinese Academy of Sciences (grant number XDB09010300), and the Key Laboratory for Radio Astronomy.

Footnotes

Please wait… references are loading.
10.3847/0067-0049/224/1/7