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

Ergosphere and shadow of a rotating regular black hole

Sushant G. Ghosh [email protected], [email protected] Centre for Theoretical Physics, Jamia Millia Islamia, New Delhi 110025, India Multidisciplinary Centre for Advanced Research and Studies (MCARS),
Jamia Millia Islamia, New Delhi 110025, India
Astrophysics and Cosmology Research Unit, School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Private Bag X54001, Durban 4000, South Africa
   Muhammed Amir [email protected] Astrophysics and Cosmology Research Unit, School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Private Bag X54001, Durban 4000, South Africa    Sunil D. Maharaj [email protected] Astrophysics and Cosmology Research Unit, School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Private Bag X54001, Durban 4000, South Africa
Abstract

The spacetime singularities in classical general relativity predicted by the celebrated singularity theorems are formed at the end of gravitational collapse. Quantum gravity is the expected theory to resolve the singularity problem, but we are now far from it. Therefore attention has shifted to models of regular black holes free from the singularities. A spherically symmetric regular toy model was obtained by Dymnikova (1992) which we demonstrate as an exact solution of Einstein’s field equations coupled to nonlinear electrodynamics for a Lagrangian with parameter bb related to magnetic charge. We construct rotating counterpart of this solution which encompasses the Kerr black hole as a special case when charge is switched off (b=0b=0). Event Horizon Telescope has released the first image of supermassive black hole M87, revealing the structure near black hole horizon. The rotating regular black hole’s shadow may be useful to determine strong field regime. We investigate ergosphere and black hole shadow of rotating regular black hole to infer that their sizes are sensitive to charge bb and have a richer chaotic structure. In particular, rotating regular black hole possess larger size, but less distorted shadows when compared with Kerr black holes. We find one to one correspondence between ergosphere and shadow of the black hole.

I Introduction

The elegant theorems of Hawking and Penrose imply that spacetime singularities are pervasive features of general relativity, so that the theory itself predicts its own failure to describe the physics of these extreme situations. The spacetime singularities that arise in gravitational collapse are always hidden inside black holes, which is the essence of the weak cosmic censorship conjecture, put forward 40 years ago by Penrose Penrose:1969pc . This signature is still one of the most important open questions in general relativity. We are far away from any robust and reliable quantum theory of gravity capable of resolving the singularities in the interior of black holes. Hence, there is significant attention towards models of black hole solutions without the central singularity. The earliest idea of regular models dates back to the pioneering work of Sakharov Sakharov:1966 and Gliner Gliner:1966 where they proposed that the spacetime in the highly dense central region of a black hole should be de Sitter-like at r0r\simeq 0. Thus spacetime filled with vacuum could provide a proper discrimination at the final stage of gravitational collapse, replacing future singularity Gliner:1966 . The prototype of these regular black holes is the Bardeen metric Bardeen:1968 , which can be formally obtained by coupling Einstein’s gravity to a nonlinear electrodynamic field AyonBeato:2000zs , thought to be an alteration of the Reissner-Nordström solution. There has been enormous attentions to obtain regular black holes Dymnikova:2004zc ; Bronnikov:2000vy ; Shankaranarayanan:2003qm ; Hayward:2005gi ; Ansoldi:2008jw ; Culetu:2014lca ; Balart:2014jia ; Balart:2014cga ; Xiang:2013sza ; most of these regular black holes were driven by the Bardeen idea Bardeen:1968 . However, these are only nonrotating black holes, which can be hardly tested by observations, as black hole spin is an important parameter in astrophysical processes.

Further, a study of rotating regular solutions is important as the astronomical observations predict that the astrophysical black hole might be a Kerr black hole Bambi:2011mj . However, the real nature of these astrophysical black holes is still to be verified Bambi:2011mj . Analyzing the iron kαk{\alpha} lines and the continuum fitting method are possible procedures, which can probe the geometry of spacetime around astrophysical black holes candidates Bambi:2014nta . The rotating solutions of the Bardeen metrics have been obtained in Ref. Bambi:2011mj . There exists a number of rotating regular black holes Bambi:2013ufa ; Neves:2014aba ; Toshmatov:2014nya ; Ghosh:2014hea ; Toshmatov:2017zpr which were discovered with the help of the Newman-Janis algorithm Newman:1965tw , and also by using other techniques Azreg-Ainou:2014aqa ; Azreg-Ainou:2014nra ; Azreg-Ainou:2014pra ; Ghosh:2014pba . The authors in Refs. Ghosh:2014mea ; Amir:2015pja ; Ghosh:2015pra ; Amir:2016nti ; Ahmed:2018fge demonstrated that the regular black holes can be considered as the particle accelerator. The quasinormal modes of test fields around the regular black holes has been discussed explicitly in Toshmatov:2015wga . An interesting study on electromagnetic perturbation of the regular black holes has been discussed thoroughly in Toshmatov:2018tyo ; Toshmatov:2018ell ; Toshmatov:2019gxg .

The black hole shadow is a dark zone in the sky and the shadow is a useful tool for measuring black hole parameters because its shape and size carry impression of the geometry surrounding the black hole. Observing the black hole shadow is a possible method to determine the spin of the black hole, and this subject is popular nowadays. This triggered theoretical investigation of the black holes for wide variety of spacetimes Chandrasekhar92 ; Falcke:1999pj ; Takahashi04 ; Takahashi:2005hy ; Bambi:2008jg ; Hioki:2009na ; Bambi:2010hf ; Amarilla:2010zq ; Amarilla:2011fxx ; Abdujabbarov:2012bnn ; Yumoto:2012kz ; Li:2013jra ; Amarilla:2013sj ; Atamurotov:2013dpa ; Atamurotov:2013sca ; Wei:2013kza ; Grenzebach:2014fha ; Papnoi:2014 ; Abdujabbarov:2015xqa ; Johannsen:2015qca ; Goddi:2016jrs ; Abdujabbarov:2016hnw ; Younsi:2016azx ; Amir:2016cen ; Grenzebach:2016 ; Zakharov ; Singh:2017vfr ; Amir:2017slq ; Kumar:2017vuh . Apart from black holes, the theoretical investigation of shadow of wormholes has also been discussed in Nedkova:2013msa ; Bambi:2013nla ; AzregAinou:2014dwa ; Ohgami:2015nra ; Ohgami:2016iqm ; Abdujabbarov:2016efm ; Shaikh:2018kfv ; Amir:2018szm ; Amir:2018pcu . Hioki and Maeda Hioki:2009na discuss a simple relation between the shape of shadow, inclination angle and spin parameter for the Kerr black hole. A coordinate-independent method was proposed to distinguish the shape of black hole shadow from other rotating black holes Abdujabbarov:2015xqa . From the observational point of view, the study of shadow is significant because of the Event Horizon Telescope111https://eventhorizontelescope.org. It was setup with an aim to image the event horizon of supermassive black holes Sgr A and M87. They have successfully released the first image of M87 which is in accordance with the Kerr black hole predicted by general relativity Akiyama:2019cqa ; Akiyama:2019fyp ; Akiyama:2019eap . The photon ring and shadow have been observed which open a new window to test general relativity and modified theories of gravity in strong field regime.

