Issue 
A&A
Volume 557, September 2013



Article Number  A48  
Number of page(s)  13  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201321874  
Published online  28 August 2013 
Spontaneously quenched γray spectra from compact sources
Department of PhysicsUniversity of Athens,
Panepistimiopolis, 15783
Zografos,
Greece
email: maroulaaki@gmail.com
Received:
10
May
2013
Accepted:
5
July
2013
Aims. We have studied a mechanism for producing intrinsic broken powerlaw γray spectra in compact sources. This is based on the principles of automatic photon quenching, according to which γrays are being absorbed on spontaneously produced soft photons whenever the injected luminosity in γrays lies above a certain critical value.
Methods. We derived an analytical expression for the critical γray compactness in the case of powerlaw injection. For the case where automatic photon quenching is relevant, we calculated analytically the emergent steadystate γray spectra. We also performed numerical calculations in order to back up our analytical results.
Results. We show that a spontaneously quenched powerlaw γray spectrum obtains a photon index 3Γ/2, where Γ is the photon index of the powerlaw at injection. Thus, large spectral breaks of the γray photon spectrum, e.g. ΔΓ ≳ 1, can be obtained by this mechanism. We also discuss additional features of this mechanism that can be tested observationally. Finally, we fit the multiwavelength spectrum of a newly discovered blazar (PKS 0447439) by using such parameters to explain the break in the γray spectrum by means of spontaneous photon quenching, under the assumption that its redshift lies in the range 0.1 < z < 0.24.
Key words: radiation mechanisms: nonthermal / gamma rays: general / BL Lacertae objects: general
© ESO, 2013
1. Introduction
The production and radiation transfer of highenergy γrays is a physical problem that has received a lot of attention over the last forty years since it can be applied on compact highenergy emitting astrophysical sources, such as active galactic nuclei (AGN) and gammaray bursts. Photonphoton absorption, in particular, turns out to be a significant physical process in compact Xray and γray emitting sources that results in electromagnetic (EM) cascades (e.g. Jelley 1966; Herterich 1974). The effects of EM cascades can be studied within either a linear or nonlinear framework. In the first approach (e.g. Protheroe 1986), the number density of target photons is assumed to be fixed, whereas in the second, the produced electron/positron pairs produce photons, which in their turn serve as targets for photonphoton absorption (Kazanas 1984; Zdziarski & Lightman 1985; Svensson 1987). The first analytical studies of EM cascades were then followed by numerical works, which aimed at computing timedependent solutions to the kinetic equations of electrons and photons taking photonphoton annihilation into account (Coppi 1992; Mastichiadis & Kirk 1995; Stern et al. 1995; Böttcher & Chiang 2002). These algorithms are now commonly used in source modelling (Mastichiadis & Kirk 1997; Kataoka et al. 2000; Konopelko et al. 2003; Katarzyński et al. 2005).
Intrinsically nonlinear effects in EM cascades initiated by photonphoton absorption, however, have only recently gained attention. First, Stawarz & Kirk (2007, hereafter SK07) studied the case where no soft (i.e. target) photons are present in a source. They investigated the necessary conditions under which γray photons can cause runaway pair production and found that these conditions can be summarized only in a single quantity, the critical γray compactness. This can be considered an upper limit of the γray compactness that depends on parameters such as the magnetic field strength and the size of the source. If the injection compactness of very high energy (VHE) photons (≳0.1 TeV) is larger than the critical compactness, the following nonlinear loop is selfsustained: γray photons are absorbed on soft photons emitted by the produced pairs through synchrotron radiation.
The work of SK07 was then expanded by Petropoulou & Mastichiadis (2011, hereafter PM11) mainly by taking into account continuous energy losses of the produced pairs. The nonlinear loop of processes called “automatic photon quenching” by SK07 can be the core of other more complex ones. In particular, Petropoulou & Mastichiadis (2012b) or PM12b from this point on, have attributed the production of VHE γrays to synchrotron emission from secondaries produced in charged pion decay, while pions were the result of photohadronic interactions between relativistic protons and soft photons. Petropoulou & Mastichiadis (2012b) have shown that the system of protons and photons resembles a preypredator system whenever automatic photon quenching operates, and it shows interesting variability patterns, such as limit cycles^{1}.
In the present work we continue the exploration of spontaneous photon quenching by studying the case of powerlaw γray injection in the source and, in that sense, it can be considered a continuation of the aforementioned works. There were two main motivations to our present study: (i) γray spectra emitted by a powerlaw distribution of relativistic particles through some radiation mechanism, e.g. synchrotron radiation and inverse compton scattering, can be modelled by a powerlaw, at least partially; and (ii) if spontaneous photon absorption affects part of the γray injection spectrum, spectral breaks are produced; we believe that this requires further investigation as it is an intrinsic mechanism for producing breaks in a γray spectrum and it could be of relevance to recent results regarding γray emitting blazars.
The present paper is structured as follows: in Sect. 2 we derive an analytical expression for the critical γray compactness in the case of powerlaw injection using certain simplifying assumptions, while we comment also on the validity range of our result. In the case where spontaneous photon quenching becomes relevant, we show that a break at the steadystate γray spectrum appears and we also calculate analytically the expected spectral change. In Sect. 3 we derive numerically the critical compactness for a wide range of parameter values and examine the effects that a primary soft photon component in the source would have on γray absorption. Possible implications of spontaneous absorption on γray emitting blazars are presented in Sect. 4. We also present a list of observationally tested characteristics that a spontaneously quenched source would, in principle, show. In the same section we also show that the spectral energy distribution (SED) of the newly discovered blazar PKS 0477439 can be explained within the framework of automatic quenching. For the required transformations between the reference systems of the blazar and the observer we have adopted a cosmology with Ω_{m} = 0.3, Ω_{Λ} = 0.7, and H_{0} = 70 km s^{1} Mpc^{1}.
2. Analytical approach
We consider a spherical region of radius R containing a magnetic field of strength B. We assume that γrays are being produced in this volume by some nonthermal emission process, e.g. proton synchrotron radiation. In our analysis however, the γray production mechanism remains unspecified, since its exact nature does not play a role in the derivation of our results. Furthermore, γrays are being injected with a luminosity that is related to the injected γray compactness as (1)where σ_{T} is the Thomson cross section. Without any substantial soft photon population inside the source, the γrays will escape without any attenuation in one crossing time. However, as SK07 and PM11 showed, the injected γray compactness cannot become arbitrarily high because, if a critical value is reached, the following loop starts operating: (i) gammarays pairproduce on soft photons, which can be initially arbitrarily low inside the source; (ii) the produced electronpositron pairs are highly relativistic since they are created with approximately half the energy of the initial γray photon, and cool mainly by synchrotron radiation, thus acting as a source of soft photons; and (iii) the emitted synchrotron photons have lower energy when compared to the γray photons and serve as targets for more γγ interactions.
There are two conditions that should be satisfied simultaneously for this network to occur:

1.
Feedback criterion This is related to the energythreshold condition for photonphoton absorption and it requiresthat the synchrotron photons emitted from the pairs havesufficient energy to pairproduce on the γrays.

