arxiv: v4 [astro-ph.sr] 25 Jun 2015

Transcription

1 The two molecular clouds in RCW 38; evidence for formation of the youngest super star cluster in the Milky Way triggered by cloud-cloud collision arxiv: v4 [astro-ph.sr] 25 Jun 2015 Y. Fukui 1, K. Torii 1, A. Ohama 1, K. Hasegawa 1, Y. Hattori 1, H. Sano 1,2, S. Ohashi 3,4, K. Fujii 3,4, S. Kuwahara 3,4, N. Mizuno 4, J. R. Dawson 5,6, H. Yamamoto 1, K. Tachihara 1, T. Okuda 4, T. Onishi 7, and A. Mizuno 8 1 Department of Physics and Astrophysics, Nagoya University, Chikusa-ku, Nagoya, Aichi , Japan 2 Institute for Advanced Research, Nagoya University, Chikusa-ku, Nagoya, Aichi , Japan 3 Department of Astronomy, School of Science, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo , Japan 4 National Astronomical Observatory of Japan, Mitaka, Tokyo , Japan 5 Department of Physics and Astronomy and MQ Research Centre in Astronomy, Astrophysics and Astrophotonics, Macquarie University, NSW 2109, Australia 6 Australia Telescope National Facility, CSIRO Astronomy and Space Science, P.O. Box 76, Epping, NSW 1710, Australia 7 Department of Astrophysics, Graduate School of Science, Osaka Prefecture University, 1-1 Gakuen-cho, Nakaku, Sakai, Osaka , Japan 8 Solar-Terrestrial Environment Laboratory, Nagoya University, Chikusa-ku, Nagoya , Japan torii@a.phys.nagoya-u.ac.jp Received ; accepted

2 2 ABSTRACT We present distributions of two molecular clouds having velocities of 2kms 1 and 14kms 1 toward RCW38, the youngest super star cluster in the Milky Way, inthe 12 CO(J =1 0and3 2)and 13 CO(J =1 0)transitions. Thetwo cloudsare likely physically associated with the cluster as verified by the high intensity ratio of the J=3 2 emission to the J=1 0 emission, the bridging feature connecting the two clouds in velocity and their morphological correspondence with the infrared dust emission. The total mass of the clouds and the cluster is too small to gravitationally bind the velocity difference. We frame a hypothesis that the two clouds are colliding with each other by chance to trigger formation of the 20 candidate O stars which are localized within 0.3pc of the cluster center in the 2kms 1 cloud. We suggest that the collision is currently continuing toward part of the 2kms 1 cloud where the bridging feature is localized. This is the third super star cluster alongside of Westerlund 2 and NGC 3603 where cloud-cloud collision triggered the cluster formation. RCW 38 is the most remarkable and youngest cluster, holding a possible sign of on-going O star formation, and is the most promising site where we may be able to witness the moment of O-star formation. Subject headings: ISM: clouds radiolines: ISM

3 3 1. Introduction 1.1. Important role of high-mass stars High-mass stars are very energetic and influential on dynamics of the interstellar medium via stellar winds, ultraviolet radiation and supernova explosions. Injection of the heavy elements at the end of their lives is another important effect of high-mass stars in galactic evolution. It is therefore one of the most important issues in astrophysics to understand formation of high-mass stars and considerable efforts have been made by a number of works published so far (for recent reviews see Zinnecker & Yorke 2007; Tan et al. 2014) Theoretical attempts for understanding high-mass star formation It is now generally accepted that most of the high-mass stars are formed by mass accretion instead of collisional stellar merging after some years of debates. In the current paradigm, monolithic collapse (Yorke & Sonnhalter 2002) and competitive accretion (Bonnell et al. 2004) are the two scenarios of the mass accretion processes. In the monolithic collapse a massive dense molecular cloud is the initial condition for the star formation, which is followed by disk accretion, whereas in the competitive accretion a massive cluster of stars/pre-stellar cores is assumed to grow in mass by accretion of the ambient gas. The two scenarios still lack detailed confrontation with the initial conditions provided by observations. This is in part due to the lack of well-resolved observational views of high-mass star formation at small distances; even the Orion nebular cluster, the most outstanding high-mass star forming region at the smallest distance 400 pc, is three times more distant than typical low-mass star-forming regions like Taurus, and the second nearests are at a distance around 2 kpc. Therefore, our understanding of high-mass star

4 4 formation yet remains to be significantly corroborated. It is interesting to note that insight-full discussion is made on the mechanisms of high-mass star formation in the review article by Zinnecker & Yorke (2007), where we find the following sentences; Rapid external shock compression (i.e., supersonic gas motions) generating high-column densities in less than a local free-fall time rather than slow quasi-static build-up of massive cores may be the recipe to set up the initial conditions for local and global bursts of massive star formation. Such bursts can both enhance or quench further massive star formation (positive and negative feedback), depending on environment (gas density and pressure), but such a discussion is beyond the scope of this review. This idea of supersonic gas motions to collect gas was not explored further in the literature in the context of high-mass star formation except for the several papers on cloud-cloud collision discussed in Section Increasing evidence for cloud-cloud collision as a triggering of high-mass star formation Remarkably, molecular observations with NANTEN2 have shown that O-type star formation in the Milky Way disk is triggered by supersonic collisions between two clouds at velocity of 10 30kms 1 in several places. The first discovery was made in the super star cluster Westerlund2 which is ionizing an Hii region RCW49 (Furukawa et al. 2009; Ohama et al. 2010), and the authors suggested that two giant molecular clouds collided to trigger formation of the cluster some 2 Myrs ago. These observations were followed by discovery of cloud-cloud collisions which triggered formation of the super star cluster/mini-starburst in NGC 3603 (Fukui et al. 2014), and formation of the single O star in Hii regions, M20 (Torii et al. 2011) and RCW120 (Torii et al. 2015). Most recently, in the Large Magellanic Cloud ALMA observations led to a discovery of formation of a 40M

5 5 star N159W-S triggered by collision between two filamentary clouds at a relative velocity of 10kms 1 (Fukui et al. 2015). These results raise a novel possibility that formation of O stars may be due to, or at least facilitated, by the supersonic compression in cloud-cloud collisions, which was not considered in the theories for high-mass star formation, either in the monolithic collapse or in the competitive accretion scenarios. On a global scale of galaxies, numerical simulations show that collisions between clouds are fairly frequent at every 10Myrs (Tasker & Tan 2009; Dobbs et al. 2015), which is similar to an evolutionary timescale of giant molecular clouds (Fukui & Kawamura 2010). It is thus becoming possible to consider cloud-cloud collision as one of the modes of high-mass star formation. Cloud-cloud collision was studied using hydrodynamics simulations by Habe & Ohta (1992) and Anathpindika (2010). These authors found that cloud-cloud collision induces the formation of cloud cores by enhanced self-gravity as a consequence of shock compression, whereas the connection of cloud-cloud collision to high-mass star formation was not discussed into depth by these authors. Recent magnetohydrodynamical (MHD) numerical simulations on cloud-cloud collision show that colliding molecular gas can indeed create dense and massive cloud cores, precursors of high-mass stars, in the shock-compressed interface, giving a theoretical explanation for the O-star formation in the above five cases (Inoue & Fukui 2013). The supersonic collision amplifies turbulence and magnetic field significantly and creates high-mass-accretion rate two orders of magnitude larger than that prior to the collision; in the simulations after the collision the density increases from cm 3 to 10 4 cm 3 and the turbulent velocity and magnetic field increase by about a factor of five, leading to a 100-times higher mass accretion rate to form stars. This indicates that the conditions for high-mass accretion rate postulated in high-mass star formation theories (McKee & Tan 2003) are created by cloud-cloud collision. Following these findings, as a next step we need to have a better understanding on the

