arXiv:1803.07237v1 [astro-ph.GA] 20 Mar 2018 smvn rmGlci ett atwt epc oteitrcigm interacting the to innermost respect the with along de east moving longitude to is Galactic west cluster the Galactic as from increases moving shells is the counte of infrared an time of expansion absence The the although cluster, star massive nryo 10 of f10 of n ae lr awl eecp.Tekntceeg n expansio and energy kinetic The 10 Telescope. be Clerk James and J it.Ti ae eot h eeto ffieepnigsel nth in shells expanding five of detection the reports paper This width. 1 2 10 CO =1–0, aaaa2382,Japan 223-8522, Kanagawa colo udmna cec n ehooy rdaeSch Graduate Technology, and Science Fundamental of School eateto hsc,Isiueo cec n Technology and Science of Institute Physics, of Department hh Tsujimoto, Shiho yee sn L using 2018 Typeset 21, March version Draft The − 48 5 l . –10 = 3 Keywords: – − 50 − 1 . 4 8 51 J ◦ . eini h aatccne a ihCO high a has center Galactic the in region 2 32 n SiO and =3–2, EETO FAOHRMLCLRBBL NTEGLCI CENTE GALACTIC THE IN BUBBLE MOLECULAR ANOTHER OF DETECTION yr r n 10 and erg r n h yia xaso ieof time expansion typical the and erg A T − E 1 X hsidctsta the that indicates This . aay etr—IM lus—IM molecules ISM: — clouds ISM: — center Galaxy: 1 twocolumn ooauOka, Tomoharu 4 . 7 – 5 . J 0 Rcie a 9 07 eie;Acpe ac 1 2018) 21, March Accepted Revised; 2017; 29, May (Received 87ln aast bandwt h oeaaRdoOsraoy4 Observatory Radio Nobeyama the with obtained sets data line =8–7 tl nAASTeX61 in style r epciey h rgno hs xadn hlsi icse.T discussed. is shells expanding these of origin The respectively. yr, ,2 1, x 1 hnaTakekawa, Shunya ri n h neatn a oldswt tfo h aatceastern Galactic the from it with collides gas interacting the and orbit l = − n utnA Roll A. Justin and 1 umte oApJ to Submitted ◦ . einmyb oeua ubeascae iha embedded an with associated bubble molecular a be may region 2 ABSTRACT eoUiest,31- ioh,Ykhm,Kngw 223- Kanagawa Yokohama, Hiyoshi, 3-14-1 University, Keio , ∼ o fSineadTcnlg,Ki nvriy -41Hiyo 3-14-1 University, Keio Technology, and Science of ool 10 5 J 1 rcrepn omlil uenv xlsosa rate a at explosions supernova multiple to correspond yr =3–2/ aaaYamada, Masaya pr ae hsitrrtto oehtcontroversial. somewhat interpretation this makes rpart J 1 10itniyrtoadeteeybodvelocity broad extremely and ratio intensity =1–0 e rae,sgetn httemsiesa cluster star massive the that suggesting creases, ieo h xadn hlsaeetmtdto estimated are shells expanding the of time n lclrgs epooeamdlweenthe wherein model a propose We gas. olecular l = − 1 ◦ . einbsdo h CO the on based region 2 1 eioTokuyama, Sekito 1 ue Iwata, Yuhei R ettlkinetic total he telescope m 5 J =1–0, 52 Japan 8522, h,Yokohama, shi, 13 side. CO 1 2 Tsujimoto et al.

1. INTRODUCTION out this paper, we adopt 8.3 kpc as the distance to the The central molecular zone (CMZ), the region within Galactic center. ∼ 200 pc from the center of our Galaxy, is character- 2. DATA 4 −3 ized by a large amount of dense [n (H2)≥10 cm ] and The mapping data were obtained by using the NRO warm (Tk =30–60 K) molecular gas (Morris et al. 1983; Paglione et al. 1998), having highly turbulent and com- 45 m radio telescope and the Tele- scope (JCMT). The observed lines are 12CO J=1–0, plex kinematics, with a large velocity width (∆V ≥ 20 13 12 − J J J km s 1). Our group surveyed the CMZ in the CO J=1– CO =1–0, CO =3–2, and SiO =8–7. ◦ ◦ ◦ 0 line (115 GHz) covering −1.5 ≤ l ≤ +3.4, |b| ≤ 0.6 2.1. NRO 45 m Data by using the Nobeyama Radio Observatory (NRO) 45 The on-the-fly (OTF) mapping observation of 12CO m telescope (Oka et al. 1998) and the CO J=3–2 line J=1–0 (115.27120 GHz) and 13CO J=1–0 (110.20135 (346 GHz) covering −1◦.5 ≤ l ≤ +1◦.0, |b|≤ 0◦.2 by us- GHz) lines was performed using the NRO 45 m ra- ing the Atacama Submillimeter Telescope Experiment dio telescope in 2016 January–March and May as (ASTE: Oka et al. 2007, 2012). By carefully inspecting parts of the large-scale surveys of the Galactic CMZ these data sets, we found four high CO J=3–2/J=1–0 (Tokuyama et al. in preparation). To investigate the intensity ratio (R3–2/1–0 ≥ 1.5) regions, at the Galac- ◦ ′ ′ l =−1.2 region, we mainly analyzed a 9.6 × 9.6 area tic longitudes; l = +1◦.3, 0◦, –0◦. 4, and–1◦.2 (Figure 1). (−1◦.32 ≤ l ≤ −1◦.16 and −0◦.20 ≤ b ≤ −0◦.04). We All of the high-ratio regions show extraordinary broad employed the FOREST receiver operated in the two- velocity widths (∆V ≥ 50 km s−1), containing several sideband mode. At 110 GHz, the half-power beamwidth high-velocity compact clouds (HVCCs), which are a (HPBW) of the telescope was 14′′ and the main beam peculiar population of molecular clouds with compact efficiency (η ) was 0.4. We used the SAM45 spec- appearances (d < 10 pc) and broad velocity widths MB −1 trometers in the 1000 MHz bandwidth (244.14 kHz (∆VLSR ≥ 50 kms ) identified in the CMZ (Oka et al. resolution) mode. The system noise temperature (Tsys) 1999, 2001, 2007; Tanaka et al. 2007; Oka et al. 2012). 13 ◦ ranged from 300 to 800 K during the CO observations, The l =+1.3 region, one of the high R3–2/1–0 regions, and from 1000 to 2000 K during the 12CO observations. contains nine expanding shells, which typically have very Antenna temperatures were calibrated using the stan- high kinetic energy reaching E ∼1053 erg, and an ex- kin dard chopper-wheel method to obtain the T ∗ scale. We pansion time t ∼ 105 yr. Most importantly, two of A exp scaled T ∗ by multiplying with 1/η to obtain the the expanding shells are associated with SiO clumps at A MB main-beam temperature (T ). All spectra were ob- their high-velocity ends. As gas-phase SiO is a well- MB tained through position switching between target scans established probe of strong interstellar shocks, a series and the clean reference position (l, b)=(+2◦, +1◦). The of supernova (SN) explosions may be responsible for the pointing errors were corrected every 1.5 hr by observing formation of the expanding shells. The kinetic energy the SiO maser source VX Sgr at 43 GHz. The point- requires tens or hundreds of SN explosions. Combined ing accuracy was better than 3′′ in both azimuth and with the typical expansion time of the shells, 105 yr, the elevation. All the NRO 45 m data were reduced using SN rate should be 10−4–10−3 yr−1. Thus we suggest the NOSTAR1 reduction package. We subtracted the that the l =+1◦.3 region is in an early stage of super- 5 6 baselines of all spectra through linear fitting. bubble formation and that a very massive (∼10 – M ⊙) star cluster is embedded in it (Tanaka et al. 2007). 2.2. JCMT Data Another high R region, that is, the l =−1◦.2 re- 3–2/1–0 The 12CO J=3–2 line (345.79599 GHz) and the SiO gion, shares the principal properties with the l =+1◦.3 J=8–7 line (347.33058 GHz) data were obtained by the region. The locations of these two high-ratio regions JCMT. The 12CO J=3–2 line observations were per- in the sky plane are roughly symmetric with respect to formed on 2013 September by the JCMT Plane Survey the Galactic center. The l =−1◦.2 region contains an ex- (JPS) team (Parsons et al. 2017) as part of the Galac- panding shell which has a size of ∼ 10 pc, showing an tic plane survey. To investigate the l =−1◦.2 region, we extraordinary broad velocity width (∆V ≥ 100 km s−1) mainly analyzed a 9.′6×9.′6 area (−1◦.32≤l ≤−1◦.16 and (Oka et al. 2012). In this paper, we report the detection −0◦.20≤b≤−0◦.04). of four additional expanding shells in the l =−1◦.2 region We also observed the SiO J=8–7 line on 2017 June. based on the newly obtained high-resolution CO survey We seta 9.′6×9.′6 area centered at (l, b)=(−1◦.24, −0◦.12) data. The l =−1◦.2 region is another candidate for a molecular bubble in the CMZ of our Galaxy. Through- 1 http://www.nro.nao.ac.jp/~nro45mrt/html/obs/otf/export-e.html Another Bubble in the Central Molecular Zone of Our Galaxy 3

0.5

2.5 3.54.5 5.5 6.57.5 8.5 L=0° 0.0 6.0 9.0 12.0 15.0

L L =+1.3° L =–1.2° –0.5 =–0.4° CO J=3–2/ J=1–0

GalacticLatitude [degrees] 2.0 1.0 0.0 –1.0 Galactic Longitude [degrees]

Figure 1. Composite CO image of the Galactic CMZ. Red and green indicate the velocity-integrated intensities of the CO J=1–0 and CO J=3–2 lines, respectively and blue is the CO J=3–2 emission integrated for data with high CO J=3–2/J=1–0 intensity ratios (R3–2/1–0 ≥ 1.5). The two zoomed-up maps show the velocity-integrated CO J=3–2 line emission around the l =+1◦. 3 and l =−1◦. 2 regions. as the mapping area, and (l, b)=(−1◦.24, −0◦.62) as the most 19 K at (l, b) = (−1◦.20, −0◦.11), while the typical clean reference position. The pointing errors were cor- intensity of the shell is 5–10 K (Figure 2). rected by observing the standard pointing sources of the In the velocity-integrated maps, especially in the CO JCMT. J=3–2 map, we also noticed four small shells inside We used a single-sideband receiver, the Heterodyne or at the edge of the LS. We named them the “small Array Receiver Program (HARP; Buckle et al. 2009). shells (S1–4)”, which are traced by the white dashed- At 345 GHz, the HPBW of the telescope is 14′′, and the and-dotted ellipses, while the LS is traced by the red- ηMB is 0.64. The ACSIS spectrometer was operated in dashed ellipse (Figure 2). The most intense position of the 1000 MHz bandwidth (976.56 kHz resolution) mode. the CO lines is at a location where the LS and S1 over- The Tsys ranged from 100 to 200 K during the HARP lap, and corresponds to the only position where the SiO observations. The JCMT data were reduced with the J=8–7 line was detected. The centers of these shells are Starlink2software package. aligned roughly parallel to the Galactic plane, ranging over 11.6 pc in projected distance. The positions and 3. RESULTS sizes of these shells (l, b, ∆l, ∆b) are summarized in Ta- 3.1. Spatial Distribution ble 1. The top four panels of Figure 2 show the velocity- 3.2. Kinematics integrated l-b maps of CO J=1–0, CO J=3–2, 13CO J=1–0, and SiO J=8–7 line emissions. The veloc- The bottom four panels of Figure 2 show the ity range for integration is from VLSR = −130 to −60 longitude–velocity (l-V) maps of the observed lines in- km s−1. First, all l-b maps of CO lines show an ellip- tegrated over latitude. The latitude range for the inte- ◦ ◦ tical morphology with a less intense central cavity; we gration was from b = −0.14 to −0.10. The straight −1 named this as the “large shell (LS)”. The major axis horizontal lines at VLSR ≃ −60km s on the l-V of the LS is parallel to the Galactic plane. Velocity- maps are absorption features due to the spiral arm integrated intensities are asymmetric, that is, they are in the Galactic disk (3 kpc arm). The emission at the −1 very intense at the Galactic eastern side, and fainter at VLSR ≃ −200km s represents the expanding molec- the western side. These are several intense clumps, in ular ring (EMR; Kaifu et al. 1972; Scoville 1972). In which the peak intensity of CO J=3–2 line reaches al- addition, another wavy emission feature appears at −1 VLSR ≃−150 km s . The LS and S1–4 can be clearly traced as parts of 2 http://starlink.eao.hawaii.edu/starlink ellipses in these l-V maps, especially in the CO J=1–0 4 Tsujimoto et al.

CO J = 3–2 CO J = 1–0 13 CO J = 1–0 SiO J = 8–7

200 300400 500 600 300 400 500 600 700 40 60 80 100 0.5 1.0 1.5 2.0

LS S2 S4

–0.12 –0.08 –0.04 S1 S3 –0.16 GalacticLatitude [deg] –0.20 0 0.1 0.2 0.3 0.4 0.5 0.20.3 0.4 3.0 ×10 –2 5.0 ×10 –2 7.0 ×10 –2 0.5 ×10 –3 1.0 ×10 –3 ] –1

