Circular geodesics in the field of double-charged dilatonic black holes

A non-extreme dilatonic charged (by two ``color electric'' charges) black hole solution is examined within a four-dimensional gravity model that incorporates two scalar (dilaton) fields and two Abelian vector fields. The scalar and vector fields interact through exponential terms containing two dilatonic coupling vectors. The solution is characterized by a dimensionless parameter $a$ $(0<a<2)$, which is a specific function of dilatonic coupling vectors. The paper presents solutions for timelike and null circular geodesics that may play a crucial role in different astrophysical scenarios, including quasinormal modes of various test fields in the eikonal approximation. For $a = 1/2, 1, 3/2, 2$, the radii of the innermost stable circular orbit are presented and analyzed.


Introduction
Since the discovery of gravitational waves in 2016, interest in the physics of stellar-mass black holes has rapidly increased [1].Gravitational waves were initially postulated by Einstein over a century ago.The extensive search for them resulted in the establishment of large, supersensitive laser interferometers such as LIGO, VIRGO, and others, along with the development of highly precise methods for detection and data analysis [2].
Observations and studies of the motion of stars at the heart of the Milky Way galaxy have confirmed the exa e-mail: kuantay@mail.rub e-mail: sulieva.gulnara0899@gmail.comc e-mail: ivashchuk@mail.rud e-mail: y.a.a.707@mail.ruistence of a supermassive black hole (BH) [3][4][5].These surveys provide information about the star clusters near the galactic center and the BH mass, which researchers use to test the predictions of general relativity (GR) [6].
In addition, the imaging process of the shadow of the supermassive black holes at the center of the M87 galaxy and in the Milky Way galaxy has required remarkable ingenuity and cohesion among scientists around the world [7,8].For this purpose, our entire planet has been used as one giant radio telescope, consisting of separate groups of telescopes scattered across all continents.The presence of the image of a BH shadow and its analysis verify the correctness and reliability of GR.All discoveries over the past 60 years testify to the significance and relevance of research in this direction.
At the same time, in addition to GR, various modified and extended theories of gravity exist, giving rise to alternative black holes with additional parameters [9][10][11][12][13][14][15].However, within the range of observational error bars, these black holes cannot be completely distinguished from ordinary Schwarzschild, Reissner-Nordström and Kerr black holes [16].
Observational data commonly serves as an important tool for constraining BH parameters in these theories [17][18][19][20].This fact allows for the examination of black holes with "colored charges" in the presence of scalar fields.This study considers classical Schwarzschild and Reissner-Nordström black holes, as well as various models of dilatonic black holes with two colored electric charges.
The dilaton is a hypothetical scalar field particle that appears in the metric through the introduction of coupling vectors (constants).This, in turn, defines a new particular subclass of black holes.In limiting cases, dilatonic solutions reduce to the Schwarzschild and Reissner-Nordström solutions, depending on the values of the coupling vectors.Given that the dilaton has not yet been discovered, it would be interesting to investigate the effects associated with its presence.
In this paper, we explore the possibility of distinguishing ordinary astrophysical black holes from dilatonic ones.Accordingly, the main objective of the paper is to study the motion of neutral test particles and photons in circular orbits in the gravitational field of astrophysical and dilatonic black holes with different color charges and coupling vectors.
To achieve this goal, the following problems are posed: -Deriving geodesic equations in the field of dilatonic black holes using the Lagrange formalism.-Calculating the angular momentum, energy, effective potential of neutral test particles and photons, and the radii of the innermost stable circular orbits (ISCO) in the field of dilatonic black holes.
The solutions of dilatonic dyonic black holes with an arbitrary coupling constant and a canonical scalar field were considered in Refs.[21][22][23][24].These solutions were obtained by solving two master equations for moduli functions.Physical parameters of the solutions were also derived, including gravitational mass, scalar charge, Hawking temperature, black hole area entropy, and parameterized post-Newtonian (PPN) parameters.While Refs.[21][22][23] focused solely on the physical characteristics of black holes, Ref. [24] explored quasi-normal modes of a massless test scalar field in the background of a gravitational field for a non-extremal dilatonic dyonic black hole.
Timelike and null geodesics 1 , including circular ones, play a crucial role in various astrophysical contexts, such as accretion disks, quasiperiodic oscillations, quasinormal modes of various test fields in the eikonal approximation [26], and shadows of supermassive black holes [27,28].These results and studies aim to distinguish ordinary black holes from dilatonic ones [17].
The novelty of the work lies in the exploration of geodesics for test particles in the gravitational field of dilatonic black holes with two scalar fields and two-color electric charges, marking the first study.Previous research conducted by some of us has focused on geodesics in the context of standard black holes, naked singularities, and 1 The geodesic motion in Euclidean Schwarzschild geometry was explored in [25].The explicit form of geodesic motion was derived using incomplete elliptic integrals of the first, second, and third kind.
In Ref. [38], the study focuses on geodesics around rotating black hole mimickers represented by the exact solution of the stationary and axially symmetric field equations in vacuum, known as the δ -Kerr metric.The authors study its optical properties using a ray-tracing code for photon motion, analyze the apparent shape of the shadow of a compact object, and compare it with a Kerr black hole.In order to provide qualitative estimates related to the observed shadow of a supermassive compact object in the M87 galaxy, the authors consider values of the object's spin a and observation inclination angle close to the measured values.The study demonstrates that, based on only one set of shadow edge observations, it is not possible to rule out the δ -Kerr solution as a viable source of the geometry outside the compact object.
It was shown in Ref. [39] that from a geometric perspective, distinguishing electrically and magnetically charged Reissner-Nordström black holes is impossible.To elucidate the distinctions between these solutions, one approach is to examine the dynamic motion of charged test particles in the vicinity of a charged black hole and scrutinize the impact of charge coupling parameters on the stability of circular orbits.The authors delve into the synchrotron radiation emitted by charged particles accelerated by a charged black hole and provide estimates for the intensity of the relativistic radiation emitted by these particles.
Recently, in Ref. [40], the authors investigated closed photon orbits in spherically symmetric static solutions of supergravity theories, a Horndeski theory, and a theory of quintessence.These orbits lie in what is called a photon sphere (or anti-photon sphere) if the orbit is unstable (stable).It was shown that in all the asymptotically flat solutions examined, which admit a regular event horizon and have an energy-momentum tensor satisfying the strong energy condition, there exists one and only one photon sphere outside the event horizon.
Note that while Refs.[41][42][43][44] explore the features of dilatonic black holes, Refs.[38,39,45] consider various solutions to the field equations that, in the limiting case, describe ordinary black holes.Refs.[38][39][40]45] also explore geodesics, which may differ in the field of ordinary black holes.Some works explore the motion of test particles and photons around static [46][47][48][49][50] and rotating [51,52] dilatonic black holes in different theories of gravity.To the best of the authors' knowledge, geodesics around dilatonic black holes with two scalar fields and two vector fields have not been studied elsewhere in the literature.
It should be mentioned that in classical physics, gravity is associated with the mass of an object, and it does not influence charge.As a result, the effects of gravitational and electromagnetic fields do not intertwine.However, in GR, gravity, in addition to mass, can be generated by rotation due to the rotational kinetic energy of the source (the Lense-Thirring effect), the presence of electric or magnetic (if they exist) charges, and the energy of their electromagnetic fields, as well as by other types of fields.Consequently, in GR, gravity affects the motion of both neutral and charged particles.These effects become more pronounced in strong field regimes near black holes and neutron stars, where the curvature of spacetime becomes substantial.
Moreover, neutral and charged particles may exhibit distinct behavior in their motion.For instance, charged particles deviating from geodesic motion may undergo significant changes in their trajectories due to the curvature of spacetime itself, caused by massive and compact objects.In addition, charged particles experience the Lorentz force when exposed to strong magnetic fields, resulting in deviations from geodesic motion.However, the influence of the magnetic field on the motion of neutral test particles may not be as strong.Similarly, in the presence of strong electric fields, charged particles encounter the Coulomb force, leading to deviations from geodesic motion.This becomes particularly relevant in astrophysics, when studying the motion of charged particles in the intense magnetic and electric fields around neutron stars or in astrophysical jets.
An exhaustive review on the influence of rotation, cosmological constant, and magnetic field on the motion of neutral and charged particles in the accretion disks around Kerr black hole candidates is provided in Ref. [53].Realistic astrophysical scenarios around compact objects such as the properties of thin and thick, neutral and charged accretion disks, quasiperiodic oscillations, and relativistic jets, are summarized there.
In this paper, a solution for a so-called double-charged dilatonic BH and particular solutions for null and timelike geodesics are considered.The BH solution under consideration takes place in a gravitational model with two ("neutral," real-valued) scalar fields and two Abelian gauge fields.(For more general dilatonic BH solutions with n electric color charges in the gravitational model with n scalar fields and n Abelian gauge fields, see Ref. [54] and references therein) .
In fact, our solution is S-dual to the dyonic-like solution from Ref. [24].The metric of the solution may be obtained directly from that of Ref. [24] by a replacement: ⃗ λ 2 → − ⃗ λ 2 , where ⃗ λ 2 is the two-dimensional dilatonic coupling vector corresponding to a second gauge field.Thus, here, we have an Abelian gauge group (U(1)) 1 × (U(1)) 2 and our double-charged black hole carries a pair of charges (Q 1 , Q 2 ).The first electric (color) charge Q 1 (of color "number 1") corresponds to the first subgroup (U(1)) 1 , while the second electric (color) charge Q 2 (of color "number 2") corresponds to the second subgroup (U(1)) 2 .As to "physical relevance" of the model and our solution under consideration, they are on the same footing with those from Ref. [24] and numerous other dyon and dyon-like dilatonic BH solutions discussed in physical journals [54].
It is worth noting that there is considerable interest in spherically symmetric solutions, as evidenced by studies such as [55][56][57][58] and others.These solutions appear in gravitational models involving scalar fields and antisymmetric forms.In our analyses, we follow the methodology of [59], where the motion of neutral test particles is studied in the field of a Reissner-Nordström BH.
The paper is organized as follows: in Sect. 2 we present charged BH solutions characterized by two scalar (dilaton) fields and two Abelian vector fields, discussing their key features.Section 3 is dedicated to the examination of geodesics followed by neutral test particles and photons.Here, we derive expressions for energy, angular momentum, and effective potential, scrutinize their behavior, and analyze the stability of circular geodesics.Finally, in Sect.4, we summarize our findings and discuss potential avenues for future research.

