Hard MeV–GeV spectra of blazars
Toruń Centre for Astronomy, Nicolaus Copernicus University, ul. Gagarina 11, 87100 Toruń, Poland
email: kat@astro.uni.torun.pl
Received: 7 March 2011
Accepted: 11 October 2011
Aims. Very high energy (VHE) gammaray emission from a distant source (z ≳ 0.2) can be efficiently absorbed by means of the electronpositron pair creation process. Analyses of the unabsorbed spectra imply that the intrinsic TeV emission of some blazars is hard, with spectral indices of 0.5 < α < 1. The absorption depends on the level of extragalactic background light (EBL) that is difficult to measure directly. This implies that it is difficult to estimate the slope of the intrinsic TeV emission. To test our blazar emission scenario that is capable of reproducing the hard spectra, we therefore used the observations made by the Fermi Gammaray Space Telescope in the unabsorbed MeV−GeV energy range.
Methods. We assume that the Xray and gammaray emission of TeV blazars is produced in a compact region of a jet uniformly filled by particles of relatively high energy (γ ≳ 10^{3}, E = γm_{e}c^{2}). In other words, we assume a low energy cutoff in the particle energy distribution. The emission produced by the particles with this energy spectrum can explain hard intrinsic spectra in the energy range from MeV up to TeV. We demonstrate how to estimate the basic physical parameters of a source in this case and how to explain the observed spectra by a precise simulation of the particle energy evolution.
Results. To test our estimation methods, we use the observations of two blazars with exceptionally hard spectral indices (α ≲ 0.5) in the MeV − GeV range and known redshifts: RGB J0710+591 and 1ES 0502+675. The estimated values of the Doppler factor and magnetic field are compared with our numerical simulations, which confirm that the particle energy distribution with a low energy cutoff can explain the observed hard spectra well. In addition, we demonstrate that the radiative cooling caused by the inverseCompton emission in the KleinNishina regime may help us to explain the hard spectra.
Key words: galaxies: active / BL Lacertae objects: individual: RGB J0710+591 / BL Lacertae objects: individual: 1ES 0502+675 / radiation mechanisms: nonthermal
© ESO, 2012
1. Introduction
The emission of some blazars is observed from radio frequencies up to VHE gamma rays. Collating observations from different energy ranges, we can show that the spectra of these objects contain two characteristic peaks (in νF(ν) plots). The first peak appears in the Xrays from a few keVs to a few hundred keVs and the second peak is observed around an energy of a few TeVs (e.g., Ghisellini et al. 1998, 2010; Massaro et al. 2004; Nieppola et al. 2006; Abdo et al. 2010b). Since the TeV emission is the most prominent feature in those sources we often call them TeV blazars, which constitutes a relatively small group of sources within the blazar family. The gammaray emission of most blazars peaks in the MeV − GeV range. Moreover, TeV gamma rays are efficiently absorbed by the EBL. This in addition limits the number of observed TeV blazars.
The high energy emission of blazars is believed to originate inside their jets. However, the observed variability on timescales from days (e.g., Catanese et al. 1997; Fossati et al. 2008) down to a few minutes (e.g., Aharonian et al. 2007a) indicates that only a small part of a jet is radiating at high energies. This kind of emission requires high energy particles, which can gain energy by acceleration at the front of a shock wave inside the jet. This is the simplest explanation of the acceleration process, which assumes that a small fraction of the jet bulk kinetic energy is transferred to the particles. In the presence of a magnetic field, the particles can radiate this energy through synchrotron emission, generating the first peak in the spectrum. Some fraction of the synchrotron photons can then be upscattered to higher energies. This is the inverseCompton (IC) scattering that gives the second peak in the TeV energies. This simple scenario has been proposed many times as part of models of the VHE emission of blazars (e.g., Dermer et al. 1997; Bloom & Marscher 1996; Inoue & Takahara 1996; Mastichiadis & Kirk 1997; Katarzyński et al. 2001).
A TeV gamma ray photon travelling through the intergalactic medium can interact with an infrared photon producing an electronpositron pair. This causes absorption of the emission above a few hundred GeVs. To calculate this absorption, we have to determine the level of the intergalactic radiation field, which is difficult to measure directly but can be estimated from simulations of star light production in evolving galaxies. There are many different solutions for the level of the absorption (e.g., De Jager & Stecker 2002; Keneiste et al. 2004; Franceschini et al. 2008; Kneiske & Dole 2010). Thus the slope of the intrinsic TeV emission cannot be accurately calculated. Nevertheless, the TeV emission of some relatively distant sources (e.g., 1ES 1101232 z = 0.186, Aharonian et al. 2006; or 1ES 0229+200, z = 0.14, Aharonian et al. 2007b) indicates that the intrinsic spectra are hard (F ∝ ν^{ − α} with 0.5 < α < 1), even if we assume as low as possible absorption.
It is difficult to explain how these hard spectra could be created. Most of the standard emission models assume a powerlaw or broken powerlaw particle energy distribution with a slope n ≃ 2 (where the number of particles N ∝ γ^{ − n}) and a minimum energy equivalent to γ_{m} ≃ 1. This gives the spectral index of the synchrotron emission α = (n − 1)/2 = 0.5. This emission is upscattered to higher energies by the same population of the particles, hence the index of the inverseCompton spectrum should be similar to the index of the synchrotron emission. Note that high energy particles (γ ≳ 10^{5}) do not scatter efficiently because of the KleinNishina (KN) restrictions. The broken powerlaw particle distribution (n_{1} and n_{2} below and above the break, respectively) implies that there is a broken synchrotron spectrum (α_{1} and α_{2}), for which the spectral index in the KN regime is approximately α ≃ 2α_{2} − α_{1} = 2.5 for n_{2} = 4 (Tavecchio et al. 1998).
A question to ask is how we can explain the hard spectra from MeVs up to TeV gamma rays. One of the simplest solution is to assume a cutoff in the low energy part of the electron spectrum, which basically means γ_{m} ≫ 1. This solution was proposed for the first time by Katarzynski et al. (2006a) to explain the observations of 1ES 1101232 in the TeV range. This cutoff leads to an additional break in the synchrotron spectrum between IR and Xray energies. The spectral index of the first part of the spectrum, from radio up to Xray energies, becomes a constant α_{0} = −1/3. This is a “tail” of the synchrotron emission produced by the low energy particles. This part of the radiation field with such a hard spectral index can also be upscattered. Therefore, the IC spectrum can also be very hard up to limiting value α = −1/3. Simultaneous optical and Xray observations of 1ES 0229+200 performed by the Swift satellite (Tavecchio et al. 2009) shows an abrupt break between the optical and the Xray range. This strongly implies that there is a low energy cutoff to the particle energy distribution.
In the present work, we explore this idea of a low energy cutoff. We focus on the unabsorbed MeV − GeV range where the spectral index of the intrinsic emission is observed directly. We demonstrate how to estimate the basic parameters of a source and how to explain the observed spectra by a simulation of the particle energy evolution. We apply our estimations and modelling to the Xray and gammaray observations of two blazars with known redshifts.
2. Observations
The spectra of TeV blazars in the MeV−GeV range are always hard with spectral indices α < 1. The average expected value of the spectral index in this particular range is α ~ 0.5. In this work we focus on the sources with α < 0.5 that we call exceptionally hard.
For all blazars observed by the Fermi Gammaray Space Telescope (Atwood et al. 2009) from 20 MeV up to 300 GeV, only in some cases were spectra published in the Fermi Large Area Telescope First Source Catalog (Abdo et al. 2010a) exceptionally hard with the photon index Γ ≲ 1.5 (where the photon index is greater by unity in comparison to the spectral index). Moreover, in many cases we either do not know the distance to the observed sources (e.g. RBS 158 Γ = 1.34 ± 0.33, RBS 621 Γ = 1.42 ± 0.32) or the spectra were obtained with a relatively large uncertainty (e.g. 1ES 1101232 Γ = 1.36 ± 0.58). Therefore, we selected only two objects for our test.
The observational parameters of the first object RGB J0710+591 (z = 0.125) in the Fermi catalogue are Γ = 1.28 ± 0.22, and F = (1.46 ± 0.49) × 10^{11} erg cm^{2} s^{1}, where F is the energy flux in the 100 MeV to 100 GeV range. This source was also discovered by the VERITAS gammaray observatory in the TeV range (Acciari et al. 2010). The photon index obtained by VERITAS above ~300 GeV was Γ = 2.69 ± 46. In our analysis, we also use Xray observations obtained by the Swift satellite in the energy range from 0.2 to 10 keV (Acciari et al. 2010). A notable characteristic of this object is that the Xray peak is observed at the level a few times higher than the GeVTeV peak. This indicates that the synchrotron cooling is dominating in this source.
The photon index of the second source 1ES 0502+675 (z = 0.341) in the Fermi catalogue is rather soft Γ = 1.75 ± 0.11. However, the MeV − GeV spectrum of this source is quite complex. Owing the limited photon statistics, the observations published by Abdo et al. (2010b,c) seem to display a broken powerlaw spectrum of IC emission. This is rather unexpected and not predicted by any model of blazar emission. However, this may have a simple explanation if we assume for example that two independent jet components were observed at the same time. The break in the spectrum appears around an energy of 1 GeV. Fitting a single powerlaw function to the observations above this energy, we obtained Γ = 1.35 ± 0.20 and Γ = 1.46 ± 0.18 for the data sets above quoted, respectively. The VERITAS team reported the detection of this source in the TeV energy range (Ong 2009). However, no public information about the observed flux and slope has been publicly available so far. In addition we use quasisimultaneous observations from Swift (Abdo et al. 2010b). The comparison between the Xray peak level and the GeV observations shows that in this particular source the radiative cooling generated by the IC emission is the predominant cooling process.
3. Estimations
To estimate the basic parameters of a source of VHE emission, we adapted approach proposed by Tavecchio et al. (1998). We assume a spherical source of radius R, filled uniformly by a constant magnetic field (B) and relativistic particles (K – density). In addition, we assume that the particle energy distribution is a broken powerlaw (n_{1},n_{2},γ_{b} – break energy) and that the source travels with relativistic velocity (δ – Doppler factor) inside a jet. Finally, we included a low energy cutoff in the particle energy spectrum (γ_{m} – minimum energy). This is our main extension to the initial approach that let us estimate the basic physical parameters from the exceptionally hard spectrum. Even this simple scenario has eight free parameters. Therefore, using observational quantities we can only estimate the relation between the main parameters B and δ. However, we can estimate this relation in four different ways.
Fig. 1 Approximate synchrotron and IC spectra calculated for different values of γ_{m}, from almost unity up to γ_{b}. This test demonstrates how the spectral index in the MeV − GeV range depends on the γ_{m} value. The dashed line shows the difference between the estimated (α_{1} = 0.5) and calculated (α ≃ 0.65) spectral slope below the peak. 

