Implications of automatic photon quenching on compact gammaray sources
Department of PhysicsUniversity of Athens, Panepistimiopolis, 15783 Zografos, Greece
Received: 22 February 2011
Accepted: 26 April 2011
Aims. We investigate photon quenching in compact nonthermal sources. This involves photonphoton annihilation and lepton synchrotron radiation in a network that can become nonlinear. As a result the γray luminosity of a source cannot exceed a critical limit that depends only on the radius of the source and on the magnetic field.
Methods. We perform analytic and numerical calculations that verify previous results and extend them so that the basic properties of photon quenching are investigated.
Results. We apply the above to the 2006 TeV observations of quasar 3C 279 and obtain the parameter space of allowed values for the radius of the emitting source, its magnetic field strength and the Doppler factor of the flow. We argue that the TeV observations favour either a modest Doppler factor and a low magnetic field or a high Doppler factor and a high magnetic field.
Key words: infrared: general / acceleration of particles / radiation mechnisms: nonthermal
© ESO, 2011
1. Introduction
Astrophysical γray sources are often subject to the socalled “compactness” problem: when the flux at longer wavelengths is combined with the dimensions of the source, which are inferred from lightcrossing time arguments, a high optical thickness to photonphoton annihilation and effective absorption of γrays are implied.
The possibility that highenergy photons could pairproduce on soft target photons instead of escape in compact sources was first discussed by Jelley (1966), while the ratio of the luminosity to the size of the source, L/R, emerges as a determining factor of whether or not a highenergy photon will actually be absorbed (Herterich 1974). This was followed by work that took into account photonphoton annihilation not only as a sink of γrays, but also as a source of electronpositron pairs inside nonthermal compact sources (Bonometto & Rees 1971; Guilbert et al. 1983; Kazanas 1984; Zdziarski & Lightman 1985; Svensson 1987). The aim of these models was to calculate selfconsistently the photon flux escaping from the sources by taking into account the energy redistribution caused by the γannihilationinduced pair cascades. The assumption was that highenergy particles or γrays were injected uniformly in a source that also contains soft photons, and the system was followed to its final steadystate through the solution of a set of kinetic equations describing the physical processes at work. Several groups have also developed numerical codes for the computation of timedependent solutions to the kinetic equations taking photonphoton annihilation into account (Coppi 1992; Mastichiadis & Kirk 1995; Stern et al. 1995; Böttcher & Chiang 2002). These algorithms are commonly used in source modelling (Mastichiadis & Kirk 1997; Kataoka et al. 2000; Konopelko et al. 2003; Katarzyński et al. 2005).
One of the assumptions of these models, as stated above, is the presence of soft photons in the source, which serve as targets for the γray annihilation. A different approach has been recently presented by Stawarz & Kirk (2007, henceforth SK). These authors focused on the nonlinear effects induced by photonphoton annihilation and investigated the necessary conditions under which γray photons can cause runaway pair production^{1}. Stawarz & Kirk (2007) showed that there is a limit to the γray luminosity escaping from a source, which does not rely on the existing soft photon population but is instead a theoretical limit depending only on parameters such as the source size and its magnetic field strength. Violation of this limit leads to automatic quenching of the γrays. This involves a network of processes, namely photonphoton annihilation and lepton synchrotron radiation, which can become nonlinear once certain criteria are satisfied. In this case electronpositron pairs grow spontaneously in the system and the “excessive” γrays are absorbed on the synchrotron photons emitted by the pairs. As a result the system reaches a final steady state where the γrays, soft photons, and electronpositron pairs have all reached equilibrium. Therefore this can occur even in the hypothetical case when there are no soft photons in the source, at least initially.
Indeed, photon quenching has some interesting implications for emission models of γrays because it gives a robust upper limit on the γray luminosity. The aim of the present paper is to investigate the implications of this network for the highenergy compact astrophysical sources. In Sect. 2 we will expand the analytical approach of SK, who studied the dynamical system of photons and relativistic pairs using δfunctions for the photonphoton annihilation cross section and for the synchrotron emissivity, while they treated synchrotron losses as catastrophic. In Sect. 3 we will use a numerical approach that will allow us to study the properties of quenching using the full cross section for γγ absorption and the full synchrotron emissivity. As an example, in Sect. 4 the above will be applied to the 2006 MAGIC TeV observations of quasar 3C 279 to extract a parameter space of allowed values for the source parameters, i.e., the radius R and the magnetic field strength B, as well as for the Doppler factor δ of the flow. Finally in Sect. 5 we conclude and give a brief discussion of the main points of the present work.
2. Analytical approach
2.1. First principles
We begin by assuming a spherical source of radius R containing a tangled magnetic field of strength B. We also assume that monoenergetic γrays of energy ϵ_{γ} (in units of m_{e}c^{2}) are uniformly produced by some unspecified mechanism throughout the volume of the source. If these are injected with a luminosity , one can define 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 SK showed, the injected γray compactness cannot become arbitrarily high because if a critical value is reached, the following loop starts operating:

1.
gammarays pairproduce on soft photons, which can bearbitrarily low inside the source;

2.
the produced electronpositron pairs cool by emitting synchrotron photons, thus acting as a source of soft photons;

