This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

\draft\Received

\langlereception date\rangle \Accepted\langleacception date\rangle \Published\langlepublication date\rangle 11affiliationtext: Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency,
3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 252-5210, Japan
22affiliationtext: Department of Astronomy, The University of Tokyo, Bunkyo, Tokyo 113-0033, Japan 33affiliationtext: College of Arts and Sciences, J.F. Oberlin University, Machida, Tokyo 194-0294, Japan 44affiliationtext: National Radio Astronomy Observatory, P.O.Box O, Socorro, NM 87801-0387, USA 55affiliationtext: Japan Space Forum, Kanda-surugadai, Chiyoda-ku,Tokyo,101-0062, Japan 66affiliationtext: National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan \KeyWordsGalaxy: center1 — stars: formation2 — ISM: molecules3

Cloud-Cloud Collision in the Galactic Center Arc

Masato Tsuboi1    Yoshimi Kitamura1    Kenta Uehara2    Ryosuke Miyawaki3    Takahiro Tsutsumi4    Atsushi Miyazaki5    and Makoto Miyoshi6 tsuboi@vsop.isas.jaxa.jp
Abstract

We performed a search of cloud-cloud collision (CCC) sites in the Sagittarius A molecular cloud (SgrAMC) based on the survey observations using the Nobeyama 45-m telescope in the C32S J=10J=1-0 and SiO v=0J=21v=0~{}J=2-1 emission lines. We found candidates being abundant in shocked molecular gas in the Galactic Center Arc (GCA). One of them, M0.014-0.054, is located in the mapping area of our previous ALMA mosaic observation. We explored the structure and kinematics of M0.014-0.054 in the C32S J=21J=2-1, C34S J=21J=2-1, SiO v=0J=21v=0~{}J=2-1, H13COJ+=10{}^{+}J=1-0, and SO N,J=2,21,1N,J=2,2-1,1 emission lines and fainter emission lines. M0.014-0.054 is likely formed by the CCC between the vertical molecular filaments (VP) of the GCA, and other molecular filaments along Galactic longitude. The bridging features between these colliding filaments on the PV diagram are found, which are the characteristics expected in CCC sites. We also found continuum compact objects in M0.014-0.054, which have no counterpart in the H42α\alpha recombination line. They are detected in the SO emission line, and would be “Hot Molecular Core (HMC)”s. Because the LTE mass of one HMC is larger than the virial mass, it is bound gravitationally. This is also detected in the CCS emission line. The embedded star would be too young to ionize the surrounding molecular cloud. The VP is traced by poloidal magnetic field. Because the strength of the magnetic field is estimated to be m\sim mGauss using the CF method, the VP is supported against fragmentation. The star formation in the HMC of M0.014-0.054 is likely induced by the CCC between the stable filaments, which may be a common mechanism in the SgrAMC.

1 Introduction

The Galactic center region is the nucleus of the nearest spiral galaxy. The Central Molecular Zone (CMZ) (Morris & Serabyn, 1996) is a large molecular cloud reservoir in the galaxy, which extends along the galactic plane up to l±1l\sim\pm 1^{\circ}. The mass of the CMZ is estimated to be 5×107\sim 5\times 10^{7} M (e.g. Morris & Serabyn (1996); Tsuboi, Handa & Ukita (1999)). The molecular clouds in the CMZ are much denser, warmer, and more turbulent than those in the Galactic disk region. The CMZ is recognized to be a laboratory for peculiar phenomena, which will be found in central molecular cloud reservoirs of normal external galaxies by future telescopes. Young and luminous star clusters which contain over several ten OB stars, for example Arches cluster and Quintuplet cluster, have been found in the CMZ by IR observations (e.g. Genzel et al. (1996); Figer et al. (1999); Figer et al. (2002)). They are as luminous as those which are nearly hard to be found in the Galactic disk region. These star clusters presumably have been formed in the cradle molecular clouds in the CMZ. The star formation would be influenced by external factors, such as interactions with SNRs and/or cloud-cloud collisions(CCC) because they are crowded in the region (e.g. Morris (1993); Hasegawa et al. (1994); Hasegawa et al. (2008)). However, it is difficult to demonstrate observationally how the cradle molecular clouds produce such massive clusters because almost these clusters have already lost the surrounding molecular materials. The Galactic Center 50 kms-1 molecular cloud (50MC) is an exception, which has still abundant molecular gas and several compact HII regions. In the previous observations, we have found a half-shell structure filled with shocked molecular gas, which would be made by CCC (Tsuboi et al. (2011), Tsuboi et al. (2015), Tsuboi et al. (2019)). As known widely, the dense molecular clouds in the CMZ seem to exist as molecular ridges along the galactic plane, which are bundles of the molecular filaments (Bally et al. (1987); Oka et al. (1998); Tsuboi, Handa & Ukita (1999)). The molecular filaments would collide each other, where star formation would be activated. It is an open issue whether CCCs usually induce star formation in the CMZ or not.

First, we searched CCC sites based on the survey observations for the Sagittarius A molecular cloud complex (SgrAMC) with the Nobeyama 45-m telescope (Tsuboi, Handa & Ukita, 1999; Tsuboi et al., 2011). The SgrAMC is one of most conspicuous molecular cloud complexes in the CMZ. One of the candidates, M0.014-0.054 near the 50MC, is located serendipitously in the mosaic observation area using ALMA of the 50MC although the observation itself had other science objective (Uehara et al., 2019). Then we looked for the signs of CCCs and induced star formation in the ALMA data. Throughout this paper, we adopt 8 kpc as the distance to the Galactic center (e.g. Boehle et al. (2016)). Then, 1\arcsec1\arcsec corresponds to about 0.04 pc at the distance.

Refer to caption
Figure 1: a Integrated intensity map of the Galactic Center Arc (GCA) in the C32S J=10J=1-0 emission line with the Nobeyama 45-m telescope (Tsuboi, Handa & Ukita (1999)). The velocity range is VLSR=30V_{\mathrm{LSR}}=-30 to 0 km s-1. The angular resolution is 45\arcsec45\arcsec in FWHM, which indicates as a white circle in the lower left corner. A red arrow indicates a candidate of the CCC site; M0.014-0.054. The contours show the 1.44 GHz (20-cm) continuum emission for comparison (Yusef-Zadeh, Morris, & Chance (1984)). The contour levels are (1,2,3,4,5,6,7,8,9,10,15,20,30,and50)×95(1,2,3,4,5,6,7,8,9,10,15,20,30,\mathrm{and}~{}50)\times 95 K in TBT_{\mathrm{B}}. The angular resolution is 19\arcsec×15\arcsec,PA=10519\arcsec\times 15\arcsec,PA=105^{\circ} in FWHM. b Position velocity diagram along “Vertical Part (VP)” of the GCA extending roughly north and south. The sampling area is shown as a white rectangle in a. A red arrow indicates M0.014-0.054. c Integrated intensity map of the GCA in the SiO v=0,J=21v=0,J=2-1 emission line with the Nobeyama 45-m telescope (Tsuboi et al. (2011)). The velocity range is the same as that of a. The angular resolution is 26\arcsec26\arcsec in FWHM, which indicates as a white circle in the lower left corner. A red arrow indicates M0.014-0.054. The red broken line square indicates the mapping area of the ALMA Cycle 1 observation (2012.1.00080.S. PI M.Tsuboi). d The 330 MHz (90-cm) continuum emission for comparison (LaRosa et al., 2000). The angular resolution is 43\arcsec×24\arcsec,PA=6543\arcsec\times 24\arcsec,PA=65^{\circ} in FWHM. The yellow arrows show “threads”.

2 Search for the Cloud-Cloud Collision Candidates

We searched the CCC candidates in the SgrAMC based on the existing survey observations with the Nobeyama 45-m telescope in the C32S J=10J=1-0 (48.99095748.990957 GHz) and SiO v=0,J=21v=0,J=2-1 (86.84699586.846995 GHz) emission lines (C32S; Tsuboi, Handa & Ukita (1999), SiO; Tsuboi et al. (2011)). These surveys have the highest angular resolutions among single dish observations. The C32S emission line is a tracer of medium dense molecular cloud, n(H2)cl104n({\mathrm{H}}_{2})_{\mathrm{cl}}\sim 10^{4} cm-3. Because this emission line is moderately optically thick, τ0.53\tau\sim 0.5-3, in the SgrAMC (e.g. Tsuboi, Handa & Ukita (1999)), this emission shows mainly the location of the medium dense molecular cloud. On the other hand, the SiO emission line is a famous tracer of strong C-shock wave (ΔV>30\Delta V>30 km s-1) which propagated in the region within 10510^{5} yr (e.g. Gusdorf et al. (2008), Jiménez-Serra et al. (2008)). Because the emission line is often detected in the CCC sites, the detection indicates the CCC candidates.

Figure 1a is the integrated intensity map of the SgrAMC in the C32S J=10J=1-0 emission line with the velocity range of VLSR=30V_{\mathrm{LSR}}=-30 to 0 km s-1. The contours in the figure show the Galactic Center Arc (GCA) in the 20-cm continuum emission for comparison (Yusef-Zadeh, Morris, & Chance (1984)). The molecular cloud is identified as two curved ridges along the GCA (e.g. Serabyn & Güsten (1987)). Moreover, a large molecular cloud connecting with the curved ridges is located apparently in the “continuum gap” between the GCA and Sagittarius A east SNR. This cloud extends vertically north, l0.02,b0.05l\sim-0.02^{\circ},~{}b\sim 0.05^{\circ}, and south, l0.03,b0.12l\sim 0.03^{\circ},~{}b\sim-0.12^{\circ}. We call it the “Vertical Part (VP)” here. M0.014-0.054 is identified as a strong compact feature which is apparently located on the VP in the C32S integrated intensity map. Figure 1b is the position-velocity (PV) diagram in the C32S J=10J=1-0 emission line along the VP. We identified the VP as a ridge-like feature with no noticeable velocity gradient on the position-velocity diagram. M0.014-0.054 is also seen as a compact feature on the ridge.

While Figure 1c is the integrated intensity map in the SiO v=0,J=21v=0,J=2-1 emission line with the same area and velocity range as in Figure 1a. Almost all molecular clouds in the area become faint or disappear in the SiO map. This shows that no strong C-shock wave propagated in the area except for several compact components including M0.014-0.054. M0.014-0.054 extends roughly east and west. The east-west extent of M0.014-0.054 seems to correspond with the width of the VP seen in the C32S emission line. M0.014-0.054 is more prominent in the SiO emission line than in the C32S emission line. This shows that shocked molecular gas originated by the C-shock wave within 10510^{5} yr is abundant in M0.014-0.054. This suggests that the object would be made by some historical event including cloud-cloud collision.

In addition, Figure 1d shows the continuum emission map at 330 MHz (LaRosa et al., 2000). Well-known continuum features, “threads” are identified (yellow arrows). There is also a emission extending from the Sagittarius A east SNR to east, which crosses “threads” and reaches up to l0.10,b0.06l\sim 0.10^{\circ},~{}b\sim-0.06^{\circ}. Because M0.014-0.054 is adjacent to these continuum features, this would be associated physically with them.

3 ALMA Observation

We have performed the observation of a 330\arcsec×330\arcsec330\arcsec\times 330\arcsec area covering the 50MC in the C32S J=21J=2-1 (97.98095397.980953 GHz), C34S J=21J=2-1 (96.41295096.412950 GHz), H13COJ+=10{}^{+}J=1-0 (86.75428886.754288 GHz), SiO v=0J=21v=0~{}J=2-1 (86.84699586.846995 GHz), CH3OH (96.739363,96.741377,96.739363,96.741377, and 96.74454996.744549 GHz), C2H (87.31692587.316925 and 87.32862487.328624 GHz), c-C3H2 JKa,Kc=21,210,1J_{K_{a},K_{c}}=2_{1,2}-1_{0,1} (85.33890685.338906 GHz), and SO N,J=2,21,1N,J=2,2-1,1 (86.09398386.093983 GHz) emission lines and many fainter emission lines as ALMA Cycle1 observation (2012.1.00080.S. PI M.Tsuboi). We also have observed the continuum emission at 86 GHz, simultaneously. The entire ALMA observation consists of a 137 pointing mosaic of the 12-m array and a 52 pointing mosaic of the 7-m array (ACA). Additionally, the single-dish data has been obtained by Total Power Array. The red broken line square in Figure 1c indicates the mapping area. It is fortunate that M0.014-0.054 is located in the mapping area of the ALMA observation by chance. The molecular emission lines used for imaging here are summarized in Table 1.

