This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

11institutetext: Institute for Fusion Studies, The University of Texas at Austin, Austin, TX 78712, USA; 11email: [email protected] 22institutetext: School of Physics, Free University of Tbilisi, 0183, Tbilisi, Georgia; 22email: [email protected] 33institutetext: E. Kharadze Georgian National Astrophysical Observatory, Abastumani, 0301, Georgia

Resonant energization of particles by radio AGN

S. M. Mahajan 11    Z. N. Osmanov 2233
(Received ; accepted )
Abstract

Aims. A new mechanism of particle acceleration, based on the resonant interaction of a classical electromagnetic wave (EM) with a quantum wave (associated with a relativistic particle), is explored.

Methods. In a model illustrative calculation, we study the fate of a Klein Gordon wave subjected to the intense radio frequency waves generated in the vicinity of an active galactic nuclei (AGN). In the framework of the paper we examine a quantum wave associated with a relativistic particle, and it is shown that the group velocity of the wave approaches the speed of light, implying that the particles resonantly exchange energy with EM waves, eventually leading to acceleration of particles to very high energies.

Results. For typical parameters of under accreting Eddington radio AGN, it is shown that the resonant energization could catapult particles to extreme energies 101620\sim 10^{16-20}eV.

Key Words.:
acceleration of particles - galaxies: active - galaxies: jets - plasmas

1 Introduction

One of the fundamental problems of modern astrophysics is to figure the origin of particle acceleration to extremely high energies. It is very helpful to note that observations in the last decade have established/confirmed a strong correlation of very high energy (VHE) protons with active galactic nuclei (Kim & Kim, 2013). In this paper, we advance a possibly acceleration mechanism operating on particles comprising the magnetospheres in nearby regions of radio AGNs. The mechanism, based on resonant energy transfer from the intense electromagnetic (EM) waves to relativistic quantum particle waves (modeled here by as Klein-Gordon Waves), is an interesting manifestation of the strong affinity between two waves with a similar mathematical structure; the dominant behavior of both the classical electromagnetic wave (the radio waves, for example) and the relativistic Klein-Gordon wave is controlled by essentially the same hyperbolic wave operator.

The literature is full of plausible acceleration mechanisms- the so-called Fermi mechanism (Fermi, 1949) and its several variants (Bell, 1978; Catanese & Weeks, 1999) might account for the origin of extremely high energy cosmic rays; this process, though, has limitations and is efficient when the particles are initially accelerated (Rieger & Mannheim, 2000). One of the papers relying on extreme conditions in the vicinity of compact objects (Blandford & Znajek, 1977) proposes that the poloidal magnetic field close to a black hole induces electrostatic potential, which might efficiently accelerate charged particles; the process of acceleration, however, is limited by a significant screening effect. The so-called turbulent acceleration (Lemoine, 2021) might account for very high energies (Pezzi et al., 2022). In the framework of this mechanism the turbulent plasma is described as a combination of Alfvén waves and magnetosonic modes, where the particle acceleration comes from the phase locking of trajectories of particles with the mentioned waves. As a result, this mechanism might provide energies in the PeV domain (Pezzi et al., 2022). Magneto-centrifugal mechanism of energy pumping has been considered for AGN magnetospheres in a series of papers (Rieger & Mannheim, 2000; Osmanov et al., 2007; Osmanov, 2010) to explain energization of particles up to TeV energies. Taking into account generation of centrifugally driven electrostatic field, similar mechanisms have been invoked for similar class of objects (Osmanov et al., 2014) - the central black hole of Milky Way, for instance (Osmanov et al., 2017).

Harnessing of the strong wave-wave resonant interaction (Mahajan & Asenjo, 2022) for energy transfer is what distinguishes this effort from the aforementioned schemes that also operate in extreme conditions in the vicinity of highly compact objects. In particular, if the phase velocity of the EM wave coincides with the phase velocity of a quantum wave associated with a particle, the latter will be energized to extremely high energies. The resonant enhancement is particularly pronounced when the plasma frequency is much less than the frequency of the EM radiation(Sec.2).

It is, perhaps, important to emphasize that the the manifest wave- wave interaction of this paper is possible only because of the quantum nature of the particle.

