<<

arXiv:2010.03615v1 [gr-qc] 7 Oct 2020 [ n ae eeaie oa rirr ttc peial ymti m symmetric spherically static, arbitrary an to generalized later and § ‡ † ∗ shado the [ of rotat size and for shape shadow The BH BH. of of [ shape parameter The rotation sha the situated. in with [ is observed Synge observer deformation by the no calculated where was is BH there Schwarzschild and a disk for spacetimes circular given a the of generally since properties investigated is been important have several also determine shadows would BH of contour the to [ limit f field GL strong of in to aspects [ background tool Various important contexts . an diverse is of GL in theories The T alternative phenomen eclipse. other astronomical solar and interesting a shadow. most during 1919 the BH in of fact, Eddington one the In indeed e GL. is is distributed GL sources of light The spot virtue are as there the the universe if by dark and our left GR is M87 in which of BHs sky galaxy predictions server’s of the fundamental existence the in for of orbiting studies Telescope discover mathematical Horizon gas exciting earlier Event a [ be the of M87 to Recently, galaxy disk out Equa [ of shadow. emerged (4D) Field center hole been the Einstein’s have black four of which and usual fields solutions (GL) in the the GR GR, by as of predicted emerge objects which compact (BHs) holes Black xeietltssi h ekfil eiewihmksG standa a propos GR gravitation makes which of regime theory weak-field geometric the in a tests is experimental (GR) Relativity General 24 21 20 lcrncaddress: Electronic address: Electronic address: Electronic lcrncaddress: Electronic ftecnoro h hdwo oaigBswr investigated. were BHs rotating of wh [ shadow the the et.al for of Perlic shadow contour BH by the of by size of performed and ] shape successfully the is of ] study analytical The ]. eeto fLgtadSao atb ulCagdSrnyB Stringy Charged Dual a by Cast Shadow and Light of Deflection a ASnmes 45.h 76.f 04.70. 97.60.Lf, 04.50.+h, numbers: PACS space cha hole significant co black not have Schwarzschild are s We well-known there to also it. that compared by We observed hole cast and shadow studied. it of around is radius disk charges the obtain magnetic to and spacetime s angle charged electric deflection dual of of a values by variation light The of compac gravity. deflection l particular the dilaton-Maxwell of study the deflection gravi we of as paper, The well nature this as physical object universe. compact about gravitatin a a information the by of cast s field in shadow gravitational the the holes the of of by black light one of of are deflection existence shadows hole the black prove and lensing Gravitational eateto hsc,Grkl agiVswvdaaa H Vishwavidyalaya, Kangri Gurukula Physics, of Department b eateto hsc,Da ig olg,Uiest fDe of University College, Singh Dyal Physics, of Department [email protected] [email protected] [email protected] [email protected] hba Kala Shubham c etrfrSaeRsac,NrhWs nvriy Mahiken University, North-West Research, Space for Center 11 – 16 .Bzae l rtivsiae nltclyteG nteSchwarzsch the in GL the analytically investigated first al. et Bozza ]. 4 – 9 hc hw rgtrn ihadr,cnrlso.Tern sactu is ring The spot. central dark, a with ring bright a shows which ] 17 a .Tetcnqehssbeunl enapidtm n gi o diffe for again and time applied been subsequently has technique The ]. , ∗ Saurabh .INTRODUCTION I. b − , 19 † 22 s ewt Nandan Hemwati safnto ftems fteB n fterdu coordinate radius the of and BH the of mass the of function a as ] .Sm-nltccluain [ calculations Semi-analytic ]. rB pctmsaedsusdb ieetgroups different by discussed are spacetimes BH or h hp faB hdwfrannrttn BH non-rotating a for shadow BH a of shape The . oysc sabakhl.Investigation hole. black a as such body g eyhr,icuigfo eidteosre [ observer the behind from including verywhere, ulgoeisna H n te opc objects compact other and BHs near geodesics null n Hi oee olne iclradi varies it and circular longer no however is BH ing ihteipc aaee o different for parameter impact the with ti [ etric frttn Hwsfis acltdb Bardeen by calculated first was BH rotating of w 3 hc a rttse n ofimdb Arthur by confirmed and tested first was which a betadohrrltdapcs In aspects. related other and object t iei GR. in time a,c .Aogoeo h acntn consequences fascinating the of one Among ]. rnybakhl pctm eie in derived spacetime hole black tringy gsi h hdwcs yti black this by cast shadow the in nges h hdwi enda h eino h ob- the of region the as defined is shadow the ik n ad [ Maeda and Hioki gtaon tmypoieteuseful the provide may it around ight rda 4 0,Utrkad India Uttarakhand, 404, 249 aridwar dter fgaiaini u nvre[ universe our in gravitation of theory rd o ie h nua imtro h shadow the of diameter angular The size. dow ela rvdsapsil a ovrf one verify to way possible a provides as well db isenadhstuypse various passed truly has and Einstein by ed uytesao fti lc hole black this of shadow the tudy rnetosrainleiecsto evidences observational trongest , ET rvddtefis mg faB in BH a of image first the provided (EHT) h,NwDli100,Idaand India 110003, Delhi New lhi, sdrda pial hnemission thin optically an nsidered ‡ ainllnigaie u othe to due arises lensing tational hrceievrospoete fBsi GR in BHs of properties various characterize e natohsc r rvttoa lensing gravitational are astrophysics in ies l ls fPeasiDmasispacetimes Plebanski-Demianski of class ole 18 n rte Sharma Prateek and .Frhr h hoeia set related aspects theoretical the Further, ]. 75 ot Africa South 2745, g in EE)aeoeo h peculiar the of one are (EFEs) tions 23 i td hrfr uprsthe supports therefore study his n e ueia method numerical new a and ] 25 lodsusdarelation a discussed also ] a § akHole lack l Hspacetime BH ild lyabright a ally etBHs rent 1 , 10 2 ]. ]. 2 between the shape of contour of the shadow and the inclination angle and the spin parameter of the Kerr BH. Various other attempts to study the shadow in various theories of GR and other alternative theories of gravity have also made in recent times [26–32]. The general formula for the radius of the shadow of spherically symmetric BHs was considered by Perlic et al. as well as more recently by Konoplya [33, 34]. Some recent studies related to the shadows of the BHs are also focused towards the accretion disk properties and the images of the shadow cast by these ra- diated emissions [35, 36]. Simulations of those accretions has thrown quite a bit light on the central host of the source.