The purpose of this paper is to construct the rotating counterpart or Kerr-like regular black hole from the spherically symmetric regular black hole proposed by Dymnikova Dymnikova:1992ux . The black hole is singularity free or regular black hole with an additional parameter bb, and for definiteness, we can call it the rotating regular black hole. We investigate several properties including the shadow of this black hole to explicitly bring out the effect of parameter bb, and compare with the Kerr black hole. It turns out that the parameter bb leaves a significant imprint on the black hole shadow, and also affects the other properties.

The paper is organized as follows. In Sec. II, we show that the spherically symmetric regular black hole obtained by Dymnikova Dymnikova:1992ux is an exact solution of general relativity coupled to nonlinear electrodynamics. We introduce a line element of the rotating regular black hole spacetime and discuss associated sources with it in Sec. III. The weak energy conditions is the subject of Sec. IV. In Sec. V, we discuss some basic properties of the rotating regular black hole. Shadow of rotating regular black hole is the subject of Sec. VI where we derive analytical formulae for the shadow and we conclude in the Sec. VII. We consider the signature convention (,+,+,+-,+,+,+) for the spacetime metric and use the geometric unit G=c=1G=c=1 throughout the paper.

II Nonlinear electrodynamics and exact regular black hole solution

The spherically symmetric metric proposed by Dymnikova Dymnikova:1992ux represents a regular black hole with de Sitter core instead of a singularity. The stress-energy tensor responsible for the geometry describes a smooth transition from the standard vacuum state at infinity to isotropic vacuum state through anisotropic vacuum state in the intermediate region. The spacetime metric of the Dymnikova solution reads

ds2\displaystyle ds^{2} =\displaystyle= [12M(1er3/b3)r]dt2+[12M(1er3/b3)r]1dr2+r2dΩ2,\displaystyle-\left[1-\frac{2M(1-e^{-r^{3}/b^{3}})}{r}\right]dt^{2}+\left[1-\frac{2M(1-e^{-r^{3}/b^{3}})}{r}\right]^{-1}dr^{2}+r^{2}d\Omega^{2}, (1)

where dΩ2=dθ2+sin2θdϕ2d\Omega^{2}=d\theta^{2}+\sin^{2}\theta d\phi^{2}, MM is the black hole mass, and the parameter bb is given by b3=r02Mb^{3}=r_{0}^{2}M. The corresponding energy-momentum tensor of metric (1) is Dymnikova:1992ux

T00\displaystyle T_{0}^{0} =\displaystyle= T11=3r02er3/b3,\displaystyle T_{1}^{1}=\frac{3}{r_{0}^{2}}e^{-r^{3}/b^{3}},
T22\displaystyle T_{2}^{2} =\displaystyle= T33=3r02(13r32b3)er3/b3.\displaystyle T_{3}^{3}=\frac{3}{r_{0}^{2}}\left(1-\frac{3r^{3}}{2b^{3}}\right)e^{-r^{3}/b^{3}}. (2)

The metric (1) admits two horizons, i.e., the event horizon and the Cauchy horizon, but there is no singularity. We know that the horizons are zeros of grr=0g^{rr}=0, which are given Dymnikova:1992ux by

r+M[1𝒪(exp(M2r02))],r_{+}\approx M\left[1-\mathcal{O}\left(\exp(-\frac{M^{2}}{r_{0}^{2}})\right)\right],
rr0[1𝒪(exp(r04M))].r_{-}\approx r_{0}\left[1-\mathcal{O}\left(\exp(-\frac{r_{0}}{4M})\right)\right].

The metric (1) is regular, everywhere including at r=0r=0, which is evident from the behavior of the curvature invariants. Thus the metric (1) that for large rr coincides with the Schwarzschild solution, for small rr behaves like the de Sitter spacetime and describes a spherically symmetric regular black hole. In what follows, we show the metric (1) as an exact solution of general relativity coupled to an appropriate nonlinear electrodynamics.

We start with action of general relativity coupled to the nonlinear electrodynamics which is given by

S=116πd4xg[R(F)],S=\frac{1}{16\pi}\int d^{4}x\sqrt{-g}\left[R-\mathcal{L}(F)\right], (3)

where RR is the Ricci scalar and (F)\mathcal{L}(F) is an arbitrary function of F=14FμνFμνF=\frac{1}{4}F_{\mu\nu}F^{\mu\nu} with Fμν=2[μAν]F_{\mu\nu}=2\nabla_{[\mu}A_{\nu]}, and AμA_{\mu} denotes the electromagnetic potential. The Einstein field equations are derived from the action (3) which read simply

Gμν\displaystyle G_{\mu\nu} =\displaystyle= 2(FFμαFναgμν),\displaystyle 2\left(\frac{\partial\mathcal{L}}{\partial F}F_{\mu\alpha}F^{\alpha}_{\nu}-g_{\mu\nu}\mathcal{L}\right), (5)
μ(FFμν)=0.\displaystyle\nabla_{\mu}\left(\frac{\partial\mathcal{L}}{\partial F}F^{\mu\nu}\right)=0.

We further consider the spherically symmetric line element of the following form

ds2\displaystyle ds^{2} =\displaystyle= [12M(r)r]dt2+[12M(r)r]1dr2+r2dΩ2.\displaystyle-\left[1-\frac{2M(r)}{r}\right]dt^{2}+\left[1-\frac{2M(r)}{r}\right]^{-1}dr^{2}+r^{2}d\Omega^{2}. (6)

For the spherically symmetric spacetime the only nonzero components of FμνF_{\mu\nu} appropriate to magnetic charge are FθϕF_{\theta\phi} and FϕθF_{\phi\theta} with Fθϕ=FϕθF_{\theta\phi}=-F_{\phi\theta} such that Fθϕ=2bsinθF_{\theta\phi}=2b\sin\theta Fernando:2016ksb . Now it is straight forward to compute

F=2b2r4=2B2,F=\frac{2b^{2}}{r^{4}}=2B^{2}, (7)

where BB is the absolute value of magnetic field. The regular solution Dymnikova:1992ux in which we are much interested comes from particular nonlinear electrodynamics source

(F)=3sb2exp[(2b2F)3/4],\mathcal{L}(F)=\frac{3}{sb^{2}}\exp\left[-\left(\frac{2}{b^{2}F}\right)^{3/4}\right], (8)

where the parameter ss is given by s=|b|/Ms={|b|}/{M} with bb and MM are arbitrary constant can be related to the magnetic charge and the black hole mass, respectively. The (t,t)(t,t) component of the Einstein’s field equations (5) gives

M(r)=r2(F).M^{\prime}(r)=r^{2}\mathcal{L}(F). (9)

On substituting (8) into (9) and then integrating it results in

rM(r)𝑑r=MM(r)=Mer3/b3,\int_{r}^{\infty}M^{\prime}(r)dr=M-M(r)=Me^{-r^{3}/b^{3}}, (10)

where the integration constant is purposely chosen as MM, thus the above equation yields

M(r)=M(1er3/b3).M(r)=M(1-e^{-r^{3}/b^{3}}). (11)

When we substitute (11) in (6) eventually we get the spherically symmetric Dymnikova regular spacetime (1). Thus, we can say that the Dymnikova spacetime (1) can be obtained exactly with source from a particular nonlinear electrodynamics where the source is given by the Lagrangian (8).

