Issue |
A&A
Volume 526, February 2011
|
|
---|---|---|
Article Number | A61 | |
Number of page(s) | 6 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/200913267 | |
Published online | 22 December 2010 |
TeV gamma ray opacity in PKS 2155-304
1
Research School of Astronomy and Astrophysics, Australian National
University,
Cotter Rd. Weston,
ACT
2611,
Australia
e-mail: geoff@mso.anu.edu.au
2
Landessternwarte Heidelberg,
Königstuhl 12,
69117
Heidelberg,
Germany
e-mail: S.Wagner@lsw.uni-heidelberg.de
Received:
8
September
2009
Accepted:
4
November
2010
Context. Recently observed rapid time variability in the gamma ray emission from the blazar PKS 2155-304 indicates that either the site of gamma ray emission is close to the black hole or the jet has a very high Lorentz factor. Consideration of the opacity of gamma rays close to the black hole provides additional information related to these two possibilities.
Aims. We investigate the TeV gamma ray opacity resulting from pair production on soft photons radiated by the black hole accretion disk, considering situations where the radiation produced by the disk is close to that predicted by the Shakura-Sunyaev model and situations where the disk radiation is much less.
Methods. We utilise expressions for the pair opacity of very high energy gamma rays developed by Gould and Schréder and Donea and Protheroe and use expressions for the disk radiation field produced by an accretion disk external to a Poynting flux dominated jet. A lower level of disk emission is modelled by using the spectrum corresponding to a lower accretion rate. We also consider changes in the shape of the gamma ray sepctrum as a flare emerges from the optically thick region.
Results. If the gravitational power resulting from accretion is dissipated by radiation in the region of the disk external to the jet, then TeV gamma rays are opaque out to of order 40 gravitational radii ≈6 × 1015 cm from the black hole. However, the opacity is considerably reduced if the disk radiates a negligible amount of gravitational power as would be the case if there is a significant disk wind external to the jet. In either case the absence of spectral changes during observed flares excludes scenarios in which the TeV gamma rays are emitted even in regions of modest pair opacity.
Key words: black hole physics / accretion, accretion disks / galaxies: jets / gamma rays: theory / BL Lacertae objects: general
© ESO, 2010
1. Introduction
TeV gamma ray observations of the blazar PKS 2155-304 have revealed rapid time variability
on a scale of Δt ≈ five minutes (Aharonian
et al. 2007). For a Doppler factor of δ this implies a spatial
scale, D, in the frame of the emitting plasma given by:
(1)With the gravitational
radius
rg = GM/c2
one finds
(2)assuming
a black hole mass ≈109 M⊙ (Aharonian et al. 2007). Thus, there are two immediate
possibilities: (1) if the Doppler factor is typical of that estimated for a number of other
gamma ray blazars, or typical of that inferred on the VLBI scale,
i.e. δ ≈ 10, then the size of the emitting region is of order the
gravitational radius of the black hole. This naturally places the emission region close to
the black hole. (2) If the diameter of the emitting region is of order
1015−16 cm, i.e. typical of estimates based upon the variability of other
blazars and modelling of their emission regions, then one requires δ ≳ 100
if the size of the emitting region is not to be too small a fraction of the jet diameter.
These broadly indicative estimates highlight the interest that is raised by the rapid
variability detected in very high energy gamma ray observations of PKS 2155-304.
Gamma rays are subject to absorption by pair creation on soft photons and this process significantly influences the region from which gamma ray emission can be detected. The presence of an intense radiation field prohibits the escape of gamma rays. Therefore, it is in principle possible to constrain the emission site of the TeV emission in blazars, and in PKS 2155-304 in particular, by taking pair opacity into account.
The threshold energy ϵth of a soft photon that pair produces
with a gamma ray of energy E when the angle between the directions of the
two photons is θ, is (3)(Gould & Schréder 1979). The minimum
threshold energy is
(4)for
a head-on collision (θ = π) so that TeV opacity is
primarily influenced by the infrared radiation field. It is well known that absorption of
TeV gamma rays by the diffuse infrared background limits the redshift from which
TeV emission can be detected. Here we are concerned with the infrared radiation density in
the vicinity of the black hole.
Gamma ray pair opacity in AGN has been considered previously by Blandford & Levinson (1995), who introduced the concept of the gamma-sphere and by Donea & Protheroe (2003). Blandford & Levinson (1995) primarily considered the effect of absorption in the EGRET energy band (≲10 GeV) in which case the soft photon energy exceeds ≈30 eV. Hence, their input soft photon spectrum is dominated by hard and soft X-rays, either produced within the jet or by scattering of coronal and disk photons in the broad line region of the AGN. They adopted empirical forms for the input spectra. They also considered the cascading production of lower energy gamma rays and electron-positron pairs. Donea & Protheroe (2003), on the other hand, were concerned with absorption in the TeV energy range and considered the effects of pair opacity resulting from infrared photons originating from the accretion disk, the broad line region and a possible torus. Since they considered much less energetic soft photon inputs, their disk photon spectrum was determined theoretically using a version of the “jet-disk symbiosis” model (Donea & Biermann 1996). Both Blandford & Levinson (1995) and Donea & Protheroe (2003) developed models for a range of black hole masses and input spectra and consider the effect of absorption on a steady emitter of gamma rays.
In this paper we consider in detail variations in the spectral index of the gamma rays as an emitting region in the jet emerges from the optically thick region. We follow the approach of Donea & Protheroe (2003) but primarily concentrate on the effect of disk emission on pair opacity. Our disk emission model differs from Donea & Protheroe (2003) and is based on a Poynting flux-dominated jet with a jet radius of a few gravitational radii together with a self-consistent Shakura-Sunyaev spectrum external to the jet radius (see Sect. 3). The parameters in our calculations are most directly relevant to PKS 2155-304 but are also relevant to other gamma ray emitting jets since we consider a range of accretion rates.
2. Gamma ray absorption in the vicinity of accretion disks
We consider the opacity of gamma rays escaping along the jet perpendicular to the disk
towards an observer in pole-on geometry, as is generally assumed for Blazars. Referring to
Fig. 1, the rate of change of optical
depth τ with height z above the accretion disk
(the inverse absorption length) is given by Eq. (8) below. We first define the symbols and terms in this expression:
E is the gamma ray energy, θ is the angle between the
direction of the soft photon and that of the gamma ray photon,
μ = cosθ,
μmin,max = cosθmax,min
where θmax,min are respectively the maximum and
minimum values of θ,
ϵmin(E) = ϵ0 is
the minimum threshold energy (see Eq. (3)),
dn/dϵdΩ is the number density of soft
photons per unit energy per unit solid angle. The square of the electron-positron Lorentz
factor in the centre of mass frame is s and
β = (1−1/s)1/2
is the corresponding velocity; r0 is the electron radius. The
pair production cross
section σγγ(β) is:
(Gould & Schréder 1979). The rate of change of
optical depth with distance is:
(8)(Gould & Schréder 1979; see also Donea & Protheroe 2003).
![]() |
Fig. 1 The geometry of gamma ray absorption resulting from pair production in the vicinity of an accretion disk. |
Approximating the accretion disk spectrum at each radius by a local black body, the number
density of soft photons per unit energy per steradian is: (9)where
T(r) is the accretion disk temperature at a radial
distance r = ztanθ.
Following Gould & Schréder (1979) we
evaluate the integral (8) utilising the
normalized energy variable
η = ϵ/ϵ0
(where ϵ0 is given by Eq. (4)) and the variable
s = s(μ) defined by Eq. (7). The result is (10)where
α is the fine structure constant,
λe = h/mec
is the electron Compton wavelength, and the limits on the integrals are given by:
(11)Once
dτ/dz has been numerically evaluated,
the optical depth between any two points along the z axis may be evaluated
by a further numerical integration with respect to z.
Relativistic leptons generated via pair-absorption upscatter ambient radiation up to gamma-ray energies resulting in a second-order contribution to the observed gamma-ray spectra. This contribution is negligible for the considerations in this study and is not taken into account here.
3. Spectral energy distribution of soft disk photons
An accretion disk model is required in order to determine the temperature appearing in Eq. (10). The standard model is the Shakura & Sunyaev (1973) model in which the gravitational power is entirely dissipated in the disk and appears as (approximately) black body radiation from the surface. However, the mechanical power in winds and jets reduces the flux of radiation from the surface (Kuncic & Bicknell 2004). Moreover, in the case of AGN some of the gravitational power not channeled through outflows is not dissipated in the disk but in the corona in a way that is not currently understood, notwithstanding the insights provided by partly phenomenological models (Haardt & Maraschi 1991, 1993; Haardt et al. 1994). In view of our imperfect understanding of the physics of coronal emission, it is not clear for blazars such as PKS 2155-304, how much of the gravitational power is channelled into disk and coronal emission, respectively.
In view of the uncertainties in the role of coronal emission and hence in the accretion disk spectral energy distribution we ignore potential dissipation in the corona and adopt the following approach:
-
1.
We assume that the base of the jet occupies a radial rangebetween the minimally stable circularorbit rms and the outer jet radius rjet. We adopt the limiting case (Bicknell & Li 2007) in which all of the angular momentum transported inside r = rjet is transported vertically by a Poynting flux-dominated jet. Consistent with this assumption, most of the accretion power liberated between rjet and rms is channeled into the kinetic power of two jets and most of the accreted mass goes into the black hole. In this limit, the power of each jet and the accretion rate, Ma, are related by the following equation (for a Newtonian accretion disk; see Bicknell & Li 2007):
(12)In the models presented below we estimate the gamma ray opacity for specific values of the parameters rms and rjet; we discuss the effect of varying these parameters further in Sect. 8.
-
2.
Equation (12) may be used to determine the mass accretion rate given the jet power. We adopt a fiducial power ≈1044 erg s-1 as implied by the analysis of Begelman et al. (2008). Adopting a black hole mass M ≈ 109 M⊙ for PKS 2155-304 (Aharonian et al. 2007), an Eddington luminosity, LEdd ≈ 1.3 × 1047 erg s-1 and defining the normalized accretion rate
gives:
(13)Thus, for this Newtonian accretion disk model, we estimate ṁ ≈ 10-2 for rms/rg = 2 and rjet/rg = 5.
-
3.
Previous models of PKS 2155-304 (e.g. Aharonian et al. 2005) imply a lower jet power ≈1043 erg s-1 so that we also consider lower values of ṁ in our models.
-
4.
If, in the region of the disk external to the jet, there is no energetically significant wind, then the accretion disk temperature for r > rjet is given by a Shakura-Sunyaev model in which the turbulent stress is zero at r = rjet. Hence, the photon number density per unit energy per unit solid angle is as given by Eq. (9) with the accretion disk temperature for r > rjet given by:
(14)where σ is the Stefan-Boltzmann constant.
-
5.
As well as allowing for lower jet powers, models with lower ṁ give us some insight into the situation where there is an energetically significant wind, which reduces the radiative flux from the disk (Kuncic & Bicknell 2004). We therefore consider a range of values 10-5 ≤ ṁ ≤ 10-2.
-
6.
We neglect gravitational effects on the energy of photons and their trajectories close to the black hole. This is discussed further below in Sect. 7
In this model the ratio of the disk power to total power is rms/rjet. Hence, for jet radii of a few gravitational radii the radiation field implied by Eq. (14) may be a factor of a few lower in intensity than expected from the standard Shakura & Sunyaev (1973) model. Nevertheless, it is not realistic to expect that the radiative intensity is 2−3 orders of magnitude less than that from a Shakura-Sunyaev disk on the basis of the kinetic luminosity of the jet alone. The most direct way in which the disk radiative intensity could be reduced would be if most of the gravitational power outside the jet were to be channeled into a slower disk wind. (See Sect. 4.)
The above ansatz for calculating the accretion disk spectrum differs from Donea & Protheroe (2003) who based their calculation on a non-magnetic disk and jet model (“jet symbiosis model”) in which the flux of angular momentum is proportional to the jet mass flux. This aspect is essential since there is no magnetic field to transport angular momentum. The mass flux in their model is small in order that the jet have a high asymptotic velocity. Hence, little angular momentum is transported by the jet and their disk model is close to the Shakura-Sunyaev model. In our model there is a relative deficit of radiation from the region of the disk interior to rjet, as illustrated in Fig. 4.
4. Estimation of opacity
![]() |
Fig. 2 The pair opacity of gamma rays for different energies and different values of the
normalised accretion |
Consider first the case of an accretion disk corresponding to a normalized accretion rate ṁ = 10-2. The gamma ray opacity between a height z above the accretion disk and 500 gravitational radii is shown in the panel (a) of Fig. 2 for a number of different energies below and above a TeV. For TeV gamma rays the region is optically thick out to ≈100 gravitational radii ≈1.5 × 1016 cm. Moreover the height of the gammasphere increases with energy; 10 TeV emission remains optically thick out to approximately 500 gravitational radii.
The optical depths for the lower jet power of 1043 erg s-1 corresponding to ṁ ≈ 10-3 (panel (b) of Fig. 2, display a similar trend, although there is a noticeable decrease in optical depth. At a TeV, for example, the emission starts to become optically thin at around 30 gravitational radii ≈4.5 × 1015 cm.
The most direct way to decrease the opacity is to reduce the level of disk radiation. This occurs if there is a disk wind in the region outside the jet thereby reducing the gravitational power that is dissipated in disk radiation. The details depend on the disk spectrum, which in turn depends on the dependence of wind mass-loss rate on radius. Here, we obtain a semi-quantitative estimate of the reduction of the opacity by simply reducing ṁ. Thus, the lower two panels of Fig. 2 indicate the opacities corresponding to ṁ = 10-4 and 10-5. Nevertheless, for ṁ = 10-5 the 10 TeV emission is still optically thick out to about 40 gravitational radii ≈6 × 1015 cm.
Note that for the two highest accretion rates, the 1012 eV and 1012.5 eV opacity curves cross at a height of about 20 rg. This occurs for the following reasons: at small heights above the accretion disk, the photon field is dominated by the inner part of the accretion disk. For the higher accretion rates the threshold energy for disk photons that pair produce with 1012.5 eV gamma rays is considerably less than the mean energy of photons emitted from that region of the disk. The decrease of the pair cross section with energy means that these disk photons are less effective in producing pairs than a comparable number density of lower energy photons. Moreover the threshold energies of photons that pair-produce with 1012 eV photons are closer to the mean disk photon energy. This is why their opacity in this region is slightly larger. As the emission height increases a larger region of the disk contributes to the opacity and the effect of the energetic photons from the inner disk is minimised. This effect disappears at the lower accretion rates since the inner disk is not as hot for those photons.
5. Spectral variations
The intrinsic spectrum is likely to exhibit a cutoff at high energies. Hence, it is not easy to unambiguously identify pair-absorption in single epoch spectra even for significant opacities.
Figure 2 shows that at a given energy the opacity is dependent on the height above the accretion disk. Therefore, a gamma-ray emitting region emerging from an optically thick region becomes optically thin at progressively higher energies. That is, its spectral index changes with time. For an emitting region that is close enough to the central black hole at the onset of a flare to be tranversing the regime of moderate absorption while the flare is ongoing, such spectral changes can be observed. We consider the case of the flares in PKS 2155-304 reported by (Aharonian et al. 2007). No spectral changes were observed during the 120 min of continuous observations.
In order to indicate how the spectrum may evolve observationally, we proceed as follows:
let Nobs(E,z0)
(in erg cm-2 s-1 erg-1) be the observed spectrum at a
height z0 above which absorption is negligible,
Nabs(E,z) be the observed spectrum at a lower
height z
and Δτ(E,z,z0) the
difference in optical depth between the two heights. Then, if the intrinsic spectrum of the
flare does not change during the time of propagation from z
to z0, (15)The differential
optical depth between the two heights is given by
(16)and
dτ/dz is given by Eq. (10).
![]() |
Fig. 3 The observed high energy gamma ray spectrum of PKS 2155-304 modified by the optical depth at three different heights, assuming for illustration that the observed spectrum is emitted at a height of 300 rg above the accretion disk. The filled circles represent the observed spectrum. The various symbols represent the absorbed spectrum at the following heights: open circles z = 50 rg; Diamonds: z = 100 rg; Crosses: z = 150 rg. |
By way of illustration we select the optical depth derived for the normalised accretion rate ṁ = 10-3. The optical depth is unity at 10 TeV for a height z = 304.3 so that the spectrum of a flare observed at this height could be assumed to be close to intrinsic. Thus, we take z0 = 300 rg and evaluate the absorbed spectrum according to Eq. (15) for heights z = 50,100 and 150 rg. The template spectrum is the integrated spectrum presented by (Aharonian et al. 2007), chosen to indicate realistic error bars and deviations from a power-law slope. The results of applying the relevant differential optical depths to the observed spectrum are shown in Fig. 3. As can be seen in Fig. 3 the major temporal changes in the spectrum occur for energies greater than a TeV and naturally the greatest changes are observed for the radiation emitted from closest to the black hole. The observed time for a flare to propagate from z = 50 rg to 100 rg is 2.9/δ days. The smallest changes in the observed spectrum in this example result from the propagation from 150 to 300 rg; this would take 8.65/δ days. Therefore, the absence of an observable spectral change rules out stratified absorption if the flare propagates without much change in its intrinsic spectrum over this period. It would be contrived if intrinsic spectral changes would compensate the spectral variations resulting from varying amounts of pair-absorption.
6. Discussion of some physical aspects of the models
![]() |
Fig. 4 Comparison of threshold energy and mean emitted energy for two accretion disk models. The mean energies of photons emitted by the disk at each radius are indicated by the solid curves. The curve labelled “a” represents the mean photon energy from a Shakura-Sunyaev accretion disk with parameters Mb = 109 M⊙, ṁ = 10-2, rjet/rg = 2; the curve labeled “b” represents the mean energy from a modified disk with the same parameters except that rjet/rg = 5. The dot-dash curves to the upper right show the threshold energy for disk photons emitted from a given radius interacting with 100 GeV gamma-rays. The dashed curves represent the threshold energy for disk photons interacting with 1 TeV gamma rays. The different threshold curves are derived from gamma-ray emission heights of 10, 25, 50 and 100 gravitational radii respectively, as indicated on the curves. |
It is instructive to compare our disk model with the Shakura-Sunyaev model by comparing the threshold energy for absorption by photons emitted from a given disk radius with the mean energy emitted at that radius. Since the threshold energy depends inversely on the angle between the two photons (see Eq. (3)), it increases as the radius from which the low energy photons are emitted, decreases. Figure 4 shows the comparison of mean emitted energy with threshold energy as a function of disk radius for two accretion disk models. The first is a Shakura-Sunyaev model with black hole mass, Mbh = 109 M⊙, normalized accretion rate, ṁ = 10-2, no jet and rms/rg = 2; the second model has the same Mbh, ṁ and rms but has a jet with rjet/rg = 5. The decrease in disk emission and the deficit of high energy photons of the second model are apparent. We also show in this plot the threshold energy as a function of disk radius for photons interacting with (a) 1 TeV and (b) 100 GeV photons at heights of 10−100 gravitational radii. Clearly, the interaction of disk photons with 100 GeV gamma-rays is almost insignificant irrespective of the particular disk model. However, it is clear from the comparison of threshold energies with the mean disk photon energy that there is substantial interaction between disk photons and 1 TeV gamma-rays. However, the modified disk model decreases the interaction from photons emitted close to the jet since the emission is truncated for r < 5 rg, the mean disk photon energy decreases substantially near r = 5 rg and the threshold energy also rises sharply in this region.
The comparison of mean emitted and threshold energy in Fig. 4 provides additional justification for the neglect of gravitational redshift. For heights of the gamma ray emitting region in excess of 25 rg the mean disk photon energy is less than the threshold inward of about 15 rg indicating a reduction in the contribution to the opacity of photons emitted close to the black hole. For gamma ray emission heights ≈10 rg from the black hole there will be some effect from the gravitational redshift. However, our estimates of opacity in this region are much better than “factor of order unity” and the trends discussed in the following section are robust.
7. Relativistic effects
We have used a Newtonian disk model, assuming that the details of disk emission external to ≈5 rg would not make a large difference to our estimates of opacity. In this section we estimate effects related to general relativity.
There are three principal effects to consider: (1) the redshift of the emitted gamma rays. (2) The change in energy of a photon along its trajectory from the disk surface to the gamma ray emitting region along the axis of symmetry of the black hole. (3) The change in the locations of the footpoints (in the disk) of photons which intercept the symmetry axis. (4) The disk spectrum modified by relativistic effects. Here we consider the first three of these with (3) being particularly important.
We consider relativistic effects in the following way. We first consider the equations for
particle trajectories in the Kerr metric given by Bardeen
et al. (1972) specifically for the case of zero rest mass. Moroever, we consider
photons which intersect the symmetry axis; these have zero angular
momentum, L. The equations for the photon trajectories in Boyer-Lindquist
coordinates, (r,Θ,φ,t), where r is in
units of gravitational radii and t is in units
of GM/c3 are:
where
In
order to determine the constant of the motion discovered by Carter (1968), which enters into the above equations, we consider photons which
intercept the axis of symmetry at the location
r = r0 at the angle θ
indicated in Fig. 1. In this case,
(26)
![]() |
Fig. 5 Trajectories in Boyer-Lindquist coordinates of photons intercepting the symmetry axis of a maximally rotating black hole. The event horizon and static limit surface of the black hole are indicated. |
We integrate the above equations backwards from the axis of symmetry until they either
intercept the accretion disk (θ = π/2)
or the static limit surface at .
A typical set of trajectories for a maximally rotating black hole is shown in Fig. 5. In the figure it is apparent that photons are affected
by the gravitational field: all trajectories beginning at the accretion disk curve downwards
as they approach the vicinity of the black hole; trajectories beginning near the black hole
(within about 5 rg) and initially traveling outwards are
pulled back towards the axis. The effect of this is to focus the photons propagating towards
the axis, increasing the photon density there. Another way of looking at this is that the
footpoints of photons intercepting the axis at a given angle are moved inwards compared to
where they would be if there were no gravitational field. The relative movement is greatest
for the trajectories which begin furthest out from the black hole. Relativistic effects will
therefore increase pair-absorption.
In order to semi-quantitatively estimate the effect of the enhancement of photon density
resulting from the gravitational field, we consider the integrand appearing in the
expression (8) for the absorption length.
The absorption cross-section peaks at s ≈ 2 so that for a specific angle
θ = cos-1μ the dominant contribution from the
photon density occurs at an energy
ϵ = 4(mec2)2E-1(1 − μ)-1,
i.e. at twice the threshold energy. We have evaluated the photon density contribution at
this energy using the Newtonian model adopted in this paper and also by taking into account
the shifted footpoints of the photon trajectories discussed above but still using the
Newtonian disk model. The disk temperature is higher at each shifted footpoint with a
consequent shift on the blackbody spectrum. In addition, if ϵ is the photon
energy on axis and ϵ0 is the the photon energy at the disk
surface, then the photon density on axis is related to the photon density at the disk by:
(27)All energies are
calculated in the local non-rotating frame (LNRF) of Bardeen
(1970).
This procedure gives us an approximate idea of the differences between our Newtonian model and a fully relativistic treatment. As one expects the errors are largest for gamma ray emission regions closest to the black hole – as large as 30 percent for z = 10 rg and decreasing slowly for points at larger distances. At z = 100 rg the estimated error is about 15 percent; at z = 500 rg there is about a 5 percent effect. These estimates confirm that the optical depths calculated here are reasonably accurate and that the qualitative conclusions concerning the dependence of optical depth on distance from the black hole and the effect of the radiation field of the disk are robust.
8. Discussion and conclusions
For the parameters that we have used in this paper, we have shown that if the disk external to the jet in PKS 2155-304 dissipates a substantial portion of the gravitational power resulting from accretion in radiation, as indicated by our truncated Shakura & Sunyaev (1973) model, then the pair opacity close to the black hole is significant at ≳TeV gamma-ray energies and possibly lower, depending on the precise value of the normalized accretion rate. For a modest value of the normalised accretion rate, ṁ = 10-3 corresponding to a jet power ≈1043 erg s-1, the opaque region at an energy of 1 TeV extends to approximately 400 gravitational radii ≈6 × 1016 cm. Within the region pair-absorption is expected to be detectable especially through spectral monitoring. The absence of spectral changes during outbursts provides constraints against even modest amounts of pair-absorption.
How does this conclusion depend upon parameters such as the jet radius, which we have not varied in this paper and also the assumption that the jet is Poynting flux dominated? The TeV optical depth is high for rjet = 5 rg and ṁ ≳ 10-3 and any reduction in jet radius would make the optical depth even larger since there would be more radiation illuminating the jet. An increase in jet radius decreases the ambient radiation density. However, even for rjet = 10 rg the optical depth for energies of 1012 and 1012.5 eV is unity at heights of 37 and 124 rg respectively. A jet radius of 10 rg is larger than is usually contemplated for jets in active galaxies. Hence our conclusion that for ṁ ≳ 10-3 the very high energy gamma ray emission (say E ≳ 1012 eV) is optically thick out to approximately 100 gravitational radii is secure.
The major source of uncertainty in the calculations is the spectrum of the accretion disk. This is not only influenced by relativistic effects but also by the effect of potential disk winds outside the jet region (Kuncic & Bicknell 2007). We have assumed that the jet is disk-driven as proposed by Blandford & Payne (1982) and investigated by numerous workers since. If the jet is driven by the process proposed by Blandford & Znajek (1977), then in that case the disk exists at small radii (certainly within 5 rg) and again the ambient radiation field exceeds what we have estimated here.
We have also investigated the alternative case in this paper, namely that the disk luminosity is well below that implied by the Shakura & Sunyaev (1973) model. We have noted that this would be the case if a disk wind existed which advected energy away from the disk, reducing the dissipation of gravitational power. A comparison with M 87 is useful: taking into account the revised mass estimate of M 87’s black hole, Mbh ≈ 6.4 × 109 M⊙ (Gebhardt & Thomas 2009), the ratio of disk luminosity to Eddington luminosity is approximately 9 × 10-7 (Whysong & Antonucci 2004; Jolley & Kuncic 2007). In the models we have used in this paper, the ratio of disk luminosity to Eddington luminosity is (rjet/rg)ṁ = 0.2 ṁ. For the lowest values of ṁ = 10-5, the ratio is comparable to that inferred for M 87. Thus there is good reason to suppose that the disk in PKS 2155-304 may not radiate at anything like the Shakura-Sunyaev rate and that the gamma ray emission originates from close to the accretion disk. For ṁ = 10-5 spectral changes at 1 TeV would only be detectable if the emitting region is closer than ≈30 rg.
Acknowledgments
We gratefully acknowledge support from the Go8-DAAD scheme which enabled this research to be carried out. We also acknowledge valuable support from Sonderforschungbereich 439.
References
- Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2005, A&A, 442, 895 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2007, ApJ, 664, L71 [Google Scholar]
- Aharonian, F., Akhperjanian, A. G., Anton, G., et al. 2009, ApJ, 696, L150 [NASA ADS] [CrossRef] [Google Scholar]
- Bardeen, J. M. 1970, ApJ, 162, 71 [NASA ADS] [CrossRef] [Google Scholar]
- Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347 [NASA ADS] [CrossRef] [Google Scholar]
- Begelman, M. C., Fabian, A. C., & Rees, M. J. 2008, MNRAS, 384, L19 [NASA ADS] [CrossRef] [Google Scholar]
- Bicknell, G. V., & Li, J. 2007, ApSS, 311, 275 [NASA ADS] [Google Scholar]
- Blandford, R. D., & Levinson, A. 1995, ApJ, 441, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
- Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 [NASA ADS] [CrossRef] [Google Scholar]
- Carter, B. 1968, Phys. Rev., 174, 1559 [NASA ADS] [CrossRef] [Google Scholar]
- Donea, A.-C., & Biermann, P. L. 1996, A&A, 316, 43 [NASA ADS] [Google Scholar]
- Donea, A.-C., & Protheroe, R. J. 2003, Astrop. Phys., 18, 377 [NASA ADS] [CrossRef] [Google Scholar]
- Gebhardt, K., & Thomas, J. 2009, ApJ, 700, 1690 [NASA ADS] [CrossRef] [Google Scholar]
- Gould, R. J., & Schréder, G. P. 1979, Phys. Rev., 155, 306 [Google Scholar]
- Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 [NASA ADS] [CrossRef] [Google Scholar]
- Haardt, F., & Maraschi, L. 1993, ApJ, 413, 507 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Haardt, F., Maraschi, L., & Ghisellini, G. 1994, ApJ, 432, L95 [NASA ADS] [CrossRef] [Google Scholar]
- Jolley, E. J. D., & Kuncic, Z. 2007, ApSS, 311, 257 [NASA ADS] [Google Scholar]
- Kuncic, Z., & Bicknell, G. V. 2004, ApJ, 616, 669 [NASA ADS] [CrossRef] [Google Scholar]
- Kuncic, Z., & Bicknell, G. V. 2007, Mod. Phys. Lett. A, 22, 1685 [NASA ADS] [CrossRef] [Google Scholar]
- Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
- Whysong, D., & Antonucci, R. 2004, ApJ, 602, 116 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
![]() |
Fig. 1 The geometry of gamma ray absorption resulting from pair production in the vicinity of an accretion disk. |
In the text |
![]() |
Fig. 2 The pair opacity of gamma rays for different energies and different values of the
normalised accretion |
In the text |
![]() |
Fig. 3 The observed high energy gamma ray spectrum of PKS 2155-304 modified by the optical depth at three different heights, assuming for illustration that the observed spectrum is emitted at a height of 300 rg above the accretion disk. The filled circles represent the observed spectrum. The various symbols represent the absorbed spectrum at the following heights: open circles z = 50 rg; Diamonds: z = 100 rg; Crosses: z = 150 rg. |
In the text |
![]() |
Fig. 4 Comparison of threshold energy and mean emitted energy for two accretion disk models. The mean energies of photons emitted by the disk at each radius are indicated by the solid curves. The curve labelled “a” represents the mean photon energy from a Shakura-Sunyaev accretion disk with parameters Mb = 109 M⊙, ṁ = 10-2, rjet/rg = 2; the curve labeled “b” represents the mean energy from a modified disk with the same parameters except that rjet/rg = 5. The dot-dash curves to the upper right show the threshold energy for disk photons emitted from a given radius interacting with 100 GeV gamma-rays. The dashed curves represent the threshold energy for disk photons interacting with 1 TeV gamma rays. The different threshold curves are derived from gamma-ray emission heights of 10, 25, 50 and 100 gravitational radii respectively, as indicated on the curves. |
In the text |
![]() |
Fig. 5 Trajectories in Boyer-Lindquist coordinates of photons intercepting the symmetry axis of a maximally rotating black hole. The event horizon and static limit surface of the black hole are indicated. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.