S4

[km s S3 S2 LSR

V S1 –150 –100 –50 LS –200 –1.16 –1.20 –1.24 –1.28 –1.32 Galactic Longitude [deg]

−1 Figure 2. The top four maps are the velocity-integrated (VLSR = −130 to −60 km s ) l-b maps and the bottom four are the l-V maps integrated over latitudes b = −0◦. 14 to − 0◦. 10. From left to right, we see the CO J=3–2, CO J=1–0, 13CO J=1–0 and SiO J=8–7 maps. The red-dashed ellipse shows the main large shell (LS), the white dashed-and-dotted ellipses show the small shells (S1–4).

Table 1. Properties of Expanding Shells

Name l b ∆l ∆b Vexp M log10 Ekin texp log10 Pkin −1 3 4 −1 (degree) (degree) (pc) (pc) (km s ) (10 M ⊙) (erg) (10 year) (erg s ) LS −1.24 −0.12 13 9.1 60 19.0 50.8 8.6 38.4 S1 −1.21 −0.12 4.3 5.8 40 3.81 49.8 5.9 37.5 S2 −1.24 −0.12 5.6 6.0 35 2.06 49.4 7.8 37.0 S3 −1.26 −0.13 3.6 4.4 20 1.49 48.8 9.5 36.3 S4 −1.29 −0.12 4.0 3.1 15 0.99 48.3 11.2 35.8 Total 27.3 50.9 38.5 and J=3–2 maps. The elliptical shapes in the position– 4. DISCUSSION velocity maps indicate expanding motion. The expan- 4.1. Driving Mechanisms sion velocities (Vexp) of the small shells are slower than those of the LS. The expansion velocities of the shells The l-b and l-V maps show incomplete elliptical shapes are listed in Table 1. Line intensities are enhanced at with broad velocity widths (Figure 2). Several processes are proposed to explain the morphology and kinemat- locations where the shells overlap. The LS/S1 overlap- ◦ ping area shows especially intense emission with clumpy ics of the l =−1.2 region; (1) cloud-cloud collision, (2) morphology in the l-b space. This corresponds to the magnetically floated loop, and (3) multiple supernova Galactic eastern edge of the LS and the S1, and the explosions. high-velocity ends of the S1. Another Bubble in the Central Molecular Zone of Our Galaxy 5

Cloud-Cloud Collision —A collision of clouds may gen- belong to shells redundantly. We suppose that all gas in erate a broad velocity width and elliptical spatial mor- the l =−1◦.2 region belongs to the five shells. The gas phology (Habe & Ohta 1992). In this case, position- column density is calculated from the CO J=1–0 data velocity maps should show less intense emission bridges by assuming a local thermodynamic equilibrium (LTE). which connect two intense components with different The optical depth (τ =2.5) and the excitation tempera- velocities (Matsumura et al. 2012; Haworth et al. 2015; ture (Tex = 10 K) are determined based on the CO J=1– Torii et al. 2017). The l =−1◦.2 region has two velocity 0 and 13CO J=1–0 data. We adopt [CO]/[13CO] = 24, −1 −1 4 components at VLSR ∼ −50km s and −150km s , [H2]/[CO]=10 (Langer & Penzias 1990). ◦ and a broad velocity width feature at l ∼−1.20 seems The kinetic energy (Ekin) of each shell is calculated to connect these components. Cloud-cloud collisions by may frequently occur at the interface between x1 and 1 2 x2 orbits, as well as at the inner ends of the dust lanes E = MV . (1) kin 2 exp where the ‘spray’ of incoming gas collides with the inner- most x1 orbit (Binney et al. 1991). These collisions defi- The expansion time of a shell is derived according to its nitely generate broad velocity width features in position- size and expansion velocity. By considering the asym- velocity maps (e.g., Liszt 2006). metric shape of the shells, we define the size of a shell as Magnetically Floated Loop —We may expect broad veloc- ity width features at two footpoints of the molecular ∆l ∆b R = , (2) loops. There are two molecular loops, loop 1 and loop r 2 2 2, discovered in the Galactic center region (Fukui et al. and thus t is calculated as 2006). They are interpreted as formed by magnetic exp floatation caused by Parker instability. Large velocity R texp = . (3) dispersion due to the collision of gas at each footpoint is Vexp expected by the magnetic loop model. We note that the ◦ western edge of LS is at l ∼−1.3 and the eastern foot- The calculated physical parameters, M, Ekin, texp, and point of the loop 1 is at l ∼−2◦, indicating no physical kinetic power (Pkin = Ekin/texp) are listed in Table 1. connection between them, and no molecular loop struc- Such a large kinetic power implies the close association ture has been found above/below the l =−1◦.2 region. of an energy source to this cloud.

◦ Multiple Supernova Explosions —The interaction with 4.3. Energy Source of the l =−1.2 Region a supernova explosion naturally explains the ellipti- The total kinetic energy of the five expanding shells cal morphology, broad velocity width, and expanding amounts to 1050.9 erg. By considering the typical bary- 51 shell kinematics of the molecular gas. The multiple onic energy supplied by an SN explosion, Ebar ∼ 10 shells may indicate multiple supernova explosions in the erg (Chevalier 1974), the energy conversion efficiency ◦ l =−1.2 region within their expansion times. Although (η = Ekin/Ebar) to the gas kinetic energy is as low there remains some subjectivity in the shell identifica- as ∼ 0.1–0.3 (Chevalier 1974; Salpeter 1976), and the tion, this interpretation can explain the observational total kinetic energy corresponds to ∼ 100.4–0.9 SNe. facts most naturally. The detection of the SiO J=8–7 Combined with the typical expansion time of the shells line, which is a good tracer of shocked molecular gas, at (∼ 105 yr), the SN rate is estimated to be ∼ 10−5–10−4 the high-velocity end of the S1 may support the origin yr−1. This high SN rate within the small volume of related to a local explosive event. Hereafter we discuss the l =−1◦.2 region indicates that a young massive star ◦ the l =−1.2 region basically along with the multiple cluster is embedded here. The absence of HII regions supernova explosion model. indicates that the cluster age is older than 10 Myr, the main sequence lifetime of a 15 M ⊙ star. The presence of 4.2. Physical Parameters a type-II SN indicates a cluster age of less than 30 Myr. Here, we estimate the physical parameters of the If we assume a uniform stellar age (instantaneous shells. In principle, the gas mass (M) can be derived star formation), Scalo initial mass function (IMF; Scalo by summing up the column density of gas in each shell. 1986) with lower/higher mass cutoffs of 0.08 M ⊙/100 M ⊙, Some shells overlap in the line-of-site direction, and it is and mass of heaviest stars of 8–15 M ⊙, the cluster mass 5.6 −1 difficult to separate individual shells in the data cube. is estimated to be Mcl ∼ 10 η M ⊙. This is two or- Thus, we simply define each shell as a finite thickness el- ders of magnitude larger than those of the Arches and 4 lipsoid in the l-b-V space. Pixels in an overlapping area Quintuplet clusters, ∼ 10 M ⊙ (Figer et al. 1999), and 6 Tsujimoto et al.

Table 2. Properties of Gases and Embedded Star Cluster 15.0

Region log10 Ekin texp log10 Pkin Mcl

− yr] S4 5 1 4 (erg) (10 year) (erg s ) (M⊙) l =+1◦. 3 52.6–53.0 4.5–5.3 39.5 106.6 η−1 S3 l =−1◦. 2 50.9 0.6–1.1 38.5 105.9 η−1 LS