6 6 physical details of the collision, which include direct observations of the shocked interface as well as the statistical properties of cloud-cloud collisions. We are carrying out a systematic study of high-mass star formation toward Hii regions in order to elucidate the role of cloud-cloud collision with NANTEN2, Mopra, and ASTE telescopes. These Hii regions include Spitzer bubbles, compact Hii regions, and mini-starbursts including NGC 6334, NGC6357, W33, etc. In order to place cloud-cloud collision in a wide context of star formation, we have the following specific questions to be answered; Q.I How the mass and number of the formed O stars are determined in the collision; how is the dependence of the stellar properties on the collision parameters such as the relative velocity between the clouds, cloud mass, cloud density, and cloud morphology. Q. II What is the origin of the mass segregation/stellar distribution in super star clusters/mini-starbursts? Q. III What are the density, velocity field and magnetic field in the shocked interface layer. This is most relevant to our understanding of details of high-mass protostellatr formation/evolution. We note that in the above four objects in the Milky Way molecular gas is already ionized in the central pc of the O star(s), hindering one to see directly the initial conditions of the natal gas. Q. IV Can cloud-cloud collision also trigger formation of low- to intermediate-mass stars? If yes, what is the mass function of the formed stars. Q.V What is the fraction of O stars formed by cloud-cloud collisions in the Milky Way? Are cloud-cloud collisions one of the major modes of O-star formation?

7 RCW38 the youngest super star cluster in the Milky Way Hoping to answer the above five questions, we in this paper focus ourselves on the super star cluster, RCW38 (Rodgers et al. 1960), youngest in the Milky Way (see Table 2 of super star clusters in Portegies Zwart et al. 2010). The age of RCW38 is assumed to be 0.5Myrs so that it does not conflict with the near infrared colors (Wolk et al. 2006), while the age may be even smaller as discussed in the present paper (Section 4.3). The extreme youth will benefit study of the molecular gas close to the central O stars, which is particularly important to address Q. III which depends on the ionization/dispersal of the molecular gas. RCW38 has large total mass, similar to that of the ONC, and is the active high-mass star forming region, second nearest to the sun, at a distance of 1.7kpc. An overview of the previous studies on RCW38 is summarized in an article in the Handbook of Star Formation (Wolk et al. 2008). RCW38 is known as one of the closest high-mass star forming regions containing more than thousand stars (Wolk et al. 2006; Lada & Lada 2003) (Figure1). At the central part of RCW 38, two remarkable infrared peaks have been identified (Frogel & Persson 1974). The brightest at 2µm is labeled IRS2, which corresponds to the O5.5 binary located at the center of the RCW38 cluster (DeRose et al. 2009). Furniss et al. (1975) derived a total infrared luminosity toward IRS2 to be L. On the other hand, the brightest feature in 10µm is found 0.1pc west of IRS2, and is labelled IRS1. IRS1 is a dust ridge extending by pc in the north-south direction (Figure1b), which has a color temperature of about 175 K and includes several dust condensations (Smith et al. 1999). Infrared and millimeter radio observations has indicated a ring-like or horseshoe structure around IRS2 about 1 2 across (Figure1; Huchtmeier 1974; Vigil 2004; Wolk et al. 2006, 2008). Inside the ring-like shape two regions are cleared of dust and forms cavities, one is the region centered on IRS2 with a diameter of 0.1pc, and another is just west of IRS1

8 8 with a similar size (Smith et al. 1999; Wolk et al. 2008), suggesting that they were formed by feedbacks from high-mass stars. A large-area infrared image shows that RCW 38 includes numerous filaments and bubbles (see Figure 1a). Many observational studies to investigate the cluster members of RCW 38 have been carried out so far. Near infrared observations by using the Very Large Telescope (VLT) identified more than 300 young stars for the 0.5pc 2 area centered on IRS2 (DeRose et al. 2009). Chandra observations reported by Wolk et al. (2006) was used to identify 345 X-ray sources likely associated with RCW 38. Completeness arguments suggest a total estimated cluster size of between 1500 and 2400 stars. Winston et al. (2011) identified 624 young stellar objects (YSOs) by the Spitzer/IRAC observations toward a region of RCW38. In addition, about sixty O-star candidates have been identified by Wolk et al. (2006) and Winston et al. (2011) in a large area of RCW38, and assumed an age of the RCW38 cluster to be 0.5Myrs. Out of them, about twenty O-star candidates are found in the central 1pc (Wolk et al. 2006). Molecular line observations toward RCW 38 are very limited. Yamaguchi et al. (1999) observed a large area of RCW38 in 13 CO J=1 0 emission and found a molecular cloud at 2 3kms 1 which coincides with IRS2, although detailed structures including the ring-like component were not resolved due to a relatively coarse resolution. On the other hand, Zinchenko et al. (1995) observed the central region of RCW38 at HPBW of 1 in the CS emission to resolve the ring shape surrounding IRS2. 12 CO J=1 0 observations with SEST at 45 resolution for a small region performed by Gyulbudaghian & May (2008) show two molecular clouds likely associated with RCW 38 at two different velocity ranges of 3 +2kms 1 and +3 +8kms 1. In this paper, we present new molecular line observations toward RCW 38 obtained with NANTEN2, Mopra, and ASTE telescopes. The NANTEN2 CO J=1 0 observations covered a large area of RCW38 at an angular resolution of 4 to describe the large

9 9 scale molecular morphologies and velocities, while the Mopra CO J=1 0 and ASTE CO J=3 2 observations at angular resolutions of revealed the detailed molecular gas distributions and physical conditions such as temperature and density around IRS 1 and IRS2. 2. Observations 2.1. NANTEN2 CO J=1 0 observations The NANTEN2 4-m mm/sub-mm telescope situated in Chile was used to observe a large area of RCW38 in the 12 CO and 13 CO J=1 0 transitions during from May 2012 to December The 4K cooled SIS mixer receiver equipped with NANTEN2 provided a typical system temperature T sys of 250K in DSB. The backend was a digital spectrometer which provides channels at 1 GHz bandwidth and 61 khz resolution, corresponding to 2600kms 1 and 0.17kms 1, respectively, at 110GHz. The obtained data were smoothed to a velocity resolution of 0.8kms 1 and angular resolution of 240. The pointing accuracy was checked to be better than 15 with daily observations toward the Sun and IRC10216 (R.A., Dec.)=(9 h 47 m 57 ṣ 406, ). The absolute intensity calibration was done with daily observations of the ρ-ophiuchus molecular cloud (Ra., Dec.)=(16 h 32 m 23 ṣ 3, ) and Perseus molecular cloud (Ra., Dec.)=(3 h 29 m 19 ṣ 0, ) Mopra CO J=1 0 observations Detailed CO J=1 0 distributions around RCW 38 were obtained by using the 22-m ATNF (Australia Telescope National Facility) Mopra mm telescope in Australia at a high spatial resolution of 33 during July We simultaneously observed 12 CO J=1 0, 13 CO