Open with DEXTER 
3.1. Peak positions
The simplest estimation of the B − δ relation involves the comparison of the synchrotron peak position (ν_{s}) with the inverseCompton peak frequency (ν_{c}). The positions of the peaks do not depend on γ_{m} thus we can directly use the formula derived by Tavecchio et al. (1998) (1)In Fig. 1, we illustrate how an approximate IC spectrum evolve with time for different values of γ_{m}, from γ_{m} ≳ 1 up to γ_{b}. The spectra are calculated according to the prescription given in the Appendix. This simple test shows that the IC peak remains at the same position for the wide range of γ_{m} values with the exception for γ_{m} ≃ γ_{b} and this gives the limit to the above formula. The IC spectrum with α = −1/3 that could be obtained for γ_{m} = γ_{b} was never observed directly. The observed spectral indices described in the previous section are in range from 0.3 to 0.5. Therefore, we can safely use this B and δ relationship for our estimations.
The spectra in Fig. 1 were calculated for n_{1} = 2 and n_{2} = 4, which directly gives the spectral indices of the synchrotron emission below (α_{1} = 0.5) and above (α_{2} = 3/2) the peak, respectively. The value of n_{1} = 2 is rather typical for many different astrophysical sources. This may be the result of the firstorder Fermi acceleration at a shock front or the synchrotron or the inverseCompton cooling in the Thompson regime. This index is crucial for the emission of TeV blazars because this part of the particle spectrum provides synchrotron photons for the IC scattering. The IC spectrum in the MeV − GeV range comes mostly from the scattering in the Thompson regime. Therefore, the spectral index of this emission should be similar to the index of the scattered radiation field. However, the complexity of the scattering (see Fig. A.1) produces a curved spectrum below the peak. In the MeV − GeV range, just below the peak, the spectrum can be approximated well by a powerlaw function with the index α ≃ 0.65. This differs significantly from the estimated value α = 0.5 that is clearly illustrated in Fig. 1. Note that this concerns only the case where γ_{m} ~ 1 because γ_{m} controls the slope in the other cases. Since the somewhat classical value of n_{1} = 2, postulated by many different particle evolution scenarios, gives α = 0.65 all spectra with α < 0.65 should be classified as exceptionally hard. Using this revised criterion, all spectra considered here may be regarded as exceptionally hard.
3.2. Peak levels
The second relationship between B and δ derived by Tavecchio et al. (1998) was obtained from the wellknown formula (2)where U_{syn} and U_{B} are the energy densities of the synchrotron radiation field and the magnetic field, respectively, and L_{s} and L_{c} are the total luminosities of the synchrotron and the IC emission, respectively. The luminosities were calculated from the observed spectral slopes and the peak emission levels. Unfortunately, we cannot adopt this approach in our estimations because usually we do not observe the break in the synchrotron emission (ν_{m}) that is related to γ_{m}. We instead compare (3)the emissivities (Eqs. (A.2), (B.1)) with the observed emission levels at the peaks. This gives the particle density (Eq. (B.2)) for the first part of the particle energy distribution (4)where (5)is the simplified version of the IC emissivity calculated as a sum of two dominant components j_{c,sim} = j_{(1,1)} + j_{(1,2)} (Eqs. (B.8), (B.9)). Note that in the component j_{(1,2)} the lower integration boundary was neglected. Finally, using the synchrotron emissivity (Eq. (A.2)) and the transformation given by Eq. (A.5) we can derive another relation between B and δ(6)where the luminosity distance (D_{L}) is calculated in a ΛCDM universe for h = 0.72, Ω_{Λ} = 0.7, and Ω_{m} = 0.3 for all estimations and the spectra calculated in this work.
3.3. Spectral index in MeV–GeV range
Another constraint on the basic physical parameters can be obtained from the slope of the gammaray emission in the MeV–GeV range. This constraint assumes that the particle minimal energy (γ_{m}) is directly related to the radiative cooling inside the source. This is one of the main assumptions in our model that we discuss in detail in the next section. According to this assumption (7)where t_{cool} is the characteristic cooling time equivalent to the source evolution time in this particular case. This and the synchrotron energy density (U_{syn}) are unknown parameters. Therefore, for our estimation we assume that U_{syn} = x_{B}U_{B} is energy independent and can be parametrized as some fraction (x_{B}) of U_{B}. In the same way, we parametrize t_{cool} = x_{R}R/c as the multiplied crossing time. Since the spectral peaks are observed at almost the same level we may conclude that x_{B} ~ 1. In addition, we know that x_{R} must be greater than unity, although, in our estimations we assume that x_{R} ≥ 2. This gives the shortest reasonable time for the source evolution. With the above assumptions, we can derive an upper limit to the magnetic field value (8)where γ_{m} value required by the above formula can be obtained from the spectral slope in the MeV − GeV range. However, there is no simple analytic formula that gives γ_{m}.
We derive the the minimum energy (γ_{m}) generating spectra from γ_{m} = 1 up to γ_{m} = γ_{b} and searching for the required spectral index (9)where and are the corresponding frequencies to the energies 100 MeV and 100 GeV, respectively, in the comoving frame. The frequency transformation, which is required for this calculation, introduces the Doppler factor dependency into the estimation. Note that it is necessary to use the full emissivity (Eq. (B.1)) for this particular estimation. The emissivity can be calculated for any nonzero value of R,B, and K because the absolute level of the emission is unimportant in this calculation. Finally, to calculate the emissivity we need the energy spectrum break that is given by a simple formula (10)derived by Tavecchio et al. (1998).
3.4. Pair absorption
The optical depth for the pair absorption inside a spherical source can be approximated as (11)where is the number density of soft photons per energy interval (e.g., Coppi & Blandford 1990). The above formula incorporates a simple relation between the soft photon and the gammaray photon energy . Using the flux transformation (Eq. (A.5)), we can write (12)which should be smaller than unity for an optically thin source. This gives a lower limit to the value of the Doppler factor (13)where ν_{0} = (m_{e}c^{2}/h)^{2}/ν_{c}. Note that this formula is more accurate for a spherical source than the relationship derived by Dondi & Ghisellini (1995) used in the paper by Tavecchio et al. (1998). Our formula infers a lower limit that is almost five times higher than the above quoted relationship.
3.5. Estimation results
It is difficult to determine the position of the peaks in RGB J0710+591 from the available observations. The Swift energy range is quite narrow and the gammaray observations made by Fermi and VERITAS are not simultaneous. We analysed therefore two alternative solutions focusing first on the Swift and Fermi observations and then using the Swift and VERITAS data with different assumption about the peak positions.
To estimate the basic physical parameters of RGB J0710+591 we assume both the position of the synchrotron peak (ν_{s} = 7 × 10^{17} Hz) and the emission level at the peak (ν_{s}F_{s}(ν_{s}) = 3 × 10^{11} erg cm^{2} s^{1}). It is also necessary to assume the frequency of the IC peak (ν_{c} = 9.5 × 10^{25} Hz) and the emission level at this peak (ν_{c}F_{c}(ν_{c}) = 10^{11} erg cm^{2} s^{1}). However, for this peak we assume additional discrepancy of about ± 30% in both the position and the emission level as well. This discrepancy is introduced because of the unknown extragalactic absorption. The spectral index of the synchrotron emission below the peak is assumed to be α_{1} = 0.5 (n_{1} = 2), whereas above the peak we use α_{2} = 3/2 (n_{2} = 4).
Fig. 2 The upper panel shows the high energy emission of RGB J0710+591 and approximated spectra used for our estimations. In the lower panel, we illustrate the constraints to the basic parameters obtained in four different ways and the values selected for the modelling. 

