Blazar origin of some IceCube events

Recently ANTARES collaboration presented a time dependent analysis to a selected number of flaring blazars to look for upward going muon events produced from the charge current interaction of the muon neutrinos. We use the same list of flaring blazars to look for possible positional correlation with the IceCube neutrino events. In the context of photohadronic model we propose that the neutrinos are produced within the nuclear region of the blazar where Fermi accelerated high energy protons interact with the background synchrotron/SSC photons. Although we found that some objects from the ANTARES list are within the error circles of few IceCube events, the statistical analysis shows that none of these sources have a significant correlation.


I. INTRODUCTION
Interaction of ultra high energy cosmic rays (UHECRs) with the background medium (photons and protons) produce high energy γ-rays and neutrinos. On their way to Earth the UHECRs can be deflected in the magnetic field and the high energy γ-rays can be absorbed. So both of these heavenly messengers will lose their directionality. On the other hand neutrinos will be directly pointing to the source, that is why neutrinos are considered as ideal cosmic messengers.
The IceCube detector located at South Pole in Antarctic ice is precisely built to look for high energy neutrinos (above few TeV) by measuring the Cherenkov radiation of the secondary particles created in each neutrino event. The energy deposited by each event, their direction and topology can be calculated from the trail of the observed Cherenkov light. In 2012 the IceCube Collaboration published two years of data (2010-2012) in which 28 neutrino events with energies between 30 to 1200 TeV were observed [1]. Twenty one of these events are shower-like and the rest are muon tracks. In this analysis two events were PeV neutrino shower events. Adding a third year of analysis in a total 988-days data revealed a total of 37 events, of which 9 are track events and the rest are shower events [2]. The shower events have large angular errors (an average of 15 • ) than the track events (about 1 • ). These events have flavors, directions and energies inconsistent with those expected from the atmospheric muon and neutrino backgrounds. So the study of arrival directions are helpful to find sources of high energy neutrinos and the relevant acceleration mechanism acting within the source.
The isotropic distribution of these IceCube neutrino events suggest contribution from at least some extragalactic sources. There exist different types of potential astrophysical sources to produce UHECRs and hence high energy neutrinos and γ-rays. The list includes: Gamma-ray bursts (GRBs) [3], core of active galactic nuclei (AGN) [4], high energy peaked blazars (HBLs) [5][6][7], starburst galaxies [8] and sources from Galactic center [9]. In Ref. [5] many positional correlations of BL Lac objects and galactic pulsar wind nebulae with the IceCube events are shown. There are also nonstandard physics interpretations of these IceCube events from the decay of superheavy dark matter particles, leptoquark interaction and decay of exotic neutrinos [10] (See [11] for a recent review).
Recently ANTARES collaboration presented a time dependent analysis [12] to look for upward going muon tracks by charge current interaction of ν µ from flaring blazars selected from the Fermi-LAT and TeV γ-ray observed by ground based telescopes H.E.S.S, MAGIC and VERITAS respectively. In this analysis the most significant correlation was found with a GeV flaring blazar from the Fermi-LAT catalog. However, the post-trial probability estimate shows that the event was compatible with background fluctuations. In this work we would like to analyse the above list of Fermi-LAT flaring blazars to see if there is any correlation with the IceCube neutrino events. We use the unbinned maximum likelihood method (MLM) with two different values of the spectral index for our analysis of the positional correlation of these objects.