10 10 J=1 0, and C 18 O J=1 0 toward a area of RCW38 with the OTF mode with a unit field of 4 4. T sys was K in the SSB. The Mopra backend system MOPS which provides 4096 channels across MHz in each of the two orthogonal polarizations was used in the observations. The velocity resolution was 0.088kms 1 and the velocity coverage was 360kms 1 at 115GHz. The obtained data was smoothed to be HPBW of 40 with a 2D Gaussian function and to a 0.6kms 1 velocity resolution. The pointing accuracy was checked every 1 hour to keep within 7 with observations of 86GHz SiO masers. We made daily observations of Orion-KL (R.A., Dec.)=( 5 h 35 m 14 ṣ 5, ) to estimate extended beam efficiency by making comparisons with the peak temperature of 100 K in Orion-KL shown in Figure 6 of Ladd et al. (2005). We finally obtained an extend beam efficiency of ASTE 12 CO J=3 2 observations 12 CO J=3 2 line observations of RCW38 were performed by using the ASTE 10-m telescope located in Chile in 2013 (Ezawa et al. 2004, 2008). The waveguide-type sideband-separating SIS mixer receiver for the single sideband (SSB) CAT345 having system temperature of 250K (Inoue et al. 2008) and the digital spectrometer MAC with the narrow-band mode providing 128 MHz bandwidth and MHz resolution (Sorai et al. 2000), which corresponds to 111kms 1 velocity coverage and 0.11kms 1 velocity resolution at 345 GHz, was used. The observations were made with the OTF mode at a grid spacing of 7.5, and the HPBW was 22 at the CO J=3 2 frequency. The pointing accuracy was checked every 1.5hours to keep within 5 with observations of RAFGL5254 (R.A., Dec.)=(09 h 13 m 53 ṣ 94, ). The absolute intensity calibration was done with observations of Orion-KL (R.A., Dec.)=( 5 h 35 m 14 ṣ 5, ), and the day-to-day fluctuations of the peak intensity were within 10%.

11 11 3. Results 3.1. Large-scale molecular distributions with NANTEN2 We first present large-scale molecular distributions of RCW 38 using the NANTEN2 12 CO J=1 0 dataset. Compared with the NANTEN observations given by Yamaguchi et al. (1999), although the beam size of 180 is the same, the grid spacing of 60 in the present NANTEN2 observations is quarter, which enables us to reveal the gas distributions into more detail. As shown in the CO distributions in Figure2, toward the direction of RCW38, two CO clouds are distributed at radial velocities v LSR of 4 +8kms 1 and kms 1. In Figure 2(a), the blue-shifted cloud has a remarkably strong peak just coinciding with the center part of RCW38. The compact peak with a size of 0. 1 ( 3pc at RCW38) has a feature elongated toward the southwest of RCW 38, showing a characteristic head-tail structure, which was also detected in Yamaguchi et al. (1999). The head-tail structure is surrounded by the diffuse CO emission distributed around v LSR of 0 +8kms 1, covering the present observed region. Contrary to the blue-shifted cloud, as shown in Figure 2(b), the red-shifted cloud has weak and dispersed CO distribution. A CO peak at (l,b , 1. 1) with relatively strong emission is located just at the southeast of RCW 38, pointing toward the center of RCW38. This CO peak is well covered with the Mopra and ASTE observations as indicated by a box with solid lines in Figure2(b). There are several other CO components around the CO peak, and they all seem to be continuously distributed along the east-west direction roughly at a range of (l,b) (267. 5, 1. 2) (268. 4, 0. 9). In the position-velocity diagram in Figure 2(c), the red-shifted cloud can be identified separately from, although it looks being connected with, the blue-shifted cloud and the diffuse emission, suggesting a

12 12 possible physical relationship of the red-shifted cloud with RCW Detailed molecular distributions with Mopra and ASTE We present the detailed molecular distributions around RCW 38 using the Mopra and ASTE datasets. In Figure 3, the integrated intensity maps (left and center panels) and the position-velocity maps (right panels) are shown with 12 CO J=3 2 (upper panels) and 13 CO J=1 0 (lower panels). In Figures 3(a) and 3(d), the outstanding compact peak seen in Figure2(a) is resolved. It shows a ring-like structure having a diameter of 1 1.5pc, forming a cavity centered on IRS 2, and several filamentary and bubble-like structures are radially extended outside the ring-like structure by 1 2pc. In the 12 CO position-velocity diagram in Figure3(c), intensity decrease is seen around v LSR of 0 2kms 1. The 13 CO in Figure3(f) shows a single peak at 2kms 1, indicating that the dip at the same velocity in the 12 CO is due to self-absorption 1. The blue-shifted cloud has a peak at 2kms 1. The ring-like structure in 13 CO consists of four clumps (clumps I IV), and the ring-like structure itself has an elongated shape along the northeast-southwest direction (Figure 4). The sizes of the clumps range from 0.3pc to 0.5pc, and v from 4kms 1 to 5kms 1 (Table 1). In contrast to the complicated distributions of the 2kms 1 cloud, the 14kms 1 cloud has a simple ridge elongated from the north to the south (Figures3b and 3e). We shall hereafter call the 2kms 1 cloud and the 14kms 1 cloud as the ring cloud and the finger cloud, respectively, by their shapes. The finger cloud is not likely affected by self-absorption in 12 CO, since 12 CO and 13 CO show similar distributions in the position-velocity map in Figures 3(c) and 3(f). At its bottom, the ridge intersects another elongated structure 1 The two velocity features reported by Gyulbudaghian & May (2008) in 12 CO J=1 0 correspond to the two peaks caused by self-absorption.

13 13 running from the east to the west, and these two orthogonal elongated structures comprise the finger cloud in the present observed area. As shown in Figures3(c) and 3(f), the finger cloud shows a uniform velocity gradient of 1kms 1 pc 1 along the ridge. We find three connecting features in velocity between the CO emission of the ring cloud and the finger cloud at Dec , and , respectively, as indicated by arrows in Figure 3(c) and they are not seen beyond the velocity of the blue cloud. Their positions are shown in Figure Intensity ratios of 12 CO J=3 2 to 12 CO J=1 0 In order to investigate the physical properties of the ring cloud and the finger cloud, we present distributions of the intensity ratio of 12 CO J=3 2 to 12 CO J=1 0 (hereafter R 3 2/1 0 ). Taking intensity ratios between different J-levels of CO is a useful diagnose of the molecular gas properties (e.g., Ohama et al. 2010; Torii et al. 2011; Fukui et al. 2014). As seen in the R 3 2/1 0 position-velocity diagram in Figure6, both of the ring cloud and the finger cloud have typical ratios of , up to over 1.0. Since the ring cloud is strongly affected by the self-absorption in 12 CO spectra around its central velocity range, 0 4kms 1, the ratio around this velocity range is not reliable. On the other hand, low intensity ratios are seen at the northern part (Dec.> ) and southern part (Dec.< ) of the ring cloud, where the diffuse CO emission is widely distributed, while the finger cloud retains its relatively high intensity ratios of throughout the cloud.