Table 1: Molecular emission lines used for imaging.
Name transition Rest frequency Remarks
νrest\nu_{\mathrm{rest}}[GHz]
C32S J=21J=2-1 97.98095397.980953 medium dense molecular cloud
moderately optically thick
34SO N,J=2,31,2N,J=2,3-1,2 97.71540197.715401 “Hot Molecular Core”
CH3OH JKa,Kc=21,111,0AJ_{K_{a},K_{c}}=2_{1,1}-1_{1,0}A_{--} 97.58280897.582808 mild C-shock
CH3OH JKa,Kc=21,211,1EJ_{K_{a},K_{c}}=2_{-1,2}-1_{-1,1}E 96.73936396.739363 mild C-shock
CH3OH JKa,Kc=20,210,1A++J_{K_{a},K_{c}}=2_{0,2}-1_{0,1}A_{++} 96.74137796.741377 mild C-shock
CH3OH JKa,Kc=20,210,1EJ_{K_{a},K_{c}}=2_{0,2}-1_{0,1}E 96.74454996.744549 mild C-shock
CH3CHO JKa,Kc=52,342,2EJ_{K_{a},K_{c}}=5_{2,3}-4_{2,2}E 96.47552396.475523
C34S J=21J=2-1 96.41295096.412950 dense molecular cloud (opt. thin)
C2H N=10J=3/21/2F=21N=1-0~{}J=3/2-1/2~{}F=2-1 87.31692587.316925 dense regions exposed to UV radiation
C2H N=10J=3/21/2F=10N=1-0~{}J=3/2-1/2~{}F=1-0 87.32862487.328624 dense regions exposed to UV radiation
HN13C J=10F=01J=1-0~{}F=0-1 87.09073587.090735 dense molecular cloud (opt. thin)
HN13C J=10F=21J=1-0~{}F=2-1 87.09085987.090859 dense molecular cloud (opt. thin)
HN13C J=10F=11J=1-0~{}F=1-1 87.09094287.090942 dense molecular cloud (opt. thin)
28SiO v=0J=21v=0~{}J=2-1 86.84699586.846995 strong C-shock
H13CO+ J=10J=1-0 86.754288 dense molecular cloud (opt. thin)
CCS N,J=7,66,5N,J=7,6-6,5 86.18141386.181413 early stage of star formation
32SO N,J=2,21,1N,J=2,2-1,1 86.09398386.093983 “Hot Molecular Core”
HCOOH JKa,Kc=41,431,3J_{K_{a},K_{c}}=4_{1,4}-3_{1,3} 86.5461886.54618
HC15N J=10J=1-0 86.05496786.054967 dense molecular cloud (opt. thin)
29SiO v=0J=21v=0~{}J=2-1 85.75918885.759188 strong C-shock
HOCO+ JKa,Kc=40,430,3J_{K_{a},K_{c}}=4_{0,4}-3_{0,3} 85.53148085.531480 dense molecular cloud (opt. thin)
c-C3H2 JKa,Kc=21,210,1J_{K_{a},K_{c}}=2_{1,2}-1_{0,1} ortho 85.33890685.338906 dense regions exposed to UV radiation
Hydrogen atom H42α\alpha 85.688485.6884 ionized gas

https://physics.nist.gov/cgi-bin/micro/table5/start.pl. e.g. Nagy et al. (2015)

The C34S and H13CO+ emission lines are tracers of dense molecular cloud, n(H2)cl105n({\mathrm{H}}_{2})_{\mathrm{cl}}\sim 10^{5} cm-3. These emission lines are thought to be optically thin even in the SgrAMC. The optical thickness of the C34S emission line will be demonstrated to be thin for the typical case in the subsection 4.4. That of the H13CO+ emission line is estimated to be τ<0.2\tau<0.2 based on the observed mean temperature ratio, TBT_{\mathrm{B}}(H12CO+)/TBT_{\mathrm{B}}( H13CO+) 7\sim 7 (e.g. Armijos-Abendaño et al. (2015)). As mentioned previously, the SiO emission line is a well-known tracer of strong C-shock wave, while the CH3OH molecules are enhanced even by mild C-shock (ΔV10\Delta V\sim 10 km s-1 e.g. Hartquist et al. (1995)). The C2H and c-C3H2 emission lines are well-known tracers of photodissociation regions (PDRs). The abundance of both radicals was found to increase close to the dissociation front (e.g. Jansen et al. (1995); Fuente, Rodriguez-Franco, & Martin-Pintado (1996)). The HN13C and HC15N emission lines have critical densities in the environment of the SgrAMC as high as n(H2)cl106n({\mathrm{H}}_{2})_{\mathrm{cl}}\sim 10^{6} cm-3 and n(H2)cl107n({\mathrm{H}}_{2})_{\mathrm{cl}}\sim 10^{7} cm-3, respectively. These intensities can vary widely by their path length, abundance, and isotope ratio. The SO and 34SO emission lines are usually emitted from “Hot Molecular Core (HMC)”s where the molecular cloud is heated up to 100 K by newborn stars. The H42α\alpha recombination line and continuum emission at 86 GHz are simultaneously observed. The H42α\alpha recombination line traces ionized gas. The continuum emission at 86 GHz is usually emitted by ionized gas. However, there is a possibility that the emission is originated by warm dust or non-thermal mechanism in the Galactic center region. In addition, the CCS emission line is partly in the observation frequency range. The CCS molecule is abundant in the early stage of star formation but decreases with increasing time (e.g. Hirahara et al. (1992)).

The resultant maps have FWHM angular resolutions of 2\farcs5×1\farcs8,PA30\sim 2\farcs 5\times\sim 1\farcs 8,PA\sim-30^{\circ} and 1\farcs9×1\farcs3,PA37\sim 1\farcs 9\times\sim 1\farcs 3,PA\sim-37^{\circ} using “natural weighting” and “briggs weighting” as u-v sampling, respectively. The original velocity resolution is 1.71.7 km s-1(488 kHz). The typical rms noise levels in the emission free areas of the resultant maps using “natural weighting” and “briggs weighting” are 0.003\sim 0.003 Jy/beam/1.71.7 km s-1 or 0.11\sim 0.11 K in TBT_{\mathrm{B}} and 0.005\sim 0.005 Jy/beam/1.71.7 km s-1 or 0.17\sim 0.17 K in TBT_{\mathrm{B}}, respectively. J0006-0623, J1517-2422, J717-3342, J1733-1304, J1743-3058, J1744-3116 and J2148+0657 were used as phase calibrators. The flux density scale was determined using Titan, Neptune and Mars. The calibration and imaging of the data were done by CASA (McMullin et al., 2007).

Although the area around M0.014-0.054 is less dense comparing to the 50MC (Uehara et al., 2019), the C32S and C34S emission lines are significantly resolved out only by the interferometer observations when they have been processed in the same procedure mentioned above. This is because the molecular cloud is widely extended in the area. The combining with the single-dish data, Total Power Array data, is required to recover the missing flux. We performed the procedure using the CASA task ”FEATHER”. We present the detailed description of the procedure and the full results in other papers (Uehara et al. (2019), Uehara et al. (2020)).

4 Results

Refer to caption
Figure 2: Integrated intensity maps of M0.014-0.054 and the “Vertical Part” in the C32S J=21J=2-1(a), C34S J=21J=2-1(b), H13COJ+=10{}^{+}J=1-0(c), SiO v=0J=21v=0~{}J=2-1 (d), CH3OH (96.741 GHz)(e), C2H (87.317 GHz)(f), c-C3H2 (g), HN13C (h), HC15N (i), SO N,J=2,21,1N,J=2,2-1,1 (j), and H42α\alpha(k) emission lines with ALMA. The C32S map is before the combining with the single-dish data. The velocity range is VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1. The angular resolutions are 2\farcs5×1\farcs8\sim 2\farcs 5\times\sim 1\farcs 8 in FWHM and PA30PA\sim-30^{\circ}, which indicate as white ovals in the lower left corners. The rms noise of these maps is 0.036\sim 0.036 K or 0.90.9 K km s-1. The panel (l) shows the map of 850 μ\mum continuum with JCMT is also shown for comparison (pseudo color: Pierce-Price et al. (2000)). The overlaid contours show the continuum emission at 86 GHz with ALMA. The contour levels are (1,2,4,8,16)×0.023(1,2,4,8,16)\times 0.023 K. The angular resolution of JCMT is 17\arcsec\sim 17\arcsec in FWHM.
Refer to caption
Figure 2: Continued.

4.1 Distributions of Molecular Line Emissions around M0.014-0.054 and the “Vertical Part”

Figures 2a, 2b, 2c, 2d, and 2e show the integrated intensity maps with ALMA of M0.014-0.054 and the VP in the C32S, C34S, H13CO+, SiO, and CH3OH (blended lines around 96.741 GHz) emission lines, respectively. Note that the data for the C32S and C34S maps are not combined with the single-dish data. The velocity range of each map is VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1. The velocity range was determined from Figures 1b and 1c.

M0.014-0.054 and the VP are clearly detected in the map of the C32S emission line. The VP is resolved into the bundles of the molecular filaments extending roughly from north to south. M0.014-0.054 is extended roughly from west to east and over the filaments. The filaments line up nearly at a right angle to M0.014-0.054 in the far distance area. However, they become disordered and tangled in the vicinity of M0.014-0.054. The positions of the peaks in the C32S emission line do not correspond to those in the maps of the other emission lines which will be mentioned in the following part. This is probably because the C32S emission line is partially optically thick in the region. M0.014-0.054 is also clearly detected in the maps of the C34S and H13CO+ emission lines. They resemble each other very well. These show that dense molecular cloud with n(H2)cl105n({\mathrm{H}}_{2})_{\mathrm{cl}}\sim 10^{5} cm-3 is abundant in M0.014-0.054. In contrast, the molecular filaments in the C34S and H13CO+ emission lines are fainter than those in the C32S emission line, suggesting that this shows that the molecular cloud in the filaments is less dense than that in M0.014-0.054.

M0.014-0.054 is clearly detected in the maps of the SiO and CH3OH emission lines. The detection in the SiO emission line shows that a strong C-shock wave (ΔV>30\Delta V>30 km s-1) propagated within 10510^{5} yr over the entire M0.014-0.054. Although the molecular filaments of the VP are also detected in the maps of these emission lines, they are faint on the whole. For example, the prominent parts are located around l0.010,b0.060l\sim 0^{\circ}.010,~{}b\sim-0^{\circ}.060, l0.015,b0.087l\sim 0^{\circ}.015,~{}b\sim-0^{\circ}.087, and l0.025,b0.035l\sim 0^{\circ}.025,~{}b\sim-0^{\circ}.035 in the maps of the C32S and H13CO+ emission lines. However, they almost disappear in the SiO emission line.

The C2H (87.317 GHz) and c-C3H2 emission lines show a similar distribution in M0.014-0.054 (see Figures 2f and 2g). The other C2H (87.329 GHz) emission line is also detected although this line is fainter. The VP almost disappears in the c-C3H2 emission line. It has been reported that the NN[C2H]/NN[c-C3H2] column density ratio increases by the UV radiation (32\sim 32 at the PDR edge; Cuadrado et al. (2015), up to 80\sim 80 at planetary nebulae; Schmidt, Zack, & Ziurys (2018)). Because the VP is less dense than M0.014-0.054, even the interior of the VP should have the high ratio. The c-C3H2 and HN13C emission lines show a similar distribution in M0.014-0.054 but the latter is faint (see Figures 2g and 2h). The HC15N and SO emission lines show a similar distribution (see Figures 2i and 2j). M0.014-0.054 is identified as a faint extended feature with a single compact peak in the maps of the HC15N and SO emission lines. The peak of the feature is located in the vicinity of strongest peaks in the other emission lines except for the C32S emission line. The filaments of the VP disappears in the HC15N and SO emission lines.

The structures mentioned above including M0.014-0.054 and the filaments of the VP have no ionized gas counterparts which appears in the H42α\alpha recombination line (see Figure 2k). However, two distinct components in M0.014-0.054, are identified in the continuum emission at 86 GHz (see Figure 2l). The east one corresponds to the peak in the SO emission line. As mentioned previously, the SO emission line is a tracer of HMCs. Then the peaks would involve HMCs. The map of the 850 μ\mum continuum with JCMT (Pierce-Price et al., 2000) is shown for comparison. The whole of M0.014-0.054 is also detected in the 850 μ\mum continuum although the VP is not identified. Two peaks detected at 850 μ\mum are associated with the peaks detected at 86 GHz.

Refer to caption
Figure 3: a Integrated intensity map with VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1 in the C32S J=21J=2-1 emission line. This is after the combining with the single-dish data. The angular resolution is 1\farcs9×1\farcs3\sim 1\farcs 9\times\sim 1\farcs 3 in FWHM and PA37PA\sim-37^{\circ}, which indicates as an oval in the lower right corner. The rectangle shows the sampling area of b. b Position-velocity diagram in the C32S J=21J=2-1 emission line along the “Vertical Part” (white rectangle in a and c). c Integrated intensity map with VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1 in the C32S J=21J=2-1 emission line.
Refer to caption
Figure 4: a Integrated intensity map with VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1 in the C34S J=21J=2-1 emission line. This is after the combining with the single-dish data. The angular resolution is 1\farcs9×1\farcs3\sim 1\farcs 9\times\sim 1\farcs 3 in FWHM and PA37PA\sim-37^{\circ}, which indicates as an oval in the lower right corner. The rectangle shows the sampling area of b. b Position-velocity diagram in the C34S J=21J=2-1 emission line along the “Vertical Part”. c Integrated intensity map with VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1 in the C34S J=21J=2-1 emission line.