2.
Marginal stability criterion This is related to the optical depth for photonphoton absorption, which must be above unity in order to establish the growth of the instability.
By making suitable simplifying assumptions, one can derive an analytic relation for the first condition (see also SK07). Thus, combining (i) the threshold condition for γγ absorption ϵx = 2, where ϵ and x are the γray and soft photon energies in units of m_{e}c^{2}, respectively (this normalization will be used for all photon energies throughout the text unless stated otherwise); (ii) the fact that there is equipartition of energy among the created electronpositron pairs γ_{e} = γ = ϵ/2; and (iii) the assumption that the required soft photons are the synchrotron photons that the electrons/positrons radiate, i.e. x_{s} = bγ^{2}, where b = B/B_{crit} and G, one finds the relation (2)that defines, for a certain magnetic field strength, the γray photon energy above which automatic photon quenching becomes relevant.
In what follows, we will concentrate on the second condition, since our aim is to determine the value of the injected γray compactness that ensures the growth of the instability. The corresponding calculations in the case of monoenergetic γray injection can be found in SK07 and PM11. Here we focus on the more astrophysically relevant case of a powerlaw γray injection. To treat this problem analytically we discretize the powerlaw of γrays. In particular, we begin by calculating the critical compactness in the case where two monoenergetic γrays are injected. We repeat the calculation for the injection of three monoenergetic γrays and, finally, we generalize our result for N monoenergetic injection functions. Furthermore, our treatment is built upon the following assumptions and approximations:

1.
Only two physical processes are taken into account, i.e.photonphoton absorption and synchrotron radiation of theproduced pairs. Inverse Compton scattering of pairs on thesynchrotron produced photons can be safely neglected becauseof the strong magnetic field that is typically required for theautomatic photon quenching loop to function, and because of thelarge^{2}Lorentz factors of the produced pairs.

2.
Only the equations describing the evolution of γrays and synchrotron (soft) photons are taken into account. The equation for the pairs is neglected, since these have synchrotron cooling timescales much smaller than the crossing time of the source. Thus, all the injected energy into pairs is transformed into synchrotron radiation.

3.
Synchrotron emissivity is approximated by a δfunction, i.e. j_{s}(x) = j_{0}δ(x − x_{s}), where x_{s} = bγ^{2} is the synchrotron critical energy.

4.
The synchrotron energy losses of pairs are treated as catastrophic escape from the considered energy range. In other words, an electron with Lorentz factor γ loses its energy by radiating synchrotron photons at energy bϵ^{2}/4.

