arxiv: v1 [astro-ph.he] 13 Jul 2018

Draft version July 16, 218 Typeset using L A TEX manuscript style in AASTeX62 Active galactic nuclei with GeV activities and the PeV neutrino source candidate TXS 56+56 arxiv:187.521v1 [astro-ph.he] 13 Jul 218 Neng-Hui Liao, 1 Yu-Liang Xin, 1 Yun-Feng Liang, 1 Xiao-Lei Guo, 1, 2 Shang Li, 1 Hao-Ning He, 1, 3 Qiang Yuan, 1 and Yi-Zhong Fan 1 1 Key Laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 2134, China 2 School of Astronomy and Space Science, University of Science and Technology of China, Hefei 2326, Anhui, China 3 Astrophysical Big Bang Laboratory, RIKEN, Wako, Saitama, Japan ABSTRACT On 217 September 22 the IceCube neutrino observatory detected a track-like, veryhigh-energy event (IceCube-17922A) that is spatially associated with TXS 56+56, a quasar at a redshift of z =.3365. This source is characterized by the increased acitivies in a very wide energy range (from radio to TeV) during these days. To investigate the possible connection of the PeV neutrino emission with the GeV activity of blazars, in this work we select 116 bright sources and analyze their lightcurves and spectra. We focus on the sources displaying GeV activities. Among these blazars, TXS 56+56 seems to be typical in many aspects but is distinguished by the very strong GeV activties. We suggest to search for neutrino outburst in the historical data of IceCube, as recently done for TXS 56+56, from the directions of these more energetic and harder blazars with strong GeV activities. Keywords: galaxies: active gamma rays: galaxies blazars Corresponding author: Neng-Hui Liao, Hao-Ning He, Yi-Zhong Fan liaonh@pmo.ac.cn, hnhe@pmo.ac.cn, yzfan@pmo.ac.cn

2 1. INTRODUCTION The discovery of astrophysical neutrino flux around PeV energies by IceCube (Aartsen et al. 217) is an important milestone in high energy astronomy (IceCube Collaboration 213). The neutrino is an ideal astronomical messenger since they travel undistorted from the sources and are therefore a valuable probe of the innermost regions of the energetic and enigmatic objects in the cosmos. Since then, the accumulating observations of neutrino events from IceCube suggest that a significant fraction of the observed neutrinos are of extragalactic origin due to their isotropic distribution and reveal a flux of neutrinos with a total energy density comparable with that of the extragalactic rays observed by Fermi-LAT (Aartsen et al. 214, 215; Ackermann et al. 216). Blazars are an extreme subtype of Active Galactic Nuclei (AGNs), dominating the extragalactic -ray sky (Acero et al. 215; Madejski & Sikora 216). Emissions from their strong collimated jets are overwhelming due to relativistic beaming effects and hence blazars are characterized by the luminous and highly variable broadband continuum emissions (e.g., Blandford & Rees 1978; Ulrich et al. 1997). Blazars are traditionally divided into flat spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs) based on their optical spectra. Spectral energy distribution of these jetted AGNs generally exhibits a two-bump structure in logνf ν logν representation and the low energy bump is widely believed to be contributed by synchrotron emission. Therefore, blazars are also classified as low-synchrotron-peaked sources (LSPs, ν syn peak < 114 Hz), intermediate-synchrotronpeaked sources (ISPs, 1 14 Hz < ν syn peak < 115 Hz), and high-synchrotron-peaked sources (ν syn peak > 1 15 Hz)(Ackermann et al. 215). However, the origin of the other bump extending to -ray regime is still under debate. On one hand, the leptonic radiation model (e.g., Maraschi et al. 1992; Dermer & Schlickeiser 1993; Sikora et al. 1994; B lażejowski et al. 2) in which the rays are from inverse Compton scattering of soft photons by the same population of electrons that emit the synchrotron emission, can naturally describe the tightly connected multiwavelength variability of blazars (e.g., Abdo et al. 21; Liao et al. 214). On the other hand, observational phenomena, like the orphan -ray flare (e.g., Böttcher 25), support the hadronic scenario that rays and neutrinos are produced via interactions of high-energy protons with gas (i.e., the pp-interactions) in the jets (Schuster et al.

3 22) or in interactions of protons with internal (Mannheim 1995) or external (Atoyan & Dermer 21) photon fields (p-interactions). Among the various possibilities of potential extragalactic sources of neutrinos (see Ahlers & Halzen (215) for a review), including star-forming galaxies (e.g., Loeb & Waxman 26), gamma-ray bursts (e.g., Waxman & Bahcall 1997), galaxy clusters (e.g., Berezinsky et al. 1997) and so on, blazars are believed to be promising sources(e.g., Tavecchio & Ghisellini 215; Halzen & Kheirandish 216). This kind of sources also contribute to the majority of extragalactic -ray background which is consistent with the measured neutrino flux level (Murase et al. 213). Meanwhile, a correlation between cosmic neutrinos and blazar catalogs has been argued (Padovani et al. 216). Moreover, since strong -ray flares of blazars have been frequently detected, a spatial association combined with a coincidence in time with a flaring blazar may represent a smoking gun for the origin of the IceCube flux. A coincidence between a 2 PeV neutrino event and the blazar PKS B1424 418 provides an interesting hint in this context (Kadler et al. 216). On 217 Spetmber 22, the IceCube neutrino observatory (hereafter IceCube) detected a track-like, very-high-energy event with a high probability of being of astrophysical origin (Kopper & Blaufuss 217). Inside the error region of the neutrino event, there is a Fermi-LAT source, TXS 56+56 (Ajello et al. 217). Interestingly, the GeV emission is found to be at high state (Tanaka et al. 217). Enhanced emission are also found in radio, optical, X-ray and TeV bands (Tanaka et al. 217; Fox et al. 217; Franckowiak et al. 217; Mirzoyan 217). Therefore, TXS 56+56 is a promising PeV neutrino source candidate (IceCube collaboration 218; IceCube collaboration et al. 218). Motivated by the coincidence of the GeV activity with the PeV neutrino emission, in this work we analyze the Fermi-LAT data (Atwood et al. 29, 213) of some bright AGNs and identify the GeV flares. We then compare the properties of these bright AGNs displaying strong GeV activities with TXS 56+56 to look for possible indication of the high energy neutrino sources. This work is organized as follows: in Section 2 we introduce our sample and the data analysis. The results are reported in Section 3. We summarize our results with some discussions in Section 4. 2. THE SAMPLE AND DATA ANALYSES