3.
the soft photons serve as targets for more γγ interactions.
There are two conditions that should be satisfied simultaneously for this network to occur: the first, which is a feedback condition, requires that the synchrotron photons emitted from the pairs have sufficient energy to pairproduce on the γrays. By making suitable simplifying assumptions, one can derive an analytic relation for it – see also SK. Thus, combining (1) the threshold condition for γγabsorption ϵ_{γ}ϵ_{0} = 2; (2) the fact that there is equipartition of energy among the created electronpositron pairs γ_{p} = γ_{e} = γ = ϵ_{γ}/2; and (3) the assumption that the required soft photons of energy ϵ_{0} are the synchrotron photons that the electrons/positrons radiate, i.e., ϵ_{0} = bγ^{2} where b = B/B_{crit} and G is the critical value of the magnetic field, one derives the minimum value of the magnetic field required for quenching to become relevant (2)Thus for B ≥ B_{q} the feedback criterion is satisfied. Note that this is a condition that only contains the emitted γray energy and the magnetic field strength; moreover, it can easily be satisfied, at least if one is to use the values inferred from typical modelling of the sources (Böttcher 2007; Böttcher et al. 2009).
Along with the feedback criterion, a second criterion must be imposed for the quenching to be fully operative. A simple way to see this is with the following consideration. Assume that the γrays pairproduce on some soft photon and that the created electronpositron pairs cool by emitting synchrotron photons. Because an electron emits several such photons before cooling the critical condition occurs if the number density of the γrays is such that at least one of the synchrotron photons pairproduces on a γray instead of escaping from the source. Thus the condition for criticality can be written as (3)where n(ϵ_{γ}) is the number density of γrays, σ_{γγ} is the cross section for γγ interactions and is the number of synchrotron photons emitted by an electron with Lorentz factor γ before it cools. Assuming that the pairproducing collisions occur close to threshold, approximating the cross section there by σ_{γγ} ≃ σ_{T}/3 and using and , the critical condition can be written as (4)where the feedback condition (2) and Eq. (1) were also used. Taken as equalities, relations (2) and (4) define the marginal stability criterion, which essentially is a condition for the γray luminosity.
2.2. Kinetic equations
We begin by writing the kinetic equations that describe the distributions of γray photons, soft photons, and electrons in the source. These are respectively and (7)where n, n_{0} and, n_{e} are the differential γray, soft photon, and electron number densities, respectively, and ϵ_{γ}, x, and γ are the corresponding energies normalized in m_{e}c^{2} units. The densities refer to the number of particles contained in a volume element σ_{T}R. In other words, if expresses the number of particles of species i per ergs per cm^{3}, then . Time is normalized with respect to the photon crossing/escape time from the source, t_{cr} = R/c. Thus, τ, which appears in the kinetic equations, is dimensionless and equals . The operators and ℒ denote injection and losses, respectively. The only processes that we take into account are γγ annihilation and synchrotron cooling of the produced pairs. The synchrotron emissivity is approximated by a δfunction, i.e., j_{s}(x) = j_{0}δ(x − ϵ_{0}), where ϵ_{0} = bγ^{2} is the synchrotron critical energy. In all cases we will assume a monoenergetic γray injection at energy ϵ_{γ}.
2.2.1. δfunction approximation for the σ_{γγ} cross section
First we will examine the stability of the system using the simplest form for the photonphoton annihilation cross section (Zdziarski & Lightman 1985): (8)where x,ϵ_{γ} are the energies of the soft and γray photons, respectively. The particle injection and loss operators take the following forms: and (13)because each photonphoton annihilation results in a pair of leptons, each one with approximately half of the initial γray energy. In the above equations the “magnetic compactness” was introduced (14)The system of Eqs. (5)–(7) can now be written as (15)In this section we will examine the stability of the trivial stationary solution of the system (15): (16)which corresponds to the free propagation of γray photons through the source. To investigate the stability of the system, we assume that initially arbitrarily small perturbations of the soft photon and electron densities are present, which lead to the perturbation of the γray photon density n′. After linearization, the system (15) becomes (17)The stability analysis of the system above can be simplified when one works with the Laplace transformed number densities: (18)with s being the solution to the eigenvalue problem: For large τ both the soft photon and electron distributions behave as e^{sτ}. If s = −1, then and the γray photons freely escape from the source. If s > 0, one finds solutions that grow with elapsing time. The marginally stable solution, which is obtained when one sets s = 0, will be examined below.
Combining Eqs. (19)–(21) leads to an ordinary differential equation for the perturbed electron distribution with the following solution: (22)where . Setting γ = γ_{cr} in the above equation one finds (23)Then the perturbed photon distribution is found using Eq. (20).
As already mentioned in Sect. 2.1, the γray compactness cannot become arbitrarily high but instead saturates at a critical value . When only a negligible fraction of the hard photon luminosity will be absorbed, i.e., , or equivalently . We call critical the injected luminosity required for to occur. After this point any further increase of the injected luminosity will be seen as an increase of the luminosity in soft photons. Using Eqs. (19), (20), and (23) one derives (24)where . Solving the above equation with respect to , we find (25)To compare our results with the value found by SK, who have used catastrophic synchrotron losses, we have to take the following steps:

1.
the case of synchrotron cooling degenerates to the“catastrophic losses case” when the threshold energy forabsorption 2/ϵ_{γ} equals the maximum soft photon energy x_{max};

2.
in SK the volume element is defined as V = πR^{3}, whereas we used V = 4/3πR^{3}. Thus, the replacement should be made;

3.
finally, in SK the luminosity compactness is not given by the conventional definition, which we adopted, but is multiplied by the energy ϵ_{γ}.
Thus, our result should be compared to . The critical compactness as a function of ϵ_{γ} for two values of the magnetic field is shown in Fig. 1. The two crosses mark the values of SK, which completely agree with ours in the limit of catastrophic losses.
Fig. 1 Loglog plot of the critical luminosity compactness as a function of the injected γray photon energy ϵ_{γ}, when the δfunction approximation for the cross section σ_{γγ} was used for two values of the magnetic field strength: B = 40 G (solid line) and B = 4 G (dashed line). The values calculated by SK using only catastrophic losses are marked with crosses. 

Open with DEXTER 
We note also that Eq. (4), which was obtained by making robust calculations, also agrees with Eq. (25) in the limit of catastrophic losses.
2.2.2. Step function approximation for the σ_{γγ} cross section
In this paragraph we analytically find the critical luminosity compactness after adopting a more realistic expression for the annihilation cross section: (26)where the normalization used is σ_{0} = 4/3. The initial Eqs. (5)–(7) keep the same form. Only the operator of photon losses becomes (27)Following the same steps as those described in Sect. 2.2.1, we find for s = 0 (28)where and γ_{max} = ϵ_{γ}/2. Equation (28) can be solved iteratively when writing as a sum of approximations, i.e., . Assuming for the first approximation the form n^{(0)} = Cγ^{2}, we find the other terms of the series expansion: where and λ < 1 for typical values of the parameters. Thus, the series of approximations converges and can be found in closed form: (31)The luminosity in soft and hard photons can then be found (32)(33)The critical γray compactness is given by (34)Figure 2 shows as a function of ϵ_{γ} for two values of the magnetic field. When the threshold energy for absorption 2/ϵ_{γ} is much lower than the maximum energy of the soft synchrotron photons , the critical compactness as a function of ϵ_{γ} behaves as for the δfunction cross section. However, when the two energies become comparable, abruptly increases, because the number of soft photons serving as targets for the absorption significantly decreases. This effect can be seen in Fig. 2 as a change in the curvature of the curves.
Fig. 2 Loglog plot of the critical luminosity compactness as a function of the injected γray photon energy ϵ_{γ}, when the Θfunction approximation for the cross section σ_{γγ} was used for two values of the magnetic field strength: B = 40 G (solid line) and B = 4 G (dashed line). 

Open with DEXTER 
3. Numerical approach
To numerically investigate the properties of quenching, one needs to solve again the system of Eqs. (5)–(7) augmented to include more physical processes. 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: (35)and (36)where n_{e} and n_{γ} are the differential electron and photon number density, respectively, normalized as in Sect. 2. Here we considered the following processes: (i) photonphoton pair production, which acts as a source term for electrons () and a sink term for photons (); (ii) synchrotron radiation, which acts as a loss term for electrons () and a source term for photons (); (iii) synchrotron selfabsorption, which acts as a loss term for photons (); and (iv) inverse Compton scattering, which 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 (1997) and Petropoulou & Mastichiadis (2009). The photons are assumed to escape the source in one crossing time, therefore t_{γ,esc} = R/c.
The parameters of the problem are completely specified once the values of radius R and of the magnetic field B are set and the injection rate is specified. Because is readily related to the injected luminosity through the relation (37)where V is the volume of the source, the injection rate is specified once the injection compactness and the functional dependence of on ϵ_{γ} is set. The electron physical escape timescale from the source t_{e,esc} is another free parameter which, 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 synchrotron photons, we assume that at t = 0 there are no electrons or photons in the source, so we set n_{e}(γ,0) = n_{γ}(x,0) = 0.
3.1. Monoenergetic γray injection
As a first example we study the case where the γrays injected are monoenergetic at an energy ϵ_{γ}. The expected trivial solution of the kinetic equations is that because there are no soft photons in the source, the injected γrays will freely escape and their spectrum upon escape will simply be their production one. Electrons and photons at energies different from ϵ_{γ} will remain at their initial values, i.e., n_{e}(γ,t) = n_{γ}(x,t) = 0.
Fig. 3 Critical γray compactness as a function of ϵ_{γ} when (i) σ_{γγ} and synchrotron emissivity are approximated by a Θ and a δfunction respectively; and (ii) when for both quantities the full expressions are used. For the first case the numerical result (points) and the analytical one given by Eq. (34) (solid line) are shown. The full problem can only be treated numerically and the corresponding result is shown with a dashed line. The size of the source is R = 3 × 10^{16} cm and the magnetic field strength B = 3.57 G. 