5.
The cross section of photonphoton absorption (in units of σ_{T}) is approximated as (3)which is the same as the cross section given by Coppi & Blandford (1990) apart from the logarithmic term ln(x).
2.1. Marginal stability criterion for injection of two monoenergetic γrays
If we assume that γrays are being injected into the source at energies ϵ_{1} and ϵ_{2} (ϵ_{1} < ϵ_{2}) with compactnesses , where i = 1,2, we can use suitable parameters to ensure that the γray photon with the minimum energy satisfies the feedback criterion, i.e. ϵ_{1} > ϵ_{q} = 2/b^{1/3}. Then, all higher energy γray photons also satisfy the feedback criterion and the corresponding emitted synchrotron photons have energies . Gammaray photons with energy ϵ_{i} can, therefore, be absorbed on both soft photon distributions because the energy threshold criterion is satisfied for all the four possible photonphoton interactions, i.e. ϵ_{i}x_{s,j} > 2 for i,j = 1,2. We note also that we do not consider absorption of γrays on less energetic γrays. Moreover, the number densities of γray photons and of the corresponding soft photons are denoted as and , respectively, where the symbols imply the relations and ; the densities refer to the number of photons contained in a volume element σ_{T}R. In other words, if expresses the number of photons per erg per cm^{3}, then . The dimensionless photon number densities are also related to the compactnesses through the relation (4)for discrete monoenergetic injection^{3}. Using the notations introduced above along with the assumptions 1–5 the system can be described by four equations where time is normalized with respect to the photon crossing/escape time from the source, i.e. τ = ct/R, and the operators L and Q denote losses and injection respectively; the loss term in Eq. (6) due to photonphoton absorption is omitted since it is negligible. We note that number densities and rates are equivalent in the dimensionless form of the above equations. The explicit expressions of the operators in Eqs. (5) and (6) are where the normalization constant is calculated by equating the total γray energy loss rate with the total energy injection rate into soft photons, i.e. − ^{∫}dϵ ϵL_{γγ} = ^{∫}dx xQ_{γγ}, and it is given by (10)We also note that in the case of continuous powerlaw injection, Eq. (9) should be replaced by Q_{inj}(ϵ) = 3ℓ_{inj}(ϵ)/ϵ. The trivial stationary solution of Eqs. (5) and (6) is , where and it corresponds to the case where the injection rate of γrays equals the photon escape rate from the source. Following the methodology described in SK07 and PM11, we introduce perturbations to all photon number densities and linearize the set of Eqs. (5) and (6) around the trivial solution. The linearized system can be written in the form dY/dτ = AY where (11)is the vector of the perturbed number densities and A is the matrix of the linearized system of equations (12)In order to build a finite number of soft photons in the source, the perturbations must grow with time. This is ensured if at least one of the eigenvalues of matrix A is positive. After some algebraic manipulation we find that, indeed, one eigenvalue can become positive if the following condition holds: (13)The same methodology can be applied in the case where γrays are being injected as a δ function at three energies ϵ_{1} < ϵ_{2} < ϵ_{3}. This leads to an analogous critical condition (14)
2.2. Critical compactness for powerlaw γray injection
The above can be generalized for the case of N monoenergetic γrays with energies ϵ_{1} < ϵ_{i} < ϵ_{N} to find the marginal stability criterion (15)We have verified that the above relation also applies to cases where the feedback criterion is not satisfied for γrays of the minimum energy but for higher energy photons, i.e. ϵ_{k} = 2/b^{1/3} with k > 1, with only a slight modification: the summation starts from i = k. Moreover, if we assume that the injection rate of γray photons can be modelled as a powerlaw, e.g. (16)the criterion of Eq. (15) takes the form (17)If N → ∞ and (ϵ_{i + 1}−ϵ_{i}) → 0, the discrete sum of Eq. (17) can be trasformed into an integral which leads to (18)Finally, if we replace the normalization constant Q_{0} by the integrated γray compactness over all photon energies (see also comment on Eq. (9)) using (19)we find that ℓ_{inj} ≥ ℓ_{γ,cr} where ℓ_{γ,cr} is the critical compactness for a powerlaw γray injection. This has a compact form (20)where ϵ_{min} ≡ ϵ_{1}, ϵ_{max} ≡ ϵ_{N}, and ϵ_{M} = max [ ϵ_{min},ϵ_{q} ] (see Eq. (2) for the definition of ϵ_{q}).
We were able to derive an analytical and rather simple expression of the critical compactness in the case of powerlaw injection at the cost, however, of a series of approximations/assumptions that may limit the validity range of our result. It is reasonable therefore, before closing the present section to check the range of validity of Eq. (20). For this, we made a comparison between this expression and the numerically derived values^{4}, which is illustrated in Fig. 1. Both panels show the dependence of ℓ_{γ,cr} on ϵ_{min} for two different pairs of photon indices marked on the plots. Lines and symbols are used for plotting the analytically and numerically derived values respectively, while different types of lines/symbols correspond to different values of the photon index. For values of ϵ_{min} below ϵ_{q}, which for the values used in this example equals 2 × 10^{4}, we find a good agreement between our analytical and numerical results, apart from a numerical factor of ≃2; we note that the dependence of ℓ_{γ,cr} given by Eq. (20) on ϵ_{min} coincides with the numerically determined dependence. However, the analytical solution of ℓ_{γ,cr} fails in the energy range of ϵ_{min} ≥ ϵ_{q}, since in this regime the approximation 4 listed at the beginning of Sect. 2 proves to be crude.
Fig. 1
Top panel: critical compactness ℓ_{γ,cr} as a function of the minimum energy of the γray spectrum ϵ_{min} for Γ = 1.6 (solid line) and Γ = 2 (dashed line). The numerically derived values for the two cases are shown with circles and triangles, respectively. Bottom panel: same as in top panel except for different photon indices, which are marked on the plot. Other parameters used are: ϵ_{max} = 2.3 × 10^{5} (in m_{e}c^{2} units), B = 40 G, and R = 3 × 10^{16} cm. In both panels the grey area denotes the region where ϵ_{min} > ϵ_{q}; see text for the definition of ϵ_{q}. 
2.3. Steady state γray quenched spectrum
Assuming that the injected γray compactness exceeds the critical value derived in the previous section, we search for steadystate solutions of the system. There is a main difference between the analytical approach that we follow in this section and the one presented in Sect. 2.1: here we take into account the cooling of pairs due to synchrotron radiation, i.e. synchrotron losses are treated as a continuous energy loss mechanism. For this reason, in addition to the equations of γrays and soft photons, we include in our analysis a third equation for electron/positron pairs or simply electrons from this point on. Dropping the time derivatives, the kinetic equations for γrays, soft photons, and electrons are written respectively as and (23)where n_{e} is the dimensionless electron number density, and an electron escape timescale t_{e,esc} = t_{cr} was assumed. All normalizations and approximations, apart from the fourth in our list, are the same as in Sect. 2.1. The synchrotron emissivity is approximated by a δfunction (see Sect. 2.1). We note that in what follows we can safely neglect electron escape from the source since for magnetic field strengths relevant to automatic quenching, synchrotron cooling is the dominant term in the electron equation; thus, the lefthand side of Eq. (23) is essentially equal to zero. Then, the injection and loss operators take the forms where the form of implies the use of a δfunction for the synchrotron emissivity (see e.g. Mastichiadis & Kirk 1995) and α_{1} = (2/3)ℓ_{B}b^{− 3/2}. The loss operators are given by where is the maximum energy of the produced soft photon distribution^{5}, α_{2} = 4ℓ_{B}/3 with ℓ_{B} = σ_{T}RU_{B}/m_{e}c^{2} being the magnetic compactness, and σ_{γγ} is the dimensionless (in units of σ_{T}) cross section for photonphoton absorption (see point (5) of Sect. 2). Within this approximation Eq. (27) takes the form (29)which further implies that γrays with energy ϵ < ϵ_{br} ≡ 2/x_{max} cannot be absorbed. Thus, pairs with γ < γ_{br} ≡ ϵ_{br}/2 cannot be produced, i.e. the injection term in Eq. (23) vanishes. The above can be summarized by inserting the step function H(ϵ − ϵ_{br}) in the expression of given by Eq. (29). On the one hand, for γ < γ_{br} the electron distribution has the trivial form n_{e} ∝ γ^{2}. On the other hand, for γ ≥ γ_{br}, the distribution is determined by synchrotron losses and pair injection.
We now proceed to calculate the electron distribution for γ > γ_{br}. After having inserted the above expressions for the operators into Eqs. (21)–(23), we combine them into one nonlinear integrodifferential equation, where the unknown function is n_{e}(γ): (30)where (31)Integration of Eq. (30) leads to (32)where we have used the notation and . If γ_{br} > γ_{min} then , whereas if γ_{br} < γ_{min} then can be equal either to γ_{min} or γ. The condition γ_{br} < γ_{min} corresponds to the physical case where the entire γray powerlaw spectrum is affected by automatic quenching. Therefore, the steadystate γray spectrum will show, in general, no break (see also numerical example in Fig. 8 of PM11). Since in the present work we are interested in producing broken powerlaw γray spectra, we will examine only the case where γ_{br} > γ_{min} and therefore . Since an exact solution of Eq. (32) cannot be obtained analytically, we assume a certain form for the solution we are searching for. A reasonable guess is that the electron distribution is a powerlaw, i.e. n_{e} = N_{0}γ^{− p}, where N_{0} and p are the parameters to be defined. Inserting this function into Eq. (32) we obtain (33)where and A(N_{0},p) = N_{0}b^{p/2}α_{1}σ_{0}/(p + 1) is a function only of the unknown parameters for fixed values of the magnetic field and source size. We also note that we have neglected the contribution of terms calculated at the upper limit (x_{max}) while performing the integral I(γ′), so that the solution of the final integral appearing at the righthand side of Eq. (33) can be expressed in closed form as (34)where 2F1(a,b,c;z) is the Gauss hypergeometric function. Without making any further approximations, the above equation cannot be solved for p and N_{0} analytically. However, if the term Aγ^{′β} inside the integral is much larger than unity, this is, then, greatly simplified and Eq. (33) takes the form (35)which is satisfied for The above solution is valid as long as Γ ≠ 1. In the opposite case, one should search for more general forms of the electron distribution, i.e. n_{e}(γ) ∝ N_{0}(γ)γ^{− p}, and follow the same procedure. However, for the purposes of the present work, the solution for Γ ≠ 1 is sufficient. Summarizing, the electron distribution is given by (38)where we have demanded the solution to be continuous at γ = γ_{br}. We note that the slope of the electron distribution above γ_{br} is Γ + 1, i.e. steeper than that at injection (see term γ^{′−Γ} in Eq. (32)) by one. This denotes the efficient synchrotron cooling of the whole distribution. We emphasize that the results for γ ≳ γ_{br} should be taken cautiously into consideration, since the above solution is not valid for all γ above γ_{br}; we note that it was derived under the approximation Aγ^{β} > 1 or equivalently (39)After having determined the form of the electron distribution we can then calculate the steadystate solutions for γrays and soft photons. For ϵ ≤ ϵ_{br} the injected spectrum remains unaffected by automatic quenching, whereas for ϵ > ϵ_{br} the solution n_{γ}(ϵ) can be found by inserting first the second branch of n_{e} into Eq. (22) and then by using the derived expression for n_{s}(x) in Eq. (21). The results are summarized below (40)where α = Γ/2 + 1 and (41)As the solution for n_{e} is not valid for γ ≳ γ_{br} (see Eq. (39)), the behaviour of n_{γ} close to the transition, i.e. for ϵ ≳ ϵ_{br}, should also be considered with caution. The validity range set by Eq. (39) can be translated into terms of photon energies, i.e. ϵ > ϵ_{⋆}, where (42)It can be easily verified that for ϵ > ϵ_{⋆} the second term in the denominator of n_{γ} is larger than unity^{6}, which simplifies the functional dependence of n_{γ} on energy, i.e. n_{γ} ∝ ϵ^{− 3Γ/2}. Summarizing, we find that the asymptotic photon index of a spontaneously quenched γray spectrum is well defined and it is given by 3Γ/2; this result is in complete agreement with the one derived numerically in PM11 (see Fig. 9 therein). The spectral break in the case of automatic quenching depends, therefore, linearly on the photon index at injection, i.e. ΔΓ = Γ/2.
Figure 2 shows the γray spectra given by Eq. (40) for Γ = 2 along with ϵ_{⋆} (crosses) for different values of the injection parameter Q_{0} that ensure the operation of spontaneous photon quenching. Our solution for ϵ > ϵ_{br} becomes progressively valid over a larger range of energies as Q_{0} increases. Moreover, for the highest value of Q_{0}, the photon index of the quenched part of the spectrum is 3, i.e. it has obtained its asymptotic value defined by 3Γ/2. The analytical solutions presented in Fig. 2 are compared to those derived using the numerical code described in the following section and PM11. Solid and dashed lines in Fig. 3 correspond to the analytical and numerical solutions respectively. The agreement between the two is better for larger values of the injection compactness, i.e. when the absorption term in the equation of γray photons becomes larger.
Before closing the present section and for reasons of completeness we will make a short comment on our choice of assuming continuous instead of catastrophic energy losses of electrons. We have also derived the steadystate solution in the case of monoenergetic injection of Nγrays using catastrophic losses. However, the obtained results, when compared to the numerically derived results, were not reasonable, in the sense that no large spectral breaks were produced because γray absorption was underestimated. The main reason behind this disagreement with the numerically derived results is the approximation 4. By assuming catastrophic losses of the produced pairs we neglect soft photons with x < bϵ^{2}/4, i.e. we artificially decrease the optical depth for absorption of a γray photon with certain energy ϵ. Thus, although the assumption of catastrophic losses proves to be suffiecient for the derivation of ℓ_{γ,cr} (Sect. 2.1) it proves to be too crude for more quantitative results.
3. Numerical approach
In this section we will present