Open with DEXTER 
Fig. 3 The upper panel shows the multifrequency emission of 1ES 0502+675 observed by the Swift & Fermi instruments (Abdo et al. 2010b – squares; Abdo et al. 2010c – circles) and the approximate spectra used for the estimation presented in the lower panel. The values of B and δ selected for our modelling are indicated by the dot. 

Open with DEXTER 
With the above assumptions, we can estimate the values of B and δ in the four different ways described in the previous subsections. The spectra used for the estimation and the estimation result are presented in Fig. 2. On the basis of the constraints obtained from the peak positions (Eq. (1)) and the emission levels (Eq. (6)), we found two crossing areas because of the discrepancy in the IC peak. The two other methods provide upper limit to the magnetic field value and the lower limit to the Doppler factor.
To calculate some of the limiting curves, we have to determine the radius of the emitting region. This parameter in principle can be constrained from the observed variability timescales R ≲ ct_{var}δ/(1 + z). However, there is no information about the variability of RGB J0710+591. We therefore selected as large as possible a value of the radius R = 2 × 10^{16} cm, which still provides a good agreement for the four different estimation methods.
In the second approach, we focus on the VERITAS observations of RGB J0710+591 assuming that the IC peak is at the frequency ν_{c} = 2 × 10^{26} Hz and the emission level at this peak is ν_{c}F_{c}(ν_{c}) = 9.5 × 10^{12} erg cm^{2} s^{1}. In addition me must assume that the synchrotron peak is placed above the observed range (ν_{s} = 3 × 10^{18} Hz and ν_{s}F_{s}(ν_{s}) = 4 × 10^{11} erg cm^{2} s^{1}). This assumption is similar to the solution proposed by Acciari et al. (2010) that helps us to explain a “flat” (Γ ~ 2) intrinsic TeV spectrum.
In the case of 1ES 0502+675, we assumed the same parameters for the synchrotron peak as for RGB J0710+591 and that ν_{c} = 2.5 × 10^{26} Hz, and ν_{c}F_{c}(ν_{c}) = 1.5 × 10^{10} erg cm^{2} s^{1} for the IC peak. Moreover, we used α_{1} = 0.375 (n_{1} = 1.75) and α_{2} = 2 (n_{2} = 5). We explain this particular choice in the next section.
The main difference between this estimation and the calculations performed for RGB J0710+591 is in the upper limit to the magnetic field obtained from the spectral slope in the MeV–GeV range. This estimate was calculated under the assumption that x_{B} ≃ 1 (U_{syn} ≃ U_{B}). However, the observations indicate that the emission at the IC peak, in this particular source, might even be one order of magnitude higher than the emission at the synchrotron peak. This would indicate that x_{B} ≫ 1 but a more precise value of this parameter cannot be estimated from the observations because the IC scattering at the peak is already in the KleinNishina regime. This illustrates a limitation of this particular estimation method. Finally, we selected as large as possible radius (R = 10^{16} cm), which provided a good agreement between our four different methods of the estimation.
4. Modelling
The information that we can obtain from the analysed observations is significantly limited. We therefore decided to apply quite simple model of the emission that is directly related to our estimation methods. The basic assumptions of the model that we adopted are identical to the assumptions made for the estimations. The VHE originates in a spherical source uniformly filled by relativistic particles and an entangled magnetic field. The source evolves in time, “fresh” particles are continuously injected into the source where they lose energy by means of the synchrotron and the IC emission. This is a simple adaptation of the jet internal shock model, where the particles efficiently accelerated at a shock front fill the downstream region of the shock.
The evolution of the particle energy spectrum inside the source is described by the kinetic equation (14)where is the synchrotron cooling rate and the IC cooling rate () is calculated according to the prescription of Sauge & Henri (2004) for an isotropic distribution of soft photons, using the full KleinNishina crosssection in the headon approximation (Jones 1968). The injection rate is described by a powerlaw function (15)where γ_{i} ≫ 1 describes the minimum energy of the injected particles. Finally, we simulate the pair creation process within the source, where the pair injection rate P(γ,t) is calculated using the approach described by Sauge & Henri (2004). This process in principle should be negligible when we use the parameters obtained from our estimations. However, in the estimation our intention was to use as large as possible a radius that gives a relatively small value of the Doppler factor. The value is only slightly larger than the lower limit obtained from our simple estimation (Eq. (13)). Therefore, the pair creation process may have a small impact on the evolution of the particle energy spectrum in our more precise simulation. To solve the kinetic equation, we used the numerical method described by Chiaberge & Ghisellini (1999). Calculating spectra for each time step, we are carefully checked the energetic balance between the injected and radiated energy.
Such simple timedependent modelling has been already proposed many times to explain evolution of the TeV blazars emission (e.g., Dermer et al. 1999; Kusunose et al. 2000; Böttcher & Chiang 2002; Sauge & Henri 2004). The main difference between our approach and the other models is in the minimum injected energy, where in the other approaches γ_{i} = 1 is usually assumed. In other words, we assume that the acceleration process at the shock front is very efficient and “pushes” almost all particles towards high relativistic energies. This is somehow similar to the formation of either a thermal or quasithermal particle energy distribution (e.g., Schlickeiser 1984; Katarzyński et al. 2006b).
To obtain an exceptionally hard spectrum in the MeV − GeV range, we needed to apply a cutoff to the low energy part of the spectrum. This cutoff can be produced in a natural way by the radiative cooling. Injecting only high energy particles, we may obtain a broken powerlaw particle distribution with the break at the γ_{i} energy. The minimum energy (γ_{m}) that the particles can reach in the source evolution time is given by Eq. (7) (e.g., Kardashev 1962). This shows that in principle, the particles might cool down to γ ~ 1, if the source evolution is long enough. However, the minimum energy is inversely proportional to time, which means that the high energy particles (γ ~ 10^{5}) lose energy more rapidly than the medium energy particles (γ ~ 10^{3}). Moreover, in the downstream region of the shock there may be another acceleration process, for example a turbulent secondorder Fermi acceleration. This relatively weak process may not produce any high energy activity but it can be strong enough to compensate at some point for the radiative cooling and to keep the minimum energy significantly above the particle rest energy. This is, however, a much more complex approach that requires an analysis of possible activity. With no information about timedependent flux variations, we simply assume that the source evolution time is just 2R/c.
The spectral index of the particle energy spectrum below the break will be constant n_{1} = 2, whereas above the peak the index depends on the injection slope n_{2} = n_{i} + 1. This, however, is true only for the radiative cooling caused by either the synchrotron emission or the IC emission in the Thompson regime. In the KleinNishina regime, the scattering efficiency of the high energy particles is significantly lower. This modifies the high energy part of the energy distribution (e.g., Moderski et al. 2005). It is difficult to describe this effect analytically in our particular scenario. We may perform only a simple estimation assuming that the energy density of the synchrotron radiation field is given by (16)where is the minimum frequency of the scattered photons and introduces the KN restriction in a simple form (e.g., Chiaberge & Ghisellini 1999). Neglecting the lower boundary in the above integral and assuming that and I_{s} ∝ (ν′)^{ − α1}, we may approximate the IC cooling ratio as . Hence, the stationary (Ṅ = 0) solution of the kinetic equation is given by N ∝ γ^{ − 2 − α1}, where the slope is significantly lower than the value N ∝ γ^{2} obtained for a constant classical synchrotron or IC cooling. This agrees with our more precise numerical tests, which indicate that the radiative cooling dominated by the IC scattering in the KN regime can reduce the particle energy slope by a factor ≲ 1/2. This appears to be true for the cooled (below the break) and injected (above the break) part of the energy spectrum. Since the emission of 1ES 0502+675 seems to be significantly affected by the scattering in the KN regime, we assume that n_{1} = 1.75 for the estimation made in the previous section for this source. We note that the energy spectrum index less than two helps in addition to explain exceptionally hard MeV − GeV spectra.
Fig. 4 The modelling of the RGB J0710+591 spectral emission. The upper panel shows the first approach where we considered the Swift and Fermi observations only. In the lower panel, we demonstrate an alternative approach that can explain the origin of the TeV emission observed by VERITAS. Note that the Fermi and the VERITAS observations were not simultaneous. To obtain the observed IC spectra, we used the EBL lowerlimit model proposed by Kneiske & Dole (2010), where the dashed lines show the intrinsic unabsorbed emission. 