III Rotating regular black holes

In this section, we aim to study the rotating counterpart of the Dymnikova spacetime which is a generalization of the Kerr black hole. We start by constructing a rotating counterpart of the regular spacetime (1) which in the Boyer-Lindquist coordinates reads

ds2\displaystyle ds^{2} =\displaystyle= (12M(r)rΣ)dt24aM(r)rsin2θΣdtdϕ+ΣΔdr2+Σdθ2\displaystyle-\left(1-\frac{2M(r)r}{\Sigma}\right)dt^{2}-\frac{4aM(r)r\sin^{2}\theta}{\Sigma}dtd\phi+\frac{\Sigma}{\Delta}dr^{2}+\Sigma d\theta^{2} (12)
+(r2+a2+2a2M(r)rsin2θΣ)sin2θdϕ2,\displaystyle+\left(r^{2}+a^{2}+\frac{2a^{2}M(r)r\sin^{2}\theta}{\Sigma}\right)\sin^{2}\theta d\phi^{2},

where

Σ=r2+a2cos2θ,Δ=r2+a22rM(r).\Sigma=r^{2}+a^{2}\cos^{2}\theta,\quad\Delta=r^{2}+a^{2}-2rM(r).

Here aa represents the rotation parameter and M(r)M(r) is the mass function given in (11). Henceforth, we use the term rotating regular black hole for the spacetime metric (12). It is noticeable that in the absence of charge (b=0b=0), we obtain the Kerr spacetime. The spacetime (12) has independence on tt and ϕ\phi coordinates that means tt and ϕ\phi coordinates are cyclic coordinates in the Lagrangian which corresponding to two Killing vectors, namely, the time translation Killing vector χa\chi^{a}, and the azimuthal Killing vector ζa\zeta^{a}.

Now we are going to check the validity of the rotating solution by computing the source nonlinear electrodynamics equations. The action and corresponding field equations are explicitly expressed in the previous section. The vector potential in case of spherically symmetric spacetime is given by Aμ=bcosθδμϕA_{\mu}=-b\cos\theta\delta^{\phi}_{\mu}. The vector potential for the rotating spacetime gets modify Toshmatov:2017zpr and given as follows

Aμ=[bacosθΣ,0,0,b(r2+a2)cosθΣ].A_{\mu}=\left[\frac{ba\cos\theta}{\Sigma},0,0,-\frac{b(r^{2}+a^{2})\cos\theta}{\Sigma}\right]. (13)

It is easy to compute the field strength (F)(F) for the rotating spacetime, which turns out to be in following form

F=2b2[(r2a2cos2θ)24a2r2cos2θ]Σ4.F=\frac{2b^{2}[(r^{2}-a^{2}\cos^{2}\theta)^{2}-4a^{2}r^{2}\cos^{2}\theta]}{\Sigma^{4}}. (14)

We can immediately recover F=2b2/r4F=2b^{2}/r^{4} while substitute a=0a=0 in (14). In order to compute the source for the rotating spacetime, we solve the Einstein field equations, Gμν=TμνG_{\mu\nu}=T_{\mu\nu} with respect to \mathcal{L} and /F\partial\mathcal{L}/\partial F. Consequently, the Lagrangian density reads simply

=3Mr4er3/b3[r4b3+2a2r2(b33r3)cos2θ+3a4(3b32r3)cos4θ]b6Σ4,\displaystyle\mathcal{L}=\frac{3Mr^{4}e^{-r^{3}/b^{3}}[r^{4}b^{3}+2a^{2}r^{2}(b^{3}-3r^{3})\cos^{2}\theta+3a^{4}(3b^{3}-2r^{3})\cos^{4}\theta]}{b^{6}\Sigma^{4}}, (15)

and /F\partial\mathcal{L}/\partial F has the form

F=6Mr2er3/b3[3r5a2(4b33r3)cos2θ]8b8.\displaystyle\frac{\partial\mathcal{L}}{\partial F}=\frac{6Mr^{2}e^{-r^{3}/b^{3}}[3r^{5}-a^{2}(4b^{3}-3r^{3})\cos^{2}\theta]}{8b^{8}}. (16)

Interestingly, when we substitute a=0a=0 into Eqs. (15) and (16) as a result we obtain

=3Mer3/b3b3,F=9Mr7er3/b38b8.\mathcal{L}=\frac{3Me^{-r^{3}/b^{3}}}{b^{3}},\quad\frac{\partial\mathcal{L}}{\partial F}=\frac{9Mr^{7}e^{-r^{3}/b^{3}}}{8b^{8}}. (17)

These expressions are exactly similar to that of the spherically symmetric spacetime.

IV Energy conditions

Now we are much interested in computing the orthonormal basis also known as tetrads of the spacetime (12) so that we can analyze the matter associated with the rotating regular spacetime Neves:2014aba ; Bardeen:1972fi . The computation of the orthonormal basis gives the following form Neves:2014aba ; Bardeen:1972fi :

eμ(a)=((gttΩgtϕ)0000±grr0000gθθ0gtϕ/gϕϕ00gϕϕ),e^{(a)}_{\mu}=\left(\begin{array}[]{cccc}\sqrt{\mp(g_{tt}-\Omega g_{t\phi})}&0&0&0\\ 0&\sqrt{\pm g_{rr}}&0&0\\ 0&0&\sqrt{g_{\theta\theta}}&0\\ {g_{t\phi}}/{\sqrt{g_{\phi\phi}}}&0&0&\sqrt{g_{\phi\phi}}\end{array}\right), (18)

where Ω=gtϕ/gϕϕ\Omega=g_{t\phi}/{g_{\phi\phi}} is the angular velocity of the rotating regular black hole. We further determine the components of the energy-momentum tensor by using the following relation Neves:2014aba ; Bardeen:1972fi

T(a)(b)=eμ(a)eν(b)Gμν.T^{(a)(b)}=e^{(a)}_{\mu}e^{(b)}_{\nu}G^{\mu\nu}.

When we compute the components of the energy-momentum tensor T(a)(b)T^{(a)(b)}, then we realize that it has only diagonal components Neves:2014aba ; Bardeen:1972fi , i.e.,

T(a)(b)=diag(ρ,P1,P2,P3),T^{(a)(b)}=\mbox{diag}(\rho,P_{1},P_{2},P_{3}), (19)

where ρ\rho and PP are the matter density and the pressure, respectively. In case of the rotating regular black hole, the quantities ρ\rho and PP that appearing in (19) have the forms as follows

ρ\displaystyle\rho =\displaystyle= 6Mr4er3/b3b3Σ2=P1,\displaystyle\frac{6Mr^{4}e^{-r^{3}/b^{3}}}{b^{3}\Sigma^{2}}=-P_{1},
P2\displaystyle P_{2} =\displaystyle= 3Mr2er3/b3[3r3Σ+2b3(r2+2a2cos2θ)]b6Σ2=P3.\displaystyle-\frac{3Mr^{2}e^{-r^{3}/b^{3}}\left[-3r^{3}\Sigma+2b^{3}(r^{2}+2a^{2}\cos^{2}\theta)\right]}{b^{6}\Sigma^{2}}=P_{3}. (20)

When we compare them with the components of static counterpart (2), we find that the components of (IV) are different because of the black hole rotation. This also happens when we consider rotating Vaidya solution Carmeli:1975kg , it has in addition null radiation some another source. Having the expressions of matter density ρ\rho and pressure PiP_{i}, we can verify the weak energy condition for the spacetime (12) that requires the inequalities ρ0\rho\geq 0 and ρ+Pi0\rho+P_{i}\geq 0, (i=1,2,3i=1,2,3) to be satisfied.

Refer to captionRefer to caption
Refer to captionRefer to caption
Figure 1: Radial and angular dependence of matter density ρ\rho for the rotating regular black hole, for the given values of the rotation parameter a=0.4,0.6a=0.4,0.6, the parameter b=0.4,0.5b=0.4,0.5, and θ=π/4\theta=\pi/4 (x=cosθx=\cos\theta).
Refer to captionRefer to caption
Refer to captionRefer to caption
Figure 2: Radial and angular dependence of ρ+P2\rho+P_{2} for the rotating regular black hole, for the given values of the rotation parameter a=0.4,0.6a=0.4,0.6, the parameter b=0.4,0.5b=0.4,0.5, and θ=π/4\theta=\pi/4 (x=cosθx=\cos\theta).

Moreover, a straight forward calculation gives

ρ+P2=ρ+P3=3Mr2er3/b3[3r3Σ4a2b3cos2θ]b6Σ2,\rho+P_{2}=\rho+P_{3}=\frac{3Mr^{2}e^{-r^{3}/b^{3}}\left[3r^{3}\Sigma-4a^{2}b^{3}\cos^{2}\theta\right]}{b^{6}\Sigma^{2}}, (21)

which shows that the energy condition may be violated and it can be seen in the plots (cf. Fig. 1 and 2). We find that the matter density is always positive for all values of aa and bb as can be seen from the Fig. 1. It turns out that the rotating regular black hole solution may violate energy conditions; however, this happens for all the rotating regular black holes in the literature (see, e.g., Neves:2014aba ; Bambi:2013ufa ). Despite small violation of such solutions are important from phenomenologically and also they are important as astrophysical black holes are rotating.

V Properties of rotating regular spacetime

In this section, we will discuss some important properties of the rotating regular black hole solution, for instance, the curvature scalars, the horizons, and the ergospheres etc. It is important to study these basic properties of the spacetime from physical point of view.

V.1 Curvature invariants

In order to check the curvature singularity within the spacetime, we compute the curvature invariants or scalars of the rotating regular black hole spacetime. We find that the curvature invariants of the spacetime (12) have cumbersome forms, but in the limit θπ/2\theta\rightarrow\pi/2, they reduce to the following expressions

limθπ/2\displaystyle\lim_{\theta\rightarrow\pi/2}\mathcal{R} =\displaystyle= 18M2(9r612r3b3+8b6)e2r3/b3b12,\displaystyle\frac{18M^{2}(9r^{6}-12r^{3}b^{3}+8b^{6})e^{-2r^{3}/b^{3}}}{b^{12}},
limθπ/2K\displaystyle\lim_{\theta\rightarrow\pi/2}K =\displaystyle= 12M2e2r3/b3b12r6[4b12(1+e2r3/b3)4(3r6b6+2r3b9+2b12)er3/b3\displaystyle\frac{12M^{2}e^{-2r^{3}/b^{3}}}{b^{12}r^{6}}\Big{[}4b^{12}(1+e^{2r^{3}/b^{3}})-4(3r^{6}b^{6}+2r^{3}b^{9}+2b^{12})e^{r^{3}/b^{3}}
+27r12+24r6b6+8r3b9],\displaystyle+27r^{12}+24r^{6}b^{6}+8r^{3}b^{9}\Big{]},
limθπ/2R\displaystyle\lim_{\theta\rightarrow\pi/2}R =\displaystyle= 6M(3r34b3)er3/b3b6.\displaystyle-\frac{6M(3r^{3}-4b^{3})e^{-r^{3}/b^{3}}}{b^{6}}. (22)

We further consider the limit r0r\rightarrow 0 of (22), in this limit they turn out to be

limr0limθπ/2\displaystyle\lim_{r\rightarrow 0}\lim_{\theta\rightarrow\pi/2}\mathcal{R} =\displaystyle= 144M2b6,limr0limθπ/2K=96M2b6,limr0limθπ/2R=24Mb3.\displaystyle\frac{144M^{2}}{b^{6}},\quad\lim_{r\rightarrow 0}\lim_{\theta\rightarrow\pi/2}K=\frac{96M^{2}}{b^{6}},\quad\lim_{r\rightarrow 0}\lim_{\theta\rightarrow\pi/2}R=\frac{24M}{b^{3}}. (23)
Refer to captionRefer to caption
Figure 3: The behavior of Ricci square and Riemann square for the rotating regular black hole.

The curvature invariants have finite values for M,b0M,b\neq 0, therefore, we can say that the metric (12) is regular everywhere as can be seen from the Fig. 3. The figure shows that the Kretschmann scalar and the square of Ricci tensor are well behaving for different values of the parameters bb and aa. It is noticeable that the presence of the exponential factor in the mass of the black hole (12), i.e., er3/b3e^{-r^{3}/b^{3}}, makes the black hole singularity free.

V.2 Horizons

Now we wish to discuss the effect of charge on the structure of horizons and ergosphere. It turns out that as like the Kerr black hole, the spacetime has two surfaces, viz., static limit surfaces and horizons. The static limit surface is a surface on which the time translation Killing vectors (χa\chi^{a}) becomes null or χaχa=0\chi^{a}\chi_{a}=0. It requires that the gttg_{tt} component to be vanished,

r2+a2cos2θ2M(1er3/b3)r=0.r^{2}+a^{2}\cos^{2}\theta-2M(1-e^{-r^{3}/b^{3}})r=0. (24)

The typical behavior of the static limit surface is depicted in Fig. 4 for various values of the parameters (see also Table 1). It turns out that the behavior is similar to that of the Kerr black holes. However, the radii of the static limit surface shrink with increasing values of charge bb (Table 1).

Refer to caption
Refer to caption
Figure 4: Plots showing the variation of static limit surface with radius rr for different values of the parameter bb and the rotation parameter aa.
Table 1: Radius of EHs, SLSs and δa=rSLS+rEH+\delta^{a}=r^{+}_{SLS}-r^{+}_{EH} for different values of parameter bb
a=0.40a=0.40 a=0.45a=0.45 a=0.50a=0.50 a=0.55a=0.55
bb rH+r^{+}_{H} rsls+r^{+}_{sls} δ0.40\delta^{0.40} rH+r^{+}_{H} rsls+r^{+}_{sls} δ0.45\delta^{0.45} rH+r^{+}_{H} rsls+r^{+}_{sls} δ0.50\delta^{0.50} rH+r^{+}_{H} rsls+r^{+}_{sls} δ0.55\delta^{0.55}
0 1.91652 1.95917 0.04265 1.89303 1.94802 0.05499 1.86603 1.93541 0.06938 1.83516 1.92128 0.08612
0.20 1.91652 1.95917 0.04265 1.89303 1.94802 0.05499 1.86603 1.93541 0.06938 1.83516 1.92128 0.08612
0.40 1.91652 1.95917 0.04265 1.89303 1.94802 0.05499 1.86603 1.93541 0.06938 1.83516 1.92128 0.08612
0.80 1.91648 1.95915 0.04267 1.89298 1.94800 0.05502 1.86594 1.93539 0.06945 1.83502 1.92124 0.08622
1.10 1.90991 1.95491 0.04500 1.88464 1.94325 0.05861 1.85500 1.92997 0.07497 1.82009 1.91499 0.0949
1.20 1.89844 1.94703 0.04859 1.87049 1.93456 0.06407 1.83687 1.92027 0.08340 1.79566 1.90399 0.10833
1.30 1.81955 1.89605 0.07650 1.76301 1.87791 0.11490 1.79581 1.89955 0.10374 1.73422 1.88046 0.14624

Since the spacetime (12) has a coordinate singularity at Δ=0\Delta=0, which corresponds to the horizon of the rotating regular black holes. Actually the event horizon is located at the larger root (r+r_{+}) of

Δ=r2+a22M(1er3/b3)r,\Delta=r^{2}+a^{2}-2M(1-e^{-r^{3}/b^{3}})r, (25)

where (24) coincides with (25) when either θ=0\theta=0 or π\pi. Clearly, the radii of horizons depend on the parameter bb and they are different from the Kerr black hole. We solve (25) for horizons numerically as well as plot them in Fig. 5 for different values of parameters aa and bb. It turns out that the horizons of the rotating regular black hole have similar behavior like the Kerr black hole, but there exist several extremal black holes corresponding to the various values of charge bb.

Refer to caption
Refer to caption
Figure 5: Plots showing the variation variation of Δ\Delta with radius rr for different values of the parameter bb and the rotation parameter aa.
Table 2: Radius of inner horizon and outer horizon for different values of parameter bb and aa
a=0.40a=0.40 a=0.45a=0.45 a=0.50a=0.50 a=0.55a=0.55
bb rHr^{-}_{H} rH+r^{+}_{H} rHr^{-}_{H} rH+r^{+}_{H} rHr^{-}_{H} rH+r^{+}_{H} rHr^{-}_{H} rH+r^{+}_{H}
0 0.28085 1.91652 0.29917 1.89303 0.31791 1.86603 0.33734 1.83516
0.40 0.43462 1.91652 0.45736 1.89303 0.48063 1.86603 0.50460 1.83516
0.80 0.74894 1.91648 0.77631 1.89298 0.80536 1.86594 0.83625 1.83502
1.10 1.02687 1.90991 1.05876 1.88464 1.09424 1.85500 1.13402 1.82009
1.20 1.13547 1.89844 1.17071 1.87049 1.21110 1.83687 1.25840 1.79566
1.30 1.41374 1.81955 1.47841 1.76301 1.35248 1.79581 1.42108 1.73422

The Fig. 5 and Table 2 shows the existence of two roots, for a set of values of parameters aa and bb, which corresponds to the Cauchy horizon (smaller root) and event horizon (larger root). We find that for a given value of aa, there exists a critical value of parameter b=bcb=b_{c}, where the two horizons coincide (r=r±r=r_{\pm}) corresponding to the extremal black holes (Fig. 5 and 6). When b<bcb<b_{c}, we have rotating black hole with two horizons, and for b>bcb>b_{c} no black hole will form (Fig. 5 and 6).

V.3 Ergosphere

Refer to caption Refer to caption Refer to caption Refer to caption
Refer to caption Refer to caption Refer to caption Refer to caption
Refer to caption Refer to caption Refer to caption Refer to caption
Refer to caption Refer to caption Refer to caption Refer to caption
Figure 6: Plots showing the behavior of the ergoregion in the xz-plane for different values of the parameters bb and aa.

An ergosphere is bounded by the two above discussed surfaces, namely, the static limit surface and the event horizon. It lies outside the black hole. The ergoregion of the rotating black holes is a region of spacetime where the time translation Killing vector χa\chi^{a} becomes spacelike, and also that every timelike vector acquires rotational (ϕ\phi) counterpart. Hence, a particle can enter into ergoregion and can leave again, but it cannot remain stationary without the ergosphere. The ergospheres for the rotating regular black hole are shown in the Fig. 6 for different values of parameters bb, which are polar plots of (24) and (25).

Interestingly, we note that ergosphere becomes more prolate with increasing values of charge bb (see Fig. 6 horizontally) and with increasing spin aa (see Fig. 6 from top to bottom). It means the faster rotating and charged black holes are more prolate thereby increasing the area of the ergosphere. Furthermore, we observe from the Fig. 6 that we can find a critical parameter b=bcb=b_{c} where two surface shrinks to one, and for b>bcb>b_{c}, we have no ergoregion. It turns out that energy can be extracted from ergoregion via Penrose process Penrose:1971uk and the magnetic charge shall affect the Penrose process.

VI Black hole shadow

In this section, we are going to study one of the most interesting astrophysical phenomena for the rotating regular black hole known as black hole shadow. This can be understood as follows. When a black hole is placed in between an observer and the bright background source, the photons with a small angular momentum fall into the black hole without reaching the observer create a dark spot in the sky, which is called the black hole shadow. Thus, to obtain the apparent shape or shadow of the black hole and discuss the properties of shadow, we must discuss the geodesic of the photon in the background of the rotating regular black hole. These geodesics can be determined by the Hamilton-Jacobi formulation Carter:1968rr given by

Sσ=12gμνSxμSxν,\frac{\partial S}{\partial\sigma}=-\frac{1}{2}g^{\mu\nu}\frac{\partial S}{\partial x^{\mu}}\frac{\partial S}{\partial x^{\nu}}, (26)

which has a solution in the separable form as follows

S=12m02σEt+Lzϕ+Sr(r)+Sθ(θ),S=\frac{1}{2}m_{0}^{2}\sigma-Et+L_{z}\phi+S_{r}(r)+S_{\theta}(\theta), (27)

where m0m_{0} corresponds to the rest mass of the particle, and Sr(r)S_{r}(r) and Sθ(θ)S_{\theta}(\theta) are respectively functions of rr and θ\theta only. Note that in case of photon the rest mass, m0=0m_{0}=0. By following the standard procedure of the Hamilton-Jacobi method, we obtain the following forms of the geodesic equations,

Σdtdσ=a(aEsin2θLz)+(r2+a2)[(r2+a2)EaLz]r2+a22M(1er3/b3)r,\Sigma\frac{dt}{d\sigma}=-a\left(aE\sin^{2}\theta-L_{z}\right)+\frac{(r^{2}+a^{2})\left[(r^{2}+a^{2})E-aL_{z}\right]}{r^{2}+a^{2}-2M(1-e^{-r^{3}/b^{3}})r}, (28)
Σdϕdσ=(aELzcsc2θ)+a[(r2+a2)EaLz]r2+a22M(1er3/b3)r,\Sigma\frac{d\phi}{d\sigma}=-\left(aE-L_{z}\csc^{2}\theta\right)+\frac{a\left[(r^{2}+a^{2})E-aL_{z}\right]}{r^{2}+a^{2}-2M(1-e^{-r^{3}/b^{3}})r}, (29)
Σdrdσ=±,\Sigma\frac{dr}{d\sigma}=\pm\sqrt{\mathcal{R}}, (30)
Σdθdσ=±Θ,\Sigma\frac{d\theta}{d\sigma}=\pm\sqrt{\Theta}, (31)