The paper is organised as follows: in Sec. 2 we outline a theory of our mechanism, consider radio-load AGN applying the theory, discuss obtained results and in Sec. 3 we summarise them.

2 Basic Theory of Resonant energization

In this section we will briefly outline the mechanism of resonant energization developed in (Mahajan & Asenjo, 2022). Enough detail is given to make this paper self contained.

The group velocity of the quantum wave associated with a relativistic particle of energy (momentum) =E (P),

υg=EP=P(P2+m2)1/2,\upsilon_{g}=\frac{\partial E}{\partial P}=\frac{P}{\left(P^{2}+m^{2}\right)^{1/2}}, (1)

approaches the speed of light (c=1c=1) when PmP\gg m, where P=γmυP=\gamma m\upsilon is the momentum of the particle, γ\gamma its relativistic factor and υ\upsilon - the velocity. Thus for extremely high values of the momentum, the particle could, resonantly, exchange energy with an EM wave that propagates in unison. The resonant energization phenomena is well illustrated in a model calculation (Mahajan & Asenjo, 2022) in which a Klein-Gordon (KG) wave is subjected to a circularly polarized EM wave; the latter is described by the EM four potentials AμA^{\mu} (A0=Az=0,A^{0}=A^{z}=0,),

Ax=Acos(ωtkz),Ay=Acos(ωtkz),A^{x}=A\cos\left(\omega t-kz\right),\;A^{y}=-A\cos\left(\omega t-kz\right), (2)

where ω\omega and kk are, respectively, the frequency and the wavenumber of the EM wave propagating along the zz axis. Because, the directions perpendicular to z are ignorable, the KG/EM system obeys the Klein-Gordon equation (see Mahajan & Asenjo (2016) for details),

(t2z2+2qAKcos(ωtkz)+(K2+m2+q2A2))Ψ=0,\left(\partial^{2}_{t}-\partial_{z}^{2}+2qAK_{\perp}\cos\left(\omega t-kz\right)+\left(K_{\perp}^{2}+m^{2}+q^{2}A^{2}\right)\right)\Psi=0, (3)

where qq is the particle charge, and KK_{\perp} is the wave-vector’s perpendicular component, which on the other hand is the measure of the (conserved) perpendicular momentum and labels (suppressed) the wave function. The explicit tt and zz dependence in Eq. (3) implies that the energy E and PzP_{z} will be functions of time. The resonant solutions of Eq. (3) will emerge when we demand Ψ=Ψ(ωtkz)Ψ(ξ)\Psi=\Psi(\omega t-kz)\equiv\Psi(\xi) (ξωtkz\xi\equiv\omega t-kz); the result is a simple Mathew equation

(ω2k2)d2ψdξ2+(μ+λcosξ)=0,\left(\omega^{2}-k^{2}\right)\frac{d^{2}\psi}{d\xi^{2}}+\left(\mu+\lambda\cos\xi\right)=0, (4)