(i)
the dependence of the critical injection compactness onvarious parameters, e.g. on the minimum and maximum energyof injected γrays, as well as on on the γray photon index, for a widerange of values.

(ii)
the effects that the presence of lowenergy photons has on the automatic absorption of γrays.
For a detailed study of the above a numerical treatment is required; as far as the first point is concerned, we have already shown that the analytical approach breaks down (for example, see Sect. 2.1 for sufficiently high values of the minimum energy of injected γrays).
To numerically investigate the properties of quenching one needs to solve again the system of Eqs. (5) and (6), where the discretized photon number densities should be replaced by their continuous functional form. For completeness we have augmented it to include more physical processes.
Fig. 2
Analytical solution of steadystate spontaneously quenched γray spectra for different values of the injection parameter Q_{0} starting from 1.25 × 10^{3}. Each value is increased by a factor of 5 over its previous value. Crosses denote in each case the value of ϵ_{⋆}. Other parameters used are: Γ = 2, ϵ_{min} = 33, ϵ_{max} = 2.3 × 10^{5}, B = 40 G, and R = 3 × 10^{16} cm. The break energy is then ϵ_{br} = 160. 
Fig. 3
Numerical (dashed lines) and analytical (solid lines) solutions for the steadystate γray spectra. The injection rate increases from bottom to top starting from Q_{0} = 1.25 × 10^{3}. The parameters used are the same as in Fig. 2. 
As in the numerical code, there is no need to treat the timeevolution of soft photons and γrays through separate equations, the system can be written (43)and (44)where n_{e} and n_{γ} are the differential electron and photon number densities, respectively, normalized as in Sect. 2. Here we considered the following processes: (i) photonphoton pair production that acts as a source term for electrons () and a sink term for photons (); (ii) synchrotron radiation that acts as a loss term for electrons () and a source term for photons (); (iii) synchrotron selfabsorption that acts as a loss term for photons (); and (iv) inverse Compton scattering that acts as a loss term for electrons () and a source term for photons (). In addition to the above, we assume that γrays are injected into the source through the term (). The functional forms of the various rates have been presented elsewhere (Mastichiadis & Kirk 1995, 1997; Petropoulou & Mastichiadis 2009). The photons are assumed to escape the source in one crossing time, therefore t_{γ,esc} = R/c. The electron physical escape timescale from the source t_{e,esc} is another free parameter that, however, is not important in our case. Thus, we will fix it at value t_{e,esc} = t_{γ,esc} = R/c. The final settings are the initial conditions for the electron and photon number densities. Because we are investigating the spontaneous growth of pairs and their emitted synchrotron photons, we assume that at t = 0 there are no electrons in the source, so we set n_{e}(γ,0) = η → 0. Moreover, during the injection of photons in a certain γray energy range it is important to keep the background photons used in the numerical code at a level as low as possible in order to avoid artificial growth of the instability.
3.1. Critical γray compactness for powerlaw injection
The procedure we follow for the numerical determination of the critical compactness is as follows. We start by injecting γrays at a low rate in a specific energy range (ϵ_{min},ϵ_{max}), e.g. ℓ_{inj} ≃ 10^{5}, and then we increase ℓ_{inj} over its previous value by a factor of 0.2 in logarithm. The increase of ℓ_{inj} is directly related to the increase of the normalization Q_{0} of the injection γray spectrum; we note that for Γ ≠ 2. For each value we allow the system to reach a steady state and then we examine the shape of the γray spectrum. We define as ℓ_{γ,cr} that value of ℓ_{inj} that causes the first deviation from the spectral shape at injection; for this reason, we consider it to be a strict limit.
Fig. 4
Log–log diagram of the critical γray compactness as a function of ϵ_{min} for constant ϵ_{max} = 2.3 × 10^{5} (top panel) and as a function of ϵ_{max} for constant ϵ_{min} = 1.4 × 10^{5} (bottom panel) for three photon indices marked on the plot. Other parameters used are B = 40 G and R = 3 × 10^{16} cm. 
Figure 4 shows ℓ_{γ,cr} as a function of ϵ_{min} (top panel) or ϵ_{max} (bottom panel) for three different slopes of the injection spectrum marked on the plot. Other parameters used are B = 40 G and R = 3 × 10^{16} cm. In the top panel, the maximum energy of the powerlaw spectrum is fixed at ϵ_{max} = 2.3 × 10^{5}, whereas in the bottom panel the minimum energy is taken to be constant and equal to ϵ_{min} = 1.4 × 10^{5}. For soft injection spectra, e.g. Γ = 2.4, we find that ℓ_{γ,cr} strongly depends on ϵ_{min}. The situation is exactly the opposite for hard γray spectra (see solid line in bottom panel of the same figure). The critical compactness that we derive for ϵ_{min} ≈ ϵ_{max} ≈ 2.5 × 10^{5} is ℓ_{γ,cr} ≈ 2 × 10^{3} (see bottom panel of Fig. 4) and it corresponds to monoenergetic γray injection. Thus, it should be compared to ℓ_{γ,cr} ≈ 5 × 10^{4}, which was derived analytically in PM11 for δfunction injection at ϵ = 2.5 × 10^{5} (see Fig. 2 therein). The difference between the two results is not worrying, since the analytical values in PM11 were about a factor of four lower than the accurate values that were derived numerically.
3.2. Effects of a primary soft photon component
The main difference between automatic γray absorption and the widely used photonphoton absorption on a preexisting photon field (primary photons), is that in the first case no target field is initially present in the source. It is as if the system finds its own equilibrium by selfproducing the soft photons required for quenching the extra γray luminosity. In many physical scenarios, however, primary photons are present in the source, e.g. synchrotron radiation from primary electrons, and therefore γrays are more likely to be absorbed on both primary and secondary photon fields. Here arises the question whether or not the effects of spontaneous photon quenching can be disentagled from those of linear absorption^{7}.
In the limiting cases where γrays are being absorbed on either primary or secondary soft photons there is a straighforward relation between the photon index of the absorbed γray spectrum and that of photon target field where s is the photon index of the primary photon distribution (see Appendix A for the derivation of Γ_{abs} in the first case).
Using the numerical code presented in Sect. 3 we will first verify the above analytical relations and then study intermediate cases. For this, we assume a spherical region with size R = 3 × 10^{16} cm containing a magnetic field B = 40 G. Very high energy γrays with a photon index Γ = 2.4 are injected into this volume between energies ϵ_{min} = 23 and ϵ_{max} = 2.3 × 10^{5} with compactness ℓ_{inj} being a free parameter. Primary photons are produced via synchrotron radiation of a powerlaw electron distribution with index p = 1.5. The maximum Lorentz factor γ_{max} and the electron injection compactness , which is defined as , are also free parameters; here is the total electron injection luminosity. The only processes we consider in the following examples are synchrotron emission, photonphoton absorption, and escape from the source. We note that the effect of inverse Compton scattering is negligible for the magnetic field and electron energies assumed here. The results regarding the two regimes are summarized in Fig. 5. In the top panel we have used (solid line) and (dashed line), while we kept γ_{max} = 3.6 × 10^{5} fixed. For these parameter values the maximum energies of primary and secondary synchrotron photons are 10^{1} and 1.3 × 10^{2}, respectively. The slopes of different powerlaw segments in a x^{2}n_{γ}(x) plot are also shown. The spectral break of the γray spectrum differs between the two regimes and the numerical results are in agreement with those given by Eqs. (45) and (46). In particular, we find that the absorbed γray spectrum has a photon index Γ_{abs} = 3.5 which should be compared to 3Γ/2 = 3.6 for the nonlinear quenching case. To estimate the spectral change for the linear absorption case the photon index of the soft photon distribution is required, which in this example is s = (p/2) + 1 = 1.75. The expected value of Γ_{abs} is 3.15 while the derived value is 3. We also note that the spectral shape of the synchrotron component emitted by the produced pairs (solid line) is in agreement with that expected from our solution for the electron distribution (see Eq. (38)); emission from lower energy electrons (n_{e} ∝ γ^{2}) results in n_{γ} ∝ x^{− 3/2}, while from higher energy electrons (n_{e} ∝ γ^{−Γ−1}) corresponds to n_{s} ∝ x^{−Γ/2−1} = x^{2.2}. In the bottom panel, along with the examples shown in the top panel (solid and dashdotted lines), we plot two intermediate cases with a progressively higher γray compactness, i.e. ℓ_{inj} = 8 × 10^{2} (dashed line) and ℓ_{inj} = 8 × 10^{1} (dotted line), while we used the same . Already from ℓ_{inj} = 8 × 10^{2}, which is still above the critical value, the contribution of photons produced via automatic γray quenching is evident as a bump in the primary soft photon component. We note also that the presence of primary soft photons enhances the automatic absorption of γrays for the same ℓ_{inj} (compare the dotted and dashdotted lines).
Fig. 5
Top panel: comparison between two limiting cases where γray photons are absorbed only on synchrotron photons emitted by secondary pairs (solid line) or by primary electrons (dashed line). The injection compactnesses of γrays and primary electrons are (solid line) and (dashed line), respectively. Different parts of the spectra have different powerlaw dependencies that are marked on the plot. Bottom panel: example of an indermediate case where γrays are being partially absorbed on synchrotron photons from secondaries. The injection compactness of primary electrons is , while ℓ_{inj} takes the following values: 8 × 10^{3} (solid line), 8 × 10^{2} (dashed line), and 8 × 10^{1} (dotted line). The dashdotted line that is obtained for (same as solid line in top panel) is plotted for comparison reasons. Other parameters used for the plot are: B = 40 G and R = 3 × 10^{16} cm. 
4. Implications on γray emitting blazars
4.1. General remarks
The mechanism of automatic photon quenching sets an upper limit to the intrinsic γray luminosity of a compact source and, therefore, can be applied on γray emitting blazars for constraining physical quantities of their VHE emission region (for more details see PM11 and Petropoulou & Mastichiadis 2012a, hereafter PM12a). This can be relevant to recent observations that have revealed the presence of very hard intrinsic TeV γray spectra of blazars, even when corrected with low extragalactic background light (EBL) flux levels, e.g. 1ES 1101232 (Aharonian et al. 2006) and 1ES 0229+200 (Aharonian et al. 2007). The SEDs of some hard γray sources are very difficult to explain within onezone emission models, and therefore, they are often attributed to a second component, whose emission in longer wavelengths is hidden by that of the first component (e.g. Costamante 2012). Thus, the physical conditions of the component emitting in the TeV energy range can be chosen quite arbitrarily only by demanding to fit the VHE part of the spectrum. The presence of very hard TeV γray spectra (typically Γ_{int} ≃ 1.5) implies that these sources cannot be spontaneously quenched, i.e. their intrinsic γray luminosity should be less than the critical one. We note, however, that there is an alternative scenario that employs the photonphoton absorption on internal soft photon fields that explains the formation of very hard intrinsic γray spectra (Aharonian et al. 2008).
By now there are several simultaneous and quasisimultaneous observations of blazars in the GeV and TeV energy range, which clearly show that the γray spectrum cannot be fitted by a simple powerlaw over the whole GeVTeV energy range, but rather by a broken powerlaw. The change of the photon index in many cases is large (ΔΓ ≳ 1) and it cannot be explained using simple arguments such as cooling breaks (ΔΓ = 0.5). In particular, for highpeaked blazars that are bright in the Fermi energy band (e.g. Mrk421, PKS 2155304, PKS 0447439, etc.), the peak of their SED seems to fall in the highGeV energy part of the spectrum (≃100 GeV). In these cases GeV and TeV emission correspond to parts of the spectrum below and above the highenergy hump respectively. Thus, it is commonly considered that the VHE γray spectra are intrinsically much softer (e.g. exponential cutoff effects) than the GeV spectra (Costamante 2012).
Here, however, we investigate another explanation of γray spectral breaks. In our framework, the injection γray spectrum is described by a single powerlaw from GeV up to TeV energies, and the spectrum of the escaping γray radiation is modified due to internal spontaneous photon absorption. We note that the instability of automatic photon quenching offers an alternative mechanism for producing intrinsic broken powerlaw spectra (see Sect. 2).
A plot of the photon index in the GeV energy band (as measured by the Fermi satellite) versus the one in the TeV energy band (as measured by MAGIC and H.E.S.S. telescopes) is shown in Fig. 6. The sources used for this plot are listed in Table 1 along with the observed values of the photon indices and the reference papers. Filled and open symbols show the photon indices of the observed and of the corrected for EBL absorption VHE spectra. We note that in all cases we used the model C by Finke et al. (2010) for the EBL correction; in this model, the EBL flux from UV to the near – IR is also similar to that of Domínguez et al. (2011). Different symbols, in particular circles and squares, denote TeV observations made by MAGIC and H.E.S.S., respectively. The solid line represents the relation between the photon indices that is expected by spontaneous photon quenching, i.e. Γ_{TeV} = 3Γ_{GeV}/2, whereas the dashed line (Γ_{TeV} = Γ_{GeV}) is plotted only to guide the eye. Some of the data points within their error bars are compatible with the theoretical predictions. For the purposes of the present work, the choice of a particular EBL model does not significantly affect our results. For example, the EBL model used here predicts slightly higher optical depth than the one of Franceschini et al. (2008) for E_{γ} ≲ 5 TeV (10 TeV) for z = 0.6 (0.1), respectively. Thus, correction of the VHE spectra with the EBL model of Franceschini et al. (2008) would result in slightly softer intrinsic VHE spectra, i.e. some of the open symbols in Fig. 6 would move upwards in the vertical direction. Moreover, the fact that the analysis of the Fermi and MAGIC/H.E.S.S./VERITAS data has not been made over the same energy intervals for all sources listed in Table 1 makes difficult to draw any conclusions regarding the statistical properties of this sample. Nonsimultaneity of GeV and TeV observations in some cases, e.g. NGC 1275, also makes any coherent comparison difficult. For these reasons, this type of plot could be used only as a first indicator for searching among sources that could be explained by the mechanism of automatic quenching.
Fig. 6
Photon index of the TeV spectrum versus the one of the GeV spectrum for several γray emitting blazars. 
In these cases, there are five additional features that can be tested observationally:

(i)
Automatic photon quenching is a radiative instability thatredistributes the energy within a photon population. Theabsorbed γray luminosity appears, therefore, in the lower part ofthe multiwavelength (MW) spectrum, usually in the Xrayregime.

(ii)
The spectral break in the γray spectrum of blazars that have low Xray emission compared to that of VHE γrays cannot, in general, be attributed to automatic photon absorption.

(iii)
If the VHE γray spectrum is spontaneously absorbed, there is a straighforward relation between the photon indices of the absorbed part of the γray spectrum and that of the soft photon component. In Sect. 2.3 we have shown that Γ_{abs} = 3Γ/2, where Γ is the photon index of the γray spectrum at injection. The steady state electron distribution due to pair production is n_{e}(γ) ∝ γ^{−Γ−1} (see Eq. (38)) and the corresponding photon index of the synchrotron spectrum is given by Γ_{soft} = Γ/2 + 1. Thus, the relation between the two photon indices is Γ_{abs} = 3(Γ_{soft}−1) (see also the numerical example in the top panel of Fig. 5).

(iv)
Strong correlation between the soft component of the MW spectrum and the unabsorbed part of the γray spectrum is to be expected in cases where the intrinsic γray luminosity varies.

(v)
An increase of the intrinsic γray compactness is accompanied by a shift of the break energy towards lower energies.
As far as the first observational prediction is concerned, one can show that the maximum energy of the soft component produced by automatic photon quenching falls for reasonable parameter values, in the Xray regime. First, the feedback criterion for automatic quenching must be satisfied, at least, by γray photons having the maximum energy (see also Sect. 2). This is written as (47)where we have used the observed quantities instead of those measured in the comoving frame of the blob that has a Doppler factor δ. From this point on and in what follows the convention E_{X} ≡ E/10^{X} in eV will be adopted for photon energies, unless stated otherwise^{8}. In Sect. 2 we have shown that a spontaneously quenched γray spectrum shows a break at the energy . Combining this expression with the fact that the observed break energy usually lies in the GeV energy band, we derive a second relation between the magnetic field and the Doppler factor of the blob (48)Since the break energy is by definition smaller than the maximum energy, the magnetic field B_{br} always satisfies the inequality in Eq. (47). We note that the magnetic field required is, generally, strong. Even for a small value of the Doppler factor, e.g. δ = 10, one needs B ≃ 40 G (see also PM12a for a related discussion). Thus, spontaneously quenched γray spectra cannot operate in the context of onezone leptonic models, such as synchrotron selfCompton (SSC) that usually requires weak magnetic fields (e.g. Böttcher et al. 2009). Finally, the observed maximum energy of the produced soft photons is given by (49)where we have used the magnetic field strength given by Eq. (48).
If a γray emitting blazar happens to be a spontaneously quenched source, then one can make a strong prediction about the flux correlation between soft (usually Xray) photons, the unabsorbed part of the γray spectrum (GeV energy band) and the absorbed one (typically TeV energy band). Increase of the intrinsic γray compactness amplifies, in general, the absorption of VHE γrays, which leads to an increase of the soft photon component. The number of photon targets for the γrays is then increased, which further sustains the nonlinear loop of photonphoton absorption. The part of the γray spectrum that is not affected by automatic quenching follows exactly the variations of the injection compactness, whereas the spontaneously quenched varies in the inverse way. The above are shown in Fig. 7. First, we allowed the system to reach a steady state (for the parameters used see Fig. 7). Then, we imposed a Lorentzian variation on the injected γray compactness: (50)where τ is the comoving time in units of t_{cr} and (51)where τ_{c} and G are free parameters that control the position of the maximum and the width at half maximum, respectively. Then, we calculated the photon compactness of the soft component () of the unabsorbed (ϵ_{min} ≤ ϵ ≤ 10ϵ_{min} ≈ ϵ_{br}) and spontaneously quenched (0.1ϵ_{max} ≤ ϵ ≤ ϵ_{max}) γrays and plotted them as a function of time in Fig. 7. The evolution of the break energy of the spectrum from high to lower energies can be seen in Fig. 8, where snapshots of MW spectra are plotted.
Fig. 7
Plot of photon compactness as a function of time for various components. Absorbed and unabsorbed γrays are plotted with solid and dashed lines, respectively, while soft photons are shown with a dotted line. Parameters used for this plot: ϵ_{max} = 2.3 × 10^{5}, ϵ_{min} = 23, Γ = 1.5, B = 40 G, R = 3 × 10^{16} cm, ℓ_{inj} = 1.7 × 10^{3} ≳ ℓ_{γ,cr}, G = 10, and τ_{c} = 20. 
Fig. 8
Snapshots of MW spectra obtained for a variable γray injection compactness. The type of variation and the parameter values are the same as in Fig. 7. The arrow shows the socalled hardtosoft evolution of the γray spectrum break energy. 
4.2. Application to blazar PKS 0447439
The bright blazar PKS 0447439 has recently been detected at high energy (Abdo et al. 2009a) and VHE (Abramowski et al. 2013) γrays. Although the redshift of the source is still disputable (see e.g. Prandini et al. 2012; Pita et al. 2012, for different estimations and discussions), in the following we will adopt the value z = 0.2. Furthermore, a higher redshift, that is equivalent to a larger optical depth for γγ absorption, would imply a harder intrinsic VHE spectrum and, therefore, a smaller value of ΔΓ ≡ Γ_{TeV} − Γ_{GeV}. However, taking into account the errorbars of the photon indices, a fit can still be achieved for a range of redshift values (0.1 < z < 0.24). We have focused on PKS 0447439 since it satisfies most of the conditions that were presented in the previous section. In particular:

The photon indices in the GeV and TeV energy ranges areΓ_{GeV} = 1.85 ± 0.05 and Γ_{TeV} ≈ 2.5 ± 0.37, respectively^{9}. This is in agreement with what we have shown, i.e. that a γray spectrum with Γ = 1.8 can steepen up to Γ = 2.7 owing to spontaneous quenching.

The Xray luminosity is less than the γray luminosity but of the same order of magnitude.

An anticorrelated variability between VHE γrays and Xrays may be suggested, although the number of data points is small and therefore still inconclusive (see Fig. 5 in Abramowski et al. 2013).
As a first attempt, we did not specify the production mechanism of γrays. Instead, we assumed that they are being injected into the emission region with a rate given by (52)where Q_{0} is a normalization constant that is related to the injection γray compactness as (53)We also included the injection of primary electrons at a rate (54)where Q_{0,e} is related to the electron injection compactness in the same way as in Eq. (53). The parameters of the fit shown in Fig. 9 can be found in Table 2.
Fig. 9
Multiwavelength spectrum of PKS 0447439 during the period November 2009–January 2010. Filled squares represent the Swift/UVOT, Swift/XRT, Fermi, and H.E.S.S. data from low to high energies respectively. Solid and dashed lines show the SEDs with and without the injection of primary electrons. Our model SEDs are corrected for EBL absorption asumming z = 0.2 and using model C of Finke et al. (2010). 
Figure 9 shows that the synchrotron emission of secondary electron/positron pairs explains the Xray emission, while the synchrotron emission from primary leptons is required to fit the optical data. In this example, primary soft photons, i.e. photons that are not produced by automatic quenching of γrays, are not targets for photonphoton absorption because of their low energies, and therefore, linear absorption does not interfere with the nonlinear one. We note that inverse Compton scattering was also taken into account; both primary and secondary soft photon fields were used in the scatterings. As a second step, one could attribute γray emission to a particular production mechanism, e.g. to synchrotron radiation from relativistic protons, since the choice of a relatively large value of the magnetic field makes the application of automatic photon quenching more relevant to hadronic emission models; a detailed hadronic modelling of the source, however, lies outside the scope of the present paper.
5. Discussion
A very interesting, yet largely unnoticed property of γray radiation transfer is the presence of an upper limit at the production rate per volume of γray photons. If the γray compactness at injection exceeds this critical value, then soft photons are produced spontaneously in the source, serve as targets for highenergy photons, and absorb the excessive γray luminosity. Thus, soft photons act as a thermostat and appear irrespective of the γray production mechanism. These ideas were put forward in SK07, who coined the term “automatic photon quenching” to describe this nonlinear mechanism, and were expanded by PM11 and PM12b. The present paper continues the exploration of automatic quenching in the case where γrays are injected with a powerlaw distribution. Therefore, this can be considered a continuation of an earlier work (PM11), in which the quenching in the case of monoenergetic γray injection was studied.
In Sect. 2 we have derived an analytical expression of the critical compactness that is required for an injected powerlaw γray spectrum to be quenched. A series of approximations/assumptions, which were presented in detail in the same section, were necessary for the above derivation. In particular, the catastrophic losses approximation for electrons proved to be crude enough and made our analytical results valid in a particular parameter range; we have commented on that through an indicative example, where the analytically and numerically derived values of the critical compactness were compared. In cases where automatic photon quenching applies, we have calculated the steadystate γray spectra and shown that spontaneous photon absorption produces a break of ΔΓ = Γ/2 between the unabsorbed and the absorbed parts of the injected powerlaw; here Γ is the slope of the injected γ −rays.
In Sect. 3 we have implemented a numerical code that solves the full radiative transfer problem inside a spherical volume, in order to derive the critical γray compactness for a wide range of parameter values, e.g. for various photon indices, as well as for different minimum and maximum energies of the γray spectrum. We have also examined the effects that a primary soft photon component would have on the absorption of γrays. We found that these depend on both the compactness and the spectral shape of the external component.
In Sect. 4 we have examined the implications of quenching on γray emitting blazars. We have also given a set of crireria that can be used in order to deduce, observationally, whether a γray blazar is spontaneously quenched or not. We note that the relevant information is imprinted both on the SED and on the variability patterns of the source. We also point out that blazar PKS 0477439 meets several of the criteria and there is a distinct possibility that its highenergy spectrum is quenched, while the Xrays are produced from the reprocession of γrays.
In the present paper we have intentionally avoided pinpointing a specific mechanism for γray production. However, the range of parameters used and especially the choice of a rather large magnetic field value (of the order of a few Gauss), imply that these ideas are more effectively applied to hadronic models. In this case γrays could be produced either by protonsynchrotron radiation or by pion production (Mastichiadis et al. 2013).
Automatic quenching might have some farreaching implications for γray blazars. If there is evidence that it operates at some level, then part (or all) of the UV/Xray component should be reprocessed γray emission, i.e. there is no need for a primary component to produce all of the observed soft radiation. It also predicts that, for AGN related parameters, breaks in the γray spectra should appear in the high GeV – low TeV regime. Moreover, a hardtosoft evolution of the spectral break is expected whenever the injected γray flux increases. Therefore, future observations, especially with CTA (Sol et al. 2013), could prove decisive in detecting the presence of such spectral breaks. If, on the other hand, the sources do not show signs of spontaneous quenching, then some interesting constraints apply to the source parameters, as this relates the γray luminosity to the size of the source, the magnetic field strength and the Doppler factor (see PM12a for an application of this kind on quasar 3C 279). These constraints can be quite severe, especially if the source undergoes strong flaring episodes and the absence of quenching could only mean either a very large value of the Doppler factor or a low magnetic field in the production region. Both aspects have strong implications for the physical conditions prevailing in the emitting regions of γray blazars.
Similar results are presented in Mastichiadis et al. (2005) but they are caused by a different intrinsic nonlinear process known as the “PPSloop” (Kirk & Mastichiadis 1992).
The VHE photon index is calculated after correcting for EBL absortpion for z = 0.2 and using the model by Finke et al. (2010).
Acknowledgments
We would like to thank the anonymous referee for useful comments/suggestions on the manuscript and Dimitra Lingri for helping us compile Table 1. This research has been cofinanced by the European Union (European Social Fund – ESF) and Greek national funds through the Operational Program “Education and Lifelong Learning” of the National Strategic Reference Framework (NSRF) – Research Funding Program: Heracleitus II. Investing in knowledge society through the European Social Fund.
References
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009a, ApJS, 183, 46 [NASA ADS] [CrossRef] [Google Scholar]
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009b, ApJ, 700, 597 [NASA ADS] [CrossRef] [Google Scholar]
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011a, ApJ, 727, 129 [NASA ADS] [CrossRef] [Google Scholar]
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011b, ApJ, 736, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Abramowski, A., Acero, F., Aharonian, F., et al. (H.E.S.S. Collaboration) 2011, A&A, 533, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Abramowski, A., Acero, F., Aharonian, F., et al. (H.E.S.S. Collaboration) 2012, A&A, 538, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Abramowski, A., Acero, F., Akhperjanian, A. G., et al. (H.E.S.S. Collaboration) 2013, A&A, 552, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Acciari, V. A., Aliu, E., Arlen, T., et al. 2010a, ApJ, 708, L100 [NASA ADS] [CrossRef] [Google Scholar]
 Acciari, V. A., Aliu, E., Arlen, T., et al. 2010b, 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., Barres de Almeida, U., et al. 2007, A&A, 475, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aharonian, F. A., Khangulyan, D., & Costamante, L. 2008, MNRAS, 387, 1206 [NASA ADS] [CrossRef] [Google Scholar]
 Aleksić, J., Antonelli, L. A., Antoranz, P., et al. 2010, ApJ, 723, L207 [NASA ADS] [CrossRef] [Google Scholar]
 Aleksić, J., Antonelli, L. A., Antoranz, P., et al. 2011a, ApJ, 730, L8 [NASA ADS] [CrossRef] [Google Scholar]
 Aleksić, J., Antonelli, L. A., Antoranz, P., et al. 2011b, A&A, 530, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aleksić, J., Alvarez, E. A., Antonelli, L. A., et al. 2012a, A&A, 544, A142 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aleksić, J., Alvarez, E. A., Antonelli, L. A., et al. 2012b, A&A, 539, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aliu, E., Aune, T., Beilicke, M., et al. 2011, ApJ, 742, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Aliu, E., Archambault, S., Arlen, T., et al. 2012, ApJ, 750, 94 [NASA ADS] [CrossRef] [Google Scholar]
 Böttcher, M., & Chiang, J. 2002, ApJ, 581, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Böttcher, M., Reimer, A., & Marscher, A. P. 2009, ApJ, 703, 1168 [NASA ADS] [CrossRef] [Google Scholar]
 Coppi, P. S. 1992, MNRAS, 258, 657 [NASA ADS] [CrossRef] [Google Scholar]
 Coppi, P. S., & Blandford, R. D. 1990, MNRAS, 245, 453 [NASA ADS] [Google Scholar]
 Costamante, L. 2012 [arXiv:1208.0810] [Google Scholar]
 Domínguez, A., Primack, J. R., Rosario, D. J., et al. 2011, MNRAS, 410, 2556 [NASA ADS] [CrossRef] [Google Scholar]
 Finke, J. D., Razzaque, S., & Dermer, C. D. 2010, ApJ, 712, 238 [NASA ADS] [CrossRef] [Google Scholar]
 Franceschini, A., Rodighiero, G., & Vaccari, M. 2008, A&A, 487, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Herterich, K. 1974, Nature, 250, 311 [NASA ADS] [CrossRef] [Google Scholar]
 Jelley, J. V. 1966, Nature, 211, 472 [NASA ADS] [CrossRef] [Google Scholar]
 Kataoka, J., Takahashi, T., Makino, F., et al. 2000, ApJ, 528, 243 [NASA ADS] [CrossRef] [Google Scholar]
 Katarzyński, K., Ghisellini, G., Tavecchio, F., et al. 2005, A&A, 433, 479 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kazanas, D. 1984, ApJ, 287, 112 [NASA ADS] [CrossRef] [Google Scholar]
 Kirk, J. G., & Mastichiadis, A. 1992, Nature, 360, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Konopelko, A., Mastichiadis, A., Kirk, J., de Jager, O. C., & Stecker, F. W. 2003, ApJ, 597, 851 [NASA ADS] [CrossRef] [Google Scholar]
 Mastichiadis, A., & Kirk, J. G. 1995, A&A, 295, 613 [NASA ADS] [Google Scholar]
 Mastichiadis, A., & Kirk, J. G. 1997, A&A, 320, 19 [NASA ADS] [Google Scholar]
 Mastichiadis, A., Protheroe, R. J., & Kirk, J. G. 2005, A&A, 433, 765 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mastichiadis, A., Petropoulou, M., & Dimitrakoudis, S. 2013, MNRAS, submitted [arXiv:1304.2957] [Google Scholar]
 Petropoulou, M., & Mastichiadis, A. 2009, A&A, 507, 599 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Petropoulou, M., & Mastichiadis, A. 2011, A&A, 532, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Petropoulou, M., & Mastichiadis, A. 2012a, MNRAS, 426, 462 [NASA ADS] [CrossRef] [Google Scholar]
 Petropoulou, M., & Mastichiadis, A. 2012b, MNRAS, 421, 2325 [Google Scholar]
 Pita, S., Goldoni, P., Boisson, C., et al. 2012, in AIP Conf. Ser., 1505, eds. F. A. Aharonian, W. Hofmann, & F. M. Rieger, 566 [Google Scholar]
 Prandini, E., Bonnoli, G., & Tavecchio, F. 2012, A&A, 543, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Protheroe, R. J. 1986, MNRAS, 221, 769 [NASA ADS] [Google Scholar]
 Sol, H., Zech, A., Boisson, C., et al. 2013, Astropart. Phys., 43, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Stawarz, Ł., & Kirk, J. G. 2007, ApJ, 661, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Stern, B. E., Begelman, M. C., Sikora, M., & Svensson, R. 1995, MNRAS, 272, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Svensson, R. 1987, MNRAS, 227, 403 [NASA ADS] [Google Scholar]
 Zdziarski, A. A., & Lightman, A. P. 1985, ApJ, 294, L79 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: γray spectral break due to absorption on primary soft photons
