Properties of Particle Trajectory Around a Weakly Magnetized Black Hole
Haldar Amritendu1, †, Biswas Ritabrata2, ‡
Department of Physics, Sripat Singh College, Jiaganj, Murshidabad 742123, India
Department of Mathematics, University of Burdwan, Golapbag Academic Complex, Burdwan 713 104, WestBengal, India

 

† Corresponding author. E-mail: amritendu.h@gmail.com biswas.ritabrata@gmail.com

Project grant of Goverment of West Bengal, Department of Higher Education, Science and Technology and Biotechnology, File no:- ST/P/S&T/16G-19/2017, AH wishes to thank the Department of Mathematics, the University of Burdwan for the research facilities provided during the work. RB thanks IUCAA, Pune, India for providing Visiting Associateship

Abstract
Abstract

In this paper, we consider charged accelerating AdS black holes with nonlinear electromagnetic source. The metric chosen by us is of a regular black hole, which shows regular nature at poles and a conical effect, which corresponds to a cosmic string. In such a space time construction of the Lagrangian for a charged particle is done. Cyclic coordinates as well as the corresponding symmetry generators, i.e., the Killing vectors are found. Conservation laws corresponding to the symmetries are counted. Euler-Lagrange equations are found. The orbit is mainly taken to be a circular one and effective potential is found. The minimum velocity obtained by a particle to escape from innermost stable circular orbit is found. The value of this escape velocity is plotted with respect to the radius of the event horizon of the central black hole for different parametric values. The nature of the escape velocity is studied when the central object is working with gravitational force and charge simultaneously. Effective potential and effective force are also plotted. The range of radius of event horizon for which the effective force turns to be positive is found out. A pathway of future studies of accretion disc around such black holes is made.

1 Introduction

Black Hole (BH hereafter) is the important prediction of General Relativity (GR hereafter). BHs are specified by the curved space-time geometry and bounded by the event horizons. Nowadays, the study of dynamics of particles around the BHs in the astrophysical background is found in many literature. Both the particles and photons in the vicinity of the BHs are highly attracted by the strong gravitational pull and as a result the accreted particles cause the gain of mass of the BHs. Again it is also possible that the BHs throw the particles away with high relativistic velocity due to angular momentum barrier. Furthermore the highly energised charged particles may also be escaped from their stable orbits around the BHs due to the collisions with other particles and are moved under the influence of the Lorentz force in the electromagnetic fields.

If the gravitational and electromagnetic fields are strong, the motions of the charged particles become, in general, chaotic and undetermined and so the corresponding orbits become unstable except the inner most stable circular orbit. The weak magnetic field does not affect the geometry of the concerned central BH but it may affect the motion of charged particles.[12] In the experimental point of view,[34] the magnetic field around the BHs occurs due to existence of plasma in the form of an accretion disc or a charged gas cloud[56] surrounding the BHs.

Mechanically, when the particles of a system collide with each other, the Center of Mass Energy (CME hereafter) occurs. Near the event horizon of a BH, the high CME is produced due to the collision of two particles. The mechanism of collision of two particles falling towards the Kerr BH has been proposed by Banados Silk and West (BSW mechanism hereafter).[7] Furthermore, the same authors have also shown that the CME, in the equatorial plane, may be highest for a fastly rotating BH. The BSW mechanism has been investigated for different BHs.[818] A general review of the mechanism of collision is viewed in Ref. [19]. The CME of the particles at the inner horizon of Kerr BH,[20] the CME of the collision of the particles around Kerr-Newmann BH,[21] near the horizon(s) of Kerr-Taub-NUT BH,[22] cylindrical BH,[23] Plebanski-Demianski BH,[24] Kerr-Newmann-Taub-NUT BH,[25] charged dilaton BH[26] have been studied.