Charged black hole solution
The action of a model containing two scalar fields, 2form and dilatonic coupling vectors is given by where is the vector of scalar fields belonging to and R[g] is the Ricci scalar.Here, and in what follows, we set G = c = 1 (where c is the speed of light in vacuum.) We consider a so-called double-charged black hole solution to the field equations corresponding to the action (2.1) which is defined on the (oriented) manifold and has the following form with the 2-form defined by where Q 1 and Q 2 are the (color) electric charges, µ > 0 is the extremality parameter, dΩ2 = dθ 2 + sin 2 θ dφ 2 is the canonical metric on the unit sphere S 2 (0 < θ < π, 0 < φ < 2π), τ = sin θ dθ ∧ dφ is the standard volume form on S 2 and the moduli function is adopted as with P > 0 obeying or equivalently (2.9) All the remaining parameters of the solution are defined as follows 2 ) i = 1, 2 and Here, the following additional restrictions on dilatonic coupling vectors are imposed We note that Because of relations (2.14) and (2.15), the Q 2 s are well defined.Note that the restrictions (2.14) imply relations ⃗ λ s ̸ = ⃗ 0, s = 1, 2, and (2.2).
Indeed, in this case we have the sum of two nonnegative terms in (2.12): (2.17) For non-collinear vectors ⃗ λ 1 and ⃗ λ 2 we get 0 < a < 2 , and for collinear ones we get a = 2 .This solution may be verified just by a straightforward substitution into the equations of motion.It may also be extracted as a special dilatonic BH solution from Ref. [54].
In addition, in Ref. [24], the definition of gravitational mass was obtained in relation to µ and P parameters: (2.18) The following relations can be found/verified from Eqs. (2.10)-(2.13) and The calculation of the scalar curvature for the metric ds 2 = g µν dx µ dx ν in (2.4) yields [24] For the Schwarzschild (a = 0) and Reissner-Nordström (a = 2) solutions, one immediately obtains R[g] = 0.