Open with DEXTER 
Before numerically examining the stability of the above system using the full expressions for the photonphoton annihilation cross section and the emissivities, we perform a test by comparing the analytical result given by Eq. (34) with the respective numerical. The two solutions completely agree as shown in Fig. 3. In the same figure we also plot the critical γray compactness found numerically when all quantities used were given by their full expressions. The qualitative behaviour of the full problem solution is quite well reproduced by the analytical one obtained while using the Θ and δapproximations for the cross section and synchrotron emissivity respectively. This allows us to use the analytical expression given by Eq. (34) for the full problem as well, when investigating the properties of quenching at least qualitatively. From this figure it becomes apparent that γray quenching can affect all photons with energies that satisfy the feedback condition. Starting around this energy and moving progressively to higher ones, we find that the critical compactness first sharply decreases until it reaches a minimum and then starts increasing again, as asymptotically. The exact value of the critical compactness also depends on the strength of the magnetic field. For a higher Bfield the required minimum energy ϵ_{γ} decreases, which can be easily deduced by Eq. (2). Note also that because of the marginal relation (34), the required critical injected compactness becomes higher. Thus, by increasing B, the corresponding curve is shifted towards lower γray energies and higher critical compactnesses, keeping its shape (see e.g. Fig. 2).
Figure 4 depicts the results of a series of runs where we kept ϵ_{γ} constant and varied . For sufficiently low values of the latter parameter, all γrays indeed escape from the system without any impedance. If we call the compactness of the outgoing radiation at energy ϵ_{γ}, then obviously, in this case, , where is the compactness of the outgoing luminosity. However, above a critical value of , which we will call , the loop described above starts operating, transferring the γray luminosity to softer radiation. Thus, for , we get because saturates – by we denote the compactness of the reprocessed luminosity, which appears at lower energies. This can be seen in Fig. 4, where the dashed line curve represents , the dotted line and the full line , which is always equal to because there are no sinks of photons present. For the particular example shown in Fig. 4, the values of the other parameters used are, R = 3 × 10^{16} cm, B = 40 G and ϵ_{γ} = 2.3 × 10^{4}.
Fig. 4 Plot of the compactness of emerging hard (dashed), soft (dotted), and total radiated luminosity (full line) as a function of the injected compactness . The γrays are assumed to be monoenergetic, injected at energy ϵ_{γ} = 2.3 × 10^{4}. The size of the source is R = 3 × 10^{16} cm, while the magnetic field strength is B = 40 G. The emerging soft radiation is a result of the selfquenching of the system. 

Open with DEXTER 
Fig. 5 Multiwavelength spectrum for a choice of such that . The other parameters are the same as in Fig. 4. Full line corresponds to the case where all relevant radiative processes were taken into account, while the dashed line corresponds to the case where only the two “core” processes were considered, i.e., synchrotron and γγ interactions. 

Open with DEXTER 
Figure 5 shows the multiwavelength (MW) spectrum for the particular value of , which results in equipartition between soft and hard luminosities, i.e., . While the hard luminosity is simply the escaping monoenergetic radiation at ϵ_{γ}, the soft spectrum mainly consists of the synchrotron component produced from the secondary pairs, which were injected at energies γ_{p} = γ_{e} = ϵ_{γ}/2. The peak of this component is at , while it has a spectral index of − 1/2 (or + 1/2 in νF_{ν} units) as a result of the electron cooling to a γ^{2} distribution below the energy of injection. There is also a synchrotron selfCompton (ssc) component, which is greatly suppressed by KleinNishina effects. As can be seen from the figure, the corrections introduced from the inclusion of inverse Compton scattering and synchrotron selfabsorption are negligible, at least for our set of parameters.
Fig. 6 Time evolution of the soft (dotted line), hard (dashed line) and total l_{γ} (solid line) compactness for the same set of parameters as in Fig. 5. 

Open with DEXTER 
Figure 6 shows the time evolution of the different luminosity compactnesses for the example case presented in Fig. 5. Although the dynamical, nonlinear system of equations describing the system (see for example system of Eq. (15)) is of third order, which in principle can lead to unstable solutions or/and “chaotic” behaviour (SK), the system quickly reaches a stable state after a few crossing times. Figure 6 shows that although initially absent, a soft photon population builds up, which leads to a decrease of , because γrays are being absorbed by the soft photons.
3.2. Powerlaw injection
As a next step we investigate the effects of quenching for a powerlaw injection of γrays, i.e., when the injection rate can be written as (38)where Q_{0} is a normalization constant depending on the value of α, the low and high energy cutoffs of the distribution.
Before proceeding to more quantitative calculations, we remark that obviously, for a given B the feedback criterion for quenching to start will be that of Eq. (2) where ϵ_{γ} should be replaced by ϵ_{γmax}. Assuming that this condition holds and that the compactness of the γrays is above the critical one, the production of synchrotron photons with energy is guaranteed. This means that γrays with energy ϵ_{γ} ≳ ϵ_{γmin,eff} should interact in γγ interactions, where (39)Thus, depending on the relation between ϵ_{γmin,eff} and ϵ_{γmin} we can distinguish two cases:

(i)
ϵ_{γmin} < ϵ_{γmin,eff}: in this case quenching will affect γrays from ϵ_{γmax} down to ϵ_{γmin,eff}. For ϵ_{γ} ≤ ϵ_{γmin,eff} the spectrum will be unaffected, i.e., it will be the same as the one at injection.

(ii)
ϵ_{γmin,eff} < ϵ_{γmin}: in this case quenching will affect the whole γray spectrum.
The above two cases are exemplified in Figs. 7 and 8. In both figures we assumed that γrays with a powerlaw form of slope α = 2 are injected in a source of radius R = 3 × 10^{16} cm and magnetic field strength B = 40 G. The only difference between the two cases is in the lower and upper limits, which are ϵ_{γmin} = 23 and ϵ_{γmax} = 2.3 × 10^{4} for the case shown in Fig. 7 and ϵ_{γmin} = 230 and ϵ_{γmax} = 2.3 × 10^{5} for the case shown in Fig. 8.
Fig. 7 MW spectra for a powerlaw injection of γrays of the form for ϵ_{γmin} < ϵ_{γ} < ϵ_{γmax} with ϵ_{γmin} = 23 and ϵ_{γmax} = 2.3 × 10^{4} for various values of the normalization Q_{0}. This starts from Q_{0} = 4 × 10^{3} and then exceeds its previous value by a factor of 2. The dashed line curves represent the runs without quenching. The source parameters are B = 40 G and R = 3 × 10^{16} cm. 

Open with DEXTER 
In the case of Fig. 7 we have ϵ_{γmin} < ϵ_{γmin,eff} = 1.7 × 10^{4} and, as can be seen from the figure, the powerlaw breaks above ϵ_{γmin,eff} because of quenching – indeed, we find that the value of ϵ_{γmin,eff} is in reality lower by a factor of 2 than the one given in Eq. (39), because in the numerical calculations shown here we use the full expressions for the synchrotron emissivity and the crosssection for γγ and not the δfunction approximations implied in the derivation of the above relation. It is interesting to note that as the γray luminosity increases and thus the quenching becomes more effective, the break between the unabsorbed and absorbed part of the powerlaw becomes more abrupt.
In the case of Fig. 8, on the other hand, the choice of the parameters were such as to make ϵ_{γmin,eff} comparable to ϵ_{γmin}. Thus, according to the qualitative analysis presented above, quenching should affect the whole γray distribution. Indeed, as can be deduced from the figure, as the luminosity increases, absorption is moving progressively to lower γray energies, until it affects the whole distribution.
Fig. 8 Same as in Fig. 7 with ϵ_{γmin} = 2.3 × 10^{2} and ϵ_{γmax} = 2.3 × 10^{5} with Q_{0} starting from Q_{0} = 2.5 × 10^{4}. 

Open with DEXTER 
We also investigated the relation between the slope of the injected (α_{inj}) and the “quenched” (α_{que}) γray spectrum. For that reason, we ran the code for different α_{inj}, choosing each time suitable values of the source parameters to enable quenching. Figure 9 shows the result of the runs (solid line) and the slope α_{inj} + 1 (dashed line) for comparison reasons. We find that . Thus, quenching becomes more pronounced as the injected spectrum gets steeper.
Fig. 9 Slope of the quenched γray spectrum as a function of the injected spectrum’s slope (solid line). To facilitate a comparison the line α_{inj} + 1 is also shown (dashed line). 

Open with DEXTER 
4. Application to the onezone homogeneous models: the 3C 279 case
The above has a straightforward application to the socalled onezone homogeneous models, which are customarily used to fit the MW spectra of Active Galactic Nuclei (AGN) (Böttcher et al. 2009; Aleksić et al. 2011). According to these models, the “source” contains electrons and/or protons, which radiate through various processes and form a spectrum that can then be fitted to the observed data by varying certain key parameters – for a detailed description see Mastichiadis & Kirk (1997) and Konopelko et al. (2003).
Alternatively, one can assume that the γrays are produced by some mechanism and investigate the parameter space that allows the escape of these γrays. To perform this analysis, one does not require the existence of other softer photons that act as targets, but one simply has to look for the threshold conditions for photon quenching. In other words, the relevant question one has to ask is: assume a spherical source of radius R containing a magnetic field B and producing a γray spectrum that is measured at Earth with a flux F_{γ}(ϵ_{γ}). If the source is moving with a Doppler factor δ with respect to us, what is the parameter space of R, B, and δ that allows the γrays to escape? As we will show, this can give very robust lower limits for the Doppler factor δ.
The luminous quasar 3C 279 is a good example. A recent comprehensive review of observations can be found in Böttcher et al. (2009). Here we will mainly focus on the 2006 campaign, which discovered the source at VHE γrays, showing a high TeV flux (MAGIC Collaboration et al. 2008), while the Xrays were at a much lower level.
To apply the ideas discussed in the previous section, we follow the procedure below.
We assume a spherical source of radius R moving with Doppler factor δ with respect to us and containing a magnetic field of strength B. We assume that γrays are produced by some unspecified mechanism inside the source and, upon escape, they form the observed γray spectrum. If is the compactness of the γrays (cf. Eq. (1)), we can relate it to the observed integrated flux by the relation (40)where D_{L} = 3.08 Gpc is the luminosity distance to the source. Here we used a cosmology with Ω_{m} = 0.3, Ω_{Λ} = 0.7 and H_{0} = 70 km s^{1} Mpc^{1}.
To model the TeV emission, we performed a χ^{2} fit to the MAGIC data, corrected for intergalactic γγ absorption (Böttcher et al. 2009), of the form νF_{ν} = Aν^{ − α} for where A = 1.3 × 10^{31} Jy Hz, α = 0.7, Hz and Hz. Then we injected γrays in the code of the form for ϵ_{γmin} < ϵ_{γ} < ϵ_{γmax} where the normalization is set by condition (37) with the integration taking place between and .
Fig. 10 MW spectra of 3C 279 for R = 10^{16} cm, B = 40 G and δ = 10^{0.8} (dotted line), 10^{1.0} (dashed line), 10^{1.2} (dotdashed line) and 10^{1.4} (solid line). In all cases the γray injection would have produced the solid line in the absence of quenching. Black squares and bowtie represent the observational data from February 2006. 

Open with DEXTER 
Fig. 11 MW spectra of 3C 279 for R = 10^{16} cm, δ = 10^{1.2} and B = .16 G (solid line), 1.6 G (dashed line) and 16 G (dotted line). In all cases the γray injection would have produced the solid line in the absence of quenching. Black squares and bowtie represent the observational data from February 2006. 

Open with DEXTER 
Figure 10 depicts the resulting MW spectra for R = 10^{16} cm, B = 40 G and δ varying between 10^{0.8} to 10^{1.4} with increments of 0.2 in the exponent. Only the run with the highest value of δ can fit the TeV data. The other three, even though they have parameters that could also fit the data, they do not because of quenching. It is important to repeat at this point that the γray luminosity produced in these cases is the maximum possible that the system can radiate for the particular choice of R, B, and δ. Any effort to increase the γray flux by e.g. increasing would only result in an increase of the reprocessed radiation.
Figure 11 depicts a series of similar runs with R = 10^{16} cm, δ = 10^{1.2} and B = .16 G (solid line), B = 1.6 G (dashed line) and B = 16 G (dotted line). While for low B values there is no quenching, this begins to appear at higher values of the magnetic field. Note, however, that the value of B = 1.6 G produces an acceptable fit, because the soft reprocessed radiation does not violate any observable limit (most notably the Xrays) and it can produce a marginally acceptable fit to TeV γrays. Based on this, one can define a limiting value of the magnetic field B_{q,mx}(R,δ) with the property that if B > B_{q,mx}, the quenching does not allow the internally produced γray luminosity to escape. Therefore B_{q,mx} is the highest allowed value that the Bfield can take for the observed flux of TeV γrays.
Figure 12 depicts the locus B_{q,mx} as a function of the size of the source R for various Doppler factors δ. Each line separates the parameter space into two regions: because of the quenching, a choice of R and B above the locus produces a far too low γray flux that cannot fit the TeV observations. For example, if one assumes that δ = 10 and R = 2.5 × 10^{16} cm – so that it corresponds to a variability timescale of 1 day, the magnetic field cannot exceed 0.3 G. However, owing to the steep dependence of the flux on δ (cf. Eq. (40)), the above condition is greatly relaxed for higher values of δ.
We should also mention that B_{q,mx} not only depends on δ and R, but also on the upper cutoff ϵ_{γmax} of the injected powerlaw in γrays. For the calculations performed here, we chose ϵ_{γmax} such as to match the highest observed frequency in the data set, i.e., 10^{26} Hz. If one assumes that the internally produced γrays can extend to even higher energies, then the value of B_{q,mx} is affected as well. Accordingly we found that if we were to repeat the same calculations with ϵ_{γmax} ten times higher than the value we used, then quenching would affect the flux more and the lines of Fig. 12 would have to shift upwards by a factor of 2. Therefore, Fig. 12 can be taken as a conservative lower limit on B_{q,mx}.
Fig. 12 Parameter space of the allowed values of the magnetic field B as a function of the size of the source R for Doppler factors of δ = 10 (dotted line), 15.8 (dashed line), and 25 (solid line). The space above a line for a certain δ corresponds to combinations of R and B that cannot successfully fit the TeV observations of 3C 279. 