The study of the dynamics of particles around the BHs in the presence of magnetic field has a special significance in recent time. The authors of some literatures[2729] have investigated the motion of a charged particle near weakly magnetized Schwarzschild BH. In some other literatures[3035] the authors have analyzed the chaotic motion of a charged particles around Kerr BH near magnetic field. Moreover in another paper[36] the authors have studied the circular motion of charged particles around Reissner-Nordstrom BH. Reference [37] has investigated the dynamics of charged particles around slowly rotating Kerr BH with magnetic field and also Ref. [38] has discussed about the dynamics of particles around Schwarzschild BH in the presence of quintessence and magnetic field. Recently, Ref. [39] has investigated the dynamics of particles around a regular BH surrounded by external magnetic field. Zhou et al.[40] has studied the geodesic structure of the Janis-Newman-Winicour space time. This space time contains a strong curvature naked singularity. Solving the geodesic equation and analysing the behavior of effective potential, they have shown all geodesic types of the test particles and photons. In 2009 Horava Lifshitz gravity was introduced. Timelike geodesic motion in Horava Lifshitz space time was studied in Ref. [41]. Differently energised particles’ motions are studied.

Motivated by the previous works, our objective in this paper is to apply the Euler-Lagrange equations of motion to study the properties of particle trajectories around the charged accelerating AdS (Anti-de-Sitter) BH surrounded by weak magnetic field.

This paper is organized as follows: in next section, we study a charged accelerating AdS BH metric with non-linear electromagnetic source and establish the Euler-Lagrange equation of motion radially. In Sec. 3, we analyze the dynamics of charged particles with and without magnetic field around a charged accelerating AdS BH. In Sec 4, we investigate the CME of two colliding charged particles in the aforesaid conditions and in Sec. 5, we calculate the effective force and explain it with graphical representation. Finally in the last section, we present a brief conclusion of the work.

2 Charged Accelerating AdS (Anti-de-Sitter) Black Holes with Non-Linear Electromagnetic Source

A charged accelerating AdS BHs can be expressed by the metric as[4243]

with the electromagnetic field tensor F, which is related to the gauge potential B as:

and where

The factor Ω determines the conformal infinity or boundary of the AdS space-time. In the above expression (3), m, q, (where Λ is cosmological constant) and A represents the mass, electric charge of BHs, AdS radius and magnitude of acceleration of the BHs respectively. By looking at the angular part of the metric and the behavior of g(θ) at the both poles θ+ = 0 and θ = π, one can discover the presence of “cosmic string”. The regularity of the metric at a pole requires,

We assume that the BH is regular on the North Pole (θ = 0) with K = 1 + 2Am+q2A2 and on the South Pole (θ = π) with K = 1 − 2Am + q2A2 but there is a “conical defect”, which is δ = 2π (1 − g/K+) = 8πmA/(1 − 2Am + q2A2), which corresponds to a “cosmic string” with tension μ = δ/8π = πmA/(1 − 2Am + q2A2).

3 Dynamics of Charged Particles
3.1 In the Absence of Magnetic Field

The Lagrangian of a charged particle of mass m and charge q in the absence of magnetic field is expressed as[29,38]

By employing the Lagrangian dynamics, we now examine the motion of that particle in the background of the BH (1). The Lagrangian of that particle is as follows:

The over dot indicates the derivative with respect to proper time. It is evident from the Lagrangian (6) that t, θ, and ϕ are cyclic coordinates and hence it leads to corresponding symmetry generators, which are also known as Killing vectors. The metric (1) remains invariant under the Killing vector fields X and we obtain:

There are three commuting integral of motion corresponding Eq. (1) in which two of them are generated by the Killing vectors

where and , reflecting invariance under the time translation and rotation around the symmetry axis. The corresponding conserved quantities, viz., specific energy E and specific azimuthal angular momentum Lz are given by

where pμ = muμ is the linear momentum of the particle in the absence of magnetic field. And the other conservation laws corresponding to the symmetries are depicted as: angular momenta

It is obvious from Eqs. (9) and (10) that if the particle, which moves in the equatorial plane, X2 and X3 become irrelevant.

Here we apply the Euler-Lagrange equation of motion for r only and we obtain

where f(r), , and are obtained from Eqs. (3) and (12).