S2 is an order of magnitude lower than that of the cluster S1 expansion time [10 embedded in the l =+1◦.3 region (Table 2). The absence of an infrared counterpart seems to be at 5.0 10.0 variance with the cluster-origin model. In fact, the to- –1.20 –1.25 –1.30 5.6 −1 tal luminosity of the young Mcl ∼ 10 η M ⊙ cluster Galactic Longitude [deg] 7.1 −1 amounts to Lcl ∼10 η L⊙ (Williams & Perry 1994), while the far-infrared luminosity of the l =−1◦.2 region is Figure 3. Calculated expansion times plotted against the 5.9 Galactic longitude of the center of each shell. The error ∼10 L⊙. The order of magnitude discrepancy may in- bars were calculated by substituting R with the major/minor dicate that the IMF in the hypothetical cluster is abnor- radius of each shell. mal. For instance, a shallower slope, as well as a higher low-mass cutoff of the IMF, significantly decrease the estimated cluster mass. Alternatively, an abundance of The discovery of expanding shells at the longitude of dark stellar remnants in the cluster could suggest other ∼−1◦.2 adds another molecular bubble in our Galaxy’s kinds of explosions, such as Type Ia SNe and neutron CMZ. The age of the embedded cluster is similar to that star mergers (Rosswog et al. 2013). of the cluster in the l =+1◦.3 region, which may have ◦ 4.4. Kinematics of Gas and Cluster been formed by the same mechanism as the l =−1.2 region. This suggests that active star formation may Figure 3 shows a plot of texp of the shells against the have occurred at multiple places in the innermost x1 Galactic longitude (l). The apparent gradient in texp orbit ∼107 yr ago. along l indicates a relative velocity of ∼ 200 km s−1 between the cluster and interacting gas. This relative velocity is similar to that of the Galactic rotation. The We are grateful to the Nobeyama Radio Observatory shorter t in the larger longitude implies that the clus- exp (NRO) and James Clerk Maxwell Telescope (JCMT) ter is relatively moving toward a positive Galactic lon- staff for their excellent support of the observations. The gitude relative to the interacting gas. NRO is a branch of the National Astronomical Observa- This situation can be understood by considering the tory of Japan, National Institutes of Natural Sciences. cluster and interacting gas on different orbits. First, The JCMT is operated by the East Asian Observatory we assume Sofue’s configuration of dense molecular gas on behalf of The National Astronomical Observatory (Sofue 1995) and Binney’s concept of gas kinematics in of Japan, Academia Sinica Institute of Astronomy and the CMZ (Binney et al. 1991). Large-scale CO J=1– Astrophysics, the Korea Astronomy and Space Science 0 data (Oka et al. 1998) show that the l =−1◦.2 and Institute, the National Astronomical Observatories of l =+1◦.3 regions correspond to the extensions of Sofue’s China, and the Chinese Academy of Sciences (Grant arms I and II, respectively. Figure 4 shows a schematic No. XDB09000000), with additional funding support face-on view of the CMZ. Accorging to Binney’s con- from the Science and Technology Facilities Council of cept, “sprayed gas” arises from locations where the in- the United Kingdom and participating universities in nermost x orbit and the outermost x orbit intersect. 1 2 the United Kingdom and Canada. This research was The sprayed gas flows into the inner orbits, and then supported by a grant from the Hayakawa Satio Fund collides with the other side of the innermost x orbit to 1 awarded by the Astronomical Society of Japan. T.O. generate a pair of shock fronts. If the hypothetical clus- acknowledges support from JSPS Grant-in-Aid for Sci- ter exists on the innermost x orbit, the presence of a 1 entific Research (B) No. 15H03643. S.T. acknowledges high relative velocity between the cluster and interacting support from JSPS Grant-in-Aid for Research Fellow gas can be understood. A parcel of sprayed gas might No. 15J04405. encounter the massive cluster, then receive expanding motions by SNe. Software: NOSTAR, Starlink(Currie et al. 2014) Another Bubble in the Central Molecular Zone of Our Galaxy 7

orbit 1 L innermost x =–1.2°

Arm II Sgr C

Sgr A*

Sgr B orbit x 2 Arm I

orbit 1

L=+1.3° innermost x To Sun

Figure 4. Schematic face-on view of the CMZ. The black dashed line represents the innermost x1 orbit, and the black cross shows the location of Sgr A*. The bottom left cyan object is the l =+1◦.3 complex, and the two cyan objects on arms I and II are the Sgr B and Sgr C regions, respectively. The orange arrow shows the direction along which the cluster is moving. The arrows show the direction along which the “sprayed” molecular gas is moving. The sprayed gas, shown as dotted and filled green objects, reaches to the other side of the innermost x1 orbit and encounters a cluster.

REFERENCES

Binney, J., Gerhard, O., E., Stark, A., A., Bally, J., and Oka, T., Hasegawa, T., Sato, F., Tsuboi, M., and Miyazaki, Uchida, K., I., 1991, MNRAS, 252, 210 A., 1998, ApJS, 118, 455 Buckle, J., V., Hills, R., E., Smith, H., et al., 2009, MNRAS Oka, T., White, G., J., Hasegawa, T., et al., 1999, ApJ, 399, 1026 515, 249 Chevalier, R., A., 1974, ApJ, 188, 501 Oka, T., Hasegawa, T., Sato, F., Tsuboi, M., and Miyazaki, Currie, M., J., Berry, D., S., Jennes, T., et al. 2014, ASPC, A., 2001, PASJ, 53, 787 485, 391 Oka, T., Nagai, M., Kamegai, K., Tanaka, K., and Kuboi, Figer, D., F., Kim, S., S., Morris, M., et al., 1999, ApJ, 525, N., 2007 PASJ, 59, 15 750 Oka, T., Onodera, Y., Nagai, M., et al., 2012. ApJS, 201, 14 Fukui, Y., Yamamoto, H., Fujishita, M., et al. 2006, Sci, Paglione, T., A., D., Jackson, J., M., Bolatto, A., D., and 314, 106 Heyer, M., H., 1998, ApJ, 493, 680 Habe, A., and Ohta, K., 1992, PASJ, 44, 203 Parsons, H., Dempsey, J., T., Thomas, H., S., et al. 2017, Haworth, T., J., Tasker, E., J., Fukui, Y., et al., 2015, ApJS, 234, 22 MNRAS, 450, 10 Rosswog, S., Piran, T., and Nakar, E., 2013, MNRAS, 430, Kaifu, N., Kato, T., and Iguchi, T., 1972, NPhS, 238, 105 2585 Langer, W., D., and Penzias, A., A., 1990, ApJ, 357, 477 Salpeter, E., E., 1976, ApJ, 206, 673 Liszt, H., S., 2006, A&A, 447, 533 Scalo J.M., 1986, FCPh, 11, 1 Matsumura, S., Oka, T., Tanaka, K., et al., 2012, ApJ, 756, Scoville, N., Z., 1972, ApJL, 175, L127 87 Sofue, Y., 1995, PASJ, 47, 527 Morris, M., Polish, N., Zuckerman, B., and Kaifu, N., 1983, Strickland, D. K., and Heckman, T. M. 2009, ApJ, 697, AJ, 88, 1228 2030 8 Tsujimoto et al.

Tanaka, K., Kamegai, K., Nagai, M., and Oka, T., 2007, Williams, R., J., R., Perry, J., J., 1994, MNRAS, 269, 538 PASJ, 59, 323 Tokuyama, S., Oka, T., Takekawa, S., et al., in prep. Torii, K., Hattori, Y., Hasegawa, K., et al., 2017, ApJ, 835, 142