EDP Sciences
Free Access
Issue
A&A
Volume 530, June 2011
Article Number A75
Number of page(s) 5
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201014725
Published online 12 May 2011

© ESO, 2011

1. Introduction

During the last decades the interest in very high energy (VHE) radiation of active galactic nuclei (AGNs) has considerably increased because of the detailed data from the new telescopes (MAGIC, HESS, EGRET, Fermi). Active galactic nuclei, and especially blazars reveal a strong connection between high-energy γ-ray and radio emission. Since the EGRET era, a direct connection between radio and γ-ray flux density has been intensively discussed. In particular, Bloom (2008) performed a statistical analysis of the broadband properties of EGRET blazars and showed a strong positive correlation between the total radio luminosity at 8.4 GHz and the γ-ray luminosity in the EGRET band (). The comparison of the radio flux density at 8.4 GHz and the γ-ray (>100 MeV) photon flux measured by the Fermi γ-ray Space Telescope has been presented by Giroletti et al. (2010) and the Spearman’s rank correlation coefficient was estimated to be 0.57. Kovalev et al. (2009) found a positive correlation between the pc-scale radio flux at 15 GHz and the γ-ray photon flux in a sub-band of the Fermi/LAT detector (100 MeV–1 GeV), and a time separation of a few months between γ-ray and radio flares.

Most commonly the VHE emission of AGNs is explained in terms of the inverse Compton upscattering (Blandford et al. 1990) and curvature radiation (Gangadhara 1996; Thomas & Gangadhara 2005). The major reason for this is that if the magnetic field is strong, the synchrotron cooling timescale is very short and relativistic particles very soon transit to their ground Landau level, which prevents the subsequent emission process (Osmanov & Machabeli 2010).

Recently, it has been shown that the quasi-linear diffusion (QLD) may create necessary conditions in the magnetospheres of AGNs and pulsars for avoiding the synchrotron emission damping despite the strong magnetic field. According to the model of Kazbegi et al. (1992), the cyclotron instability appears in the magnetospheres of pulsars. These unstable modes lead to the feedback of the excited waves on particles by means of the process of the diffusion. The pitch angles are no longer negligible and the emission process continues (e.g. Machabeli & Usov 1979; Lominadze et al. 1979; Malov & Machabeli 2001; Chkheidze & Machabeli 2007; Chkheidze et al. 2011). This model was applied to the Crab pulsar (Machabeli & Osmanov 2009, 2010; Chkheidze et al. 2011) to explain the recent results from the MAGIC Cherenkov telescope between 2007 October and 2008 February (Albert et al. 2002).

In the context of works presented by Bloom (2008) and Giroletti et al. (2010) (where the observationally evident strong correlation of high-energy γ-rays and low-frequency radiation is discussed), it is of fundamental importance that according to the QLD, two different radiation domains are excited. As was studied by Osmanov & Machabeli (2010) and Osmanov (2010), the cyclotron resonance efficiently generates transverse momenta, which leads to non-zero pitch angles and the subsequent synchrotron radiation. This means that a connection between two radiation domains is naturally introduced into the quasi-linear diffusion process. In particular, examining the QLD in AGNs, Osmanov & Machabeli (2010) studied the generation of X-rays associated with the radio emission, and the work of Osmanov (2010) was related to high-energy emission associated with submillimeter/infrared radiation.

In this paper we study the production of high-energy radiation in the GeV-TeV interval in the magnetospheres of supermassive black holes. The paper is organized as follows. In Sect. 2 we present the model, in Sect. 3 we apply the mechanism to AGNs and in Sect. 4 we summarize our results.

2. Theory

In general, Lorentz factors of magnetospheric plasma particles in AGNs1 lie in a broad interval ranging from  ~1 up to  ~108 (Osmanov et al. 2007; Rieger & Aharonian 2008), therefore, for simplicity one can consider the magnetosphere to be composed of the electron-positron plasma component with relatively low Lorentz factors, γp, and highly relativistic electrons, the so-called beam component electrons with the Lorentz factor γb (γb ≫ γp).

