arXiv:astro-ph/0108031v1 2 Aug 2001 eli h vrapdeeg region. agre energy spectra overlapped two TeV the The in 1 conference. well around this spectrum in presented preliminary is the region and tele- CANGAROO-II 2000, 10m in the scope using position same the from h bandfil teghadmxmmeeto nryare energy electron 4 maximum and en- strength field maximum accura obtained considerable the the the and with electrons shock accelerated of of the ergy strength around absolute field the magnetic estimate the to Background us Microwave enables Cosmic This 2.7k (CMB). with process (IC) Com Inverse ton the using calculated model the matches spectrum N06wsotie rmtedt bevdi 96and of 1996 rim in observed north-east CANGAROO 3.8m data the the the using from 1997 from TeV obtained 20 was and SN1006 TeV 1.5 between Abstract. naee l 96,adtesbeun eeto fTeV of detection u.ac.jp) subsequent the and to: (Will-Correspondence 1996), ROSAT and 1995) al. al. et et ingale (Koyama (G327.6 SN1006 ASCA SNR X- by Ia Type non-thermal +14.6) of intense rims the the from is emission observations ray such of set First to up accelerated ∼ are which particles (SNR), Remnants non-thermal Supernova of the existence dra- the rays gamma revealed TeV and matically rays X both in observations Recent Introduction 1 Tanimori T. NE the around 1006 SN of Rim spectrum gamma-ray TeV the of Study 1–4 2001: ICRC of Proceedings 11 Sakurazawa Matsubara 1 12 2 3 13 4 5 14 6 7 16 8 9 10 aut fMngmn nomto,YmnsiGki Unive Gakuin Yamanashi Information, Management of Faculty aut fSine brk nvriy io brk 310- Ibaraki Mito, University, Ibaraki Science, of Faculty nttt fSaeadAtoatclSine Sagamihara, Science, Astronautical and Space of Institute eateto hsc,Ymgt nvriy aaaa Yam Yamagata, University, Yamagata Physics, of Department eateto hsc,TkoIsiueo ehooy Megu Technology, of Institute Tokyo Physics, of Department niuefrCsi a eerh nvriyo oy,Kash Tokyo, of University Research, Ray Cosmic for Insitute aut fEgneig hnh nvriy aao Nagan Nagano, University, Shinshu Engineering, of Faculty T aoaoy aoaUiest,Ngy,Aci464-8601 Aichi Nagoya, University, Nagoya Laboratory, STE eateto hsc,Gaut colo cec,KooUn Kyoto Science, of School Graduate Physics, of Department ± brk rfcua nvriyo elhSine,Ai Ib Ami, Sciences, Health of University Prefectural Ibaraki eateto hsc,TkiUiest,Hrtua Kanag Hiratsuka, University, Tokai Physics, of Department eateto hsc n ah hsc,Uiest fAdel of University Physics, Math. and Physics of Department eateto hsc,Kngw nvriy ooaa Kan Yokohama, University, Kanagawa Physics, of Department eateto hsc,OaaCt nvriy sk,Osaka Osaka, University, City Osaka Physics, of Department ainlAtooia bevtr fJpn iaa Tokyo Mitaka, Japan, of Observatory Astronomical National low aedtce gi h e am-a emission gamma-ray TeV the again detected have we Also 0 e rsmbyb h hc ceeainprocess. acceleration shock the by presumably TeV 100 1 µ n 0TVrespectively. TeV 50 and G 9 h ifrnilsetu fTVgmarays gamma TeV of spectrum differential The 1 .Mizumoto Y. , 6 .Naito T. , .Suzuki R. , .aioi([email protected] T.Tanimori 2 .Yoshida T. , 7 T.Tamura , 10 .Mori M. ,

c eekvtlsoe This telescope. Cerenkov ˇ oencsGslshf 2001 Gesellschaft Copernicus 3 ..Edwards P.G. , 14 7 .Yanagita S. , .Moriya M. , 4 6 3 .Gunji S. , .Muraishi H. , n .Yoshikoshi T. and , cy: in 52 Japan 8512, e aaaa2981,Japan 229-8510, Kanagawa rk 0-34 Japan 300-0394, araki 8-53 Japan 380-8553, o p- w 5-22 Japan 259-1292, awa ok,Tko1285,Japan 152-8551, Tokyo ro-ku, w,2788 hb,Japan Chiba, 277-8582 iwa, gt 9-50 Japan 990-8560, agata Japan , 8-58 Japan 181-8588, ie A50,Australia 5005, SA aide, 5-55 Japan 558-8585, st,Kf,Ymnsi4087,Japan 400-8575, Yamanashi Kofu, rsity, vriy yt,6680,Japan 606-8502, Kyoto, iversity, gw 2-66 Japan 221-8686, agawa 5 n e am asfo RXJ1713.7 from emission the rays X-ray in gamma synchrotron TeV detected strong and been Follow- bands; have energy SNRs two 1997). two same Yanagita other and success, Yoshida this ing 1996; Jager de (Mastichiadis background and cosmic K with 2.7 electrons of energy photons energy high low of In- (IC) the scattering on Compton based verse predictions with fluxes consistent integral were detected TeV The above CAN- 1998). by al. rim et (NE) (Tanimori East GAROO Norh the from emission gamma-ray oeaaefrom considered High emanate be could to SNRs protons). from detected (mainly rays rays gamma energy cosmic galactic of generation Beppo by 1999). al., (P¨uhlhofer detected et HEGRA and were 1997) al., A et Cassiopeia (Allen SAX from those and also respectively, 1997) 2001) al., and al., et et (Muraishi (Koyama 2001 in 1997 CANGAROO in by ASCA by detected were iheeg ii fprn lcrn,a neeg u peak flux energy an as electrons, a parent by explained of be limit can energy steepening high TeV this the scheme, around this spectra In steeper region. a to turnover a with band TeV hc rn nS 06i o euu ( tenuous too is 1006 SN th in around electrons density front energy matter shock the high because very 1998), by al. et radiation (Tanimori to IC the likely by On are explained SN1006 be from SNRs. detected in gamma-rays hand, protons other accelerated and matter swept-up as hsmdl nadto,peit atrpwrlaw power flatter a predicts addition, ( in gamma spectrum TeV model, and X-ray This radio, of data rays. the all fit can process IC naee l 96 ognrt h bevdTVgamma-ray TeV observed the from generate flux to 1996) al. et ingale .Hara S. , tpeet Nsaegnrlyfvrda h iefrthe for site the as favored generally are SNRs present, At 11 .Muraki Y. , 15 6 π .Hara T. , E ◦ eas h oe ae ntesynchrotron- the on based model The decays. s s dE CC2001 ICRC iha ne of index an with ) 9 π .Nishijima K. , ◦ 2 .Kawachi A. , easidcdb olsosbetween collisions by induced decays 12 7 ..Patterson J.R. , .Kifune T. , s − = 96(G347.3 3946 − ∼ ∼ . cm 0.4 1 7,8 . 6 Y. , eo the below − 13 3 Will- : K. , − 0.5) e 2 of the synchrotron emission appears in the soft X-ray region. At present there still remains the possibility that an addi- tional component of gamma rays from π◦ decays with the index of s ∼ −2.2 would gradually dominate instead of IC 900 400 radiation in the sub-TeV energy region. Detailed study of the 800 350 700 300 spectral shape is quite important to answer this question. 600 250 500 In April 2000, the new 10m CANGAROO-II telescope 200 400 was completed (Tanimori, et al., 2001 and Kawachi et al., 150 300 2001). Both RXJ1713.7−3946 and SN1006 have been de- 200 100 tected again with high statistics by this new telescope, and 100 (a) 50 (b) 0 0 their results are presented in this conference (Enomoto et al, 0 20 40 60 80 0 20 40 60 80 2001; Hara et al., 2001). Thus we have firmly established the Alpha(deg.) Alpha(deg.) TeV gamma-ray emission from SNRs. Here we report on a further analysis of the same data set 500 used in the previous paper (Tanimori et al., 1998) in order to 400 extract spectral information for the reasons mentioned above. 300

200 2 Analysis and Results 100 (c) 0 Observations were done with the 3.8m Cerenkovˇ imaging 0 20 40 60 80 telescope of the CANGAROO Collaboration (Patterson & Alpha(deg.) Kifune 1992; Hara et al. 1993) near Woomera, South Aus- ◦ ◦ tralia (136 47’ E and 31 06’S) in 1996 and 1997. Fig. 1. (a) Plot of the event distribution as a function of the α for the An imaging analysis using the conventional parameteriza- 1997 data. (b) is the similar plot after applying the timing selection, tion was applied for the data, where the parameter cuts were and (c) is also the similar plot for rejected events of (a) by this cut. varied as a function of the energy of the shower. These vari- ations of parameter cuts were estimated from the simulation study. Here we add the additional analysis to enhance the of the gamma-ray events remained. gamma-ray signals using the fine timing information of each In order to obtain the differential energy spectrum, sev- hit photomultiplier (about 1 ns time resolution). An image eral α plots were made by varying the minimum and maxi- obtained from a single telescope gives only one angle: the mum numbers of detected Cerenkovˇ photons. The collecting projection of the direction on a plane normal to the mirror area, trigger efficiency and threshold energy corresponding axis. Stereo observations are required to obtain the two an- to each α plot were independently estimated from simula- gles necessary to determine the direction of a gamma-ray in tions, where events were generated between 1 TeV and 30 the sky. For a single telescope, the unknown angle of the TeV using an initial power-law index of ∼−2.2. This initial shower direction lies on the extended long axis of the ellip- index was estimated from two observed integral fluxes 1996 tic image. However if we know both the the arrival times and (threshold energy ≥ 3.5 TeV) and 1997 (threshold energy ≥ the height of the emitting point of Cerenkovˇ light, the shower 1.7 TeV). At first the spectrum was calculated using this in- direction can be reconstructed. In other words, the arrival dex. Using the new index obtained from this spectrum, the timings of Cerenkovˇ photons emitted from a shower can be above simulations and calculations were iterated in a similar calculated if both the height and the direction of a shower are way until the spectral index converged. Other initial values known. Therefore, by assuming the height of shower gen- of the index were tried, but the resultant index converged to eration (∼ 10, 000m is used here), the unknown angle can the similar value. ˇ be calculated using the observed arrival timings of Cerenkov The resultant differential spectrum, J(E), between 1.5 TeV photons. The simulation indicated that the obtained angle is and 20 TeV is plotted in Fig.2 in which only statistical errors quite insensitive to the assumption of the shower height. We are shown. It can be written as: actually obtained this angle event by event, and used it only to estimate the direction of the shower development along the −11 E −2.3±0.2 −1 −2 −1 J(E)=(1.1±0.4)×10 ( ) TeV cm s ,(1) long axis of the image. By selecting the shower images for 1T eV events from the target direction, gamma-ray events should be enhanced, and indeed α peaks in boththe 1996and 1997data where errors quoted are statistical. Systematic errors are esti- became more significant. The result of the 1997 data is pre- mated to be ∼30%. Differential fluxes of 1996 and 1997 data sented in Figs. 1a and b and the significance of the α peaks were also individually calculated, and those fluxes were con- was increased from 7.5 σ to 8.8 σ by this method. After this sistent within one sigma errors. Furthermore the differential cut, about 70% of background events that survived through flux without applying the timing reconstruction method was imaging selections can be further removed, while about 75% obtained to be quite similar to this flux. 3

-22 the convex shape. The resultant field strength means the 10 Emax ∼ 51 TeV from the above formula. 3 µG Figure 3 shows the wide band energy spectrum from the 2 µG 1 µG ] radio to TeV regions at the north rim of SN1006 and also -1 -23 4 µG 10 6 µG the above fitting (Naito et al., 1999). All data are fitted very eV

-1 8 µG well all over the wide band. Here an allowable spectrum due

s ◦

-2 to π decays generated by high energy proton is also plotted -24 10 taking into account the upper limit of GeV gamma-rayfluxes, which obviously conflicts with the observed TeV gamma-ray flux. Those upperlimits in the GeV region of this figure were -25 10 calculated from the EGRET archive data using a maximum

[photons cm likelihood method (Mattox et al., 1996). These results seems

ε very reasonable consequence considering the tenuous shell -26 −3 10 of (≤∼0.4 cm ; Willingale et al. 1996).

dF/d Thus the identification of the parent particles of the TeV gamma-rays (electron or proton) will be possible by observ- -27 ing the wide spectrum from sub- to multi-TeV regionas shown 10 12 13 in Fig.3. A Gamma-ray spectrum flatter than −2.0 in this re- 10 10 gion is surely due to the I.C. process, while that due to π◦ photon energy ε [eV] decay generated by collision between ISM and high energy proton is expected to be steeper than −2.0. Fig. 2. Observed differential spectrum from the NE rim, where only From this fitting, the energy flux of TeV gamma rays due statistical errors are shown. In addition, spectra of TeV gamma rays to the IC process at 300 GeV is estimated to be ∼ 4 eV for several magnetic field strengths B calculated by the IC model cm−2s−1, whichis lowerthantheupperlimit (8 eVcm−2s−1 are presented. at 300GeV)reportedin Chadwick et al. (2000). Also we have observed again TeV gamma rays from the NE rim of SN1006 observed by the 10m telescope as presented in this confer- 3 Discussion ence (Hara et al., 2001). Although the new data is still very The resultant differential spectrum shows a spectral index of preliminary, the spectra are consistent in the overlapping en- −2.3 between 1.5 TeV and 20 TeV. The profile of the ob- ergy region. served differential spectrum may resolve which IC or π◦ de- In this conference, we present the interesting spectrum of cay mechanism dominates this TeV emission. The model TeV gamma rays of PSR1706−44 showing an obvious break used here originated in Yoshida and Yanagita (1997), simi- of the spectrum around 1 TeV (Kushida et al., 2001). This lar to the previous paper, while we improved the treatment breaking can be fit by the IC model quite well similar to the of the approximation for the synchrotron flux (Naito et al., case of SN1006. These two results indicate the major role of 1999): a simple δ function centered on the characteristic IC process for the production of high energy gamma rays in frequency was changed to the full function describing syn- the universe. chrotron emission (Rybicki and Lightman 1979 and refer- However, our result showing the IC process dominance for ences therein). The electron spectrum was assumed to be as the production of TeV gamma rays in SN1006 does not im- ply the nonexistence of the plenty of high energy protons ac- s dNe E E celerated by a diffusive shock. Conventional diffusive shock = N0 2 exp − . (2) dE mec   Emax  theories naturally predict efficient acceleration of protons in SNRs. In SN1006, a deficiency of target protons surround- In addition, the latest data of both radio (Reynolds and El- ing the SNR may conceal accelerated protons from being re- lison 1992) and X-rays from ROSAT and ASCA (Willingale vealed. et al. 1996; Ozaki 1998) were used in fitting the synchrotron RXJ1713.7−3946 is the second SNR emitting both syn- flux. In particular, we used the partial fluxes emitted from chrotron X-rays and TeV gamma rays. However, the en- the NE rim for all data: for the ASCA data the partial flux vironment of RXJ1713.7−3946 is obviously different from has been directly obtained by Ozaki (1998), and for the radio SN1006; for example, association of RXJ1713.7−3946 with and the ROSAT data the partial fluxes were estimated using a molecularcloud is reported from the CO observation( Slane the morphological distributions of radio (Winkler and Long et al., 1999). Our new result on RXJ1713.7−3946 indicates a 1997) and ROSAT (Willingale et al. 1996) data. Comparing more intense flux below 1 TeV,while the flux above 1 TeV of the calculated synchrotron flux to the observational data, we it is similar to that of SN1006 (Enomoto et al., 2001). Those obtain s = −2.2 and (Emax/TeV ) (B/µG) = 101. In features looks inconsistent with the scenario of the IC model. Fig.2 several expected IC spectra arep plotted by varying the Future analyses of SN1006 and RXJ1713.7−3946 will soon strength of the ambient magnetic field, where B =4 ± 1µG provide both the differential spectra between ∼ 300 GeV and is most probable, and this model matches successfully in ∼ 10 TeV and the image of TeV gamma ray emission respec- 4

2 10 IRAS EGRET upper limit upper limit ROSAT CANGAROO 1

] 10 -1

s ASCA -2

0 [ eV cm 10 ε Inverse Compton dF/d 2 B = 4µ G ε Synchrotron -1 radio or 10 0 ν π decay S

ν parent proton spectrum α=2.2 Emax =5e15 -2 10 no*E50=2.5

-5 0 5 10 15 10 10 10 10 10 photon energyε [eV]

Fig. 3. Multi-band spectrum of energy fluxes observed from the NE rim, where observed fluxes or upper limits of radio (Reynolds 1998), infrared, soft X-ray (estimated from Willingale et al. 1996), hard X-ray (Ozaki 1998), GeV gamma rays (calculated from the EGRET archive data), and TeV gamma rays are presented. Solid lines are the fits based on the model of IC model and neutral pion decay. tively, which will tell us what kind of particle is accelerated Naito, T., et al., Aston. Nachr.,320, 205-206, 1999. and generates TeV gamma rays in each SNR. Ozaki, M. 1998a, private communication Ozaki, M. 1998b, ph. D thesis. Acknowledgements. Patterson, J. R., & Kifune, T., Australian and New Zealand Physi- cist, 29, 58, 1992. This work is supported by a Grant-in-Aid in Scientific Re- P¨uhlhofer, G., et al. Proc. of the 26th ICRC (Salt Lake City), 3, search of the Japan Ministry of Education, Culture, Science, 492-495, 1999. Reynolds, S.P., & Ellison, D.C., Astrophys. J., 399, L75-L78, 1992. Sports and Technology, and the Australian Research Council. Rybicki, G.B., & Lightman, A.P., 1979, Radiative Process in Astro- physics (Wiley, New York), 1979. Slane, P., et al., Astrophys. J., 525, 357-367, 1999. References Susukita, K., Memoirs of the Faculty of Science, Kyoto University, Series A of Physics, Astrophysics, Geophysics and Chemistry Allen, G.E., et al., Astrophys. J., 487, L97-L100, 1997. Vol.39 No.5, 1998. Chadwick, P.M. et al., in Proc. of GeV-TeV Gamma Ray Astro- Tanimori, T., et al., Astrophys. J., 497, L25-L28, 1998. physics Workshop edited by Dingus et al., AIP Conf. Proc. 512, Tanimori, T. to appear in Prog. Theor. Phys., 143, 2001. 210-214, 2000. Weekes, T.C., et al., Astrophys. J., 342, 379, 1989. Enomoto, R., et al., in these Proceedings, 2001. Willingale, R., et al., Mon. Not. R. Astron. Soc., 278, 749, 1996. Hara, S., et al., in these Proceedings, 2001. Winkler, P., F., & Long, K., S., Astrophys. J., 491, 829, 1997 Hara, T., et al. Nucl. Instr. and Meth., A332, 300, 1993. Yoshida, T., & Yanagita, S., in Proc. 2nd INTEGRAL Workshop Kawachi, A., et al., Astropart. Phys. 14, 261-269, 2001. ‘Transparent Universe’, ESA SP-382, 85-88, 1997. Koyama, K., et al. Nature, 378, 255-258, 1995 K. Koyama et al, Publ. Astron. Soc. Japan,49, L7, 1997. Kushida, J. et al., in these Proceedings, 2001. Mastichiadis, A., & de Jager, O.C., Astron. Astrophys., 311, L5-L8, 1996. Mattox, J.R., et al., Astrophys. J., 461, 396, 1996. Muraishi, H., et al., Astron. Astrophys., 354, L57-L60, 2001.