The non-trivial geometry of different types of BHs is of crucial interest in GR and other alternative theories of gravity [37]. All the four dimensional geometries, we choose to work with generally have a Schwarzschild limit (obtainable by setting the relevant parameter zero in the line element) and one such BH arise in context of theory [38, 39]. The formulation as mentioned in [33, 34] can therefore be applied to a dual charged stringy BH spacetime to observe the difference of the shape and size of the shadow cast by it along with its comparison with most general BH spacetime i.e. Schwarzschild BH in GR. It may further describe new examples of the contours of the shadows of such BHs in view of optically thin emission disk in a much finer way.

The outline of our paper is as follows. In section 2, we introduce the spacetime metric and obtained the equations of motion for light rays. Section 3 includes the bending angle and circular light orbits which are of crucial relevance for the formation of the shadow. In section 4, we have calculated the radius of sphere and then obtained the angular radius of the shadow. In section 5, we consider a simple model of accretion flow by considering an optically thin emission region surrounding the dual charge stringy BH to see the images of the shadows cast by such accretion flow. Finally, conclusions draw are summarised in section 6.

II. THE STRINGY BH SPACETIMES AND PHOTON MOTION

The well-known asymptotically flat solutions which represent BHs with the electric and dual magnetic charge in dilaton-Maxwell gravity are due to Garfinkle, Horowitz and Strominger (GHS) [38, 39].The 4D effective action to describe the above solutions reads as,

S = d4x√ g R + 2( φ)2 + e−2φF 2 , (1) − − ∇ Z   here, F is Maxwell field associated with U(1) sub group of E E and φ is dilaton, it is a scalar field which couples µν 8 × 8 to Maxwell field. The field equations yielded by extremizing the action with respect to U(1) potential Aµ, φ and gµν are

(e−2φF µν )=0, (2) ∇µ

1 2φ + e−2φF 2 =0, (3) ∇ 2

1 R =2 φ φ +2e−2φF F ρ g e−2φF 2. (4) µν ∇µ ∇ν µρ ν − 2 µν

For Fµν = 0, it reduces to standard Einstein scalar field action as in GR. The spacetime geometry arises from the field equations (2 - 4) have line elements those are causally similar to Schwarzschild geometry in GR as evident from the metric given below. The metric for the BH with electric charge is given as [40],

m 2 2 1 r 2 dr 2 2 dsE = − 2 dt + + r dΩ2, (5) − m sinh α 1 m 1+ r  − r    2 2 2 2 where dΩ2 = (dθ + sin θ dφ ) is the metric on a two dimensional unit sphere and α is a parameter related to the electric charge. The magnetically charged BH can be obtained from the electrically charged solution by electromagnetic duality transformation [39] and the metric with a magnetic charge parameter is given as follows,

m 2 2 1 r 2 dr 2 2 dsM = − 2 dt + 2 + r dΩ2. (6) − 1 Q  1 m 1 Q − mr − r − mr      3

Here Q is the magnetic charge of the BH. In the respective limits i.e. α = 0 and Q = 0, the Schwarzschild BH geometry as in GR in both the cases can be reproduced easily. Since the metric is spherically symmetric so each plane can be considered as an equatorial one, so that one can choose π ˙ θ = 2 and consequently θ = 0. Using the Hamilton-Jacobi equation, the equation of motion can be derived easily. The Hamiltonian for light rays given as,

2 2 2 1 ik 1 pt pr pφ H = g pipk = + + . (7) 2 2 −gtt grr gφφ !

The light rays are the solutions to the equation of motion with,

∂H ˙i ∂H p˙i = i , x = . (8) −∂x ∂pi We will now discuss the photon momentum and effective potentials for the cases corresponding to the line elements (5) and (6) respectively.

Case I [Corresponding to line element (5)]-

The momenta for electric charge stringy BH given as follows,

p m sinh2 α m −1 t˙ = t = E 1+ 1 , (9) −gtt − r − r    

˙ pφ L φ = = 2 , (10) −gφφ r

pr m r˙ = = pr 1 . (11) grr − r   Here E and L are represent the energy and angular momentum of test respectively. The effective potential for this case reads accordingly as below,

2 m L 1 r VE (r)= − . (12) r2 m sinh2 α 1+ r    Case II [Corresponding to line element (6)]-

The momenta for magnetic charge stringy BH is given as follows,

p Q2 m −1 t˙ = t = E 1 1 , (13) −gtt − − mr − r    

˙ pφ L φ = = 2 , (14) −gφφ r

2 pr m Q r˙ = = pr 1 1 . (15) grr − r − mr     The effective potential for this case reads as,

L2 1 m − r VM (r)= 2 2 . (16) r 1 Q  − mr   4

For H = 0, the Hamiltonian for light rays as described by equation (7) leads to,

p2 p2 p2 t + r + φ =0. (17) −gtt grr gφφ ! such that the radial momentum is given as follows,

2 2 grr 2 grr pr = pt pφ . (18) gtt − gφφ

Here the dot represents the derivatives with respect to an affine parameter and a prime denotes the derivatives with respect to r. The momenta pt and pφ are constants of motion and thus represent the energy and angular momentum of light rays respectively. From equations (15) and (14), one can obtain,

dr r˙ (g )p = = φφ r , (19) dφ φ˙ (grr)L and using pr from (18),

g 2 dr √ φφ E 2 = 2 h(r) 1 (20) dφ ± √grr r L − where, the function h(r) for electric and magnetic charge BH is calculated respectively as given below,

2 r2 1+ m sinh α 2 gφφ r hE(r) = = m , (21) gtt  1  − r and 

2 r2 1 Q 2 gφφ − mr hM (r) = = m . (22) gtt 1  − r  The effective potentials as given in equations (12) and (16) for Case I and Case II have been already discussed in [40] in detail. In order to observe the effect of dual charge, we present effective potential of the same BH for different values of electric charge(α) and magnetic charge(Q). The effective potentials for different values of α and Q and their comparison with Schwarzschild BH is shown in Fig.1. It can easily observed from effective potential that as the value of electric charge parameter increases the value of effective potential decreases (see Fig.1 Case I) and vice-versa for magnetic charge parameter (see Fig.1 Case II ). There is no minima present in each curve therefore there exist no stable orbit for the and only unstable orbits may exist in each case.

III. BENDING ANGLE AND CIRCULAR LIGHT ORBITS

Let us consider a light ray that comes in from infinity, reaches a minimum at a radius R, and goes out to infinity again. With such consideration, the integration over the orbit then leads to the following formula [33] for the bending angle δ,

∞ √grr dr δ = π +2 . (23) − g 2 R √ φφ pt 2 Z p2 h(r) 1 φ − r  Here R corresponds to the turning point of the trajectory and the condition dr needs to hold necessarily. This dφ |R=0 equation relates R to the constant of motion pφ/pt as below,

2 2 2 2 pφ L hM (R) = hE(R) = 2 = 2 . (24) pt E 5

FIG. 1: Effective potentials of the Stringy BH for different values of electric (Case I) and magnetic charge (Case II) respectively. The cases with α = 0, Q = 0 correspond to the Schwarzschild BH geometry in GR.

The deflection angle may then be rewritten as a function of R, ∞ √grr dr δ = π +2 . (25) − R √gφφ h(r)2 Z 2 1 h(R) − r  The expression of bending angle given by equation (25) takes the form for the Case I as, ∞ dr δE = π +2 . (26) 2 − R r2 1+ m sinh α Z  r  r(r m) 2 1 − b (1−m/r) − s  p However, the deflection angle for Case II reads as, ∞ dr δM = π +2 . (27) − R Q2 r2(1−Q2/mr) Z r(r m)(1 ) 2 1 − − mr b (1−m/r) − q r  where b = L/E is the impact parameter. The expression obtained in equation (26) and (27) with prescribed limits reduce to Schwarzschild BH case in GR. The plots represent the variation of deflection angle with impact parameter of BH for magnetic and electric charged stringy BH (see Fig.2). The cases with vanishing of charge parameters correspond to the Schwarzschild case as studied in [41]. In Fig.(2), every single curve for both cases indicate that, by increasing the value of impact parameter, the bending angle decreases for different values of Q and α. However, with an increase in the value of charge parameters, one can observe that the value of closest approach show the lateral inversion behaviour with one another. For magnetic charge (Q), the critical value of the closest approach decreases since the light goes closer to the BH however vice-versa for electric charge (α). A circular light orbits corresponds to zero radial velocity and acceleration, so thatr ˙ = 0 andr ¨ = 0. We will now discuss the circular light orbits for the above mentioned cases. Case I-

Using pr = 0 in equation (17) along with equation (11), we obtain

2 2 m sinh α E 1+ r L2 = , (28)  (1 m )  r2 − r and the angular momentum term can be separated as follows,

2 E2r2(1 + m sinh α ) L2 = r . (29) E (1 m/r) − 6

Differentiating equation (17) w.r.t affine parameter along with pr = 0, we have

mE2(1 + sinh2 α)r3 2L2(m r)2 − − =0. (30) (m 2)2r3 − and separating the angular momentum term from above expression, as

mE2(1 + sinh2 α)r3 L2 = . (31) E (m r)2 − Case II-

Using pr = 0 in equation (17) along with equation (15), we obtain

Q2 m −1 L2 E2 1 1 = , (32) − mr − r r2     and separating the angular term leads to,

r2E2(1 Q2/mr) L2 = − . (33) M (1 m/r) −

Differentiating equation 17 w.r.t affine parameter along with pr = 0, we have

2mL2(m r)2 + (m2 Q2)E2r3 − − − =0, (34) m(m r)2r3 − and once again separating the angular momentum term, we have

r3E2(m2 Q2) L2 = − . (35) M 2m(m r)2 − Subtracting the set of these two equations i.e. (29), (31) and (33), (35) from each other yields, the radius of a circular light orbit in the form d 0= h(r)2. (36) dr

The solution of equation (36) at r = rph gives the radius of photon sphere for both the cases as discussed in detail in the forthcoming section.

IV. RADIUS OF THE SHADOW

The motion of light rays near the photon sphere are crucial to understand the properties of BH shadows. The radius rph of the outermost photon sphere is the critical rph value of the minimal radius R mentioned above. If a light ray comes in from infinity and reaches a minimum radius R bigger than rph, it will go out to infinity again. The case R = rph corresponds to a light ray that spirals asymptotically towards a circular photon orbit in the sphere of radius rph. We consider a light ray that is sent from the observer’s position at r0 into the past under an angle β [33] with respect to the radial direction such that

√grr dr cot β = R=0. (37) √gφφ dφ |

The equation (20) may then be written in terms of minimal radius R as follows

2 dr √grr h (r) = 2 1. (38) dφ ±√gφφ sh (R) − 7

FIG. 2: Deflection angle with the impact parameter of stringy BH for different values of α (Case I) and Q (Case II) with m = 1. Here, the impact parameter is a function of the closest distance of approach.

Using equations (37) and (38),