14 Temperature and density of the molecular gas We here utilize the large velocity gradient (LVG) analysis (e.g., Goldreich & Kwan 1974) to estimate kinetic temperature T k and molecular number density n(h 2 ) of the finger cloud and the ring cloud to interpret their intensity ratio distributions. We first present curves of R 3 2/1 0 as a function of T k in Figure7 for various densities from 10 2 cm 3 to 10 4 cm 3, where two values of X(CO)/(dv/dr), 10 5 (kms 1 pc 1 ) 1 and (kms 1 pc 1 ) 1, are assumed. We here adopt an abundance ratio X(CO) = [ 12 CO]/[H 2 ] = 10 4 (e.g., Frerking et al. 1982; Leung et al. 1984), and typical dv/dr is estimated for clumps A D as 4kms 1 /0.4pc=10kms 1 pc 1 and 10kms 1 /0.3pc 33kms 1 pc 1, respectively. Figure7 provides a useful guide for the R 3 2/1 0 distributions in Figures6. If R 3 2/1 0 is larger than 0.7, which is depicted by a dashed line in Figure7, T k is always higher than 10K for every n(h 2 ). 10K is a typical value of the galactic molecular gas without any star formation, suggesting that the gas with T k higher than 10K is heated by some additional heating. On the other hand, if R 3 2/1 0 is as small as 0.4, T k cannot be determined, since it has a large variation from <10K to >100K depending on n(h 2 ). Both of the ring cloud and the finger cloud show R 3 2/1 0 higher than 0.7 up to over 1.0 at many points, suggesting that these two clouds have T k higher than 10K. It is reasonable that the heating is due to the O stars in RCW38. On the other hand, the physical condition of the diffuse CO emission surrounding the ring cloud with R 3 2/1 0 of 0.4 is unclear. In order to provide more quantitative constraints for T k and n(h 2 ), we here add the ratio of 13 CO J=1 0 to 12 CO J=1 0 (hereafter R 13/12 ) in the LVG analysis, with an assumption of the abundance ratio [ 12 CO]/[ 13 CO] = 77 (Wilson & Rood 1994). Ten target regions are finally chosen (regions A J) in the ring cloud and the finger cloud for T k and n(h 2 ) estimates. The CO spectra and the LVG results are presented in Figures8 and 9

15 15 for the individual target regions. In the CO spectra, the velocity ranges used for the LVG analysis are shown by shade. For the ring cloud, the five target regions are chosen to have apparent self-absorption in their 12 CO spectra. For the ring cloud, regions A C are taken around the rim of the dense part, while regions D and E are at the low R 3 2/1 0 regions widely distributed in the north and the south. On the other hand, regions F J in Figure9 are taken to cover a large area of the finger cloud. In the diagrams of the LVG results, R 3 2/1 0 and R 13/12 distributions are shown with the blue lines and red lines, respectively, and their errors are shown with the colored area. Here the errors are estimated with 1σ fluctuations of the spectra and 10% relative calibration error for each CO transition. Since 12 CO J=1 0 and 13 CO J=1 0 were taken simultaneously with the same receiver and the same backend at Mopra, we here assume that the 10% relative calibration error is canceled for R 13/12, and is adopted only for R 3 2/1 0. In addition, dv/dr is estimated for individual regions, and X(CO) is assumed to be 10 4 same as in Figure7. Finally, T k and n(h 2 ) are given at the region where R 3 2/1 0 curve and R 13/12 curve overlap. The results are summarized as follows: Regions A C in the ring cloud and all regions in the finger cloud (regions F J) show T k higher than 10K, typically 30 40K, up to more than 50K. Typical n(h 2 ) for these regions are about cm 3, where only lower limits are given for regions A, B, and F2 and no solutions are given for regions F1 and I. On the other hand, regions D and E located at the south of the ring cloud show significantly low T k less than 10K and n(h 2 ) of cm 3. As a summary, both of the ring cloud and the finger cloud show remarkably high temperatures of T k > 20 30K. Some heating source/mechanism is necessary to understand the results, and radiative heating by the O stars in RCW38 is a reasonable candidate. The present results suggest that, despite of a large velocity separation of 12kms 1, both of the ring cloud and the finger cloud are associated with RCW38.

16 Molecular masses of the ring cloud and the finger cloud By adopting the distance of RCW38, 1.7kpc, to both the ring cloud and the finger cloud, their molecular masses can be estimated. We first estimate the molecular mass by using the X(CO) factor, which is an empirical conversion factor from 12 CO J=1 0 integrated intensity into H 2 column density. We hereafter refer to the molecular mass estimated with an X(CO) factor as M X(CO). As already stated, since 12 CO is strongly affected by the self-absorption in particular in the central 1 2pc of RCW38, M X(CO) may be underestimated. It is, however, still quite useful to estimate M X(CO) to give the lower limit of the molecular mass. We here assume an X(CO) factor of cm 2 (Kkms 1 ) 1 (Strong et al. 1988). Using the NANTEN CO dataset, M X(CO) in the central peak of the ring cloud seen in Figure2(a) is estimated to be M, for the points with the 12 CO J=1 0 integrated intensity W( 12 CO1 0) equal to or larger than 65Kkms 1, and M X(CO) of the head-tail structure depicted by the dashed line in Figure2(a) is M at W( 12 CO1 0) 30Kkms 1. On the other hand, M X(CO) of the finger cloud at the region depicted by the dashed lines in Figure2(b) is M, an order of magnitude smaller than that of the ring cloud. M X(CO) of the central peak of the ring cloud is also estimated using the Mopra 12 CO J=1 0 dataset (Figure3) to be M at W( 12 CO1 0) 70Kkms 1, and that in the ring structure in the inner 1 1.5pc is M (W( 12 CO1 0) 140Kkms 1 ). On the other hand, M X(CO) of the ring cloud in Figure3(b) is calculated to be M at W( 12 CO1 0) 10Kkms 1. We next estimate the molecular mass of the ring cloud using the Mopra 13 CO J=1 0 data with the local thermodynamic equilibrium (LTE) assumption. In estimate of the molecular mass with an LTE assumption (hereafter M LTE ), excitation temperature T ex must

17 17 be provided, and the peak intensity of the 12 CO J=1 0 emission T( 12 CO) is usually used to estimate T ex with the following equation, T ex = 5.53 ln{1+5.53/(t( 12 CO)+0.819)}. (1) Although the 12 CO J=1 0 emission strongly suffers from the self-absorption as shown in the spectra for clumps A D, their observed brightness temperature is still remarkably high up to 45K in clump C, which corresponds to T ex of 49K with equation(1). This figure is roughly consistent with the LVG analysis resulting in T k > 30 50K (Figure8) in regions A C which are located in the outskirts of the ring cloud. We therefore adopt an uniform T ex of 49K in the present M LTE estimates. M LTE of the ring-like structure shown in Figure3(d) is finally estimated to be M at the 13 CO J=1 0 integrated intensity W( 13 CO1 0) 22.5Kkms 1, which is by a factor of 2.7 larger than M X(CO) in the ring-like structure. The present estimate also provides a remarkably high H 2 column density N(H 2 ) of cm 2 as a typical value in the ring-like structure. On the other hand, M LTE and typical N(H 2 ) of the central cavity surrounded by the ring-like structure of 0.3pc 0.4pc is estimated to be 120M and cm 2. M LTE of the four 13 CO clumps shown in Figure4(a) is estimated at W( 13 CO1 0) 32.5Kkms 1, and the results are given as follows: M for clump A, M for clump B, M for clump C, and M for clump D. These results are also summarized in Table Comparisons with infrared images Figure 10 shows a comparison between the CO J=1 0 distribution and the Spitzer 3.6 µm image which consists of several filamentary features around the central peak. Only panel (c) is for 13 CO. Figure10(a), a large-scale view, shows that the ring cloud

