Multi-TeV flaring from blazars: Markarian 421 a case study

The TeV blazar Markarian 421 underwent multi-TeV flaring during April 2004 and simultaneously observed in x-ray and TeV energies. It was observed that the TeV outbursts had no counterparts in the lower energies, which implies that this might be an orphan flare. In the context of hadronic model, we have shown that this multi-TeV flaring can be produced due to the interaction of Fermi-accelerated protons of energy $\lesssim 168$ TeV with the background photons in the low energy tail of the synchrotron self-Compton spectrum of the blazar jet. We fit very well the flaring spectrum with this model. Based on this study, we speculate that Mrk 501 and PG 1553+113 are possible candidates for orphan flaring in the future.


Introduction
Active galactic nuclei (AGN) emit electromagnetic radiation the spectrum of which stretches from radio to gamma-rays and exhibit large luminosity variations on time scales ranging from less than an hour up to several years. A super massive black hole is believed to sit at the center of the AGN surrounded by an accretion disk in the inner region and a torus of gas cloud in the outer region. Oppositely directed relativistic jets are ejected from the AGN which are perpendicular to the accretion disk and the torus. In the framework of the unification scheme of AGN, blazars and radio galaxies are intrinsically the same objects, viewed at different angles with respect to the jet axis. When the angle between the jet and the line of sight is small they are called blazars and in contrast, for radio galaxies, the angle between the jet and the light of sight is large. Almost all AGN detected at very high a e-mail: sarira@nucleares.unam.mx b e-mail: luis.miranda@correo.nucleares.unam.mx c e-mail: Subhash.Rajpoot@csulb.edu energy (VHE) (>100 GeV) are blazars with the exception of the three objects, Centaurus A (Cen A) [1,2], M87 and NGC 1275 which are radio galaxies [3,4]. The spectral energy distribution (SED) of these AGN have a double peak structure in the ν−ν F ν plane. In the context of the leptonic model, the low energy peak corresponds to the synchrotron radiation from a population of relativistic electrons in the jet. Although the general consensus is that the high energy peak corresponds to the synchrotron self Compton (SSC) scattering of the high energy electrons with their self-produced synchrotron photons, this result remains inconclusive for various reasons [2]. However, this leptonic model is very successful in explaining the multi wavelength emission from blazars and FR I galaxies [1][2][3][4].
Although the SSC scenario seems to work very well to explain the SED of AGN up to the second peak [5,6], difficulties arise in explaining the multi-TeV emission detected in Cen A [7], flares from the radio galaxy M87 [8], the flares from blazars 1ES 1959 + 650 [9,10] and Markarian 421 (Mrk 421) [11]. Also the inevitable outcome of the leptonic scenario is that emission in multi-TeV energy has to be accompanied by a simultaneously enhanced emission in the synchrotron peak. Unfortunately the enhanced synchrotron emission was not observed in the flaring of 1ES 1959+650 in June 2002 [9] and also probably in the flaring of Mrk 421 in April 2004 [11], which implies that the SSC scenario may not be efficient enough to contribute in the multi-TeV regime.
To explain the orphan flaring of 1ES1959 + 650 a hadronic synchrotron mirror model was proposed by Böttcher [12]. In this model, the high energy protons from the jet interact with the primary synchrotron photons that have been reflected off clouds located at a few pc above the accretion disk. These photons are blue shifted in the jet frame so that there is a substantial decrease in the high energy protons to overcome the threshold for -resonance. Similarly a structured leptonic jet model is proposed to explain the orphan TeV flare [13].
In the hadronic models, the second bump is believed to occur due to the synchrotron emission from the ultra-high energy protons, or emission from the secondary electronpositron pairs or muons produced from the charged pion decay as a result of the interactions of the high energy protons with the background low energy photons. However, the efficiency of the process depends on two factors. One is the strength of the magnetic field in the jet and the other is the proton synchrotron cooling time which has to be smaller than the expansion time scale of the jet. Here the decay of neutral pions to γ -ray pairs can explain the high energy peak in the multi-TeV range [14,15].
Previously it was shown that the above hadronic processes are inefficient to explain the multi-TeV emission. In order to explain the high energy peaks efficient acceleration of relativistic protons to ultra-high energies within the inner part of the jet outflow is required. At the same time the jet kinetic power has to exceed the Eddington luminosity, L Edd , by several orders of magnitude. The hadronic model is employed by Zdziarski et al. [16] to explain the broad band spectra of radio-loud active galactic nuclei and by Cao et al. [17] to explain the TeV spectrum of the blazar 1ES 1101-232. The authors of Refs. [16,17] have shown that super Eddington luminosity in protons is required to explain the high energy peaks. This situation arises because the photon density in the jet is low. However, this problem can be circumvented if one assumes that within the jet there is a dense inner jet region [18]. The multi-TeV peak can also be produced by pp interactions. It was shown earlier that normally the pp process is inefficient to produce γ -rays in the blazar environment unless the jet-cloud interaction is taken into account [19][20][21]. However, the pp interactions are employed to explain the emission of gamma-rays from a population of massive stars surrounding the AGN jets [22,23]. Also emission of high energy gamma-rays by the interaction of the AGN jets with the tidally disrupted atmosphere of Red Giant is studied in [24,25].