2 2 h(R) sin β = 2 . (39) h(r0)

The boundary of the shadow βs is determined by light rays that spiral asymptotically towards a circular light orbit at radius rph. We will now first discuss the photon orbits for the Case I and Case II to obtain βs.

Case I-

From equation (36), the radius of photon sphere for BH with electric charge is calculated as,

E 1 2 m cosh α 17 + cosh(2α) rph = 3m m sinh α + . (40) 4 " − p√2 # The critical value of minimal radius R is however given as,

E E 2 rph(rph + m sinh α) RE = . (41) m 1 rE − ph   The angular radius of the shadow of BH with electric charge as described by equation (39) is then obtained as follows,

2 2 m sinh α m 2 r 1+ 1 2 E h(rph) ph rph − r0 sin βs = 2 = 2 , (42) h(r0) r2 1 m 1+m sinh α  0 − rph r0    where rph has to be determined from equation (40).

Case II-

The radius of photon sphere for BH with magnetic charge is calculated as,

2 2 4 2 2 2 M (3m/2+ Q /2m)+ 9m /4+ Q /4m +3Q /2 4Q rph = − . (43) p 2 8

The critical value of the minimal radius R thus reads as,

2 rM rM Q ph ph − m RM = . (44)  m  1 rM − ph   Therefore the angular radius of the shadow of BH with magnetic charge is then obtained as follows,

2 2 m Q 2 r 1 1 2 M h(rph) ph − r0 − mrph sin βs = 2 = 2 , (45) h(r0) r2 1 m  1 Q  0 − rph − mr0    where rph has to be determined from equation (43).

The expressions obtained in equations (42) and (45) with prescribed limits reduce to Schwarzschild BH case in GR [19].

V. OPTICALLY THIN EMISSION DISK SURROUNDED BY A STRINGY BH

Here, we consider a simple model of optically thin, radiating accretion flow around the stringy BH. For the emission mechanisms, certain assumptions shall be made for the calculation of the intensity from the radiating accretion flow. The observed specific intensity Iν0 at the observed photon frequency νobs at the point (X, Y ) of the observer’s image (usually measured in ergs−1cm−2str−1Hz−1) is given by[36, 42],

3 Iobs(νobs,X,Y )= g j(νe)dlprop, (46) Zγ where g = νobs/νe is the redshift factor, νe is the photon frequency as measured in the rest-frame of the emitter, j(νe) ρ is the emissivity per unit volume in the rest-frame of the emitter, and dlprop = kρuedλ is the infinitesimal proper length as measured in the rest-frame of the emitter . The redshift factor as in equation (46) is calculated from,

ρ kρuobs g = σ , (47) kσue

µ ρ µ where k is the 4-velocity of the photons, ue 4-velocity of the accreting gas emitting the radiation and uobs is 4-velocity of the observer with λ is the affine parameter along the photon path γ. The whole integral given by equation (46) should be evaluated over the path γ of the photons i.e. (null geodesics). In order to generalise the formalism for both the cases, we will restrict to the equatorial plane only and assume that the gas is radially free falling with a four-velocity whose components are described below,

t 1 ue = , gtt(r)

r 1 gtt(r) ue = − , −sgtt(r)grr(r) θ ue = 0, φ ue = 0, (48) where the metric components gtt, grr and gφφ are corresponding to the line elements for Case I and Case II. In previous section, the components of four-velocity were already calculated. In order to ease further calculations, we obtain an equation between the radial and time component of the four-velocity as follows,

k 1 b2 r = g , (49) k ± rr g − g t s  tt φφ  9

(a)

(b)

(c)

FIG. 3: (Case I)-The first column represents the intensity distribution with respect to the impact parameter and the second column represents images of optically thin emission region surrounding the electrically charge stringy BH with m = 1, for (a) α = 0 (i.e. Schwarzschild BH), (b) α = 0.25 and (c) α = 0.50. 10

(a)

(b)

FIG. 4: (Case-II)-The first column represents the intensity distribution with respect to the impact parameter and the second column represents images of optically thin emission region surrounding the magnetically charge stringy BH with m = 1, for (a) Q = 0.25 and (b) Q = 0.50. The Schwarzschild BH case is similar as in Fig.3(a)(Case I). where the sign +( ) corresponds to photon approaches (goes away) from the massive object. The redshift function g is then given by,− 1 g = . (50)