In this section we consider the QLD mechanism and derive an average value of the pitch angles driven by the diffusion. In general the process of energy dissipation that takes place in strong magnetic fields is dynamically controlled by two forces: one is the synchrotron radiative force (Landau & Lifshitz 1971) (1)where , ψ is the pitch angle, ωB ≡ eB/mc is the cyclotron frequency, B is the magnetic induction, c is the speed of light, e and m are electron’s charge and the rest mass, respectively; the other is the force G, responsible for the conservation of the adiabatic invariant, (Landau & Lifshitz 1971) (2)where ρ is the curvature radius of the magnetic field lines. In general, in the magnetospheres of AGNs two effects take place: the energy dissipation provoked by the forces, F,G, which lead to the decrease of the pitch angles, the QLD, which opposes the first effect. The dynamical process saturates when the effects of the above-mentioned forces are balanced by the diffusion, which leads to non-zero pitch angles and the subsequent radiation.

As discussed by Machabeli & Usov (1979), Machabeli & Osmanov (2010), Osmanov & Machabeli (2010), for an activation of the QLD, the cyclotron modes must be excited, which leads to the diffusion in terms of feedback of the cyclotron waves on particles. On the other hand, according to the work of Kazbegi et al. (1992), the anomalous Doppler effect can efficiently induce unstable cyclotron waves with the corresponding frequency (Machabeli & Usov 1979; Osmanov & Machabeli 2010) (3)where , is the plasma frequency and np is the plasma density. Here we assume that energy in plasma is uniformly distributed, therefore, the plasma density can be approximated to be nbγb/γp, where nb is the beam density (Machabeli & Osmanov 2009). Although the aforementioned expression is a direct consequence of a resonance character of the cyclotron instability, in exciting waves all resonance particles (with broad energy spectra) participate, therefore the range of spectral frequencies is wide and the waves are not characterized by a monochromatic signature. In the framework of the present approach, one has two simultaneously generated connected VHE γ-rays and low-frequency cyclotron modes. This result seems to be promising in the context of recent investigations, because it shows the strong connection between the γ-ray and radio fluxes (Bloom 2008; Giroletti et al. 2010).

A similar problem was studied by Osmanov & Machabeli (2010) who examined the physical regime |G ⊥ | ≫  |F ⊥ |  and  |G ∥ |  ≪ |F ∥ | and studied the production of X-rays (0.13–100 keV) connected with the radio emission (0.04–35 MHz). A different regime, |G ⊥ | ≪ |F ⊥ | and |G ∥ | ≪ |F ∥ |, was examined by Osmanov (2010). This regime shows the possibility of production of MeV–GeV γ-rays, which are strongly correlated with submillimeter/infrared radiation domains. In the present paper we examine the same regime, but for different parameters. Unlike the previous work, where the QLD was studied in the light cylinder area where the magnetic field is of the order of  ~ 300 G, in the present paper the physics of the QLD is examined for a region close to the supermassive black hole event horizon where B ~ 104 G (see Thorne et al. 1988). Another difference is the fact that here we study the cyclotron waves generated by plasma particles with Lorentz factors, γp ~ 2 − 4, whereas values of γp considered by Osmanov (2010) were higher by two orders of magnitude. As we will see, the aforementioned differences significantly change the results. By comparing the transverse and longitudinal components of the forces, one obtains Evidently, the relations |G ⊥ | ≪ |F ⊥ | and |G ∥ | ≪ |F ∥ | are satisfied if ψ ≫ 2.2 × 10-6 rad if we consider the physical parameters γb ~ 108, B ~ 104 G, Rg ~ ρ.

By assuming a quasi-stationary scenario (/∂t = 0), the corresponding kinetic equation will be different from that of Machabeli & Osmanov (2010) (6)where f = f(ψ,γb) is the distribution function of particles, (7)are the diffusion coefficients and |Ek|2 is the energy density per unit wavelength.

