Arxiv:0712.2902V2

arXiv:0712.2902v2 [astro-ph] 27 May 2008 rls,osre nsvrlAN mle eycmatVHE most compact at very of a implies region AGNs emission several in days of observed scales less, time on or flux, gamma-ray VHE varying rapidly The Introduction 1. h H hrceitc,icuigvraiiyo iesca time on variability reproduce to including ∆ enough characteristics, flexible VHE be may the single models conventional, whether SSC shown be mogeneous to remains it and ing, enwdl icse,bt ntecneto usremissio pulsar of context has the magnetospheres in both rotating discussed, by widely Centrifuga flows been pro- emission. plasma Compton gamma-ray of inverse VHE acceleration variable via to could lead system, cesses sub-Edd hole a of black vicinity the accreting centrifugal in that occurring electrons, possibility of the eration explore we Here 2006). al. l 07 lete l 07,aegnrlydi et generally Aharonian are 1996; 2007), al. al. et et (Gaidos Albert minutes variation 2007; of blazar al. VHE scales observed time fastest the Mastichi on Yet, (e.g., 2002). and photons a Kirk electrons soft & successfully of of be upscattering acceleration Compton shock principle, verse internal in by may, for hours counted several of scales in eg,Slaie l 98 igr20) ntecs of case the In 2004). with expl Rieger AGNs 1998; their aligned al. likely for et are Salvati set-ups (e.g., and geometrical tions models non-standard emission require VHE to jet standard within stand m where cm, ein o lzrsucswt hi espitn lotd almost pointing jets their (i.e., us with towards sources rectly blazar For region. edo Send coe 1 2018 31, October srnm Astrophysics & Astronomy t ∼ < a,a bevdi h ai aayM7(hrna et (Aharonian M87 galaxy radio the in observed as day, 2 eevd..07 cetd... accepted ...2007; Received 3 2 1 e words. Key ( accreting advection-dominated some that suggest results Conclusions. black central the onto Results. accretion assuming M87, det galaxy are radio scale, time variability characteristic associated Methods. Aims. Context. to t conventional for challenging proves M87, from reported as opoie ikpoosrslsi a e spectrum TeV flat a in results photons disk Comptonized cl so h order the of is scale nlz h nrycntansipsdfrhgl underlum highly magneto for jet imposed constraints rotating energy by the analyze particles of acceleration trifugal ff rn eussto requests print ulnIsiuefrAvne tde,3 izila Plac Fitzwilliam 31 Studies, Advanced for Institute Dublin uoenAscae aoaoyfrGmaRyAtooy j Astronomy, Gamma-Ray for Laboratory Associated European a-lnkIsiu ¨rKrpyi,Sufrhcwg1, Saupfercheckweg f¨ur Kernphysik, Max-Planck-Institut γ aibeVEgmaryeiso rmnnbaa AGNs non-blazar from emission gamma-ray VHE Variable ∼ oiae yrcn oko usrtp atceaccelerat particle pulsar-type on work recent by Motivated δ (10 eso htfrahgl nelmnu orelk 8,cent M87, like source underluminous highly a for that show We stebl ope atro h H emitting VHE the of factor Doppler bulk the is h bevto frpdyvral eyhg nry(H)g (VHE) energy high very variable rapidly of observation The h aiu oet atrfrcnrfglyaclrtde accelerated centrifugally for factor Lorentz maximum The 7 − aais cie–glxe:jt aito ehns:no mechanism: radiation – jets galaxies: – active galaxies: etiua ceeaincudpoieantrlexplanat natural a provide could acceleration Centrifugal 10 δ 8 ,alwn nes opo Tosn pcteigo sub- of upscattering (Thomson) Compton inverse allowing ), ..Rieger F.M. : ∼ δ r L ,ltl a egie rmjtboost- jet from gained be can little 1, ∼ / c hc ntecs fM7crepnsto corresponds M87 of case the in which , Γ R aucitn.aal no. manuscript b ≤ ∼ c 5,VEvraiiyo time on variability VHE 15), ∆ t δ = 2 . ..Rieger F.M. 6 × ffi utt under- to cult 10 15 L rie.Apiain r rsne o odtosexpect conditions for presented are Applications ermined. ( ∆ te othe to etter t / pee ntevcnt faceigsprasv lc hol black supermassive accreting of vicinity the in spheres ington o-lzr Gscudwl eosral H mtigsour emitting VHE observable be well could AGNs non-blazar) accel- L 1 d) 1 , non- ν ABSTRACT nu systems. inous ana- 2 91 edleg emn;e-mail: Germany; Heidelberg, 69117 ertclaclrto n msinmodels. emission and acceleration heoretical adis ho- ,Dbi ,Ireland 4, Dublin e, oet cu na deto-oiae AA)mode. (ADAF) advection-dominated an in occur to hole ∝ les in- n ..Aharonian F.A. and c- s, i- n δ ν l − o nM7(eoo hrna 07,w eeaietecen- the re-examine we 2007), Aharonian & (Neronov M87 in ion ityspotdb NSadMPG and CNRS by supported ointly α ≃ c l 96 agdaa&Lsh19;R0) h ailparticl radial et The RM00). Machabeli 1997; 1994; Lesch Rogava & Gangadhara & o 1996; (Machabeli moving al. t field while the experience energy along rotational then ward gain will and motion), corotating force plasma and (bead-on-wire centrifugal jet mag- base line the magnetospheric the field at the injected the to by particle, parallel test screened charged component is A line field field electric netic the Fend (e.g., where hole fr black topology a the and with of velocity rigidly model rotational rotate the to efficiency of two-dimensional assumed tion is field idealized magnetic the acceleration where an and consider We energization Particle 2.1. httegnrtdMDjttplg ol loalwe allow also could topology jet MHD generated the that htcnrfglydie uflw rmAN ol ec bu reach indica (e.g. could factors example, AGNs Lorentz formation for from simulations, outflows jet centrifugally-driven MHD Gangadhara Meier that Detailed relativistic 1999; & 2001). Camenzind of al. 1997; Thomas Fendt et context 2005; 1982; Payne the al. & Blandford in et and Contopolous Machabeli 1997; 2007) 1999; Chedia Bogovalov 1994; al. 1996; Rogava Gangadhara et & 1996; Machabeli al. 1969; et Gold (e.g., models .Cnrfglaclrto fts particles test of acceleration Centrifugal 2. the to relevant parameters for depth M87. more ex- galaxy in we radio Here 2000, issue 2007). Mannheim this al. et & amine Osmanov Rieger 2002; sy Xu 1997; hole RM00; Lesch henceforth black & accreting (Gangadhara sub-Eddington par tems in test energies charged high relativistic very of acceleration centrifugal ihseta index spectral with 1 . E teml–gmary:ter niiul M87 individual: – theory rays: gamma – nthermal o yia ih yidrrdu of radius cylinder light typical a for d 7 o o h hlegn H msinfaue nM7 Our M87. in features emission VHE challenging the for ion etosi h rsneo nes opo oss n the and losses, Compton inverse of presence the in lectrons mary rmnnaindatv aatcnce (AGNs), nuclei galactic active non-aligned from amma-rays iual ceeae lcrn a ec oet factors Lorentz reach may electrons accelerated rifugally ditor mds htn oteTVrgm.Usatrn of Upscattering regime. TeV the to photons disk mm Γ 1 , b 3 α ∼ c 0 thsbe elzdframs decade a almost for realized been has It 10. ≃ 1 . .Tecaatrsi aiblt time variability characteristic The 2. [email protected] dt epeeti the in present be to ed r L ≃ 5 r ytm and systems e s . ces.

c S 2018 ESO ilsto ticles up 1997), t ffi cient with ac- ut- he lk s- te e . , 2 F.M. Rieger and F.A. Aharonian: Variable VHE gamma-ray emission from non-blazar AGNs motion can be most conveniently analyzed in the framework of provided that the corotation condition can be satisfied for such Hamiltonian dynamics (RM00) by noting that the Hamiltonian a range of Lorentz factors. The latter qualification seems impor- H for a particle (rest mass m0) on a relativistically rotating wire tant, as it may well happen that the particle motion becomes so 2 2 2 2 is a constantofmotion,and is givenby H = γ m0 c (1 Ω r /c ), perturbed by radiation recoil that the bead-on-the-wire approxi- where γ = 1/(1 Ω2r2/c2 r˙2/c2)1/2 is the Lorentz− factor, mation is no longer a useful concept. For highly underluminous − − < 3 15 Ω = c/rL the angular velocity of the field line and rL the light AGN sources with, e.g. Uph 0.01 erg/cm and rL 5 10 ∼ cylinder radius. The resultant equation of motion cm (see below, 3), for which Eq. (5) would otherwise≃ × imply IC > 9§ 1/2 ∂2r ∂r ∂γ γmax 1.3 10 m˜ , this may in fact be the case. Indeed, from γ + = γ Ω2 , ∼ × 2 r (1) the bead-on-wire requirement that the radiation reaction force ∂t ∂t ∂t that results from inverse Compton scattering, i.e., F P /c rad ≃ IC can be solved analytically yielding (Machabeli & Rogava 1994; (with PIC the single particle Thomson power) should be (much) RM00) smaller than the Lorentz force FL, one finds that achievable elec- r(t) = r cn(λ Ω t) (2) tron Lorentz factors should be (much) smaller than L 0 − for the time-dependence of the radial coordinate r, assuming a 1/2 1/2 RR 8 B(rL) 0.01 particle to be injected at r0 with Lorentz factor γ0, where cn is γmax 7.3 10 , (6) ≃ × 10G ! Uph ! the Jacobian elliptic cosine and λ0 is a Legendre elliptic integral of the first kind, and where B(rL) is the magnetic field strength at the light cylinder 1 radius. Yet, even for cases where radiative losses might be ne- γ(r(t)) = , (3) √m˜ (1 Ω2r2(t)/c2) glected, the breakdown of the bead-on-the-wire approximation − (roughly occurring when the Coriolis force exceeds the Lorentz 2 2 2 2 when expressed in terms of r(t),withm ˜ = 1/(γ0[1 Ω r /c ]) force) will prevent a particle from achieving infinite energies − 0 ≤ 1. If we neglect, for a moment, radiative losses and the break- (RM00). In the simplest case, this restricts achievable particle down of the bead-on-the-wire approximation, a particle would energies to Lorentz factors below 1 reach the light cylinder rL within a time tv = λ0/Ω rL/c, where, as a consequence of the reversal of the centrifugal≤ ac- 2/3 2/3 2/3 BB 8 1/6 B(rL) me rL celeration, it would change direction and move inward again γ 2.0 10 m˜ − , (7) max ≃ × 10G ! m ! 1015cm (Machabeli & Rogava 1994). Knowing the dependence of γ 0   on r, we can easily determine the local acceleration time scale which implies lower Lorentz factors for protons than for elec- (RM00) trons. In any case, once the Lorentz factors become too high, the γ c 1 Ω2r2/c2 inertia of the particles overcomes the tension in the field line, t = = − , (4) so that the line is swept back opposite to the sense of rotation, acc p γ˙ 2 Ω2 r 1 m˜ [1 Ω2r2/c2] slowing down acceleration and introducing curvature radiative − − p asy 1/4 1/2 losses, thus ultimately preventing infinite energy growth. In what which for large γ approaches tacc = 1/(2 Ω m˜ γ ) 16.6m ˜ 1/4 (106/γ)1/2 (r /1015cm) [s]. The (local) acceleration≃ follows, it is assumed that achievable Lorentz factors γ always − L satisfy the relation γ< min γIC , γRR , γBB . time scale thus decreases with increasing γ. Equation (3) implies { max max max} that the particle Lorentz factor increases dramatically as a parti- cle approaches the light cylinder. Therefore, in many cases the 3. Application to the radio galaxy M87 characteristic linear size of our region of interest is much smaller than the light cylinder and the curvature radius of the field, a fact 3.1. Phenomenological background that may qualify the presumed straight field line approach. The nearby (distance 16 Mpc) giant elliptical galaxy M87 ∼ 9 hosts one of the most massive black holes MBH 3 10 M 2.2. Efficiency constraints and maximum particle energy ≃ × ⊙ (e.g., Marconi et al. 1997), with Schwarzschild radius rs = 2GM /c2 8.9 1014 cm, and a prominent one-sided (kpc- In reality, unlimited growth will be prohibited by radiative en- BH ≃ × ergy losses, the breakdown of the bead-on-the-wire approxi- scale) jet visible from radio to X-ray wavelengths (e.g., Marshall mation or the bending of the field line with increasing inertia et al. 2002; see Ly et al. 2007 for a possible radio counter-jet de- (Gangadhara & Lesch 1997; RM00). Inverse Compton upscat- tection). HST observations have revealed superluminal motion of jet components at 0.5 kpc from the central black hole, in- tering of accretion disk photons, for example, leads to an en- ∼ dicative of bulk flow Lorentz factors Γb 6 and a jet orienta- ergy loss of the particle characterized by a time scale tcool 1/γ ∼ ∝ tion of θ 19◦ to the line of sight (Biretta et al. 1999; but see that decreases faster than tacc and so introduces a natural lim- ∼ itation. Indeed, only for highly underluminous AGN sources also Ly et al. 2007 for larger radio θ), suggesting that M87 is a will centrifugal acceleration be sufficiently efficient to accelerate non-blazar jet source, characterized by only moderate Doppler electrons well beyond Lorentz factors of one hundred (RM00; factors. Superluminal radio features have also been detected in Xu 2002). Assuming the inverse Compton scattering process to HST-1 located at around 100 pc (Cheung et al. 2007), although be approximately describable by the Thomson limit for a quasi- no superluminalmotion has been found on small scales (Kovalev 3 et al. 2007). HESS observations have recently shown M87 to be isotropic photon distribution with energy density Uph [erg/cm ], we have t 3 107/[γ U ] s. Balancing acceleration by cool- cool ≃ · ph 1 We are grateful to Osmanov et al. (2007) for making us aware of ing (tacc = tcool) for electrons thus gives a maximum electron an incorrect conclusion in RM00. Although formula (19) in RM00 is Lorentz factor (see also Osmanov et al. 2007) BB correct, an inconsistent set of parameters has been used to estimate γmax 2 for the applications presented. This becomes relevant for highly under- √m˜ 1015cm γIC 3.2 106 . (5) luminous sources, where achievable Lorentz factors can be much higher max 2 r ≃ × Uph L ! than previously concluded in RM00. F.M. Rieger and F.A. Aharonian: Variable VHE gamma-ray emission from non-blazar AGNs 3 a rapidly variable (observed time scale of 2 days) TeV emit- is generally produced in the innermost region of the accretion ting source, yet with a relatively low (isotropic)∼ TeV luminosity flow. For M87, the characteristic synchrotron (peak) frequency of several times 1040 erg/s (Aharonian et al. 2006). The total becomes nuclear (disk and jet) bolometric luminosity of M87 has been 5/4 2 42 1/2 r − Te estimated to be of order Lbol 10 erg/s or less (Reynolds et al. ν (r) 4 1013 m˙ Hz , (8) ∼ s b 9 1996; Owen et al. 2000), indicating that M87 is a highly under- ≃ × rs ! 5 10 K 6 × luminous source with le 3 10− , where le = Lbol/LEdd and s ≤ × and the associated radio luminosity LR νsLν is given by LEdd is the Eddington luminosity. This has led to the proposal ∼ that M87 is a prototype galaxy, where accretion occurs in a two- x 8/5 T 21/5 ν 7/5 L 1042 m˙ 4/5 M e erg/s , (9) temperature, advective-dominated (ADAF) mode characterized R b 3 9 11 ≃ 10  5 10 K 10 Hz by an intrinsically low radiative efficiency (Reynolds et al. 1996; × where (x /103) 1 denotes the dimensionless synchrotron self- Camenzind 1999; Di Matteo et al. 2003). Magnetic flux dragged M ∼ inwards may then build up a rigidly rotating, dipolar magneto- absorption frequency (Yi & Boughn 1998). On the light cylinder 11 sphere, along which disk plasma can be centrifugally accelerated scale (rL 5rs) this implies a peak frequency νs(rL) 2 10 ≃ 40 ≃ × to (bulk) outflow Lorentz factors of Γ (5 10). The generated Hz and a luminosity LR 2 10 erg/s. Thomson upscattering b 2 ∼ × light cylinder scale is likely to be of order≃ r− 5 r (Camenzind (ν γ νs) of these mm soft photons by centrifugally acceler- L s ∼ 7 8 & Krockenberger 1992; Fendt 1997; Camenzind∼ 1999; cf. also ated electrons with γ up to (10 10 ) will thus, in princi- ple, lead to VHE photons with∼ energies− up to (0.1 10) TeV. Fendt & Memola 2001 for higher rL if differential rotation is Comptonization of cyclosynchrotron soft photons∼ adds− a power important). The magnetic field lines in global (quasi force-free) α law tail to the disk spectrum above ν , i.e., Lν Lν (ν/ν )− c MHD wind solutions are radial near the black hole horizon, but s ≃ s s asymptotically collimated into a cylindrical structure beyond the with power index αc = ln τes/ ln A (Mahadevan 1997). For M87 withm ˙ = 1.6 10−3, viscosity parameter α = 0.3 and light cylinder, typically on radial scales 10 rL. This seems con- b − ∼ T = 5 109 K, we obtain× an electron scattering depth τ 0.04 sistent with high frequency VLBI observations in M87, indicat- e × es ≃ > and an amplification factor A 15.7, so that αc 1.2. ing a jet that forms with opening angles 60◦ at the jet base ∼ (3) Suppose that during an≃ active state, test≃ particles are in- (Ly et al. 2007), as expected in MHD models, and a jet radius of > jected with γ0 2 at a constant rate Q and accelerated up to a 50 rs close to the origin (Krichbaum et al. 2006). It has been ∼RR IC ∼proposed recently that efficient pulsar-type particle acceleration threshold γb < γmax < γmax, above which they are considered to may occur in such an environment(Neronov & Aharonian 2007, leave the centrifugal acceleration process due to the breakdown henceforth NA07). of corotation. The differential particle energy distribution n(γ) along a field line would satisfy the simplified transport equation ∂n ∂ γ γ n 3.2. Implications for particle acceleration + n δ(γb γ) = Q δ(γ γ0) , (10) ∂t ∂γ "tacc − tcool # ! − τesc − − Let us consider the implications of these findings for the cen- 1 1/2 trifugal acceleration of particles in M87: with tcool γ− and tacc γ− [cf. Eq. (4)]. Above injection, the steady-state∝ distribution∝ in the acceleration region thus be- (1) Firstly, even for the most limiting case where all of the comes observed bolometric luminosity is assumed to originate within a 1 region rL, so that the energy density of the radiation field may γ − 3/2 2 n(γ) 1 γ− H(γb γ) , (11) be approximated by Uph = Lbol/(4πr c), Eq. (5) would imply IC L ∝ − rγmax ! − that Lorentz factors γIC 107m˜ 1/2 can be achieved, allowing max IC Thomson upscattering of≃ infrared ( 0.01 eV) photons to the i.e., a power law distribution with index 3/2for γb γ . The − ≪ max TeV regime. In reality, this case is certainly≥ over-restrictive, as it emergent (singly scattered, Thomson) inverse Compton spec- assumes that neither the observed jet nor the disk regions beyond trum jIC(ν) from such a hard electron distribution would fol- α rL make a significant contribution to the bolometric luminosity low a power law jIC(ν) ν− with index α = 0.25 for ν 2 ∝ ≪ output, which we consider unlikely. Indeed, if the relevant lu- 4 γbνs. Integrating Eq. (11) over γ gives the number of parti- minosity is an order of magnitude or more smaller, as expected cles along a field line N Q/Ω. Electrons, escaping quasi- ≃ 7 in the ADAF scenario (see estimate L below), the maximum monoenergetically with γb 10 from the acceleration mech- R ∼ Lorentz factors implied by Eq. (5) will be at least two orders of anism and encountering the Comptonized disk photons (ν > νs), magnitude higher. This suggests that electron Lorentz factors up can Thomson upscatter them to the TeV regime, producing a 7 8 2 to γ (10 10 ) may be well possible (cf. Eqs. [5]-[6]). power law-like energy distribution above 4γbνs with index αc ∼ − 1.2, consistent with the value 1.22 0.15 derived for the HESS∼ (2) If accretion in M87 indeed occurs in an ADAF mode, ± the emitted disk spectrum will range from the radio up to the X- 2005 observations of M87 (Aharonian et al. 2006). (4) The number of escaping particles per unit time is n ray regime and beyond: the radio part is produced by synchro- e ∼ cyclotron emission of thermal electrons (T 5 109 K), the n(γb)γb/τesc. Thus, within some time smaller than the cooling e ≃ × time (∆t = ρ t ,ρ< 1), we accumulate N n ∆t particles that optical/UV/X-rays arise via inverse Compton scattering of ra- cool b ∼ e dio soft photons, and the hard X-rays are due to bremsstrahlung can IC upscatter Comptonized disk photons. We can roughly es- and multiple Compton scattering (Mahadevan 1997; Narayan et timate the associated TeV luminosity from LIC Nb PIC, where = 2 ∼ al. 1998; Yi 1999). The ADAF equipartition magnetic field for PIC 1.3σT cγbUph is the single particle Compton power per 4 1/2 5/4 unit volume. This gives M87 is of order B 2.5 10 m˙ (r/r )− G, wherem ˙ is eq ∼ × s the accretion rate in units of the Eddington rate (cf. Yi 1999). γ 2 N 5r 3 40 ρ b s / . For the inferred Bondi accretion ratem ˙ = m˙ b 1.6 10− (Di LIC 10 7 36 erg s (12) 3 ∼ × ∼ 5 10  10  rL ! Matteo et al. 2003) this gives B(rs) 10 G, suggesting possible × (radial) jet magnetic field strengths∼ close to the light cylinder of To achieve a Compton luminosity comparable to the observed B(r ) (10 50) G. The highest radio emission in an ADAF TeV luminosity of L 3 1040 erg/s (Aharonian et al. 2006), L ∼ − TeV ≃ × 4 F.M. Rieger and F.A. Aharonian: Variable VHE gamma-ray emission from non-blazar AGNs we thus need N 3 1036/ρ particles along field lines. Denoting strong relativistic beaming effects, while for luminous quasars the relevant acceleration∼ × volume by ∆V η r2 ∆r, with charac- internal absorption of gamma-rays becomes dominant. Based ∼ L teristic length scale ∆r = γ/(∂γ/∂r) (γ0/γb) rL and η < 1, on a simple toy model we have shown that efficient centrifu- | | ∼ 2 2 the corresponding kinetic energy density n(γb)γbme c /∆V (for gal acceleration of electrons in the vicinity of the light cylinder > 5 2 ρ η 10− ) is still well below the energy density B /(8π) of could provide a natural explanation for variable (time scale of ∼ the magnetic field, suggesting that the presumed (quasi force- one day) VHE emission with a hard inverse Compton spectrum free) MHD field structure is still a valid approximation (cf. also as observed in M87. Our models fits well with other evidence Osmanov et al. 2007). for advection-dominated accretion in M87 and may indeed be (5) In principle, TeV gamma-rays can be strongly attenuated regarded as providing some further corroboration for the pres- due to photon-photon pair production in the background disk ence of such modes in highly underluminous AGNs. photon field. The narrow dependence of the cross-section σγγ As always, there are a number of subtleties whose impact on the on the product of photon energies implies that VHE photons of presented results need to be explored in more details including energy E interact most efficiently with infrared background pho- general relativistic effects, anisotropic scattering modifications, tons of energy ǫIR (1TeV/E) eV. The optical depth τ for a γ- quasi rigid rotation and plasma instabilities. The extent to which ≃ ff ray photon in a background field of infrared luminosity LIR and our conclusions might be a ected may require fully relativistic size RIR thus becomes (cf. NA07) modelling. Yet, given the demonstrated potential of centrifugal acceleration and our current understanding of relativistic jet for- LIRσγγ L r E τ , = . IR l , mation, this may represent a program worth pursuing. (E RIR) 0 2 41 (13) 4πRIRǫIR ≃ 10 erg/s! RIR ! 1TeV Acknowledgements. Discussion with and comments by John Kirk, Christian indicating that due to its low bolometric luminosity M87 could Fendt and Karl Mannheim are gratefully acknowledged. be well transparent to VHE gamma-rays, even if almost all of 41 the observed infrared luminosity LIR 10 erg/s (Whysong & References Antonucci 2004) is (somewhat unrealistically)≃ taken to be pro- Aharonian, F., et al. (HESS collaboration) 2006, Science 314, 1424 duced on a scale R r 60 r . Note that even if τ would IR ∼ l ∼ s Aharonian, F., et al. (HESS collaboration) 2007, ApJL 664, L71 become larger than one, γ-rays from the last transparent layer Albert, J., et al. (MAGIC collaboration) 2007, ApJ 669, 862 are still able to escape, so that the VHE flux would not simply Biretta, J.A., Sparks, W.B., Macchetto, F. 1999, ApJ 520, 621 decrease exponentially by exp( τ), but only by a factor of τ Blandford, R.D., Payne, D.G. 1982, MNRAS 199,883 − ∼ Bogovalov, S.V. 1997, A&A 327, 662 (NA07). Camenzind, M., Krockenberger, M. 1992, A&A 255, 59 (6) The number of electrons escaping quasi- Camenzind, M. 1999, in: The radio galaxy Messier 87, eds. H.-J. R¨oser & K. monoenergetically from the acceleration mechanism in the Meisenheimer, LNP 530, p. 252 vicinity of the light cylinder is of order Nb. Once these Chedia, O.V., Kahniashvili, T.A., et al. 1996, Ap&SS 239, 57 energetic particles encounter non-vanishing perpendicu- Cheung, C.C., Harris, D.E., Stawarz, L. 2007, ApJL 663, L65 Contopoulos, I., Kazanas, D., Fendt, C. 1999, ApJ 511, 351 lar and/or the turbulent plasma magnetic fields, they can Di Matteo, T. et al. 2003, ApJ 582, 133 1/3 produce synchrotron emission arising as Lν ν be- Fendt, C. 1997, A&A 319, 1025 low, and decaying exponentially above the peak∝ frequency Fendt, C., Memola, E. 2001, A&A 365, 631 7 2 Gaidos, J.A., et al. (Whipple collaboration) 1996, Nature 383, 319 νsyn 50(γb/5 10 ) (B sin α/1G) MeV with a total lumi- ∼ × Gangadhara, R.T. 1996, A&A 314, 853 nosity of order L P N 0.06 L (B sin α)2, where P syn ∼ syn b ∼ IC syn Gangadhara, R.T., Lesch, H. 1997, A&A 323, L45 is the single particle synchrotron power. In order to satisfy the Georganopoulos, M., Perlman, E.S., Kazanas, D. 2003, ApJL 634, L33 restrictions imposed by the existing (yet non-contemporaneous) Gold, T. 1969, Nature 221, 25 upper limit on the M87 flux in the EGRET energy band above Krichbaum, T.P. et al. 2006, J. Phys. Conf. Ser. 54, 328 Kovalev Y.Y., et al. 2007, ApJL 668, L27 100 MeV (e.g., Reimer et al. 2003), the effectively encountered Ly, C., Walker, R.C., Junor, W. 2007, ApJ 660, 200 fields should be smaller than 1 Gauss. This seems consistent Machabeli, G.Z., Rogava, A.D. 1994, Phys. Rev. A 50, 98 with independent estimates suggesting∼ a strength of the random Machabeli, G.Z., Nanobashvili, I.S., Rogava, A.D. 1996, Radiophysics and field component close to the black hole of below one Gauss Quantum Electronics 39, 26 (NA07). The overall spectral energy distribution in M87 is then Machabeli, G.Z., Osmanov, Z.N., Mahajan, S.M. 2005, PhPL 12, 2901 Mahadevan, R. 1997, ApJ 477, 585 likely to consist of a number of different contributions, involving Marconi, A. et al. 1997, MNRAS 289, L21 also other leptonic (Georganopoulos et al. 2003; NA07) and Marshall, H.L., et al. 2002, MNRAS 564, 683 perhaps even hadronic (Reimer et al. 2004) processes. If so, then Mastichiadis, A., Kirk, J.G. 2002, PASA 19, 138 no straighforward X-ray–TeV correlation might be expected. Meier, D.L., Koide, S., Uchida, Y. 2001, Science 291, 84 Narayan, R., Mahadevan R., Quataert, E. 1998, in: Theory of Black Hole (7) As shown above, accelerating particles up to the light Accretion Disks, eds. M.A. Abramowicz et al., Cambridge, p. 148 cylinder typically takes a time rL/c, suggesting a characteristic Neronov, A., Aharonian, F.A. 2007, ApJ 671, 85 (NA07) rL 5 rs Osmanov, Z., Rogava, A., Bodo, G. 2007, A&A 470, 395 variability time scale for M87 of tv c c 1.7 days, well consistent with the observed TeV time≃ scale∼ of≃∆t 2 days, a Owen, F.N., Eilek, J.A., Kassim, N.E. 2000, ApJ 543, 611 ∼ Reimer, A., Protheroe, R.J., Donea, A.-C. 2004, A&A 419, 89 fact that may further validate the assumptions of the presented Reimer, O., Pohl, M., Sreekumar, P., Mattox, J.R. 2003, ApJ 588, 155 model. Reynolds, C.S. et al. 1996, MNRAS 283, L111 Rieger, F.M. 2004, ApJL 615, L5 Rieger, F.M., Mannheim, K. 2000, A&A 353, 473 (RM00) 4. Conclusions Salvati, M., Spada, M., Pacini, F. 1998, ApJL 495, L19 Thomas, R.M.C., Gangadhara, R.T. 2007, A&A 467, 911 VHE radiation from low-luminous, non-blazar AGN jet sources Whysong, D., Antonucci, R. 2004, ApJ 602, 116 like M87 could providean ideal test laboratoryforthe analysis of Xu, Y.D. 2002, A&A 381, 357 particle acceleration processes close to the supermassive black Yi, I. 1999, in: Astrophysical Disks, ASP Conf. Ser. 160, eds. J.A. Sellwood, J. Goodman, p. 279 hole event horizon. In blazars with their relativistic jets point- Yi, I., Boughn, S.P. 1998, ApJ 499, 198 ing towards us, most of these traces are likely to be masked by