1 kr 1−gtt gtt ± kt gtt grr s  The specific emissivity for a simple model in which the emission is monochromatic with emitter’s-rest frame frequency 2 ν⋆, and the emission has a 1/r radial profile is calculated as, δ (ν ν ) j(ν ) D e − ⋆ , (51) e ∝ r2 where δD is the Dirac delta function. The proper length can be written as, k dl = k uρdλ = t dr, (52) prop ρ e −g kr | | Integrating the Intensity over all the observed frequencies, we obtain the observed photon intensity as given below, g3k I (X, Y ) t dr. (53) obs ∝− r2kr Zγ 11

Now we have the desired intensity equation so we can proceed to compute the images. The numerical routine was implemented in a python package EinsteinPY [43] and modified for the use of stringy BH spacetimes as Case I and Case II accordingly.

VI. SUMMARY, CONCLUSION AND FUTURE DIRECTIONS

In this paper, we have studied deflection angle and BH shadow of a dual charged stringy BH. First, we have analysed the null geodesics and the motion of a photon around this BH spacetime in brief. Using Hamilton-Jacobi formalism, we find out the equations of motion and the deflection angle of light around the above stringy BH spacetimes accordingly. The maximum bending angle takes place at a critical value of impact parameter and this critical value increases with an increase in the electric charge. In contrast, in case of magnetically charge stringy BH, the critical value for maximum deflection angle shifts towards a lower value as magnetic charge increases. However, with an increase in the value of electric and magnetic charge parameters as compared to the Schwarzschild BH case, it is observed that the value of closest approach shows the lateral inversion behaviour with one another. Further, we obtained the circular orbits for photons which gives us the necessary condition that whether it forms an unstable photon sphere or not. An unstable photon sphere constitutes for the shadow of the BH and the radius of the photon spheres is calculated accordingly for both the cases of dual charged stringy BH spacetime. The boundaries of the shadow were then found along with the angular radius. It is observed that the obtained expressions of radius of photon sphere and angular size of BH shadow with prescribed limits reduce to the Schwarzschild BH geometry in GR. Finally, we consider an optically thin emission disk around the BH (simple inverse square radial profile) and see the images produced by it. The images (see Figs.3, 4) produced in it show that there is not much significant changes for different values of electric and magnetic charge. The image while compared to the case of Schwarzschild BH seems indistinguishable however difference in deflection angle of a stringy BH for both cases as compared to Schwarzschild BH is clearly distinguishable. We therefore need a careful attention to distinguish among the geometries of stringy BHs and Schwarzschild BH from the view point of their optical and other properties.The stability of circular geodesics of these stringy charged BHs needs careful attention and we intended to report on this issue in near future. Further, there will also be detail investigation of lensing and shadow phenomenon for a Kerr-like BHs in various alternative theories of gravity in our future studies.

VII. ACKNOWLEDGMENTS

The authors are grateful to the anonymous referee for the valuable comments and suggestions which helped us to improve the quality of this paper. The author SK is thankful to the Uttarakhand State Council of Science and Technology (UCOST), Dehradun for financial assistance through R&D grant number UCS&T/RD-18/18-19/16038/4. The author HN and PS acknowledge the financial support provided by Science and Engineering Research Board (SERB), New Delhi through the grant number EMR/2017/000339. The authors (SK, HN and PS) also acknowledge the facilities at ICARD, Gurukula Kangri Vishwavidyalaya, Haridwar. The authors also thank to Dr. Rajibul Shaikh for useful suggestions.

[1] J. B. Hartle, Gravity: An introduction to einstein’s (2003). [2] P. S. Joshi, Global aspects in gravitation and cosmology., Int. Ser. Monogr. Phys. 87. [3] S. Chandrasekhar, The mathematical theory of black holes, Vol. 69, Oxford university press, 1998. [4] K. Akiyama, et al., First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive , Astrophys. J. 875 (1) (2019) L1. arXiv:1906.11238, doi:10.3847/2041-8213/ab0ec7. [5] K. Akiyama, et al., First M87 Event Horizon Telescope Results. II. Array and Instrumentation, Astrophys. J. 875 (1) (2019) L2. arXiv:1906.11239, doi:10.3847/2041-8213/ab0c96. [6] K. Akiyama, et al., First M87 Event Horizon Telescope Results. III. Data Processing and Calibration, Astrophys. J. 875 (1) (2019) L3. arXiv:1906.11240, doi:10.3847/2041-8213/ab0c57. [7] K. Akiyama, et al., First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole, Astrophys. J. 875 (1) (2019) L4. arXiv:1906.11241, doi:10.3847/2041-8213/ab0e85. [8] K. Akiyama, et al., First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring, Astrophys. J. 875 (1) (2019) L5. arXiv:1906.11242, doi:10.3847/2041-8213/ab0f43. [9] K. Akiyama, et al., First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole, Astrophys. J. 875 (1) (2019) L6. arXiv:1906.11243, doi:10.3847/2041-8213/ab1141. [10] S. E. Gralla, D. E. Holz, R. M. Wald, Black Hole Shadows, Photon Rings, and Lensing Rings, Phys. Rev. D 100 (2) (2019) 024018. arXiv:1906.00873, doi:10.1103/PhysRevD.100.024018. 12