The constant of motion for that particle is given by Eq. (12). The particle is associated with the total specific angular momentum as[3,39]

Applying the normalization condition for four-velocity, we calculate from Eq. (6) as

Since we assume the system is spherically symmetric, all θ = constant planes will be equivalent to the equatorial plane for which θ = π/2 and hence from Eq. (16) we have

Again if we consider the orbit of the particle through, which it is moving is circular, then Eq. (17) reduces to

and it is equivalent to the effective potential Ueff(r) of that particle. Hence Eq. (18) shows that the total energy and/or the effective potential will vanish at the horizon(s).

The critical azimuthal angular momentum of a particle that follows a particular orbit where the effective potential is extremum (i.e., maximum or minimum) and it is expressed as

Hence the energy of the particle would be

After collision, the energy of the particle takes the form as:

So comparing between Eqs. (18) and (21) we infer that the energy after collision is grater than that of the total energy before collision as the extra term rv is present in Eq. (21) and, which is obvious due to the collision. In this expression v is the minimum velocity of the particle required to escape from innermost stable circular orbit (ISCO hereafter) and is given as:

3.2 In the Presence of Magnetic Field

We now consider the case of weakly magnetized BH and investigate the motion of a charged particle having charge q in presence of magnetic field in the BH exterior. The general Killing vector equation is[44]

where ξμ is a killing vector and this equation coincides with the Maxwell equation for 4-potential Aμ in the Lorentz gauge . The special choice[45]

corresponds to the test magnetic field, where is magnetic field strength. The 4-potential discussed in Eq. (24) is invariant under the symmetries corresponding to the Killing vectors as discussed above, i.e.,

A magnetic field vector with respect to an observer whose 4-velocity is uμ, defined as:

where , , and g = det(gμν). ϵμνλσ is the Levi symbol and Fλσ is the Maxwell tensor, which is given by

The Lagrangian of a particle of mass m and charge q in the presence of magnetic field is expressed by[29,38]

where Aμ is the 4-vector potential for the electromagnetic field.

By employing the Lagrangian dynamics, we now examine the motion of a charged particle in the background of the BH (1). Here the Lagrangian of that particle is as follows:

The generalized 4-momentum of the particle is denoted as,

Then the new conservation laws corresponding to the symmetries are defined below

Here we denote . Hence the new constants of motion will be given by

The equation of motion of the charged particle in this case is obtained by applying by the Euler-Lagrange equation, which will be of the form as:

where f(r), and are obtained from Eq. (32). Applying the normalisation condition for four-velocity in Eq. (29) and using Eq. (32), we have the energy and the corresponding effective potential as:

For the condition as applied in Eq. (16), we have from Eqs. (34) and (35) that

3.3 Dimensionless Form of the Equations

In order to integrate the dynamical equations, we need to make these equations dimensionless. We use the following transformation relations[38,4647] as:

Using these relations (38), Eqs. (33), (34), and (35) acquire the form as:

where

After collision, for θ = π/2 and constant ρ the energy given in Eq. (39) reduces to

Hence the escape velocity of the charged particle is expressed as:

We have plotted the variation of escape velocity v of the particle with respect to the radius r of the BH for different values of q, l, K, b, Lz, and Ec in Figs. 1(a)1(e)

Fig. 1 The variation of escape velocity v with respect to r (a) K, Lz and b are fixed, with varying q and l. (b) q, l, Lz and b are fixed with varying K. (c) q, l, K and b are fixed with varying Lz. (d) q, l, Lz and K are fixed with varying b. (e) K, Lz q, l and b are fixed, with varying Ec.

In Fig. 1(a), we vary q and l simultaneously keeping K, b, Lz, and Ec fixed and we notice that with increment of q and l, for same r, the escape velocity v decreases. If r is small enough, v has a negative value, which depicts no physical particle can escape from so near points of a compact object considered here. But once distance from the center of the gravitating object is made higher, we can have an escape velocity, which firstly will increase steeply with the increment of r and latter the rate of increment will be reduced. However, the value is ever increasing with r. If we increase q and l keeping all the other parameters constant, the curves stay of the same nature only except the fact that they are amplified. This means, for a higher charge, escape velocity is lower. So is for the higher space time curvature, i.e., higher the curvature lower is the escape velocity. This phenomenon is quite obvious as whenever the central engine is attracting gravitationally and electrically it is more hard for the particle to escape from such an object. Besides, if the curvature is high, higher velocity is required.