For estimating |Ek|2, we assume that half of the plasma energy density, mc2nbγb/2, converts to the energy density of the waves |Ek|2k, then an expression of |Ek|2 writes as follows (8)For solving the aforementioned kinetic equation with respect to the pitch angles we separate the variables by expressing the distribution function as χ(ψ)f(γb). Then the equation governing the distribution by pitch angles is given by (9)and the corresponding solution writes as (10)where (11)For simplicity we use the average value of the pitch angle, which after taking Eqs. (10, 11) into account leads to (12)Clearly, despite very efficient synchrotron losses the quasi-linear diffusion works against the dissipative forces and maintains non-zero pitch angles, which prevent the synchrotron mechanism from damping. Under these conditions, the relativistic particles with Lorentz factors, γ, radiate in the synchrotron regime and emit photons with energies (Rybicki & Lightman 1979) (13)

3. Results

According to the observational evidence, AGNs are characterized by emission in the VHE domain, which cannot be explained in the framework of thermal processes. This can also be explained through the existence of ultra-relativistic particles in the magnetosphere. In general, there are several acceleration mechanisms that might provide very high Lorentz factors. We show that processes such as Fermi-type acceleration (Catanese & Weekes 1999) and acceleration owing the the black hole dynamo mechanism (Levinson 2000) can provide Lorentz factors of the order of  ~107−9 close to the supermassive black holes, i.e. at r ~ Rg, where Rg ≡ 2GMBH/c2 ≈ 3M9 × 1014 cm is the gravitational radius of the black hole, M9 ≡ M/(109   M) and M is the solar mass.

Osmanov & Machabeli (2010) have shown that the synchrotron cooling timescale for ultra-relativistic electrons in the magnetospheres of AGNs close to the black hole is of the order of  ~10-5–10-4 s. To estimate the efficiency of the energy losses it is reasonable to calculate the kinematic timescale tkin ~ Rg/c. Considering the supermassive black hole with M9 = 1, one can see that tkin ~ 104 s, which by many orders of magnitude exceeds the dissipation timescale, therefore, the cooling process is extremely efficient and the electrons stop emitting very soon.

Anomalous Doppler effect significantly changes the overall radiative pattern of the system. Equation (3) clearly shows that the mentioned effect efficiently induces unstable cyclotron waves with the frequency (Osmanov & Machabeli 2010) (14)which in turn lead to the feedback of these modes to the particles and switch on the QLD, which generates the non-vanishing pitch angles and the subsequent synchrotron process.

thumbnail Fig. 1

Behaviour of ϵTeV with respect to γb. The set of parameters is M9 = 1, B = 104 G, nb = 20 cm-3 and γp ∈  {2;3;4} .

Open with DEXTER

As a first example, we examine a supermassive black hole with M9 = 1 and the typical value of the magnetic induction for a nearby zone of the black hole B = 104 G. Since we assume that energy in plasma is uniformly distributed, the value of plasma density nbγb/γp, for nb = 20 cm-3 and B = 104 G is of the order of 109 cm-3, which is compatible with observations of narrow and broad line regions and indicates that the plasma density in the environment of AGNs lies in the interval  ~ (105 − 1014) cm-3 (Xu et al. 2007; Arav et al. 2002). We consider the beam component’s Lorentz factors to be of the order of 108 and the plasma component with the Lorentz factor, γp = 2. From Eq. (14) one can show that the cyclotron frequency of the order of 87 MHz is excited, which in turn, by means of the QLD, leads to the average value of the pitch angle  ~ 0.2 rad (see Eqs. (11), (12)) and the subsequent synchrotron emission in the VHE domain  ~ 240 GeV (see Eq. (13)).

In Fig. 1 we show the dependence of on γb. The set of parameters is M9 = 1, B = 104 G, nb = 20 cm-3 and γp ∈  { 2;3;4 } . Clearly, physically realistic parameters might guarantee the synchrotron emission from  ~ 70 GeV (γb = 108, γp = 4) to 0.9 TeV (γb = 2 × 108, γp = 2). These plots show a continuously increasing character of with respect to γb, which is a natural result of the fact that more energetic beam particles radiate photons with higher energies. Unlike the beam Lorentz factors, the behaviour of with respect to γp is a continuously decreasing function. This can be seen from Eqs. (11), (12): , which by combining with Eqs. (7), (8) reduces to . Therefore, by considering higher values of γp, the resulting synchrotron emission energy will be lower.