18 18 and the intense 3.6 µm emission corresponds well to each other, confirming the results by Kaneda et al. (2013). The finger cloud in Figre 10(b) shows possible correspondence with the southern part of the 3.6µm distribution at (R.A., Dec.)=(08 h 59 m 5 s, ). Figure10(c) compares the YSOs and O star candidates in the ring cloud. The YSO distribution corresponds to the overall distribution of the ring cloud, and the O stars with the northwester part of the cavity. Figure10(d) indicates that nearly 20 of the O stars are located toward the tip of the finger cloud overlapping with the bridge feature, while, to be more into detail, the distribution of the O stars shows an offset slightly toward the south of the peak of the bridge feature. The other YSOs do not show strong correlation with the finger cloud in the south (Figure10(c)). Figure 11 shows velocity channel distributions of CO overlayed on the infrared image. We label the nine infrared filamentary features by a i in the first panel of Figure11. The CO features generally show good correspondence with CO as listed in Table 3. Correspondence of the infrared image with the finger cloud is not so obvious in Figure11. Figure12, an overlay of the finger cloud with the VLT image at a smaller scale, however, shows a good agreement with IRS1 which is elongated in the north to south while the angular resolution is limited for spatially resolving the bridging feature. We interpret that the dark lane toward IRS 1 along with another dust feature comprising the western cavity corresponds to the finger cloud which probably lies on the near side of the cluster. IRS1 is clearly irradiated by the cluster and the physical association of the bridging feature with the cluster is supported. To summarize the comparison, we have clear indications of the physical association of the ring cloud and the finger cloud with RCW38.

19 19 4. Discussion 4.1. RCW38 as the youngest super star cluster in the Milky Way The number of known super star clusters in the Milky Way is small, only 20 (see Table 2 in Portegies Zwart et al. 2010). Among them RCW38 is the youngest one with an estimated age of Myrs (Wolk et al. 2006, Section 4.3 of the present work) while the ages of the others are 2 3Myrs. According to eye inspection of the infrared images, only four clusters, NGC 3603, RCW 38, Westerlund 2 and DBS[2003]179, are associated with dust nebulosity. The remaining clusters have no associated nebulosity; the most typical cases are seen in the Arches and Quintuplet clusters where the nebulosity is completely dispersed by the stellar feedback including ionization and stellar winds. Figure 13 shows the averaged CO intensity in the three clusters RCW38, NGC3603, and Westerlund2, indicating that RCW38 has the largest molecular mass within a few pc of the cluster. On the contrary, in Westerlund2 and NGC3603 the central part of the clouds are already strongly ionized/dispersed, where the radii of the cavities are at least a few pc. The rich remaining gas is a unique feature of RCW38 and may provide a best opportunity to study the initial conditions for the super star cluster formation The two clouds; physical association and evidence for collision The present observations have revealed the distribution of the two molecular clouds with 12kms 1 velocity separation toward RCW38 at a spatial resolution of 0.1 pc. The following pieces of evidence offer robust verification of the association between the two clouds and the cluster. 1. The ring cloud shows a central cavity created by the feedback of the cluster. The ring

20 20 cloud is also associated with the infrared filamentary/bubble-like features heated by the cluster. The finger cloud shows a good correspondence with IRS1 and the western cavity. These morphological correspondences support the physical association of the clouds with the cluster. 2. The line intensity ratio of the 12 CO J=3 2 emission to 13 CO J=1 0 emission shows significantly high values, indicating high gas temperatures toward the cluster. LVG calculations show that the temperature there is as high as 30 K, consistent with local heating by the O stars in the cluster. It is also possible that shock heating may provide additional heating, which also supports the physical interaction between the two clouds with the cluster. 3. Another sign of the interaction between the two clouds is seen as the bridging features between them in velocity at least in three places including the direction of the cluster center. The cloud velocity separation is large 12kms 1 as seen in the previous two cases Westerlund2 and NGC3603 (Furukawa et al. 2009; Fukui et al. 2014). The velocity is too large to be gravitationally bound by the total mass of the clouds and cluster, less than 10 5 M. This indicates that the association of the two clouds is by chance; it is possible that the relatively rich molecular gas in the Vela ridge region (Yamaguchi et al. 1999, 2001) favors more frequent encounters than the average over the Galaxy. Accordingly, we frame a hypothesis that the two clouds collided with each other recently and that the collision triggered formation of stars, primarily O stars, as have been suggested for Westerlund 2 and NGC The implications of the bridging feature on cloud-cloud collision are discussed in detail by Haworth et al. (2015). We note that the fourth super star cluster with nebulosity DBS[2003]179 is also associated with two molecular clouds which appear colliding to form the cluster (Kuwahara et al. in preparation).

21 21 Stellar wind acceleration may be an alternative interpretation to cloud-cloud collision. We see velocity gradient in the finger cloud (Section 3), which may suggest some acceleration by the O stars via the stellar winds. This is however not likely in the present case because the gradient is extended outside the cluster and rather uniform over the finger cloud; the velocity shift should be most prominent toward the cluster if the winds play a role. We also note that the finger cloud is probably located on the nearside of the cluster as suggested by the coincidence with IRS 1 (Figure 12), implying that the wind acceleration will cause blue-shift instead of the observed red-shift. In addition, the inferred time scale from O star formation, Myrs, may be too short for the winds to affect the clouds. The collision between the two clouds is therefore plausible Details of the collisional interaction The remaining rich molecular gas in RCW38 allows us a unique opportunity to study the interaction between the clouds in the early stage of the cloud-cloud collision and triggering. In Figure 10(d) we find the distribution of the finger cloud and the bridging feature show a good spatial correlation with the O star candidates. This suggests that the cloud-cloud collision preferentially triggered localized O-star formation but not the extended low-mass YSO s other than the O stars (Figure 1). It is likely that the extended low-mass star formation occurred prior to the cloud-cloud collision in the ring cloud. By assuming that the relative cloud motion is inclined to the line of sight by 45, we estimate the actual relative velocity of the two clouds to be 17kms 1, where the velocity may be decelerated by momentum conversation in collision. The radius of the cluster s influence is seen to be about 1.0pc from the radius of the ring (Figure10). By adopting an age spread of youngest star formation, 0.1 Myrs, derived for the youngest stars in NGC 3603 (Kudryavtseva et al. 2012), we assume that the duration of O star formation is 0.1Myrs.

22 22 On the other hand, the star formation by collision was initiated 0.1Myrs ago as roughly estimated by the ratio [the cloud size]/ [relative velocity]=(1pc)/(10kms 1 ), where we consider the uncertainty due to the projection effect and possible deceleration of the shocked interface layer which makes the shock velocity less than the relative velocity of the two clouds (e.g., see Takahira et al. 2014). We infer that RCW38 consists of the localized young O stars of 0.1Myrs old and the older low-mass YSOs. Here we note that the quoted age 0.5Myr (Wolk et al. 2006) is an assumed value not conflicting with the stellar near infrared colors and does not exclude 0.1Myrs as the stellar age. Duality in stellar age is suggested in RCW38 as found in NGC3603 (Harayama et al. 2008), and R136 in the LMC (Crowther et al. 2010). It is thus plausible that collision between the ring and finger clouds triggered formation of O stars in RCW38. The collision may be still continuing at present. Figure12 shows that parts of the molecular clouds are ionized; the ring cloud has a cavity of 1pc radius where the mass of the ionized gas is roughly estimated to be 500M (see Table 1), and the finger cloud seems to be eroded on the eastern side toward IRS2. A possible scenario is as follows; the collision created shocked interface layer which formed a few tens of proto O stars within 1 pc of the cluster center. MHD numerical simulations of supersonically colliding molecular flows indicate indeed a top-heavy core mass function which favors O star formation (Inoue & Fukui 2013). The protostars grow by mass accretion and the most luminous member becomes bright in the ultraviolet to ionize the accretion disks of the surrounding proto O stars. Then their mass growth is halted and the O stars become exposed as observed now. The existence of IRS 1 with associated molecular gas indicates that there is still molecular gas interacting with the ring cloud. This collisional interaction may lead to formation of more O stars in another 0.1Myrs in the C 18 O clump IV whose mass is 500M.