Photohadronic model
Sahu et al. [18,26,27] have employed the hadronic model to address the orphan TeV flaring of 1ES 1959 + 650 and the multi-TeV emission from Cen A and M87. Recently this model is also used to make estimates of the neutrinos flux from the TeV-blazars and their spatial correlation with the IceCube events [28]. Here we use the same hadronic model to explain the multi-TeV orphan flaring of Mrk 421. In this model the Fermi-accelerated high energy protons interact with the SSC photons in the core region of the jet (few times the Schwarzschild radius R S ) to produce the -resonance. Subsequently the -resonance decays to charged and neutral pions as follows: The decay of neutral pions to TeV photons gives the multi-TeV SED. Throughout our work we use natural units c = h = 1. The π 0 -decay TeV photon energy E γ and the target SSC photon energy ε γ in the observer frame are related as follows: where D is the Doppler factor and z is the redshift. The observed TeV γ -ray energy and the proton energy E p are related through where is the bulk Lorentz factor of the relativistic jet.
Multi-TeV emission is observed from many blazars and FR I galaxies by Cherenkov telescope arrays. These emissions are classified into non-flaring and flaring events. For the non-flaring events (e.g. emission from Cen A [1]), the injected proton spectrum is a power-law spectrum given by These protons will interact with the background SSC photons having comoving number density n γ (henceforth implies the jet comoving frame) and satisfying the kinematical conditions in Eqs. (2) and (3). As discussed in Ref. [18], the flaring occurs within a compact and confined volume of a smaller inner cone which is enclosed in a bigger outer cone as shown in Fig. 1. Similar internal compact jet models are also proposed to explain the fast variability of the blazars Mrk 501 and PKS 2155-304 [29,30]. In our case, the injected proton spectrum is a powerlaw spectrum supplemented with an exponential decay factor. We can express the injected proton spectrum for both nonflaring and flaring processes in a unified manner as [18,31] where E p,c is the break energy for the high energy protons. Above this break energy the proton spectrum falls off very rapidly. The high energy protons will interact in the flaring region where the comoving photon number density is n γ, f to produce the -resonance. The photon density in the flaring region is much higher than the rest of the blob (non-flaring region with the photon density n γ ) probably due to the copious annihilation of electron-positron pairs, the splitting of photons in the magnetic field, the production of enhanced IC photons in this region and the Poynting flux dominated flow which can form from the magnetic reconnection in the strongly magnetized plasma around the base of the jet [29,32,33]. Here we assume that the spectrum of high energy protons has a exponential decay above the E p,c only for the flaring case. Since the inner jet is more compact than the outer jet, the number density of photons, their energy density and the magnetic field are comparatively higher than the outer jet values. Due to higher photon density in the inner region, the multi-TeV γ -rays and neutrinos can be produced through intermediate -resonance as shown in Eq. (1). The optical depth of the -resonance process in the inner jet region is given by where the resonant cross section is σ ∼ 5 × 10 −28 cm 2 .
The efficiency of the pγ process depends on the physical conditions of the interaction region, such as the size, the distance from the base of the jet, the photon density and their distribution in the region of interest.
In the inner region we compare the dynamical time scale t d = R f with the pγ interaction time scale t pγ = (n γ, f σ K pγ ) −1 to constraint the seed photon density so that multi-TeV photons can be produced. For a moderate efficiency of this process, we can assume t pγ > t d and this gives τ pγ < 2, where the inelasticity parameter is assigned the usual value of K pγ = 0.5. Also by assuming the Eddington luminosity is equally shared by the jet and the counter jet, the luminosity within the inner region for a seed photon energy ε γ will satisfy (4π n γ, f R f ε γ ) L Edd /2. This puts an upper limit on the seed photon density as From Eq. (6) we can estimate the photon density in this region and this is discussed in Sect. 3. For a self consistent treatment, in principle, we should use the photon density n γ, f in the hidden internal jet and solve the coupled transport equations for leptons and photons along the jet axis by taking into account their respective cooling mechanisms as well as the injection spectrum of the primary particles [34]. To avoid this complication we assume the scaling behavior of the photon densities in different background energies as follows: This implies that the ratio of photon densities at two different background energies ε γ 1 and ε γ 2 in the flaring and the nonflaring states remains the same. The photon density in the outer region is calculated from the observed flux in the usual way. In the photohadronic scenario, the number of π 0 -decay photons at a given energy are proportional to both the number of high energy protons and the density of the SSC background photons in the jet, i.e. N (E γ ) ∝ N (E p )n γ . For the flaring case n γ is replaced by the photon density in the flaring region (n γ, f ). The γ -ray flux from the π 0 decay is deduced to be Using the scaling behavior of Eq. (7), the observed multi-TeV photon flux from π 0 -decay at two different observed photon energies E γ 1 and E γ 2 can be expressed as where E γ 1,2 correspond to the proton energy E p 1,2 . In this derivation we have used the relations where E c is the γ -ray cut-off energy corresponding to E p,c . By using the known flux at a particular energy in the flaring/non-flaring state, we can calculate the flux at other energies using Eq. (9). The normalization constant can also be calculated.
In terms of SSC photon energy and its luminosity, the photon number density n γ is expressed as where η is the efficiency of SSC process and κ describes whether the jet is continuous (κ = 0) or discrete (κ = 1). In this work we take η = 1 for 100 % efficiency. The SSC photon luminosity is expressed in terms of the observed flux ( SSC (ε γ ) = ε 2 γ dN γ /dε γ ) and is given by Furthermore, by using Eq. (2), we can simplify the ratio of photon densities given in Eq. (7) to Here we take κ = 0 as representative value for the calculation of n γ . We note that κ = 1 will only lead to the suppression of the photon number density for the outer enveloping jet without changing the main conclusions of our work. This follows from Eq. (12) where proportionality factors depending on κ cancel out. Equation (9) in terms of observed SSC flux and E γ takes the following form: and it uses the SED of SSC photon calculated using the leptonic model. Here the multi-TeV flux is proportional to E −α+3 γ and SSC (ε γ ), while in Eq. (9) it is proportional to E −α+2 γ and to the photon number density n γ (ε γ ). In the photohadronic process ( pγ ), the multi-TeV photon flux is expressed as Both ε γ and E γ satisfy the condition given in Eq.
(2) and the dimensionless constant A γ is given by Equation (14) will be used to calculate the multi-TeV flux from both non-flaring (without exponential decay term) and flaring events from AGN and their subclasses if the emission is due to photohadronic process from the core region. For different blazars/AGN, the value of A γ will be different. This is discussed in the next section. We can calculate the Fermiaccelerated high energy proton flux F p from the TeV γ -ray flux through the relation The optical depth τ pγ is given in Eq. (5). For the observed highest energy proton energy E p , F p (E p ) will always be smaller than the Eddington flux F Edd . This condition puts a lower limit on the optical depth of the process and is given by From the comparison of different time scales and from Eq. (17) we will be able to constraint the seed photon density in the inner jet region.