4.2 Kinematics along the “Vertical Part”

Figures 3a and 4a show the integrated intensity maps with VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1 in the C32S J=21J=2-1 and C34S J=21J=2-1 emission lines, respectively. Because the data of these maps are combined with the single-dish data, the mapping areas become narrow slightly. The appearances of the VP in these emission lines resemble each other very well. Figures 3b and 4b show the PV diagrams along the VP in the C32S J=21J=2-1 and C34S J=21J=2-1 emission lines, respectively. The sampling areas are shown as the rectangles in Figures 3a and 4a. In the PV diagrams, the VP is identified as a long feature with VLSR30V_{\mathrm{LSR}}\sim-30 to 0 km s-1. There is no clear velocity gradient in the feature. The appearance and velocity gradient are consistent with those of the corresponding feature shown in Figure 1b. On the other hand, M0.014-0.054 is identified as a compact curved feature on the filaments of the VP (Figures 3b and 4b). This feature suggests that M0.014-0.054 had been affected by any external interaction (e.g. Haworth et al. (2015), Haworth et al. (2015b)). Additionally, a large extended feature with VLSR30V_{\mathrm{LSR}}\sim 30 to 8080 km s-1 is also identified. The appearances and kinematics of these features are consistent with those of the corresponding features shown in Figure 1b.

In the PV diagrams, there are two compact components with VLSR20V_{\mathrm{LSR}}\sim 20 km s-1 adjoining the extended feature mentioned above. In order to clarify the angular extension of these components, the integrated intensity maps with VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1 are shown in Figures 3c and 4c. These features are resolved into two molecular filaments which are almost parallel to each other in the maps. The two components in the PV diagrams are thought to be a part of the molecular filament bundle, which will be discussed in the following subsection.

In addition, a faint feature is also identified around the angular offset of M0.014-0.054 and velocity of VLSR50V_{\mathrm{LSR}}\sim-50 km s-1 in the C34S J=21J=2-1 emission line (see Figure 4b). The feature is the contamination by the CH3CHO JKa,Kc=52,442,3EJ_{K_{a},K_{c}}=5_{-2,4}-4_{-2,3}E emission line (96.42562096.425620 GHz) because this feature is not detected in Figure 3b.

Refer to caption
Figure 5: a Integrated intensity maps around the ”Molecular Filament Bundle” with VLSR=40V_{\mathrm{LSR}}=-40 to 1010 km s-1 (pseudo color) and VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1 (contours) in the C32S J=21J=2-1 emission line. These are after the combining with the single-dish data. The angular resolution is 1\farcs9×1\farcs3\sim 1\farcs 9\times\sim 1\farcs 3 in FWHM and PA37PA\sim-37^{\circ}, which indicates as an oval in the lower left corner. The contour levels are 32,48, 64, and 80 K km s-1. b Integrated intensity maps around the ”Molecular Filaments Bundle” with VLSR=40V_{\mathrm{LSR}}=-40 to 1010 km s-1 (pseudo color) and VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1 (contours) in the C34S J=21J=2-1 emission line. These are after the combining with the single-dish data. The angular resolution is 1\farcs9×1\farcs3\sim 1\farcs 9\times\sim 1\farcs 3 in FWHM and PA37PA\sim-37^{\circ}, which indicates as an oval in the lower left corner. The contour levels are 4, 6, 8, and 10 K km s-1.

4.3 “Molecular Filament Bundle”

Another molecular cloud ridge extends roughly east and west over this area (see Figures 5a and 5b). The western and eastern parts of this structure are seen in the negative and positiveLSR velocities, respectively. We call this structure “Molecular Filament Bundle” (MFB). Figure 5a and 5b show the integrated intensity maps around the MFB with VLSR=40V_{\mathrm{LSR}}=-40 to 1010 km s-1 (pseudo color) and VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1 (contours) in the C32S J=21J=2-1 and C34S J=21J=2-1 emission lines, respectively. The MFB would be extended over the mapping area. The MFB is almost along the Galactic longitude but curved slightly. The MFB is resolved into the bundle of a few thin molecular filaments in these maps. The width of the MFB is changed periodically from 40\arcsec\sim 40\arcsec around l0.035l\sim-0^{\circ}.035 and l0.005l\sim 0^{\circ}.005 to 20\arcsec\sim 20\arcsec around l0.050l\sim-0^{\circ}.050 and l0.025l\sim-0^{\circ}.025. These filaments are seem to be twisted each other. The appearance of the MFB is similar to each other in both the emission lines. M0.014-0.054 is located around the apparent intersection of the MFB and VP. There is shocked molecular gas abundantly in M0.014-0.054 as shown in Section 4.1. Although the velocity of M0.014-0.054 is similar to that of the VP, it is VLSR40V_{\mathrm{LSR}}\sim 40 km s-1 different from that of the MFB. As will be discussed in detail later, M0.014-0.054 is connected to the MFB by the ”Bridge” on the PV diagram and a collisionally excited maser spot is also located in M0.014-0.054. These features demonstrate that the MFB and VP collided each other within 10510^{5} years and M0.014-0.054 was made in the interaction by the CCC consequently. In addition, M0.014-0.054 is resolved into a series of compact features in these figures. The strongest feature (labeled A) and second one (labeled B) in the CS emission line maps are centered at l0.021,b0.051l\sim 0.021^{\circ},b\sim-0.051^{\circ} (the left circle) and l0.014,b0.052l\sim 0.014^{\circ},b\sim-0.052^{\circ} (the right circle), respectively.

4.4 Molecular Could Masses of M0.014-0.054, “Vertical Part”, and “Molecular Filament Bundle”

The derivation of molecular cloud mass should suffer from the optical thickness of the observed emission line. The ratio of the observed brightness temperatures of the C32S and C34S emission lines is given by

R=TB(C32S)TB(C34S)=f(C32S)T(C32S)ex(1eτ)f(C34S)T(C34S)ex(1eτ22.35),R=\frac{T_{\mathrm{B}}\mathrm{(C^{32}S)}}{T_{\mathrm{B}}\mathrm{(C^{34}S)}}=\frac{f\mathrm{(C^{32}S)}T\mathrm{{}_{ex}(C^{32}S)}(1-e^{-\tau})}{f\mathrm{(C^{34}S)}T\mathrm{{}_{ex}(C^{34}S)}(1-e^{\frac{-\tau}{22.35}})}, (1)

where TexT\mathrm{{}_{ex}} and ff are the excitation temperature and beam-filling factor of the emission line, respectively. The isotope abundance ratio of 32S and 34S in the molecular cloud is assumed to be equal to the natural abundance ratio of 22.3522.35111https://physics.nist.gov/cgibin/Compositions/stand_alone.plhttps://physics.nist.gov/cgi-bin/Compositions/stand\_alone.pl.

The mean line intensity ratio in M0.014-0.054 is calculated to be R9.4R\sim 9.4. Therefore the mean optical thickness of the C32S emission line toward M0.014-0.054 is estimated to be τ2.1\tau\sim 2.1 assuming that T(C32S)ex=T(C34S)exT\mathrm{{}_{ex}(C^{32}S)}=T\mathrm{{}_{ex}(C^{34}S)} and f(C32S)=f(C34S)f(\mathrm{C^{32}S})=f(\mathrm{C^{34}S}). The correction factor for the optical thickness of the C32S emission line is τ1eτ2.4\frac{\tau}{1-e^{-\tau}}\sim 2.4. The corrected integrated intensity of the object in the CS J=21J=2-1 emission line is given by

I(CS)corr.=τ1eτobjectT(C32S)Bdvds.I\mathrm{{}_{corr.}(CS)}=\frac{\tau}{1-e^{-\tau}}\int_{\mathrm{object}}\int T\mathrm{{}_{B}(C^{32}S)}dvds. (2)

The total number of the H2 molecules in the object is given by

objectN(H2)LTEds=4.56×1011e2.35TexI(CS)corr.X(CS)(1e4.7Tex),\int_{\mathrm{object}}N{\mathrm{{}_{LTE}(H_{2})}}ds=\frac{4.56\times 10^{11}e^{\frac{2.35}{T_{\mathrm{ex}}}}I\mathrm{{}_{corr.}(CS)}}{X\mathrm{(CS)}(1-e^{\frac{-4.7}{T_{\mathrm{ex}}}})}, (3)

where N(H2)LTEN{\mathrm{{}_{LTE}(H_{2})}} is the LTE molecular column density. Here the Einstein AA coefficient of the CS J=21J=2-1 emission line is assumed to be A21=2.2×105A_{21}=2.2\times 10^{-5} s-1.

The excitation temperature of the CS J=21J=2-1 emission line is assumed to be Tex=80T_{\mathrm{ex}}=80 K (TK=80T_{\mathrm{K}}=80 K in Ao et al. (2013)). The fractional abundance of the CS molecule is assumed to be X(CS)=N(CS)N(H2)=1×108X\mathrm{(CS)}=\frac{N\mathrm{(CS)}}{N\mathrm{(H_{2})}}=1\times 10^{-8}, which is usually used for molecular clouds in the disk region. The LTE molecular cloud mass is given by

MLTE[M]=μobjectN(H2)LTEds2.4×1038TexI(CS)corr.,M_{\mathrm{LTE}}[M_{\odot}]=\mu\int_{\mathrm{object}}N{\mathrm{{}_{LTE}(H_{2})}}ds\simeq 2.4\times 10^{-38}T_{\mathrm{ex}}I\mathrm{{}_{corr.}(CS)}, (4)

where μ\mu is the mean molecular weight per H2 molecule: μ=2.8\mu=2.8 in amu =2.4×1057M=2.4\times 10^{-57}M_{\odot}. The integrated intensity of the whole of M0.014-0.054 in the C32S emission line is object4010TB(C32S)𝑑v𝑑s=1.15×1040[Kkms1cm2]\int_{\mathrm{object}}\int_{-40}^{10}T_{\mathrm{B}}\mathrm{(C^{32}S)}dvds=1.15\times 10^{40}\mathrm{[K~{}kms^{-1}cm^{2}]}. The LTE molecular cloud mass is estimated to be MC32S,LTE5.3×104(Tex80)[M]M_{\mathrm{C32S,LTE}}\simeq 5.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})[M_{\odot}].

Although the mean optical thickness of the C32S emission line is fairly high in M0.014-0.054, that of the C34S emission line is estimated to be as small as τ2.122.350.09\tau\sim\frac{2.1}{22.35}\sim 0.09. We can estimate the molecular cloud mass using the C34S emission line data without the correction for the optical thickness. The integrated intensity of the whole of M0.014-0.054 in the C34S emission line is 4010TB(C34S)𝑑v𝑑s=1.22×1039[Kkms1cm2]\int\int_{-40}^{10}T_{\mathrm{B}}\mathrm{(C^{34}S)}dvds=1.22\times 10^{39}\mathrm{[K~{}kms^{-1}cm^{2}]}. The LTE molecular cloud mass is estimated to be MC34S,LTE5.3×104(Tex80)[M]M_{\mathrm{C34S,LTE}}\simeq 5.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})[M_{\odot}]. The LTE molecular cloud mass is consistent with that from the C32S emission line observation. These are summarized in Table 2.

Using the same procedure and the C32S emission line data, the LTE molecular cloud masses of the VP and MFB are MC32S,LTE4.3×104(Tex80)[M]M_{\mathrm{C32S,LTE}}\simeq 4.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})[M_{\odot}] and MC32S,LTE8.4×104(Tex80)[M]M_{\mathrm{C32S,LTE}}\simeq 8.4\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})[M_{\odot}], respectively. Meanwhile, using the C34S emission line data, the LTE molecular cloud masses of the VP and MFB are MC34S,LTE4.2×104(Tex80)[M]M_{\mathrm{C34S,LTE}}\simeq 4.2\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})[M_{\odot}] and MC34S,LTE8.1×104(Tex80)[M]M_{\mathrm{C34S,LTE}}\simeq 8.1\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})[M_{\odot}], respectively. These are also summarized in Table 2.

Table 2: Molecular Could Masses
Object τ¯\bar{\tau}(C32S) MC32S,LTEM_{\mathrm{C32S,LTE}} MC34S,LTEM_{\mathrm{C34S,LTE}} MvirM_{\mathrm{vir}} MvirMLTE\frac{M_{\mathrm{vir}}}{M_{\mathrm{LTE}}}
[MM_{\mathrm{\odot}}] [MM_{\mathrm{\odot}}] [MM_{\mathrm{\odot}}]
M0.014-0.054 2.12.1 5.3×104(Tex80)5.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) 5.3×104(Tex80)5.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) - -
Vertical Part 1.31.3 4.3×104(Tex80)4.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) 4.2×104(Tex80)4.2\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) - -
Molecular Filament Bundle 1.81.8 8.4×104(Tex80)8.4\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) 8.1×104(Tex80)8.1\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) - -
Object A 3.83.8 - 1.3×104(Tex80)1.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) 5.0×1035.0\times 10^{3} 0.4
Object B 2.72.7 - 1.2×104(Tex80)1.2\times 10^{4}(\frac{T_{\mathrm{ex}}}{80}) - -

