Resonant energization of particles by radio AGN
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 eV.
Key Words.:
acceleration of particles - galaxies: active - galaxies: jets - plasmas1 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),
(1) |
approaches the speed of light () when , where is the momentum of the particle, its relativistic factor and - 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 (),
(2) |
where and are, respectively, the frequency and the wavenumber of the EM wave propagating along the axis. Because, the directions perpendicular to z are ignorable, the KG/EM system obeys the Klein-Gordon equation (see Mahajan & Asenjo (2016) for details),
(3) |
where is the particle charge, and 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 and dependence in Eq. (3) implies that the energy E and will be functions of time. The resonant solutions of Eq. (3) will emerge when we demand (); the result is a simple Mathew equation
(4) |
with and . Notice that the mathematical operator in Eq. (3) translates to (in Eq. (4) when we impose a solution whose phase factor ()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 tends to zero for EM waves traveling in a tenuous medium. The implication, of course, is that nontrivial solutions must demand derivative of to become commensurately large. Consequently, the energy and z momentum of the KG state must also become large since both are proportional to the derivative (as we will soon see).
Following simple rules of quantum mechanics, the expectation value of the energy is
(7) |
and the acceleration takes place in the range . Note that, for this resonant state, the expectation value of the z momentum .
In order to understand the accessibility of the high energy states, let us calculate the root mean square of the rate of energy gain
(8) |
From Eqs. (6-8), it is obvious that both the energy and the r.m.s rate of energy increase (acceleration) are resonantly enhanced by the factor . If we now invoke the standard relativistic dispersion relation for EM propagation in highly under dense plasmas (see for example Mahajan & Asenjo (2016) )
(9) |
where denotes the plasma frequency and is the number density of plasma particles, we may derive an explicit expression for the relativistic factor (associated with the high energy state) ,
(10) |
the first two factors reflect resonant enhancement. If , very high energy states become accessible even for moderate value of . In terms of these parameters, the acceleration rate becomes
(11) |
which, in physical units, may be approximated as ()
(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)
(13) |
where is the Thomson cross section, is the energy density of EM emission, denotes the bolometric luminosity of the radio source and 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 - .
As the particle relativistic factor (energy) increases during resonant energization, so would the energy loss due to IC . Further acceleration will, therefore be terminated as soon as the rate of energy loss balances the rate of energy gain. Equating , gives us the maximum allowed Lorentz factor
(14) |
where is the frequency of EM radiation and we have used an expression of the Poynting flux in terms of the vector potential . We normalize the frequency by Hz (one should note that the frequency interval of radio AGN is 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)


(15) |
were cm-3 is the number density of protons in the interstellar gas (Shapiro & Teukolsky, 2004), is the gravitational constant; the black hole is assumed to be supermassive with , where g denotes the Solar mass, and
(16) |
is the speed of sound of the interstellar gas and is its adiabatic constant, is the Boltzmann’s constant and the temperature of the interstellar gas is assumed to be 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, , and energy density of emission, are in equipartition: . The estimated magnetic induction, then, will be
(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
(18) |
the loss rate pertaining to the relativistic regime (Rybicki & Lightman, 1979). We find
(19) |
Since , 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 writes as
(20) |
From this expression, after comparing to it becomes evident that inverse Compton process might significantly restrict the maximum attainable energies only for higher values of . 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 , the achieved energy 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 Hz, , pc, cm-3, K and . From the plots it is evident that is a continuously increasing function of luminosity, which is a natural result of the relativistic factor’s behaviour, , (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 erg/sec; all other parameters are the same as in Fig. 1. It is clear form the plots that , as predicted by Eq. (20): , is a continuously decreasing function of distance. Energies achieved by particles is of the order of 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)