Open with DEXTER 
Figure 4 shows the modelling of the RGB J0710+591 emission. We present the results obtained for two different sets of physical parameters derived from the estimations performed in the previous section. Note that the IC spectrum obtained in the first case is unable to explain the intrinsic shape of the emission. In the second case, the intrinsic spectrum is reproduced well but the level of the MeV − GeV emission is significantly lower than the observed level. However, the Fermi and the VERITAS observations are not simultaneous, hence there may be a difference in the emission levels. The detailed values of the physical parameters used for the modelling are given in Table 1. There are only six important parameters illustrating the simplicity of the model. In addition, we assumed that the maximum Lorentz factor of the injected particles is γ > 10^{7} and the source evolution time t_{evo} = 2R/c.
For 1ES 0502+675, we do not attempt to reproduce the Swift UVOT observations. This part of the emission might be dominated by either the host galaxy or more extended jet structures as demonstrated for example in the case of Mrk 501 (Katarzynski et al. 2001). Moreover, no direct correlation between Xray and optical variability in other TeV blazars also implies that the opticalUV emission has a different origin. Our modelling is also unable to explain the inverted spectral slope in the 100 MeV to 1 GeV energy range. This is an unexpected and puzzling feature of the emission that will require additional investigation.
However, this may be simply the IC emission produced by the extended jet structures, which also dominate the emission in the opticalUV range.
5. Summary
Several blazars observed by the Fermi GammaRay Space Telescope have exceptionally hard spectra in the MeV − GeV range. We have demonstrated that in principle all sources with a spectral index α ≲ 0.65 should be classified as sources with an exceptionally hard MeV − GeV spectrum. Therefore, we have attempted to reproduce these spectral slopes by assuming a low energy cutoff in the particle energy distribution. This approach had been previously developed to explain hard spectra in the TeV range and at least in one case (1ES 0229+200 Tavecchio et al. 2009) it was found to be correct.
Parameters used for the modelling.
Fig. 5 The spectra obtained from the modelling of 1ES 0502+675 for the physical parameters obtained from our estimation. 