Typical error of these masses is estimated to be 30\sim 30 %.

5 The Cloud-Cloud Collision between the “Vertical Part” and “Molecular Filament Bundle”

As mentioned in the previous section, the MFB, VP, and M0.014-0.054 are clearly seen in Figures 5a and 5b. M0.014-0.054 is located positionally in the interacting area between the MFB and VP. Figures 6a, 6c, and 6e show the enlarged integrated intensity maps around the interacting area in the C32S J=21J=2-1, CH3OH (blended lines around 96.741 GHz), and 28SiO J=21J=2-1 emission lines with the velocity ranges of VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1, respectively. The velocity ranges include that of the MFB. The C32S data includes the single-dish data. On the other hand, Figures 7a, 7c, and 7e show the enlarged integrated intensity maps of the same area in the C32S J=21J=2-1, CH3OH, and 28SiO J=21J=2-1 emission lines with the velocity ranges of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1, respectively. The velocity ranges include those of the VP and M0.014-0.054. As mentioned previously, the SiO and CH3OH emission lines are good tracers of hard (ΔV30\Delta V\simeq 30 km s-1) and mild (ΔV10\Delta V\simeq 10 km s-1) C-shock waves, respectively. Although the filaments of the MFB are clearly seen in Figures 6a, they almost disappear in Figures 6c and 6e except for several faint features. The appearances in both shock tracer lines quite resemble each other. On the other hand, M0.014-0.054 is prominent in Figures 7a, 7c and 7e. The VP is faint in Figures 7c and 7e although it is identified clearly in Figures 7a. These indicate that the whole of the MFB and VP did not suffer from any C-shock wave with ΔV10\Delta V\gtrsim 10 km s-1 but a hard C-shock wave is propagating in M0.014-0.054.

Figures 6b, 6d and 6f show the PV diagrams along the MFB in the C32S J=21J=2-1, CH3OH, and 28SiO J=21J=2-1 emission lines, respectively. The sampling area is shown as the rectangles in Figures 6a, 6c and 6e. In Figures 6b, the MFB is clearly identified as a long inclined ridge with VLSR20V_{\mathrm{LSR}}\sim-20 at the western end to 305030-50 km s-1 at the eastern end. The MFB almost disappears in the SiO and CH3OH emission lines except for the following several faint features. In these PV diagrams, M0.014-0.054 is identified around the angular offset of 80\arcsec-80\arcsec. Figures 7b, 7d, and 7f shows the PV diagrams along M0.014-0.054 in the C32S J=21J=2-1, CH3OH, and 28SiO J=21J=2-1 emission lines, respectively. The sampling area is along M0.014-0.054, which is shown as the rectangles in Figures 7a ,7c, and 7e. In these PV diagrams, M0.014-0.054 is clearly identified as a long feature with no clear velocity gradient although the velocity width becomes large around the east end. The radial velocity difference between M0.014-0.054 on the VP and the MFB is ΔVrad=ΔVcosϕ40\Delta V_{\mathrm{rad}}=\Delta V\cos\phi\sim 40 km s-1, where ΔV\Delta V is the collision velocity in 3D space and ϕ\phi is the angle between the line-of-sight and the axis of the collision. The shock wave induced by the CCC between the VP and MFB is as strong as bringing the enhancement of the SiO molecule in M0.014-0.054.

The “Bridge” features are identified in these PV diagrams of Figures 6 and 7, which are connecting between M0.014-0.054 on the VP and the MFB (arrows, also see a red arrow in Figure 3b). The “Bridge”s also correspond to the faint features seen in Figures 6c and 6e. They contain shocked molecular gas made by hard C-shock wave. Recent hydrodynamical simulation studies have shown that such connecting features in PV diagrams are reproduced as a characteristic feature of a CCC (e.g. Haworth et al. (2015), Haworth et al. (2015b)). Therefore the “Bridge” features support the scenario that the MFB and VP collided each other and made the “Bridge” with the intermediate velocity between those of the MFB and VP, which should be caused by the momentum exchange of the molecular cloud in them. Especially, the feature around the angular offset of 60\arcsec-60\arcsec is more prominent than that around the angular offset of 80\arcsec-80\arcsec around in these emission lines. SiO molecules have been known to be taken in the mantle of interstellar dusts within 10510^{5} years and disappear from interstellar gas as mentioned previously. If the CCC was occurs once in such past, the difference between the“Bridge”s may be made by the elapsed time from the collision.

However, if the MFB and VP collided with the proper motion velocity of V=40×tanϕV=40\times\tan\phi km/s before 10510^{5} yr, one molecular cloud is now separating from another molecular cloud by Δd4.1×tanϕ\Delta d\sim 4.1\times\tan\phi pc. The separation is corresponding to the angular separation of Δθ0.03×tanϕ\Delta\theta\sim 0.03^{\circ}\times\tan\phi at the distance of 8 kpc, which can be easily distinguished in this observation although ϕ\phi is unknown. However, such large angular separations are not observed in the “Bridge”s connecting two clouds. There are at least two possibilities to explain it. As will be mentioned in Section 7, mGauss magnetic fields run along these filaments. In this case, the molecular cloud can move along the filaments but cannot move crossing them. Therefore, the molecular cloud in the Bridge would be fixed around the collision area by the magnetic field which are probably tangled by the collision. On the other hand, 10510^{5} yr is the upper limit of the SiO molecules surviving in the interstellar space before re-absorption into dust and not the elapsed time from the CCC. Therefore, this may suggest the possibility that the HMCs are formed in a shorter time than the surviving upper limit of the SiO molecule although it is required that the angle ϕ\phi is somewhat small.

Refer to caption
Figure 6: a Enlarged integrated intensity map around the ”Molecular Filament Bundle” in the C32S J=21J=2-1 emission line with the velocity ranges of VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1. b PV diagram along the ”Molecular Filament Bundle” in the C32S J=21J=2-1 emission line. The sampling area is shown as the rectangles in a. The C32S maps are after the combining with the single-dish data. c Enlarged integrated intensity map in the CH3OH emission of VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1. d PV diagram along the ”Molecular Filament Bundle” in the CH3OH emission line. The sampling area is shown as the rectangles in c. There is a faint component with ΔV50\Delta V\sim-50 km s-1 from M0.014-0.054. This is the contamination of the CH3OHJKa,Kc=21,111,0J_{K_{a},K_{c}}=2_{1,1}-1_{1,0}E emission line at 96.75550796.755507 GHz. e Enlarged integrated intensity map in the 28SiO J=21J=2-1 emission of VLSR=10V_{\mathrm{LSR}}=10 to 3535 km s-1. f PV diagram along along the ”Molecular Filament Bundle” in the 28SiO J=21J=2-1 emission line. The sampling area is shown as the rectangles in e.
Refer to caption
Figure 7: a Enlarged integrated intensity map around M0.014-0.054 in the C32S J=21J=2-1 emission line with the velocity ranges of 𝐚VLSR=40{\bf a}~{}V_{\mathrm{LSR}}=-40 to 1010 km s-1. The cross shows the position of the collisionally excited CH3OH maser emission at 36.2 GHz in the same velocity range (Yusef-Zadeh et al., 2013). b PV diagram along M0.014-0.054 in the C32S J=21J=2-1 emission line. The sampling area is shown as the rectangles in a. The C32S maps are after the combining with the single-dish data. The cross shows the position and velocity of the collisionally excited CH3OH maser emission at 36.2 GHz (Yusef-Zadeh et al., 2013). c Enlarged integrated intensity map in the CH3OH emission of 𝐚VLSR=40{\bf a}~{}V_{\mathrm{LSR}}=-40 to 1010 km s-1. d PV diagram along M0.014-0.054 in the CH3OH emission line. The sampling area is shown as the rectangles in c. There is a faint component with ΔV50\Delta V\sim-50 km s-1 from M0.014-0.054. This is the contamination of the CH3OHJKa,Kc=21,111,0J_{K_{a},K_{c}}=2_{1,1}-1_{1,0}E emission line at 96.75550796.755507 GHz. e Enlarged integrated intensity map in the 28SiO J=21J=2-1 emission of 𝐚VLSR=40{\bf a}~{}V_{\mathrm{LSR}}=-40 to 1010 km s-1. f PV diagram along M0.014-0.054 in the 28SiO J=21J=2-1 emission line. The sampling area is shown as the rectangles in e.

The crosses in Figures 7a, 7c, and 7e show the position, l=0.02055b=0.05013l=0^{\circ}.02055~{}b=-0^{\circ}.05013, of the CH3OH J=4130J=4_{-1}-3_{0} maser emission spot at 36.2 GHz in the velocity range of V=0.0±8.3V=0.0\pm 8.3 km s-1 (Yusef-Zadeh et al. (2013)). The maser emission spot is located in the vicinity of the strongest intensity peak of the maps. While the crosses in Figures 7b, 7d, and 7f show the position and velocity of the CH3OH J=4130J=4_{-1}-3_{0} maser emission spot in the PV diagrams. The velocity range of the maser emission spot overlaps with the positive velocity edge of M0.014-0.054. Because the CH3OH maser emission is known to be excited collisionally by shock wave propagating in molecular cloud, these facts indicate that the CCC between the MFB and VP is on-going.

Figures 6b and 7b also shows the large feature with VLSR20V_{\mathrm{LSR}}\sim 20 to 8080 km s-1, which has been mentioned in Section 4.2. This feature seems to be made of several ridge-like features. Although the ridge-like features are also identified in Figures 6d, 6f, 7d, and 7f, the spines of the ridges seem to be enhanced. In addition, a faint feature is also identified around the angular extent of M0.014-0.054 and velocity range of VLSR=50V_{\mathrm{LSR}}=-50 to70-70 km s-1 in the PV diagrams of the CH3OH emission line (see Figures 6d and 7d). The feature is the contamination by the CH3OH JKa,Kc=21,111,0EJ_{K_{a},K_{c}}=2_{1,1}-1_{1,0}E emission line (96.75550796.755507 GHz) because this feature is not detected in the C32S J=21J=2-1 emission line (see Figures 6b and 7b).

6 “Hot Molecular Core”s in the Interacting Area

6.1 Identification of Hot Molecular Cores in M0.014-0.054

Figure 8a shows the continuum map at 86 GHz of M0.014-0.054. The noise level of the map is 1σ=0.011\sigma=0.01 K in TBT_{\mathrm{B}}. M0.014-0.054 is detected as a series of the compact objects, A, B, C, D, and E. The eastern two objects are prominent, which correspond to “A” and “B” shown in Figures 5a and 5b. The cross in the map indicates the position, l=0.02055b=0.05013l=0^{\circ}.02055~{}b=-0^{\circ}.05013, of the CH3OH J=4130J=4_{-1}-3_{0} maser emission spot at 36.2 GHz in the velocity range of V=0.0±8.3V=0.0\pm 8.3 km s-1 (Yusef-Zadeh et al. (2013)). This is located in the object A. The mean brightness temperatures of the continuum emission in the objects A and B are T¯B,cont0.030±0.002\bar{T}_{\mathrm{B,cont}}\sim 0.030\pm 0.002 K and T¯B,cont0.020±0.001\bar{T}_{\mathrm{B,cont}}\sim 0.020\pm 0.001 K, respectively. The integration areas are shown as the circles in the map. The brightness temperature ratio of the recombination line to continuum emission is given by the formula, TB,line/TB,cont3×104[Te/K]1.65[ν/GHz]1.1T_{\mathrm{B,line}}/T_{\mathrm{B,cont}}\sim 3\times 10^{4}[T_{\mathrm{e}}/\mathrm{K}]^{-1.65}[\nu/\mathrm{GHz}]^{1.1} in the assumption of LTE and optically thin condition (e.g. Mezger & Henderson (1967)). Therefore, the ratio at 86 GHz of ionized gas with up to Te1×104T_{\mathrm{e}}\lesssim 1\times 10^{4} K is estimated to be larger than unity, T(H42α)B/TB,cont>1T\mathrm{{}_{B}(H42\alpha)}/T_{\mathrm{B,cont}}>1. The mean brightness temperatures of the H42α\alpha recombination line are expected to be T¯(H42α)B>0.030\bar{T}\mathrm{{}_{B}(H42\alpha)}>0.030 K in the object A and T¯(H42α)B>0.020\bar{T}\mathrm{{}_{B}(H42\alpha)}>0.020 K in the object B, respectively. Figure 8b shows the integrated intensity map with the velocity range of 25-25 to 0 km s-1 of the H42α\alpha recombination line. The velocity width is as large as that of the recombination line from the ionized gas at Te1×104T_{\mathrm{e}}\sim 1\times 10^{4} K. The upper limits of the mean brightness temperatures in the H42α\alpha recombination line are 5σ=0.0195\sigma=0.019 K at the object A and 5σ=0.0195\sigma=0.019 K at the object B, respectively. The integration areas are the same as those mentioned above. The objects A and B are not detected in the map. There is no ionized gas in the objects A and B of M0.014-0.054 above the detection limit of our observation. In the cases of the objects C, D, and E, we cannot exclude the existence of ionized gas by this procedure because of their weak intensities. However, the low TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2] ratios suggest that the objects C, D, and E have no ionized gas as will be discussed in Subsection 6.4.