Open with DEXTER 
5. Summary and discussion
We have examined the effects of photon quenching on compact γray sources. Automatic photon quenching, first proposed by SK, is a nonlinear network of two processes that are very common in highenergy astrophysics, namely photonphoton annihilation and lepton synchrotron radiation. We verified and expanded both analytically and numerically the results of SK and we showed that this network, when applied to compact γray sources, can give limits on the Doppler factor, which depend on the source radius and on the magnetic field strength only. This is a theoretically derived limit and as such does not depend on the soft photon populations that might be present inside the source.
It is important to note that photon quenching can be investigated only when the radiation transfer problem is solved selfconsistently because it is a nonlinear process that involves redistribution of energy from the γrays to the lower energy parts of the photon spectrum. For this we used a system of kinetic equations that treats in the simplest case γray annihilation and electron synchrotron radiation. Using suitable δfunction approximations for the photonphoton annihilation cross section and for the photon synchrotron emissivity, we were able to verify the results of SK in the limit where synchrotron losses can be considered as catastrophic (see Fig. 1). On the other hand, using a more realistic step function for the cross section, we were able to find a very good agreement with numerical results, which utilized the full photonphoton annihilation cross section. Furthermore, we extended our numerical study in the case where γrays are injected in the form of a powerlaw. We found that quenching can be responsible for wide spectral breaks, which become more pronounced as the injected γrays become steeper (see Fig. 9).
Quenching has some interesting consequences when applied to existing γray sources. For example, for the 2006 MAGIC observations of the quasar 3C 279, modelling with combinations of low values of the Doppler factor δ and high values of the magnetic field strength B lead to strong absorption of the injected γrays. The only way to avoid quenching is either to assume typical values of δ (i.e., around 10) and low values of B or, alternatively, high values for both δ and B – see Fig. 12.
We restricted our analysis to cases where γrays are injected into the source without being too specific in assigning a particular radiation mechanism for producing them. We preferred this approach because our aim was to calculate the effects of quenching using simple functional forms for the γrays and not to model the source in detail. For example, we could have used proton synchrotron radiation to fit the 3C 279 TeV observations. However, this would have involved extra parameters for the radiating relativistic protons and would have complicated the analysis beyond the level at which we would like to present it here.
In our analytical treatment we ignored inverse Compton scattering, because its inclusion would have complicated the analysis. On the other hand, it was taken into account in all our numerical calculations, although it is greatly suppressed by KleinNishina effects (see Fig. 5). This process can still make an impact in cases where the soft photon compactness greatly exceeds the magnetic one (cf. Eq. (14)). However, this has to be seen in conjunction with the overall fitting of the MW spectrum of the source and we avoided presenting it here, although we took it into account in Sect. 4 when performing the analysis on quasar 3C 279.
Concluding, quenching has many farreaching implications for the modelling of compact γray sources such as Active Galactic Nuclei and Gamma Ray Bursts. Because quenching operates in an autoregulatory manner by coupling electrons to photons, its effects on the photon spectra cannot be manifested if the modelling is done by assuming an adhoc radiating electron and/or proton distribution. This could lead to parameter choices that lie within the quenching “forbidden” parameter space, making them in essence invalid.
On the other hand, if taken fully into account, quenching can give robust limits for important source parameters such as its size, magnetic field, and the Doppler factor for a given γray flux. In the case of 3C 279 for example, we found that to avoid quenching, the source should have a large Doppler factor (δ ≳ 10) whenever high B values are adopted. Similar constraints could be found for other AGNs with spectra extending to even higher energies (~1–10 TeV). In these cases the conditions for quenching are more stringent, because they critically depend on the maximum energy of the γray photons and are certainly worth investigating further.
We note that in the same framework an analogous study where ultrarelativistic protons are the primary injected particles has earlier been presented by Kirk & Mastichiadis (1992).
Acknowledgments
We would like to thank Dr. John Kirk for useful comments and discussion. 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 NSRF – Research Funding Program: Heracleitus II.
References
 Aleksić, J., Antonelli, L. A., Antoranz, P., et al. 2011, A&A, 530, A4 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Bonometto, S., & Rees, M. J. 1971, MNRAS, 152, 21 [NASA ADS] (In the text)
 Böttcher, M. 2007, Ap&SS, 309, 95 [NASA ADS] [CrossRef] (In the text)
 Böttcher, M., & Chiang, J. 2002, ApJ, 581, 127 [NASA ADS] [CrossRef] (In the text)
 Böttcher, M., Reimer, A., & Marscher, A. P. 2009, ApJ, 703, 1168 [NASA ADS] [CrossRef] (In the text)
 Coppi, P. S. 1992, MNRAS, 258, 657 [NASA ADS] [CrossRef] (In the text)
 Guilbert, P. W., Fabian, A. C., & Rees, M. J. 1983, MNRAS, 205, 593 [NASA ADS] (In the text)
 Herterich, K. 1974, Nature, 250, 311 [NASA ADS] [CrossRef] (In the text)
 Jelley, J. V. 1966, Nature, 211, 472 [NASA ADS] [CrossRef] (In the text)
 Kataoka, J., Takahashi, T., Makino, F., et al. 2000, ApJ, 528, 243 [NASA ADS] [CrossRef] (In the text)
 Katarzyński, K., Ghisellini, G., Tavecchio, F., et al. 2005, A&A, 433, 479 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Kazanas, D. 1984, ApJ, 287, 112 [NASA ADS] [CrossRef] (In the text)
 Kirk, J. G., & Mastichiadis, A. 1992, Nature, 360, 135 [NASA ADS] [CrossRef] (In the text)
 Konopelko, A., Mastichiadis, A., Kirk, J., de Jager, O. C., & Stecker, F. W. 2003, ApJ, 597, 851 [NASA ADS] [CrossRef] (In the text)
 MAGIC Collaboration, Albert, J., Aliu, E., Anderhub, H., et al. 2008, Science, 320, 1752 [NASA ADS] [CrossRef] [PubMed] (In the text)
 Mastichiadis, A., & Kirk, J. G. 1995, A&A, 295, 613 [NASA ADS] (In the text)
 Mastichiadis, A., & Kirk, J. G. 1997, A&A, 320, 19 [NASA ADS] (In the text)
 Petropoulou, M., & Mastichiadis, A. 2009, A&A, 507, 599 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Stawarz, Ł., & Kirk, J. G. 2007, ApJ, 661, L17 [NASA ADS] [CrossRef] (In the text)
 Stern, B. E., Begelman, M. C., Sikora, M., & Svensson, R. 1995, MNRAS, 272, 291 [NASA ADS] [CrossRef] (In the text)
 Svensson, R. 1987, MNRAS, 227, 403 [NASA ADS] (In the text)
 Zdziarski, A. A., & Lightman, A. P. 1985, ApJ, 294, L79 [NASA ADS] [CrossRef] (In the text)