Mrk 421
Mrk 421 is a high synchrotron peaked BL Lac object (HBL) and was the first extragalactic source with a redshift of z = 0.031 to be established as a TeV emitter [35]. It has a luminosity distance d L of about 129.8 Mpc. Its central supermassive black hole is assumed to have a mass M B H (2−9)×10 8 M corresponding to a Schwarzschild radius of (0.6 − 2.7) × 10 14 cm and the Eddington luminosity L Edd = (2.5 − 11.3) × 10 46 erg s −1 . The synchrotron peak of its SED is in the soft to medium X-ray range and the SSC peak is in the GeV range. It is one of the fastest varying γ -ray sources. In the past, through dedicated multi wavelength observations, the source has been studied intensively. These studies show a correlation between X-rays and very high energy (VHE) γ -rays. A one-zone SSC model explains the observed SED reasonably well [36]. Several major flares were observed in the 2003/2004 season. During April 2004, a large flare took place both in the X-rays and the TeV energy bands. The flare lasted for more than two weeks (from MJD 53,104 to roughly MJD 53,120). But due to a large data gap between MJD 53,093 and 53,104, it is difficult to exactly quantify the duration. The source was observed simultaneously at TeV energies with the Whipple 10 m telescope and at X-ray energies with the Rossi X-ray timing explorer (RXTE). It was also observed simultaneously at radio and optical wavelengths. During the flaring it was observed that the TeV flares had no coincident counterparts at longer wavelengths. Also it was observed that the X-ray flux reached its peak 1.5 days before the TeV flux did during this outburst.
Generally it is believed that the TeV flare might not be a true orphan flare, like the one observed in 1ES 1959+650. On the other hand remarkable similarities between the orphan TeV flare in 1ES 1959+650 and Mrk 421 were observed, including similar variation patterns in the X-ray spectrum. By using the one-zone SSC model, the average SED of Mrk 421 is fitted in Fig. 11 of Ref. [11]. In Fig. 2  Metsähovi radio telescope at Helsinki. The single cross (red) in the optical range is measured by Fred Lawrence Observatory (FLWO) 1.2 m telescope. The points of the first peak in the X-ray range are from Rossi X-ray Explorer (RXTE) [11].
In the one-zone leptonic model, the blob of size R b ∼ 0.7 × 10 16 cm moves down the conical jet with a Lorentz factor 14 and a Doppler factor of D = 14. The emitting region is filled with an isotropic electron population and a randomly oriented magnetic field B = 0.26 G. In the present work we study the flaring of Mrk 421 during April 2004. We use the parameters of the one-zone leptonic model of ref. [11]. The parameters of the one-zone synchrotron model are summarized in Table 1. In principle the Lorentz factor in the inner jet should be larger than the outer jet. But here for simplicity we assume that out in .

