Abstract
We present a new analysis of the interstellar protons toward the TeV γ-ray SNR RX J0852.0−4622 (G266.2−1.2, Vela Jr.). We used the NANTEN2 12CO(J = 1–0) and Australia Telescope Compact Array and Parkes H i data sets in order to derive the molecular and atomic gas associated with the TeV γ-ray shell of the SNR. We find that atomic gas over a velocity range from VLSR = −4 to 50 km s−1 or 60 km s−1 is associated with the entire SNR, while molecular gas is associated with a limited portion of the SNR. The large velocity dispersion of the H i is ascribed to the expanding motion of a few H i shells overlapping toward the SNR but is not due to the Galactic rotation. The total masses of the associated H i and molecular gases are estimated to be M⊙ and ∼103 M⊙, respectively. A comparison with the High Energy Stereoscopic System TeV γ-rays indicates that the interstellar protons have an average density around 100 cm−3 and shows a good spatial correspondence with the TeV γ-rays. The total cosmic-ray proton energy is estimated to be ∼1048 erg for the hadronic γ-ray production, which may still be an underestimate by a factor of a few due to a small filling factor of the SNR volume by the interstellar protons. This result presents a third case, after RX J1713.7−3946 and HESS J1731−347, of the good spatial correspondence between the TeV γ-rays and the interstellar protons, lending further support for a hadronic component in the γ-rays from young TeV γ-ray SNRs.
Export citation and abstract BibTeX RIS
1. Introduction
The origin of the cosmic rays (CRs) is one of the most fundamental issues in modern astrophysics since the discovery of the CRs by Victor Franz Hess in 1912. It is particularly important to understand the origin of the CR protons, the dominant constituent of the CRs. There have been a number of studies to address the acceleration sites of CR protons in the Galaxy (e.g., Hess & Steinmaurer 1935). It is widely accepted that CR particles are accelerated up to the energies below 1015 eV via the diffusive shock acceleration (DSA), which takes place for instance between upstream and downstream of the high velocity shock front of the SNR (e.g., Bell 1978; Blandford & Ostriker 1978). It is crucial to verify observational signatures for the hadronic γ-rays of the CR proton origin.
The recent advent of high-resolution γ-ray observations have allowed us to image some 20 SNRs, offering a new opportunity to identify the hadronic γ-rays. Among them, four SNRs, RX J1713.7−3946, RX J0852.0−4622, RCW 86, and HESS J1731−347, show shell-type TeV γ-rays morphology as revealed by High Energy Stereoscopic System (H.E.S.S.; Aharonian et al. 2004, 2005, 2006, 2007a, 2007b, 2009; H.E.S.S. Collaboration et al. 2011, 2016a, 2016b, 2016c). Potential hadronic γ-rays are produced by the proton–proton collisions followed by neutral pion decay, while any leptonic γ-rays come from CR electrons via the inverse-Compton effect and possibly Bremsstrahlung.
Previous interpretations of the γ-ray observations usually used mainly the γ-ray spectra in order to discern the above two processes, while some attempts were made to consider explicitly the target protons in the interstellar medium (ISM; e.g., Aharonian et al. 2006, 2007b). A conventional assumption is that the target protons are uniform, having densities of 1 cm−3 (e.g., Aharonian et al. 2006). Such a low density is preferred because the DSA works efficiently at such low density. It is, however, becoming recognized that the γ-ray spectrum alone may not be sufficient to settle the origin of the γ-rays because the penetration of the CR protons into the target ISM may significantly depend on the density of the ISM (Gabici et al. 2007); the higher energy CR protons can penetrate deeply into the surrounding molecular cloud cores, whereas the lower energy CR protons cannot, making the γ-ray spectrum significantly harder than the fully interacting case (Inoue et al. 2012, hereafter I12: see also references therein). I12 therefore suggest that the γ-ray spectra are not usable to discern the γ-ray production mechanisms, but that the spatial correspondence between the γ-rays and ISM distribution is a key element in testing for a hadronic γ-ray component. The effect of the clumpy ISM was also discussed by Zirakashvili & Aharonian (2010) and Gabici & Aharonian (2014). It is often noted that the acceleration via DSA must happen in a low-density space, which has too low a density for the hadronic origin to be effective (e.g., Ellison et al. 2010). This is, however, not a difficulty if one takes into account the highly inhomogeneous ISM distribution as is commonly the case in a stellar wind evacuated cavity with a dense surrounding ISM shell-like in RX J1713.7−3946 (Sano et al. 2010).
Several of the previous studies explicitly analyzed the distribution of the candidate target ISM protons observed by CO and compared them with the γ-ray distribution in the most typical TeV γ-ray SNRs including RX J1713.7−3946 and RX J0852.0−4622. Aharonian et al. (2006) compared the TeV γ-ray distribution with the millimeter-wave rotational transition of CO, the tracer of H2, obtained with the NANTEN 4 m telescope. These authors found some similarity between CO and γ-rays in RX J1713.7−3946, whereas part of the γ-ray shell was completely missed in CO. The authors thus did not reach a firm conclusion on the target ISM protons in the hadronic scenario. In the case of RX J0852.0−4622, Aharonian et al. (2007b) made a similar comparison in a velocity range of 0–20 km s−1 but found little sign of the target protons in CO. So, the ISM associated with RX J0852.0−4622 remained ambiguous.
Fukui et al. (2012; hereafter F12) carried out a detailed analysis of the ISM protons toward the RX J1713.7−3946 by employing both the molecular and atomic protons, and have shown for the first time that these ISM protons have a good spatial correspondence with the TeV γ-rays. This new study showed that the atomic protons are equally important as molecular protons as the target in the hadronic process, which was not previously taken into account. The study by F12 provides a necessary condition for the hadronic process, lending support for the hadronic scenario, although this correspondence alone does not exclude the leptonic process. F12 and I12 further considered the other relevant aspects including magneto-hydro-dynamical numerical simulations of the SN shocks and X-ray observations (see also Sano et al. 2010, 2013, 2015) and argued that the TeV γ-rays in RX J1713.7−3946 is likely emitted via the hadronic process.
Another young SNR RX J0852.0−4622 was discovered by Aschenbach (1998), and RX J0852.0−4622 shows a hard X-ray spectrum toward part of the more extended Vela SNR in the ROSAT All-Sky Survey image. TeV γ-rays were detected and imaged toward RX J0852.0−4622 by H.E.S.S. RX J0852.0−4622 shows similar properties with RX J1713.7−3946; they are both young with ages of 2400–5100 years for RX J0852.0−4622 (Allen et al. 2015): ∼1600 year for RX J1713.7−3946 (e.g., Fukui et al. 2003; Moriguchi et al. 2005) having synchrotron X-ray emission without thermal features and share shell-like X-ray/TeV γ-ray morphology. RX J0852.0−4622 has an apparently large diameter of about 2°. The size will allow us to test spatial correspondence between the γ-rays and the ISM at 012 angular resolution in FWHM of H.E.S.S. We have two more shell-like TeV γ-ray SNRs RCW 86 and HESS J1731−347 (Aharonian et al. 2009; H.E.S.S. Collaboration et al. 2011). In HESS J1731−347, Fukuda et al. (2014) carried out a comparative analysis of the ISM and γ-rays and have shown that the ISM protons show that their spatial distributions are similar with each other, being consistent with a dominant hadronic component of g-rays plus minor contribution of a leptonic component. In RCW 86, a similar comparative analysis is being carried out by H. Sano et al. (2017, in preparation).
The distance of RX J0852.0−4622 was not well determined in previous works (e.g., Slane et al. 2001a, 2001b; Iyudin et al. 2007; Katsuda et al. 2008; Pannuti et al. 2010; Allen et al. 2015; Maxted et al. 2017). A possible distance of RX J0852.0−4622 was 250 ± 30 pc similar to the Vela SNR (Cha et al. 1999), while another possible distance was larger than that of the Vela SNR. Slane et al. (2001a, 2001b) argued that RX J0852.0−4622 is physically associated with the giant molecular cloud, the Vela Molecular Ridge (VMR, May et al. 1988; Yamaguchi et al. 1999b), and the distance of the VMR was estimated to be 700–200 pc (Liseau et al. 1992). It is, however, not established if the VMR is physically associated with RX J0852.0−4622 (see, e.g., Pannuti et al. 2010). Toward the northwestern rim of RX J0852.0−4622, two observations were made with the XMM-Newton and an expansion velocity was derived and an age of (1.7–4.3) × 103 years was estimated by Katsuda et al. (2008), assuming the free expansion phase. Recently, Allen et al. (2015) improved the expansion measurement by using two Chandra data sets separated by 4.5 years toward the northwestern rim of RX J0852.0−4622. They derived an age of (2.4–5.1) × 103 year old and the range of distance from 700 to 800 pc, which are roughly consistent with the previous study by Katsuda et al. (2008). We will therefore adopt the distance ∼750 pc to RX J0852.0−4622 in the present paper. Analyses of a central X-ray source in this SNR using data from multiple observatories suggest that the progenitor of the SNR is a high-mass star, which led to a core-collapse type SNe (Aschenbach 1998; Mereghetti 2001; Pavlov et al. 2001; Slane et al. 2001a). This conclusion is also consistent with the estimate of the current SNR shell expansion speed (Chen & Gehrels 1999).
In this paper, we present the results of an analysis of the ISM protons toward RX J0852.0−4622 by using both the CO and H i data. Section 2 gives the observations of CO and H i. Section 3 describes the high-energy γ-ray and X-ray data. Section 4 gives the results of the CO and H i analyses and Section 5 provides a discussion. We conclude the paper in Section 6.
2. Observations
2.1. CO Observations
Observations of the 12CO(J = 1–0) transition were carried out by the NANTEN 4 m telescope of the Nagoya University at the Las Campanas Observatory (2400 m above the sea level) in Chile in 1999 May–July (Moriguchi et al. 2001). The half-power beam width at the frequency of the 12CO(J = 1–0), 115.290 GHz, was . The observations were made by position switching mode and grid spacing was . SIS (superconductor-insulator-superconductor) mixer receiver provided Tsys of ∼250 K in the single side band including the atmosphere in the direction of the zenith. The spectrometer was an AOS (acousto-optical spectrometer) with 40 MHz bandwidth and 40 kHz resolution, providing 100 km s−1 velocity coverage and a velocity resolution of 0.1 km s−1. The velocity always refers to that of the local standard of rest. The rms noise per channel is ∼0.5 K.
2.2. H i Observations
We used the high resolution 21 cm H i data obtained with Australia Telescope Compact Array (ATCA) in Narrabri, New South Wales in Australia consisting of six 22 m dishes. The region of RX J0852.0−4622 in l = 266° was observed as part of the Southern Galactic Plane Survey (SGPS; McClure-Griffiths et al. 2005) and was combined with the H i data taken with the 64 m Parkes telescope. The Parkes H i data were taken for b = −10° to +10°, while the ATCA covers b = −15 to . In order to supplement the southwestern half of the SNR new H i observations were conducted during 24 hr on 2011 February 26–27 and March 29–30, with the ATCA in the EW352 and EW367 configurations (Project ID: C2449, PI: Y. Fukui). We employed the mosaicking technique, with 43 pointings arranged in a hexagonal grid at the Nyquist separation of . The absolute flux density was scaled by observing PKS B1934−638, which was used as the primary bandpass and amplitude calibrator. We also periodically observed PKS 0823−500 for phase and gain calibration. The MIRIAD software package was used for the data reduction (Sault et al. 1995). We combined the ATCA data set with single dish observations taken with the Parkes 64 m telescope and the SGPS data set. The final beam size of H i is 245'' × 130'' with a position angle of 117°. The typical rms noise level was 1.4 K per channel for a velocity resolution of 0.82 km s−1.
3. Distributions of the SNR
Figure 1(a) shows the TeV γ-ray and X-ray distributions of RX J0852.0−4622. These high energy features are thin and shell-like and are enhanced toward the northern half.
The γ-ray distribution (Figure 1 of Aharonian et al. 2007b) is obtained in the energy above 0.3 TeV with H.E.S.S. installed in Namibia, which utilized the first four 12 m diameter Cherenkov telescopes of the current five-telescope H.E.S.S. array. The H.E.S.S. image has a point-spread function of 006 at the 68% radius (hereafter termed "r68") or the half-power full width of 0144 for 33 hr observations.11 The shell is well resolved for the 2° diameter, making RX J0852.0−4622 suitable for testing the spatial correspondence between the ISM and the γ-rays. The γ-rays peaked toward (l, b) ∼ (26697, −100) correspond to a pulsar wind nebula PSR J0855−4644 at a distance of below 900 pc and is not related to the RX J0852.0−4622 (Acero et al. 2013). The total γ-ray flux in the energy range from 0.3 to 30 TeV is estimated to be () × 10−12 erg cm−2 s−1 (H.E.S.S. Collaboration et al. 2016b).
ROSAT observations (Aschenbach 1998; Aschenbach et al. 1999), ASCA observations (Tsunemi et al. 2000; Slane et al. 2001a, 2001b), and Suzaku observations (Takeda et al. 2016) show that the object is shell-like, with luminous northern and southern rims and less luminous northeastern and southwestern rims. The X-rays are superposed on the γ-ray image in Figure 1(a). The most prominent X-ray peak is seen toward the northwestern shell, (l, b) = (2654, −12) (Tsunemi et al. 2000; Slane et al. 2001a, 2001b). The X-rays are nonthermal with a photon index around 2.7 and the absorbing column density of (2–4) × 1021 cm−2 (e.g., Tsunemi et al. 2000; Slane et al. 2001a, 2001b; Bamba et al. 2005; Iyudin et al. 2005). ASCA and Chandra observations of the SNR show that the nonthermal X-rays are dominant with an absorption column density (1–4) × 1021 cm−2 toward the X-ray peak. A spatially resolved spectroscopic X-ray study revealed details of the fine structure in the luminous northwestern rim complex of G266.2−1.2 by observations made with Chandra (Bamba et al. 2005; Iyudin et al. 2005; Pannuti et al. 2010; Allen et al. 2015).
4. ISM Distributions
4.1. CO and H i Distributions
CO is a tracer of molecular gas with density of a few times 100 cm−3 or higher and the H i traces the atomic gas with lower density less than several 100 cm−3. The density range between them may be probed by the cool H i gas as seen in self-absorption if the background H i is bright enough. Such H i self-absorption is in fact observed in RX J1713.7−3946 and is identified as part of the target ISM protons in the hadronic γ-ray production (F12).
Figure 1(b) shows a three-color image of RX J0852.0−4622, consisting of 12CO(J = 1–0) in red, H i in green, and X-ray (2–5.7 keV) in blue. H i emission delineates the outer boundary of the X-ray shell, except for the southeast and northeast. In the western region, the X-ray shell is well correlated with the CO clumps.
Figures 2 and 3 show the velocity channel distributions of 12CO(J = 1–0) and H i every 5 km s−1 in the range from −10 to 80 km s−1, where the H.E.S.S. TeV γ-ray contours are superposed. The giant molecular cloud at 0–15 km s−1 is the VMR (Murphy & May 1991; Yamaguchi et al. 1999a, 1999b; Moriguchi et al. 2001). The other CO features are all small and clumpy. We list the observational parameters of the small CO clouds in Table 1; the position of peak intensity within the cloud, velocity, line intensity, size, and mass.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageTable 1. Properties of CO Clouds toward the SNR RX J0852.0−4622
Name | l | b | Size | Mass | |||
---|---|---|---|---|---|---|---|
(degree) | (degree) | (K) | (km s−1) | (km s−1) | (pc) | (104 ) | |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) |
CO0S | 267.07 | −1.63 | 2.86 | 3.47 | 3.48 | 3.2 | 220 |
CO20E | 266.87 | −0.67 | 11.46 | 19.80 | 2.06 | 6.8 | 1330 |
CO25W | 266.47 | −2.04 | 7.48 | 24.00 | 1.72 | 6.2 | 650 |
CO25C | 266.13 | −1.00 | 2.04 | 24.45 | 2.68 | 1.9 | 40 |
CO30E | 266.67 | −0.87 | 6.28 | 31.09 | 2.78 | 3.0 | 180 |
CO45NW | 265.33 | −1.10 | 4.42 | 44.68 | 4.07 | 4.0 | 330 |
CO60NW | 265.53 | −1.47 | 3.31 | 63.98 | 2.25 | 6.3 | 410 |
Note. Column (1): name of CO cloud. Columns (2–7): physical properties of the CO cloud obtained by a single Gaussian fitting. Columns (2)–(3): position of the peak intensity in the Galactic coordinate. Column (4): radiation temperature. Column (5): center velocity. Column (6): full width at half maximum (FWHM) of line width. Column (7): size of CO cloud defined as , where A is the area of cloud surface surrounded by the CO contours in Figure 5. Column (8): mass of CO cloud is defined as mH μ Ω , where mH is the mass of the atomic hydrogen, μ is the mean molecular weight, D is the distance to RX J0852.0−4622, Ω is the angular size in pixel, and N() is column density of molecular hydrogen for each pixel. We used by taking into account the helium abundance of 20% relative to the molecular hydrogen in mass, and N() = 2.0 × 1020 [W(12CO) (K km s−1)] () (Bertsch et al. 1993).
Download table as: ASCIITypeset image
4.2. Associated Clouds
We shall first identify seven candidate CO clouds named CO 0 S, CO 20 E, CO 25 W, CO 25 C, CO 30 E, CO 45 NW, and CO 60 NW, where their typical velocity is indicated by Figure 5. In addition, we show five H i features, H i 0 S, H i 25 W, H i 30 E, H i 45 N, and H i 60 NW in Figure 6, which are candidates for the H i counterparts of the CO. By combining these CO with the H i and the shell distribution, we select the plausible candidates for the clouds interacting with the SNR.
Figure 4(a) shows an overlay between CO and H i. The CO corresponds well the southwestern rim of the shell. This shows a clear association of CO cloud CO 25 W with the SNR. H i 25 W is also associated with this CO showing a shift in position from the center to the northwest in Figure 4(a). The shift is consistent with part of an expanding shell on the near side at an expanding velocity around 5–10 km s−1. We also note that H i 25 W is likely located toward the tangential edge of the shell, suggesting that the central velocity of the ISM associated with the SNR is around 25 km s−1.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageFigures 4(b) and (c) show two additional cases of the association. Figure 4(b) shows CO 30 E with a H i tail extending outward from the center, H i 30 E. The H i is V-shaped pointing toward the center of the SNR. CO 30 E is located at the tip of the H i and is elongated in the radial direction. We suggest that the distribution presents a blown-off cloud overtaken by the stellar wind by the SN progenitor, where the dense head of CO survived against the stellar wind with a H i tail. A similar morphology is seen toward CO 25 C in Figure 4(c). The CO is also elongated toward the center and the H i, part of H i 45 N, shows an elongation in the same direction with the CO at the inner H i tip (Figure 4(c)). We discuss more quantitative details of the interaction in the next section. The three CO features, CO 25 W and CO 25 C, in a velocity range from 20 to 30 km s−1 suggest the physical association of the CO and H i with RX J0852.0−4622. We also note that the small double feature CO 0 S is seen toward the southeast of the γ-ray feature. These features and the H i 0 S (Figure 6) show correspondence with the γ-rays. To summarize, the CO in a velocity range from −3 to 30 km s−1 with a central velocity 25 km s−1 shows signs of the association with the SNR.
Based on these candidates for the association, we extended a search for CO and H i in a velocity from −5 to 65 km s−1 in Figures 2 and 3, and have summarized the possible candidate features in Figures 5 and 6. These candidate CO and H i features are CO 20 E, CO 45 NW, and CO 60 NW, and H i 45 N and H i 60 NW. It is notable that CO 25 W is found toward the γ-ray shell (panel 25–30 km s−1).
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution image4.3. Distance of the ISM: H i Supershells
The central velocity of the ISM associated with RX J0852.0−4622 is around 25 km s−1 with a velocity span of ∼50 km s−1. 25 km s−1 corresponds to a kinematic distance of 4.3 kpc if the galactic rotation model is adopted (Brand & Blitz 1993), significantly different from the adopted distance value of 750 pc (Katsuda et al. 2008; Allen et al. 2015). We shall look at the H i distribution around the SNR in order to clarify this discrepancy.
Figures 7(a), (b), and (c) show the H i distributions in the velocities of 66 km s−1, 45 km s−1, and 35 km s−1, respectively. By a conventional Galactic rotation model, these velocities correspond to distances of 8.7, 6.4, and 5.4 kpc (Brand & Blitz 1993). We argue here that the Galactic rotation is not the dominant cause of the H i velocity but that the expanding motion of several H i supershells is mainly responsible for the line-of-sight velocity. In fact, Suad et al. (2014) and Arnal & Corti (2007) identified the H i supershells, GS 268−01+066 and GS 263−02+45 (see Figures 7(a) and (b)) toward the SNR. In addition, we newly identified the H i supershell "GS 265−04+35" as shown in Figure 7(c). The supershell has also an expanding motion and shows the front and rear wall in the H i spectrum (see more details in the Appendix). Figure 7(d) shows a schematic view of the supershell boundaries toward the SNR RX J0852.0−4622. The TeV γ-ray contours overlap three supershells, indicating that the expanding motion of the supershells dominates the H i velocity field around the SNR. Their detailed physical parameters are shown in Table 2.
Download figure:
Standard image High-resolution imageTable 2. Observed Properties of H i Supershells toward RX J0852.0−4622
ID | Name | l | b | vc | r | a | b | ϕ | vexp | Reference |
---|---|---|---|---|---|---|---|---|---|---|
(degree) | (degree) | (km s−1) | (degree) | (degree) | (degree) | (degree) | (km s−1) | |||
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | (11) |
1 | GS 268−01+066 | 267.8 | −1.1 | 66 | ⋯ | 2.2 ± 0.2 | 1.7 ± 0.2 | −36.9 | 10.3 | (1) |
2 | GS 263−02+45 | 262.6 | −1.8 | 45 | 3.6 ± 0.4 | ⋯ | ⋯ | ⋯ | 14 | (2) |
3 | GS 265−04+35 | 264.9 | −3.6 | 35 | 2.47 ± 0.05 | ⋯ | ⋯ | ⋯ | 9 | This work |
Note. Column (1): supershell ID. Column (2): supershell name. Columns (3–4): central position of the the supershell. Column (5): central radial velocity of the supershell. Column (6): radius of the fitted circle. Columns (7–8): major and minor semi-axis of the fitted ellipse. Column (9): major axis inclination ϕ relative to the Galactic longitude, measured counterclockwise from the Galactic plane. Column (10): velocity extension of the supershell. Column (11): (1) Suad et al. (2014); (2) Arnal & Corti (2007).
Download table as: ASCIITypeset image
The H i position–velocity diagram on a large scale is shown in Figure 8. The H i gas toward the SNR is not ordered in the spiral arm at 60–80 km s−1 and the H i has several velocity features at 0 km s−1, 20–30 km s−1, and 40–50 km s−1 toward RX J0852.0−4622 in l = 263°–268°, all of which show kinematic properties consistent with part of an expanding shell (Figure 7). We therefore suggest that the apparent velocity shifts of the H i in the region are due to the expansion of supershells, which are driven by several stellar clusters located in the center of each shell. The known molecular supershells, including the Carina flare supershell, show an expanding velocity of 10 km s−1 with a ∼10 Myr age (Fukui et al. 1999; Matsunaga et al. 2001; Dawson et al. 2008a, 2008b, 2011a, 2011b, 2015). We generally find no direct observational hints for any stellar clusters although it is not unusual that we see no clear signs of clusters, which should become faint, in 10 Myr due to the evolution of high mass stars in the cluster.
Download figure:
Standard image High-resolution image4.4. Total ISM Protons
We here derive the total ISM proton density. The molecular column density is calculated by canonical conversion factors and the error corresponds to 3σ noise fluctuations. The total intensity of 12CO(J = 1–0) can be converted into the molecular column density N(H2) (cm−2) by the following relation:
where the factor N(H2) (cm−2)/W(12CO) (K km s−1) is adopted as (cm−2/K km s−1) (Bertsch et al. 1993). Then, the total proton density is given as Np(H2) = 2N(H2).
Usually, the atomic proton column density is estimated by assuming that the 21 cm H i line is optically thin. If this approximation is valid, the H i column density Np(H i) (cm−2) is estimated as follows (Dickey & Lockman 1990):
where Tb (K) and V (km s−1) are the 21 cm intensity and the peak velocity. Fukui et al. (2014, 2015) made a new analysis by using the submillimeter dust optical depth derived by the Planck satellite (Planck Collaboration et al. 2014) and concluded that the H i emission is generally optically thick with an optical depth of around 1 in the local ISM within 200 pc of the Sun at the Galactic latitude higher than 15°. This optical depth correction increases the H i density by a factor of ∼2 as compared with the optically thin case. We adopt this correction by using the relationship between the H i integrated intensity and the 353 GHz dust optical depth given in Figure 9 (see also Planck Collaboration et al. 2014). Since RX J0852.0−4622 is close to the Galactic plane, the 353 GHz dust optical depth is not available for a single velocity component as in the local space. Instead, we adopt the empirical relationship between and in Figure 9 in order to estimate the from . The ratio of the and is determined as follows in the optically thin limit (Fukui et al. 2015, Equation (3)):
hence,
where the nonlinear dust property, derived by Roy et al. (2013) and Okamoto et al. (2017) is assumed, and this nonlinearity does not alter significantly in the present range.
Download figure:
Standard image High-resolution imageThe H i velocity range is taken to be 20–50 km s−1 and −4–1 km s−1, where VMR is dominant in the velocity range from −5 to 15 km s−1. The CO 0 S, CO 25 W, and CO 25 W locations are also taken into account as the molecular components. The total ISM proton column density is given by the sum as follows:
Figures 10(a), (b), and (c) present the total molecular protons, atomic protons and the sum of the molecular and atomic protons, respectively. The resolution is adjusted to the major axis of the beam size ().
Download figure:
Standard image High-resolution imageOur assumed distance around 1 kpc for most of the associated ISM features is independently confirmed by the visual extinction AV seen in the southern half of the SNR, where the foreground contamination by the VMR is not significant. Figure 11 shows the Av distribution from the Digitized Sky Survey I (DSS I; Dobashi et al. 2005) where AV is estimated by star counting of the 2MASS data. The direction (l, b) = (266°, −1°) is where the number of stars is small as compared with the Galactic center and AV is not very accurately estimated in general. Nevertheless, we see a clear indication of enhanced AV features in Figure 11, where AV is typically 0.7–1.7 mag as is clearly seen toward the southern part of the γ-ray shell, and the CO 0 S, H i 0 S, and CO 25 W clouds. The peaks have a total proton column density of (5–10) × 1021 cm−2 in Figure 10(c), corresponding to AV = 3–5 mag (AV = Np(H)/1.7 × 1021 cm−2 mag). The fact that the extinction is visible toward the SNR indicates that the cloud is relatively close to the Sun, and the image in Figure 11 is consistent with a distance around 1 kpc.
4.5. TeV γ-rays
The histogram in Figure 12 indicates the average TeV γ-ray count taken every 01 as a function of radius from the shell center (l, b) = (26628, −124) obtained by H.E.S.S. (Aharonian et al. 2007b). The error is conservatively estimated as the (oversampling-corrected total smoothed count)0.5. We shall here analyze the γ-ray distribution by using a spherical shell model. The distribution is fit by a spherically symmetric shell, following the method of F12 and a Gaussian intensity distribution in radius of the γ-ray counts as follows:
where A is a factor for normalization, the radius r0 at the peak, and the width of the shell σ. The green line in Figure 12 shows the results. , . They correspond to ∼12 pc and ∼2.4 pc, at 750 pc.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageThe azimuthal γ-ray distributions are obtained by taking the two circles centered at (l, b) = (26628, −124), where the radii of these circles are determined to be 064 and 091, respectively, at the one-third count level of the peak of the Gaussian shell shown in Figure 12. The projected distribution is shown by the orange color. The range of the radius is the same as that adopted by Aharonian et al. (2007b). The fitting seems adequate in a range from 00 to 12.
5. Discussion
Spatial correspondence between γ-rays and ISM protons is a key to discerning the γ-ray production mechanism (I12, F12). Figure 10 shows a comparison between γ-rays and ISM in the azimuthal distribution (the angle is measured clockwise with 0° in the southeast). The TeV γ-ray counts were normalized.
The correspondence between the two is remarkably good, while some small deviations between them are seen at angles of −60° to 0° and 120° to 150°, and at the inside of the shell. This offers a third case after RX J1713.7−3946 (F12), which supports a hadronic origin of the γ-rays via spatial correspondence. The total H i and H2 masses involved in the shell are estimated to be and within the radius of ∼15 pc, respectively.
In Figure 10(a), there is a point that shows significant deviation at an angle of 0° to 30° and this corresponds to the position of the PWN unrelated to the SNR shell, because the pulsar characteristic age is 140 kyr and the large velocity ∼3000 km s−1 needed if the pulsar is a remnant of the SNR (e.g., Acero et al. 2013; H.E.S.S. Collaboration et al. 2016b). The positions between 30° and 150° show a trend that the ISM proton column density underestimates the γ-rays by less than 10%. These positions show no CO emission and the H i is only responsible for the ISM protons.
In Figure 10(c), the shell seems more intense in the γ-rays and may suggest that CR is enhanced as suggested by the significantly enhanced X-rays. We estimate the total energy of CR protons above 1 GeV by using the equation below (H.E.S.S. Collaboration et al. 2016b):
where d is the distance, 750 pc, and n the proton density, ∼100 cm−3, by adopting the shell radius ∼15 pc and the thickness ∼9 pc. is estimated to be erg, corresponding to 0.1% of the total kinetic energy of an SNe ∼1051 erg. The H i distribution seems to be fairly uneven in space and velocity as suggested by the nonuniform distribution in Figure 6. The coupling between CR protons and the target protons may not be complete. According to Inoue et al. (2012), the effective mean target density for CR protons ntg can be written as ntg ≃ n f, where f is the volume filling factor of the interstellar protons. The value is therefore to be regarded as a lower limit.
In RX J1713.7−3946, the total molecular mass and atomic mass are ∼104 for each and the total energy of the CR protons of 1048 erg (F12). In RX J0852.0−4622, the total atomic mass is , while the total molecular mass is . Given the radii of the two SNRs, the average density over the whole volume of RX J1713.7−3946 is similar to that in RX J0852.0−4622. We find that these two young TeV γ-ray SNR have total CR proton energy at an order of 1048 erg. The energy is significantly smaller than those discussed before (e.g., Aharonian et al. 2006) and suggests that the fraction of the explosion energy converted into CRs appears fairly low 0.1% in such a young stage coupled with the effect of the volume filling factor. The efficiency may possibly grow in time and this can be tested by exploring middle-aged SNRs with an age of 104 years. In fact, for the middle-aged SNRs W44, W28, and IC443, the total CR proton energy is ∼1049 erg (Giuliani et al. 2010, 2011; Yoshiike et al. 2013, 2017).
Nonthermal X-rays indicate particle acceleration in the SN blast waves. Slane et al. (2001a, 2001b) analyzed the X-rays with ASCA observations. They showed that the nonthermal X-rays are dominant similarly to RX J1713.7−3946. The absorption column density in front of the X-rays is cm−2 toward the western rim, which is consistent with the present ISM proton distribution in Figure 10(c). ISM inhomogeneity is suggested by the CO clumps in Figure 5 (CO 30 E and CO 25 C) overtaken by the blast waves. The situation is similar to the inside of RX J1713.7−3946. DSA works for CR acceleration in the low-density cavity and the CR protons can collide with the ISM protons in the SNR. The diffusion length is estimated to be 1 pc for the magnetic field of 10 μG and CR proton energy of 100 TeV, large enough for the CR protons to interact with the CO and H i.
We note that there is a correlation between the γ-rays, X-rays, and the ISM. All of these components are enhanced on the western rim of the SNR but the eastern rim is quite weak. This distribution is interpreted in terms of the enhanced turbulence in the western rim such as suggested by magneto-hydro-dynamical numerical simulations by I12. These authors showed that the shock interacting with dense clumps create turbulence and that the turbulence amplifies the magnetic field up to 1 mG. Thus enhanced magnetic field can intensify the synchrotron X-rays and efficiently accelerate CRs, leading to the increase of the γ- and X-rays.
RX J0852.0−4622 is suggested to be a core-collapse SNR (Aschenbach 1998). The parent cloud of the SN progenitor is likely the one at 25 km s−1 elongated along the plane by 100 pc to the west. It is the most massive H i complex in the region. The total H i mass of the cloud is roughly estimated to be 105 and the H2 mass 104 . The H i cloud is found to be located in the intersection among at least three H i supershells, which are expanding at 10–15 km s−1. Collisions among the shells may be responsible for the formation of the parent cloud in this relatively diffuse environment outside the solar circle (e.g., Dawson et al. 2011a, 2011b).
6. Summary
We have carried out a combined analysis of CO and H i toward the youngest TeV γ-ray SNR RX J0852.0−4622, following the first good spatial correlation found in RX J1713.7−3946. The main conclusions of the present study are summarized as follows.
- 1.The ISM in a velocity range from −4 to 50 km s−1 is likely associated with the SNR. The ISM is dominated by the H i gas and the mass of the molecular gas probed by CO corresponds to about 4% of the H i gas. This association is supported by the morphological signs of the interaction between CO/H i and the SNR, such as cometary-tailed and shell-like shapes, in 20–30 km s−1. The visual extinction in the south of the SNR corresponding to the ISM protons lends more support for the association.
- 2.The total ISM protons show a good spatial correspondence with the TeV γ-rays in the azimuthal distribution. This provides a third case of such correspondence next to the SNR RX J1713.7−3946 and HESS J1713−347, a necessary condition for a hadronic component of the γ-ray emission. The total CR proton energy is estimated to be 1048 erg from the γ-rays measured by H.E.S.S. and the average ISM proton density 100 cm−3 derived from the present study. The CR energy is a similar value obtained for RX J1713.7−3946 and is 0.1% of the total kinetic energy of an SNe.
- 3.The large velocity splitting of the H i gas is likely due to the expansion of a few supershells driven by star clusters but is not due to the Galactic rotation. The velocity range from 0 to 60 km s−1 corresponds to that of the inter-arm features between the local arm and the Perseus arm. We find a possible parent H i/CO cloud where the high mass stellar projenitor of the SNR was formed. This cloud located toward the interface of a few supershells and has a mass of 105 in H i.
The NANTEN project is based on a mutual agreement between Nagoya University and the Carnegie Institution of Washington (CIW). We greatly appreciate the hospitality of all the staff members of the Las Campanas Observatory of CIW. We are thankful to many Japanese public donors and companies who contributed to the realization of the project. This study was financially supported by Grants-in-Aid for Scientific Research (KAKENHI) of the Japanese society for the Promotion of Science (JSPS, grant Nos. 12J10082, 24224005, 25287035, and 15H05694). This work also was supported by "Building of Consortia for the Development of Human Resources in Science and Technology" of Ministry of Education, Culture, Sports, Science and Technology (MEXT, grant No. 01-M1-0305).
Software: MIRIAD (Sault et al. 1995).
Appendix: Identification of the Supershell GS 265−04+35
In order to identify the supershell GS 265−04+35, we used the following three steps.
- 1.Searching for a central velocity ofthe H i supershell—Investigating the H i profile toward the supershell and searching for the velocity with minimum H i intensity. We defined the velocity as the central velocity of an expanding motion (see McClure-Griffiths et al. 2002). The front and rear walls were also identified (Figure 13). The expanding velocity was estimated to be ∼9 km s−1 as a difference between the central and front wall velocities (see also Table 2).
- 2.Estimating a geometric center of the H i supershell—In the H i map of the central velocity, we estimated a radial distance of the supershell as defined at peak H i intensity every 15° for an assumed geometric center. We circulated only the radial distance with the peak H i intensity of 30 K or higher. We minimized the dispersion of radial distances in azimuth angles, which allows the geometric center to be (l, b) = (2649, −36) (see also Table 2).
- 3.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageFinally, we confirmed that GS 265−04+35 is the expanding H i supershell.
Footnotes
- 11
A higher statistics γ-ray image was presented by the H.E.S.S. Collaboration et al. (2016b) with an angular resolution of r68 = 008, poorer than 006 reported for the morphological analysis in Aharonian et al. (2007b). Since there is no significant difference between the two images except for the angular resolution and photon statistics, we therefore use the previous image in Aharonian et al. (2007b) to improve the morphological study in this paper.