Geodesic motion
Geodesics are a fundamental tool for understanding the motion of test particles in the field of dilatonic black holes and can provide valuable insights into the nature of these objects and their effects on surrounding spacetime.Geodesics can be derived from the Lagrangian using the Euler-Lagrange equations: where ẋα = dx α /dτ = u α is the 4-velocity vector, e.g. of a test particle, moving along the curve x α (τ), τ is the proper time for massive particles moving along timelike geodesics and affine parameter in the case of null geodesics, respectively, and α = 0, 1, 2, 3.The generalized momentum p α = g αβ (x) ẋβ for Lagrangian (3.1) is normalized as : where k = −1, 0, 1 for spacelike, null, and timelike geodesics, correspondingly.For general k, the quantity (−k/2) is the energy integral of motion for the Lagrangian (3.1).
For simplicity, we do not consider spacelike geodesics in the paper.We consider null or timelike geodesics in the equatorial plane (θ = π/2).In this case, the Lagrangian for the metric given by Eq. (2.4) reads For cyclic coordinates t and φ , one obtains the integrals of motion associated for k = 1 with the total energy E = Ẽm and angular momentum L = Lm of a test (neutral point-like) particle of mass m, which are the constants of motion.In order to find the equation for the R-coordinate for Ṙ ̸ = 0, one has to use Eq.(3.3), which has the following form for the line element from Eq. (2.4) Eq. (3.7) can be presented in terms of the effective potential: which is explicitly given by It may be readily verified that the Lagrange equation for the radial coordinate R is equivalent to the following one For Ṙ ̸ = 0, it can be obtained just by differentiating Eq. (3.8) with respect to parameter τ and then dividing the result by Ṙ.In case Ṙ = 0, the radial equation reads It does not follow from Eq. (3.8) and should be considered separately.Now we focus on timelike geodesics (k = 1).The behavior of the effective potential as a function of R/µ for the fixed value of Q/µ = 0.6 (which is roughly P/µ = 0.0863) and for different values of L/(µm) is illustrated in Fig. 1.The cases of a = 0, a = 1, and a = 2 formally correspond to the Schwarzschild, Sen [60], and Reissner-Nordström solutions, respectively.It is clearly seen that the maximum of the effective potential for a = 0 exceeds the maximums for other cases of a.
In Fig. 2, the effective potential at fixed L ISCO is illustrated as a function of the dimensionless radial coordinate R/(2µ) for different a = 0, 0.5, 1, 1.5, 2. Dots show the inflection points.As one may see, in terms of R and µ these points differ from the standard Schwarzschild and Reissner-Nordström black holes.
Moreover, a three-dimensional plot of the effective potential with fixed value of Q/µ = 0.6 is shown in Fig. 3 as a function of R/µ and L/(µm).
In Table 1 we present the radii of ISCO, corresponding orbital angular momentum, energy, and effective potential of test particles for fixed Q/µ = 0.6 as a function of parameter a.The graphical representation of these values (points) is shown in Figs. 2 -8.

