arXiv:2003.04795v1 [astro-ph.HE] 10 Mar 2020 1 3 2 4 5 6 colo srnm n pc cec,NnigUniversity, Nanjing Science, Space and Astronomy of School h ers uenv N18A a enpredicted been ( has while 1987A, long a SN for Thus, supernova nearest same. from the the particularly is govern- supernovae, from gamma-rays physics emission of gamma-ray the production amount (CSM), the large ing a medium up circumstellar swept su- already early of have the our they from in as distinct pernova (SNRs) are SNRs remnants Although supernova from Galaxy. historic detected of been have dozen energies a TeV to GeV from sion 2011 e aoaoyo akMte n pc srnm,Prl Mo Purple Astronomy, Space and Matter Dark of Laboratory Key Depar Cosmology, China and 230026, Galaxies Anhui in Research for Labrotory Key ayhve al. et Malyshev h e msinfo h N18Argo sobserved as region 1987A SN with the from emission GeV the eateto hsc n srnm,Uiest fNevada, of University Astronomy, and Physics of Department e aoaoyo oenAtooyadAtohsc (Nanji Astrophysics and Astronomy Modern of laboratory Key colo srnm n pc cec,Uiest fScienc of University Science, Space and Astronomy of School n hycnpoueGVTVgmary i the via gamma-rays ( GeV-TeV decay pion produce neutral can shocks they the by su- accelerated and of be number also large may a Protons from pernovae. observed been have en- as X-ray ergies, through emit radio-submillimeter which collisionless the energies, from The high photons to electrons forms. accelerates shock shock collisionless a medium, rpittpstuigL using typeset Preprint stesproaeet xad ntecircum-stellar the in expands ejecta supernova the As EEDPTU ICVR FGVGMARYEISO RMSU FROM EMISSION GAMMA-RAY GEV OF DISCOVERY SERENDIPITOUS A , Fermi 2015 oniec ihtescn ers uenv–N20d rmo from G 2004dj fading supernova–SN a with nearest of galaxies discovery second star-forming e serendipitous the supernova a nearest with report the conver we coincidence of is Here flux energy gamma-ray picture. kinetic low supernova this The the time. of early fraction via very small emission a gamma-ray only GeV-TeV since produce electron to the of expected radio evidence are The the constitutes protons). explosion and supernova electrons the (i.e., circum-stella the particles and accelerates ejecta which supernova a between interaction The 6 urudn ihdniysel hc eeeae h jcaadc and ejecta the decelerates which shell, high-density surrounding icvr e msin rmtosa-omn aais .. r 299 survey Arp i.e., gamma-ray galaxies, our star-forming addition, two from In emissions GeV protons. discovers relativistic to energy hoQagXi Shao-Qiang × ; LT u h oaino hssuc over- source this of location the But –LAT. age l 2019 al. et Wang 10 47 r.W nepe hsgmaryeiso rsn rmtesuperno the from arising emission gamma-ray this interpret We erg. ( 2019 1. eehoe l 2011 al. et Berezhko INTRODUCTION UVYO ERYSA-OMN AAISWT FERMI-LAT WITH GALAXIES STAR-FORMING NEARBY OF SURVEY uaee l 2011 al. et Murase A eotdarcn nacmn of enhancement recent a reported ) T E tl ATX .1.0 v. AASTeX6 style X 1,5 u-uLiu Ruo-Yu , .Ide,gmaryemis- gamma-ray Indeed, ). Fermi , ; 2015 1,5 eehoe al. et Berezhko LT h oa am-a nryrlae yS 04ji about is 2004dj SN by released energy gamma-ray total The –LAT. in-uWang Xiang-Yu , .Recently, ). n ehooyo hn,Hfi nu 306 China 230026, Anhui Hefei, China, of Technology and e ajn 103 hn;rlunueuc;[email protected] [email protected]; China; 210023, Nanjing ABSTRACT a ea,40 ayadPrwy a ea,N 89154-4002 NV Vegas, Las Parkway, Maryland 4505 Vegas, Las mn fAtooy nvriyo cec n ehooyof Technology and Science of University Astronomy, of tment gUiest) iityo dcto,Nnig202,Chi 210023, Nanjing Education, of Ministry University), ng nanOsraoy hns cdm fSine,Nnig2 Nanjing Sciences, of Academy Chinese Observatory, untain 1,5 u-h Yang Rui-Zhi , eat ihtecruselrmdu CM created (CSM) in- medium ejecta circumstellar stellar the moving with fast teracts The red-supergiants. be the of circle error there the in However, quasar iPTF14hls. a is supernova, consistent peculiar temporally a and with spatially source gamma-ray n aaee agsepce yteter.Recently, interest- theory. reach the by to expected start ranges parameter limits ing observational their but env xlddi h aayNC20 tadis- a at 2403 ( NGC 2004 Mpc 3.5 galaxy about the of tance in exploded pernova source. variable con- gamma-ray its the address with conclusively infrared nection to this difficult the of it observations makes to multi-wavelength quasar of according lack candidate The blazar data. a is which al. et Yuan ( al. et Ackermann rmsprlmnu N aebe erhdby searched been have SNe super-luminous clarified. be from to remains emission GeV the so of sources, gamma-ray nature potential the other several with laps 2018 N20d stenaetadbihetTp I su- IIP Type brightest and nearest the is 2004dj SN am-a msinfo yeInSeand SNe IIn Type from emission Gamma-ray .Tepoeioso N I r huh to thought are IIP SNe of progenitors The ). ,rsetvl.N vdnefrasga a found, was signal a for evidence No respectively. ), pp ceeain h ceeae protons accelerated The acceleration. olsos u h u suulylow usually is flux the but collisions, eimdie togsokwave shock strong a drives medium r onverts ( 2018 n -a msinosre after observed emission X-ray and rgmarysre fnearby of survey gamma-ray ur poin N18A geswith agrees 1987A, SN xplosion, e notesokeeg tthe at energy shock the into ted 2,6 Vgmarysuc nspatial in source gamma-ray eV fnab trfriggalaxies star-forming nearby of n 3,frtefis time. first the for M33, and in Yuan Qiang , eotdtedtcino variable a of detection the reported ) ( ∼ 2015 aeet neatn iha with interacting ejecta va %o h jcaskinetic ejecta’s the of 1% aaoe l 2004 al. et Nakano and ) ENV 04JI A IN 2004DJ PERNOVA 3 igZhang Bing , eal-iac tal. et Renault-Tinacci Fermi .cn LTsource, –LAT na ; aa tal. et Patat hn,Hefei, China, 03,China 10033, 4 2 Xi et al. by the stellar wind of the red-supergiants. SN 2004dj in a spatial bin of 0.025◦ per pixel. We make use of was detected in a wide range of wavelengths from recent developments of the Science Tools (v11r5p3) for radio, through infrared to X-rays during the first likelihood analysis. The background model for each ce- several years after the explosion (Meikle et al. 2011; lestial ROI contains all sources listed in the 4FGL along Chakraborti et al. 2012; Nayana et al. 2018). with the standard diffuse emission background, i.e. the In this paper, we report a serendipitous discovery of foreground for Galactic diffuse emission (gll iem v7.fits) a GeV gamma-ray source coincident in position with released and described by the Fermi–LAT collabora- SN 2004dj. This discovery is made when we search for tion through the Fermi Science Support Center (FSSC) GeV emission from nearby star-forming galaxies. By (Acero et al. 2016) and the background for spatially analysing the 11.4 years of Fermi–LAT data, we find isotropic diffuse emission with a spectral shape described a GeV source from the direction of the nearby galaxy by iso P8R3 SOURCE V2 v01.txt. In the likelihood NGC 2403. The flux of the this source, however, does analysis, we allow all sources that are separated from not obey the well-known relation between the gamma- ROI center by less than 6.5◦ to have a free normaliza- ray luminosity and infrared luminosity for star-forming tion to allow for the broad PSF at the lowest energies galaxies, disfavoring the usual cosmic ray-ISM interac- (the 68% containment radius of photons at normal inci- tion origin. We also find the flux is decaying during dence with an energy of 300 MeV is roughly 2.5◦). This the Fermi–LAT observation period from 2008 to the choice ensures that 99.9% of the predicted gamma-ray present. Motivated by this, we search for transient counts is contained within the chosen radius. The nor- sources in the error region of the gamma-ray emission malizations of both Galactic and extragalactic diffuse and find that SN 2004dj lies within the error region. emission models are left free. We employ the gttsmap tool to evaluate the 6◦ × 6◦ map of the Test Statistic 2. SAMPLE SELECTION AND DATA ANALYSIS (TS), defined as TS = −2(lnL0 − lnL), where L0 is We perform a search for possible gamma-ray emis- the maximum-likelihood value for null hypothesis and sion from galaxies in the IRAS Revised Bright Galax- L is the maximum-likelihood with the additional point ies Sample (Sanders et al. 2003), using 11.4 years of source with a power-law spectrum, and then survey the gamma-ray data taken by the Fermi–LAT telescope. new gamma-ray sources when the significance of the This is a complete flux-limited sample of all extra- gamma-ray excess above backgrounds is TS > 25. The galactic objects brighter than 5.24 Jy at 60 µm, cov- best location and uncertainty of the new sources can be ering the entire sky surveyed by IRAS at Galactic determined by maximizing the TS value and using the latitudes |b| > 5◦. In the sample, 15 infrared(IR)- distribution of Localization Test Statistic (LTS), defined bright galaxies have been detected in gamma-rays with by twice the log of the likelihood ratio of any position Fermi–LAT and listed in Fermi–LAT Fourth Source with respect to the maximum. We claim a gamma-ray Catalog (4FGL,The Fermi-LAT collaboration (2019)), source associated with a target galaxy when the target including six galaxies with AGNs (i.e., Cen A, IC galaxy lies within the 95% confidence location region 4402, NGC 3067, NGC 3683, NGC 1275, NGC of the gamma-ray source. Note that new background 3424 (Peng et al. 2019)), seven star-forming galax- sources could be found and, if so, we re-do the analysis ies (SMC, LMC, M31, NGC253, M82, NGC 2146, using an updated background model including the new Arp 220(Abdo et al. 2010a; Ackermann et al. 2012; background sources. Tang et al. 2014; Peng et al. 2016; Griffin et al. 2016)), and two star-forming galaxies with obscured AGNs 3. RESULTS (NGC 1068, NGC 4945,Ackermann et al. 2012) . These Our analysis results in the detection of two new galaxies are excluded from our sample for the gamma- gamma-ray sources that are, respectively, spatially co- ray emission search. incident with the galaxy Arp 299 and the galaxy M33 This work uses ∼ 11.4 years (MET 239557417- (see Table. 1, more details will be presented in an 595385929) of Fermi–LAT SOURCE class events with upcoming paper, Xi et al. 2020, in prep). In ad- reconstructed energies between 300 MeV and 500 GeV, dition, we find that the gamma-ray source reported excluding those with a zenith angle larger than 90◦ to in 4FGL (4FGL J0737.4+6535) is spatially coincident avoid Earth limb contamination. We implement a stan- with the galaxy NGC 2403, since the best-location of dard sequence of analysis steps. For each IR-bright 4FGL J0737.4+6535 is within the spatial extended re- galaxies, we select the events in a 17◦ × 17◦ region gion (∼ 0.1◦) of NGC 2403. For the rest galaxies in of interest (ROI) centered at the galactic center and our sample, we do not find any significant gamma-ray use gtmktime tool to select time intervals expressed by detection. (DATA QUAL > 0) && (LAT CONFIG ==1). We bin It has been found that there is an empirical cor- the data in 20 logarithmically spaced bins in energy and relation between the gamma-ray luminosity in (0.1 − GeV emission from SN 2004dj 3
Table 1.
Name dL Optical center gamma-ray location F0.1−100GeV Γ TS [Mpc] (R.A.,Dec) (R.A.,Dec) [10−12 erg cm−2 s−1 ] (1) (2) (3) (4) (5) (6) (7) SN 2004dj 3.49 (114.321◦, 65.599◦) (114.395◦, 65.573◦)±0.041◦a 2.29±0.51a 1.92±0.07a 73.1a Arp 299 47.74 (172.136◦, 58.561◦) (172.050◦, 58.526◦)±0.111◦ 1.08±0.28 2.07±0.20 27.8 M 33 0.84 (23.475◦, 30.669◦) (23.609◦, 30.784◦)±0.089◦ 1.28±0.42 2.33±0.24 25.1 aThese values are derived by the first 5.7 years analysis.
Note—(1) source name; (2) Optical Right ascension J2000; (3) Optical Declination J2000; (4) Best location of the γ − ray source detected by Fermi–LAT. The uncertainty of position correspond to 95% containment radius, which derived by fitting the distribution of LTS to a 2D Gaussian function; (5) 100 MeV-100 GeV γ − ray average fluxes; (6) Power-law spectral photon index derived by broad band spectrum fitting; (7)TS value of the gamma-ray event excess using Fermi–LAT data above 0.3 GeV at the source position.
100) GeV and total IR luminosity (8 − 1000 µm) for local group galaxies and nearby star forming galaxies (Ackermann et al. 2012; Peng et al. 2016). This cor- 1 ()⊙ ) ) ⊙ ) relation is generally interpreted as that the gamma- 1 ( ( ray emission arises from cosmic-ray (CR) protons in- ⊙ ⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙ 1 1 )⊙ ) ) teracting with the interstellar medium (ISM) via the ⊙ ) ⊙ ⊙ ⊙1 ⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙ proton–proton (pp) collisions (Ackermann et al. 2012; 1 1 ⊙ 1 ⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙ ⊙ ⊙( Zhang et al. 2019). In Figure 1, we show the rela- ⊙ ⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙ ⊙ ⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙ tion between the gamma-ray luminosity and IR lumi- 1 ⊙ ⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙ nosity for the three gamma-ray sources found in our ⊙ 1 work. While Arp 299 and M33 are well consistent ⊙ 1⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙⊙