1. Introduction
Blazars are a subclass of Active Galactic Nuclei (AGN) and the dominant extra galactic population in gamma-rays [1]. These objects show rapid variability in the entire electromagnetic spectrum and have non thermal spectra which implies that the observed photons originate within the highly relativistic jets oriented very close to the observers line of sight [2]. Due to the small viewing angle of the jet, it is possible to observe the strong relativistic effects, such as the boosting of the emitted power and a shortening of the characteristic time scales, as short as minutes [3,4]. Thus, these objects are important to study the energy extraction mechanisms from the central super-massive black hole, physical properties of the astrophysical jets, acceleration mechanisms of the charged particles in the jet and production of ultra high energy cosmic rays, very high energy γ-rays and neutrinos.
The spectral energy distribution (SED) of these blazars has a double peak structure
in the v - vFv plane.
The low energy peak corresponds to the synchrotron radiation from a population of
relativistic electrons in the jet and the high energy peak believed to be due to the
synchrotron self Compton (SSC) scattering of the high energy electrons with their
self-produced synchrotron photons [5,6]. Depending mostly on the optical spectra,
blazars can be divied into BL Lacertae objects (BL Lacs) and flat- spectrum radio
quasars (FSRQ) [7]. Based on the
location of the first peak, BL Lacs can be further classified into low energy peaked
blazars (LBLs) (
Flaring seems to be the major activity of the blazars which is unpredictable and switches between quiescent and active states involving different time scales and fluxes. While in some blazars a strong temporal correlation between X-ray and multi-TeV γ-ray has been observed, outbursts in some other have no low energy counterparts (orphan flaring) [15,16] and explanation of such extreme activity needs to be addressed through different mechanisms. It is also very important to have simultaneous multiwavelength observations of the flaring period to constrain different theoretical models of emission in different energy regimes.
Different theoretical models have been proposed to explain the flaring from AGN and its subclasses. These models are mainly clssified into two categories: the leptonic models and the hadronic models. In the leptonic model scenario, the high energy electrons upscatter the low energy photons in the jet through the SSC process and the photons can be in the GeV-TeV region. This model has severe limitations to explain VHE γ-rays and the orphan flaring observed in many blazars. The multi-zone leptonic model can explain these high energy emissions, however, one has to increase the number of parameters in the model. In the hadronic synchrotron-proton blazar model [19,21], emission of synchrotron photons from protons take place. In this scenario, the Fermi accelerated protons in the jet magnetic field, emit synchrotron radiation which will be suppressed by a factor of m-4p, where mp is the proton mass. So ultra high energy proton flux is needed to explain the VHE gamma-rays. It also needs a strong magnetic field for the synchrotron process to be effective, but a strong magnetic field in the jet is not very usual. In the jet-in-jet model of Giannios et al. [22] minijets are formed within the jet due to flow instabilities and these minijets move relativistically with respect to the main jet flow. The interaction of the daughter jets with the main jet are responsible for the production of VHE gamma rays. While the minijets are aligned with our line of sight, the VHE gamma rays are beamed with large Doppler factor. This scenario can explain the 2010 flare of the radio galaxy M87 but does not provide a quantitive prediction of the light curve of the flare. The lepto-hadronic model [23] fits to the low energy γ-ray spectrum by Large Area Telescope (LAT) and High-Energy Stereoscopic System (H.E.S.S.) low state but not the flaring state. Similarly, the magnetosphere model [24-26] can explain the hard TeV spectrum but in this case also there is no detailed quantitive predication for the VHE light curve. Also interaction of Fermi acclerated protons with the MeV photons emitted by the Wein fireball in the base of the jet can explain the orphan falring from 1ES 1959+650 and Mrk 421 [27]. It is well known that EBL plays an inportant role for the attenuation of γ-rays even from the nearby blazars and this model does not take into account the EBL effect. The VHE SED of blazars are modelled using different hadronic models [28-31]. Using the hadronic models, the high energy cosmic ray and the neutrino fluxes are also estimated.
In a series of papers Sahu et al. [35-42] have explained the GeV-TeV flaring from many blazars using photohadronic scenario. In this scenario, Fermi accelerated protons interact with the background photons in the jet environment to produce Δ-resonance. Subsequent decay of the Δ-resonance to pions can produce VHE photons and neutrinos. The produced gamma-rays can explain very well the observed spectra from the flaring blazars.
The TeV photons of the flare can interact with the background soft photons in the jet to produce e+e- pairs. However, production of the lepton pair within the jet depends on the size of the emitting region and the photon density in it. Also the required target soft photon threshold energy ∈γ ≥ 2m2e/Eγ is needed. It has been observed that the jet medium is transparent to pair production where the optical depth is very small [35,43]. Also the TeV photons on their way to Earth can interact with the extragalactic background light (EBL) to produce the lepton pair. TeV photons from the sources in the cosmologically local Universe (low redshift sources) are believed to propagate unimpeded by the EBL, although the effect is found to be non negligible [43].
In this work the goal is to use the photohadronic model of Sahu et al. [35] and different template EBL models [47,48] to explain the observed GeV-TeV flares from different HBLs. In this article, I review the photohadronic model of Sahu et al. and discuss about its applicability to explain the flaring events from high energy blazars. The plan of the papers is as follows. In Sec. 2 a short account of different EBL models is given. In Sec. 3 a detail discussion about the photohadronic model and the kinematical condition for the process ργ → Δ and its subsequent decay is given. Also the relation between the observed GeV-TeV gamma-ray flux and the gamma-ray from the photohadronic process is shown. Sec. 4 is dedicated to discuss the results of different flaring blazars and compare the results with other models. In Sec. 5 a general remark is given regarding the problems of different models. Finally I summarize the photohadronic model and discuss about the future research in this field in Sec. 6.
2. EBL Models
A distinctive feature of the high energy γ-ray astronomy is that the high energy gamma-rays undergo energy dependent attenuation en route to Earth by the intervening EBL through electron-positron pair production [44]. This interaction process not only attenuates the absolute flux but also significantly changes the spectral shape of the VHE photons. The diffuse EBL contains a record of the star formation history of the Universe. A proper understanding of the EBL SED is very important for the correct interpretation of the deabsorbed VHE spectrum from the source. The direct measurement of the EBL is very difficult with high uncertainties mainly due to the contribution of zodiacal light [45,49], and galaxy counts result in a lower limit since the number of unresolved sources (faint galaxies) is unknown [50].
Keeping in mind the observational constraints and the uncertainty associated with the direct detection of the EBL contribution, several approaches with different degrees of complexity have been developed to calculate the EBL density as a function of energy for different redshifts. A wide range of models have been developed to model the EBL SED based on our knowledge of galaxy and star formation rate and at the same time incorporating the observational inputs [46-48,51-54]. Mainly three types of EBL models exist: backward and forward evolution models and semi-analytical galaxy formation models with a combination of information about galaxy evolution and observed properties of galaxy spectra. In the backward evolution scenarios [52], one starts from the observed properties of galaxies in the local universe and evolve them from cosmological initial conditions or extrapolating backward in time using parametric models of the evolution of galaxies. This extrapolation induces uncertainties in the properties of the EBL which increases at high redshifts. However, the forward evolution models [47,51] predict the temporal evolution of galaxies forward in time starting from the cosmological initial conditions. Although, these models are successful in reproducing the general characteristics of the observed EBL, cannot account for the detailed evolution of important quantities such as the metallicity and dust content, which can significantly affect the shape of the EBL. Finally, semi-analytical models have been developed which follow the formation of large scale structures driven by cold dark matter in the universe by using the cosmological parameters from observations. This method also accounts for the merging of the dark matter halos and the emergence of galaxies which form as baryonic matter falls into the potential wells of these halos. Such models are successful in reproducing observed properties of galaxies from local universe up to z∼6.
3. Photohadronic Scenario
In general, in the leptonic one-zone synchrotron and SSC jet model the emitting region is a blob with comoving radius R'b (where ' implies the jet comoving frame and without prime are in the observer frame), moving with a velocity βc corresponding to a bulk Lorentz factor Γ and seen at an angle θob by an observer which results with a Doppler factor D = Γ-1 (1 - βc cosθ ob )-1 [10,15]. The emitting region is filled with an isotropic electron population and a randomly oriented magnetic field B’. The electrons have a power-law spectrum. The energy spectrum of the Fermi-accelerated protons in the blazar jet is also assumed to be of power-law. Due to high radiative losses, electron acceleration is limited. On the other hand, protons and heavy nuclei can reach Ultra-High Energy (UHE) through the same acceleration mechanism.
The photohadronic model discussed here relies on the above standard interpretation of the leptonic model to explain both low and high energy peaks by synchrotron and SSC photons respectively as in the case of any other AGN. Thereafter, it is assumed that the flaring occurs within a compact and confined volume of size R'f inside the blob of radius R'b. The geometrical description of the jet structure during a flare is shown in Fig. 1. In this scenario the internal and the external jets are moving with almost the same bulk Lorentz factor Γin ≃ Γext ≃Γ and the Doppler factor D as the blob (for blazars Γ≃D). Within the confined volume, the injected spectrum of the Fermi accelerated charged particles have a power-law, and for the protons with energy Ep it is given as
where the spectral index α ≤ 2 [5]. Also in this small volume, the comoving photon number density n'γ,f (flaring) is much higher than rest of the blob n'γ (non-flaring).
The dominant mechanism through which the high energy protons interact with the background photons in the inner jet region is given by
which has a cross section σΔ ∼ 5 x 10-28 cm2.
Subsequently, the charged and neutral pions will decay through
i.e. the ratio of photon densities at two different background energies ∈γ1 and ∈γ2 in the flaring (n'γ,f) and in the non-flaring (n'γ) states remains almost the same. The photon density in the outer region is known to us from the observed flux. So by using the above relation, we can express the unknown inner photon density in terms of the known outer density which can be calculated from the observed flux in the usual way using the observed/fitted spectral energy distribution (SED). Henceforth, for our calculation, we shall use n’γ and its corresponding flux rather than the one from the inner jet region which is not known.
3.1 Kinematical condition
For the above process in Eq. (2) to take place, the center-of-mass energy of the interaction has to exceed the Δ-mass 1.232 GeV which corresponds to the kinematical condition
where E'p and ∈'γ are the proton and the background photon energies in the comoving frame of the jet, respectively. Also for high energy protons we take βp ≃ 1. Since in the comoving frame the protons collide with the SSC photons from all directions, in our calculation we consider an average value (1-cosθ) ∼ 1 (θ in the range of 0 and π). In the observer frame, one can re-write the matching condition as
Here
is the observed background photon energy, while
is the energy of the proton as measured by the observer on Earth, if it could escape the source and reach earth without energy loss and z is the redshift of the object.
In the comoving frame, each pion carries ∼ 0.2 of the proton energy. Considering that each π0 decays into two γ-rays, the π0 -decay γ-ray energy in the observer frame (Eγ) can be written as
The matching condition between the π0 -decay photon energy Eγ and the target photon energy ∈γ is therefore
So from the known flare energy Eγ of a blazar the seed photon energy ∈γ can be calculated when the Doppler factor and the bulk Lorentz factor is known from the leptonic model fit to the blazar SED.
3.2 Flux calculation
The observed VHE γ-ray flux depends on the background seed photon density and the
differential power-spectrum of the Fermi accelerated protons given as
The optical depth of the Δ-resonance process in the inner jet region is given by
The efficiency of the ργ 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 ργ interaction
time scale
From Eq. (11), we can estimate the photon density in this region. 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 (
Using the Eqs. (12) and (13) we can simplify the ratio of photon densities given in Eq. (3) to
The γ-ray flux from the π0 decay is deduced to be
The EBL effect attenuates the VHE flux by a factor of e-τγγ, where τγγ is the optical depth which depends on the energy of the propagating VHE γ-ray and the redshift z of the source.
Including the EBL effect, the relation between observed flux Fγ and the intrinsic flux Fint is given as
Then the EBL corrected observed multi-TeV photon flux from π0 -decay at two different observed photon energies Eγ1 and Eγ2 can be expressed as
where we have used
The ΦSSC at different energies are
calculated using the leptonic model. Here, the multi-TeV flux is proportional to
Both ∈γ and Eγ satisfy the condition given in Eq. (9) and the dimensionless constant Aγ is given by
Comparing Eqs. (16) and (19), the intrinsic flux Fint is given as
Using Eq. (19), we can calculate the EBL corrected multi-TeV flux where Aγ can be fixed from observed flare data. We can calculate the Fermi accelerated high energy proton flux Fp from the TeV γ-ray flux through the relation [36]
The optical depth τpγ is given in Eq. (10). For the observed highest energy γ-ray, Eγ corresponding to a proton energy Ep , the proton flux Fp (Ep) will be always smaller than the Eddington flux FEdd. This condition puts a lower limit on the optical depth of the process and is given by
From the comparison of different times scales and from Eq. (23) we will be able to constrain the seed photon density in the inner jet region.
It is observed that, for the observed flare energy Eγ
the range of ∈γ is always in the
low energy tail region of the SSC band and the corresponding SSC flux in this
range of seed photon energy is exactly a power-law given by
ΦSSC
From the leptonic model fit to the observed multiwavelength data (up to second peak) during a quiescent/flaring state we can get the SED for the SSC region from which Φ0 and β can be obtained easily. By expressing the observed flux Fγ in terms of the intrinsic flux Fγ,in and the EBL correction as
where the intrinsic flux is given in Ref. [40] as,
where Aγ is a dimensionless normalization constant and can be fixed by fitting the observed VHE data. As discussed above the power index β is fixed from the tail region of the SSC SED for a given leptonic model which fits the low energy data well. So the Fermi accelerated proton spectral index α is the only free parameter to fit the intrinsic spectrum.
4. Results
From the continuous monitoring and dedicated multi-wavelength observations of the nearest HBLs Markarian 421 (Mrk 421, z=0.0308), Mrk 501 (z=0.033 [58,59]) and 1ES 1959+650 (z=0.047 [60]), several major multi-TeV flares have been observed [61-65]. Strong temporal correlation in different wavebands, particularly in X-rays and VHE γ-rays has been observed in some flaring events. However, in some other flaring events no such correlation is observed [15,16], which seems unusual for a leptonic origin [9,10,12,66] of the multi-TeV emissions and needs to be addressed through other alternative mechanisms [55,67-71,98]. Below we shall discuss the flaring of Mrk 421, Mrk 501 and 1ES1959+650 separately.
4.1 Markarian 421
Mrk 421 was the first extragalactic source detected in the multi-TeV domain [72] and also it is one of the fastest varying γ-ray sources. It has a luminosity distance d L of about 129.8 Mpc and its central supermassive black hole is assumed to have a mass MBH ≃ (2-9)x108 Mʘ corresponding to a Schwarzschild radius of (0.6-2.7) x 1014 cm and the Eddington luminosity LEdd = (2.5 - 11.3) x 1046 erg s-1. The synchrotron peak (1st peak) of its SED is in the soft to medium X-ray range and the SSC peak (2nd peak) is in the GeV range. Through dedicated multi wavelength observations, the source has been studied intensively. These studies show a correlation between X-rays and VHE γ-rays. A one-zone SSC model explains the observed SED reasonably well [73]. Several large flares were observed in 2000 - 2001 [74-76] and 2003 - 2004 [16,77]. 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 Modified Julian Date (MJD) 53,104 to roughly MJD 53,120). 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) [16]. It was also observed simultaneously at lower wavelengths (both radio and optical). 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. 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. A strong outburst in multi-TeV energy in Mrk 421 was first detected by Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory on 16th of February 2010 and follow up observations were done by the HESS telescopes during four subsequent nights [61].
A six month long multi-instrument campaign by the Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) telescopes observed VHE flaring from Mrk 421 on 25th of April 2014 and the flux (above 300 GeV) was about 16 times brighter than the usual one. This triggered a joint ToO program by X-ray Multi-mirror Mission - Newton (XMM-Newton), VERITAS, and MAGIC instruments. These three instruments individually observed approximately 3 h each day on April 29, May 1, and May 3 of 2014 [78]. The simultaneous VERITAS-XMM-Newton observation is published recently and it is shown that the observed multiwavelength spectra are consistent with one-zone synchrotron self-Compton model [78].
The observed multi-TeV flux from the above flaring events are EBL corrected and for this correction well known EBL models of Dominguez et al. [46] (EBL-D) and Inoue et al. [48] (EBL-I) are used. In Fig. 1 the attenuation factor for these two models as functions of observed gamma-ray energy E γ are shown and both are practically the same in all the energy ranges (there is a minor difference in the energy range 600 Gev ≤ E γ ≤ 1TeV.
4.1.1 The flare of April 2004
The multi-TeV flare of April 2004 was the first flare observed in multiwavelength by the Whipple telescope in the energy range 0.25TeV(6.0 x 1025 Hz) ≤ Eγ ≤ 16.85TeV(4.1 x 1027 Hz) and it was difficult to explain by one-zone leptonic model [16]. As discussed above, the photohadronic interpretation of the flare data needs the leptonic model as input and here the one-zone leptonic model of Ref. [16] (lep-1) is used. In this model the bulk Lorentz is Γ = D = 14. The above range of Eγ corresponds to the Fermi accelerated proton energy in the range 2.5 TeV ≤ Ep ≤ Ep ≤ 168 TeV and the corresponding background photon energy is in the range 23.6MeV(5.7 x 1021Hz) ≥ ϵγ ≥ 0.35MeV(8.4 x 1019Hz). This range of ∈γ is in the low energy tail region of the SSC SED and its flux is expressed as power-law given in Eq. (24) with Φ0 = 6.0 x 10-10 erg cm-2s-1 and β = 0.48 which is shown in Fig. 2.
Very good fit to the multi-TeV flare data is obtained for α = 2.7, and the
normalization constant A
γ for EBL-D is 2.3 and for EBL-I is 2.5 respectively which are
shown in Fig. 3. Below 4 TeV both these
EBL models fit the data very well, above this energy there is a slight
difference due to the change in the attenuation factor. It is observed that
EBL-D fits better than the EBL-I. To compare the photohadronic model fit
without EBL correction but with an exponential cut-off
(dN/dEγ α
E-αγ exp
(-Eγ /
Ec)) [38], and the multi-zone leptonic model fit
(magenta curve) [16] are
also shown in the same figure (red curve). It can be seen from Fig. 3 that the multi-zone fit is not so
good compared to other fits for Eγ ≤ 15 TeV.
However, for higher energy it has the same behavior as EBL-D and EBL-I. With
the exponential cut-off scenario, a good fit is obtained for the spectral
index α = 2.7 and the cut-off energy
EC = 6.2 TeV. Again
comparing this with the EBL corrected models, below 4 TeV, all these fits
are exactly the same. However, above 4 TeV we observe some discrepancy among
these fits and above 10 TeV the fits of EBL-D and EBL-I fall faster than the
exponential cut-off scenario. Comparison of EBL-I with the exponential
cut-off scenario shows that, for Eγ ≤ 10 TeV
both e-Eγ/Ec and
e-τγγ are almost the same
and above 10 TeV the attenuation factor falls faster than the exponential
cut-off. The intrinsic flux in EBL-D and EBL-I are almost the same and
having power-law behavior with
By comparing different time scales i.e. expansion time
scale, interaction time scale of ργ interaction and using the fact that the
high energy proton luminosity to be smaller than the Eddington luminosity in
the inner jet region of size
R'f ≃ 3x1015 m,
the range of optical depth for the Δ-resonance production is estimated as
0.02 < πργ < 0.13. This corresponds to a photon density in
the inner jet region as
4.1.2 The flare of February 2010
On 16th of February 2010, a strong outburst in multi-TeV gamma-rays from Mrk 421 was observed by VERITAS telescopes and follow up observations were carried out by HESS telescopes from 17th to 20th of February a total of 6.5 h. These data were taken in 11 runs with each run ∼ 28 minutes duration [61]. The HESS telescopes observed the flare in the energy range 1.67 TeV (4.0 x 1026Hz ≤ E γ ≤ 20.95 TeV (5.0 x 1027Hz). During this period there was no observation of multiwavelength SED, particularly in the SSC band. So to interpret the observed flare data, the observed SED at an earlier and later epochs were used. The first one is the lep-1, which is used for the interpretation of the April 2004 flare and the second leptonic SED is from the multiwavelength observation of Mrk 421 during January to March 2013, undertaken by GLAST-AGILE Support Program of the Whole Earth Blazar Telescope (GASp-WEBT) Swift, NuSTAR Fermi-LAT, MAGIC, VERITAS [79] and fitted with one-zone leptonic model where the bulk Lorentz factor Γ = 25 is used (lep-2). The parameters of lep-1 and lep-2 are shown in Table I.
lep-1 | lep-2 | |
D (Doppler factor) | 14 | 25 |
R’b (Blob radius) | 0.7 x 1016 cm | 0.9 x 1016 cm |
R’b (Inner blob radius) | ≈ 3 x 1015 cm | ≈ 3 x 1015 cm |
B’ (Magnetic field) | 0.26 G | 0.17 G |
EBL Model | α, Aγ | α, Aγ |
EBL-I | 2.7,2.3 | 2.6,2.4 |
EBL-D | 2.7,2.5 | 2.6,2.5 |
Again, the above range of Eγ corresponds to the proton energy in the range 16.7 TeV ≤ Ep ≤ 210 TeV. Using lep-1, where Γ=D=14, the seed photon energy lies in the range 0.28MeV(6.8x1019 Hz) ≤ ∈γ ≤ 3.53 MeV(8.5x1020 Hz) which is again in the tail region of the SSC SED as shown in Fig. 2. Very good fit to the multi-TeV spectrum is obtained by using the EBL-D and EBL-I and the parameters are respectively α = 3.1 and Aγ = 58.0 for EBL-D and α = 3.2 and Aγ = 28.0 for EBL-I which correspond to very soft spectrum and the intrinsic spectrum is also soft (between -0.68 to -0.58). In the low energy limit the spectrum shoots up very high and this behavior is not seen by HESS. So we can ignore the fit for which α > 3. This soft power-law problem arises because β = 0.48 is small and we can use the leptonic models which have β > 0.48 as a result of which we should get α > 3. The time averaged differential energy spectrum of this observation is also fitted with a power-law with exponential cut-off having four parameters [61].
For lep-2 with Γ = D = 25, the energy range 1.67TeV(4.0 x 1026 Hz) ≤ Eγ ≤ 20.95TeV(5.0 x 1027 Hz) corresponds to the seed photon energy in the range 0.90MeV(2.17 x 1020 Hz) ≤ ∈γ ≤ 11.26MeV(2.72 x 1021 Hz), which is again in the tail region of the SSC SED as shown in Fig. 5. This is fitted with a power-law with β = 1.1 and Φ0 = 4.37 x 10-9 TeV cm-2s-1. Again EBL-D and EBL-I are used to fit the 2010 flare data in the photohadronic model which are also shown in Fig. 5. The best fit parameters are α = 2.6 and Aγ = 2.5 for EBL-D and α = 2.6 and Aγ = 2.4 for EBL-I respectively also the flux decreases towards low energy regime as expected. Both the EBL-D and EBL-I corrections to the photohadronic model give practically the same result.
In the low energy regime, the lep-1 fit flux increases drastically but this behavior is absent with the lep-2 fit. It is to be noted that, lep-1 corresponds to the observation during the year 2003-2004 and lep-2 is the recent one of January 2013. So we believe that during each observation period, the photon density distribution in the jet changes and this change in the seed photon changes the spectral behavior of the observed multi-TeV gamma-rays. This implies almost simultaneous observation in multiwavelength is essential to fit the observed data. A minor difference between EBL-D and EBL-I predictions for 0.6TeV ≤ E γ ≤ 20 TeV is observed for 2010 flare data as shown in Fig. 6 and the intrinsic flux is Fγ,in α F-07γ Comparison of the intrinsic flux Fγ,in of 2004 and 2010 multi-TeV flaring shows that different spectral shape of the observed events are solely due to the diversity in the shape of the seed photon density distribution (particularly in the SSC tail region) during different epochs.
4.2 Markarian 501
Markarian 501 (RA:251.46°, DEC:39.76°) is one of the brightest extragalactic sources in X-ray/TeV sky [43] and also the second extragalactic object (after Mrk 421) identified as VHE emitter by Whipple telescope in 1996. Since its discovery, the multiwavelength correlation of Mrk 501 has been studied intensively and during this period it has undergone many major outbursts on long time scales and rapid flares on short times scales mostly in the X-rays and TeV energies [80-90]. In the year 2009, Mrk 501 was observed as a part of large scale multiwavelength campaign covering a period of 4.5 months (from March 9 to August 1, 2009) [63]. The scientific goal of this extended observation was to collect a simultaneous, complete multifrequency data set to test the current theoretical models of broadband blazar emission mechanism. Between April 17 to May 5, this HBL was observed by both space and ground based observatories, covering the entire electromagnetic spectrum even including the variation in optical polarization [63]. A very strong VHE flare was detected on May 1st first by Whipple telescope and 1.5 hours later with VERITAS. Both of these telescopes continued simultaneous observation of this VHE flare until the end of the night and the detected flux was enhanced by a factor of ∼ 10 above the average baseline flux. Also a dramatic increase in the flux by a factor ∼ 4 in 25 minutes and a falling time of ∼ 50 minutes was observed. The flux measured at lower energies before and after the VHE flare did not show any significant variation. But, Swift-XRT (in X-ray) and UVOT (in optical) did observe moderate flux variability [63]. Using the one-zone SSC model, the average SED of this multiwavelength (up to second peak) campaign of Mrk 501 is interpreted satisfactorily [63].
The very strong VHE flare data of May 1st observed by Whipple telescope and the long outburst observed by HEGRA telescopes in 1997 were modeled using the photohadronic model of Ref.[39]. The model of Dominguez et al. [46] is used to correct for the EBL effect on the observed data. Here we shall discuss about the fit to this flare data by different models. To explain the VHE flare of May 1st, 2009, we use the parameters of the one-zone leptonic model [63] whose parameters are shown in Table II.
Parameter | Description | Value |
MBH | Black hole mass | (0.9 - 3.5) x 109 Mʘ |
z | Redshift | 0.034 |
Γ | Bulk Lorentz Factor | 12 |
D | Doppler Factor | 12 |
R’b | Blob Radius | 1.2 x 1016 cm |
B’ | Magnetic Field | 0.03 G |
R’f | Inner blob Radius | 5 x 1015 cm |
α | Spectral index | 2.4 |
The observed VHE flare of May 1st was in the range ∼ 317 GeV ≤ Eγ ≤ 5 TeV. In the present context, the above range of Eγ corresponds to the Fermi accelerated proton energy in the range 3.2 TeV ≤ Ep ≤ 50 TeV which interacts with the inner jet background SSC photons in the energy range 13.6 MeV(3.29 x 1021 Hz) ≥ ϵγ ≥ 0.86 MeV(2.1 x 1021 Hz) and finally produce the observed VHE photons. The range of ∈γ lies in the beginning of the SSC spectrum. A very good fit to the data is obtained for α = 2.4 and Aγ = 89 in in Eq. (19). For comparison we have fitted the data with an exponential cut-off function (dashed dotted curve) and the best fit is obtained for α = 2.6, Ec = 30 TeV and Aγ = 66. Also we have shown the Whipple fit (dashed curve) for comparison, where it is fitted by the function dNγ /dEγ = 9.1 x 10-7 (Eγ/1TeV)-2.1 ph m-2s-1TeV-1. It is observed that the EBL correction to the VHE γ-ray is small but not insignificant (black curve in Fig. 7) which has a faster fall above 10 TeV. The Whipple data fit very well with the above three scenarios. However, above 5 TeV, both the the EBL corrected fit and the exponential fit differ from the Whipple fit. Again the EBL fit and the exponential fit differ above 10 TeV and the former one falls faster than the latter as clearly shown in Fig. 7. Even though all these fit very well with the Whipple data, deviation is obvious in the VHE limit. Again, between March 16th and October 1st, 1997, Mrk 501 was in a flaring state which was monitored in TeV γ-rays with the HEGRA stereoscopic system of imaging atmospheric Cherenkov telescopes (IACTs). During this long outburst period (a total exposure time of 110 h) more than 38,000 TeV photons were detected and the energy spectrum of the source was above 10 TeV. A time-averaged energy spectrum of this observation period was fitted with a power-law accompanied with an exponential cutoff [89]. The same data was reanalyzed with an improved energy resolution [90] and found that except for the highest energy, the two analysis were in very good agreement. At the highest energy the spectrum was found to be much steeper than the conventional analysis. In Fig. 7, along with the 2009 flare spectrum, we have also shown the conventional analysis and the improved energy resolution analysis of the 1997 outburst observed by HEGRA telescopes. The photohadronic model of Ref. [39] fits very well with the reanalysis result of 1997 flare data beyond 10 TeV. The entire SED (from low to VHE) is shown in Fig. 8 and the intrinsic flux (red curve in Fig. 7) is plotted to show the EBL contribution.
By comparing different time scales and the Eddington luminosity (as done for Mrk 421), the optical depth is in the range 0.04 < τpγ < 0.13 and this corresponds to the range of photon density in the inner jet region as 1.5 x 1010 cm-3 < n'γ,f < 5.1 x 1010 cm-3. Due to the adiabatic expansion of the inner jet, the photon density will be reduced to n'γ and also the optical depth τpγ << 1. This will drastically reduce the Δ-resonance production efficiency from the ργ process.
4.3 1ES1959+650
The AGN 1ES 1959+650 (z = 0.047) was first detected in the Einstein IPC Slew Survey [91] and classified as a HBL subclass, based on its X-ray to radio flux ratio with a luminosity distance of dL =210 Mpc and the mass of the central black hole is estimated to be ∼ 1.5 x 108 Mʘ. In 1998, VHE gamma-ray from 1ES 1959+650 was observed by the Seven Telescope Array in Utah and later on other observations were also reported. In May 2002, 1ES 1959+650 had a strong TeV outburst which was observed by Whipple [13] and HEGRA experiments [14] as well as in the X-ray range by RXTE experiments. The X-ray flux smoothly declined throughout the following month. However, during this smooth decline period, a second TeV flare was observed after few days (on 4th of June) of the initial one without a X-ray counterpart [15]. As this flare was not accompanied by low energy counterparts, it is called orphan flare. So the observation of the orphan flare in 1ES 1959+650 is in striking disagreement with the predictions of the leptonic models thus challenging the SSC interpretation of the TeV emission. Similar behavior also observed in the flare of April 2004 from Markarian 421. Non observation of a significant X-ray activity could naturally be interpreted by the suppression of electron acceleration and inverse Compton scattering as production mechanism for very high energy (VHE) gamma rays in favor of alternative scenarios. To explain the orphan flare, A hadronic synchrotron mirror model was proposed by Böttcher [93] to explain this orphan TeV flare from 1ES1959+650 [94]. In this model, the flare is explained through the decay of neutral pions to gamma rays when the former are produced due to the interaction of high energy cosmic ray (HECR) protons 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 will be substantial decrease in the HECR proton energy to overcome the threshold for Δ-resonance and, at the same time, it is an alternative to the standard scenario where HECR protons interact with the synchrotron photons, where one needs HECR protons to be Fermi accelerated to very high energy. However, how efficiently these photons will be reflected from the cloud is rather unclear.
The SED of the 1ES1959+650 is fitted quite well with the leptonic one-zone synchrotron and SSC model [95,96]. In all these models, although the blob size differ by about 1 to 2 orders of magnitudes (1.4 x 1014 cm ≤ R'b ≤ 1.4 x 1016 cm), the bulk Lorentz factor Γ and D are almost the same (18 ≤ Γ ≃ D ≤ 20). The multiwavelength observation of 1ES 1959+650 was performed in May, 2006 and the SED fitted with the above one-zone model by Tagliaferri et. al [95], for which the parameters used are Γ ≃ D = 18, R'b = 7.3 x1015 cm and B’ = 0.25 G. To explain the VHE spectrum of the orphan flare in photohadronic model [35], Γ = 18 is also used.
The observed flare energy was in the range 1.26TeV(3.05 x 1026Hz) ≃
9.4TeV(2.3 x 1027Hz). Gamma-ray in this energy range is produced when protons in
the energy range 12 TeV ≤ Ep ≤ 94
TeV collide with the background photons in the energy interval 7.5
MeV(1.8x1021Hz) ≥ ϵγ
≥ 1MeV(2.4x1020 Hz). The range of ∈γ lies exactly in
the low energy tail of the SSC photons as shown in Fig. 9, calculated using the one-zone leptonic model and flux in
this region also has a power-law behavior
ΦSSC
Best fit to the observed VHE spectrum is obtained with the values of α =2.83 and Ec = 4.2 TeV [14] for power-law with exponential cut-off. The γ-ray cut-off energy of 4.2 TeV corresponds to Ep,c = 42 TeV and above the cut-off energy the flux decreases. For power-law with EBL correction, the best fit is obtained for α = 2.8 and Aγ = 87 (black curve) using Eq. (19). For comparison, both the exponential and the EBL fits are shown in Fig. 10. The exponential fit falls faster than the EBL fit above 16 TeV. The time averaged TeV energy spectrum (above 1.4 TeV) of the flaring state of the 1ES 1959+650 was well fitted with pure power-law by the HEGRA collaboration and the power-law spectral index is α=2.83±0.14stat±0.008sys or by a power-law with an exponential cut-off at Ec = (4.2+0.8-06stat ± 0.9sys) TeV and a spectral index of 1.83±0.15stat±0.08sys [14,97].
In Fig. 9, the flare data is shown along with the fit. Normally, it is observed that the flux increases for Eγ = 1.2 TeV due to the high proton flux in this energy range. However, in order not to violate the Eddington luminosity, the proton energy spectrum must break to a harder index (e. g.α ∼ 2.3) below 12 TeV. So a break in proton spectrum is introduced at Ep,b ∼ 12 TeV below which α = 2.3 and above this energy α = 2.8 so that the gamma-ray flux falls below 12 TeV.
5. General Remark
A general discussion of shortcomings in different models is discussed here. In proton synchrotron models [55], the interaction of high energy protons with the synchrotron photons in the jet can produce γ-rays from π0 decay and can explain the multi-TeV emission from blazars. Also Zdziarski et al. [98,99] have used the hadronic model to explain the broad band spectra of radio-loud AGN. These scenarios require super Eddington luminosity in protons (∼ 106 times the Eddington flux) to explain the multi-TeV emission. Also, synchrotron emission from the ultra high energy protons in the jet magnetic field can explain the VHE γ-ray SED [68] but strong magnetic field at the emission site is a necessary requirement which is probably not there. In an alternative scenario, ultra high energy protons escaping from the jet region produce VHE photons by interacting with the cosmic microwave background (CMB) photons and/or EBL which avoids the absorption in the inner jet region [100]. This explains the transparency of the universe to VHE γ-rays due to their proximity to the Earth compared to the one produced in the source which travels a longer distance. Also the TeV spectrum is independent of the intrinsic spectrum but depends on the output of the high energy cosmic rays in the source. This model fits very well to multi-TeV spectra from many sources [69,101-103]. However, in this scenario, it is assumed that the source produces VHE protons with energies 1017 - 1019 eV (this is true with most of the hadronic models) and a weak extragalactic magnetic field in the range 10-17 G < B <10-14 G is needed. So far we have not observed either high energy cosmic ray event or high energy neutrino event from nearby HBLs (e.g. Mrk 401, Mrk 501 etc.) when they were flaring and also the extragalactic magnetic field seems to be weak. However, in the photohadronic model discussed above, ργ→Δ process takes place in the hidden inner jet region where the photon density is order of magnitude higher than the normal jet and overcomes the super-Eddington energy budget. Also, the energy of the Fermi accelerated proton is Ep = 10Eγ which can be achieved easily in the jet.
Normally it is difficult to explain the GeV-TeV emission from HBLs through one-zone leptonic model. On the other hand, multi-zone leptonic models overcome this problem and explain the multi-TeV data but one has to increase the number of parameters. A power-law with an exponential cut-off energy Ec is the conventional method used in the literature to explain the exponential fall of the multi-TeV flare from flaring blazars where the free parameter Ec depends on some unknown mechanism. The EBL correction does the same job as the exponential decay factor. So the inclusion of EBL effect eliminates the necessity of the extra parameter Ec.
The high energy protons will be accompanied by high energy electrons and these electrons will emit synchrotron photons in the magnetic field of the jet. The energy range of the photons emitted lie in between the high energy end of the synchrotron spectrum and the low energy tail of the SSC spectrum, thus may not be observed due to their low flux in this region. The secondary charged particles e±, π±, μ± produced within the jet will also emit synchrotron photons and their energy will be much smaller than the synchrotron photons produced from the primary electrons accompanying the high energy protons. So even though the charged particles emit synchrotron radiation, their contribution will be small.
The problem with the photohadronic model is that, the gamma-rays are produced from the pion decay within the inner jet region where photon density is very high compered to outer region. So γγ → e+e- should be efficient, which will attenuate the propagation of multi-TeV gamma-rays from the production site. However, observation of multi-TeV gamma-ray events from many flaring HBLs observed by MAGIC, VERITAS, HESS telescopes guarantees that gamma-rays do escape from the jet without undergoing pair production, so the pair production process can be constrained. The photohadronic model works well for high energy gamma-rays (above ∼ 100 GeV), so in the lower energy regime, it is the leptonic model which contributes to the multiwavelength SED. In principle both should coincide at some intermediate energy.
6. Summary
One-zone leptonic model is very successful to explain the synchrotron and SSC peaks
of AGN. But this model has difficulties to explain the multi-TeV emission/flaring
from many HBLs. In this article, we review the photohadronic model which is very
successful to explain the multi-TeV flare data from the nearest HBLs Mrk 421, Mrk
501 and 1ES1959+650. As high energy gamma-ray is attenuated by the EBL, it is
necessary to account for it which is considered here and compared with the
exponential cut-off scenario and other power-law fits. The main objective here is to
understand the flaring from nearby HBLs so that the mechanism to produced intrinsic
flux can be understood well and also it can be useful to understand the flaring from
far off HBLs. This model uses the leptonic model parameters which explains the first
two peaks of the AGN very well as input. In the photohadronic scenario, the observed
flux is proportional to the SSC flux
ΦSSC which is a power-law,
ΦSSC