In Fig. 1(b), we change the values of K by fixing the other parameters and we find that with increase of K, the particle escapes from higher values of r but for same r, the escape velocity v increases. As we vary K only again for low r we see unphysical escape velocity. With r escape velocity increases. For low r, low K a high v will be required. For high r, low K, the v required is low. For every two Klow and Khigh curves there is a point rcrit where the system needs same escape velocity. This is due to the quadratic nature of K in the expression of v in Eq. (42).

We vary Lz keeping the parameters q, l, K, b, and Ec are unchanged in Fig. 1(c) and we observe that its nature is similar as Fig. 1(b). Here we noticed that with increase of Lz, the particle escapes from higher values of r but for same r, the escape velocity v decreases.

Figures 1(d) and 1(e) have the similar nature. One is plotted by varying the strength of magnetic field b and other with escape energy Ec. We observe from the figures that for higher strength of magnetic field as well as escape energy the escape velocity of the particle is also higher. Hence we conclude that the magnetic field, which escape the particle from the vicinity of the BH plays an important role in transfer mechanism of energy.

We have drawn the curves of effective potential Ueff of the particle corresponding to circular orbits with respect to radius r of BH for different values of q, l, K, b, and Lz in Figs. 2(a)2(d). We observe a similar feature that there is a minimum value of Ueff at r lying between 1 and 2. Initially Ueff decreases very sharply (almost straight down) with increases of r and gets negative value when r exceeds the value about 1 and it reaches to a minimum value. Further increases of r, Ueff increases to nearly zero and then rapidly decreases.

Fig. 2 The variation of effective potential Ueff with respect to r (a) K, Lz, and b are fixed, with varying q and l. (b) q, l, Lz, and b are fixed with varying K. (c) q, l, Lz, and K are fixed with varying b. (d) q, l, K, and b are fixed with varying Lz.

In Fig. 2(a), it is found that Ueff gets nearly zero value at about r = 3 for all the variation of q and l. But Figs. 2(b)2(d) show that Ueff reaches nearly zero value for different values of r. Due to increase of corresponding parameters, the values of r where Ueff reaches nearly zero value also increase.

4 Center of Mass Energy of the Colliding Charged Particles

The CME of the colliding particles is expressed as[39]

where m0 is the mass and uμ is the 4-velocity of each particles respectively.

4.1 In the Absence of Magnetic Field

Applying Eq. (12) in Eq. (43), we obtain

4.2 In the Presence of Magnetic Field

Applying Eq. (32) in Eq. (43), we obtain

It is evident from Eqs. (44) and (45) that the CME of the charged particles with and without magnetic field do not change at horizon(s).

5 Effective Force

The effective force acting on the charged particles in the flat background normally measured by the Lorentz force but in the curved background it may be determined as follows

We have plotted here the graphs of effective force feff acting on the charged particle vs. the radius r of BH for different values of q, l, K, b, and Lz in Figs. 3(a)3(d). We observe from all the graphs that for r = 0.3 feff has very high value, further very small increases of r, feff decreases sharply and reaches a minimum value at r = 0.4. feff acquires zero value i.e., constant effective potential at about r = 0.6. Figures 3(a) and 3(c) are in similar nature. Both the graphs show for higher values of q and l (Fig. 3(a)) and b (Fig. 3(c)), feff increases more rapidly with increases of r. Again surprisingly the feature of Figs. 3(b) and 3(d) are similar. Here we vary K and Lz respectively. Both these curves show feff reaches to Positive value within the range of r about 0.6 to 1.4. For higher values of K and Lz the peak value of feff is also high. Beyond r = 1.4, feff decreases slightly and then increases rapidly with increases of r this increasing is higher for higher values of K and Lz respectively.