with μ=K2+m2+q2A2\mu=K_{\perp}^{2}+m^{2}+q^{2}A^{2} and λ=2qAK\lambda=2qAK_{\perp}. Notice that the mathematical operator t2z2\partial^{2}_{t}-\partial^{2}_{z} in Eq. (3) translates to ω2k2\omega^{2}-k^{2} (in Eq. (4) when we impose a solution whose phase factor (ξωtkz\xi\equiv\omega t-kz)is exactly that of the EM wave. In some sense this is the obvious mathematics underlying wave resonance. In particular, notice that the signature of resonance is already explicit in Eq. (4); the equation is singular since ω2k2\omega^{2}-k^{2} tends to zero for EM waves traveling in a tenuous medium. The implication, of course, is that nontrivial solutions must demand derivative of Ψ\Psi to become commensurately large. Consequently, the energy and z momentum of the KG state must also become large since both are proportional to the ξ\xi derivative (as we will soon see).

The leading order WKB solution of Eq. (4) is given by (Mahajan & Asenjo, 2022)

Ψ=Ψ0exp(iSξ+iαsinξ),\Psi=\Psi_{0}\exp(-iS\xi+i\alpha\sin\xi), (5)

where

S=K2+m2+q2A2ω2k2,α=qAKμ(ω2k2)S=\sqrt{\frac{K_{\perp}^{2}+m^{2}+q^{2}A^{2}}{\omega^{2}-k^{2}}},\quad\alpha=\frac{qAK_{\perp}}{\sqrt{\mu\left(\omega^{2}-k^{2}\right)}} (6)

Following simple rules of quantum mechanics, the expectation value of the energy is

E=iΨΨtΨΨ=12lll(S+αcosξ)𝑑z=\langle E\rangle=\frac{i\langle\Psi^{\star}\frac{\partial\Psi}{\partial t}\rangle}{\langle\Psi^{\star}\Psi\rangle}=\frac{1}{2l}\int_{-l}^{l}\left(S+\alpha\cos\xi\right)dz=
=S+αωcosωtsinklkl;=S+\alpha\omega\cos\omega t\;\frac{\sin kl}{kl}; (7)

and the acceleration takes place in the range z=[l,+l]z=[-l,+l]. Note that, for this resonant state, the expectation value of the z momentum Kz=(k/ω)E\langle K_{z}\rangle=(k/\omega)\langle E\rangle .

In order to understand the accessibility of the high energy states, let us calculate the root mean square of the rate of energy gain

dEdt¯=[12πππd(ωt)(dEdt)2]1/2=\bar{\frac{dE}{dt}}=\left[\frac{1}{2\pi}\int_{-\pi}^{\pi}d(\omega t)\left(\frac{d\langle E\rangle}{dt}\right)^{2}\right]^{1/2}=
=2ω2(ω2k2)1/2qAK(m2+K2+q2A2)1/2sinklkl.=\frac{\sqrt{2}\omega^{2}}{\left(\omega^{2}-k^{2}\right)^{1/2}}\frac{qAK_{\perp}}{\left(m^{2}+K_{\perp}^{2}+q^{2}A^{2}\right)^{1/2}}\frac{\sin kl}{kl}. (8)

From Eqs. (6-8), it is obvious that both the energy E\langle E\rangle and the r.m.s rate of energy increase (acceleration) are resonantly enhanced by the factor ω/(ω2k2)\omega/\sqrt{(\omega^{2}-k^{2})}. If we now invoke the standard relativistic dispersion relation for EM propagation in highly under dense plasmas (see for example Mahajan & Asenjo (2016) )

ω2k2=ωp21+q2A2/m2,\omega^{2}-k^{2}=\frac{\omega_{p}^{2}}{\sqrt{1+q^{2}A^{2}/m^{2}}}, (9)

where ωp=4πnq2/m\omega_{p}=\sqrt{4\pi nq^{2}/m} denotes the plasma frequency and nn is the number density of plasma particles, we may derive an explicit expression for the relativistic factor (associated with the high energy state) γ=E/mS/m\gamma=E/m\approx S/m,

γ=ωωp(1+q2A2m2)1/4(1+q2A2+K2m2)1/2;\gamma=\frac{\omega}{\omega_{p}}\left(1+\frac{q^{2}A^{2}}{m^{2}}\right)^{1/4}\left(1+\frac{q^{2}A^{2}+K_{\perp}^{2}}{m^{2}}\right)^{1/2}; (10)

the first two factors reflect resonant enhancement. If ω>>ωp\omega>>\omega_{p}, very high energy states become accessible even for moderate value of q2A2/m2q^{2}A^{2}/m^{2}. In terms of these parameters, the acceleration rate becomes

dEdt=2ω2ωpsinklklqAK(m2+q2A2)1/4(m2+K2+q2A2)1/2,\frac{dE}{dt}=\frac{\sqrt{2}\omega^{2}}{\omega_{p}}\frac{\sin kl}{kl}\frac{qAK_{\perp}\left(m^{2}+q^{2}A^{2}\right)^{1/4}}{\left(m^{2}+K_{\perp}^{2}+q^{2}A^{2}\right)^{1/2}}, (11)

which, in physical units, may be approximated as (sin(kL)/kL1,qA>mc2,Kc\sin(kL)/kL\simeq 1,qA>mc^{2},\hbar K_{\perp}c)

dEdt2ω2ωp(qAmc2)1/2Kc.\frac{dE}{dt}\simeq\frac{\sqrt{2}\hbar\omega^{2}}{\omega_{p}}\left(\frac{qA}{mc^{2}}\right)^{1/2}K_{\perp}c. (12)

This completes our review of the theory of resonant energization.

3 Radio AGNs and Processes potentially limiting energy gain

We are now ready to apply the essentials of this theory to the particular case of a plasma in the magnetospheres of a radio AGN that emits copious EM energy in the radio frequency range. Although we will, for this paper, work out the mechanism for the radio AGN, the resonant mechanism of particle energization pertains to any AGNs emitting in the entire electromagnetic spectrum- all the way to gamma rays.

To assess the effectiveness of the resonant energy transfer from the EM waves to the relativistic KG particle-waves, we must consider processes that will impede the acceleration process (summed up in Eq. (12). Perhaps, in the present context, the two most important impeding (cooling) processes will be :1) the inverse Compton (IC) scattering of the charged particles with the ambient photon field, and 2) synchrotron radiation when relatively strong magnetic fields are present. We will deal with them in the following subsections.