where \mathcal{R} and Θ\Theta can be expressed as follows

\displaystyle\mathcal{R} =\displaystyle= [(r2+a2)EaLz]2[r2+a22M(1er3/b3)r][𝒦+(LzaE)2],\displaystyle\left[(r^{2}+a^{2})E-aL_{z}\right]^{2}-\left[r^{2}+a^{2}-2M(1-e^{-r^{3}/b^{3}})r\right]\left[\mathcal{K}+(L_{z}-aE)^{2}\right],
Θ\displaystyle\Theta =\displaystyle= 𝒦+cos2θ(a2E2Lz2csc2θ).\displaystyle\mathcal{K}+\cos^{2}\theta\left(a^{2}E^{2}-L_{z}^{2}\csc^{2}\theta\right). (32)

These geodesic equations represent first-order differential equations with respect to the affine parameter σ\sigma and contain three constants of the motion: the energy EE, the angular momentum LzL_{z}, and the Carter constant 𝒦\mathcal{K} Carter:1968rr . If we substitute b=0b=0, in the geodesic equations, they reduce to the Kerr spacetime case Hioki:2009na . Now in order to get the shadow of the rotating regular black hole, we introduce the two independent dimensionless quantities or impact parameters such that ξ=Lz/E\xi=L_{z}/E and η=𝒦/E2\eta=\mathcal{K}/E^{2}. Equation =0\mathcal{R}=0 explores the radial turning points of the photons and the spherical photon orbits around the black hole can be determined by the following relations Bardeen:1972fi

=0andddr=0.\mathcal{R}=0\quad\text{and}\quad\frac{d\mathcal{R}}{dr}=0. (33)

On solving (33) for impact parameters ξ\xi and η\eta, we immediately get their forms as follows

ξ\displaystyle\xi =\displaystyle= M[a2(b33r3)3r2(b3+r3)]b3er3/b3[a2(r+M)+r2(r3M)]aM(b33r3)+ab3(rM)er3/b3,\displaystyle\frac{M\left[a^{2}(b^{3}-3r^{3})-3r^{2}(b^{3}+r^{3})\right]-b^{3}e^{r^{3}/b^{3}}\left[a^{2}(r+M)+r^{2}(r-3M)\right]}{aM(b^{3}-3r^{3})+ab^{3}(r-M)e^{r^{3}/b^{3}}},
η\displaystyle\eta =\displaystyle= 4a2b3r3Mer3/b3[b3(er3/b31)3r3]r4[b3er3/b3(r3M)+3M(b3+r3)]2[aM(b33r3)+ab3(rM)er3/b3]2.\displaystyle\frac{4a^{2}b^{3}r^{3}Me^{r^{3}/b^{3}}\left[b^{3}(e^{r^{3}/b^{3}}-1)-3r^{3}\right]-r^{4}\left[b^{3}e^{r^{3}/b^{3}}(r-3M)+3M(b^{3}+r^{3})\right]^{2}}{\left[aM(b^{3}-3r^{3})+ab^{3}(r-M)e^{r^{3}/b^{3}}\right]^{2}}. (34)

It turns out that the impact parameters have a dependency on charge bb. If charge bb is set to be zero, then Eq. (34) reduces to the Kerr black hole case Hioki:2009na . This calculation is helpful to determine the celestial coordinates for a distant observer in order to get the shadow of the rotating regular black hole. The celestial coordinates for a distant observer in (α\alpha, β\beta)-plane are given Bardeen:1973gb by

α\displaystyle\alpha =\displaystyle= limr0(r02sinθ0dϕdr),\displaystyle\lim_{r_{0}\rightarrow\infty}\left(-r_{0}^{2}\sin\theta_{0}\frac{d\phi}{dr}\right),
β\displaystyle\beta =\displaystyle= limr0(r02dθdr),\displaystyle\lim_{r_{0}\rightarrow\infty}\left(r_{0}^{2}\frac{d\theta}{dr}\right), (35)

where r0r_{0} is the distance from the black hole to the observer and θ0\theta_{0} represents the inclination angle between the direction to observer and the rotation axis of the black hole. By using (28), (29), (30), (31), and (32), the celestial coordinates (35) transform into

α\displaystyle\alpha =\displaystyle= ξcscθ0,\displaystyle-\xi\csc\theta_{0},
β\displaystyle\beta =\displaystyle= ±η+a2cos2θξ2cot2θ0.\displaystyle\pm\sqrt{\eta+a^{2}\cos^{2}\theta-\xi^{2}\cot^{2}\theta_{0}}. (36)
Refer to caption Refer to caption Refer to caption Refer to caption
Figure 7: Plots showing the shapes of shadow with variation of inclination angle θ0\theta_{0} for fixed values of bb and aa.
Refer to caption Refer to caption Refer to caption Refer to caption
Figure 8: Plots showing the shapes of shadow for different values of the parameter bb and aa.

Equation (36) represents a direct relationship between the celestial coordinates (α,β\alpha,\beta) and the impact parameters (ξ,η\xi,\eta). Every single photon that approaches the observer determines a point on the (α\alpha, β\beta)-plane of the image, which can be visible with the help of telescope. The size and shape of the shadow depends on the parameters of black hole as well as on the inclination angle θ0\theta_{0}. The pictorial representation of the rotating regular black hole’s shadow can be seen from Figs. 7, 8 and 9. We obtain several images of the rotating regular black hole’s shadow for the particular choice of the parameters aa, bb, and θ0\theta_{0}. A distortion in the shape of the black hole’s shadow arises that increases with higher values of aa as well as θ0\theta_{0}. Fig. 9 depicts the effect of parameter bb on the shadow of the rotating regular black hole, which includes the case of the Kerr black hole (b=0b=0) for a comparison.

Refer to captionRefer to caption
Refer to captionRefer to caption
Figure 9: Plots showing the shapes of shadow for different values of the parameter bb.

Moreover, the actual size and distortion in shape of the rotating regular black hole’s shadow can be determined by the observables, namely, RsR_{s} and δs\delta_{s} Hioki:2009na . Here RsR_{s} and δs\delta_{s} correspond to the actual size or radius of the shadow and distortion in the shape of the shadow, respectively. In order to compute the radius of the shadow, we consider a circle passing through the three different points, namely, top, bottom, and rightmost corresponding to (αt,βt)(\alpha_{t},\beta_{t}), (αb,βb)(\alpha_{b},\beta_{b}), and (αr,0)(\alpha_{r},0), respectively Hioki:2009na . We are considering the definition of the observables given by Hioki and Maeda Hioki:2009na :