Fig. 3 The variation of effective force feff with respect to r (a) K, Lz, and b are fixed, with varying q and l. (b) q, l, Lz and b are fixed with varying K. (c) q, l, Lz and K are fixed with varying b. (d) q, l, K and b are fixed with varying Lz.
6 Conclusion

In our study, we investigate the dynamics of charged particles around charged accelerating AdS (anti-de-Sitter) black hole in the absence and presence of magnetic field by computing Euler-Lagrange equation of motion for radial component only. First of all we choose the charged accelerating AdS (anti-de Sitter) BH with non-linear electromagnetic source and analysis the effective potential, escape energy and therefore escape velocity in the absence and presence of magnetic field. We analysis the dependence of that physical quantities on the radius of BH graphically. We observe from Figs. 1(a)1(e) that the escape velocity of the charged particle increases with increase of the radius of BH as expected. But higher values of charge, Ads radius (Fig. 1(a)), magnetic field (Fig. 1(d)) and escape energy (Fig. 1(e)) the escape velocity for same radius decreases. Whereas the escape velocity for same radius decreases initially for higher values of azimuthal angular momentum (Fig. 1(b)) and cosmic spring constant (Fig. 1(c)). But after a certain range of radius the escape velocity increases for for higher values of azimuthal angular momentum (Fig. 1(b)) and cosmic spring constant (Fig. 1(c)). We observe from Figs. 1(d) and 1(e) that for higher strength of magnetic field as well as escape energy the escape velocity of the particle is also higher. Hence we conclude that the magnetic field, which escape the particle from the vicinity of the BH plays an important role in transfer mechanism of energy.

The variation of effective potential with radius of BH is studied in Figs. 2(a)2(d). All these graphs show that there is a minimum negative value when radius is very small. Maybe this is a stable circular orbit. There also may exist many stable circular orbits. Moreover we investigate the CME of two colliding particles with and without magnetic field.

Finally, we calculate the effective force acting on the charged particles in the curved background and examine its dependence on the radius of black bole graphically. We find there is a minimum negative value of effective force at very small radius of BH and then it reaches to zero, stay within few range of the radius for different values of charge and AdS radius (Fig. 3(a)) and magnetic field Fig. 3(c)). Again for different values of cosmic spring constant (Fig. 3(b)) and azimuthal angular momentum (Fig. 3(d)) it reaches to a positive value whose peak value is higher for higher values cosmic spring constant and azimuthal angular momentum separately.

We have studied the dynamics of both charged and chargeless particles. It is seen that more the charge, lower is the escape velocity. Dependency of escape velocity on K is somewhere different. For the quadratic nature of K, we see between two KLow and Khigh curves, there is a point rcrit where the system needs same escape velocity. In a nutshell, it is shown in this paper that as we shift from the Schwarzschild black hole and incorporate charge and magnetic field, the attracting nature increases and particles near to the central engine becomes more redshifted. If we compare our result with existing studies like[39] we observe their speculation is that the particle moving around such black holes should not feel the absence of the singularity of the central regular black hole. The effective part, according to them, must be the magnetic field present in the accretion disc. However, our study also indicates the same. Increment in magnetic field strength signifies that the rotating particle needs comparatively lesser value of escape velocity. Dynamics of a neutral and a charged particle around a Schwarzschild singularity embedded in quintessence universe has been studied in Ref. [38]. Energy conditions for stable orbits has been constructed. Their result has been more clearly stated in our work. In future, if the motion of a series of particles are been counted, then we will have an idea of accretion disc around such a regular black hole.