3.1 Maximum Energy Gain limited by Inverse Compton (IC) Scattering

We begin with the well-known expression for the IC cooling power (Rybicki & Lightman, 1979)

PCompσTcγ2U,P_{{}_{Comp}}\simeq\sigma_{{}_{T}}c\gamma^{2}U, (13)

where σT\sigma_{{}_{T}} is the Thomson cross section, U=L/(4πrc2)U=L/(4\pi rc^{2}) is the energy density of EM emission, LL denotes the bolometric luminosity of the radio source and rr is the distance from the central object. As an example we consider the typical length-scale, where the generation of very high energy (VHE) particles might take place - r0.001pcr\simeq 0.001pc.

As the particle relativistic factor (energy) increases during resonant energization, so would the energy loss due to IC (γ2)(\sim\gamma^{2}). Further acceleration will, therefore be terminated as soon as the rate of energy loss balances the rate of energy gain. Equating dE/dt¯PIC\bar{dE/dt}\simeq P_{{}_{IC}}, gives us the maximum allowed Lorentz factor

γCompmax(d¯Edt×4πrc2σTL)1/2\gamma_{{}_{Comp}}^{max}\simeq\left(\frac{\bar{d}E}{dt}\times\frac{4\pi rc^{2}}{\sigma_{{}_{T}}L}\right)^{1/2}\simeq
1.4×1015×(r0.001pc)5/8×\simeq 1.4\times 10^{15}\times\left(\frac{r}{0.001\;pc}\right)^{5/8}\times
×(f109Hz)3/4(1043erg/sL)3/8,\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\times\left(\frac{f}{10^{9}\;Hz}\right)^{3/4}\left(\frac{10^{43}\;erg/s}{L}\right)^{3/8}, (14)

where f=ω/2πf=\omega/2\pi is the frequency of EM radiation and we have used an expression of the Poynting flux in terms of the vector potential S=A2ω2/(4πc3)S=A^{2}\omega^{2}/(4\pi c^{3}). We normalize the frequency by 10910^{9} Hz (one should note that the frequency interval of radio AGN is 1071110^{7-11} Hz (Beckmann & Shrader, 2012)). In Eq. (14) we have considered a spherically symmetric accretion model, leading to the following number density of particles (Shapiro & Teukolsky, 2004)

Refer to caption
Figure 1: Here we show the plots of the maximum energy E(L) for different emission frequencies, f=(1;10;100)×109f=(1;10;100)\times 10^{9} Hz. The set of parameters is M=108×MM=10^{8}\times M_{\odot}, r=0.001r=0.001 pc, n1n_{\infty}\simeq 1 cm-3, T104T_{\infty}\simeq 10^{4} K and Γ=5/3\Gamma=5/3.
Refer to caption
Figure 2: The plots of particles’ energy versus distance form the central object. The set of parameters is the same as in Fig.1, except L=1043L=10^{43} erg/s.
nn24×(GMa2r)3/2,n\simeq n_{\infty}\frac{\sqrt{2}}{4}\times\left(\frac{GM}{a_{\infty}^{2}r}\right)^{3/2}, (15)