[11] S. Liebes, Gravitational Lenses, Phys. Rev. 133 (1964) B835–B844. doi:10.1103/PhysRev.133.B835. [12] S. Refsdal, J. Surdej, Gravitational lenses, Rept. Prog. Phys. 57 (1994) 117–186. doi:10.1088/0034-4885/57/2/001. [13] A. M. Nzioki, P. K. S. Dunsby, R. Goswami, S. Carloni, A Geometrical Approach to Strong Gravitational Lensing in f(R) Gravity, Phys. Rev. D83 (2011) 024030. arXiv:1002.2056, doi:10.1103/PhysRevD.83.024030. [14] R. S. Kuniyal, H. Nandan, U. Papnoi, R. Uniyal, K. Purohit, Strong lensing and observables around 5d myers–perry black hole spacetime, Modern Physics Letters A 33 (23) (2018) 1850126. [15] R. Uniyal, H. Nandan, P. Jetzer, Bending angle of light in equatorial plane of Kerr–Sen Black Hole, Phys. Lett. B782 (2018) 185–192. arXiv:1803.04268, doi:10.1016/j.physletb.2018.05.006. [16] M. Azreg-A¨ınou, S. Bahamonde, M. Jamil, Strong Gravitational Lensing by a Charged Kiselev Black Hole, Eur. Phys. J. C77 (6) (2017) 414. arXiv:1701.02239, doi:10.1140/epjc/s10052-017-4969-4. [17] V. Bozza, S. Capozziello, G. Iovane, G. Scarpetta, Strong field limit of black hole gravitational lensing, Gen. Rel. Grav. 33 (2001) 1535–1548. arXiv:gr-qc/0102068, doi:10.1023/A:1012292927358. [18] V. Bozza, Gravitational lensing in the strong field limit, Phys. Rev. D66 (2002) 103001. arXiv:gr-qc/0208075, doi:10.1103/PhysRevD.66.103001. [19] J. L. Synge, The Escape of Photons from Gravitationally Intense Stars, Mon. Not. Roy. Astron. Soc. 131 (3) (1966) 463–466. doi:10.1093/mnras/131.3.463. [20] S. Hawking, B. Carter, J. M. Bardeen, H. Gursky, K. S. Thorne, R. Ruffini, I. D. Novikov, et al., Black Holes, Vol. 23, CRC Press, 1973. [21] A. Grenzebach, V. Perlick, C. L¨ammerzahl, Photon Regions and Shadows of Kerr-Newman-NUT Black Holes with a Cosmological Constant, Phys. Rev. D89 (12) (2014) 124004. arXiv:1403.5234, doi:10.1103/PhysRevD.89.124004. [22] A. Grenzebach, V. Perlick, C. L¨ammerzahl, Photon Regions and Shadows of Accelerated Black Holes, Int. J. Mod. Phys. D24 (09) (2015) 1542024. arXiv:1503.03036, doi:10.1142/S0218271815420249. [23] T. Johannsen, Photon Rings around Kerr and Kerr-like Black Holes, Astrophys. J. 777 (2013) 170. arXiv:1501.02814, doi:10.1088/0004-637X/777/2/170. [24] Z. Younsi, A. Zhidenko, L. Rezzolla, R. Konoplya, Y. Mizuno, New method for shadow calculations: Ap- plication to parametrized axisymmetric black holes, Phys. Rev. D94 (8) (2016) 084025. arXiv:1607.05767, doi:10.1103/PhysRevD.94.084025. [25] K. Hioki, K.-i. Maeda, Measurement of the Kerr Spin Parameter by Observation of a Compact Object’s Shadow, Phys. Rev. D80 (2009) 024042. arXiv:0904.3575, doi:10.1103/PhysRevD.80.024042. [26] R. Narayan, M. D. Johnson, C. F. Gammie, The Shadow of a Spherically Accreting Black Hole, Astrophys. J. 885 (2) (2019) L33. arXiv:1910.02957, doi:10.3847/2041-8213/ab518c. [27] S. Schneider, V. Perlick, The shadow of a collapsing dark star, Gen. Rel. Grav. 50 (6) (2018) 58. arXiv:1802.04901, doi:10.1007/s10714-018-2379-z. [28] P. V. P. Cunha, C. A. R. Herdeiro, B. Kleihaus, J. Kunz, E. Radu, Shadows of Einstein–dilaton–Gauss–Bonnet black holes, Phys. Lett. B768 (2017) 373–379. arXiv:1701.00079, doi:10.1016/j.physletb.2017.03.020. [29] S. Vagnozzi, C. Bambi, L. Visinelli, Concerns regarding the use of black hole shadows as standard rulersarXiv:2001.02986, doi:10.1088/1361-6382/ab7965. [30] O. Yu. Tsupko, Z. Fan, G. S. Bisnovatyi-Kogan, Black hole shadow as a standard ruler in cosmology, Class. Quant. Grav. 37 (6) (2020) 065016. arXiv:1905.10509, doi:10.1088/1361-6382/ab6f7d. [31] U. Papnoi, F. Atamurotov, S. G. Ghosh, B. Ahmedov, Shadow of five-dimensional rotating Myers-Perry black hole, Phys. Rev. D90 (2) (2014) 024073. arXiv:1407.0834, doi:10.1103/PhysRevD.90.024073. [32] F. Atamurotov, B. Ahmedov, Optical properties of black hole in the presence of plasma: shadow, Phys. Rev. D92 (2015) 084005. arXiv:1507.08131, doi:10.1103/PhysRevD.92.084005. [33] V. Perlick, O. Yu. Tsupko, G. S. Bisnovatyi-Kogan, Influence of a plasma on the shadow of a spherically symmetric black hole, Phys. Rev. D92 (10) (2015) 104031. arXiv:1507.04217, doi:10.1103/PhysRevD.92.104031. [34] R. A. Konoplya, Shadow of a black hole surrounded by dark , Phys. Lett. B795 (2019) 1–6. arXiv:1905.00064, doi:10.1016/j.physletb.2019.05.043. [35] R. Shaikh, P. Kocherlakota, R. Narayan, P. S. Joshi, Shadows of spherically symmetric black holes and naked singularities, MNRAS 482. arXiv:arXiv:1802.08060, doi:10.1093/mnras/sty2624. [36] M. Jaroszynski, A. Kurpiewski, Optics near kerr black holes: spectra of advection dominated accretion flows, Astron. Astrophys. 326 (1997) 419. arXiv:astro-ph/9705044. [37] G. G. L. Nashed, S. Capozziello, Charged spherically symmetric black holes in f(r) gravity and their stability analysis, Phys. Rev. D 99 (2019) 104018. doi:10.1103/PhysRevD.99.104018. URL https://link.aps.org/doi/10.1103/PhysRevD.99.104018 [38] D. Garfinkle, G. T. Horowitz, A. Strominger, Charged black holes in , Phys. Rev. D43 (1991) 3140, [Erratum: Phys. Rev.D45,3888(1992)]. doi:10.1103/PhysRevD.43.3140,10.1103/PhysRevD.45.3888. [39] G. T. Horowitz, The dark side of string theory: Black holes and black strings., in: In *Trieste 1992, Proceedings, String theory and ’92* 55-99, 1992. arXiv:hep-th/9210119. [40] A. Dasgupta, H. Nandan, S. Kar, Kinematics of geodesic flows in stringy black hole backgrounds, Phys. Rev. D79 (2009) 124004. arXiv:0809.3074, doi:10.1103/PhysRevD.79.124004. [41] J. Gerard, S. Pireaux, The Observable light deflection anglearXiv:gr-qc/9907034. [42] C. Bambi, Can the supermassive objects at the centers of galaxies be traversable wormholes?, Phys. Rev. D 87. doi:10.1103/PhysRevD.87.107501. [43] S. Bapat, et al., einsteinpy/einsteinpy: Einsteinpy 0.2.1 (2019). doi:10.5281/zenodo.3526198.