Results
The flaring of Mrk 421 in April 2004 was observed in the energy range 0.25 TeV (6.0 × 10 25 Hz) ≤ E γ ≤ 16.85 TeV(4.1 × 10 27 Hz) by the Whipple telescope. In the hadronic model alluded to previously, this corresponds to the Fermi-accelerated proton energy in the range 2.5 TeV ≤ E p ≤ 168 TeV and the corresponding background photon energy will lie in the range 23.6 MeV (5.7 × 10 21 Hz) ≥ Table 1 These parameters (up to B ) are taken from the one-zone synchrotron model of Ref. [11] which are used to fit the SED of Mrk 421. The last three parameters are obtained from the best fit to the observed flare data in our model   Fig. 2. This shows that for orphan flaring, the cutoff energy E c is due to the change from the synchrotron band to the SSC band and can be calculated from their crossover energy region. In Fig. 3 we show the observed multi-TeV SED and the predictions of our model. In our results, the presence of SSC (ε γ ) in Eq. (14) modifies the power-law with the exponential fall-off scenario. From the best fit parameters we obtain the value of the dimensionless constant A γ 20 in Eq. (15). By using the same best fit parameters, we have also calculated the value of A γ from the multi-TeV flare of 1ES 1959+650 and M87 and the multi-TeV emission from Centaurus A, which are determined to be 86, 2 and O(10 −3 ), respectively. We observe that for orphan flaring the condition A γ 1 is satisfied as is seen from 1ES 1959+650 and Mrk421. On the other hand for non-orphan flaring A γ ≤ 1.
In the flaring state, as has been alluded to before, in general, the flux of the two opposing jets can be as high as F Edd /2. However, the highest energy protons with E p = 168 TeV must have flux F p < F Edd /2. Using Eq. (16) this constraint translates into τ pγ > 0.02, which corresponds to n γ, f > 1.3 × 10 10 cm −3 in the inner jet. However, the hidden jet lies between R s and R b . As one representative value we take R f 3 × 10 15 . This value of R f corresponds to a day scale variability, which follows from the fact that the flaring lasted for more than two weeks. From Eq. (6) the seed photon density for ε γ = 0.35 MeV satisfies the inequality n γ, f < 8.9 × 10 10 cm −3 , which translates to the optical depth to be constrained as τ pγ < 0.13. The upshot of all this is that, while the value of τ pγ < 2 is a generic constraint, a more refined present analysis constrains the optical depth to lie in the range 0.02 < τ pγ < 0.13. From this the constraint on the photon number density in the inner jet region is determined to be 1.3 × 10 10 cm −3 < n γ, f < 8.9 × 10 10 cm −3 , which shows that the photon density in this region is high. Due to the adiabatic expansion of the inner blob, the photon density will be reduced to n γ and the energy will dissipate once these photons cross into the bigger outer cone. In spite of the two-zones structure, only the outer zone will be responsible for the observed synchrotron and the IC peaks. From the leptonic model fit to the SED, the magnetic field in the outer jet region is determined to be B = 0.26 G. In general, higher values of the magnetic field are expected to be present in the inner jet region. The maximum proton energy in the hidden jet region will be E p,max ∼ 10 18 (B f /1G) eV and for larger magnetic fields, E p,max can even be higher.
In comparison, in the non-flaring jet scenario case we estimate the proton flux needed to explain the observed TeV γ -rays as follows: Corresponding to the seed photon energy ε γ = 0.35 MeV, the photon density and the optical depth are, respectively, n γ ∼ 4 × 10 3 cm −3 and τ pγ ∼ 1.4×10 −8 . This corresponds to the observed TeV photon energy E γ = 16.8 TeV and the proton energy E p = 168, respectively. The high energy γ -ray flux for E γ = 16.8 is F γ ∼ 3.16 × 10 −11 erg cm −2 s −1 . Using Eq. (16) we obtain the high energy proton flux as F p 10 6 × F Edd . This shows that the normal jet model needs super Eddington power in the protons to explain the high energy peak [16], whereas our inner jet scenario eliminates this extreme energy requirement.
In passing we note that, in the energy range considered, the high energy protons will be accompanied by high energy electrons. These electrons will emit synchrotron photons in the energy range 4×10 19 Hz to 2×10 23 Hz when encountering the magnetic field of the jet. These energy range photons lie in the lower part of the SSC spectrum and thus will not be observed due to their flux being low. The high energy electrons will also emit SSC photons and their energy is given by E IC ∼ γ 2 e ε syn . By considering the electron Lorentz factor in the range 7 × 10 2 ≤ γ e ≤ 4 × 10 4 [36] and the peak energy of the synchrotron photons to be ε syn ∼ 10 18 Hz, we determine the SSC process contribution to the flux in the energy range 2 GeV ≤ E IC ≤ 6.6 TeV. The details of the SSC flux depend on the breaks in SSC spectrum and the spectral index. It is observed that during the flaring of Mrk 421, the X-ray emission reached the peak days after the TeV emission, which poses a serious challenge to the SSC model [11].
The multi-TeV photons in the energy range 0.25 TeV ≤ E γ ≤ 16.8 T eV can interact with the background photons to produce e + e − pairs in which the individual electrons or positrons have energy E γ /2. To produce the e + e − pair, the required threshold seed photon energy of ε γ ≥ 2m 2 e /E γ is needed. During the flaring, multi-TeV γ -rays will interact with soft seed photons in the energy range 0.05 eV ≤ ε γ ≤ 3.5 eV (in between the infrared and the visible range), where σ γ γ ∼ 1.7×10 −25 cm 2 is the maximum cross section. For higher ε γ , σ γ γ will be smaller. We consider the sources of these soft photons to be from the synchrotron emission of 1-10 GeV electrons interacting with the magnetic field B ∼ 1 Gauss and the flux of the ambient photons coming from the disk. On the other hand the multi-TeV γ -rays are produced beyond this region where the photons make the low energy tail (0. 35-23.5 MeV range) of the IC photons. The regions 0.05-3.5 eV and 0.35-23.5 MeV represent two distinct energy ranges. The TeV photons will mostly encounter the tail region of the SSC spectrum. The pair production cross section for ε γ ≥ 0.35 MeV is very small σ γ γ ≤ 10 −30 cm 2 , which corresponds to a mean free path of λ γ γ ≥ 10 19 cm. Hence, the TeV photons will not be attenuated much due to the e + e − pair production. Also it has been observed that, during the flaring of Mrk 421, the variation in the light curves at optical and radio wavelengths are minimal [11], which implies that the low energy photon production was sup-pressed. In addition to the above scenario, the positrons produced from the π + decay have energy E γ /2 and will radiate synchrotron photons in the energy range 2 × 10 17 to 9 × 10 20 Hz. The photon flux F e + ,syn from the synchrotron radiation of e + will be much smaller than F γ (E γ = 0.25 TeV)/8, i.e., F e + ,syn 8 × 10 −11 erg cm −2 s −1 . This flux is well below the observed flux limit in the normal case as can be seen from Fig. 2. From the above analysis our conclusion is that the photon fluxes from the synchrotron emission of the electrons and positrons are not observable during the flaring event of Mrk 421 in April 2004. This situation is very much like the one discussed in the orphan flaring of 1ES 1959+650. Hence we conclude that the flaring of Mrk 421 is orphan in nature. In principle, the flux of multi-TeV γ -rays from the extragalactic sources is attenuated due to interactions with the diffuse extragalactic background light (EBL) through the process γ TeV + γ b → e + e − due to the energy dependent optical depth [37][38][39]. But for low redshifts and the energy range of our interest, the optical depth τ γ γ does not vary much. Hence we can assume an almost constant optical depth [31] so that the spectral shape remains nearly unchanged.
Our study shows that the flaring of blazars can be explained through photohadronic interactions in a compact and confined region within the blazar jet where the photon density is high. In particular, from the study of Mrk 421, we deduce that the orphan flaring can only be possible for those blazars which have a deep valley in between the end of the synchrotron SED and the beginning of the SSC SED as shown in Fig. 2. We note that the HBLs Mrk 501 and PG 1553+113 are possible candidates for orphan flaring in the future. For the orphan flaring events we find A γ 1 and for non-orphan flaring A γ ≤ 1.
If Mrk 501 were to produce an orphan flare then the flare energy will lie in the range 1 TeV ≤ E γ ≤ 8.6 TeV and this corresponds to the background SSC photon energy in the range 4.3 MeV ≥ ε γ ≥ 0.5 MeV. We have estimated this by taking the parameters of Mrk 501 from Ref. [40] as follows: D = 12, z = 0.034. In this case the maximum Fermiaccelerated proton energy will be E p ≤ 10E γ ∼ 86 TeV. Also, another condition which must be satisfied for orphan flaring is The value of SSC (0.5 MeV) will be obtained from the SED of Mrk 501. It is to be noted that, if D changes, accordingly the values of E γ , ε γ and the E p will also change. However, in our scenario the condition of Eq. (18) is independent of these changes.

Conclusions
We conclude that the orphan flaring of Mrk 421 can be explained well by the photohadronic model. We observe that, in this model, the multi-TeV photon flux is proportional to SSC (ε γ ) E −α+3 γ supplemented with an exponential decay term as shown in Eq. (14). This implies that the Fermiaccelerated protons interact with the background photons (in the low energy tail) of the SSC spectrum. In the present work we have shown that the higher photon density in the internal jet scenario eliminates the requirement of the super Eddington luminosity to explain the multi-TeV emission. During the April 2004 flaring of Mrk 421, we have demonstrated that the flux from the synchrotron emission from the high energy e + and e − is suppressed relative to the normal flux, implying that the flaring was orphan in nature, like the one observed in 1ES 1959+650. We have also considered the type of blazar spectrum that will result in orphan flaring and predict that the Mrk 501 and the PG 1553+113 are possible candidates for orphan flaring in the future. Monitoring of these objects by the TeV gamma-ray telescopes will shed more light on the details of the orphan flaring mechanism.
We wish to thank the anonymous referee for helpful comments. The work of S.S. is partially supported by DGAPA-UNAM (Mexico) Project No. IN110815.