Particle motion around generic black holes coupled to non-linear electrodynamics

We study spherically symmetric magnetically charged generic black hole solutions of general relativity coupled to non-linear electrodynamics (NED). For characteristic values of the generic spacetime parameters we give the position of horizons in dependence on the charge parameter, demonstrating separation of the black hole and no-horizon solutions, and possibility of existence of solutions containing three horizons. We show that null, weak and strong energy conditions are violated when the outer horizon is approaching the center. We study effective potentials for photons and massive test particles and location of circular photon orbits (CPO) and innermost stable circular orbit (ISCO). We show that the unstable photon orbit can become stable, leading to the possibility of photon capture which affects on silhouette of the central object. The position of ISCO approaches the horizon with increasing charge parameter q and the energy at ISCO decreases with increasing charge parameter. We investigate this phenomenon and summarize for a variety of the generic spacetime parameters the upper estimate on the spin parameter of the Kerr black which can be mimicked by the generic charged black hole solutions.

We study spherically symmetric magnetically charged generic black hole solutions of general relativity coupled to non-linear electrodynamics (NED). For characteristic values of the generic spacetime parameters we give the position of horizons in dependence on the charge parameter, demonstrating separation of the black hole and no-horizon solutions, and possibility of existence of solutions containing three horizons. We show that null, weak and strong energy conditions are violated when the outer horizon is approaching the center. We study effective potentials for photons and massive test particles and location of circular photon orbits (CPO) and innermost stable circular orbit (ISCO). We show that the unstable photon orbit can become stable, leading to the possibility of photon capture which affects on silhouette of the central object. The position of ISCO approaches the horizon with increasing charge parameter q and the energy at ISCO decreases with increasing charge parameter. We investigate this phenomenon and summarize for a variety of the generic spacetime parameters the upper estimate on the spin parameter of the Kerr black which can be mimicked by the generic charged black hole solutions.

I. INTRODUCTION
Physical singularity at the center of the standard vacuum solutions within the Einstein general relativity exists due to the Penrose and Hawking theorems [1]. However, there are several effective methods to avoid the singularity, and get a singularity-free non-vacuum or modified gravity solution describing compact objects in astrophysics. The conformal gravity models [2][3][4][5][6], quantum corrections [7,8], dark energy star -gravastar models [9][10][11][12] can be considered as examples of models providing the singularity-free solutions.
Other interesting way of avoiding the singularity is solution of general relativity coupled to non-linear electrodynamics. Recently the general procedure for obtaining the solutions describing the compact objects with the electric and/or magnetic charges in general relativity coupled to non-linear electrodynamics has been described [13,14]. Such generic solutions could be quite regular black hole (or no-horizon) solutions, if their mass parameter is related exclusively to the NED part [13,14]. In Ref. [14] authors suggested the regular black hole models as the gravitational field of a non-linear electric or * jaroslav.vrba@fpf.slu.cz † ahmadjon@astrin.uz ‡ arman.tursunov@fpf.slu.cz § ahmedov@astrin.uz ¶ zdenek.stuchlik@fpf.slu.cz magnetic monopole. The properties of the rotating analogue of this type black hole have been considered in [15].
In the special cases of generic regular NED black holes and no horizon spacetimes (Bardeen, Hayward) the motion of particles and photons was investigated in [16][17][18][19] Here we study the properties of the singular generic black hole solutions containing also the standard selfinteracting mass parameter, analyze the energy conditions and test particle motion in their environment.
The main purpose of the studying the particle motion around compact objects is to test the model or solution of gravity theories. Particularly, the accretion disc around the astrophysical black holes is considered as a relevant model of some X-ray sources and the observation of inner edge radii of accretion disc can give constraints on parameters of alternative and modified theories of gravity [20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39]. Above presented facts are good motivation for study test particle's motion around black hole described by the general relativistic solution coupled with non-linear electrodynamics.
The paper is organized as follows. In Sect. II we describe generic black hole and consider existence of horizons and at energy conditions. Sect. III is devoted to studying the motion and circular orbits of photons and neutral massive test particles around generic black hole. In Sect. IV the possibility to distinguish the spin of the Kerr black hole from the charge parameter for the generic black hole solution coupled to NED is discussed. Finally, in Sect. V we summarize our obtained results.
Throughout the paper we use a space-like signature (−, +, +, +), a system of units in which G = c = 1 and arXiv:1909.12026v1 [gr-qc] 26 Sep 2019 we restore them when we need to compare our results with observational data. Greek indices run from 0 to 3, Latin indices from 1 to 3.

II. GENERIC BLACK HOLE HORIZONS AND ENERGY CONDITIONS
A. Generic black hole spacetime The generic black hole spacetime metric has the following Schwarzschild-like form in the standard Schwarzschild coordinates We consider the family of solutions where the mass function m(r) is defined as [13,14] where M is the pure gravitational self-interaction mass and the Arnowitt-Deser-Misner (ADM) mass is defined as being the electromagnetically induced gravitational mass, where α > 0 describes the strength of the nonlinear effect, µ > 0 characterizes the degree of nonlinearity, q is the integration constant related to the magnetic charge by relation 3 and ν > 0 is free parameter governing character of its field in concrete solution and we suppose q > 0 as in [40]. The Lagrangian density governing the electromagnetic non-linear self-interaction takes the form [40] where F = F µν F µν . The metric (1) with mass function (2) is used throughout this paper. The magnetic charge can be expressed as [13][14][15] The regular black hole (no physical singularities even at the center r = 0) appears when we assume with condition µ ≥ 3 [40]. Requirement (4) with (2) leads to mass function The electromagnetic field potential of the non-rotating magnetically charged black hole spacetime takes the form The magnetic charge in a black hole could be a remnant of the phase transitions in the early Universe where in accord with Grand-Unification theory magnetic monopoles could be created as defects of the transition [41]. Alternatively see [42][43][44].
In the following we consider the generic singular black hole solutions with the self-gravitating term M = 0, i.e., we do not assume condition (4) for the non-singular black holes in the rest of the paper. We can see immediately that for r → ∞ there is m(r) → M , and the mass parameter M governs the gravitational mass at large distances. We can thus define dimensionless radius r/M → r and dimensionless charge parameter q/M → q. Then we can use dimensionless parameter α/M 2 → α. For using of these dimensionless quantities for simplicity we can put M = 1.

B. Generic black hole horizons
Now we explore existence and behavior of the horizon(s) of generic black holes. The lapse function can be written in the form Four families can be found depending on the spacetime parameters. One can distinguish them according to the number of horizons from zero to three possible horizons which can appear in the considered generic black hole spacetime. The existence and location of horizons comes from a solution of equation g tt = 0, it follows from (7) that we have to solve equation Interesting characteristics is related to the existence of limiting value of charge parameter q L for which all horizons disappear and the object becomes a naked singularity for any bigger value q > q L . We can see that q L depends on all the other black hole parameters, the greater is parameter ν the greater is q L . The parameter  µ has evidently weaker and opposite effect. Influence of the parameter α on horizon is plotted in Fig. 3 and shows how q L increases with parameter α. The outer horizon also increases with increasing of the parameter α (for fixed q), since the negative electromagnetic contribution of mass decreases.
Very interesting result appears when the "innermost" horizon shrinks to zero. At this point the object can lose physical singularity and become non-singular (if µ ≥ 3). Eq. (4) then implies q 3 = α. For that case there may exist two, one or zero horizon. In case of non-singular compact object, for the chosen values of parameters µ and ν, the number of horizons can range from one to three, while one of the horizons is always located at r = 0, as can be seen e.g. in Fig.1. Existence of the event horizon defines the existence of black hole. Therefore the object having zero-radius horizon and two more horizons is a non-singular black hole. The outermost horizon is the event horizon, the middle horizon is a form of the Cauchy horizon. When the two horizons coincide, the object is a non-singular extremal black hole with zero-radius horizon and event horizon.

C. Energy conditions
Here we study energy conditions for spacetime described by (1) and (2), using the Einstein equations is the stressenergy tensor releted to electromagnetic field. It is also obvious (because of the symmetry of spacetime) that P 2 = P 3 . We assume the above assumptions on parameters and black hole solution for the rest of this paper. The energy conditions were tested in [4] or [45][46][47] and even more restrictive causality and unitarity conditions were studied in [48].

Null energy condition
The null energy condition (NEC) can be written as where n µ is a null-like vector. We can rewrite it and show the physical meaning in the form The symmetry of the spacetime (G t t = G r r ) helps again because for i = 1 it is identity. We study only ρ + P 2 ≥ 0 for now and the NEC takes the form First, ρ + P 2 is approaching zero as r tends to infinity. What is happening near the horizon r → r h is a more complex issue, because the horizon position cannot be determined analytically. However, Eq. (11) clearly says if µ > 3, and the outer horizon is very close to center, NEC can be violated. From Figs. 1 and 2, it follows that the outer horizon does not exist when q → q L and sufficiently large ν. In this situation, the NEC is satisfied at the horizon. An example of NEC is shown in figure 4.
Then is easy to find that ρ + P 2 has no zero in the interval r ∈ (r h , ∞). The NEC is fulfilled everywhere in the outer part of spacetime if it is fulfilled at the horizon.

Weak energy condition
The weak energy condition is defined as where u µ is a time-like vector. We can rewrite it and show the physical meaning in the form We examine only ρ ≥ 0 because the second part was studied in the previous section. Again we start with the limiting cases r → r h and r → ∞. ρ is greater than zero at the horizon and tends to zero at infinity.
Our condition ρ ≥ 0 transforms to From Eq. (14) it immediately follow that it is impossible (for our assumptions) to get negative value for ρ. This part (ρ ≥ 0) of the weak energy condition (WEC) is fulfilled everywhere between horizon and infinity. Examples of ρ profile are shown in Fig. 5. The result for the WEC is therefore the same as for the NEC.

Strong energy condition
The strong energy condition is defined in the form where u µ is a time-like vector. It can be transformed to more physically transparent form If we use the spacetime symmetry again (ρ + P 1 = 0 and P 2 = P 3 ), the condition reads The limit r → ∞ is zero. But (17) can be negative at the horizon if µ > 1 and outer horizon is close to the center. The behavior of horizon is seen in Figs. 1 and 2.
The function P 2 has no zero value at interval r ∈ (r h , ∞), it results that the strong energy condition (SEC) is fulfilled in outer spacetime if value at horizon is positive. Fig. 6 presents P 2 for several values of parameters µ and ν.

III. NEUTRAL TEST PARTICLE MOTION AROUND GENERIC BLACK HOLE
Now consider the neutral test particle motion around the generic black hole. We distinguish two situations photons and massive particles. We examine the effective potential for both cases and look at the special orbitsunstable circular photon orbit (UCPO) for photons and ISCO for massive particles. The circular geodesics are determined by the local extrema of the effective potential, dV eff dr = 0. These orbits are stable (unstable) for The local extrema of the effective potential cannot be given in an analytical form due to the general form of the lapse function of the generic spacetimes. Therefore, all the calculations are given numerically and results are presented in graphical form. The effective potential for the radial motion of test particle can be obtained from the normalization condition where δ is 0 or −1 for photons and for massive particles, respectively. Due to the spacetime symmetries we define constants of motion The motion is always fixed to a central plane that can be for simplicity chosen as the equatorial plane (θ = π/2). The radial motion is determined by where τ is the particle's proper time and these constants can be identified as energy and angular momentum of the particle at infinity. The effective potential takes the form where f (r) is the lapse function given by Eq. 7.

A. Orbits of photons
In this subsection we discuss orbits of photons (δ = 0). It is well known that photons do not follow null-geodesics of the generic spacetime as there are coupled to the NED and one had to introduce effective metric [49,50] given by the relation where L = L(F ) is Lagrangian density, F µν is electromagnetic tensor and F = F αβ F αβ is electromagnetic field invariant. Then L F denotes differentiation of Lagrangian density L with respect to invariant F and L F F is double differentiation with respect to F . The optical phenomena in the special cases of Bardeen black hole were for its effective geometry studied in [50,51]. Here we study the photons motion in the effective geometry of the generic NED black holes. The photon's orbits can be determined in relation to the impact parameter b = L E and the effective potential reads [50] V eff = f (r) The radius of circular photon orbits is given by dV eff dr = 0, their impact parameter is given by the effective potential at the given radius.
In Figs. 7 and 8 typical effective potentials for photons are shown and compared to the Schwarzschild case when the photon motion is governed by the spacetime geometry. The Schwarzschild effective potentials were scaled by constant factors for better resolution -the factors are listed in the description of figures. These factors reflect the large differences of the photons redshift in the effective geometry as compared to the spacetime (Schwarzschild) geometry demonstrated for the first time in the special case of the Bardeen spacetimes in [50]. The position of UCPO is a little shifted toward the horizon with rising parameter µ. The parameter ν has a bit different effect. The character of effective potential for small values (ν ≈ 1) is very different from that for greater values (ν ≈ 3). This peculiar behavior for small ν causes that in such specific spacetime UCPO becomes stable circular photon orbit, hence there arises a possible window for trapping the photon in potential well in finite area.
This can be seen in Fig. 9 where the influence of parameters µ and ν on UCPO is plotted. In the left box, we can identify the change of shape of the effective potential by diverging of UCPO. In this point, the shape is switched and there is no circular photon orbit. How UCPO changes with parameters µ, ν and q is seen in Figs. 10 and 11.
The influence of the parameter α is shown in the Fig. 12.