23 Conditions for O star formation The present clouds show signs of collisional interaction in the three positions in Figure 5. Only one of them shows O star formation. This may give us a hint on conditions for O-star formation. The relative velocities between the two clouds are 12kms 1, 11kms 1 and 7kms 1, respectively, in the three positions from the north to the south. The column densities toward the positions are cm 2, cm 2 and cm 2, respectively in the ring cloud and those in the finger cloud are nearly uniform at cm 2. We infer that the column density of one of the colliding clouds is an important condition, and that the column density around cm 2 may be the threshold value for O star formation in cloud-cloud collision. Alternatively, the three regions may have different time scales of collision, and the other two positions of the interaction may still be in an early stage prior to O star formation. The ring cloud has a total mass of 10 4 M, while the finger cloud has a total mass of 10 3 M. They are smaller than the giant molecular clouds whose mass is M, suggesting that the total cloud mass is not critical for O star formation. These results are compared with the triggering in the other regions of cloud-cloud collision, NGC3603, Westerlund2, M20, and RCW120 (Table 4). NGC3603 has 30 O stars. One of the colliding clouds has a high molecular column density of cm 2 at its peak while the other has a low molecular column density of cm 2. This is a similar situation to the two clouds in RCW38, although the significant ionization toward the cluster makes these column densities less certain than in RCW 38. The high column density cloud has a total mass of 10 5 M, while the other more spatially extended has a total mass of 10 4 M. Westerlund 2 has 14 O/WR stars observationally confirmed. The cluster shows even more significant ionization of the molecular gas and one cannot infer the initial conditions of the colliding spots of the clouds better than in RCW38. The two clouds are giant

24 24 molecular clouds having 10 5 M in each and their average molecular column densities are both equal to cm 2 in the parts outside the cluster. M20 has an O star with four lower-mass stars (two B stars, one A star and one F star) within 0.1pc of the O star in the collisional interaction region. Its age is very young, 0.3Myrs. So this has a similar properties to RCW38 including rich molecular gas close to the O stars. The two clouds have 10 3 M for each with molecular column densities of cm 2. RCW120 has an O star in the Spitzer bubble with an age of 0.5Myrs. There is no direct evidence for additional stars near the O star. The mass of the cloud is M with a molecular column density of cm 2. One may infer the molecular column density of the other cloud to be cm 2 from the peak outside the bubble, where the other cloud is already dissipated by the collisional interaction according to the model by Torii et al. (2015). Although the statistics is not yet significant, it still makes sense to discuss some trend of collision-triggered O star formation for the five regions in the Milky Way. By summarizing the above, we suggest the followings as possible answers to the five question raised in the introduction; Q.I How the mass and number of the formed O stars are determined in the collision; how is the dependence of the stellar properties on the collision parameters such as the relative velocity between the clouds, cloud mass, cloud density, and cloud morphology. A.I Formation of O stars requires high molecular column density of cm 2 in one of the colliding clouds, whereas a single O star can be formed at lower molecular column density of a few times cm 2. It is not necessary that both of the two clouds have the high column density. For a relative velocity of 20 30kms 1, we see

25 25 no strong dependence of O star formation on the relative velocity, whereas it seems that supersonic collision at 20kms 1 or higher is required for O-star formation. We suggest in RCW38 that most of the O stars are formed by triggering in the supersonic cloud-cloud collision, although we see some non-clustered O-star candidates which are located outside the collisional-interaction regions. It is not excluded that they were formed by collisional triggering with lower-intensity/smaller clouds not detected with the present beam size. Q. II What is the origin of the mass segregation/stellar distribution in super star clusters/mini-starbursts? A. II We infer that the O star distribution holds the distribution of the collisional interaction at least right after the formation. This is best indicated by the present work on RCW38. Q. III What are the density, velocity field and magnetic field in the shocked interface layer. This is most relevant to our understanding of details of high-mass protostellatr formation/evolution. We note that in the above four objects in the Miky Way molecular gas is already ionized in the central pc of the O star(s), hindering one to see directly the initial conditions of the natal gas. A.III We have identified the possible shocked interface in RCW38 and M20 at 0.1pc resolution. Observations at higher resolutions with ALMA are required for these youngest objects to obtain details including the cloud core mass function down to 0.01M in core mass which is predicted to be top-heavy by numerical simulations (Inoue & Fukui 2013). Q. IV Can cloud-cloud collision also trigger formation of low- to intermediate-mass stars? If yes, what is the mass function of the formed stars.

26 26 A.IV In the three objects except for M20, it seems difficult to identify the trigger-formed member lower-mass stars due to confusion. We however have strong sings for the significantly top-heavy, nearly flat, mass function toward the O stars in NGC 3603 (Harayama et al. 2008) and R136 in the LMC (Crowther et al. 2010). Collisiontriggered formation of lower-mass stars is also seen in the four stars (two B, one A and one F) in M20, certainly showing a top-heavy stellar mass function. Q.V What is the fraction of O stars formed by cloud-cloud collisions in the Milky Way? Are cloud-cloud collisions one of the major modes of O-star formation? A.V We do not yet have an answer on the question Q.V in the present work. The answer remains to be explored in future works A comparison of O-star formation under the cloud-cloud collision In the cloud-cloud collision the timescale of the collision is estimated from the relative velocity and the size of the system. Also the mass of the formed stars is estimated by near infrared observations. In RCW38 the O5.5 star IRS2 has mass of 40M, and the mass accretion rate required is roughly estimated to be M yr 1 for 0.1Myrs, if the accretion rate is assumed to be constant. Similar high mass accretion rates and short timescales are suggested in the other four regions of cloud-cloud collision. It is shown that the age spreads of young stars are remarkably short, in the order of 0.1Myrs, and that the star formation took place a few Myrs ago in the two clusters NGC3603 and Westerlund 1 by careful measurements of stellar ages with HST and VLT (Kudryavtseva et al. 2012). These observations are particularly relevant in the context of cloud-cloud collision; a possible interpretation is that the age, a few Myrs, corresponds to the timing of the collision, and the age spread, 0.1 Myrs, corresponds to the duration of the shock interaction which will

27 27 be terminated when the O star began to ionize the gas disks of the individual protostars. Since the highest-mass stars have M in these clusters, the mass accretion rate is estimated to be M yr 1 to 10 3 M yr 1 for 0.1Myrs as is consistent with RCW 38. We next discuss the mass range of stars formed by cloud-cloud collision. Zinnecker & Yorke (2007) gave the following discussion; Massive star formation is not a simple scaled-up version of low-mass star formation, particularly when it comes to ZAMS O stars (> 20M ). The formation of early-type B stars (10M ), however, may still be considered as a continuation of low-mass star formation. The difference between a 10M star and a 20M star is often not clearly discerned in the discussion of high-mass star formation, and both of them are treated as a single category high-mass stars. It is however important to appreciate their difference in the physical properties as discussed by Zinnecker & Yorke (2007); for instance, early B stars cannot ionize the accreting gas, whereas O stars ionize the surroundings. This naturally causes significant difference in their formation mechanisms, too, in the accretion processes. As a natural consequence of the observational bias toward super star clusters and Hii regions, the cloud-cloud collision found so far is associated with O stars having mass larger than 20M. This seems complimentary to the infrared dark clouds (e.g., Perault et al. 1996; Sridharan et al. 2005) and hot molecular cloud cores (e.g., Kurtz et al. 2000; Menten et al. 2005), where Hii regions are usually not seen. This may suggest that the two classes of dense clouds form mainly early B stars but not O stars. The supersonic flow caused by the collision allows strong non-gravitational accumulation of the dense gas into a small volume very quickly within a typical timescale of 0.1Myrs. This becomes efficient by collimation of the colliding flow into a small volume of a 1pc 3 scale via the guided flow in the curved magnetic field (Inoue & Fukui 2013) or by non-magnetic