There are two possibilities to explain the situation that the ionized gas does not exist although the mm-wave continuum is detected. One possibility is that these continuum emissions are made by the low frequency extension of the dust emission shown in the JCMT map of Figure 2l. As mentioned previously, the continuum emissions of the objects A and B at 86 GHz correspond to “Double peaks” of the 850μ850~{}\mum dust emission with T¯B0.3\bar{T}_{\mathrm{B}}\sim 0.3 K (see Figure 2l). If the dust β\beta is 2\sim 2, the mean brightness temperatures at 86 GHz of the objects A and B are consistent with the dust emission at the low frequency. In addition, the objects C, D, and E also have the corresponding components in the 850μ850~{}\mum dust emission (see Figure 2l). Another possibility is that this emission is an artifact by contamination from other molecular emission lines. If so, the appearance of the 86 GHz continuum emission in Figure 8a should resemble those in the other molecular emission lines. However, they do not always resemble that of the 86 GHz continuum emission. Therefore, the second possibility seems unlikely. These will be discussed in the followings.

Refer to caption
Figure 8: a Continuum map of M0.014-0.054 at 86 GHz. The first contour level and contour interval are both 0.014 K. The cross indicates the position, l=0.02055b=0.05013l=0^{\circ}.02055~{}b=-0^{\circ}.05013, of the CH3OH J=4130J=4_{-1}-3_{0} maser emission spot at 36.2 GHz in the velocity range of V=0.0±8.3V=0.0\pm 8.3 km s-1 (Yusef-Zadeh et al. (2013)). b Integrated intensity map with the velocity range of 25-25 to 0 km s-1 (1σ=0.0291\sigma=0.029 K) of the H42α\alpha recombination line. The continuum map at 86 GHz (contours) is superimposed on it. c Integrated intensity map of M0.014-0.054 in the 32SO N,J=2,21,1N,J=2,2-1,1 emission line with the velocity ranges of 𝐚VLSR=25{\bf a}~{}V_{\mathrm{LSR}}=-25 to 0 km s-1. d Integrated intensity map in the 34SO N,J=2,31,2N,J=2,3-1,2 emission line with the velocity ranges of 𝐚VLSR=25{\bf a}~{}V_{\mathrm{LSR}}=-25 to 0 km s-1. e Integrated intensity map in the HCOOH JKa,Kc=41,431,3J_{K_{a},K_{c}}=4_{1,4}-3_{1,3} emission line. f Integrated intensity map in the CH3CHO JKa,Kc=52,342,2EJ_{K_{a},K_{c}}=5_{2,3}-4_{2,2}E emission line.

Figures 8c and 8d show the enlarged integrated intensity maps of M0.014-0.054 in the SO and 34SO emission lines with the velocity ranges of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1, respectively. As mentioned previously, the SO and 34SO emission lines are good tracers of the HMC. A compact feature in the SO and 34SO emission lines is located in the object A. In contrast, these emissions are associated with the object B but not centered in it. Or rather, these emissions traces the southern limb of the object. These suggest that the object B has a hollow structure in these molecules. The objects A and B would be HMCs despite some different structures. The appearances in these emission lines resemble those in the HC15N emission line which is shown in Figure 9e. Chemical models of HMCs suggest that the HC15N emission line is enhanced in HMCs with the age of 105\gtrsim 10^{5} years (e.g. see Figure 11 in Stéphan et al. (2018)). Therefore star formation activity may have fairly advanced in the HMCs. In addition, faint features in the SO, 34SO and HC15N emission lines are also detected around the objects C, D, and E. This may suggest that these objects are in the similar evolutional stage to the objects A and B.

Figures 8e and 8f show the enlarged integrated intensity maps of M0.014-0.054 in the HCOOH and CH3CHO(96.475523GHz) emission lines with the velocity ranges of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1, respectively. Although the emissions of the HCOOH and CH3CHO emission lines are very weak, they are also centered in the object A. The association with the other objects is marginal. Figure 10c shows the wide field peak intensity map of the whole area in the HCOOH emission line. The emissions of HCOOH are identified only on several spots in the 50MC except for in M0.014-0.054. These spots seem to correspond to HMCs identified by Miyawaki et al (2020). The concentration to HMCs of the HCOOH emission should be remarkable. The HCOOH emission line may be a good tracer of HMCs. Although the distribution in the CH3CHO emission line is not clear in the 50MC because of the contamination of the C34S emission line, the similar concentration might be seen.

Refer to caption
Figure 9: Enlarged integrated intensity map of M0.014-0.054 of the C32S J=21J=2-1 emission line with the velocity ranges of 𝐚VLSR=25{\bf a}~{}V_{\mathrm{LSR}}=-25 to 0 km s-1. Contours show the continuum map at 86 GHz. The first contour level and contour interval are both 0.014 K. b Map of the C34S J=21J=2-1 emission line. c Map of the H13CO+ J=10J=1-0 emission line. d Enlarged integrated intensity map of the HN13J=10F=11J=1-0~{}F=1-1 emission line. e Map of the HC15J=10J=1-0 emission line. f Map of the HOCO+ JKa,Kc=40,430,3J_{K_{a},K_{c}}=4_{0,4}-3_{0,3} emission line.
Refer to caption
Figure 9: (Continued.) g Map of the CH3OH blended emission lines around 96.741 GHz. Contours show the continuum map at 86 GHz. The first contour level and contour interval are both 0.014 K. h Map of the CH3OH emission lines at 96.583 GHz. i Map of the 28SiO J=21J=2-1 emission line. j Map of the 29SiO J=21J=2-1 emission line.

6.2 Physical Properties of Hot Molecular Cores

Figures 9a, 9b, 9c, and 9d show the enlarged integrated intensity maps of M0.014-0.054 in the C32S, C34S, H13CO+, and HN13C emission lines with the velocity ranges of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1. The image in the C32S emission line has a ridge-like feature along the distribution of the objects A-E in the continuum map, which does not resemble the images in the other emission lines as mentioned previously. The image would suffer from large optical thickness of the C32S emission line. The images in the C34S, H13CO+ and HN13C emission lines resemble each other. These emissions have strong intensity peaks centered at the object A. Meanwhile, the emissions trace the southern limb of the object B rather than the continuum emission peak. There are also their emissions associated with the objects C, D, and E. These similar images would be caused by tracing the similar number density of Hydrogen molecules, n(H2)105n({\mathrm{H}}_{2})\sim 10^{5} cm-3. Figure 9e shows the map in the HC15N emission line with the same velocity range. The critical density of this emission line, n(H2)crit107n({\mathrm{H}}_{2})_{\mathrm{crit}}\sim 10^{7} cm-3, is higher than those of the others. This emission becomes strong at the object A, indicating that the object is filled by a very dense molecular cloud. While, although there is the emission along the south limb of the object B, it dose not centered at the object. The center part of the object B is not filled by such very dense molecular cloud.

The LTE molecular cloud masses of the objects A and B are estimated from the C32S and C34S emission line observations using the procedure shown in Section 4.4. The mean intensity ratios are calculated to be R6.3R\sim 6.3 in the object A and R8.1R\sim 8.1 in the object B, respectively. Then the optical thicknesses of the C32S emission line are τ3.8\tau\sim 3.8 in the object A and τ2.7\tau\sim 2.7 in the object B, respectively. While those of the C34S emission line are τ3.8/22.35=0.17\tau\sim 3.8/22.35=0.17 in the object A and τ2.7/22.35=0.12\tau\sim 2.7/22.35=0.12 in the object B, respectively. Because the C32S emission line of the object is fairly thick, the LTE molecular cloud masses of the objects A and B are estimated from the C34S emission line data using the optical thickness correction factors of τ1eτ1.09\frac{\tau}{1-e^{-\tau}}\sim 1.09 and 1.06\sim 1.06, respectively. The integration velocity range is VLSR=40V_{\mathrm{LSR}}=-40 to 1010 km s-1. The LTE molecular cloud masses of the objects A and B are MLTE=1.3×104(Tex80)MM_{\mathrm{LTE}}=1.3\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})M_{\odot} and MLTE=1.2×104(Tex80)MM_{\mathrm{LTE}}=1.2\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})M_{\odot}, respectively. Although these masses are consistent with the LTE masses of the HMCs detected in the 50MC, M¯LTE8×103(Tex80)M\bar{M}_{\mathrm{LTE}}\sim 8\times 10^{3}(\frac{T_{\mathrm{ex}}}{80})M_{\odot} (Miyawaki et al. (2020)), they are much larger than those of the molecular cloud cores detected in the 50MC, M¯LTE2×102(Tex50)M\bar{M}_{\mathrm{LTE}}\sim 2\times 10^{2}(\frac{T_{\mathrm{ex}}}{50})M_{\odot} (Uehara et al. (2019)).

The mean physical radius of the object A is estimated to be RHWHM=0.33R_{\mathrm{HWHM}}=0.33 pc as the half width at half maximum (HWHM) by 2-dimensional Gaussian fit for the C34S image. The FWHM velocity width is estimated to be VFWHM=8.5±0.7V_{\mathrm{FWHM}}=8.5\pm 0.7 km s-1 by 1-dimensional Gaussian fit for the the C34S line profile. Using the following formula for a spherical object,

Mvir[M]=210×ΔVFWHM[kms1]2R[pc]M_{\mathrm{vir}}[M_{\odot}]=210\times\Delta V_{\mathrm{FWHM}}[\mathrm{kms}^{-1}]^{2}R[\mathrm{pc}] (5)

the virial mass of the object A is estimated to be Mvir=5.0×103MM_{\mathrm{vir}}=5.0\times 10^{3}M_{\odot}. It is difficult to estimate the virial mass of the object B using this formula since the shape of the C34S image is not spherical (see Figure 9b). Because the CCC is ongoing around the object A as mentioned previously, the observed velocity width may be widened, and the derived viral mass should be on the upper limit. Therefore, the virial parameter of the object A is estimated to be Mvir/MLTE1M_{\mathrm{vir}}/M_{\mathrm{LTE}}\lesssim 1. Although the fractional abundance of the CS molecule and the excitation temperature of the CS emission line have large ambiguity, the object A could be bound gravitationally. These masses are also summarized in Table 2.

The flux density at 850 μ\mum in the object A is Sν240S_{\nu}\sim 240 Jy (see Figure 2l). Because the beam size of JCMT is not sufficiently small to resolve the structure seen in the map and there are other components in the line-of-sight (see Figures 3 and 4), the flux density should be the upper limit. Assuming that the dust temperature is 33 K (see Fig. 3 in Etxaluze et al. (2016)) and the gas-to-dust ratio is 100, the gas mass is estimated to be M7×103MM\sim 7\times 10^{3}M_{\odot}. This mass is consistent with the LTE molecular cloud mass and larger than the virial mass. This also indicates that the object A is bound gravitationally.

Refer to caption
Figure 10: a and b Peak intensity maps of the field including the 50MC and Sgr A in the HOCO+ emission line. c, d, and e Peak intensity maps in the HCOOH, C34S, and H13CO+ emission lines, respectively. The velocity range of the maps is VLSR=150V_{\mathrm{LSR}}=-150 to 150150 km s-1. The FWHM beam sizes are shown in the lower left corners. The contours in a, c, d and e show the continuum emission at 86 GHz for comparison. The contour levels are (1,2,8,and32)×3.5(1,2,8,\mathrm{and}~{}32)\times 3.5 K. The contours in b show the continuum emission at 1.44 GHz for comparison (Yusef-Zadeh, Morris, & Chance (1984)). The contour levels are (1,2,3,4,5,6,7,8,9,10,15,20,30,40,and50)×24(1,2,3,4,5,6,7,8,9,10,15,20,30,40,\mathrm{and}~{}50)\times 24 K.

6.3 Chemical Properties of Hot Molecular Cores

6.3.1 HOCO+ Ion

The HOCO+ emission line has been detected toward the molecular clouds in the CMZ including the 50MC (e.g. Minh et al. (1991)). Figure 9f shows the map of M0.014-0.054 in the HOCO+ emission line with the velocity range of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1. This emission is clearly centered on the object A. The distribution closely resembles those of the C34S and H13CO+ emission lines (see Figures 9b and 9c). The C34S emission line is considered to be a faithful tracer of the dense molecular cloud with n(H2)crit105n({\mathrm{H}}_{2})_{\mathrm{crit}}\sim 10^{5} cm-3. The effective critical density of the H13CO+ emission line is similar to that of the C34S emission line. Although this similarity is probably why the distribution of the C34S and H13CO+ emission lines resemble each other, it cannot be explained why that of the HOCO+ emission line resembles them.