All Figures
Fig. 1 Loglog plot of the critical luminosity compactness as a function of the injected γray photon energy ϵ_{γ}, when the δfunction approximation for the cross section σ_{γγ} was used for two values of the magnetic field strength: B = 40 G (solid line) and B = 4 G (dashed line). The values calculated by SK using only catastrophic losses are marked with crosses. 

Open with DEXTER  
In the text 
Fig. 2 Loglog plot of the critical luminosity compactness as a function of the injected γray photon energy ϵ_{γ}, when the Θfunction approximation for the cross section σ_{γγ} was used for two values of the magnetic field strength: B = 40 G (solid line) and B = 4 G (dashed line). 

Open with DEXTER  
In the text 
Fig. 3 Critical γray compactness as a function of ϵ_{γ} when (i) σ_{γγ} and synchrotron emissivity are approximated by a Θ and a δfunction respectively; and (ii) when for both quantities the full expressions are used. For the first case the numerical result (points) and the analytical one given by Eq. (34) (solid line) are shown. The full problem can only be treated numerically and the corresponding result is shown with a dashed line. The size of the source is R = 3 × 10^{16} cm and the magnetic field strength B = 3.57 G. 

Open with DEXTER  
In the text 
Fig. 4 Plot of the compactness of emerging hard (dashed), soft (dotted), and total radiated luminosity (full line) as a function of the injected compactness . The γrays are assumed to be monoenergetic, injected at energy ϵ_{γ} = 2.3 × 10^{4}. The size of the source is R = 3 × 10^{16} cm, while the magnetic field strength is B = 40 G. The emerging soft radiation is a result of the selfquenching of the system. 

Open with DEXTER  
In the text 
Fig. 5 Multiwavelength spectrum for a choice of such that . The other parameters are the same as in Fig. 4. Full line corresponds to the case where all relevant radiative processes were taken into account, while the dashed line corresponds to the case where only the two “core” processes were considered, i.e., synchrotron and γγ interactions. 

Open with DEXTER  
In the text 
Fig. 6 Time evolution of the soft (dotted line), hard (dashed line) and total l_{γ} (solid line) compactness for the same set of parameters as in Fig. 5. 

Open with DEXTER  
In the text 
Fig. 7 MW spectra for a powerlaw injection of γrays of the form for ϵ_{γmin} < ϵ_{γ} < ϵ_{γmax} with ϵ_{γmin} = 23 and ϵ_{γmax} = 2.3 × 10^{4} for various values of the normalization Q_{0}. This starts from Q_{0} = 4 × 10^{3} and then exceeds its previous value by a factor of 2. The dashed line curves represent the runs without quenching. The source parameters are B = 40 G and R = 3 × 10^{16} cm. 

Open with DEXTER  
In the text 
Fig. 8 Same as in Fig. 7 with ϵ_{γmin} = 2.3 × 10^{2} and ϵ_{γmax} = 2.3 × 10^{5} with Q_{0} starting from Q_{0} = 2.5 × 10^{4}. 

Open with DEXTER  
In the text 
Fig. 9 Slope of the quenched γray spectrum as a function of the injected spectrum’s slope (solid line). To facilitate a comparison the line α_{inj} + 1 is also shown (dashed line). 

Open with DEXTER  
In the text 
Fig. 10 MW spectra of 3C 279 for R = 10^{16} cm, B = 40 G and δ = 10^{0.8} (dotted line), 10^{1.0} (dashed line), 10^{1.2} (dotdashed line) and 10^{1.4} (solid line). In all cases the γray injection would have produced the solid line in the absence of quenching. Black squares and bowtie represent the observational data from February 2006. 

Open with DEXTER  
In the text 
Fig. 11 MW spectra of 3C 279 for R = 10^{16} cm, δ = 10^{1.2} and B = .16 G (solid line), 1.6 G (dashed line) and 16 G (dotted line). In all cases the γray injection would have produced the solid line in the absence of quenching. Black squares and bowtie represent the observational data from February 2006. 

Open with DEXTER  
In the text 
Fig. 12 Parameter space of the allowed values of the magnetic field B as a function of the size of the source R for Doppler factors of δ = 10 (dotted line), 15.8 (dashed line), and 25 (solid line). The space above a line for a certain δ corresponds to combinations of R and B that cannot successfully fit the TeV observations of 3C 279. 

Open with DEXTER  
In the text 