Circular geodesics
Here, we contemplate circular motions, which are characterized by condition: Ṙ = 0, so V = E/m.The derivative of the effective potential with respect to R is given by: The radial equation (3.11) implies which is substituted to Eq. (3.9) to obtain: . (3.14) Figures. 4 and 5 show E/m and L/(µm) as functions of R/µ for fixed Q/µ = 0.6.Here, the distinctions among all the curves with various a will be larger for larger values of Q/µ.  1 From Eqs. (3.13) and (3.14) it can be seen that for timelike geodesics, the motion is possible only for R(R − 3µ) + P(R(1 − a) + µ(2a − 3)) > 0, with limiting radius: where R γ+ = R 0 , in fact, is the radius of the photon sphere and R γ− is not physical since it is less than 2µ.The dependence of the normalized radius of photon sphere R 0 /µ on Q/µ is depicted in Fig. 6.As one can see, the larger Q/µ, the larger R 0 /µ, at least in this representation of parameters.

Innermost stable circular orbits
Innermost stable circular orbit (ISCO) is an important quantity in the study of geodesics in the field of black  holes.For example, it plays a crucial role in the physics of accretion disks, as it defines the inner radius of the disk.Any massive particle that goes beyond the ISCO will fall onto a black hole, so the disk is considered to have a radius larger than the ISCO to maintain its structure.The location of the ISCO and the properties of the accretion disk (such as its size, temperature, and brightness) can be used to study the features of black holes and their surrounding environment.Using the condition For a = 0, which corresponds to the Schwarzschild black hole, the ISCO radius is given by where In Fig. 7, the radii of ISCOs R ISCO /µ are presented versus Q/µ.As one can see, in the representation for vanishing Q/µ, all radii reduce to 6, and for increasing Q/µ, apart from a=0 case, R ISCO /µ grows nonlinearly.
For a = 1/2, the ISCO radius is where p = P/2µ.For a = 1, which is formally identical to the Sen black hole case, the result is straightforward (3.20) For a = 3/2, the ISCO radius is For a = 2, which corresponds to the Reissner-Nordström black hole case, the ISCO radius is where Using expressions (2.9) and (2.18) allows one to represent Eqs.(3.17) - (3.22) in terms of the gravitational mass M and net charge Q. Correspondingly, in Table 2 we present µ and P in terms of M and Q depending upon a.
Constraints for Q/M can be obtained from Table 2, by requiring µ > 0. Correspondingly, one finds In Ref [24] it was shown that for the a = 2 case under the radial coordinate transformation, R = r RN − P, along with M RN = µ + P; the metric (2.4) coincides with the Reissner-Nordström metric.Taking this into account, one can rewrite Eq. (3.22) in the following form: matches with the equation for r isco in Ref. [59].Here M RN and Q RN are the mass and charge, respectively, of the Reissner-Nordström solution.In addition, the relationship between the charges is given by Q 2 RN = Q 2 /2, whereas the masses are equal M RN = M.
As expected, in the limiting case Q → 0, R ISCO corresponds to the Schwarzschild value 6M.In the cases of 0 < a ≤ 2 the values of R ISCO depend upon the ratio of Q/µ.This behavior is illustrated in Fig. 8.The values of R ISCO at Q = µ, which can also be seen in Fig. 7, are equal to R ISCO /µ = 6.10, 6.22, 6.34, 6.46 for cases a=0.5, 1.0, 1.5, 2.0, respectively.
The energy and the angular momentum in the last stable circular orbit can be found numerically after plug-   In Fig. 9, the normalized photon sphere radius R 0 /(2µ) (left panel) and R ISCO /(2µ) (right panel) are given versus p = P/(2µ).As one can also see in these representations, R 0 /(2µ) and R ISCO /(2µ) increase with increasing p = P/(2µ).As one may notice, the Schwarzschild case with a = 0 possesses R 0 /M = 3, and for other cases, R 0 /M < 3. The Reissner-Nordström case with a = 2 has the smallest radius R 0 /M = 1 for the extreme case Q/M = √ 2. In Fig. 11 (left panel) we present the radii of ISCOs, normalized by gravitational mass M, as a function of net charge over gravitational mass Q/M.As one may see, for increasing charge, the radii of ISCOs decrease as one expects in analogy with the Reissner-Nordström solution [59].In this representation, all physical quantities can be measured and used to distinguish charged black holes with different a, i.e. different coupling constant vectors ⃗ λ 1 and ⃗ λ 2 .
It is interesting to note that a similar result was obtained in Fig. 3 (left panel) of Ref. [61] for a neutral test particle in the Sen spacetime.It is shown that for an increasing charge over a gravitational mass ratio, the radius r ISCO , normalized by gravitational mass M, decreases.In order to compare this result with our findings one should find the relationship between the net charge Q of the present paper and the one in Ref. [61], which we denote as Q Sen .To do so, one must take a close look at the line element Eq. ( 1) in Ref. [61] and compare it with the one considered here when a = 1.The two solutions, apart from notations, physical origin and interpretation, are identical.Thus, by comparing the two solutions one finds that Q = 2Q Sen and the value Q/M = 1.4 in Fig. 11 left panel for a = 1 is equal to Q Sen /M = 0.7 in Fig. 3 (left panel) of Ref. [61].
Furthermore, one of the quantities which is of great interest is the efficiency of converting matter into radiation (see for details page 662 of Ref. [62]) In Fig. 11 (right panel) the efficiency is shown versus Q/M.As we expected, the efficiency for different values of a = 0.0, 0.5, 1.0, 1.5, 2.0 is always larger than the a = 0 case.
In general, the circular orbits are allowed in the region R 0 < R < R ISCO ; however, all those orbits are unstable.For stability of circular orbits, the following condition must be fulfilled: The behavior of the second derivative of the effective potential is depicted in Figs. 12 and 13.From these plots, it can be seen where orbits become stable.The location of ISCO is determined by finding the radius at which the second derivative of the effective potential changes sign from negative to positive.