Figure 10a shows the wide field peak intensity map of the whole area of this observation, which includes the 50MC and Sgr A, in the HOCO+ emission line. The velocity range is VLSR=150V_{\mathrm{LSR}}=-150 to 150150 km s-1. The emission in M0.014-0.054, especially the object A, is the most prominent in the velocity range and area although some emissions are identified in the 50MC. The distribution might originate with the spatial variation of the fractional abundance of the HOCO+ ion. Figures 10d and 10e show the wide field peak intensity maps in the C34S and H13CO+ emission lines, respectively. In the C34S emission line, the 50MC is the most prominent molecular cloud in the area, although M0.014-0.054 is relatively inconspicuous (also see Figure 5 in Uehara et al. (2019)). On the other hand, the distribution of the H13CO+ emission line in M0.014-0.054 is more conspicuous than that in the 50MC. Although the effective critical densities of these emission lines are similar, the distributions have distinctive difference. The distribution of the HOCO+ emission line relatively resembles that of the H13CO+ emission line.

There are some possibilities of the formation of the HOCO+ ion which could make the distributions mentioned above. An ion-molecule reaction;

HCO++OHHOCO++H,\mathrm{HCO}^{+}+\mathrm{OH}\rightarrow\mathrm{HOCO}^{+}+\mathrm{H},

might be possible (e.g. Fontani, et al. (2018)) because OH radical should be abundant and the gas kinetic temperature is high in the Sgr A region. Because CO molecule is also abundant in the region, the reaction consuming the HOCO+ ion,

HOCO++COHCO++CO2,\mathrm{HOCO}^{+}+\mathrm{CO}\rightarrow\mathrm{HCO}^{+}+\mathrm{CO}_{2},

is also expected in the region (e.g. Sakai et al. (2008)). In this case, the distribution of the HOCO+ emission line would resemble that of the H13CO+ emission line.

The HOCO+ ion may be formed by the protonated reaction of CO2 molecule. H+3{}_{3}^{+} ion is abundant by the abundant cosmic ray in the Sgr A region (e.g. Oka et al. (2019)). CO2 molecules are evaporated extensively from the dust mantle by high gas kinetic temperature in M0.014-0.054 (TK80T_{\mathrm{K}}\simeq 80 K in Ao et al. (2013)). Therefore an ion-molecule reaction forming the HOCO+ ion,

H3++CO2HOCO++H2,\mathrm{H}_{3}^{+}+\mathrm{CO}_{2}\rightarrow\mathrm{HOCO}^{+}+\mathrm{H}_{2},

is also expected in the region (e.g. Sakai et al. (2008)). In this case, it is not necessary that the distribution of the HOCO+ emission line resemble that of the H13CO+ emission line. Note that these reactions are not mutually exclusive and other reactions forming and destroying HOCO+ ion are also possible.

Figure 10b shows the comparison between the peak intensity map of the HOCO+ emission line and 1.44 GHz continuum emission (Yusef-Zadeh, Morris, & Chance (1984)). The features in the 50MC observed in the HOCO+ emission line are located only along the southeast edge of Sgr A East SNR. This suggests that the physical interaction between the 50MC and Sgr A East SNR forms the HOCO+ ions or at least assists the formation. There are two possibilities explaining the scenario. The first one may be concerned in shock chemistry by the SNR. However, it is hard to specify the mechanism of the formation immediately. In the second one, abundant OH radicals made by cosmic ray of the SNR may form the ions through the reaction mentioned above (also see Tielens (2013)).

On the other hands, the formation scenario of the features observed in M0.014-0.054 is not concerned probably in SNRs because there is no known SNR around it. However, the low frequency continuum emission associated with the south edge of M0.014-0.054 is seen. This is identified more clearly at 330 MHz as mentioned in Sec.2 (see Figure 1d). A possibility explaining the continuum emission would be that the CCC between the VP and MFB accelerates cosmic ray. The strong magnetic field (see Section 7) and large collision velocity (see Section 5) may enable such acceleration around M0.014-0.054 although it has not been observed in the disk region. The accelerated cosmic ray would increase the abundance of OH radicals simultaneously (also see Tielens (2013)). Therefore the HOCO+ ions would be formed through second scenario mentioned above.

6.3.2 CH3OH and SiO Molecules

As mentioned previously, the C-type shock wave propagating in the molecular cloud enhances the abundances of the CH3OH and SiO molecules. Figures 9g and 9h show the enlarged integrated intensity maps of M0.014-0.054 in the CH3OH (blended lines around 96.741 GHz) and CH3OH (97.583 GHz) emission lines with the velocity ranges of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1, respectively. However, the higher JJ transition lines of CH3OH molecule, JKa,Kc=62,571,7J_{K_{a},K_{c}}=6_{-2,5}-7_{-1,7} E at 85.568 GHz , JKa,Kc=72,663,3J_{K_{a},K_{c}}=7_{2,6}-6_{3,3} A at 86.616 GHz, JKa,Kc=216,16225,17J_{K_{a},K_{c}}=21_{6,16}-22_{5,17} A and 216,15225,1821_{6,15}-22_{5,18} A++ at 97.679 GHz, JKa,Kc=246,19237,16J_{K_{a},K_{c}}=24_{6,19}-23_{7,16} A and 246,18237,1724_{6,18}-23_{7,17} A++ at 98.030 GHz, are not detected in this observation. Figures 9g and 9h do not well resemble each other although these are the maps in the CH3OH emission lines with J=2J=2.

Figures 9i and 9j show the enlarged integrated intensity maps in the 28SiO and 29SiO emission lines with the velocity ranges of VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1, respectively. These do not well resemble each other. On the other hand, Figures 9g and 9i resemble each other and Figures 9h and 9j resemble each other in spite of different molecules. The intensity ratio between the object A and the others in Figure 9g is smaller than that in Figure 9h although the intensity in the former is larger than that in the latter. The difference between the two CH3OH maps is similar to that between the maps of the 28SiO and 29SiO emission lines, which are the major and minor isotopes (see Figures 9i and 9j). These would be explained by that the CH3OH (96.741 GHz) and 28SiO emission lines are fairly optically thick around the object A but the CH3OH (97.583 GHz) and 29SiO emission lines are optically thin. Figure 9h and 9j would indicate faithfully the abundance enhancements of these molecules. Therefore, we conclude that the shocked molecular gas is enhanced strongly in the object A although it is somewhat enhanced in other components. This is consistent with the detection of the Class-I methanol maser only in the object A as mentioned above.

Refer to caption
Figure 11: a Integrated intensity map of of M0.014-0.054 in the C2N=10J=3/21/2F=10N=1-0~{}J=3/2-1/2~{}F=1-0 emission line. The velocity range is VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1. The unit of the color bar is mean brightness temperature. Contours show the continuum map at 86 GHz. The first contour level and contour interval are both 0.014 K. b Integrated intensity map of M0.014-0.054 in the c-C3H2 JKa,Kc=21,210,1J_{K_{a},K_{c}}=2_{1,2}-1_{0,1} emission line. c The ratio map of TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2] in M0.014-0.054. d Integrated intensity map of M0.014-0.054 in the CCS N,J=7,66,5N,J=7,6-6,5 (86.181413 GHz) emission line. The velocity range is VLSR=18V_{\mathrm{LSR}}=-18 to 8-8 km s-1.
Refer to caption
Figure 12: a Spectra of the C2H and c-C3H2 JKa,Kc=21,210,1J_{K_{a},K_{c}}=2_{1,2}-1_{0,1} emission lines of the object A. The sampling area is shown in Figure 11 (circles). b the C2H and c-C3H2 emission lines of the object B. The sampling area is shown in Figure 11 (circles).

6.4 Embedded Stars in the “Hot Molecular Core”s

The C2H and c-C3H2 molecules would fairly survive even in the molecular clouds exposed to UV radiation than others (e.g. Cuadrado et al. (2015)). Figures 11a and 11b show the integrated intensity maps of M0.014-0.054 in the C2H N=10J=3/21/2F=21N=1-0~{}J=3/2-1/2~{}F=2-1 and c-C3H2 JKa,Kc=21,210,1J_{K_{a},K_{c}}=2_{1,2}-1_{0,1} emission lines. The velocity ranges are both VLSR=25V_{\mathrm{LSR}}=-25 to 0 km s-1. The distributions in these emission lines resembles those in the C34S, H13CO+ and HN13C emission lines. The higher JJ emission line, c-C3H2 JKa,Kc=43,242,3J_{K_{a},K_{c}}=4_{3,2}-4_{2,3}, at 85.656485.6564 GHz is not detected. Other small hydrocarbon molecules, l-C3H (97.996GHz and 98.012 GHz), C4H (96.478 GHz), and C5H (97.863 GHz) emission lines are not also detected. The transition of c-C3H is not in our observation frequency range.

Figure 11c shows the ratio map of TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2] in M0.014-0.054. In addition, figures 12a and 12b show the spectra of the C2H N=10J=3/21/2F=21N=1-0~{}J=3/2-1/2~{}F=2-1 and c-C3H2 JKa,Kc=21,210,1J_{K_{a},K_{c}}=2_{1,2}-1_{0,1} emission lines on the objects A and B, respectively. The sampling areas are shown in Figures 11a and 11b as circles. The ratio in the object A is TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2]=1.73±0.03=1.73\pm 0.03. However, the ratio is seen to be inhomogeneous in the object. The ratio of the southeastern half is TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2]=2.06±0.05=2.06\pm 0.05 (see Figure 11c). The ratio of the northwestern half is TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2]=1.47±0.04=1.47\pm 0.04 (see Figure 11c). On the other hand, the ratio in the object B is TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2]=1.53±0.02=1.53\pm 0.02(see Figure 11c). The ratios of the southern half and northern half are TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2]=1.31±0.03=1.31\pm 0.03 and =1.96±0.05=1.96\pm 0.05. Although the integrated intensities of the emission lines in the northern half are weaker than those in the southern half, the ratio in the northern half is larger than that in the southern half. The ratios in the objects C, D, and E are TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2]=1.40±0.04=1.40\pm 0.04, 1.84±0.041.84\pm 0.04, and 1.47±0.041.47\pm 0.04, respectively(see Figure 11c). The ratios in the objects C, and E are about 1.5\sim 1.5. These are located on the intensity ridge of the molecular cloud in M0.014-0.054 (see Figure 11a and 11b). The ratio of the object D is slightly higher than the others.

Although the photodissociation potential of the C2H molecule, 4.9 eV (C2H \rightarrow C2 + H), is similar to that of the c-C3H2 molecule, 4.4 eV (c-C3H2 \rightarrow C3H + H), the ionization potential of the C2H molecule, 11.4 eV, is fairly higher than that of the c-C3H2 molecule, 9.2 eV (van Hemert & van Dishoeck (2008)). The difference between the ionization potentials may raise the viability of C2H molecule in the UV irradiated condition compared with that of the c-C3H2 molecule and increase the abundance ratio of XX[C2H]/XX[c-C3H2]. Consequently, TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2] ratio would increase.

The high ratios measured in the southeastern half of the object A and the northern half of the object B are distinct from those in the northwestern half of the object A, the southern half of the object B, the object C and the object E. The TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2] ratios have been reported to be increased in UV irradiated regions, for example, 3\sim 3 at the Horse head PDR (Teyssier et al. (2004)). The measured high ratios may be consistent with that of the Horse head PDR. Similar UV irradiated condition is expected in both M0.014-0.054 and the Horse head PDR. Because the object B is adjacent to the object A, these are expected to be in the similar UV irradiated condition. If the UV photons with softer spectrum, for example hνtyp.45h\nu_{\mathrm{typ.}}\sim 4-5 eV, are dominated in the objects, a slight difference between the photodissociation potentials of these molecules would decrease the abundance ratio between XX[C2H] and XX[c-C3H2]. Consequently, TBT_{\mathrm{B}}[C2H]/TBT_{\mathrm{B}}[c-C3H2] ratio is decreased. These suggest that the lower ratios of 1.5\sim 1.5 are mainly determined by the embedded stars or protostars.

Figure 11d shows the enlarged integrated intensity map of M0.014-0.054 in the CCS N,J=7,66,5N,J=7,6-6,5 (86.181413 GHz) emission line with the velocity range of VLSR=18V_{\mathrm{LSR}}=-18 to 8-8 km s-1. The same velocity range as in other panels cannot be set because this emission line is at the edge of the observation frequency range. The faint mission is centered in the object A. The emission in the object B is not centered in it but traces the southern limb of the object. The faint emission is probably associated with the objects C, D and E. The objects A and B are detected in the dust continuum emission but not detected in the recombination line as shown in Subsection 6.1. The star formation activities had started and the surfaces of the embedded stars had been heated at least to several 100 K. However, the existence of CCS molecules suggests that the chemical evolution by the star formation activity in the object A is still in the early stage (e.g. Hirahara et al. (1992)).