Open with DEXTER 
We derived four different methods to estimate the basic parameters of the emitting region. Three of these methods had been previously proposed and we simply adapted them to the case of the low energy cutoff. Our method that gives an upper limit to the magnetic field, derived from the spectral slope in the MeV–GeV range, was proposed here for the first time. This particular constraint can be very useful but it requires accurate measurements of the spectral index in the notabsorbed MeV–GeV range.
Finally, we used our modelling estimates that succeed to explain the exceptionally hard spectra in terms of particle energy evolution. Our model is relatively simple but more complex, and hence more accurate modelling would require far tighter observational constraints. We have attempted to explain the emission of two different TeV blazars. The particle evolution in RGB J0710+591 seems to be dominated by synchrotron cooling, whereas in 1ES 0502+675 IC cooling appears to dominate. In both cases, the estimated parameters of B and δ provide good results in our modelling. This type of modelling is crucial for understanding the nature of the blazar emission process, especially in nonstandard cases such as the exceptionally hard MeV–GeV emission. One interpretation of exceptionally hard spectra is that for some reason only high energy particles can produce the Xray and gammaray activity of TeV blazars.
Acknowledgments
We thank the anonymous referee for a careful review and a helpful report.
References
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010a, ApJS, 188, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Abdo, A. A., Ackermann, M., Agudo, I., et al. 2010b, ApJ, 716, 30 [NASA ADS] [CrossRef] [Google Scholar]
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010c, ApJ, 710, 1271 [NASA ADS] [CrossRef] [Google Scholar]
 Acciari, V. A., Aliu, E., Arlen, T., et al. 2010, ApJ, 715, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Aharonian, F., Akhperjanian, A. G., BazerBachi, A. R., et al. 2006, Nature, 440, 1018 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Aharonian, F., Akhperjanian, A. G., BazerBachi, A. R., et al. 2007a, ApJ, 664, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Aharonian, F., Akhperjanian, A. G., Barres de Almeida, U., et al. 2007b, A&A, 475, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 [NASA ADS] [CrossRef] [Google Scholar]
 Böttcher, M., & Chiang, J. 2002, ApJ, 581, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Bloom, S. D., & Marscher, A. P. 1996, ApJ, 461, 657 [NASA ADS] [CrossRef] [Google Scholar]
 Catanese, M., Bradbury, S. M., Breslin, A. C., et al. 1997, ApJ, 487, L143 [NASA ADS] [CrossRef] [Google Scholar]
 Chiaberge, M., & Ghisellini, G. 1999, MNRAS, 306, 551 [NASA ADS] [CrossRef] [Google Scholar]
 Coppi, P. S., & Blandford, R. D. 1990, MNRAS, 245, 453 [NASA ADS] [Google Scholar]
 De Jager, O. C., & Stecker, F. W. 2002, 566, 738 [Google Scholar]
 Dermer, C. D., Sturner, S. J., & Schlickeiser, R. 1997, ApJS, 109, 103 [NASA ADS] [CrossRef] [Google Scholar]
 Dermer, C. D., Li, H., & Chiang, J. 1999, ApL&C, 39, 1 [NASA ADS] [Google Scholar]
 Dondi, L., & Ghisellini, G. 1995, MNRAS, 273, 583 [NASA ADS] [Google Scholar]
 Fossati, G., Buckley, J. H., Bond, I. H., et al. 2008, ApJ, 677, 906 [NASA ADS] [CrossRef] [Google Scholar]
 Franceschini, A., Rodighiero, G., & Vaccari, M. 2008, A&A, 487, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ghisellini, G., Celotti, A., Fossati, G., et al. 1998, MNRAS, 301, 451 [NASA ADS] [CrossRef] [Google Scholar]
 Ghisellini, G., Tavecchio, F., Foschini, L., et al. 2010, MNRAS, 402, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Ginzburg, V. L., & Syrovatskii, S. I. 1965, ARA&A, 3, 297 [NASA ADS] [CrossRef] [Google Scholar]
 Inoue, S., & Takahara, F. 1996, ApJ, 463, 555 [NASA ADS] [CrossRef] [Google Scholar]
 Jones, F. C. 1968, Phys. Rev., 167, 1159 [NASA ADS] [CrossRef] [Google Scholar]
 Massaro, E., Perri, M., Giommi, P., et al. 2004, A&A, 413, 489 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kardashev, N. S. 1962, SvA, 6, 317 [NASA ADS] [Google Scholar]
 Kataoka, J., Mattox, J. R., Quinn, J., et al. 1999, ApJ, 514, 138 [NASA ADS] [CrossRef] [Google Scholar]
 Katarzyński, K., Sol, H., & Kus, A. 2001, A&A, 367, 809 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Katarzynski, K., Ghisellini, G., Tavecchio, F., et al. 2006a, MNRAS, 368, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Katarzynski, K., Ghisellini, G., Mastichiadis, A., et al. 2006b, A&A, 453, 47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kneiske, T. M., & Dole, H. 2010, A&A, 515, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kneiske, T. M., Bretz, T., Mannheim, K., & Hartmann, D. H. 2004, A&A, 413, 807 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kusunose, M., Takahara, F., & Li, H. 2000, ApJ, 536, 299 [NASA ADS] [CrossRef] [Google Scholar]
 Mastichiadis, A., & Kirk, J. 1997, A&A, 320, 19 [NASA ADS] [Google Scholar]
 Moderski, R., Sikora, M., Coppi, P. S., & Aharonian, F. 2005, MNRAS, 363, 954 [NASA ADS] [CrossRef] [Google Scholar]
 Nieppola, E., Tornikoski, M., & Valtaoja, E. 2006, A&A, 445, 441 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ong, R. A. 2009, ATel, 2301, 1O [NASA ADS] [Google Scholar]
 Sauge, L., & Henri, G. 2004, ApJ, 616, 136 [NASA ADS] [CrossRef] [Google Scholar]
 Schlickeiser, R. 1984, A&A, 136, 227 [NASA ADS] [Google Scholar]
 Tavecchio, F., Maraschi, L., & Ghisellini, G. 1998, ApJ, 509, 608 [NASA ADS] [CrossRef] [Google Scholar]
 Tavecchio, F., Ghisellini, G., Ghirlanda, G., et al. 2009, MNRAS, 399, 59 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Approximation of synchrotron emission