Conclusion
In this work, we have considered the solution for a doublecharged dilatonic black hole.We investigated circular geodesics of massive neutral test particles and photons, adopting various values of parameters a and P (the latter being directly linked to the net charge Q).We calculated the energy, angular momentum, and the radius of the ISCO for test particles, expressed in terms of a, P, and µ.We conducted a detailed analysis of their behavior across different scenarios, specifically for cases where a=0, 1/2, 1, 3/2, 2, corresponding to distinct configurations of the double-charged dilatonic black holes.
Our analysis predominantly relies on investigating the behavior of an effective potential that determines the position and stability characteristics of circular orbits.The stability of circular geodesics was assessed by examining the sign of the second derivative of the effective potential concerning the radial coordinate.Notably, stable orbits for neutral test particles are observed only within the range from the R ISCO and extending to infinity.
The ISCO radius and photon sphere radius were computed for specific values of a.It turned out that in this parameter representation of the line element, it was observed that for larger a with an increasing ratio of Q/µ, the ISCO and photon sphere radii also increase.The minimum values for the photon sphere radius and the ISCO radius, R 0 /µ = 3 and R ISCO /µ = 6, respectively, are attained in the Schwarzschild limiting case (a=0).
However, in the representation of the net charge over gravitational mass, the situation is utterly opposite.In the Schwarzschild case (a = 0), the largest values for the photon sphere and ISCO radii are obtained as 3M and 6M, respectively.For the Reissner-Nordström case (a = 2), one obtains M and 3M, respectively.Referring to Table 2, for a = 2 and Q/M = √ 2 the parameter P will be equal to M. Correspondingly, using the coordinate transformation r RN = R+P and the relation between charges Q RN = Q/ √ 2, the photon sphere and ISCO radii will be equal to 2M and 4M, respectively.These results align with those reported in Ref. [59].
The efficiency of converting matter into radiation is expected to be larger for the Reissner-Nordström case η = 8.14% and smaller for the Schwarzschild case η = 5.72% [63].All other configurations fall within the range defined by these two cases.This peculiarity can be used to distinguish ordinary or astrophysical black holes from the double-charged dilatonic black holes.
It would be interesting to extend the analyses of the paper in future studies related to the quasinormal modes at the final moments of black hole mergers in binaries, quasiperiodic oscillations in the X-ray systems, radiative flux, and spectral luminosity of accretion disks around astrophysical black holes.Another possible (and rather natural) generalization of our setup may be in considering the motion of a massive point-like particle carrying the electric (color) charge doublet (q 1 , q 2 ), corresponding to our gauge group (U(1)) 1 × (U(1)) 2 .Additionally, the inclusion of the source rotation and background test magnetic field will be fascinating, considering realistic astrophysical scenarios in analogy to Ref. [53].This and other intriguing problems may be considered in our future studies.