Rs\displaystyle R_{s} =\displaystyle= (αtαr)2+βt22|αtαr|,\displaystyle\frac{(\alpha_{t}-\alpha_{r})^{2}+\beta_{t}^{2}}{2|\alpha_{t}-\alpha_{r}|}, (37)
δs\displaystyle\delta_{s} =\displaystyle= (αp~αp)Rs,\displaystyle\frac{(\tilde{\alpha_{p}}-\alpha_{p})}{R_{s}}, (38)
Refer to caption
Figure 10: Illustration of the observable of the black hole shadow Amir:2016cen .

where (αp~,0)(\tilde{\alpha_{p}},0) and (αp,0)(\alpha_{p},0) are the points where the reference circle and the silhouette of the shadow cut the horizontal axis at the opposite side of (αr,0)(\alpha_{r},0) (cf., Fig. 10). This special characterization helps us to find the effect of charge bb on radius and distortion in the shape of the rotating regular black hole’s shadow. The typical behavior of these observables with charge bb can be seen in Fig. 11. We examine that an increase in the magnitude of charge bb increases the radius of the black hole’s shadow (RsR_{s}). On the other hand, the distortion in the shape of the black hole’s shadow (δs\delta_{s}) decreases with an increase in magnitude of charge bb (cf. Fig. 11).

Refer to caption
Refer to caption
Figure 11: Plots showing the behavior of observable RsR_{s} and δs\delta_{s} with parameter bb.

When a comparison of our results with the Kerr-Newman and the braneworld Kerr spacetimes is taking into account, we discover that the effect of nonlinear charge (our case) is opposite that of electric charge qq (Kerr-Newman) and tidal charge (braneworld Kerr). In our case, the radius of the black hole’s shadow increases with charge instead of being decreased at the same time the distortion decreases with charge instead of being increase.

VII Concluding remarks

We have shown that the regular black hole metric (1) is an exact solution of the Einstein’s field equations coupled to the nonlinear electrodynamics associated with the Lagrangian (8) corresponds to a magnetic charge bb. We have also constructed a rotating counterpart of the regular black hole (1) containing an additional parameter bb, which encompasses the Kerr black hole in the particular case when b=0b=0, and therefore belongs to a family of non-Kerr black holes. The source of the rotating regular black hole has also been computed in order to validate the solution. We have further discussed the energy conditions of the rotating regular black hole. The rotating regular spacetime describes the extremal black holes with degenerate horizons for a critical amount of charge b=bcb=b_{c}, the non-extremal black holes with two distinct horizons for b<bcb<b_{c} corresponding to the Cauchy and the event horizon. We have made a comprehensive analysis of horizon structure and ergosphere of the rotating regular black holes and explicitly brings out the effect of charge bb. The ergosphere is sensitive to the charge bb, which enlarges the ergoregion and becomes more prolate with an increasing magnitude of bb.

The immediate target of the Event Horizon Telescope is observations of the image by supermassive black hole Sgr A and M87 (shadow) which is a major attempt of understanding the nature of these black holes and explain the strong field gravity. Motivated from this, we have studied the shapes of the shadow from rotating regular black hole to discuss the effect of charge bb on the Kerr black hole shadow. We have done a detailed analysis of particle motion in the rotating regular black hole spacetimes. In order to achieve the goal, we have derived necessary analytical expressions to obtain the black hole shadow. The shadow is completely characterized by the two observables, namely, radius RsR_{s} and distortion δs\delta_{s}. It turns out that the radius RsR_{s} increases with an increase in the magnitude of charge bb, it results in a larger shadow of the rotating regular black hole than the Kerr black hole shadow. The rotating regular black hole is less distorted when compared with the Kerr black hole shadow as distortion δs\delta_{s} decreases with increasing charge bb. We have also compared our results with the Kerr-Newman and braneworld Kerr spacetimes; as a consequence, we found that in our case, the charge increases the radius of shadow and decreases the distortion in the shadow of the black hole. The results obtained in our study may be useful in the light of the observational outcome of the Event Horizon Telescope. The black hole shadow may help to conclude that whether the Kerr metric is an accurate description of astrophysical black holes or there is room for non-Kerr black holes like one we have discussed.

Acknowledgements.
We thank the DST INDO-SA bilateral project for Grant No. DST/INT/South Africa/P-06/2016. M.A. acknowledges that this research work is supported by the National Research Foundation, South Africa. S.D.M. acknowledges that this work is based upon research supported by South African Research Chair Initiative of the Department of Science and Technology and the National Research Foundation. We would like to thank IUCAA, Pune, for hospitality, where a part of this work was done.