Fig. A.1 An example of an approximation of SSC emission. The figure shows a broken powerlaw (γ_{b}) electron spectrum with a low energy cutoff (γ_{m}) that extends up to some maximal energy (γ_{x}). This distribution generates synchrotron emission with two breaks () and extends from some minimal () to some maximal frequency (). The energy spectrum is divided into two parts, whereas the synchrotron radiation is divided into three parts. This provides six different components of the inverseCompton emission. Each component must be integrated with different boundary conditions. The upper right panel demonstrates how to choose these conditions. This panel shows synchrotron frequencies as functions of the IC photon frequencies, where the horizontal lines indicate the characteristic frequencies in the synchrotron spectrum. Three slanted lines are calculated according to the wellknow relation for the three characteristic particle energies (γ_{m},γ_{b},γ_{x}). The single opposite slanted line shows the KleinNishina restriction according to the approximation . The crossing lines produce six rectangle areas that indicate the lower and upper integration boundaries for each component. Each component of the IC emission is indicated by a pair of numbers where the first number refers to the part of the electron energy spectrum and the second number shows the part of the radiation field used to create this component. 

Open with DEXTER 
For a powerlaw particle energy distribution (A.1)the synchrotron emissivity can be approximated by a simple powerlaw function (A.2)where α = (n − 1)/2 and (A.3)(e.g., Ginzburg & Syrovatskii 1965). Assuming an optically thin source, we can write the synchrotron luminosity L_{s} = 4πVj_{s} [erg], where V is the source volume. If our source is spherical then the intensity of the surface emission is given by (A.4)where R is the source radius. This gives L_{s} = 4π^{2}R^{2}I_{s} and using the standard luminositytoflux conversion, we can obtain the observed flux F_{s} = π(R/D_{L})^{2}I_{s} [erg cm^{2}], where D_{L} is the luminosity distance. For relativistically moving sources at cosmological distances, we have to apply additional transformations to the flux (A.5)as well as the frequency ν = ν′δ/(1 + z).
Appendix B: Approximation of inverseCompton emission
The inverseCompton emissivity can be approximated by a simple formula (B.1)where x′ = hν′/m_{e}c^{2} is the inverseCompton photon energy normalized to the electron rest energy and is the energy of a synchrotron photon also divided by m_{e}c^{2}. The particle energy spectrum, in our particular estimations, is assumed to be a broken powerlaw (B.2)where K_{2} = K_{1}γ^{n2 − n1} and γ_{m} ≫ 1. For such an energy spectrum, we can approximate the synchrotron intensity by a double broken powerlaw function (B.3)that extends from some minimal energy () up to the maximal energy () and contains breaks at the characteristic energies () related to the particle energy by (B.4)According to Eqs. (A.2) and (A.4), the normalizing coefficients in the synchrotron intensity are given by \arraycolsep1.75ptNote that inside the source the intensity of the synchrotron radiation field is on average I_{s} ≃ (3/4)Rj_{s} (e.g., Kataoka et al. 1999). Therefore, there is no constant 4/3 in the above formulae.
Since the energy spectrum (Eq. (B.2)) and the radiation field spectrum (Eq. (B.3)) are divided into the simple powerlaw functions, we can split the inverseCompton emissivity into six components where each component is described by (B.7)where index a (equal either to 1 or 2) describes parts of the particle energy spectrum and index b (equal 0, 1 or 2) indicates parts of the synchrotron spectrum. Depending on the a,b values, this gives two different results. For a = b, we have (B.8)and for a ≠ b(B.9)Finally, we have to define the integration boundaries (x_{1},x_{2}), which are different for each component. In general, we can define the lower integration boundary as a maximum of two energies of synchrotron photons (B.10)where is equal to , or for b = 0,1,2, respectively, and is calculated for the characteristic particle energies γ_{b},γ_{x} for a = 1 and 2, respectively. The upper boundary condition is calculated to be the minimum of three energies (B.11)where is equal to either , or for b = 0,1,2, respectively and is calculated for γ_{m},γ_{b} for a = 1 and 2, respectively. The last part of this maximum comes from the simple approximation of the KleinNishina limit where (B.12)An example spectrum of the approximate IC emission and a description of the integration boundaries and the components
of the spectrum is presented in Fig. A.1. This particular spectrum was calculated for a set of the parameters that gives all six components of the spectrum. Note that for the parameters that are typical for TeV blazars, the IC scattering occurs predominantly in the KN regime and only two components (j_{(1,1)}) and (j_{(2,1)}) creates the spectrum if 1 ≪ γ_{m} ≪ γ_{b}. In the case where 1 ≪ γ_{m} ≃ γ_{b}, the first part of the energy spectrum (a = 1) and the second part of the radiation field (b = 1) basically do not exist. In such a case, the spectrum is dominated by the j_{(2,0)} component with the spectral index in the MeV − GeV range α_{0} = −1/3. However, this is an extreme limiting case, never directly observed.
All Tables
All Figures
Fig. 1 Approximate synchrotron and IC spectra calculated for different values of γ_{m}, from almost unity up to γ_{b}. This test demonstrates how the spectral index in the MeV − GeV range depends on the γ_{m} value. The dashed line shows the difference between the estimated (α_{1} = 0.5) and calculated (α ≃ 0.65) spectral slope below the peak. 