B. Orbits of massive test particles
In this subsection, we examine the behavior of massive test particles (δ = −1). In that case, the effective potential is dependent on angular momentum L. There are three possible types of effective potential represented in Fig. 13. The most important case is the one with an inflex point that determines ISCO. This situation takes place for specific value of angular momentum, L ISCO .
The effect of parameters µ, ν and α on the effective potential is shown in Figs. 14 -16. The role of parameter ν seems to be similar and even stronger than for photons. Effective potentials with ν 1 and q 0.5 are very similar to Schwarzschild one, with a weak effect of parameter µ. The effect of parameter α mimics the effect of L.  The circular geodesics are located at r c , the extrema of the effective potential of a massive test particle that is determined by its specific angular momentum L. The position of r c thus depends on the specific angular momentum L. In other words, at a given r c , the specific energy E and the specific angular momentum L are determined by the local extremes of the effective potential. The radial profiles of L are illustrated in Figs. 17 and 18. Again one can see that for ν 1 and small values of q 0.5 the curves coincide. At the minimum of these curves ISCO is located. The specific energy of the circular geodesic at a given r c is determined by the effective potential extrema taken at r c for given L.
ISCO is very important from the astrophysical point of view. Dependence of ISCO on parameter q for different values of parameters µ and ν is shown in Figs. 19 and 20. We can see that this orbit is attracted to the horizon with growing parameter q, however, this effect is damped with increasing ν.

IV. ASTROPHYSICAL APPLICATION TO THE SMBH SGR A*
The interest in black hole parameters leads us to study accretion processes where ISCO plays a very important role. ISCO is located at r = 6 for non-rotating uncharged black hole. For uncharged rotating Kerr black hole ISCO is related to its mass (we set M = 1) and rotation as [52] (24) where minus is for pro-grade orbits, which are most interesting for us and plus is for retro-grade orbits, and a ≡ J/M is dimensionless spin parameter of Kerr black hole. ISCO for pro-grade orbits is pushed from non-rotating position (r = 6) toward the center. The specific energy of co-rotating ISCO comes when formulas for ISCO positions (24) - (26) are used in the formula for specific energy of the circular orbits [53] However, we can find out the ISCO position in practice in two ways. Either by direct observation [54] or estimation from the ISCO specific energy [55][56][57][58].
From the astronomical observations of the black hole binaries (microquasars) it is possible to determine the inner edge of the accretion disk or nearest stable orbits to   the center, i.e. the ISCO and estimate the spin parameter a. This kind of research has been done for supermassive black hole (SMBH) in the center of Milky way, Sgr A* [59]. But one cannot exclude the possibility that some other black hole parameter can mimic the role of the black hole spin. This black hole parameter can be magnetic or electric charge [60], [61] or alternatively the tidal charge of braneworld black holes [62][63][64]. Here we focus on the possibility to mimic the spin influence by the magnetic charge contained in the generic singular black holes in general relativity combined with NED.  (2) for parameters µ = 2 − 6, ν = 1 − 5 and α = 0.3 guessed from ISCO position.
We estimate here the maximal spin parameter which can be mimicked by an appropriate magnetic charge related to the black hole solutions discussed in our paper. Using previous results presented in Figs. 19, 20 and formulas (24) -(26), one can get how the magnetic parameter can mimic the spin of the black hole in the location of the ISCO. Figs. 21 and 22 show ISCO against charge parameter q in generic black hole spacetime which describes black hole (no naked singularity) and also ISCO for Kerr spacetime against spin parameter a. It is obvious that for some combination of parameters can magnetic charge very effectively mimics the rotation creating the degeneracy in black hole parameters. We see that for increasing parameter ν this phenomenon disappears and the ISCO is close to the Schwarzschild limit. The maximal values of spin parameter a which can be mimicked by maximal magnetic parameter q for variety of parameters µ and ν are summarized in table I. Values are rounded to two decimals.
However, estimates based on the ISCO position have to be complemented by estimates of the specific energy of ISCO orbit (E ISCO ), as this specific energy determines efficiency of the Keplerian accretion and is thus crucial for the black hole parameters in connection with the observational effects [55][56][57][58].
The dependence of ISCO specific energy on charge parameter q for several values of parameters ν and µ is given at figures 23 and 24 while its comparison with Kerr ISCO specific energy is illustrated at figures 25 and 26. One can see that the ISCO specific energy is decreasing with increase of the charge parameter q, similarly to the Kerr spacetime where the specific ISCO energy decreases with increasing spin parameter a. In the NED spacetimes this decrease is significant only for small values of parameter ν ≈ 1. With the increasing parameter ν this phenomenon disappears and the specific energy remains close to the Schwarzschild limit for increasing charge parameter q.