We derive the spectral break of a γray spectrum due to the absorption of a primary soft photon component (n_{ex}) that is present in the emission region. In this case, the absorption of γrays is a linear process in contrast to automatic quenching. We treat the target photon population as a photon tank, in the sense that n_{ex} does not evolve with time. Thus, the steadystate γray photon spectrum is derived by solving the same set of equations as those presented in Sect. 2.3, with a slightly different operator (and ): (A.1)where X_{M} = max [ x_{max},x_{0} ] and x_{max},x_{0} are the maximum energies of the secondary and primary soft photons, respectively. The above expression represents the most general physical case, where γrays are being absorbed by both primary and secondary soft photons. One can distinguish between two regimes:

1.
spontaneous or nonlinear γray absorption, ifn_{ex}(x) ≪ n_{s}(x), and

2.
linear γray absorption, if n_{ex}(x) ≫ n_{s}(x).
Here we focus on the second regime, where after following the same steps as those described in Sect. 2.3, we find the steadystate γray distribution (A.2)where n_{0} is the normalization of the primary soft photon distribution. Since we are interested in calculating the photon index of the absorbed part of the spectrum, it is sufficient to look at the asymptotic expression of n_{γ}(ϵ), (A.3)that is obtained for ϵ^{s−1} ≫ 2^{s}s/σ_{0}n_{0}. Thus, in this regime the spectral break is given by (A.4)
All Tables
All Figures
Fig. 1
Top panel: critical compactness ℓ_{γ,cr} as a function of the minimum energy of the γray spectrum ϵ_{min} for Γ = 1.6 (solid line) and Γ = 2 (dashed line). The numerically derived values for the two cases are shown with circles and triangles, respectively. Bottom panel: same as in top panel except for different photon indices, which are marked on the plot. Other parameters used are: ϵ_{max} = 2.3 × 10^{5} (in m_{e}c^{2} units), B = 40 G, and R = 3 × 10^{16} cm. In both panels the grey area denotes the region where ϵ_{min} > ϵ_{q}; see text for the definition of ϵ_{q}. 