Open with DEXTER  
In the text 
Fig. 2 The upper panel shows the high energy emission of RGB J0710+591 and approximated spectra used for our estimations. In the lower panel, we illustrate the constraints to the basic parameters obtained in four different ways and the values selected for the modelling. 

Open with DEXTER  
In the text 
Fig. 3 The upper panel shows the multifrequency emission of 1ES 0502+675 observed by the Swift & Fermi instruments (Abdo et al. 2010b – squares; Abdo et al. 2010c – circles) and the approximate spectra used for the estimation presented in the lower panel. The values of B and δ selected for our modelling are indicated by the dot. 

Open with DEXTER  
In the text 
Fig. 4 The modelling of the RGB J0710+591 spectral emission. The upper panel shows the first approach where we considered the Swift and Fermi observations only. In the lower panel, we demonstrate an alternative approach that can explain the origin of the TeV emission observed by VERITAS. Note that the Fermi and the VERITAS observations were not simultaneous. To obtain the observed IC spectra, we used the EBL lowerlimit model proposed by Kneiske & Dole (2010), where the dashed lines show the intrinsic unabsorbed emission. 

Open with DEXTER  
In the text 
Fig. 5 The spectra obtained from the modelling of 1ES 0502+675 for the physical parameters obtained from our estimation. 

Open with DEXTER  
In the text 
Fig. A.1 An example of an approximation of SSC emission. The figure shows a broken powerlaw (γ_{b}) electron spectrum with a low energy cutoff (γ_{m}) that extends up to some maximal energy (γ_{x}). This distribution generates synchrotron emission with two breaks () and extends from some minimal () to some maximal frequency (). The energy spectrum is divided into two parts, whereas the synchrotron radiation is divided into three parts. This provides six different components of the inverseCompton emission. Each component must be integrated with different boundary conditions. The upper right panel demonstrates how to choose these conditions. This panel shows synchrotron frequencies as functions of the IC photon frequencies, where the horizontal lines indicate the characteristic frequencies in the synchrotron spectrum. Three slanted lines are calculated according to the wellknow relation for the three characteristic particle energies (γ_{m},γ_{b},γ_{x}). The single opposite slanted line shows the KleinNishina restriction according to the approximation . The crossing lines produce six rectangle areas that indicate the lower and upper integration boundaries for each component. Each component of the IC emission is indicated by a pair of numbers where the first number refers to the part of the electron energy spectrum and the second number shows the part of the radiation field used to create this component. 

Open with DEXTER  
In the text 