4 2.1. The Sample Blazars included in the Fermi-LAT monitored source list 1 are chosen for this study, in which the sources are selected if their instantaneous weekly fluxes are above 1 6 ph cm 2 s 1. These sources represent the brightest and the most variable ones among the Fermi-LAT blazars. Another advantage is that the sources are randomly distributed in the sky since Fermi-LAT performs an all-sky survey. Several sources (e.g. Fermi J27 2518) are excluded in our analysis because there are only a few flux data points in the preliminary weekly light curves provided by the Fermi collaboration. Meanwhile, two sources (i.e. NRAO 19 and B2 1144+4) are not considered due to the ambiguous association relationship between the -ray source and its low-frequency counterpart. Therefore, our sample consists of total 116 blazars, including 81 FSRQs, 22 BL Lacs and 13 sources with unknown optical spectral property (BZU), according to BZCAT 2 (Massaro et al. 29). Or alternatively, there are 96 LSPs, 1 ISPs and 8 HSPs, according to the third Fermi-LAT AGN catalog (3LAC, Ackermann et al. 215), together with 2 sources lacking of relative information. The redshift distribution of our sample is between.31 (Mrk 421) and 2.852 (PKS 438 43), and there are seven sources (either BL Lacs or BZUs) lacking redshift imformation. The basic information of the sources in our sample are listed in Table 1. 2.2. Fermi-LAT Data reduction In the analysis, the latest Pass 8 version of the Fermi-LAT data (Atwood et al. 213) with Source event class (evclass = 128 & evtype = 3) are selected, recorded from May 1, 21 (Mission Elapsed Time 29436482; the time of the operation for the full IceCube detector) to May 1, 218 (Mission Elapsed Time 54682565) 3. We choose the events within a 1 region of interest (ROI) with energies between 1 MeV and 5 GeV in this analysis. In order to reduce the contamination from the Earth Limb, the events with zenith angles larger than 9 are excluded. In addition, the entire data set is filtered with gtmktime to obtain high-quality data in the good time intervals, with the expression 1 https://fermi.gsfc.nasa.gov/ssc/data/access/lat/msl lc/ 2 http://www.asdc.asi.it/bzcat/ 3 Note there is 2 day data gap for Fermi-LAT around March 218.

5 recommended by the LAT team, namely (DATA QUAL>)&&(LAT CONFIG==1). The data are analyzed with the standard LAT analysis software, ScienceTools version v1rp5 4, available from the Fermi Science Support Center, and the P8R2 SOURCE V6 instrumental response functions are adopted. First, we perform the unbinned likelihood analysis with gtlike to extract the global flux and spectral parameters of the target source. For the background subtraction, the Galactic diffuse emission and the isotropic diffuse emission are modeled by gll iem v6.fits and iso P8R2 SOURCE V6 v6.txt, which can be found from the Fermi Science Support Center 5. Meanwhile, all sources in the preliminary LAT 8-year Point Source List (FL8Y 6 ) within a radius of 15 from the ROI center are included in the source model. For convenient comparison between different sources, the spectral template of each target source is set to a power-law function (i.e. dn/de E Γ, where Γ is the spectral photon index). And during the fitting analysis, the normalizations and spectral parameters of all sources within a distance of 1 from the ROI center, together with the normalizations of the two diffuse backgrounds, are left free. The significance of the target source can be quantified by the test statistic (TS) value, which is defined as TS = 2 ln(l /L)(Mattox et al. 1996), where L and L are the maximum likelihood values of the null hypothesis and the tested model including the target source. The fittings are demanded to converge (i.e. fit quality = 3 ) to make sure the results are valid. The best-fit results for all target sources in the sample are summarized in Table 1. In the temporal analysis, we divide the data into 417 equal time bins and repeat the likelihood fitting for each time bin to extract a weekly light curve for each target. Considering the targets studied here are the brightest -ray sources in the extragalactic sky, we fix the spectral indexes of all background sources to the global fitting values, and only free the normalization parameters of background sources within 1 from the targets and two diffuse backgrounds. During the fitting analysis for each time bin, the weak background sources (i.e., TS < 5) are removed from the source 4 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 5 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/backgroundmodels.html 6 https://fermi.gsfc.nasa.gov/ssc/data/access/lat/fl8y/