In Fig. 2 we show the behaviour of ϵTeV with respect to ν and γp. The set of parameters is M9 = 1 and B = 104 G, nb = 20 cm-3. The 2D surface shows that for lower values of the plasma Lorentz factors, the synchrotron emission energy is higher and the corresponding cyclotron frequency relatively is lower. Figure 3 shows the dependance of on the cyclotron radio frequency for three different values of plasma Lorentz factors γp =  { 2;3;4 } . Evidently, the VHE γ-rays in the TeV domain are connected with the radio emission at  ~ 87 MHz, and the higher frequency radio emission  ~ 1.4 GHz is connected with the lower energy (~70 GeV) γ-ray synchrotron emission.

thumbnail Fig. 2

Behaviour of ϵTeV versus ν and γp. The set of parameters is M9 = 1 and B = 104 G, nb = 20 cm-3.

Open with DEXTER

thumbnail Fig. 3

Behaviour of ϵTeV versus ν. The set of parameters is M9 = 1 and B = 104 G, nb = 20 cm-3 and γp =  { 2;3;4 } . Evidentlythe curves are truncated, which is a direct result of the fact that the curves are one-dimensional sections of a two-dimensional surface.

Open with DEXTER

One of the major magnetospheric parameters is the density of the beam component. In Fig. 4 we show the behaviour of versus nb. The set of parameters is M9 = 1 and B = 104 G, γp = 2 and nb =  { 0.2;2;20 }  cm-3. Obviously, the higher the beam density, the higher the γ-ray synchrotron radiation energy. As we already mentioned, the average value of the pitch angle is proportional to , which by combining with (see Eqs. (7), (8)) confirms the results shown in the plots. For the considered interval of nb the photon energy varies in the range, 8.5 GeV–0.9 TeV.

thumbnail Fig. 4

Behaviour of ϵTeV versus γb. The set of parameters is M9 = 1 and B = 104 G, γp = 2 and nb =  { 0.2;2;20 }  cm-3.

Open with DEXTER

Figure 5 presents the dependance of on ν. From the plots we see that the TeV emission is connected to the radio band  ~22 MHz and the high-energy emission, 8.5 GeV is connected to  ~8.7 GHz radio domain. The present results are different from those of Osmanov & Machabeli (2010) and Osmanov (2010), where a different parametric space, has been considered. In particular, Osmanov & Machabeli (2010) examined M9 = 1, γp = 2, γb ~ 106 − 9, nb = 2000 cm-3 and Osmanov (2010) considered M9 = 1, γp ∈  { 200;250;300 } , γb ~ 108 and nb =  { 5;10;15 }  cm-3.

According to the quasi-linear approach, the efficiency of the QLD depends on the efficiency of the cyclotron instability. For this purpose it is essential to consider the growth rate of the instability. Kazbegi et al. (1992) showed that the increment, Γ, of the modes is given by (15)and (16)where is the plasma frequency of resonant (beam) electrons. For the parameters nb = 20 cm-3, γp = 2 (see Fig. 1), by assuming ux ~ c and ρ ~ Rg, one can show that , and the growth rate is defined by Eq. (15). Therefore, the corresponding cyclotron timescale, tcyc = 1/Γ, is of the order of 10-3 s, which is less by many orders of magnitude than the kinematic timescale, tkin ~ Rg/c ~ 104 s. The same relation between the timescales is valid for the rest of the considered cases; therefore, the studied mechanism is extremely efficient.

thumbnail Fig. 5

Behaviour of ϵTeV versus ν. The set of parameters is M9 = 1 and B = 104 G, γp = 2 and nb =  { 0.2;2;20 }  cm-3. The curves are truncated for the same reason as in Fig. 3.

Open with DEXTER