In the text 
Fig. 2
Analytical solution of steadystate spontaneously quenched γray spectra for different values of the injection parameter Q_{0} starting from 1.25 × 10^{3}. Each value is increased by a factor of 5 over its previous value. Crosses denote in each case the value of ϵ_{⋆}. Other parameters used are: Γ = 2, ϵ_{min} = 33, ϵ_{max} = 2.3 × 10^{5}, B = 40 G, and R = 3 × 10^{16} cm. The break energy is then ϵ_{br} = 160. 

In the text 
Fig. 3
Numerical (dashed lines) and analytical (solid lines) solutions for the steadystate γray spectra. The injection rate increases from bottom to top starting from Q_{0} = 1.25 × 10^{3}. The parameters used are the same as in Fig. 2. 

In the text 
Fig. 4
Log–log diagram of the critical γray compactness as a function of ϵ_{min} for constant ϵ_{max} = 2.3 × 10^{5} (top panel) and as a function of ϵ_{max} for constant ϵ_{min} = 1.4 × 10^{5} (bottom panel) for three photon indices marked on the plot. Other parameters used are B = 40 G and R = 3 × 10^{16} cm. 

In the text 
Fig. 5
Top panel: comparison between two limiting cases where γray photons are absorbed only on synchrotron photons emitted by secondary pairs (solid line) or by primary electrons (dashed line). The injection compactnesses of γrays and primary electrons are (solid line) and (dashed line), respectively. Different parts of the spectra have different powerlaw dependencies that are marked on the plot. Bottom panel: example of an indermediate case where γrays are being partially absorbed on synchrotron photons from secondaries. The injection compactness of primary electrons is , while ℓ_{inj} takes the following values: 8 × 10^{3} (solid line), 8 × 10^{2} (dashed line), and 8 × 10^{1} (dotted line). The dashdotted line that is obtained for (same as solid line in top panel) is plotted for comparison reasons. Other parameters used for the plot are: B = 40 G and R = 3 × 10^{16} cm. 

In the text 
Fig. 6
Photon index of the TeV spectrum versus the one of the GeV spectrum for several γray emitting blazars. 

In the text 
Fig. 7
Plot of photon compactness as a function of time for various components. Absorbed and unabsorbed γrays are plotted with solid and dashed lines, respectively, while soft photons are shown with a dotted line. Parameters used for this plot: ϵ_{max} = 2.3 × 10^{5}, ϵ_{min} = 23, Γ = 1.5, B = 40 G, R = 3 × 10^{16} cm, ℓ_{inj} = 1.7 × 10^{3} ≳ ℓ_{γ,cr}, G = 10, and τ_{c} = 20. 

In the text 
Fig. 8
Snapshots of MW spectra obtained for a variable γray injection compactness. The type of variation and the parameter values are the same as in Fig. 7. The arrow shows the socalled hardtosoft evolution of the γray spectrum break energy. 

In the text 
Fig. 9
Multiwavelength spectrum of PKS 0447439 during the period November 2009–January 2010. Filled squares represent the Swift/UVOT, Swift/XRT, Fermi, and H.E.S.S. data from low to high energies respectively. Solid and dashed lines show the SEDs with and without the injection of primary electrons. Our model SEDs are corrected for EBL absorption asumming z = 0.2 and using model C of Finke et al. (2010). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.