Reference
[1] Znajek R. Nature (London) 262 1976 270
[2] Blandford R. D. Znajek R. L. Mon. Not. R. Astron. Soc. 179 1977 433
[3] Frolov V. The Galactic Black Hole Falcke H. Hehl F. H. IoP London 2003
[4] Borm C. V. Spaans M. Astron. Astrophys. L 9 2013 553
[5] Mckinney J. C. Narayan R. Mon. Not. R. Astron. Soc. 375 2007 523
[6] Dobbie P. B. Kuncic Z. Bicknell G. V. Salmeron R. Proceeding of IAU Symposium 259 Galaxies Tentefire 2008
[7] Bandos M. Silk J. West S. M. Phys. Rev. Lett. 103 2009 111102
[8] Berti E. Cardoso V. Gualtieri L. et al. Phys. Rev. Lett. 103 2009 239001
[9] Patil M. Joshi P. S. Phys. Rev. D 82 2010 104049
[10] Patil M. Joshi P. S. Classical Quantum Gravity 28 2011 235012
[11] Wei S. W. Liu Y. X. Guo H. Fu C. E. Phys. Rev. D 82 2010 103005
[12] Zaslavskii O. B. JETP Lett. 92 2010 571
[13] Harada T. Kimura M. Phys. Rev. D 83 2011 024002
[14] Grib A. A. Pavlov V. Yu. Astropart. Phys. 34 2011 561
[15] Hussian I. Mod. Phys. Lett. A 27 2012 1250017
[16] Sadeghi J. Pourhassan B. Eur. Phys. J. C 72 2012 1984
[17] Bambi C. Modesto L. Phys. Lett. B 721 2013 329
[18] Tuesunov A. Kolo.s M. Abdujabbarov A. et al. Phys. Rev. D 88 2013 124001
[19] Harada T. Kimura M. Classical Quantum Gravity 31 2014 243001
[20] Lake K. Phys. Rev. Lett. 104 2010 211102
[21] Wei S. W. Liu Y. X. Li H. T. Chen F. W. J. High Energy Phys. 12 2010 066
[22] Liu C. Chen S. Ding C. Jing J. Phys. Lett. B 701 2011 285
[23] Said J. L. Adami K. Z. Phys. Rev. D 83 2011 104047
[24] Sharif M. Haider N. J. Theor. Exp. Phys. 117 2013 78
[25] Zakria A. Jamil M. J. High Energy Phys. 147 2015 1505
[26] Pradhan P. Astropart. Phys. 62 2015 217
[27] Forlov V. P. Shoom A. A. Phys. Rev. D 82 2010 084034
[28] Forlov V. P. Phys. Rev. D 85 2012 024020
[29] Zahrani A. A. Frolov V. P. Shoom A. A. Phys. Rev. D 87 2013 084043
[30] Nakamura Y. Ishizuka T. Astrophys. Space Sci. 210 1993 105
[31] Takahashi M. Koyama H. Astrophys. J 693 2009 472
[32] Kopacek O. Kovar J. Karas V. Stuchlik Z. AIP Conf. Proc. 1283 2010 278
[33] Preti G. Phys. Rev. D 81 2010 024008
[34] Forlov V. P. Krtous P. Phys. Rev. D 83 2011 024016
[35] Igata T. Koike T. Ishihara H. Phys. Rev. D 83 2011 065027
[36] Pugliese D. Quevedo H. Rufini R. Phys. Rev. D 83 2011 104052
[37] Majeed B. Hussain S. Jamil M. Advances in High Energy Physics 2015 2015 6712559
[38] Jamil M. Hussian S. Majeed B. Eur. Phys. J. C 75 2015 24
[39] Jawad A. Ali F. Jamil M. Debnath U. arXiv:1610.07411[gr-qc] 2016
[40] Zhou S. et al. Inter. J. Theor. Phys. 54 2015 2905
[41] Chen J. H. et al. Inter J. Modern Phys. A IJMPA 25 2010 1439
[42] Griffith J. b. Podolsky J. Int. J. Mod. Phys. D 15 2006 335
[43] Liu Hang Meng Xin-He arXiv:1607.00496[gr-qc] 2016
[44] Wald R. M. Phys. Rev. D 10 1974 1680
[45] Aliev A. N. Ozdemir N. MNRAS 336 1978 241
[46] Hsu S. D. H. Phys. Lett. B 594 2004 13
[47] Hussain S. Hussain I. Jamil M. Eur. Phys. J. C 74 2014 3210