arxiv: v2 [astro-ph.ga] 5 May 2016

Transcription

1 An Outer Arm in the Second Galactic Quadrant: Structure Xinyu Du 1,2,3, Ye Xu 1,3, Ji Yang 1,3, Yan Sun 1,2,3, Facheng Li 1,2,3, Shaobo Zhang 1,3, Xin Zhou 1,3 arxiv: v2 [astro-ph.ga] 5 May 2016 xydu@pmo.ac.cn; xuye@pmo.ac.cn ABSTRACT The lack of arm tracers, especially the 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 12 CO(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 M. Assuming that the spiral arm is a logarithmic spiral, the pitch angle is fitted as Besides combining both the CO data from MWISP and the 21 cm HI data from the Canadian Galactic Plane Survey (CGPS), the gas distribution, warp, and thickness of the Outer arm are also studied. Subject headings: catalogs Galaxy: structure ISM: molecules ISM: clouds 1. Introduction It is well-known that our galaxy has a spiral structure, (Oort et al. 1958; Bok 1959) but a detailed understanding of the Milky Way is still lacking. Since Georgelin & Georgelin 1 Purple Mountain Observatory, Chinese Academy of Science, Nanjing , China 2 Graduate University of the Chinese Academy of Sciences, 19A Yuquan Road, Shijingshan District, Beijing , China 3 Key Laboratory of Radio Astronomy, Chinese Academy of Science, Nanjing , China

2 2 (1976) presented a large-scale spiral pattern with HII 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 HI 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 outer galaxy with FCRAO CO survey results, providing the 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 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, the 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) project 1 provides such a chance to study the molecular gases with the no-bias high sensitive 12 CO(1 0), 13 CO(1 0) and C 18 O(1 0) observations. Due to such high-sensitive observations, a new spiral arm (hereafter the New arm) beyond the Outer arm has been discovered by Sun et al. (2015). Combining the 12 CO(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 13 CO(1 0) and C 18 O(1 0) data and further analyses will be published in our future works. Please note that the Outer arm does not have a particular name yet. Some authors have labeled it as 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 in this paper to avoid confusion. In Sect. 2 we introduce our CO observation conditions and archival data of atomic 1 or

3 3 hydrogen. In Sect. 3 we describe how we pick out the Outer arm clouds and briefly study their properties. In Sect. 4 we study the properties of the Outer arm, including the pitch angle, the gas distribution, the thickness, and the warp. The summary is given in Sect Observations and archival data 2.1. CO observations The 12 CO (1 0), 13 CO(1 0) and C 18 O(1 0) lines were observed simultaneously using the Purple Mountain Observatory Delingha (PMODLH) 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-insulatorsuperconductor (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 12 CO (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 MB = T A /η MB. The data were sampled every 30. All the data were reduced using the GILDAS/CLASS package 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 2. 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, which is comparable to our CO observations. 2

4 4 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 indentified 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 ) The parameters of all the 481 clouds are summarized in Table 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 V LSR 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 HI data from the CGPS, we plotted the longitude-velocity map of HI, 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 HI longitude-velocity map. (In other words, we converted the dashed curves shown in Fig. 3 into the ones shown in Fig. 5. And the converting method is presented in detail in Sect. 4.3.) And then we marked all the HCS01 clouds on that map. Combining the HI 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 clouds 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-

5 5 velocity map. However, it is necessary to point out that the Outer arm clouds that we found do not absolutely belong to that arm. From l 100 to 120, there exists an arm-blending region. The LSR velocities of HI gas of the Perseus arm and the Outer arm are mixed together (see Fig. 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 70 km s 1, or falsely pick the clouds that may be located at the inter-arm area 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 kinematic method. And the accuracy of the trigonometric parallax and photometry methods is better than that of 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 G 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, R 0 = 8.34 kpc, dθ = 0.2 km dr s 1 kpc 1, U = 10.7 km s 1, V = 15.6 km s 1, W = 8.9 km s 1, U s = 2.9 km s 1, 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 G The Reid model is much better than the IAU model (namely Θ 0 = 220 km s 1,R 0 = 8.5 kpc ). Take the cloud MWISP G as an example. Its parallax and luminosity distances both are 6.0 kpc, 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 the Section 4 of Reid et al. (2009) However, one should also keep in mind that there still exist the 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 G or the left panel of Fig. 3. A detailed discussion about the biases is presented in Sect. 4.2.

6 6 Knowing the distance d and latitude b we calculated the scale height Z by Z = d sin(b). And using b, Z, d and longitude l, we calculated the galactocentric radius R with R 2 = d 2 cos 2 (b) + R 2 0 2R 0 d cos(b) cos(l) + Z 2. The cloud solid angle A is defined by the 3σ limits. The cloud diameter D is obtained 4 after the beam deconvolution: D = d A π θ2 MB (Ladd et al. 1994), where θ MB is the main beam width. Adopting the CO-to-H 2 X factor cm 2 (K km s 1 ) 1 (Dame et al. 2001), cloud mass is calculated from M = 2 µ m H X π ( D 2 )2 T B d V, where µ = 1.36 (Hildebrand 1983) is the mean atomic weight per H atom in the ISM, m H is the H atomic mass, and T B d V is the integrated intensity of the T 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 (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 out small-size molecular clumps, which is not consistent with our cloud identification criterion. Fig.1 shows the results. One may notice that there are 102 clouds of HCS01 in Fig. 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. Since the cloud heliocentric distances of HCS01 are derived from the IAU model, the sizes and masses are correspondingly biased. In order to 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

7 7 ratio, we can detect a larger angle area for the same cloud, which results in larger sizes and masses. Fig. 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. 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 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 ln( R R ref ) = (β β ref ) tan(ψ) (1) where ψ is the spiral pitch angle, β is the galactocentric azimuth (namely the Source-GC-Sun angle), and R ref and β ref are the reference radius and reference azimuth, respectively. The

8 8 fitting algorithm is a minimizing chi-square error statistic, and the chi-square error statistic is computed as N χ 2 (k, b) = W i (y i b kx i ) 2 (2) i=1 where k = tan(ψ), b = β ref tan(ψ) + ln(r ref ), x i and y i are β and ln(r) of each cloud, respectively, and W i is the weight. The weight is defined as W = log 10 ( M M ). Fig.2 shows the fitting result. The pitch angle is 13.1, the R ref is 13.6 kpc, the β ref is 26.9 and the χ 2 is 4.1. Our fitting pitch angle is close to the result of Reid et al. (2014) (pitch=13.8 ), 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 Fig. 2. Fitting result of the Outer arm pitch angle. R is Galactocentric radius in the 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.

9 Plan view Since Georgelin & Georgelin (1976) presented the famous Milky Way plan 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 Fig. 3 shows the plan 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 R ref 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 Fig. 3 is the result after the move. We mentioned these distance biases in 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 Fig. 3.) Besides, the Reid model is a global model and it is fitted by the 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 Fig. 3.

10 10 Fig. 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 HII 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 Sect. 3.2 are retained.

11 11 Fig. 3. (left panel) Locations of molecular clouds and HMSFRs in the plan 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 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 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.

12 12 Fig. 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 > 10 3 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.

13 l-v, V-b and l-b map Fig. 5 shows the l-v map of HI 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 Fig. 3 into the ones shown in Fig. 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 Sect 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 Eq. 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 l-v curves of the arms. Fig. 6 shows the V-b montage of HI emission. The high-mass Outer arm clouds (> 10 3 M ), the New arm clouds, and the HMSFRs are marked on it. As mentioned in Sect. 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 has 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 HI super-bubble near l = 123, b = 6 associated with the Perseus arm may also lead to the 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. Fig. 7 shows the velocity-integrated intensity of HI emission. All the clouds and HMS- FRs 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 l V function as the integrated velocity window. We have adopted two polynomial fitting curves as the outline of integrated range instead of the spiral l V 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 Fig. 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 M. More details about gas distributions are presented in Sect. 4.4.

14 14 Those three figures present a 3D view of the arm. The HI and H 2 gases are roughly matched. The warp is obvious. Fig. 5. Longitude-velocity diagram of HI 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 Fig. 7.

15 15 Fig. 6. Velocity-latitude montage of HI 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 Fig. 5. Note that the low-mass clouds(< 10 3 M ) are not plotted, and the symbol size does NOT indicate mass.

16 16 Fig. 7. Velocity integrated intensity of HI from the CGPS corresponding to the Outer arm. The integrated velocity window is marked by the pink ribbon in Fig. 5. The meanings of the symbols are consistent with Fig. 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.

17 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. Fig. 8 shows the Outer arm HI surface density distribution (or column density distribution) along the galactic longitude. Assuming that the 21 cm line is optically thin, the HI surface density is calculated from Σ HI = m H TB d V, and the column density is calculated by the conversion factor of 1 M pc 2 = cm 2, where m H is the H atom mass, and T B d V is the integrated intensity. The velocity-integrated range is the pink ribbon shown in Fig. 5. Clearly, the HI surface density abruptly descends at l 100, and ascends a little at l 120. 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 HI surface density could arise from colder, optically thick HI 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 and Gibson 2010) Using the HI surface density and kinematic distance, the HI mass is obtained from the following process. First, since we know the HI surface density at every pixel of the HI integrated map (namely Fig. 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 Fig. 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 HI mass in every square degree by M HI = Σ HI d 2 HI, where Σ HI and d HI respectively indicate the surface density and kinematic distance at the corresponding square degree. Fig. 9 shows the HI and H 2 mass distribution along the galactic longitude. (Note that H 2 is different from molecular gas; the latter includes helium.) The bins of HI and H 2 are both 1 galactic longitude degree. The total mass of HI in l = [63, 155], b = [ 3, 5] of the Outer arm is M. And the total mass of H 2 in l = [100, 150], b = [ 3, 5] is M. Since the observation is not fully covered in b, the mass must be higher. The mean mass ratio of H 2 to HI in l = [100, 150], b = [ 3, 5] is about 0.1. Fig. 10 shows the HI and H 2 mass distributions along the Z scale in the region of l = [100, 150], b = [ 3, 5]. About 50%

18 18 of the gas mass is included in Z=[0.2,0.4] kpc. Interestingly, there seems to be a trend in the Z scale distribution (Fig. 10): as the total gas mass becomes large, the H 2 to HI ratio becomes higher. For example, in Fig. 10, at Z = 0.3 kpc, the total gas mass is the largest, and accordingly the H 2 to HI 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 (Fig. 9). One possible explanation may be that although the H 2 distribution is not as diffuse as HI, once it exists, it contributes a lot to the total mass. It is necessary to emphasize that we are reporting only the mass of H 2 traced by detected CO emission. Significant additional H 2 which CO cannot detect (or say the CO-dark H 2, e.g., Grenier et al. 2005) may also be present. Most of the sight lines plotted in Fig. 8 exceed the minimum column density for self-shielding H 2 (a few times 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 H 2 could exist in many sightlines where the total HI column is high enough. Fig. 8. Distribution of HI 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.

19 19 Fig. 9. Distribution of HI & H 2 masses of the Outer arm along the galactic longitude.

20 20 Fig. 10. Distribution of the HI & H 2 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 H 2 to HI.

21 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 Fig. 6, 7, 8 and 10. We can roughly obtain a first impression from Fig. 7: the Outer arm starts to bend upward at l 80, and at l 105 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 HI mass distributions along the Z scale every 2 galactic longitudes, and H 2 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 Fig. 11. The mean Z scales of the HI layer and the H 2 layer are about 0.31 kpc and 0.25 kpc, respectively. The mean thicknesses of the HI layer and the H 2 layer are about 550 pc and 80 pc, respectively. The Z scales of HI and H 2 are close ( 0.3 kpc), but the HI thickness is much larger than that of H 2 (about 7 times thicker, which is not obvious in Fig. 7). 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 HI thickness of the Perseus arm is about 200 pc corresponding to 5. And as it mentioned above, the Outer arm thickness is about 550 pc. Meanwhile, the thickness of the New arm is 400 to 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.

22 22 Fig. 11. Distribution of the Outer arm thickness (Upper panel) and scale height (lower panel) along the galactic longitude.

23 23 5. Summary Combining 115 GHz 12 CO(1-0) data (l = [100, 150], b = [ 3, 5] ) of MWISP and 21 cm HI data (l = [65, 115], b = [ 3, 5] ) of CGPS, we present the properties of 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 the 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 chi-square error statistic algorithm, and the result is (3) The total masses of molecular gas and H 2 from l = 100 to 150 in the Outer arm are about M and M respectively. And the total mass of HI gas from l = 63 to 155 is about M Since the observation is not fully covered in b, the mass must be higher. The mean mass ratio of H 2 to HI 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 HI and H 2 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: and ), the Strategic Priority Research Program of the Chinese Academy of Sciences (grant number: XDB ), and the Key Laboratory for Radio Astronomy.

24 24 REFERENCES Aalto, S., Hüttemeister, S., Scoville, N. Z., & Thaddeus, P. 1999, ApJ, 522, 165 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJ, 710, 133 Bobylev, V. V., & Bajkova, A. T. 2014, MNRAS, 437, 1549 Bok, B. J., 1959, The Observatory, 79, 58 Brand, J., & Wouterloot, J. G. A. 1994, A&AS, 103, 503 s Brunt, C. M., Kerton, C. R., & Pomerleau, C. 2003, ApJS, 144, 47 Burke, B. F. 1957, AJ, 62, 90 Burton, W. B., Gordon, M. A., Bania, T. M., & Lockman, F. J. 1975, ApJ, 202, 30 Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 Duarte-Cabral, A., Acreman, D. M., Dobbs, C. L., et al. 2015, MNRAS, 447, 2144 Foster, T., & MacWilliams, J. 2006, ApJ, 644, 214 García-Ruiz, I., Sancisi, R., & Kuijken, K. 2002, A&A, 394, 769 Georgelin, Y. M., & Georgelin, Y. P. 1976, A&A, 49, 57 Gibson, S. J. 2010, The Dynamic Interstellar Medium: Galactic Plane Survey, 438, 111 A Celebration of the Canadian Gordon, M. A., & Burton, W. B. 1976, ApJ, 208, 346 Grabelsky, D. A., Cohen, R. S., Bronfman, L., Thaddeus, P., & May, J. 1987, ApJ, 315, 122 Grenier, I. A., Casandjian, J.-M., & Terrier, R. 2005, Science, 307, 1292 Hachisuka, K., Choi, Y. K., Reid, M. J., et al. 2015, ApJ, 800, 2 Henderson, A. P., Jackson, P. D., & Kerr, F. J. 1982, ApJ, 263, 116 Heyer, M. H., Brunt, C., Snell, R. L., et al. 1998, ApJS, 115, 241 Heyer, M. H., Carpenter, J. M., & Snell, R. L. 2001, ApJ, 551, 852 Hildebrand, R. H. 1983, QJRAS, 24, 267

25 25 Honig, Z. N., & Reid, M. J. 2015, ApJ, 800, 53 Hou, L. G., & Han, J. L. 2014, A&A, 569, A125 Kalberla, P. M. W., & Kerp, J. 2009, ARA&A, 47, 27 Knee, L. B. G., & Brunt, C. M. 2001, Nature, 412, 308 Ladd, E. F., Myers, P. C., & Goodman, A. A. 1994, ApJ, 433, 117 López-Corredoira, M., Cabrera-Lavers, A., Hammersley, P. L., et al. 2010, Highlights of Astronomy, 15, 811 Moffat, A. F. J., & Vogt, N. 1975, A&AS, 20, 85 Momany, Y., Zaggia, S., Gilmore, G., et al. 2006, A&A, 451, 515 Nakanishi, H., & Sofue, Y. 2003, PASJ, 55, 191 Negueruela, I., & Marco, A. 2003, A&A, 406, 119 Oort, J. H., Kerr, F. J., Westerhout, G. 1958, MNRAS, 118, 379 Reed, B. C. 1996, AJ, 111, 804 Reid, M. J., Menten, K. M., Zheng, X. W., et al. 2009, ApJ, 700, 137 Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130 Reshetnikov, V., & Combes, F. 1998, A&A, 337, 9 Russeil, D. 2003, A&A, 397, 133 Russeil, D., Adami, C., & Georgelin, Y. M. 2007, A&A, 470, 161 Sanchez-Saavedra, M. L., Battaner, E., & Florido, E. 1990, MNRAS, 246, 458 Sancisi, R. 1976, A&A, 53, 159 Sato, M., Hirota, T., Honma, M., et al. 2007, PASJ, 59, 743 Shan, W.L., Yang, J., Shi, S. C., et al. 2012, IEEE Transactions on Terahertz Science and Technology, 2,593 Sheffer, Y., Rogers, M., Federman, S. R., et al. 2008, ApJ, 687, 1075 Snow, T. P., & McCall, B. J. 2006, ARA&A, 44, 367

26 26 Sodroski, T. J., Dwek, E., Hauser, M. G., & Kerr, F. J. 1987, ApJ, 322, 101 Steiman-Cameron, T. Y. 2010, Galaxies and their Masks, 45 Strasser, S. T., Dickey, J. M., Taylor, A. R., et al. 2007, AJ, 134, 2252 Sun, Y., Xu, Y., Yang, J., et al. 2015, ApJ, 798, L27 Taylor, J. H., & Cordes, J. M. 1993, ApJ, 411, 674 Taylor, A. R., Gibson, S. J., Peracaula, M., et al. 2003, AJ, 125, 3145 Vallée, J. P. 2008, AJ, 135, 1301 Vallée, J. P. 2015, MNRAS, 450, 4277 Visser, H. C. D. 1980, A&A, 88, 159 Wouterloot, J. G. A., Brand, J., Burton, W. B., & Kwee, K. K. 1990, A&A, 230, 21 Xu, Y., Reid, M. J., Zheng, X. W., & Menten, K. M. 2006, Science, 311, 54 This preprint was prepared with the AAS L A TEX macros v5.2.

27 27 Table 1. Parameters of molecular clouds derived by CO(1-0) Name V LSR T peak V I CO Area d R Size Mass Z Scale Note (km s 1 ) (K) (km s 1 ) (K km s 1 ) (arcmin 2 ) (kpc) (kpc) (pc) (10 3 M ) (kpc) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G BKP03

28 28 Table 1 Continued Name V LSR T peak V I CO Area d R Size Mass Z Scale Note (km s 1 ) (K) (km s 1 ) (K km s 1 ) (arcmin 2 ) (kpc) (kpc) (pc) (10 3 M ) (kpc) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G

29 29 Table 1 Continued Name V LSR T peak V I CO Area d R Size Mass Z Scale Note (km s 1 ) (K) (km s 1 ) (K km s 1 ) (arcmin 2 ) (kpc) (kpc) (pc) (10 3 M ) (kpc) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) MWISP G MWISP G MWISP G MWISP G MWISP G MWISP G HCS01;BKP03 MWISP G MWISP G MWISP G MWISP G HCS01 MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G MWISP G MWISP G HCS01;BKP03 MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G MWISP G HCS01;BKP03 MWISP G HCS01;BKP03 MWISP G MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G MWISP G HCS01;BKP03 MWISP G MWISP G MWISP G MWISP G BKP03 MWISP G BKP03 MWISP G MWISP G BKP03 MWISP G