6 model. Note that for any time bin in which the TS value of target source is smaller than 9, the 95% upper limit is calculated instead. 3. RESULTS 3.1. Global properties The analysis results of the entire 8-year LAT data are summarized in Table 1 as well. We also plot the photon index versus the photon flux and the -ray apparent luminosity (from.1 to 5 GeV, with handled k-correction) diagrams, see Fig. 1. The average photon indexes of FSRQs, BL Lacs and BZUs in our sample are 2.4±.17, 2.4±.18 and 2.29±.12, respectively, which are in agreement with the results in 3LAC (Ackermann et al. 215). Note that there is only one FSRQ, VER 521+211, whose photon index is smaller than 2. As shown in Fig. 1 and Fig. 2, FSRQs are generally brighter and more luminous than BL Lacs, and TXS 56+56 appears to be typical among other BL Lacs. The spectra of FSRQs are softer than that of BL Lacs, which may lower the neutrino detection possibility. 3.2. Temporal Behaviors 3.2.1. TXS 56+56 The weekly -ray light curve of TXS 56+56 is presented in the left panel of Fig. 3. Before MJD 57855, the source maintains at a relatively quiescent state, with only a few bins whose fluxes reach roughly three times of the 8-year averaged flux, 9.8 1 8 ph cm 2 s 1. However, since then, a strong -ray flare has appeared, with a peak flux of (5.8 ±.6) 1 7 ph cm 2 s 1 which is nearly 6 times of the 8-year averaged flux. The corresponding peak apparent -ray luminosity of the bin is 1.5 1 47 erg s 1, adopting a redshift of.3365 (Paiano et al. 218). Meanwhile, the photon index then is 2.5±.7, suggesting no significant spectral variability than the average value, 2.7±.1. The flare peaked on MJD 57981 (16th Aug. 217), about a month before the arrival time of the IceCube neutrino event. A further 3-day -ray light curve (the right panel of Fig. 3) shows that several sub-flares constitute this activity phase. Maybe there is a weak flare coincident with the neutrino event, but the error bars are relatively large.

7 3.2.2. Comparison between TXS 56+56 and Fermi-LAT bright blazars In 3LAC, about 69% FSRQs are found to be significantly variable, however, this fraction is down to 23% for BL Lacs (Ackermann et al. 215). Together with the fact that the former kind is generally brighter than the later one, it is not surprised that the vast majority of our sources are FSRQs. One famous case is CTA 12 (Becerra et al. 216; Xu et al. 216), whose daily LAT > 1 MeV flux is up to 1 5 cm 2 s 1, also see Fig. 4 for its weekly light curve. In addition to the brightness, its large variability amplitude is also remarkable. The weekly peak flux is more than one order of magnitude higher than the averaged flux. Examples of light curves of other subtypes are also shown in Fig. 4. Although most of our FSRQs are LSPs, there are several ISP FSRQs. No significant difference is found between their weekly -ray light curves. It is interesting to see that the variability amplitude for FSRQs is generally higher than HSP BL Lacs. However, for ISP BL Lacs and LSP BL Lacs, their -ray variability could be as violent as FSRQs. Worthy to note that there is one source, ON 246, which is also a ISP BL lac similar with TXS 56+56. There is also a strong -ray flare with peak flux of 5 1 7 ph cm 2 s 1 for ON 246. Since the neutrino event IceCube-17922A arrived when TXS 56+56 was flaring in GeV band, we adopt f accu /f aver as an indicator to qualify the comparison between TXS 56+56 and other Fermi bright blazars. f accu is defined as an accumulated photon flux of the flares, which should have fluxes brighter than the average (f aver ) by a factor of 3. The division of f aver is to eliminate the influence caused by a wide distribution of flux level for different sources. Therefore, f accu /f aver reflects the intensity of GeV activity of a source (a large f accu /f aver arises if a few flares have fluxes much higher than f aver or alternatively there are intense flare activities). If the neutrino is indeed tightly connected to the flare event of blazars, sources with high f accu /f aver value and hard -ray spectra are preferred to be identified as neutrino sources. In the f accu /f aver Γ aver plot, as shown in Fig. 5, TXS 56+56 is distinguished by the high f accu /f aver (i.e., the very strong GeV activities). Together with the declination of TXS 56+56 (as seen in Fig. 6; please note that the IceCube observatory has better sensitivity for sources with declination close to (Aartsen et al. 214)), it may help to explain why TXS 56+56 is the first blazar associated with a significant neutrino excess.

8 We have also investigated the distribution of the accumulated isotropic-equivalent energy of the flare emission (i.e., E accu Γ accu, with proper k correction). In Fig.7 we show the diagram of E accu Γ accu, where is the powerlaw index of the spectrum of the accumulated flare emission. In such a plot, TXS 56+56 is not distinct. Finally we study the change of the hardness of the blazar emission in the flare phase. We define a new parameter Γ = Γ aver Γ accu. In the plot of E accu Γ (see Fig. 8), TXS 56+56 is similar with other blazars, too. Therefore, TXS 56+56 seems to be a normal blazar in many ways except for its very strong GeV activities. Dedicated neurtrino searches in the directions of some bright blazars with strong GeV activities are encouraged, to check whether these sources are also important neutrino sources. 4. SUMMARY The sources of PeV neutrinos and ultra-high-energy cosmic rays are still in debate in the literature. Thanks to the successful performance of IceCube and Fermi-LAT and the quick follow-up observations in optical, radio and X-ray, significant progresses have been made. The most exciting findings are the possible IceCube-17922A/TXS 56+56 association (IceCube collaboration et al. 218) and the the coincidence between a 2 PeV neutrino event and the blazar PKS B1424 418 (Kadler et al. 216), which favors the hypothesis that blazars are important sources of PeV neutrinos and ultra-high-energy cosmic rays. Motivated by the possible connection between the GeV activity and the neutrino IceCube-17922A, in this work we have analyzed the Fermi-LAT data of a group of bright AGNs and identify strong GeV flares. We have compared the properties of these bright AGNs displaying strong GeV activities with TXS 56+56 to look for indication of the ultra-high-energy cosmic ray sources. It turns out that TXS 56+56 appears to be similar to other bright blazars studied in this work (see Fig.1, Fig.2, Fig.7 and Fig.8), except its very strong GeV activities (i.e., the high value of f accu /f aver shown in Fig.5). We suggest to carry out dedicated searches for (weak) neutrino outburst in the historical data of IceCube, as recently done for TXS 56+56 (IceCube collaboration 218), from the directions of these more energetic and harder sources with strong GeV activities. If null results are turned out, new indicator rather than the GeV activity should be identified for the neutrino sources.