28 28 focusing (Takahira et al. 2014). The high-mass star formation by cloud-cloud collision is different from the monolithic collapse because in the cloud-cloud collision scenario dense and massive cores are not required as the initial condition. The large relative motion allows the system to collect a large mass in the shock compressed layer in a timescale of 10 4 yrs even for a low density initial condition of molecular gas. The case of N159 West in the Large Magellanic Cloud indicates that a dense CS/dust clumps is seen only toward the 40M star, and the rest of the filamentary clouds has no dense CS/dust clumps, suggesting that the initial condition is of low density, less than 10 4 cm 3, even for O-star formation (Fukui et al. 2015). The Large Magellanic Cloud allows us to have a 10pc scale view at sub-pc resolution with little contamination by using ALMA, and a large sample of O star formation will help an better understanding of O star formation. 5. Conclusions We have carried out CO J=1 0 and J=3 2 observations toward the super star cluster RCW 38 with Mopra, ASTE and NANTEN2 mm/sub-mm telescopes. The main conclusions of the present study are summarized as follows; 1. We have observed two molecular clouds at velocities of 2kms 1 and 14kms 1 toward RCW38. The ring cloud (the 2kms 1 cloud) shows a ring-like shape with a cavity ionized by the cluster and has a high molecular column density of cm 2. The other the finger cloud (the 14kms 1 cloud) has a tip toward the cluster and has a lower molecular column density of cm 2. The total masses of the two clouds and that of the ring cloud are M and M, respectively. 2. It is likely that the two clouds are physically associated with the cluster as verified by the high ratio of the J=3 2 transition to the J=1 0 transition, R 3 2/1 0, toward

29 29 the cluster. The observed ratio indicates kinetic temperature of 30 K according to the LVG calculations, which is significantly higher than the usual molecular cloud temperature 10 K with no extra heat source. The association is further supported by the similar distribution of the molecular clouds with the infrared dust features heated by the cluster; the ring cloud shows remarkable correlation with the extend dust features, and the finger cloud corresponds to the infrared ridge IRS 1 toward the cluster. In addition, the two clouds are linked with each other by bridging features in velocity at least in three places including the direction of the cluster, suggesting the physical connection between the two clouds. 3. Since the total mass of the clouds and the cluster, less than 10 5 M, is too small to gravitationally bind the velocity, we suggest that the clouds encountered by chance. We present an interpretation that the two clouds collided 0.1Myrs ago at 10 kms 1 with each other to trigger formation of the high-mass stars in the cluster. The 20 O star candidates were formed in this timescale with the mass accretion rate M yr 1 for an assumed average stellar mass 20M. 4. The tip of the finger cloud and the bridging features connecting the two clouds well coincide with the youngest O stars within 0.3 pc of the cluster center, indicating that the triggering happened only toward the inner 0.3 pc and the other member low-mass stars outside 0.3 pc are extended over the ring cloud where collisional interaction is not taking place. RCW38 is the third super star cluster alongside of Westerlund2 and NGC 3603 where cloud-cloud collisions triggered the O star formation in the cluster, lending further support for an important role of supersonic collision in O star formation. Among the three super star clusters, RCW 38 is unique because it is the youngest cluster where the initial conditions prior to the O-star formation still hold without significant cloud dispersal by ionization.

30 30 5. The present findings suggest a possibility to write a recipe for O star formation; O star formation takes place where collision of two clouds at velocity 10 30kms 1 happens. One of the clouds should have a high column density cm 2, while the other may have lower column density. Since O star formation takes place only at the spots of the collision, the distribution of O stars, outputs of the collision, keeps the distribution of the collision as long as the rearrangement of O stars by the mutual gravity after the collision does not play a role. This picture provides a natural set up for the highly turbulent interstellar medium as the initial condition for high-mass star formation along the line of suggestion by Zinnecker & Yorke (2007) that supersonic gas motions accumulating gas set up the initial conditions for O star formation and that formation of O stars (ZAMS 20M ) is not a simple scale-up version of lower-mass stars (early B stars with 10M ). In order to better establish the crucial role of cloud-cloud collision in O star formation, we need a systematic study molecular clouds toward a number of young Hii regions and very young stars whose mass exceeds 20M. Such a study is on-going and will allow us to test how cloud-cloud collision dominates/regulates O star formation. It is also a challenge to test if mini-starbursts in galaxies is possible by supersonic collisions between molecular clouds. NANTEN2 is an international collaboration of ten universities, Nagoya University, Osaka Prefecture University, University of Cologne, University of Bonn, Seoul National University, University of Chile, University of New South Wales, Macquarie University, University of Sydney and Zurich Technical University. The ASTE telescope is operated by National Astronomical Observatory of Japan (NAOJ).

31 31 REFERENCES Anathpindika, S. V. 2010, MNRAS, 405, 1431 Bonnell, I. A., Vine, S. G., & Bate, M. R. 2004, MNRAS, 349, 735 Crowther, P. A., Schnurr, O., Hirschi, R., et al. 2010, MNRAS, 408, 731 DeRose, K. L., Bourke, T. L., Gutermuth, R. A., et al. 2009, AJ, 138, 33 Dobbs, C. L., Pringle, J. E., & Duarte-Cabral, A. 2015, MNRAS, 446, 3608 Ezawa, H., Kawabe, R., Kohno, K., & Yamamoto, S. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 5489, Ground-based Telescopes, ed. J. M. Oschmann, Jr., Ezawa, H., Kohno, K., Kawabe, R., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7012, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 8 Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590 Frogel, J. A., & Persson, S. E. 1974, ApJ, 192, 351 Fukui, Y., & Kawamura, A. 2010, ARA&A, 48, 547 Fukui, Y., Ohama, A., Hanaoka, N., et al. 2014, ApJ, 780, 36 Fukui, Y., Harada, R., Tokuda, K., et al. 2015, ArXiv e-prints, arxiv: , accepted for publication in ApJL Furniss, I., Jennings, R. E., & Moorwood, A. F. M. 1975, ApJ, 202, 400 Furukawa, N., Dawson, J. R., Ohama, A., et al. 2009, ApJ, 696, L115

32 32 Goldreich, P., & Kwan, J. 1974, ApJ, 189, 441 Gyulbudaghian, A. L., & May, J. 2008, Astrophysics, 51, 18 Habe, A., & Ohta, K. 1992, PASJ, 44, 203 Harayama, Y., Eisenhauer, F., & Martins, F. 2008, ApJ, 675, 1319 Haworth, T. J., Tasker, E. J., Fukui, Y., et al. 2015, MNRAS, 450, 10 Huchtmeier, W. 1974, A&A, 32, 335 Inoue, H., Muraoka, K., Sakai, T., et al. 2008, in Ninteenth International Symposium on Space Terahertz Technology, ed. W. Wild, 281 Inoue, T., & Fukui, Y. 2013, ApJ, 774, L31 Kaneda, H., Nakagawa, T., Ghosh, S. K., et al. 2013, A&A, 556, A92 Kudryavtseva, N., Brandner, W., Gennaro, M., et al. 2012, ApJ, 750, L44 Kurtz, S., Cesaroni, R., Churchwell, E., Hofner, P., & Walmsley, C. M. 2000, Protostars and Planets IV, 299 Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 Ladd, N., Purcell, C., Wong, T., & Robertson, S. 2005, PASA, 22, 62 Leung, C. M., Herbst, E., & Huebner, W. F. 1984, ApJS, 56, 231 McKee, C. F., & Tan, J. C. 2003, ApJ, 585, 850 Menten, K. M., Pillai, T., & Wyrowski, F. 2005, in IAU Symposium, Vol. 227, Massive Star Birth: A Crossroads of Astrophysics, ed. R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley, 23 34