References

  • (1) R. Penrose, Riv. Nuovo Cim.  1, 252 (1969) [Gen. Rel. Grav.  34, 1141 (2002)].
  • (2) A.D. Sakharov, Sov. Phys. JETP, 22, 241 (1966).
  • (3) E.B. Gliner, Sov. Phys. JETP, 22, 378 (1966).
  • (4) J. Bardeen, in Proceedings of GR5 (Tiflis, U.S.S.R., 1968)
  • (5) E. Ayón-Beato and A. García, Phys. Lett. B 493, 149 (2000)
  • (6) I. Dymnikova, Class. Quantum Gravity  21, 4417 (2004).
  • (7) K.A. Bronnikov, Phys. Rev. D 63, 044005 (2001).
  • (8) S. Shankaranarayanan and N. Dadhich, Int. J. Mod. Phys. D 13, 1095 (2004).
  • (9) S.A. Hayward, Phys. Rev. Lett.  96, 031103 (2006).
  • (10) S. Ansoldi, arXiv:0802.0330.
  • (11) H. Culetu, Int. J. Theor. Phys.  54, 2855 (2015).
  • (12) L. Balart and E.C. Vagenas, Phys. Lett. B 730, 14 (2014).
  • (13) L. Balart and E.C. Vagenas, Phys. Rev. D 90, 124045 (2014).
  • (14) L. Xiang, Y. Ling, and Y.G. Shen, Int. J. Mod. Phys. D 22, 1342016 (2013).
  • (15) C. Bambi, Mod. Phys. Lett. A 26, 2453 (2011).
  • (16) C. Bambi, Phys. Lett. B 730, 59 (2014).
  • (17) C. Bambi and L. Modesto, Phys. Lett. B 721, 329 (2013).
  • (18) B. Toshmatov, B. Ahmedov, A. Abdujabbarov, and Z. Stuchlík, Phys. Rev. D 89, 104017 (2014).
  • (19) S.G. Ghosh and S.D. Maharaj, Eur. Phys. J. C 75, 7 (2015).
  • (20) J.C.S. Neves and A. Saa, Phys. Lett. B 734, 44 (2014).
  • (21) B. Toshmatov, Z. Stuchlík, and B. Ahmedov, Phys. Rev. D 95, 084037 (2017).
  • (22) E.T. Newman and A.I. Janis, J. Math. Phys.  6, 915 (1965).
  • (23) M. Azreg-Aïnou, Eur. Phys. J. C 74, 2865 (2014).
  • (24) M. Azreg-Aïnou, Phys. Lett. B 730, 95 (2014).
  • (25) M. Azreg-Aïnou, Phys. Rev. D 90, 064041 (2014).
  • (26) S.G. Ghosh, Eur. Phys. J. C 75, 532 (2015).
  • (27) S.G. Ghosh, P. Sheoran, and M. Amir, Phys. Rev. D 90, 103006 (2014).
  • (28) M. Amir and S.G. Ghosh, J. High Energy Phys. 07 (2015) 015.
  • (29) S.G. Ghosh and M. Amir, Eur. Phys. J. C 75, 553 (2015).
  • (30) M. Amir, F. Ahmed and S.G. Ghosh, Eur. Phys. J. C 76, 532 (2016).
  • (31) F. Ahmed, M. Amir and S.G. Ghosh, Astrophys. Space Sci.  364, 10 (2019).
  • (32) B. Toshmatov, A. Abdujabbarov, Z. Stuchlík, and B. Ahmedov, Phys. Rev. D 91, 083008 (2015).
  • (33) B. Toshmatov, Z. Stuchlík, J. Schee, and B. Ahmedov, Phys. Rev. D 97, 084058 (2018).
  • (34) B. Toshmatov, Z. Stuchlík, and B. Ahmedov, Phys. Rev. D 98, 085021 (2018) .
  • (35) B. Toshmatov, Z. Stuchlík, B. Ahmedov, and D. Malafarina, Phys. Rev. D 99, 064043 (2019).
  • (36) S.  Chandrasekhar, The Mathematical Theory of Black Holes (Oxford University Press, New York, 1992).
  • (37) H. Falcke, F. Melia, and E. Agol, Astrophys. J.  528, L13 (2000).
  • (38) R. Takahashi, J. Korean Phys. Soc.  45, S1808 (2004) [Astrophys. J.  611, 996 (2004)].
  • (39) R. Takahashi, Publ. Astron. Soc. Jap.  57, 273 (2005).
  • (40) K. Hioki and K.i. Maeda, Phys. Rev. D 80, 024042 (2009).
  • (41) C. Bambi and K. Freese, Phys. Rev. D 79, 043002 (2009).
  • (42) C. Bambi and N. Yoshida, Class. Quant. Grav.  27, 205006 (2010).
  • (43) L. Amarilla, E.F. Eiroa, and G. Giribet, Phys. Rev. D 81, 124045 (2010).
  • (44) L. Amarilla and E.F. Eiroa, Phys. Rev. D 85, 064019 (2012).
  • (45) A. Yumoto, D. Nitta, T. Chiba, and N. Sugiyama, Phys. Rev. D 86, 103001 (2012).
  • (46) A.F. Zakharov, F. De Paolis, G. Ingrosso, and A.A. Nucita, New Astronomy Reviews 56 (2012) 64-73.
  • (47) L. Amarilla and E.F. Eiroa, Phys. Rev. D 87, 044057 (2013).
  • (48) F. Atamurotov, A. Abdujabbarov, and B. Ahmedov, Astrophys. Space Sci.  348, 179 (2013).
  • (49) F. Atamurotov, A. Abdujabbarov, and B. Ahmedov, Phys. Rev. D 88, 064004 (2013).
  • (50) S.W. Wei and Y.X. Liu, J. Cosmol. Astropart. Phys. 11 (2013) 063.
  • (51) A. Abdujabbarov, F. Atamurotov, Y. Kucukakca, B. Ahmedov, and U. Camci, Astrophys. Space Sci.  344, 429 (2013).
  • (52) T. Johannsen, Astrophys. J.  777, 170 (2013).
  • (53) A. Grenzebach, V. Perlick, and C. Lämmerzahl, Phys. Rev. D 89, 124004 (2014).
  • (54) Z. Li and C. Bambi, J. Cosmol. Astropart. Phys. 01 (2014) 041.
  • (55) U. Papnoi, F. Atamurotov, S. G. Ghosh, and B. Ahmedov, Phys. Rev. D 90, 024073 (2014).
  • (56) A.A. Abdujabbarov, L. Rezzolla, and B.J. Ahmedov, Mon. Not. Roy. Astron. Soc.  454, 2423 (2015).
  • (57) C. Goddi et al., Int. J. Mod. Phys. D 26, 1730001 (2016).
  • (58) A. Abdujabbarov, M. Amir, B. Ahmedov, and S.G. Ghosh, Phys. Rev. D 93, 104004 (2016).
  • (59) M. Amir and S.G. Ghosh, Phys. Rev. D 94, 024054 (2016).
  • (60) Z. Younsi, A. Zhidenko, L. Rezzolla, R. Konoplya, and Y. Mizuno, Phys. Rev. D 94, 084025 (2016).
  • (61) A. Grenzebach, The Shadow of Black Holes: An Analytic Description, Springer Briefs in Physics, Springer, Heidelberg (2016).
  • (62) M. Amir, B.P. Singh, and S.G. Ghosh, Eur. Phys. J. C 78, 399 (2018).
  • (63) R. Kumar, B.P. Singh, M.S. Ali, and S.G. Ghosh, arXiv:1712.09793.
  • (64) B.P. Singh and S.G. Ghosh, Ann. Phys.  395, 127 (2018).
  • (65) P.G. Nedkova, V.K. Tinchev, and S.S. Yazadjiev, Phys. Rev. D 88, 124019 (2013).
  • (66) C. Bambi, Phys. Rev. D 87, 107501 (2013).
  • (67) M. Azreg-Aïnou, J. Cosmol. Aastro. Phys. 1507 (2015) 037.
  • (68) T. Ohgami and N. Sakai, Phys. Rev. D 91, 124020 (2015).
  • (69) A. Abdujabbarov, B. Juraev, B. Ahmedov, and Z. Stuchlík, Astrophys. Space Sci.  361, 226 (2016).
  • (70) R. Shaikh, Phys. Rev. D 98, 024044 (2018).
  • (71) T. Ohgami and N. Sakai, Phys. Rev. D 94, 064071 (2016).
  • (72) M. Amir, A. Banerjee and S.D. Maharaj, Annals Phys.  400, 198 (2019).
  • (73) M. Amir, K. Jusufi, A. Banerjee and S. Hansraj, Class. Quantum Gravity  36, 215007 (2019).
  • (74) K. Akiyama et al. [Event Horizon Telescope Collaboration], Astrophys. J.  875, L1 (2019).
  • (75) K. Akiyama et al. [Event Horizon Telescope Collaboration], Astrophys. J.  875, L5 (2019).
  • (76) K. Akiyama et al. [Event Horizon Telescope Collaboration], Astrophys. J.  875, L6 (2019).
  • (77) I. Dymnikova, Gen. Relativ. Grav.  24, 235 (1992).
  • (78) S. Fernando, Int. J. Mod. Phys. D 26, 1750071 (2017).
  • (79) J.M. Bardeen, W.H. Press, and S.A. Teukolsky, Astrophys. J.  178, 347 (1972).
  • (80) M. Carmeli and M. Kaye, Annals Phys.  103, 97 (1977).
  • (81) R. Penrose and R.M. Floyd, Nature (London) 229, 177 (1971).
  • (82) B. Carter, Phys. Rev. 174, 1559 (1968).
  • (83) J.M. Bardeen, in Black holes, in Proceedings of the Les Houches Summer School, Session 215239, edited by C.  De Witt and B.S.  De Witt (Gordon and Breach, New York, 1973).