V. CONCLUSION
In the present work we have investigated properties of the generic singular spacetime (1) with mass function (2), found as solution of Einstein equation coupled to   non-linear electrodynamics. We started with existence and position of the horizon. There can exist one, two, three or no horizon. Besides, there always exists the zero-radius horizon for q = α 1/3 . The last possibility indicates naked singularity. We have found out that for any specific values of the parameters α, µ and ν, limiting value of magnetic parameter q L exists which separates black hole and naked singularity solution. This can be seen in Figs. 1 and 2. It can be seen from Fig. 3 that increasing parameter α leads to increasing location of the outer horizon. It is necessary to point out that when "innermost" horizon reaches the zero, the Eq. (4) is satisfied and we have at this point singularity free solution if µ ≥ 3. We can find one, two or no horizon in that case.
Then we have tested null, weak and strong energy conditions for the generic spacetime. We have found that, with few exceptions, especially when the outer horizon isn't shifted very close to the center and ν isn't small (ν 1), all conditions, NEC, WEC and SEC are met, as can be seen in Figs. 4 -6. In the section III, the behavior of uncharged test particles has been investigated. It is divided into two parts, the behavior of the photons -electromagnetic waves and massive test particles.
We have constructed effective potential to describe the main characteristics of the particle motion, Figs. 7 and 8. Interesting is that for ν 1, the effective potential can be reversed and unstable circular photon orbit becomes stable, leading to possibility of trapped areas and affect in the silhouette of the object [65]. We have localized the position of unstable/stable circular photon orbit. The influence of the spacetime parameters on the photon motion is very complex and is depicted in Figs. 9 -12.
Then we have investigated the effective potentials for massive particles, Figs. 13 -16. In this case the inner-  most stable circular orbit (ISCO) is very important as it determines efficiency of Keplerian accretion disks [66]. It was expected that with growing parameter q the ISCO approaching the horizon -as in Maxwellian case -but the main role here has again parameter ν, as shown in Figs. 19 and 20. The greater value of ν is the lower shift of the ISCO occur with growing q and ISCO stays close to its Schwarzschild position r = 6.
The shifting of ISCO toward the horizon with increasing the parameter q is similar to the effect of the rotation in the Kerr spacetime. If we admit the existence of magnetic charge in some objects, for example Sgr A* and assume the spacetime described by the generic singular solutions of general relativity combined with NED, then in same choices of its spacetime parameters and charge, the ISCO postiion (and energy) could be equal to those corresponding to the Kerr spacetime. In this sense the NED charge can mimic the spin parameter of Kerr spacetime. This degeneracy can arise if we estimate spin parameter directly from observation of ISCO position using e.g. GRAVITY or EHT, but also from indirect estimate from ISCO specific energy. The comparison of the shift of ISCO in the generic NED spacetimes with those related to the Kerr spacetime is shown in Figs. 21 and 22, and change of the specific energy in Figs. 25 and 26. In Table I we summarize values of the Kerr spin parameter a which can be mimicked by the magnetic charge of the generic NED black hole solution. The significant effect of mimicking of spin by the magnetic charge disappears with increasing parameter ν and for values ν 3 it is impossible completely.
Our results show that the ISCO mimic effect works for both ISCO position and energy in the charge and spin interval q, a ∈ (0.45 − 0.6).
Further investigation will focus on the motion of electrically charged particles in the field of generic black hole.