33 33 Ohama, A., Dawson, J. R., Furukawa, N., et al. 2010, ApJ, 709, 975 Perault, M., Omont, A., Simon, G., et al. 1996, A&A, 315, L165 Portegies Zwart, S. F., McMillan, S. L. W., & Gieles, M. 2010, ARA&A, 48, 431 Rodgers, A. W., Campbell, C. T., & Whiteoak, J. B. 1960, MNRAS, 121, 103 Smith, C. H., Bourke, T. L., Wright, C. M., et al. 1999, MNRAS, 303, 367 Sorai, K., Sunada, K., Okumura, S. K., et al. 2000, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 4015, Radio Telescopes, ed. H. R. Butcher, Sridharan, T. K., Williams, S. J., & Fuller, G. A. 2005, ApJ, 631, L73 Strong, A. W., Bloemen, J. B. G. M., Dame, T. M., et al. 1988, A&A, 207, 1 Takahira, K., Tasker, E. J., & Habe, A. 2014, ApJ, 792, 63 Tan, J. C., Beltrán, M. T., Caselli, P., et al. 2014, Protostars and Planets VI, 149 Tasker, E. J., & Tan, J. C. 2009, ApJ, 700, 358 Torii, K., Enokiya, R., Sano, H., et al. 2011, ApJ, 738, 46 Torii, K., Hasegawa, K., Hattori, Y., et al. 2015, ApJ, 806, 7 Vigil, M. 2004, Master s thesis, Massachusetts Institute of Technology Wilson, T. L., & Rood, R. 1994, ARA&A, 32, 191 Winston, E., Wolk, S. J., Bourke, T. L., et al. 2011, ApJ, 743, 166 Wolk, S. J., Bourke, T. L., & Vigil, M. 2008, The Embedded Massive Star Forming Region RCW 38, ed. B. Reipurth, 124

34 34 Wolk, S. J., Spitzbart, B. D., Bourke, T. L., & Alves, J. 2006, AJ, 132, 1100 Yamaguchi, R., Saito, H., Mizuno, N., et al. 1999, PASJ, 51, 791 Yamaguchi, R., Mizuno, N., Mizuno, A., et al. 2001, PASJ, 53, 985 Yorke, H. W., & Sonnhalter, C. 2002, ApJ, 569, 846 Zinchenko, I., Mattila, K., & Toriseva, M. 1995, A&AS, 111, 95 Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481 This manuscript was prepared with the AAS L A TEX macros v5.2.

35 35 Fig. 1. (a) The 3.6µm image of RCW38 obtained with the Spitzer/IRAC observations (Wolk et al. 2008). YSOs and O star candidates obtained with Wolk et al. (2006) and Winston et al. (2011) are plotted in (b). Red circles and white circles indicate the class 0/I and flat spectrum YSOs, respectively, while black dots indicate the class II and III YSOs. White crosses indicate the candidate O stars. (b) A close-up of the central region of RCW38 from the VLT observations (Wolk et al. 2006). Z band is shown in blue, H band is green, and K band is red. The bright infrared emission IRS1 and IRS2 are shown by arrows.

36 36 Fig CO J=1 0 integrated intensity of the red cloud (a) and the blue cloud (b). The observed area with ASTE is shown by black lines in (b). (c) Latitude-velocity diagram of the 12 CO J=1 0 emission. Dashed line indicates the position of RCW38 IRS2.

37 37 Fig. 3. Integrated intensity maps of the blue cloud (a and d) and the red cloud (b and e) and position-velocity maps (c and f) are shown. Upper three panels show 12 CO J=1 0, and lower three 13 CO J=1 0.

38 38 Fig. 4. Enlarged integrated intensity maps of the blue cloud in 13 CO J=1 0 (a) and C 18 O J=1 0 for a velocity range of kms 1.

39 39 Table 1: Parameters of clumps A D Name R.A. Dec. v LSR v size M LTE [ hms ] [ ] [kms 1 ] [kms 1 ] [pc] [M ] (1) (2) (3) (4) (5) (6) (7) clump I 8:59:9.2 47:29: clump II 8:59: :30: clump III 8:59:8.4 47:31: clump IV 8:59:1.2 47:30: Note. Column (1): Name. (2, 3): Position in J2000. (4, 5): Peak v LSR and velocity width v (FWHM) of the 13 CO J=1 0 profile, where fitting with a Gaussian function is used to derive the parameters. (6): Size of the 13 CO clump.

40 40 Fig. 5. A comparison between the blue cloud (color) and the red cloud (contours) in 12 CO J=1 0.

41 41 Fig. 6. The position-velocity map of the 12 CO J=3 2/J=1 0 ratio. The integration range in the R.A. is the same as Figure3(c). The dashed line indicates the position of IRS2.

42 42 X/dvdr = (solid line), (dashed line) 1.0 n=10 4 cm n= cm 3 R 3-2/ n=10 3 cm 3 n= cm 3 n=10 2 cm T [K] Fig. 7. Curves of the 12 CO J=3 2/J=1 0 ratio as a function of T k for various n(h 2 ) ranges estimated with LVG. Here the solid lines and dashed lines indicate X(CO)/(dv/dr) of 10 5 (kms 1 pc 1 ) 1 and (kms 1 pc 1 ) 1, respectively.

43 spectra are multiplied by three. 43 Fig. 8. LVG results on the n(h 2 )-T k plane for the regions A E in the blue cloud are shown. 12 CO J=3 2/J=1 0 is plotted in black, and 13 CO/ 12 CO J=1 0 is in red. The regions A E are plotted with black crosses in the CO integrated intensity map, where 12 CO J=3 2 is shown in color, while 13 CO J=1 0 is shown with contours. Spectra of 12 CO J=1 0 (black), 12 CO J=3 2 (black), and 12 CO J=1 0 (blue) are also shown. Intensities of the 13 CO J=1 0

44 44 Fig. 9. LVG results for the regions F J in the red cloud in the same manner as Figure 14.

45 45 Fig. 10. Comparisons between the two molecular clouds (contours) and the Spitzer 3.6 µm image. The velocity ranges for the left two panels (a and b) and for the right two panels (c and d) are kms 1 and kms 1, respectively. White contours in (d) indicate the bridge feature integrated over a velocity range kms 1. YSOs and O star candidates are plotted with the same manner as in Figure1.

46 46 Fig. 11. Velocity channel maps of the 12 CO J=3 2 emission (contours) superimposed on the Spitzer 3.6 µm image. In the left-top panel, approximate distributions of the infrared filaments(filaments a h) are depicted by the dashed lines. The parenthetic numbers indicate peak velocities of the filaments in kms 1.

47 47 Table 2: CO velocity range corresponding to the infrared filaments Name v LSR range [kms 1 ] (1) (2) a b c d e f g h i Note. Column: (1) Name of the filament. (2) The velocity range in which CO emission corresponds to the filament.

48 48 Fig CO J=3 2 distributions of the red cloud ( kms 1 ) superimposed on the VLT image.