On the other hand, the molecular emission lines including the CCS emission line are not centered in the object B. The observed features suggest that CCS molecules are still remained in the hollow structure of the molecular cloud around the embedded star. The star would emit soft UV radiation enough to dissociate the surrounding molecular cloud. This means that the chemical evolution in the object B is advanced comparing with that in the object A. If the star has already been in the main sequence stage, it may be a less massive star than B1 because it does not emit vast Lyman continuum emission enough to ionize the surrounding gas. In addition, if the faint features seen around the objects C, D, and E in the molecular emission lines including the CCS emission line are the remnants of the cradle molecular clouds, the chemical evolutions in these objects may be more advanced comparing with that in the object B.

Refer to caption
Figure 13: The vectors indicate the polarization data (direction of BB and polarization degree) at 350μ350\mum observed by the Caltech Submillimeter Observatory 10-m telescope (Table 4 in Chuss et al. (2003)), whic are overlaid on integrated intensity maps around the “Vertical Part” and M0.014-0.054 in the C32S J=21J=2-1 emission line with the velocity ranges of VLSR=150V_{\mathrm{LSR}}=-150 to 150150 km s-1 (pseudo color) and VLSR=40V_{\mathrm{LSR}}=-40 to 1010 km s-1 (contours), respectively. The first contour and contour interval are both 2525 K km s-1.

7 Magnetic Field in the “Vertical Part”

7.1 Direction of the Magnetic Field

We used polarization data at 350μ350\mum with the Caltech Submillimeter Observatory 10-m telescope (CSO) to obtain the direction of the magnetic field, ϕint\phi_{\mathrm{int}}, in the VP (Table 4 in Chuss et al. (2003)). Figure 13 shows the submillimeter polarization (BB vector) overlaid on the velocity integrated intensity map with the velocity ranges of 150-150 to 150150 km s-1 (pseudo color) and 40-40 to 1010 km s-1 (contours) in the C32S J=21J=2-1 emission line. There are two molecular cloud components on the line of sight, of which one belongs to the VP and another is seen in the velocity range of 3030 to 7070 km s-1 (see Figures 3b and 4b). The observed polarization vectors are the sum of of the intrinsic polarization vectors originated by these components. Although the 50MC is dominant in the western half area, the VP could be identified as the bundle of filaments in the eastern half area. Therefore the observed BB vectors in the eastern half area are expected to be originated mainly by the VP. We consider that the observed BB vectors well trace the magnetic filed lines through the classical Davis-Greenstein mechanism, in which the observed BB vectors are thought to be parallel to the magnetic filed lines. The well-ordered BB vectors suggest that the magnetic filed lines run along the filaments of the VP.

Such poloidal magnetic field (or perpendicular to the Galactic plane) has been usually observed in the non-thermal structures in the Galactic center region (e.g. Yusef-Zadeh, Morris, & Chance (1984), Tsuboi et al. (1986)). Meanwhile the magnetic field in the Galactic center molecular clouds has been unveiled to be mainly troidal (or parallel to the Galactic plane) by IR observations (e.g. Nishiyama et al. (2010)). The poloidal magnetic field in the molecular cloud is an unique case in the Galactic center region.

7.2 Strength of the Magnetic Field

The Chandrasekhar-Fermi method is used to estimate the magnetic field strength from the direction fluctuation of the magnetic filed line and internal gas kinematics (Chandrasekhar & Fermi (1953)). The magnetic field strength on the plane of sight is given by

B=4πρδvδϕ,B_{\bot}=\sqrt{4\pi\rho}\frac{\delta v}{\delta\phi}, (6)

where ρ\rho, δv\delta v, and δϕ\delta\phi are gas density, velocity dispersion, and direction fluctuation of the magnetic field, respectively. This is converted to the following formula for molecular clouds (Crutcher et al. (2004));

B[μGauss]9.3×n(H2)ΔvintΔϕint,B_{\bot}[\mu\mathrm{Gauss}]\simeq 9.3\times\sqrt{n\mathrm{(H_{2})}}\frac{\Delta v_{\mathrm{int}}}{\Delta\phi_{\mathrm{int}}}, (7)

where Δϕint\Delta\phi_{\mathrm{int}} [degree] is the direction fluctuation of the magnetic field, the n(H2)n\mathrm{(H_{2})}[cm-3] is the molecular gas density, which can be estimated by the critical density of the observed molecular line, and Δvint\Delta v_{\mathrm{int}}[km s-1] is the FWHM velocity width of the line profile integrated in the area. The relation between the observed fluctuation, Δϕobs\Delta\phi_{\mathrm{obs}}, and Δϕint\Delta\phi_{\mathrm{int}}, is given by

Δϕint=Δϕobs2ϕerr2,\Delta\phi_{\mathrm{int}}=\sqrt{\Delta\phi_{\mathrm{obs}}^{2}-\phi_{\mathrm{err}}^{2}}, (8)

where ϕerr\phi_{\mathrm{err}} is the mean error of the observation. The Δϕobs\Delta\phi_{\mathrm{obs}} is calculated from the observed magnetic field directions (Chuss et al. (2003)) as the standard deviation. The intrinsic fluctuation of the direction of the magnetic field is estimated to be Δϕint9\Delta\phi_{\mathrm{int}}\sim 9 deg in the VP (Chuss et al. (2003)). The FWHM velocity width of the C34S emission line is derived to be Δvint=13.8±0.8\Delta v_{\mathrm{int}}=13.8\pm 0.8 km s-1 by Gaussian-fit. The emission line is optically thin in the VP as mentioned previously. The sound velocity of C34S molecules is estimated to be cs,C34S=kBTKμC34SmH0.1c_{\mathrm{s,C34S}}=\sqrt{\frac{k_{\mathrm{B}}T_{\mathrm{K}}}{\mu_{\mathrm{C34S}}m_{\mathrm{H}}}}\simeq 0.1 km s-1 at TK=80T_{\mathrm{K}}=80 K where μC34S\mu_{\mathrm{C34S}} is the molecular weight of a C34S molecule, μC34S=46\mu_{\mathrm{C34S}}=46, and mHm_{\mathrm{H}} is the mass of a Hydrogen atom. The broadening of the FWHM velocity width by the sound velocity is negligible. The effective critical density of the C34S emission line is assumed to be n(H2)5×104n\mathrm{(H_{2})}\sim 5\times 10^{4}cm-3 at TK=80T_{\mathrm{K}}=80 K. Consequently, the magnetic field strength in the VP is estimated to be B3B_{\bot}\sim 3mGauss.

The Chandrasekhar-Fermi method is not valid when the magnetic field line is significantly fluctuated by the external perturbations such as HII regions and/or SNRs because the fluctuation of the magnetic field is assumed to be originated only by Alfvèn wave in the method. Because the CMZ is generally filled with the external perturbations, the magnetic field must suffer from such perturbations in varying degrees. The estimated value of the magnetic field strength should have large ambiguity.

7.3 Stability of Molecular Cloud Filaments

The critical mass per unit length of molecular filaments is given by

Mline,crit=2σtot2G,M_{\mathrm{line,crit}}=\frac{2\sigma_{\mathrm{tot}}^{2}}{G}, (9)

where σtot\sigma_{\mathrm{tot}} is the total velocity dispersion which is given by σtot=σnonth2+cs2+12VA2\sigma_{\mathrm{tot}}=\sqrt{\sigma_{\mathrm{nonth}}^{2}+c_{\mathrm{s}}^{2}+\frac{1}{2}V_{\mathrm{A}}^{2}} for a magnetized molecular cloud filament (e.g. Fiege & Pudritz (2000), Arzoumanian et al. (2013)). When the LTE mass per unit length of the filament is larger than this limit, the filament can fragment into molecular cores along the axis of the filament (e.g. Inutsuka & Miyama (1997)). In the case of the VP, the velocity dispersion of the nonthermal motion is calculated to be σnonth=Δvint2cs,C34S28ln2=5.9±0.3\sigma_{\mathrm{nonth}}=\sqrt{\frac{\Delta v_{\mathrm{int}}^{2}-c_{\mathrm{s,C34S}}^{2}}{8ln2}}=5.9\pm 0.3 km s-1. The sound velocity in the VP is estimated to be cs=kBTKμmH0.5c_{\mathrm{s}}=\sqrt{\frac{k_{\mathrm{B}}T_{\mathrm{K}}}{\mu m_{\mathrm{H}}}}\simeq 0.5 km s-1 at TK=80T_{\mathrm{K}}=80 K where μ\mu is the mean molecular weight, μ=2.8\mu=2.8. While the Alfvèn velocity is estimated to be VA=1300Bn(H2)17V_{\mathrm{A}}=\frac{1300B}{\sqrt{n({\mathrm{H_{2}}}})}\simeq 17 km s-1 at the magnetic field strength of B3B\sim 3 mGauss and the molecular gas density of n(H2)5×104n(\mathrm{H_{2}})\sim 5\times 10^{4} cm-3 as mentioned in the previous subsection. The sound velocity is much smaller than the Alfvèn velocity and the velocity dispersion of the nonthermal motion. Therefore, the critical mass per unit length of the VP is estimated to be Mline,crit=2σtot2G8×104M_{\mathrm{line,crit}}=\frac{2\sigma_{\mathrm{tot}}^{2}}{G}\sim 8\times 10^{4} MM_{\odot} pc-1. Because the magnetic field strength is not observed in the MFB, the critical mass per unit length cannot be derived. However, the lower limit can be estimated when the magnetic field is ignored. In the case of the MFB, the velocity dispersion of the nonthermal motion is calculated to be σnonth=5.1±0.3\sigma_{\mathrm{nonth}}=5.1\pm 0.3 km s-1. The critical mass per unit length of the MFB is Mline,crit1×104M_{\mathrm{line,crit}}\gtrsim 1\times 10^{4} MM_{\odot} pc-1.

On the other hands, the LTE mass per unit length of the VP is given by

Mline,LTEMLTEnL,M_{\mathrm{line,LTE}}\sim\frac{M_{\mathrm{LTE}}}{nL}, (10)

where MLTEM_{\mathrm{LTE}} is the LTE molecular cloud mass estimated in Subsection 4.4, MLTE4×104(Tex80)MM_{\mathrm{LTE}}\sim 4\times 10^{4}(\frac{T_{\mathrm{ex}}}{80})M_{\odot}, LL is the the length of the VP, L8L\sim 8 pc, and nn is the number of the observed filaments, n3n\sim 3 (see Figures 3a and 4a). Therefore, the LTE mass per unit length of the VP is estimated to be Mline,LTE2×103(Tex80)M_{\mathrm{line,LTE}}\sim 2\times 10^{3}(\frac{T_{\mathrm{ex}}}{80}) MM_{\odot} pc-1. Using the same procedure, the LTE mass per unit length of the MFB is estimated to be Mline,LTE3×103(Tex80)M_{\mathrm{line,LTE}}\sim 3\times 10^{3}(\frac{T_{\mathrm{ex}}}{80}) MM_{\odot} pc-1. Because the inclination angles of the VP and MFB are not known, these LTE masses per unit length are the upper limits. Both the LTE masses per unit length of the VP and MFB are much smaller than their critical masses per unit length, respectively. The filaments are stable for the fragmentation along the filaments. The star formation activity in the VP and MFB cannot start without any external trigger. In this observation, we found the evidences of the CCC including the “Bridge”s between the VP and MFB, and also found the evidences of the star formation including HMCs in M0.014-0.054, which is located at the intersection between the colliding filaments. Even for stable molecular filaments, the CCC is presumably make a role in the star formation as the external trigger (see also Figure 1 in Inoue &Fukui (2013)). Although there are still many issues in massive star formation, the star formation induced by CCC would be a promising mechanism in the SgrAMC.

Refer to caption
Figure 14: a Schematic display of the VP and MFB before the collision between them. The viewpoint is shown on the upper-left corner. The magnetic fields in the VP and MFB are perpendicular to the Galactic plane and parallel to it, respectively. Although the magnetic field in the MFB is not observed, it is expected that this is parallel to the Galactic plane because the magnetic field in the neighboring molecular clouds is usually along the Galactic plane (e.g. Nishiyama et al. (2010)). The molecular clouds in them are thought to be along the magnetic fields. The molecular clouds in the VP are considered to be confined around the Galactic plane by the gravity although the molecular clouds in the MFB are considered to be able to move along the magnetic fields. b Schematic display of the VP and MFB after the collision. The viewpoints of the left and right panels are along the MFB and VP, respectively. The collision velocity should be larger than Δ40\Delta\gtrsim 40 km s-1 because the SiO molecules seem to be abundant in the colliding area, which are confined in the tangled magnetic field there. In the VP, the molecular clouds should converge kinematically along the magnetic fields and form massive molecular cloud cores around the collision area. Successively, they form high mass protostars in their insides and evolve into HMCs. On the other hand, the molecular clouds in the MFB can move from the colliding area along the magnetic field because there is no confined molecular cloud like in the VP. Therefore the molecular clouds around the colliding area would become less dense than those in the VP and cannot form high mass protostars in their insides.

8 Spatial Structure of the CCC and It’s Possible Scenario