were n1n_{\infty}\simeq 1 cm-3 is the number density of protons in the interstellar gas (Shapiro & Teukolsky, 2004), GG is the gravitational constant; the black hole is assumed to be supermassive with M=108MM=10^{8}M_{\odot}, where M2×1033M_{\odot}\simeq 2\times 10^{33} g denotes the Solar mass, and

a=(ΓkBTm)1/2a_{\infty}=\left(\frac{\Gamma k_{B}T_{\infty}}{m}\right)^{1/2} (16)

is the speed of sound of the interstellar gas and Γ=5/3\Gamma=5/3 is its adiabatic constant, kBk_{B} is the Boltzmann’s constant and the temperature of the interstellar gas is assumed to be T104T_{\infty}\simeq 10^{4} K (Shapiro & Teukolsky, 2004).

3.2 Maximum Energy Gain limited by synchrotron emission

Another mechanism that potentially might limit the maximum attainable particle energy is the synchrotron radiation that is most effective when a relatively strong magnetic field is present. It is normally assumed that magnetic energy density, B2/(8π)B^{2}/(8\pi), and energy density of emission, UU are in equipartition: B2/(8π)UB^{2}/(8\pi)\simeq U. The estimated magnetic induction, then, will be

B2Lr2c8.4×0.001pcr×(L1043erg/s)1/2G.B\simeq\sqrt{\frac{2L}{r^{2}c}}\simeq 8.4\times\frac{0.001pc}{r}\times\left(\frac{L}{10^{43}\;erg/s}\right)^{1/2}\;G. (17)

The upper limit on the particle energy, imposed by synchrotron cooling, will be calculated ( as in the previous case) by balancing the resonant energization rate with

Psyn2e4B2γ23m2c3,P_{syn}\simeq\frac{2e^{4}B^{2}\gamma^{2}}{3m^{2}c^{3}}, (18)

the loss rate pertaining to the relativistic regime (Rybicki & Lightman, 1979). We find

γsynmax(dEdt×3m2c32e4B2)1/2\gamma_{syn}^{max}\simeq\left(\frac{dE}{dt}\times\frac{3m^{2}c^{3}}{2e^{4}B^{2}}\right)^{1/2}\simeq
1.9×1016×(f109Hz)3/4×\simeq 1.9\times 10^{16}\times\left(\frac{f}{10^{9}\;Hz}\right)^{3/4}\times
×(1043erg/sL)3/8(r0.001pc)9/8\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\times\left(\frac{10^{43}\;erg/s}{L}\right)^{3/8}\left(\frac{r}{0.001pc}\right)^{9/8} (19)

Since γsynmaxγCommax\gamma_{syn}^{max}\gg\gamma_{{}_{Com}}^{max}, the synchrotron losses will not impose ”any” meaningful restrictions on maximum energy that a particle could gain via resonant energization. The relevant upper bounds might come from the Inverse Compton cooling.

On the other hand, the maximum value of the Lorentz factor, that might be provided by the present mechanism is given by Eq. (10), which for realistic physical quantities with eA/mc2>>1eA/mc^{2}>>1 writes as

γωωp×(eAmc2)3/2\gamma\simeq\frac{\omega}{\omega_{p}}\times\left(\frac{eA}{mc^{2}}\right)^{3/2}\simeq
8.5×107(109Hzf)1/2×(L1043erg/s×0.001pcr)3/4.\;\;\;\;\;\;\;\simeq 8.5\times 10^{7}\left(\frac{10^{9}\;Hz}{f}\right)^{1/2}\times\left(\frac{L}{10^{43}\;erg/s}\times\frac{0.001pc}{r}\right)^{3/4}. (20)

From this expression, after comparing to γCompmax\gamma_{{}_{Comp}}^{max} it becomes evident that inverse Compton process might significantly restrict the maximum attainable energies only for higher values of γ\gamma. Therefore, the energization process for the mentioned physical quantities is fully governed by the resonance process of acceleration.

3.3 Conclusions-Discussion

Equation (20), giving an expression for the maximum energy constitutes the main result of this paper. Even more, since the mass dependance of the relativistic factor is of the form γ1/m\gamma\propto 1/m, the achieved energy γmc2\gamma mc^{2} does not depend on mass and therefore, the results are the same for protons and electrons.