II. CANDIDATES
Blazars are believed to be the most likely candidates to produce UHECRs and neutrinos [5][6][7]. These are extragalactic objects characterised by relativistic jets with a small viewing angle with respect to the line of sight and are powered by a supermassive black hole in the center of their respective galaxy. These objects are also efficient accelerators of particles through shock or diffusive Fermi acceleration processes with a power-law spectrum given as dN /dE ∝ E −κ , with the power index κ ≥ 2 [13]. Protons can reach ultra high energy through the above acceleration mechanisms. Fractions of these particles escaping from the source can constitute the UHECRs arriving on Earth. These objects also produce high energy γ-rays and neutrinos through pp and/or pγ interactions [14]. The classification for these sources are according to the properties of their emission lines: if a strong broad emission line in the optical spectrum is present, it is classified as Flat Spectrum Radio Quasar (FSRQ), otherwise is a BL Lacerate (BL Lac) object. Depending on the frequency of the first peak, the BL Lac objects are further classified into low (LBL), intermediate (IBL) and high energy (HBL) peaked objects.
The ANTARES collaboration searched for high energy cosmic muon neutrinos using the data taken during the period August 2008 to December 2012. The collaboration selected 41 very bright and variable Fermi-LAT blazars with significant time variability and having the flux > 10 −9 photons cm −2 s −1 for the γ-ray energy above 1 GeV. They have also selected seven TeV flaring objects reported by H.E.S.S., MAGIC and VERITAS telescopes with the expectation that the TeV γ-rays may be correlated with the neutrino events. From the 41 Fermi blazar list, 33 are FSRQs, 7 are BL Lacs and one is unknown. Similarly from the list of 7 TeV flaring blazars one is FSRQ and six are HBLs. It shows that both FSRQs and HBLs are probable sources of very high energy neutrinos and can be possible sources for some of the IceCube event. It is suggested that UHECRs are accelerated in the inner jet of FSRQ and interact with the background from the broad-line region (BLR), synchrotron radiation or the photon from accretion disk [15][16][17].
In a previous article [7] we proposed that photohadronic interactions of the Fermi accelerated high energy protons with the background photons in the nuclear region of the HBLs and AGN are responsible for some of the IceCube events. These objects were observed in multi-TeV γ-rays and some had also flaring. In this model it is assumed that the flaring of blazar in high energy γ-ray occurs within a compact and confined region with a comoving radius R ′ f inside the blob of radius R ′ b [18] (henceforth ′ implies jet comoving frame). In the inner region, the photon density n ′ γ,f is very high compared to the photon density n ′ γ in the outer region i.e. n ′ γ,f ≫ n ′ γ . Fermi accelerated high energy protons undergo photohadronic interaction with the seed photons in the inner region in the self-synchrotron Compton (SSC) regime through the intermediate ∆-resonance. On the other hand, in a normal blazar jet, the photohadronic process is not an efficient mechanism to produce multi-TeV γ-rays and neutrinos because n ′ γ is low, which makes the optical depth τ pγ ≪ 1. But the assumption of compact inner jet region overcome this problem where the optical depth of the ∆-resonance process is τ pγ = n ′ γ,f σ ∆ R ′ f and n ′ γ,f is unknown. We can estimate the photon density in this region by assuming that the Eddington luminosity is equally shared by the jet and the counter jet in the blazar. For a given comoving photon energy ǫ ′ γ in the synchrotron/SSC regime we can get the upper limit on the photon density . Also by comparing the proton energy loss time scale t ′ pγ ≃ (0.5 n ′ γ,f σ ∆ ) −1 and the dynamical time scale t ′ d = R ′ f we can estimate n ′ γ,f , so that the production of multi-TeV γ-rays and neutrinos take place. Not to have over production of neutrinos and γ-rays, we can assume a moderate efficiency (a few percents) by taking τ pγ < 1 which gives n ′ γ,f < (σ ∆ R ′ f ) −1 . In this work we assume 1% energy loss of the UHE protons in the inner region on the dynamical time scale t ′ d corresponding to a optical depth of τ pγ ∼ 0.01 and 15 cm and R ′ f,15 ∼ 1 [7]. In the photohadronic interaction, the intermediate ∆resonance produced will give both high energy neutrinos and γ-rays and relation between the seed photon and the neutrino energy is given by where E ν and ǫ γ are respectively the observed neutrino energy and the background photon energy. The source is located at a redshift z and the bulk Lorentz factor of the jet is Γ. The Doppler factor is given by δ. But for FSRQ and BL Lac objects Γ ≃ δ. So if z and Γ of a blazar are known we can estimate the ǫ γ from the given E ν . The neutrino flux is given as [19] where for all neutrino flavors α (e, µ and τ ), a power-law spectrum of the form is taken. The normalization constant A να is given by (4) where N ν is the number of neutrino events and A ef f,α is the effective area for different neutrino flavors. The energy integrals are done in the limit 25 TeV to 2.2 PeV. The time period T = 988 days is used [1] for the calculation of normalization constant.