Based on the observation results and discussion mentioned above, we would like to depict the spatial structure of the CCC between the VP and MFB and infer the possible scenario. As mentioned in the previous section, the VP and MFB, except for the colliding area, are stable molecular filaments against gravitational fragmentation, which are approximately perpendicular to the Galactic plane and parallel to it, respectively. The stabilities are secured by the strong magnetic field and/or large velocity dispersion in the filaments.

Figure 14a shows the schematic display of the VP and MFB before the collision. The viewpoint of the panel is along the MFB. The directions of the magnetic fields observed by CSO is considered to indicates the large scale structure of the magnetic fields in the VP because the effective angular resolution is as large as FWHM2+Grid227\arcsec\sqrt{\mathrm{FWHM}^{2}+\mathrm{Grid}^{2}}\sim 27\arcsec (Chuss et al. (2003)). The magnetic field in the VP is approximately perpendicular to the Galactic plane (see Figure 13). The molecular filaments depicted in the CS emission line are along the magnetic field. Note that it has not been clear how the vertical magnetic field are originated and how molecular clouds are took into the VP. On the other hand, those in the MFB are not observed although the molecular filaments are approximately parallel to the Galactic plane (see Figure 5). However, because the magnetic field along the Galactic plane is observed in the neighboring molecular clouds (e.g. Nishiyama et al. (2010)), such magnetic field may exist also in the MFB. Molecular clouds can generally move along the magnetic fields even if they have BB\simmGauss in the MFB. On the other hand, those in the VP may not be such case. As mentioned above, the molecular clouds in the VP would be distributed perpendicular to the Galactic plane. It is impossible that they orbit in Kepler motion keeping this distribution. The portions obeying the Kepler motion don’t feel the gravity although the portions departing from it feel the gravity. If the molecular cloud near the the Galactic plane moves along the magnetic field, it must push up other clouds captured around the Galactic plane by the gravity. Therefore, the molecular clouds in the VP are confined around the Galactic plane by the gravity.

Figure 14b shows the schematic display of the VP and MFB after the collision. The viewpoint of the left panel is along the MFB, which is the same as that in Figure 14a. While, that of the right panel is perpendicular to the Galactic plane, which is along the VP. The collision velocity should be larger than Δ40\Delta\gtrsim 40 km s-1 because the SiO emission line is enhanced around the expected colliding area. The molecular filaments in the VP are seen to get tangled up around M0.014-0.054 although they are lined up perpendicular to the Galactic plane in the farther area from it (see Figures 2 and 3). This probably shows that the magnetic field in the VP got tangled up in the colliding area because the filaments trace the magnetic field. In the VP, the molecular clouds should converge kinematically along the magnetic fields because the magnetic fields are deformed as shown in Figure 14b (C.f. Inoue &Fukui (2013)), and form massive molecular cloud cores around the collision area. Because the shocked molecular gas made by the collision cannot expand crossing the filaments, this should be confined around the colliding area. The molecular cloud cores were massive enough to be bound gravitationally as mentioned in Subsection 6.2. Therefore, they form high mass protostars in their insides and evolve into HMCs. On the other hand, the molecular clouds in the MFB can escape from the colliding area along the magnetic field because there is no confined molecular cloud like in the VP as shown the right panel of Figure 14b. Therefore the molecular clouds around the colliding area would become less dense than those in the VP and cannot form high mass protostars in their insides. This scenario explains that the HMCs are observed only in the VP.

{ack}

This work is supported in part by the Grant-in-Aids from the Ministry of Eduction, Sports, Science and Technology (MEXT) of Japan, No.16K05308 and No.19K03939. We are grateful to Prof. S. Takano at Nihon University for useful discussion. This paper is partly based on observations by the Nobeyama 45-m radio telescope. The telescope is operated by Nobeyama Radio Observatory, a branch of National Astronomical Observatory of Japan (NAOJ). The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation (NSF) of USA operated under cooperative agreement by Associated Universities, Inc.. This paper also makes use of the following ALMA data:ADS/JAO.ALMA#2012.1.00080.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC(Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work also has made use of the JCMT Science Archive (http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/jcmt/). The JCMT has historically been operated by the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council of the United Kingdom, the NRC, and the Netherlands Organization for Scientific Research.

References

  • Ao et al. (2013) Ao, Y. et al.  2013, A&A550, id A135
  • Armijos-Abendaño et al. (2015) Armijos-Abendaño, J., Martín-Pintado, J., Requena-Torres, M. A., Martín, S., & Rodríguez-Franco, A.  2015, MNRAS, 446, 3842
  • Arzoumanian et al. (2013) Arzoumanian, D., André, P., Peretto, N., & Könyves, V.  2013, A&A, 553, id A119
  • Bally et al. (1987) Bally, J., Stark, A. A., Wilson, R. W., & Henkel, C.,  1987, ApJS, 65, 13
  • Boehle et al. (2016) Boehle, A. et al. 2016, ApJ, 830, 17
  • Chandrasekhar & Fermi (1953) Chandrasekhar, S.& Fermi, E.  1953, ApJ, 118, 113
  • Chuss et al. (2003) Chuss, D. T., Davidson, J. A., Dotson, J. L., Dowell, C. D., Hildebrand, R. H., Novak, G., & Vaillancourt, J. E.  2003, ApJ, 599, 1116
  • Crutcher et al. (2004) Crutcher, R. M., Nutter, D. J., Ward-Thompson, D,, & Kirk J. M.,  2004, ApJ, 600, 279
  • Cuadrado et al. (2015) Cuadrado, S., Goicoechea, J. R., Pilleri, P., Cernicharo, J., Fuente, A., & Joblin, C.  2015, A&A, 575, A82
  • Etxaluze et al. (2016) Etxaluze, M., Smith, H. A., Tolls, V., Stark, A. A., & González-Alfonso, E.,  2011, AJ, 142, id 134
  • Fiege & Pudritz (2000) Fiege, J. D., & Pudritz, R. E.  2000, MNRAS, 311, 85
  • Figer et al. (1999) Figer, D. F.; McLean, I.  S., & Morris, M.,  1999, ApJ, 514, 202
  • Figer et al. (2002) Figer, D. F. Najarro, F., Gilmore, D., Morris, et al.  2002, ApJ, 581, 258
  • Fontani, et al. (2018) Fontani, F., Vagnoli, A., Padovani, M., Colzi, L., Caselli, P., & Rivilla, V. M.  2018, MNRAS, 481, 79
  • Fuente, Rodriguez-Franco, & Martin-Pintado (1996) Fuente, A., Rodriguez-Franco, A., & Martin-Pintado, J.  1996, A&A, 312, 599
  • Genzel et al. (1996) Genzel, R., Thatte, N., Krabbe, A., Kroker, H., & Tacconi-Garman, L. E.  1996, ApJ, 472, 153
  • Gusdorf et al. (2008) Gusdorf, A., Cabrit, S., Flower, D.R., & Pineau des Forêts, G.,  2008, A&A, 482, 809
  • Hasegawa et al. (1994) Hasegawa, T., Sato, F., Whiteoak, J.B., Miyawaki, R.,  1994, ApJ. 429, L77
  • Hasegawa et al. (2008) Hasegawa, T., Arai, T., Yamaguchi, N., Sato, F.,  2008, Ap&SS, 313, 91
  • Hartquist et al. (1995) Hartquist, T. W., Menten, K. M., Lepp, S., & Dalgarno, A.  1995, MNRAS, 272, 184
  • Haworth et al. (2015) Haworth, T. J., Tasker, E. J., Fukui, Y. et al.,  2015, MNRAS, 450, 10
  • Haworth et al. (2015b) Haworth, T. J., Shima, K., Tasker, E. J. et al.,  2015, MNRAS, 454, 1634
  • van Hemert & van Dishoeck (2008) van Hemert, M. C. & van Dishoeck, E. F.  2008, ChemPhys, 343, 292
  • Hirahara et al. (1992) Hirahara, Y., Suzuki, H., Yamamoto, S., Kawaguchi, K., Kaifu, N., Ohishi, M., Takano, S., Ishikawa, S., & Masuda, A.  1992, ApJ, 394, 539
  • Inoue &Fukui (2013) Inoue, T. & Fukui, Y.  2013, ApJ, 774, id. L31
  • Inutsuka & Miyama (1997) Inutsuka, S.-i., & Miyama, S. M.  1997, ApJ, 480, 681
  • Jansen et al. (1995) Jansen, D. J., Spaans, M., Hogerheijde, M. R.,&van Dishoeck, E. F.  1995, A&A, 303, 541
  • Jiménez-Serra et al. (2008) Jiménez-Serra, I., Caselli, Martín-Pintado, P.J., & Hartquist, T. W.,  2008, A&A482, 549
  • LaRosa et al. (2000) LaRosa, T. N., Kassim, Namir E., Lazio, T. Joseph W., & Hyman, S. D.  2000, AJ, 119, 207
  • McMullin et al. (2007) McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, Astronomical Data Analysis Software and Systems XVI (ASP Conf. Ser. 376), ed. R. A. Shaw, F. Hill, & D. J. Bell (San Francisco, CA: ASP), 127
  • Mezger & Henderson (1967) Mezger, P. G.& Henderson, A. P.  1967, ApJ, 147, 471
  • Minh et al. (1991) Minh, Y. C., Brewer, M. K., Irvine, W. M., Friberg, P., & Johansson, L. E. B.  1991, A&A, 244, 470
  • Miyawaki et al. (2020) Miyawaki, R., Tsuboi, M., Kitamura, Y., Uehara, K.,& Miyazaki, A.,  2020, in preparation.
  • Morris (1993) Morris, M. R.,  1993, ApJ, 408, 496
  • Morris & Serabyn (1996) Morris, M.,& Serabyn, E.  1996, AAR&A, 34, 645
  • Nagy et al. (2015) Nagy, Z., Ossenkopf, V., Van der Tak, F. F. S., et al.  2015, A&A, 578, A124
  • Nishiyama et al. (2010) Nishiyama, S., Hatano, H., Tamura, M., et al.  2010, ApJ, 722, L23
  • Oka et al. (1998) Oka, T., Hasegawa, T., Sato, F., Tsuboi, M., & Miyazaki, A.,  1998, ApJS, 118, 455
  • Oka et al. (2019) Oka, T., Geballe, T. R., Goto, M., Usuda, T., Benjamin, J. M., & Indriolo, N.  2019, ApJ, 883, id 54
  • Pierce-Price et al. (2000) Pierce-Price, D. et al.,  2000, ApJ, 545, L121
  • Sakai et al. (2008) Sakai, N., akai, T., Aikawa, Y., Yamamoto, S.  2008, ApJ, 675, L89
  • Serabyn & Güsten (1987) Serabyn, E. & Güsten, R.   1987, A&A, 184,133
  • Schmidt, Zack, & Ziurys (2018) Schmidt, D. R., Zack, L. N. & Ziurys, L. M.  2018, ApJ, 864, id. L31
  • Stéphan et al. (2018) Stéphan, G., S., Schilke, P., Le Bourlot, J., Schmiedeke, A., Choudhury, R., Godard, B., & Sánchez-Monge, Á.  2018, A&A, 617, 60
  • Teyssier et al. (2004) Teyssier, D., Fossé, D., Gerin, M., Pety, J., Abergel, A.& Roueff, E.  2004, A&A, 417, 135
  • Tielens (2013) Tielens, A. G. G. M.,  2013, Rev. Mod. Phys., 85, 1021
  • Tsuboi et al. (1986) Tsuboi, M., Inoue, M., Handa, T., Tabara, H., Kato, T., Sofue, Y., & Kaifu, N.  1986, AJ, 92, 818
  • Tsuboi, Handa & Ukita (1999) Tsuboi, M., Handa, T., & Ukita, N.  1999, AJ, 120, 1
  • Tsuboi et al. (2011) Tsuboi, M., Tadaki, K-I., Miyazaki, A., & Handa, T.  2011, PASJ, 63, 763
  • Tsuboi et al. (2015) Tsuboi, M., Miyazaki, A., & Uehara, K.  2015, PASJ, 67, id. 109
  • Tsuboi et al. (2019) Tsuboi, M., Kitamura, Y., Uehara, K., Miyazaki, A., Miyawaki, R., Tsutsumi, T., & Miyoshi, M.  2019, PASJ, 71, id.128
  • Uehara et al. (2019) Uehara, K., Tsuboi, M., Kitamura, Y., Miyawaki, R., & Miyazaki, A.,  2019, ApJ, 872, id. 121
  • Uehara et al. (2020) Uehara, K., Tsuboi, M., Kitamura, Y., Miyawaki, R., & Miyazaki, A.,  2020, in preparation.
  • Yusef-Zadeh, Morris, & Chance (1984) Yusef-Zadeh, F., Morris, M., & Chance, D.  1984, Nature, 310, 557
  • Yusef-Zadeh et al. (2013) Yusef-Zadeh, F., Cotton, W., Viti, S., Wardle, M., & Royster, M.,  2013, ApJ, 764, id. L19