9 This research has made use of the NASA/IPAC Extragalactic Database which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research also makes use of the SIMBAD database, operated at CDS, Strasbourg, France. This work was supported in part by NSFC under grants 11525313 (i.e., Funds for Distinguished Young Scholars), 117393 and 113398. H.N.H. is also supported by the Special Postdoctoral Researchers (SPDR) Program in RIKEN. Facilities: Fermi (LAT)

1 REFERENCES Aartsen, M. G., Ackermann, M., Adams, J., et al. 214, Physical Review Letters, 113, 1111 Aartsen, M. G., Abraham, K., Ackermann, M., et al. 215, Physical Review Letters, 115, 8112 Aartsen, M. G., Ackermann, M., Adams, J., et al. 217, Journal of Instrumentation, 12, P312 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 216, Physical Review Letters, 116, 13113 Abdo, A. A., Ackermann, M., Ajello, M., et al. 21, Nature, 463, 919 Acero, F., Ackermann, M., Ajello, M., et al. 215, ApJS, 218, 23 Ackermann, M., Ajello, M., Albert, A., et al. 216, Physical Review Letters, 116, 15115 Ackermann, M., Ajello, M., Atwood, W. B., et al. 215, ApJ, 81, 14 Ahlers, M., & Halzen, F. 215, Reports on Progress in Physics, 78, 12691 Ajello et al., arxiv:172.664 Atoyan, A., & Dermer, C. D. 21, Physical Review Letters, 87, 22112 Becerra. J., et al., 216, ATel No.8722 (http://www.astronomerstelegram.org/?read=8722) Berezinsky, V. S., Blasi, P., & Ptuskin, V. S. 1997, ApJ, 487, 529 Blandford, R. D., & Rees, M. J. 1978, in Pittsburgh Conference on BL Lac Objects, ed. A. M. Wolfe (Pittsburgh, PA: Univ. Pittsburgh Press), 328 B lażejowski, M., Sikora, M., Moderski, R., & Madejski, G. M. 2, ApJ, 545, 17 Böttcher, M. 25, ApJ, 621, 176 Dermer, C. D., & Schlickeiser, R. 1993, ApJ, 416, 458 Fox, D. B., et al. 217, ATel No.1845 (http://www.astronomerstelegram.org/?read=1845) Franckowiak, A., et al., 217, ATel No.1794 (http://www.astronomerstelegram.org/?read=1794) Halzen, F., & Kheirandish, A. 216, ApJ, 831, 12 Kadler, M., Krauß, F., Mannheim, K., et al. 216, Nature Physics, 12, 87 Kopper, C., & Blaufuss, E. 217, GCN Circ. 21916 (https://gcn.gsfc.nasa.gov/gcn3/21916.gcn3) Madejski, G., & Sikora, M. 216, ARA&A, 54, 725 Mannheim, K. 1995, Astroparticle Physics, 3, 295 Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJL, 397, L5 Massaro, E., Giommi, P., Leto, C., et al. 29, A&A, 495, 691 Mattox, J. R., Bertsch, D. L., Chiang, J., et al. 1996, ApJ, 461, 396 Mirzoyan, R. 217, ATel No.1817 (http://www.astronomerstelegram.org/?read=1817) Murase, K., Ahlers, M., & Lacki, B. C. 213, PhRvD, 88, 12131 IceCube Collaboration 213, Science, 342, 1242856 Liao, N. H., Bai, J. M., Liu, H. T., et al. 214, ApJ, 783, 83 Loeb, A., & Waxman, E. 26, JCAP, 5, 3 Padovani, P., Resconi, E., Giommi, P., Arsioli, B., & Chang, Y. L. 216, MNRAS, 457, 3582

11 Paiano, S., Falomo, R., Treves, A., & Scarpa, R. 218, ApJL, 854, L32 Schuster, C., Pohl, M., & Schlickeiser, R. 22, A&A, 382, 829 Sikora, M., Begelman, M. C., & Rees, M. J. 1994, ApJ, 421, 153 Tanaka, Y. T., Buson, S., & Kocevski, D. 217, ATel No.1791 (http://www.astronomerstelegram.org/?read=1791) Tavecchio, F., & Ghisellini, G. 215, MNRAS, 451, 152 The IceCube Collaboration. 218, Science, 361, 147 The IceCube Collaboration et al. 218, Science, 361, eaat1378 Ulrich, M.-H., Maraschi, L., & Urry, C. M. 1997, ARA&A, 35, 445 Waxman, E., & Bahcall, J. 1997, Physical Review Letters, 78, 2292 Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 29, ApJ, 697, 171 Atwood, W., Albert, A., Baldini, L., et al. 213, ArXiv e-prints, arxiv:133.3514 Xu, Z. L. et al., 216, ATel No.991 (http://www.astronomerstelegram.org/?read=991) Aartsen, M. G., Ackermann, M., Adams, J., et al. 214, ApJ, 796, 19

12 2.8 2.6 TXS56 FSRQ BL Lac BZU 2.8 2.6 aver 2.4 2.2 aver 2.4 2.2 2. 2. 1.8 1.8 TXS56 FSRQ BL Lac BZU 1 7 1 6 f aver [ph cm 2 s 1 ] 1 44 1 45 1 46 1 47 1 48 L aver [erg s 1 ] Figure 1. TXS 56+56 in the L aver f aver (left) and L aver Γ aver (right) diagrams. The superscript aver means that they are average quantities for the whole 8-year data. The red circles, green squares and blue triangles represent the FSRQs, BL Lacs and the blazars of unknown type (BZU), respectively. TXS 56+56 is marked as a green pentagram.

13 1 6 TXS56 FSRQ BL Lac BZU f aver [ph cm 2 s 1 ] 1 7 1 44 1 45 1 46 1 47 1 48 L aver [erg s 1 ] Figure 2. TXS 56+56, marked as the green pentagram, in the L aver as the same of Fig 1. f aver diagram. The legends are

14 TXS 56+56 (BL Lac, ISP) 5 TXS 56+56 (BL Lac, ISP) 4 6 35 Photon Flux [1 7 photon cm 2 s 1 ] 5 4 3 2 4 3 2 1 TS Value Photon Flux [1 7 photon cm 2 s 1 ] 8 6 4 2 3 25 2 15 1 TS Value 1 5 555 56 565 57 575 58 MJD [day] 5785 579 5795 58 585 581 5815 582 5825 MJD [day] Figure 3. The weekly -ray light curve from MJD 55317 to MJD 58239 (left) and the 3-day time bin -ray light curve from MJD 57888 to MJD 58239 (right) of TXS 56+56. For any time bin in which the TS value of TXS 56+56 is larger than 25, the photon flux is derived by free the spectral index of TXS 56+56, which is shown as the green dot. The blue dots are for the time bin with 9 < TS < 25, and the fluxes in these bin are derived by keeping spectral indexes of all sources including TXS 56+56 fixed. The red triangles are the 95% upper limits for the time bins with TS values of TXS 56+56 smaller than 9. The gray dashed line is the average flux of TXS 56+56 and the three-time average flux is shown as the black dotted line. The black solid line denotes the arrival time of the neutrino event detected by IceCube.

15 12 3 1 25 8 2 6 15 4 1 2 5 555 56 565 57 MJD [day] 575 58 BL LAC (BL Lac, LSP) 9 8 8 7 6 6 5 4 4 TS Value 35 Photon Flux [1 7 photon cm 2 s 1 ] 14 PKS 2155-34 (BL Lac, HSP) 4 TS Value Photon Flux [1 6 photon cm 2 s 1 ] CTA 12 (FSRQ, LSP) 3 2 2 1 555 56 565 57 MJD [day] 575 58 AO 235+164 (BL Lac, LSP) 14 9 1 8 1 6 4 5 7 8 6 6 5 4 4 3 2 2 2 555 56 565 57 MJD [day] 575 58 PKS 232+17 (FSRQ, ISP) TS Value 8 Photon Flux [1 7 photon cm 2 s 1 ] 1 TS Value Photon Flux [1 7 photon cm 2 s 1 ] 12 15 1 555 56 565 57 MJD [day] 575 58 ON 246 (BL Lac, ISP) 12 16 8 6 7 1 6 8 6 4 4 2 2 555 56 565 57 MJD [day] 575 58 5 6 4 5 4 3 TS Value 8 TS Value Photon Flux [1 7 photon cm 2 s 1 ] 12 Photon Flux [1 7 photon cm 2 s 1 ] 1 14 3 2 2 1 1 555 56 565 57 MJD [day] 575 58 Figure 4. Same as Fig. 3 but for different target sources. The name of the source and its type are marked as the title of each panel.

16 2.8 TXS56 FSRQ BL Lac BZU 2.6 2.4 aver 2.2 2. 1.8 1 1 1 2 1 3 f accu /f aver Figure 5. The f accu /f aver Γ aver diagram, where f accu represents the photon flux accumulated over the time bins with flux at least three times larger than f aver. TXS 56+56 is marked as the green pentagram (same legend as for Fig. 1).

17 8 6 TXS56 FSRQ BL Lac BZU Declination [degree] 4 2 2 4 6 8 1 7 1 6 f aver [ph cm 2 s 1 ] Figure 6. The Declination f aver legend as for Fig. 1). diagram. TXS 56+56 is marked as the green pentagram (same

18 2.8 2.6 2.4 2.2 accu 2. 1.8 1.6 1.4 TXS56 FSRQ BL Lac BZU 1 52 1 53 1 54 1 55 1 56 E accu [erg] Figure 7. The E accu Γ accu diagram, where E accu and Γ accu represent the accumulative equivalent energy and the averaged spectral index for the time bins in which the energy flux is three times larger than f aver, respectively. The legends are as the same of Fig. 1.

19.5 TXS56 FSRQ BL Lac.4 BZU.3.2.1..1 52 53 54 55 56 loge accu [erg] Figure 8. TXS 56+56, marked as the green pentagram, in loge accu Γ diagram (same legend as for Fig. 1).

2 Table 1. Model parameters of all target sources in the sample FL8Y Name Alias Name R.A. a Dec. a Photon Flux Spectral Energy Flux TS Redshift [deg] [deg] [1 8 ph cm 2 s 1 ] Index [1 5 MeV cm 2 s 1 ] Value FL8Y J1.4+2112 TXS 2358+29.385 21.227 8.747 ±.237 2.715 ±.24 2.9 ±.48 3238 1.16 FL8Y J18.6+135 4C +1.2 17.162 1.583 28.651 ±.326 2.369 ±.8 1.163 ±.113 37684 2.99 FL8Y J112.8+328 4C 31.3 18.21 32.138 4.876 ±.183 2.392 ±.25 1.667 ±.5 2526.63 FL8Y J133.1-521 PKS 131 522 23.274 52.1 3.652 ±.155 2.287 ±.26 1.495 ±.53 275.2 FL8Y J21.7 511 28 512 32.693 51.17 1.342 ±.196 2.343 ±.14 3.825 ±.71 1226 1.3 FL8Y J211.2+151 CGRaBS J211+151 32.85 1.86 7.334 ±.21 2.13 ±.15 4.267 ±.121 78.2 FL8Y J221.1+3556 S3 218+35 35.273 35.937 14.846 ±.283 2.28 ±.11 6.157 ±.1 18551.96 FL8Y J222.6+432 3C 66A 35.665 43.35 11.196 ±.224 1.954 ±.1 11.119 ±.26 2458.34 FL8Y J237.9+2848 4C +28.7 39.468 28.82 21.721 ±.198 2.38 ±.7 8.554 ±.97 27262 1.26 FL8Y J238.6+1637 235+164 39.662 16.616 11.963 ±.272 2.175 ±.13 6.214 ±.129 1323.94 FL8Y J246. 465 PKS 244 47 41.5 46.855 5.92 ±.167 2.419 ±.21 1.941 ±.48 4355 1.385 FL8Y J252.8 2219 PKS 25 225 43.2 22.324 8.124 ±.183 2.353 ±.16 2.959 ±.65 795 1.419 FL8Y J33.4 247 PKS 31 243 45.86 24.12 3.78 ±.119 1.879 ±.16 4.948 ±.214 7982.26 FL8Y J324.7+3411 1H 323+342 51.172 34.179 6.654 ±.332 2.871 ±.37 1.426 ±.61 1169.629 FL8Y J339.5 146 PKS 336 1 54.879 1.777 15.226 ±.262 2.346 ±.12 5.611 ±.91 1382.852 FL8Y J348.5 275 PKS 346 27 57.159 27.82 4.332 ±.158 2.261 ±.23 1.865 ±.61 3685.987496 FL8Y J359.6+557 4C +5.11 59.874 5.964 8.821 ±.342 2.593 ±.27 2.351 ±.76 1246 1.51 FL8Y J43.9 365 PKS 42 362 6.974 36.84 18.47 ±.235 2.539 ±.11 5.97 ±.61 2118 1.42284 FL8Y J428.6 3756 PKS 426 38 67.168 37.939 26.287 ±.276 2.9 ±.6 17.24 ±.231 7329 1.111 NONE b PKS 438 43 7.72 43.552 2.221 ±.237 2.626 ±.6.573 ±.44 426 2.852 FL8Y J457. 2324 PKS 454 234 74.263 23.414 24.695 ±.238 2.189 ±.7 12.48 ±.162 4981 1.3 Table 1 continued on next page

21 Table 1 (continued) FL8Y J51.2 158 PKS 458 2 75.33 1.987 11.826 ±.198 2.383 ±.12 4.15 ±.69 8876 2.285998 FL8Y J55.3+459 PKS 52+49 76.347 4.995 15.895 ±.39 2.36 ±.14 5.72 ±.13 11563.954 FL8Y J51.+18 PKS 57+17 77.51 18.12 6.96 ±.36 2.184 ±.19 3.2 ±.118 4128.416 FL8Y J515.5 4556 PKS 514 459 78.939 45.945 3.155 ±.227 2.41 ±.38 1.51 ±.47 1358.194 FL8Y J521.7+2112 VER 521+211 8.442 21.214 9.696 ±.361 1.945 ±.18 9.594 ±.373 12723 2.218 FL8Y J522.9 3628 PKS 521 36 8.742 36.459 12.144 ±.241 2.461 ±.14 3.572 ±.64 1155.56546 FL8Y J53.9+1331 PKS 528+134 82.735 13.532 3.921 ±.321 2.513 ±.35 1.183 ±.64 76 2.6 FL8Y J532. 4827 CRATES J531 4827 82.994 48.46 4.598 ±.242 2.134 ±.18 2.51 ±.94 6265 None c FL8Y J532.7+733 OG 5 83.162 7.545 1.355 ±.337 2.34 ±.16 3.729 ±.9 64 1.254 FL8Y J538.8 445 PKS 537 441 84.71 44.86 15.656 ±.225 2.72 ±.8 11.445 ±.195 3911.894 FL8Y J622.9+3326 B2 619+33 95.718 33.436 5.735 ±.26 2.244 ±.2 2.475 ±.93 3162 None c FL8Y J713.9+1935 MG2 J71354+1934 18.482 19.583 2.252 ±.239 2.595 ±.57.613 ±.32 348.54 FL8Y J721.9+712 716+714 11.43 71.35 22.819 ±.215 2.44 ±.5 14.51 ±.242 9491.3 FL8Y J725.3+1425 4C 14.23 111.32 14.42 4.887 ±.217 2.252 ±.2 2.39 ±.72 3421 1.38 FL8Y J73.3 1141 PKS 727 11 112.58 11.687 17.45 ±.36 2.33 ±.11 6.586 ±.124 12549 1.591 FL8Y J739.2+137 PKS 736+1 114.825 1.618 12.472 ±.331 2.41 ±.14 3.997 ±.81 9398.18941 FL8Y J742.6+5443 BZU J742+5444 115.666 54.74 3.954 ±.181 2.324 ±.23 1.464 ±.55 3314.72 FL8Y J752.+3318 B2 748+33 117.974 33.222 3.561 ±.167 2.412 ±.3 1.183 ±.43 1552 1.935716 FL8Y J88.2 751 PKS 85 7 122.65 7.853 7.687 ±.229 2.34 ±.18 2.443 ±.67 4439 1.837 FL8Y J83.7+241 827+243 127.49 24.22 3.856 ±.27 2.516 ±.32 1.16 ±.53 1747.9414 FL8Y J841.5+753 S5 836+71 13.352 7.895 15.342 ±.253 2.859 ±.14 3.286 ±.46 15186 2.172 FL8Y J852.5 5755 PMN J852 5755 133.161 57.925 3.874 ±.28 2.212 ±.27 1.851 ±.75 1637 None c FL8Y J854.8+26 OJ 287 133.74 2.19 8.843 ±.251 2.214 ±.14 3.943 ±.115 1848.356 FL8Y J94.9 5735 PKS 93 57 136.222 57.585 7.483 ±.249 2.266 ±.19 3.187 ±.86 3829.695 FL8Y J99.1+121 PKS B96+15 137.292 1.36 5.29 ±.32 2.515 ±.3 1.494 ±.65 199 1.2495 FL8Y J921.6+6216 OK +63 14.41 62.264 5.385 ±.192 2.334 ±.19 1.954 ±.61 583 1.446 Table 1 continued on next page

22 Table 1 (continued) FL8Y J948.9+22 PMN J948+22 147.239.374 11.425 ±.274 2.63 ±.18 2.938 ±.61 56.58384 FL8Y J958.7+6533 S4 954+65 149.697 65.565 6.34 ±.14 2.219 ±.14 2.962 ±.68 8226.367 FL8Y J16.7 2159 PKS 14 217 151.693 21.989 4.57 ±.179 2.32 ±.24 1.759 ±.58 2556.331 FL8Y J133.9+65 S4 13+61 158.464 6.852 7.498 ±.234 2.197 ±.16 3.699 ±.9 11691 1.495 FL8Y J148.4+7143 S5 144+71 162.115 71.727 17.398 ±.24 2.246 ±.9 7.715 ±.15 37461 1.15 FL8Y J158.4+133 4C +1.28 164.623 1.566 9.918 ±.266 2.229 ±.15 4.558 ±.14 9457.18516 FL8Y J114.4+3812 Mrk 421 166.114 38.29 21.273 ±.226 1.776 ±.5 42.327 ±.699 117751.31 FL8Y J1153.4+493 115+497 178.352 49.519 3.33 ±.166 2.392 ±.31 1.139 ±.4 247.33364 FL8Y J1159.5+2914 Ton 599 179.883 29.246 18.36 ±.223 2.172 ±.8 9.452 ±.148 3653.72475 FL8Y J1224.9+2122 PKS B1222+216 186.227 21.38 36.327 ±.297 2.39 ±.7 12.475 ±.113 65284.43383 FL8Y J1229.+22 3C 273 187.278 2.52 21.798 ±.364 2.864 ±.16 4.692 ±.65 122.158339 FL8Y J123.2+2517 ON 246 187.559 25.32 5.45 ±.162 2.76 ±.18 3.44 ±.117 668.135 FL8Y J1239.5+443 J123939+4449 189.9 4.7 14.912 ±.273 2.421 ±.13 4.887 ±.77 12355 1.761 FL8Y J1246.7 2547 PKS 1244 255 191.695 25.797 13.74 ±.263 2.323 ±.13 5.264 ±.89 12886.638 FL8Y J1256.1 547 3C 279 194.47 5.789 8.799 ±.428 2.33 ±.4 32.19 ±.29 196753.5362 FL8Y J1312.6+4828 GB6 B131+4844 198.181 48.475 1.938 ±.129 2.261 ±.38.834 ±.41 125.51 FL8Y J1316.1 3338 PKS 1313 333 199.33 33.65 4.87 ±.238 2.344 ±.27 1.799 ±.61 266 1.21 FL8Y J1332. 59 PKS 1329 49 23.19 5.162 7.418 ±.327 2.538 ±.26 2.96 ±.66 315 2.15 FL8Y J1345.5+4453 B3 1343+451 26.388 44.883 17.4 ±.194 2.252 ±.8 7.449 ±.17 34551 2.534 FL8Y J1427.9 426 PKS 1424 41 216.985 42.15 64.429 ±.369 2.172 ±.8 33.733 ±.556 168531 1.522 FL8Y J154.3+129 PKS 152+16 226.14 1.494 24.391 ±.25 2.254 ±.8 1.654 ±.135 4353 1.83928 FL8Y J156.1+3731 B2 154+37 226.54 37.514 6.272 ±.179 2.482 ±.21 1.893 ±.45 4379.674 FL8Y J1512.8 96 151 89 228.17 8.83 79.2 ±.455 2.354 ±.4 28.76 ±.187 15184.36 FL8Y J1517.7 2422 AP Lib 229.42 24.37 9.333 ±.47 2.161 ±.14 5.15 ±.92 8952.49 FL8Y J1522.1+3144 B2 152+31 23.54 31.74 27.55 ±.162 2.459 ±.2 8.57 ±.37 45252 1.4886 FL8Y J1532.7 1319 TXS 153 131 233.16 13.35 8.117 ±.83 2.213 ±.4 3.87 ±.24 619 None c Table 1 continued on next page

23 Table 1 (continued) FL8Y J1555.7+1111 PG 1553+113 238.93 11.19 7.17 ±.72 1.689 ±.5 2.46 ±.316 28725.36 FL8Y J1625.7 2527 PKS 1622 253 246.44 25.46 1.31 ±.14 2.31 ±.3 4.112 ±.42 3581.786 FL8Y J1626. 295 PKS B 1622 297 246.53 29.86 11.32 ±.274 2.549 ±.14 3.161 ±.61 3426.815 FL8Y J1635.2+388 1633+382 248.81 38.13 31.6 ±.334 2.36 ±.7 11.37 ±.117 56527 1.8131 FL8Y J164.5+3945 FGL J1641.4+3939 25.35 39.67 5.175 ±.94 2.441 ±.7 1.65 ±.21 193.5948 FL8Y J1653.8+3945 Mrk 51 253.47 39.76 6.222 ±.123 1.75 ±.9 13.84 ±.412 24831.33 FL8Y J17.+683 GB6 J17+683 255.4 68.5 5.295 ±.83 2.399 ±.7 1.792 ±.2 426.31 FL8Y J179.7+4318 B3 178+433 257.42 43.31 5.185 ±.67 2.35 ±.5 1.898 ±.16 3997 1.27 FL8Y J1733. 134 173 13 263.26 13.8 14.39 ±.164 2.359 ±.4 5.184 ±.39 613.92 FL8Y J1734.3+3858 OT 355 263.59 38.96 5.558 ±.188 2.378 ±.23 1.943 ±.51 3761.975 FL8Y J1748.6+75 S4 1749+7 267.14 7.1 3.986 ±.57 1.951 ±.6 3.999 ±.46 877.77 FL8Y J1751.4+938 OT 81 267.89 9.65 8.73 ±.57 2.264 ±.2 3.727 ±.17 5388.322 FL8Y J18.7+7828 S5 183+78 27.19 78.47 8.541 ±.57 2.24 ±.3 3.839 ±.17 13318.68 FL8Y J181.4+444 S4 18+44 27.38 44.7 4.391 ±.78 2.389 ±.7 1.59 ±.19 2457.663 FL8Y J1829.1 5814 PKS 1824 582 277.3 58.23 1.16 ±.49 2.593 ±.2 2.77 ±.1 4883 1.531 FL8Y J1833.6 213 PKS 183 211 278.42 21.6 45.95 ±.1 2.543 ±.3 12.91 ±.1 27416 2.57 FL8Y J1848.4+3217 CGRaBS J1848+3219 282.9 32.32 5.727 ±.19 2.41 ±.1 1.933 ±.5 1841.8 FL8Y J1848.5+3242 B2 1846+32B 282.14 32.73 7.474 ±.23 2.541 ±.1 2.14 ±.5 2339 None c FL8Y J1849.2+675 CGRaBS J1849+675 282.32 67.9 7.127 ±.57 2.32 ±.3 2.747 ±.14 7111.657 FL8Y J1911.2 26 PKS B198 21 287.79 2.11 1.44 ±.24 2.537 ±.8 2.953 ±.44 3856 1.119 FL8Y J1958. 3845 PKS 1954 388 299.499 38.752 7.18 ±.24 2.369 ±.19 2.522 ±.64 3783.63 FL8Y J2.+658 1ES 1959+65 3. 65.149 6.514 ±.159 1.81 ±.11 11.66 ±.345 19361.47 FL8Y J225.6 735 PKS 223 7 36.419 7.598 14.49 ±.234 2.257 ±.11 6.293 ±.19 15463 1.388 FL8Y J235.4+156 PKS 232+17 38.843 1.935 13.9 ±.318 2.362 ±.14 4.693 ±.91 841.61 FL8Y J2126.3 466 PKS 2123 463 321.628 46.97 3.336 ±.185 2.535 ±.38.946 ±.4 953 1.67 FL8Y J2141.6 641 PMN J2141 6411 325.3 64.26 7.776 ±.197 2.381 ±.18 2.78 ±.63 5834 None c Table 1 continued on next page

24 Table 1 (continued) FL8Y J2143.5+1743 OX 169 325.898 17.73 14.44 ±.261 2.45 ±.14 4.543 ±.75 9425.2174 FL8Y J2147.2 7536 PKS 2142 75 326.83 75.64 13.43 ±.214 2.439 ±.12 4.293 ±.62 1384 1.139 FL8Y J2151.9 327 PKS 2149 36 327.981 3.465 1.27 ±.276 2.896 ±.25 2.169 ±.51 3662 2.345 FL8Y J2158.8 313 PKS 2155 34 329.717 3.226 13.19 ±.194 1.842 ±.8 19.94 ±.445 45844.116 FL8Y J221.8+548 NRAO 676 33.431 5.816 15.26 ±.31 2.647 ±.17 3.864 ±.68 589 1.899 FL8Y J222.7+4216 BL Lac 33.68 42.278 36.25 ±.297 2.161 ±.6 19.47 ±.216 64653.686 FL8Y J2232.5+1143 CTA 12 338.152 11.731 14.3 ±.427 2.255 ±.3 45.47 ±.267 29278 1.37 FL8Y J2236.5 1433 PKS 2233-148 339.142 14.556 8.4 ±.196 2.115 ±.14 5.88 ±.136 9629.325 FL8Y J2244.2+457 TXS 2241+46 341.53 4.954 9.563 ±.232 2.91 ±.12 6.178 ±.151 11734 1.171 NONE b PKS 2247 131 342.498 12.855 15.48 ±.451 2.244 ±.12 6.92 ±.175 15991 None c FL8Y J225.7 286 PMN J225 286 342.685 28.11 3.813 ±.188 2.147 ±.25 2.124 ±.86 298.525 FL8Y J2253.9+168 3C 454.3 343.491 16.148 182.8 ±.554 2.34 ±.3 67.98 ±.273 563863.859 FL8Y J2258. 2759 PKS 2255 282 344.525 27.973 9.355 ±.325 2.422 ±.19 3.62 ±.79 613.92584 FL8Y J2311.+3425 B2 238+34 347.772 34.42 8.535 ±.23 2.352 ±.16 3.113 ±.68 6329 1.817 FL8Y J2323.5 317 PKS 232 35 35.883 3.285 8.531 ±.255 2.234 ±.17 3.88 ±.98 6797 1.41 FL8Y J2329.3 4956 PKS 2326 52 352.337 49.928 23.75 ±.231 2.232 ±.7 1.85 ±.141 4878.518 FL8Y J2345.2 1555 PMN J2345 1555 356.32 15.919 14.29 ±.254 2.12 ±.1 8.955 ±.18 23111.621 Note a: The coordinate (R.A., Dec.) of each source is adopted from the Fermi-LAT monitored source list 7. b: Two sources have no counterparts in the FL8Y catalog which marked as NONE due to they are in a relatively quiescent state during the first eight years. c: Several source marked as None have no redshift values due to the non-detection of optical spectral lines for them. 7 https://fermi.gsfc.nasa.gov/ssc/data/access/lat/msl lc/