The interest in the QLD is twofold: on the one hand, the quasi-linear diffusion is an efficient mechanism that drives the synchrotron process despite very short synchrotron cooling timescales. On the other hand, as we have seen, the QLD provides the possibility of generating connected VHE γ-rays and radio emission on mpc scales. In particular, Bloom (2008) has found a high correlation of γ-ray luminosity and radio luminosity. The work of Giroletti et. al (2010) is particularly interesting, because they show that radio emission at 8.4 GHz is connected to the VHE radiation (> 100 MeV). As we show in Fig. 5, the radio emission of the order of 8.4 GHz is associated with γ-rays in the GeV band, therefore, the study of the QLD is very promising.

4. Summary

The main aspects of the present work can be summarized as follows:

  • 1.

    We studied the role of the quasi-linear diffusion with magnetospheric plasma particles in AGNs. For this purpose we considered the kinetic equation that describes the balance between the synchrotron reaction force and the diffusion.

  • 2.

    We examined our model for different values of physical parameters. The QLD was studied versus the Lorentz factors of plasma and beam components respectively and versus the density of the beam component. We found that the higher the beam Lorentz factor, the higher the synchrotron emission energy and moreover, that the higher the plasma Lorentz factor, the lower the radio frequency. We also showed that lower beam densities lead to lower energies in the γ-ray synchrotron, and to higher energies in the radio domain.

  • 3.

    We also found that under favourable conditions the cyclotron instability generates radio emission from  ~22 MHz to  ~9 GHz and the QLD provides synchrotron energies in the interval 900 down to 9 GeV. These results show the existence of connected radio and γ-ray emissions on mpc scales.

The QLD is a mechanism able to explain a connection between radio emission and VHE γ-rays. On the other hand, a complete study requires an investigation of the spectral pattern provided by the mentioned mechanism. According to the standard theory of the synchrotron emission (Bekefi & Barrett 1977; Ginzburg 1981), it is assumed that along the line of sight the magnetic field is chaotic; correspondingly, the pitch angles vary in a broad interval (from 0 to π) and the particles have uniformly distributed

pitch angles. Contrary to this scenario, in the framework of the present model the particles are nontrivially distributed (see Eq. (10)), therefore, the spectral picture will differ from those of Bekefi & Barrett (1977), Ginzburg (1981). Another problem that has to be addressed is the variability of the corresponding VHE radiation. As the investigation shows, for the QLD to be efficient, two major requirements have to be fulfilled: the particles must be highly relativistic and the timescale of the cyclotron instability should be short enough. Therefore, the study of the variability implies the corresponding study of the acceleration timescale and the aforementioned instability timescale. We will investigate this problem in future studies.


1

The magnetosphere of AGNs with their typical lengthscales  cm is a region around the supermassive black hole, where the magnetic field (G – close to the black hole and 10–100 G – in the light cylinder zone) is dominant.

Acknowledgments

The research was supported by the Georgian National Science Foundation grant GNSF/ST07/4-193. I also thank an anonymous referee for helpful suggestions.

References

All Figures

thumbnail Fig. 1

Behaviour of ϵTeV with respect to γb. The set of parameters is M9 = 1, B = 104 G, nb = 20 cm-3 and γp ∈  {2;3;4} .

Open with DEXTER
In the text
thumbnail Fig. 2

Behaviour of ϵTeV versus ν and γp. The set of parameters is M9 = 1 and B = 104 G, nb = 20 cm-3.

Open with DEXTER
In the text
thumbnail Fig. 3

Behaviour of ϵTeV versus ν. The set of parameters is M9 = 1 and B = 104 G, nb = 20 cm-3 and γp =  { 2;3;4 } . Evidentlythe curves are truncated, which is a direct result of the fact that the curves are one-dimensional sections of a two-dimensional surface.

Open with DEXTER
In the text
thumbnail Fig. 4

Behaviour of ϵTeV versus γb. The set of parameters is M9 = 1 and B = 104 G, γp = 2 and nb =  { 0.2;2;20 }  cm-3.

Open with DEXTER
In the text
thumbnail Fig. 5

Behaviour of ϵTeV versus ν. The set of parameters is M9 = 1 and B = 104 G, γp = 2 and nb =  { 0.2;2;20 }  cm-3. The curves are truncated for the same reason as in Fig. 3.

Open with DEXTER
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.