For several representative frequencies (in the range where the AGN spectrum peaks), the total particle energy is plotted as a function of the AGN luminosity (Fig. 1). The set of parameters is f=(1;10;100)×109f=(1;10;100)\times 10^{9} Hz, M=108×MM=10^{8}\times M_{\odot}, r=0.001r=0.001 pc, n1n_{\infty}\simeq 1 cm-3, T104T_{\infty}\simeq 10^{4} K and Γ=5/3\Gamma=5/3. From the plots it is evident that EE is a continuously increasing function of luminosity, which is a natural result of the relativistic factor’s behaviour, γL3/4\gamma\propto L^{3/4}, (see Eq. (20)). One reads from Fig. 1 that the newly proposed mechanism of resonant energization, can accelerate particles to several EeVs.

In general, the number density of particles might be different for different regions of accretion matter (see Eq. (15)). In Fig.2 we plot the energy of particles versus the distance from the AGN. The luminosity is set L=1043L=10^{43} erg/sec; all other parameters are the same as in Fig. 1. It is clear form the plots that EE, as predicted by Eq. (20): Er3/4E\propto r^{-3/4}, is a continuously decreasing function of distance. Energies achieved by particles is of the order of 10162010^{16-20} eV, but for higher values of luminosity the achieved energies will reach several EeVs (see Fig. 1).

As it turned out, for physical quantities considered in this manuscript the Inverse Compton process does not impose significant restrictions on particle acceleration. But

We expect that this mechanism will operate in a variety of cosmic settings (in the vicinity of highly compact objects) and can account for some of the more energetic cosmic rays. In fact, there is a sister mechanism of resonant interaction operating directly between the KG and Gravitational waves that could result in a similar energy transfer to fast particles (Asenjo & Mahajan, 2020).

Acknowledgments

We are grateful to an anonymous referee for very interesting and valuable comments.

Data Availability

Data are available in the article and can be accessed via a DOI link.

References

  • Beckmann & Shrader (2012) Beckmann, V. & Shrader, C., 2012, Active Galactic Nuclei. Wiley-VCH
  • Bell (1978) Bell, A.R., 1978, MNRAS, 182, 147
  • Blandford & Znajek (1977) Blandford, R.D. & Znajek, R.L., 1977, MNRAS, 179, 433
  • Catanese & Weeks (1999) Catanese, M. & Weeks, T. C., 1999, PASP, 111, 1193
  • Fermi (1949) Fermi, E., 1949, Phys. Rev. 75, 1169
  • Kim & Kim (2013) Kim H. B. & Kim J., 2013, Int. J. Mod. Phys. D, 22, 1350045
  • Lemoine (2021) Lemoine, M., 2021, Phys. Rev. D, 104, 063020
  • Mahajan & Asenjo (2022) Mahajan, S.M.& Asenjo, F.A., Phys. Plasmas, 2022, 29, 022107
  • Mahajan & Asenjo (2016) Mahajan, S.M.& Asenjo, F.A., 2016, Phys. Plasmas, 23, 056301
  • Osmanov (2010) Osmanov, Z., 2010, New Astronomy, 15, 351
  • Osmanov et al. (2017) Osmanov, Z., Mahajan, S. & Machabeli, G., 2017, ApJ, 835, 4
  • Osmanov et al. (2014) Osmanov, Z., Mahajan, S., Machabeli, G. & Chkheidze, N., 2014, MNRAS, 445, 4155
  • Osmanov et al. (2007) Osmanov, Z., Rogava, A. & Bodo, G., 2007, A&A, 470, 395
  • Pezzi et al. (2022) Pezzi, O., Pasquale, B. & Mannheim, K., 2000, A&A, 353, 473
  • Rieger & Mannheim (2000) Rieger, F.M. & Matthaeus, W.H., 2022, ApJ, 928, 25
  • Rybicki & Lightman (1979) Rybicki G.B. & Lightman A.P., 1979, Radiative Processes in Astrophysics. Wiley, New York
  • Shapiro & Teukolsky (2004) Shapiro S.L. & Teukolsky, S.A., 2004, Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects. Wiley-VCH, Weinheim
  • Asenjo & Mahajan (2020) Asenjo, F.A.& Mahajan, S.M., 2020, Phys. Rev.D,101, 063010 (2020)