Fig. 1 :
Fig.1: The effective potential of test particles for L/(µm)=0, 3, 4, 5, 6 with different a=0, 0.5, 1, 1.5, 2 as a function of the dimensionless/normalized radial coordinate R/(2µ).As one may notice that for decreasing L, the minimum of the effective potential shifts from right to left

Fig. 3 :
Fig.3: Three-dimensional plot of the effective potential V as a function of L/(µm) and R/µ according to Eq. (3.9).The black solid curve indicates V when L/(µm) is defined from Eq. (3.13) for circular geodesics.Dot shows the radius of ISCO with corresponding L/(µm) and V .Here we only consider the case with a = 1 and Q/µ = 0.6

Fig. 8 :Fig. 9 :
Fig. 8: The radius of ISCO, normalized by 2µ, as a function of Q/(2µ) for different a=0.5, 1, 1.5, 2. Numbers near the points imply the values of energy E/m and angular momentum L/(µm) (underlined numbers) corresponding to a specific ISCO

Fig. 10 :
Fig. 10: The location of photon sphere radii, normalized by gravitational mass M, as a function of net charge over gravitational mass Q/M. a

Fig. 10
Fig. 10 the photon sphere radius, normalized by gravitational mass M, is shown versus Q/M for different a.As one may notice, the Schwarzschild case with a = 0 possesses R 0 /M = 3, and for other cases, R 0 /M < 3. The

Fig. 11 :
Fig. 11: Left: The radii of ISCOs, normalized by gravitational mass M, as a function of Q/M for different values of a. Right: The efficiency as a function of Q/M for different values of a.

2 V/∂R 2 Fig. 12 :
Fig. 12: Representation of the second derivative of the effective potential as a function of R/µ for a = 0.

Table 2 :
Values of parameters µ and P in terms of M and Q depending on a.Note that for vanishing Q one obtains P = 0 and µ = M.At the same time, for vanishing a, one recovers the Schwarzschild metric, and µ = M and