III. UNBINNED MAXIMUM LIKELIHOOD METHOD
To identify the possible sources of IceCube events we employ the Unbinned Maximum Likelihood Method (MLM) [20] to find spatial correlation between the blazar sample under consideration and the IceCube events. The signal and the background weights are not separable for an object and both contribute to the likelihood function, which is given by the product of the individual probability densities for the IceCube events as [21] L(n s , where N is the number of IceCube events we take into account, n s /N is the weight associate with the signal probability density function (PDF) and its values vary between 0 and 1. The background PDF depends on the the neutrino energy and the declination which is expressed as The background is constructed from the integrated effective areas of the IceCube 79 strings configuration [22]. The neutrino effective area depends on the detector geometry and the absorption of the neutrinos by the Earth. The background PDF takes into account the contribution from the atmospheric muon neutrinos. Above ∼ 100 TeV, neutrinos from the decay of charm hadrons D ± , D 0 contribute to the background neutrino flux known as prompt flux. Equal number of neutrinos and antineutrinos of electron and muon flavors are produced in this process. However, the prompt flux is poorly understood in the high energy limit. For the background calculation we also include the contribution from the prompt background [23,24]. The signal PDF is defined as the product of a spatial term and the energy term as shown below where we have defined which is a Gaussian function [25]. In the above Eq.(8), |x i − x s | 2 is the space angle difference between the source and the reconstructed event direction and σ i is the standard deviation of the i th IceCube angular error distribution. We also define The value of δχ 2 ≤ 1 signifies that the object is inside the median angular error ∆ of the IceCube event. The signal energy PDF E i depends on the event energy, spectral index κ and the declination. Here we use κ = 2 and 2.5 for our analysis. The ANTARES analysis takes into account both the temporal and energy dependence of the flaring events whereas our analysis is independent of the time. The observed IceCube events can be modelled by taking into account two hypothesis: (1) the events could be produced by atmospheric muons and the muon neutrinos (background), or (2) from an astrophysical source which also includes the background contribution. A good test of compatibility is the ratio of these two hypothesis. We can take the ratio of the likelihood with the background of unique weight (n s = 0) and the maximized likelihood of the second hypothesis with the corresponding n s values defined as n s = n * s . Now to evaluate each point source .
For this procedure we use a full-sky IceCube events. For our present analysis, we take into account 36 events out of reported 37 events (event 32 is excluded in the present analysis because its energy and direction are not reported). We calculate the significance of each source location, running 10,000 simulations in which the declination of each IceCube sample event is fixed but the right ascension is randomized. The p-value is calculated as the number of simulations with T S (sim) ≥ T S divided by the total number of simulations for a given source, where T S (sim) is the TS value obtained from the simulation. Also, the posteriori p-value for each object is estimated as the fraction of the randomized simulations that yields an equal or higher TS value for at least one of the 41 ANTARES sources. The compatibility of the second hypothesis depends on the estimate of the posteriori p-value. If the posteriori p-value is close to unity then it is consistent with the background. Below each object we also put their coordinates, Right Ascension and Declination ( R.A., Dec.) in degrees (this table is given in equatorial coordinates). The second column gives the type of object and below this we also give its redshift (z) and the bulk Lorentz factor (Γ). In the fourth column, the δχ 2 of the object is given. In the fifth and the sixth columns the deposited neutrino energy Eν /T eV and the corresponding seed photon energy ǫγ /keV are given. In columns seventh and eighth the values of the n * s and TS are given from the Maximum Likelihood Method. In columns ninth and tenth the p-value and the posteriori p-value (post p-value) are also shown. The last three objects are without (upper value) and with (lower value) the prompt contribution to the background PDF.

IV. RESULTS
In the context of recent IceCube results, we analysed the 41 flaring blazars taken from the Fermi-LAT catalog which are previously studied by the ANTARES collaboration to look for possible temporal and spatial correlation [12]. We have also analysed the 7 TeV flaring objects as discussed by ANTARES collaboration for the possible spatial correlation with the IceCube events. In fact all these 7 objects are there in the TeVCat[26] which we had already analysed in Ref. [7] and found that the only HBL, PG 1553+113 has the positional correlation with the IceCube event 17. So we don't discuss about these 7 flaring objects here any more. For our analysis of the possible correlation of IceCube events with the ANTARES sources we use the unbinned MLM and two different values of spectral index κ = 2 and 2.5. We also do the separate analysis with and without the contribution from the prompt flux coming from the charm hadron decay. Our results are summarised in Table I. All the 28 shower events with their individual errors and the 8 track events are shown in the sky map with equatorial coordinates in Fig. 1. The positions of ten FSRQs and two BL Lac objects are also shown in the sky map.

A.
Spectral index κ = 2 From the 41 Fermi blazars of ANTARES list, 32 objects have TS > 0 for the spectral index κ = 2 without the prompt contribution to the background. However, this number reduces to 19 when we include the charm contribution.
From the above 32 objects 12 are within the median angular error of at least one IceCube event having δχ 2 < 1. The FSRQ, PKS 0235-618 is the only object associated with two IceCube events (7,20). The FSRQs, PKS 2326-502, PKS 0208-512, PKS 0235-618 and PKS 0244-470 are within the error circle of event 7, while the FS-RQs, 3C454.3, B22308+34, CTA102 and PKS 2227-08 are within the error circle of event 34. Another two FS-RQs, PMNJ2345-1555 and PMNJ 2331-2148 are within the error circle of the IceCube event 21. The BL Lac objects, OJ287 and PKS0805-07 are coincident with the events 26 and 27 respectively. All the relevant parameters of the above objects are shown in Table I.
The posteriori p-value of all the above 12 objects are ≥ 99%. This shows that our result (without the prompt contribution to the atmospheric background) is consistent with the background fluctuation.
By including the prompt contribution to the background we found that 19 objects have TS > 0 of which only three objects two BL Lac objects (OJ287, PKS 0805-07) and one FSRQ (3C454.3) are within the median angular error of three IceCube events (26, 27, 34). These three objects are shown in the table.
We observed that the background photon energy ǫ γ for most of the events are below < 40 keV which shows that the photon density n ′ γ,f can be large in the inner region of the jet. By assuming a conservative 1% energy loss by the UHE protons we get the photon density in the inner region n ′ γ,f ∼ 2 × 10 10 cm −3 which has a radius R ′ f ∼ 10 15 cm. Estimate of R ′ f value depends on the outer blob radius R ′ b , while the later parameter is adjusted to fit the spectral energy distribution (SED) in the leptonic model of the objects. However, for most of the objects R ′ b > 10 15 cm is taken to fit the SED [7]. So, here we take R ′ f ∼ 10 15 cm for the estimation of n ′ γ,f . The simulation shows that the 0 < TS < 1 for all the objects.
The diffuse neutrino flux F ν for all these objects is 2.31 × 10 −9 GeV cm −2 s −1 sr −1 because for κ = 2, the integral in Eq. (2) is independent of the redshift, so the F ν is the same for all the objects. On the other hand it varies slightly for κ = 2. The high posteriori p-value for all these objects shows that our result is consistent with the background fluctuation. We also repeated the simulation for κ = 2.5. As κ changes the flux reduces and found that non of the 41 objects satisfy the condition TS > 0.

V. DISCUSSION
ANTARES collaboration looked for possible temporal and spatial correlation of 41 flaring objects selected from the Fermi-LAT catalog. We analysed the same objects for the possible spatial correlation with the IceCube events. For our analysis, we take into account the energy dependence of both the background and the signal constructed from the data of the 79 IceCube string configuration. We consider two different values of the spectral index 2 and 2.5 and also analyse our results with and without the prompt contribution to the atmospheric neutrino flux. We observed that, from the 41 flaring objects, for κ = 2, the MLM gives twelve objects (without prompt flux contribution) and three objects (with prompt flux contribution) within the error circle of some IceCube events. For these objects we have also estimated the neutrino flux. However, for all these possible candidates, the TS value is very small which leads to very high posterior p-values ≥ 99% and is consistent with the background fluctuation. It is possible that the high energy neutrino flux from these objects are much below the IceCube limit or blazars may not have powerful central engine to produce very high energy cosmic rays. So most of the events in IceCube can be from some other type of sources. We have to wait for more data to look for possible correlation of FSRQs and BL Lac objects with the IceCube events.