Front. Phys. Frontiers in Physics Front. Phys. 2296-424X Frontiers Media S.A. 1113909 10.3389/fphy.2023.1113909 Physics Review A review on analytical studies in gravitational lensing Chowdhuri et al. 10.3389/fphy.2023.1113909 Chowdhuri Abhishek * Ghosh Saptaswa * Bhattacharyya Arpan * Indian Institute of Technology Gandhinagar, Gandhinagar, India

Edited by: Farruh Atamurotov, Inha University in Tashkent, Uzbekistan

Reviewed by: Hernando Quevedo, National Autonomous University of Mexico, Mexico

Jan Schee, Silesian University in Opava, Czechia

 *Correspondence: Abhishek Chowdhuri, chowdhuri_abhishek@iitgn.ac.in; Saptaswa Ghosh, saptaswaghosh@iitgn.ac.in; Arpan Bhattacharyya, abhattacharyya@iitgn.ac.in

This article was submitted to High-Energy and Astroparticle Physics, a section of the journal Frontiers in Physics

 21 03 2023 2023 11 1113909 01 12 2022 20 02 2023 Copyright © 2023 Chowdhuri, Ghosh and Bhattacharyya. 2023 Chowdhuri, Ghosh and Bhattacharyya

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

In this study, we review some current studies on gravitational lensing for black holes, mainly in the context of general relativity. We mainly focus on the analytical studies related to lensing with references to observational results. We start with reviewing lensing in spherically symmetric Schwarzschild spacetime, showing how to calculate deflection angles before moving to the rotating counterpart, the Kerr metric. Furthermore, we extend our studies for a particular class of newly proposed solutions called black-bounce spacetimes and discuss throughout the review how to explore lensing in these spacetimes and how the various parameters can be constrained using available astrophysical and cosmological data.

 gravitational lensing strong lensing Einstein ring black holes deflection angle Department of Atomic Energy, Government of India10.13039/501100001502 1 Introduction

The importance of gravitational lensing began with Eddington’s observation of the light deflection of the Sun [1]. The experiment proved to be one of the major milestones in favor of Einstein’s theory of general relativity [2]. The theory has proved itself to be the most successful theory of classical gravity. However, GR has some hiccups in the form of singularities, not being able to understand the quantum theory of it, along with the ubiquitous nature of dark matter and dark energy, forcing us to conclude that maybe there is more to the story. Following this line of thought, efforts have been made to introduce modifications to the Einstein–Hilbert Lagrangian [3–14]. However, taking in a cue from Eddington, gravitational lensing still remains one of the most indispensable tools to check and put constraints on our theory parameters.

The theory behind gravitational lensing started when O. Lodge asked a simple question in [15] on the behavior of light in a gravitational field. Interestingly, he concluded that, unlike a convex lens, the gravitational field does not focus light rays into a focal point. However, a logarithmically shaped concave lens can indeed mimic the effect of a spherically symmetric gravitational field onto the light. Inspired by such insights, O. Chwolson came up with the possibility of observing ring-like images in cases of axial symmetry. Nowadays, these ring-like images are called Einstein rings [16, 17]. Remarkably, such images were indeed observed by J. Hewitt et al. using the Very Large Array, a collection of radio telescopes in the United States. The radio source was MG1131 + 0456, and it was found when a background (radio) Galaxy is distorted into an almost closed ring [18–20]. Later, such rings were also found in the infrared and even in the optical spectrum [21] and more recently by the James Webb Telescope [22]. Typically, the diameter ranges between a few arcseconds or less.

The discovery of quasars and such rings lead to a rapid change in the field of gravitational lensing. We are now well aware of lensing techniques such as weak lensing, where the small deformations of many background galaxies are statistically evaluated for determining the (dark) matter in a foreground Galaxy cluster, and microlensing, where the light curve of a star is registered that moves transversely to the line of sight behind a (dark and compact) mass. For a thorough exposition of gravitational lensing, including an overview of all pre-1992 observations, the reader may consult the monograph [23]. More up-to-date information can be found in [24].

In all the lensing observations mentioned previously, the bending angles are so small that a weak-field approximation for the gravitational field is applicable. A much simpler quasi-Newtonian formalism has been used. This approximation centers around a so-called lens map or lens equation, which is very intuitive. This lens map of the weak-field formalism has proven extremely useful for evaluating the aforementioned lensing phenomena. It is discussed in detail, for example, in [23]. Additional material can be found, for example, in the Living Review by Wambsganss [25], which is completely based on this approximation formalism. However, there are astrophysical scenarios where the bending angles are not small. In such cases, the weak lensing theory is not valid anymore. The objective of this review is to give an overview of the various analytical techniques used to study strong lensing and to give some insights into the relation of lensing parameters to the gravity theory itself.

This expectation is nurtured mainly by the increasing evidence of a black hole at the center of our Galaxy. Black holes arise as solutions of GR and play an important role in testing various aspects of it. Recent observations coming from the “Event Horizon Telescope” (EHT) [26–38] collaboration and LIGO, Virgo, and KAGRA collaboration [39–44] have provided us with enough observational evidence that ensures the existence of supermassive black holes. Furthermore, these observations provide a unique opportunity to test various aspects of strong gravity [45–55]. This necessitates us to go beyond weak-field approximation and consider lensing by the strong gravitational field, thereby motivating us to consider the strong-field limit of the bending angle. It provides important information about the intrinsic parameters such as spin, mass, and charge of the black holes and the parameters of the underlying gravity theory. Hence, its observational implications have been investigated in recent times 1 . An analytical description of the deflection angle for the spherically symmetric black hole in the strong gravitational field has been proposed in [76]. This was further based on [77, 78]. Later, it was extended for rotating black holes in [79]. A lens equation has been derived in [77], which has been further generalized in [80] where the underlying black hole serves as a lens. Numerous investigation of Einstein rings [81] and strong gravitational lensing based on the work has been made for various black hole spacetimes and compact objects 2 .

It should be noted that in this review, we will be focusing on the analytical computation of the deflection angle of light rays using null geodesic equations. However, there are other methods of calculating the deflection angle analytically, for instance, the material medium approach. In this approach, one maps the problem effectively to a problem of light propagation through a medium with a particular refractive index determined by the strength of the gravitational field of the original spacetime for which one wants to calculate the deflection angle. For more details, interested readers are referred to [121–133] and citations of these references.

Furthermore, in recent times, several static metrics known as the black-bounce metric have been proposed in [134–137]. What makes this spacetime interesting is that they have an extra parameter which regularizes the central singularity, unlike the usual black hole spacetime. The black-bounce spacetime interpolates between a black hole and (traversable) wormhole metric depending on the choice of underlying regularization parameters. When the solution does not admit any horizon, it corresponds to a wormhole solution. Recently, the gravitational lensing in the strong deflection limit for these black-bounce spacetimes has been studied [120, 138–141]. From the lensing perspective, this kind of spacetime provides us with extra tunable parameters; hence, it is interesting to investigate the effect of this extra parameter from the theoretical and observational points of view. In this review, besides discussing some analytical results about the computation of the strong deflection angle in the background of Schwarzschild and Kerr–Newman black holes, we also discuss the effect of this black-bounce regularization parameter on the computation of the deflection angle and, finally, on the radius of the Einstein ring.

This article is organized as follows. In Section 2, we first review the deflection angle, the angular radius of Einstein rings, and Shapiro time delay for a Schwarzschild spacetime. In Section 3, we review a general formalism of computation of the equatorial deflection angle for Kerr–Newman spacetime. Then, in Section 4, we discuss the computation of the deflection angle for the black-bounce metric. We also review the computation of the Einstein ring radius of this case and comment on its dependence on charge and regularization parameters and its observational implications. Finally, we give concluding remarks in Section 5. Some necessary details are given in Appendix A. Also, we have set the value of the speed of light c and Newton’s gravitational constant G to unity.

 2 Lensing for the spherically symmetric Schwarzschild black hole

In this section, we start by reviewing gravitational lensing in the simplest setup, that is, for the Schwarzschild solution. The metric reads d s 2 = − 1 − 2 M r d t 2 + d r 2 1 − 2 M r + r 2 d θ 2 + sin 2 ⁡ θ d ϕ 2 .

The Lagrangian for the particle L ̄ = 1 2 g μ ν ( x ) x ̇ μ x ̇ ν reads L ̄ = 1 2 g μ ν x x ̇ μ x ̇ ν = 1 2 − 1 − 2 M r t ̇ 2 + r ̇ 2 1 − 2 M r + r 2 ϕ ̇ 2 .

We are only considering geodesics on the equatorial plane, where θ = π 2 . Hence, sin θ = 1 along these geodesics and since θ is fixed, every θ ̇ component is 0. Owing to the symmetry of the solution, we have two constants of motion which we can read off from the Euler–Lagrange equations. The t-component of the equation of motion gives 1 − 2 M r t ̇ = E = c o n s t a n t , while the ϕ component gives r 2 ϕ ̇ = L = c o n s t a n t .

In this article, we will mainly focus on the lensing of light rays 3 . Hence, we will be considering the light-like trajectories with the null ray condition defined as follows: g μ ν x ̇ μ x ̇ ν = 0 .

Hence, − 1 − 2 M r t ̇ 2 + r ̇ 2 1 − 2 M r + r 2 ϕ ̇ 2 = 0 .

Then, we write the r component of the equation of motion in terms of the ϕ component. To obtain that, we divide (2.3) by (2.4) and also define the impact parameter as λ = L E .

Then, we obtain d t d ϕ = t ̇ ϕ ̇ = r 2 λ 1 − 2 M r .

Then, using (2.6), we obtain d r d ϕ 2 = r 4 λ 2 − r 2 1 − 2 M r .

Eqs 2.8, 2.9 provide all the necessary information regarding the behavior of light-like geodesics. Considering Eq. 2.9 and taking the ϕ derivative of this equation, we obtain 2 d r d ϕ d 2 r d ϕ 2 = 4 r 3 λ 2 − 2 r + 2 M d r d ϕ .

For circular light-like geodesics, we must have d r d ϕ = 0 and d 2 r d ϕ 2 = 0 which gives r 4 λ 2 − r 2 1 − 2 M r = 0 , 4 r 3 λ 2 − 2 r + 2 M = 0 .

By eliminating b from the aforementioned equations, we obtain r = 3 M. We have, thus, shown that there is a circular light-like geodesic (or photon ring) r = 3 M. As we can choose any plane through the origin as our equatorial plane, there is actually a photon ring at this radius in the sense that every great circle on this sphere is a light-like geodesic. However, the photon rings at r = 3 M are unstable in the following sense: a light-like geodesic with an initial condition that deviates slightly from that of a photon ring at r = 3 M will spiral away from r = 3 M and either go to infinity or to the horizon.

Now that we have understood how the null geodesics behave in this geometry, it is time to focus on deriving actual observable quantities that one can measure. For this, we go on to study them in the upcoming sections one by one starting with calculating deflection angles.

2.1 Formula for the deflection angle

First, we set up the problem that we want to address. We consider a light ray that comes in from infinity and then goes through a minimum radius value at r = r min and then escapes back to infinity. Due to the geometry of spacetime around the central black hole, which for us is a Schwarzchild one, there will be a deflection. This deflection is simply because the rectilinear propagation of light will not be observed in this non-trivial geometry. The deflection angle measures the degree of this deviation from rectilinear propagation. In the following discussion, we will express this in terms of r min and the mass of the central object.

We start out with (2.7) and determine the ratio 1 λ 2 at r = r min . This gives 1 λ 2 = 1 r min 2 − 2 M r min 3 .

Hence, replacing this in (2.7), we obtain d ϕ = ± d r 1 r min 2 − 2 M r min 3 r 4 − r 2 + 2 M r , which, on integrating over the coordinates of the light ray, leads to π + α ̂ = 2 ∫ r min ∞ r min d r 1 − 2 M r min r 4 − r min 2 r 2 + 2 M r min 2 r = 2 ∫ r min ∞ r min d r r 2 r 2 − r min 2 − 2 M r min r 3 − r min 3 = 2 ∫ r min ∞ r min d r 1 − 2 M r min r 3 − r min 3 r r 2 − r min 2 r r 2 − r min 2 = 2 ∫ r min ∞ 1 + 1 2 2 M r min r 3 − r min 3 r r 2 − r min 2 + O 2 M r min 2 r min d r r r 2 − r min 2 = 2 ∫ r min ∞ r min d r r r 2 − r min 2 + 2 M r min ∫ r min ∞ r 3 − r min 3 r min d r r 2 r 2 − r min 2 3 / 2 + O 2 M r min 2 = π + 4 M r min + O 2 M r min 2 .

By neglecting higher-order terms, we obtain α ̂ = 4 M r min .

Now, we will look at a point or two about this derivation:

• From the derivation, it is clear that the integrand has a singularity at the lower bound r = r min . A more detailed analysis shows that the integral is finite for all values of r = r min that are bigger than 3 2 r s , where r s = 2 M. If we consider a sequence of light rays with r = r min approaching 3 2 r s from earlier, the deflection angle δ becomes bigger and bigger, which means that the light rays make more and more turns around the center. In the limit r min → 3 2 r s , the integral goes to infinity, and the limiting light ray spirals asymptotically toward a circle at r = 3 2 r s . If an unstable photon ring is approached, the deflection angle goes to infinity, and the singularity is logarithmic.

• The second point comes during the discussion on taking the Taylor series expansion in r s r min in the fourth step. The asymptotic behavior of light rays for r min approaching 3 2 r s is relevant only for black holes and for ultracompact stars. For an ordinary star, like our Sun, r min takes a value much bigger than r s . It is under this consideration that we take a Taylor series expansion in r s r min .

2.2 Shapiro time

Combining Eqs. 2.8 and 2.9 allows us to calculate the time taken by the light ray in this particular geometry. Here, we consider a light ray that starts at a radius r L , passes through a minimum of radius r min , and terminates at a radius r o . From Eqs. 2.8 and 2.9, we obtain d r d t 2 = d r d ϕ 2 d ϕ d t 2 = r 4 λ 2 − r 2 1 − 2 M r 1 − 2 M r 2 λ 2 r 4 = r 3 λ 2 − r 1 − 2 M r r − 2 M 2 λ 2 r 5 .

At r = r min , (2.16) reads r min 4 λ 2 − r min 2 1 − 2 M r min = 0 ⇒ 1 λ 2 = r min − 2 M r min 3 .

Finally, (2.16) becomes d r d t 2 = r min − 2 M r 3 r min 3 − r 1 − 2 M r r − 2 M r min 3 r min − 2 M r 5 , which gives d t = ± r min − 2 M r 5 / 2 d r r − 2 M r min 3 / 2 r min − 2 M r 3 r min 3 − r 1 − 2 M r .

By integrating, we obtain the travel time for the light ray, Δ t = − ∫ r L r min + ∫ r min r O r min − 2 M r 5 / 2 d r r − 2 M r min 3 / 2 r min − 2 M r 3 r min 3 − r 1 − 2 M r , where the signs on the right-hand side had to be chosen in such a way that the time coordinate is always increasing along the light ray. One can perform this integral exactly in terms of an elliptic integral. However, when r min ≫ r s , we can make a Taylor approximation, in exactly the same way as we did it for the deflection formula, and obtain Δ t = ∫ r min r L + ∫ r min r O r d r r 2 − r min 2 + 2 M ∫ r min r L + ∫ r min r O d r r 2 − r min 2 + 2 M r min 2 ∫ r min r L + ∫ r min r O d r r − r min r + r min 3 / 2 + ⋯

The zeroth-order term is, of course, the Euclidean travel time for a light ray with speed c along a straight line. The deviation of the general-relativistic calculation from this zeroth-order term is known as the Shapiro time delay [142].

I. Shapiro suggested using this effect as the fourth test of general relativity (after perihelion precession, light deflection, and gravitational redshift). In the first experiment, a strong radio signal was sent to Venus when it was in opposition to Earth, and the time was measured until the signal arrived back on Earth after being reflected in Venus’s atmosphere. Later experiments were performed with transponders on spacecraft, which sent the signal back with increased intensity. The best measurement to date was performed with the Cassini spacecraft in 2002. The general-relativistic time delay was verified to be within an accuracy of 0.001% [142].

 2.3 Angular radius of Einstein rings

An Einstein ring can be observed when the light source and observer are perfectly aligned (directly opposite to each other). We want to determine the angular radius θ E of the Einstein ring in dependence on the radius coordinate r L of the light source, the radius coordinate r O of the observer, and the Schwarzschild radius (= 2M). We use the formula d ϕ = ± d r r min − 2 M c 2 r 4 r min 3 − r 2 + 2 M r .

Integrating over the light ray gives π = ∫ r min r L + ∫ r min r O d r r min − 2 M c 2 r 4 r min 3 − r 2 + 2 M r .

The integration gives r min = f r L , r O , 2 M .

With r min determined, the angular radius of the Einstein ring can be obtained by tan θ E = r d ϕ 1 − 2 M r − 1 / 2 d r | r = r O .

3 Lensing for the Kerr–Newman black hole

In Section 2, we discussed the gravitational lensing in Schwarzschild spacetime. We expand the integrand around the turning point to obtain a simplified result, but in principle, we can calculate the exact deflection angle. In this section, first and foremost, we will generalize the result of lensing for rotating spacetime. We will calculate the deflection angle for the light rays in Kerr–Newman spacetime which is an axisymmetric spacetime. One can reproduce the result for the Kerr and Schwarzschild case as special limits. We will also discuss the strong deflection angle’s analytical form. We will closely follow the notation in [120].

In Boyer–Lindquist, the line element is given by d s 2 = − Δ Σ d t − a sin 2 ⁡ θ d ϕ 2 + Σ d r 2 Δ + d θ 2 + sin 2 ⁡ θ Σ a d t − r 2 + a 2 d ϕ 2 , where Σ = r 2 + a 2 ⁡ cos 2 ⁡ θ , Δ r = Δ r = r 2 + a 2 − 2 m r + Q 2 .

In (3.3) and (3.2), m ≥ 0, Q, and a are, respectively, the ADM mass, charge, and angular momentum of the black hole. As we know, the particle Lagrangian, specifically for the photons, is given by (2.2), with x ̇ μ = d x μ d λ ̃ for some convenient parameter λ ̃ . The metric coefficients are independent of t and ϕ, implying that we will have two conserved quantities along the photon trajectory: energy E and angular momentum L. Using the constants of motion and the null ray condition in (2.5), the null geodesic equations can be written as follows: ρ 4 r ̇ 2 = r 2 + a 2 E − a L 2 − Δ r L − a E 2 + K ̄ ≔ R 2 r , Σ ϕ ̇ = − a E − L sin 2 ⁡ θ + a E r 2 + a 2 − a L Δ , Σ t ̇ = − a a E sin 2 ⁡ θ − L + r 2 + a 2 E r 2 + a 2 − a L Δ , ρ 4 θ ̇ 2 = K ̄ + cos 2 ⁡ θ a 2 E 2 − L 2 sin 2 ⁡ θ ≔ Θ θ 2 where the Carter constant, which results from the separability of the Hamilton–Jacobi equation, is K ̄ . For our purpose, we will look upon into the equatorial case only, that is, θ = π 2 , θ ̇ = 0 which automatically implies K ̄ = 0 , and then will investigate the non-equatorial case also. We can use the geodesic equations to calculate the equatorial deflection angle. Rather in the next section, we will provide the deflection angle calculations and reproduce the result for Kerr by taking the limit Q = 0.

3.1 Deflection angle for Kerr–Newman spacetime

From (3.3), we can write the radial geodesic equations as follows: r ̇ 2 L 2 + V r = 1 λ 2 , with the effective potential V r = 1 r 2 1 − a 2 λ 2 + 1 − a λ 2 − 2 m r + Q 2 r 2 , and λ is the impact parameter defined in (2.7).

Now, we think of light rays that originate at infinity, pass through the black hole, and then, return to infinity to reach the observer. The closest approach to the black hole, r 0, will be the radial turning point for these light rays, determined by r ̇ 2 L 2 r = r 0 = 1 λ 2 − V r 0 = 0 .

From (3.6), we obtain r 0 4 − λ 2 1 − a 2 λ 2 r 0 2 + 2 m λ 2 1 − a λ 2 r 0 = Q 2 λ 2 1 − a λ 2 .

Solving (3.7), we obtain r 0 λ = λ 6 1 − ω 2 1 + γ ̂ + 2 − γ ̂ − 3 6 m 1 − ω 2 λ 1 − ω 2 3 / 2 1 + γ ̂ , where ω , χ , ρ , γ ̂ are given by ω = a λ , ρ = 12 Q 2 λ 2 1 + ω 2 , γ ̂ = 1 − ρ cos 2 χ 3 , χ = arccos 3 3 m 1 − ω 2 λ 1 − ω 2 3 / 2 1 − ρ 3 / 4 1 − λ 2 1 + ω 3 54 m 2 1 − ω 1 + 3 ρ − 1 − ρ 3 / 2 .

One can convince that the non-zero charge of the black holes gives a repulsive effect on the light rays, but the effect of spin attracts the light rays toward the black hole.

 3.2 Photon sphere radius and critical impact parameter

In this subsection, we will define the photon sphere and the critical impact parameter which will be useful in the subsequent sections. The radius of the photon sphere is the value of r where the effective potential attains its maximum value. ∂ V r ∂ r r = r c = 0 .

From (3.10), we can find out r c λ c = 3 m ζ ̂ 2 1 + 1 − ζ , where ζ ̂ = 1 − a λ c 1 + a λ c , ζ = 8 Q 2 9 m 2 ζ ̂ .

It should be noted that the turning point of the photon is r 0 defined in (3.8) and attains its minimum value at r c with the corresponding impact parameter λ c . If for some λ, r 0 becomes less than r c , then the photon will fall into the black hole. λ c is called the critical impact parameter. Finally, substituting Eq. 3.8 into (3.11), we can write λ c = λ c (a, Q, l) and r c = r c (a, Q, l).

 3.3 An exact analytical computation of the deflection angle

In this subsection, we will compute the exact photon deflection angle near a Kerr–Newman black hole. To perform this, we can choose any polar plane, but to keep things simple, we choose the equatorial plane θ = π 2 and θ ̇ = 0 . Furthermore, to investigate the observational signature, we need to employ the strong deflection limit. To perform this, we have to evaluate the deflection angle in the mentioned plane and take the strong limit. We can rather perform the reverse procedure also; that is, from the beginning, we take the strong limit and then calculate the deflection angle at this limit. The second one is simpler to perform. We adopt this track while discussing the observational signature.

The analysis has been performed in [143, 144]. We review the analysis here. Before proceeding with the computation, we define the following coordinates: u ≔ 1 r .

By combining the first two equations in (3.3), we obtain d u d ϕ 2 = d u d r d r d ϕ 2 = r 2 r 6 r ̇ 2 ϕ ̇ 2 = u 4 r ̇ 2 ϕ ̇ 2 .

Now, our goal is to rewrite r ̇ and ϕ ̇ as a function of u. We have r ̇ 2 = L 2 1 λ 2 − V r = L 2 1 λ 2 − 1 r 2 1 − a 2 λ 2 + 1 − a λ 2 ( − 2 m r + Q 2 r 2 = L 2 1 λ 2 − u 2 1 − a 2 λ 2 − Q 2 1 − a λ 2 u 4 + 2 m 1 − a λ 2 u 3 ≔ L 2 B u .

Combining Eqs 3.14 and 3.3 and using Eq. 3.13, we obtain d u d ϕ 2 = 1 − 2 m u + a 2 + Q 2 u 2 1 − 2 m u − Q 2 u 2 1 − a λ 2 B u .

The photon deflection angle α ̂ can be calculated by integrating Eq. 3.15 over u from 0 to 1 r , where r 0 is the turning point, and then evaluating the resulting expression at the critical value r c [79], α ̂ = − π + 2 ∫ 0 1 r 0 d u 1 − 2 m u − Q 2 u 2 1 − a λ 1 − 2 m u + a 2 + Q 2 u 2 1 B u , with ω = a λ . This integral can be computed exactly. We give the final result here. The details of the computation of this integral are given in Appendix A. α ̂ = − π + 4 1 − ω Q 2 u 4 − u 2 u 3 − u 1 × G + + K Q + u 1 u + − u 1 Π n + , k − Π n + , ϕ , k + G − + K Q − u 1 u − − u 1 Π n + , k − Π n + , ϕ , k − G + + K Q + u 4 u + − u 4 Π n + , k − Π n + , ϕ , k − F π 2 , k + F ϕ , k − G − + K Q − u 4 u − − u 4 Π n − , k − Π n − , ϕ , k − F π 2 , k + F ϕ , k , where u 1 = X 1 − 2 m − X 2 4 m r 0 , u 2 = 1 r 0 , u 3 = X 1 − 2 m + X 2 4 m r 0 , u 4 = 2 m Q 2 − X 1 2 m r 0 , n ± = u 2 − u 1 u ± − u 1 1 + 2 m Q 2 1 − r 0 u ± 4 m 2 r 0 − Q 2 X 1 + 2 m , k 2 = X 2 + 6 m − X 1 8 m 2 r 0 − Q 2 X 2 − 2 m + 3 X 1 4 X 2 4 m 2 r 0 − Q 2 X 2 + 2 m , ψ 0 = arcsin X 2 + 2 m − X 1 4 m 2 r 0 − Q 2 X 2 + 2 m X 2 + 6 m − X 1 4 m 2 r 0 − Q 2 X 1 , G + = 2 m 1 − ω m + m 2 − a 2 + Q 2 − a 2 + Q 2 2 a 2 + Q 2 m 2 − a 2 + Q 2 , G − = a 2 + Q 2 − 2 m 1 − ω m − m 2 − a 2 + Q 2 2 a 2 + Q 2 m 2 − a 2 + Q 2 , G Q + = − Q 2 1 − ω m + m 2 − a 2 + Q 2 2 a 2 + Q 2 m 2 − a 2 + Q 2 , G Q − = Q 2 1 − ω m − m 2 − a 2 + Q 2 2 a 2 + Q 2 m 2 − a 2 + Q 2 .

It should be noted that Π(n +, ψ 0, k) and Π(n +, k) are the incomplete and complete elliptic integral of the third kind, respectively. Also, F(ψ 0, k) and F ( π 2 , k ) are the incomplete and complete elliptic integral of the first kind, respectively. Also, X 1 = 2 m Q 2 + 4 m r 0 3 Q 2 + 8 m 2 r 0 3 Q 2 1 + Q 2 2 m 2 m r 0 − 3 1 + ω 2 1 − ω − Q 2 r 0 2 cos δ 3 + 2 π 3 , where δ = arccos − 8 m 3 r 0 3 − 3 m Q 2 r 0 2 2 m − 3 r 0 1 + ω 1 − ω − 3 Q 4 r 0 5 m − 3 r 0 1 + ω 1 − ω + 10 Q 6 4 m 2 r 0 2 + Q 2 r 0 2 m − 3 r 0 1 + ω 1 − ω − 2 Q 4 3 / 2 , and X 2 can be obtained from the following equations after inserting X 1 from (3.20). X 2 2 − X 1 2 − 2 m 2 1 8 m r 0 3 − Q 2 X 1 32 m 3 r 0 4 = 1 λ 2 1 − ω 2 .

So, the exact deflection angle can be derived as discussed. One can reproduce the results for the Schwarzschild black hole in the limit a → 0 and Q → 0 and the Kerr black hole in the limit a → 0.

 3.4 Strong deflection analysis for axisymmetric spacetime

In this section, we inspect the strong limit of the equatorial deflection angle ( θ = π 2 ) of the light rays discussed in the previous section in order to gain further understanding of the deflection angle for our context and make touch with the feasible observational signature. The strong deflection limit of (4.14) is challenging to take directly. It will be simpler to perform the integration after taking the strong-field limit of the integrand of (4.4). We shall only take into account the photons having the turning point very close to the photon sphere’s radius, as was previously specified.

The metric on the equatorial plane has the following structure: d s 2 = − A ̄ r d t 2 + B ̄ r d r 2 + C ̄ r d ϕ 2 − D ̄ r d t d ϕ , with A ̄ r = Δ r − a 2 Σ , B ̄ r = Σ Δ , C ̄ r = 1 Σ r 2 + a 2 + l 2 2 − Δ r a 2 , D ̄ r = 2 Σ r 2 + a 2 + l 2 a − Δ r a .

It should be noted that all metric components are evaluated in the equatorial plane. We have already seen that the spacetime admits two conserved quantities E and L due to the existing symmetries. To keep things simple, we set E = 1. As a result, the impact parameter is λ = L. Using the fact that at the distance of closest approach, r = r 0 and r ̇ = 0 , we obtain the following results from (3.4): L = − D ̄ 0 + 4 A ̄ 0 C ̄ 0 + D ̄ 0 2 2 A ̄ 0 = r 0 r 0 2 a 2 − 2 m r 0 + Q 2 + r 0 2 + a Q 2 − 2 m r 0 Q 2 r 0 + l 2 r 0 − 2 m + r 0 2 r 0 − 2 m .

The subscript “0” denotes that functions are evaluated at r = r 0. From the equation of motion of ϕ [the second equation of (3.3)], we obtain ϕ r 0 = 2 ∫ r 0 ∞ B ̄ A ̄ 0 D ̄ + 2 L A ̄ 4 A ̄ C + D ̄ 2 C ̄ A 0 − A ̄ C 0 + L A ̄ D 0 − D ̄ A 0 d r .

In the strong limit, we only consider the photons having closest approach r c near to the radius of the photon sphere. To implement the limit, one can expand the deflection angle α ̂ around r c or λ c and perform the radial integral function. When the turning point r 0 is greater than the radius of the photon sphere r c , then we obtain a finite deflection angle; otherwise, the photon will be caught by the black hole and the deflection angle diverges. Following the method developed here [76, 79], it is easy to find out the nature of the divergence in the deflection angle when the photons are very close to the radius of the photon sphere r = r c . One can define two variables y, z as z 1 = A ̄ r , z 2 = z 1 − z 1,0 1 − z 1,0 .

Now, the azimuthal angle defined in (3.25) can be expressed in terms of these two new variables, ϕ r 0 = ∫ 0 1 R ̄ z 2 , r 0 F ̄ z 2 , r 0 d z 2 , where R ̄ z 2 , r 0 = 2 1 − z 1,0 A ′ B ̄ A ̄ 0 D ̄ + 2 L A ̄ 4 A ̄ C 2 + C ̄ D 2 , F ̄ z 2 , r 0 = 1 1 C ̄ C ̄ A 0 − A ̄ C 0 + L A ̄ D 0 − D ̄ A 0 = 1 H ̄ , H ̄ = 1 C ̄ C ̄ A 0 − A ̄ C 0 + L A ̄ D 0 − D ̄ A 0 .

The function R ̄ ( z , r 0 ) is regular for any z and r 0 values, whereas the function F ̄ ( z , r 0 ) is divergent for z = 0, that is, at r = r 0. As a result, after extracting the divergent part, one can rewrite (3.25). ϕ r 0 = ϕ R ̄ r 0 + ϕ F ̄ r 0 , where the divergent part can be written as follows: ϕ F ̄ r 0 = ∫ 0 1 R ̄ z 2 = 0 , r c F ̄ 0 z 2 , r 0 d z .

We know that the deflection angle should diverge at r 0 = r c , indicating that the photon has been captured by the black hole. Next, we want to determine the nature of the divergence. Examining the denominator will allow us to determine the nature of the divergence (4.18). In order to achieve this, we Taylor expand the denominator of F ̄ 0 ( r 0 , z 2 ) (3.33) around z 2 = 0. F ̄ 0 z 2 , r 0 ≈ 1 σ 1 r 0 z + σ 2 r 0 z 2 2 + O z 3 .

It is worth noting that if σ 1(r 0) = 0 (this occurs when r 0 coincides with the radius of the photon sphere [145]), then it is clear from (3.33) that the leading term is 1 z in the small z limit. As a result, after integration, we obtain a logarithmic divergence, as shown in (3.42). To find σ 1 and σ 2, we first Taylor expand H ̄ , which is defined in (4.21). H ̄ z , r 0 = H ̄ 0 , r 0 + ∂ H ̄ ∂ z 2 z 2 = 0 z 2 + 1 2 ! ∂ 2 H ̄ ∂ z 2 2 z 2 = 0 z 2 2 + O z 2 3 , with H ̄ 0 , r 0 = 0 .

Therefore, using (4.20), we can identify σ 1 and σ 2 as σ 1 ≔ ∂ H ̄ ∂ z 2 z 2 = 0 = 1 − A ̄ 0 A ̄ 0 ′ C ̄ 0 A ̄ 0 C ̄ 0 ′ − A ̄ 0 ′ C ̄ 0 − L A ̄ 0 D ̄ 0 ′ − A ̄ 0 ′ D ̄ 0 , and σ 2 ≔ 1 2 ! ∂ 2 H ̄ ∂ z 2 2 z 2 = 0 = 1 − A ̄ 0 2 2 C ̄ 0 2 A ̄ ′ 0 3 × 2 C ̄ 0 C ̄ 0 ′ A ̄ ′ 0 2 + C ̄ 0 C ̄ 0 ″ − 2 C ̄ ′ 0 2 A ̄ 0 A ̄ 0 ′ − C ̄ 0 C ̄ 0 ′ A ̄ 0 A ̄ 0 ″ + L A ̄ 0 C ̄ 0 A ̄ ″ 0 D ̄ 0 ′ − A ̄ 0 ′ D ̄ ″ 0 + 2 A ̄ 0 ′ C ̄ 0 ′ A ̄ 0 D ̄ 0 ′ − A ̄ 0 ′ D ̄ 0 ] .

We can write the regular part as follows: ϕ R ̄ r 0 = ∫ 0 1 G ̄ z 2 , r 0 d z 2 , where G ̄ ( z 2 , r 0 ) = R ̄ ( z 2 , r 0 ) F ̄ ( z 2 , r 0 ) − R ̄ ( z 2 = 0 , r c ) F ̄ 0 ( z 2 , r 0 ) . The coefficient σ 1 = 0, in the strong deflection limit. This implies A ̄ 0 C ̄ 0 ′ − A ̄ 0 ′ C ̄ 0 − L A ̄ 0 D ̄ 0 ′ − A ̄ 0 ′ D ̄ 0 r 0 = r c = 0 .

From (3.38), we obtain − 2 a Q 2 − m r c a 2 − 2 m r c + Q 2 + r c 2 + 2 a 2 Q 2 − m r c + r c 2 6 m 2 − 5 m r c + r c 2 + Q 2 r c ( 3 r c − 7 m ) + 2 Q 4 = 0

It is simple to verify that the expression for the critical impact parameter λ c found in (3.24) yields the identical expression for the radius of the photon sphere found in (3.11). The photon sphere is defined yet again by (4.25), according to this. Readers who are interested in learning more about the geometry of photon spheres and several complementary definitions of photon spheres are directed to [145].

For fixed values of Q, a, we can compute the radius of the photon sphere in Kerr–Newman spacetime. In the next section, we will discuss the more general case which is the Kerr–Newman black-bounce spacetime, and the solutions of the photon sphere equation will be discussed there in a more general setting.

We may now assess the divergent integral (4.23) with the help of the following equation: ϕ F ̄ r 0 ≈ r c = R ̄ z 2 = 0 , r c ∫ 0 1 1 σ 1 z 2 + σ 2 z 2 2 d z 2 = R ̄ z 2 = 0 , r 0 ≈ r c 2 σ 2 log σ 2 + σ 2 + σ 1 σ 1 .

We are aware that the function ϕ F ̄ ( r 0 ) diverges when r 0 = r c , or when the coefficient σ 1 = 0. In order to obtain the nature of ϕ F ̄ ( r 0 ≈ r c ) , the idea is to expand the σ 1(r 0) around r 0 = r c up to the first order and substitute it into (3.40). σ 1 r 0 = ∂ σ 1 ∂ r 0 r 0 = r c r 0 − r c + O r 0 − r c 2 , σ 2 r 0 = σ 2 r c + ∂ σ 2 ∂ r 0 r 0 = r c r 0 − r c + O r 0 − r c 2 .

Substituting (3.41) into (3.40) and using condition (3.38), we obtain ϕ F ̄ r 0 ≈ r c = − a ̃ log r 0 r c − 1 + b ̃ + O r 0 − r c .

Alternatively, Eq. 3.42 can be written in terms of the impact parameter as [79] α ̂ λ = − a ̄ log λ λ c − 1 + b ̄ + O λ − λ c , where the coefficients are a ̄ = 2 A ̄ c B ̄ c A ̄ c C ̄ c ″ − A ̄ c ″ C ̄ c + λ c A ̄ c ″ D ̄ c − A ̄ c D ̄ c ″ , b ̄ = − π + b R + a ̄ log 4 σ 2 c C ̄ c λ c A ̄ c D ̄ c + 2 λ c A ̄ c , λ c = L c . In the next section, we will show how the deflection angle (3.43) varies with respect to the impact parameter. The black hole can be viewed as a lens, with its gravitational field curving the path of photons. Let θ = λ D O L be the angular separation between the image and lens, and D 0l be the distance between the lens and the observer. The deflection angle (3.43) can then be expressed as follows: α ̂ θ = − a ̄ log θ D o l λ c − 1 + b ̄ + O λ − λ c , where b R = ϕ R ̄ ( r c ) = ∫ 0 1 G ̄ ( z , r c ) d z .

In order to perform the integral function, we can expand R ̄ ( z , r c ) around z = 0 and then put it in (4.23). Formally, we obtain the following expression: b R = 1 σ 2 ∫ 0 1 R ̄ z z − R ̄ 0 z d z r c = 1 σ 2 ∑ n = 1 ∞ 1 n ∂ n R ̄ ∂ z n z = 0 →  finite .

Now, we will have multiple images of the source if the deflection angle is greater than 2π. The angular radius of the Einstein ring, which is created due to the symmetric lensing of light rays coming from some distant source, can be calculated using the strong-field deflection angle formula given in (3.47) and the lens equation [146, (147). The observational signatures of the deflection angle will be discussed in the following section.

 3.5 Non-equatorial lensing

So far, we have discussed equatorial lensing. In this section, we will briefly discuss about the non-equatorial lensing in Kerr–Newmann spacetime for small inclination. We assume that the inclination is θ = π 2 − ψ , with ψ being very small. To conduct the analysis, we will closely follow the analysis presented in [79].

For the non-equatorial plane, the carter constant is K ̄ ≠ 0 . Instead, for the small inward inclination angle, we can write down the constants in terms of the inclination angle ψ as follows: L ≈ λ , K ̄ ≈ h 2 + λ 2 − a 2 ϕ 2 , with ϕ ≈ h λ .

In principle, one can parameterize the light ray coming from infinity by three parameters (ϕ, h, λ). If there is no gravitational field, the projection of the photon line on the equatorial plane has a minimum distance from the origin which is λ. Now, for given λ, the vertical distance of the light ray from the plane is h, and finally, ϕ is the inclination angle formed by the light ray with the equatorial plane.

Now, using the θ and ϕ geodesic equation in (3.3) and requiring ψ to be small, we will have d ψ d ϕ = ω r ϕ ψ ̂ 2 − ψ 2 , with ψ ̂ = h 2 λ 2 − a 2 + ϕ 2 .

Now, we are interested in computing the deflection angle. For that, we first write down the following equation: ϕ ̄ f = ∫ 0 ϕ f d ϕ ω ϕ , where ϕ f is the total azimuthal shift. Then, the deflection angle can be written as follows: α ̂ = − π + 2 ∫ r 0 ∞ d r ω r d ϕ d r = ∫ 0 1 d z ω r z R ̄ z , r 0 F ̄ z , r 0 , where ω r = λ ̄ a 2 + r r − 2 r 2 a + λ r − 2 + Q 2 λ − a , with λ ̄ = λ 2 − a 2 . We can follow the same method as mentioned in Section 3.4 to extract the divergent part as well as the finite part of the deflection angle as the function ω(r). They are given by α ̂ = − a ̂ log λ λ c − 1 + b ̂ , where a ̂ = ω z = 0 , r c R ̄ z = 0 , r c 2 σ 2 c , b ̂ = − π + b ̂ R + a ̂ log 4 σ 2 c C ̄ c λ c A ̄ c D ̄ c + 2 λ c A ̄ c , and b ̂ R = ∫ 0 1 d z ω z , r c R ̄ z , r c F ̄ z , r c − ω z = 0 , r c R ̄ z = 0 , r c F ̄ z , r c .

Our interest is to find the position of the caustics, where the magnification diverges, and is given in [79], γ ̄ k = − b ̄ + a ̄ ( b ̂ − k π ) . Here, k is a positive integer. We have one caustic for the direct photon and one for the retrograde photon for each value of k. We will find the caustic point in the weak-field limit for k = 1 and the strong-field limit for k ≥ 1, which is the regime of interest to us. In the next section, we give the results for caustic points in a more generic setup, that is, for the black-bounce metric. By taking the limit, l = 0, one can reproduce the results for the Kerr–Newman case. Furthermore from Eq. 3.54, one can see ω(r) = 1 implying that for Schwarzschild metric, there is no such difference between the equatorial and non-equatorial case.

 4 Lensing for Kerr–Newman black-bounce spacetime

Before proceeding further, we will discuss how to calculate the exact deflection angle for Kerr–Newman black-bounce spacetime [137] and then go to the observational signatures gradually. It is interesting to study because it has one more parameter (apart from mass, rotation, and charge) that regularizes the central singularity. One can reproduce the results for Schwarzschild, Kerr, and Kerr–Newman spacetimes by taking appropriate limits. We start by applying the general formalism of lensing in this special kind of axisymmetric non-singular spacetime. We will closely follow the notation of [120] throughout this section.

4.1 Brief review of Kerr–Newman black-bounce spacetime

We will begin with a brief discussion of the null geodesics in Kerr–Newman black-bounce spacetime. To perform that, we first write the corresponding metric in the Boyer–Lindquist coordinate [137]. d s 2 = g μ ν d x μ d x ν = − Δ Σ a sin 2 ⁡ θ d ϕ − d t 2 + sin 2 ⁡ θ Σ r 2 + a 2 + l 2 d ϕ − a d t 2 + Σ d r 2 Δ + d θ 2 , where Σ = r 2 + l 2 + a 2 ⁡ cos 2 ⁡ θ , Δ r = r 2 + a 2 + l 2 − 2 m r 2 + l 2 + Q 2 .

Here, m ≥ 0 is the ADM mass, Q is the black hole charge parameter, and a = J m corresponds to the angular momentum per mass. The non-vanishing regularizing parameter l > 0 accounts for the absence of the central singularity. Keeping in mind that the radial coordinate’s range in this instance is −∞ < r < ∞, in the Q = 0, l = 0 limit, we recover the Kerr metric. The location of the event horizon can be determined by equating Δ(r) = 0. R H = m + m 2 − a 2 + Q 2 2 − l 2 .

We also need to impose the reality condition. That gives m 2 − a 2 + Q 2 > 0 a n d m + m 2 − a 2 + Q 2 > l .

Following the procedure mentioned in Section 3, one can obtain the turning point and radius of the photon sphere.

 4.2 Perturbative computation of the deflection angle: Analytical results

Following the analysis mentioned in Section 3.3, one can write down the integral form of the deflection angle which is given by α ̂ = − π + 2 ∫ 0 1 r 0 2 + l 2 1 1 − l 2 u 2 1 − 2 m u − Q 2 u 2 1 − a λ 1 − 2 m u + a 2 + Q 2 u 2 1 B u .

The polynomial B(u)(1 − l 2 u 2) in (4.4) has degree six. As a result, we cannot write this integral as an elliptic integral in its entirety. However, in order to make some analytical headway, we will make the following assumption: l 2 u 2 ≪ 1 .

Then, we can Taylor expand 1 1 − l 2 u 2 = 1 + l 2 u 2 2 + O l 4 u 4 .

Finally, keeping terms upto O ( l 4 ) in (4.4), we obtain α ̂ = − π + 2 ∫ 0 1 r 0 2 + l 2 d u 1 + l 2 u 2 2 1 − 2 m u − Q 2 u 2 1 − ω 1 − 2 m u + a 2 + Q 2 u 2 1 B u + O l 4 = α ̂ K N r 0 2 + l 2 + l 2 ξ m , a , λ , Q + O l 4 , where α ̂ K N r 0 2 + l 2 = − π + 2 ∫ 0 1 r 0 2 + l 2 d u 1 − 2 m u − Q 2 u 2 1 − ω 1 − 2 m u + a 2 + Q 2 u 2 1 B u .

In the l = 0 limit, α ̂ KN reduces to the deflection angle for the Kerr–Newman black hole as mentioned in (3.16) and ξ m , a , λ , Q = ∫ 0 1 r 0 2 + l 2 u 2 1 − 2 m u − Q 2 u 2 1 − ω 1 − 2 m u + a 2 + Q 2 u 2 1 B u d u , with ω = a λ . We now rewrite B(u) as B u = − Q 2 1 − w 2 u − u 1 u − u 2 u − u 3 u − u 4 , where the roots are defined as u 1 = X 1 − 2 m − X 2 4 m r 0 2 + l 2 , u 2 = 1 r 0 2 + l 2 , u 3 = X 1 − 2 m + X 2 4 m r 0 2 + l 2 , u 4 = 2 m Q 2 − X 1 2 m r 0 2 + l 2 .

Again, we apply the same strategy as sketched in Appendix A for the Kerr–Newman case. We choose the constants X 1 and X 2 in such a way so that we can write down the roots in the following order u 1 < u 2 < u 3 < u 4. Similar to the Kerr–Newman case, u 2, u 3, u 4 turn out to be the positive roots, while u 1 turns to be a negative root. To find out the roots, we need to substitute Eqs 4.9.12.–.4.4.12 into (4.8) and compare it with the coefficients of u 0, u 2, u 3, u 4 in (3.14), and then, we will obtain [144] Q 2 X 2 2 − X 1 − 2 m X 1 + 6 m + 4 X 1 2 = 16 m 2 r 0 2 + l 2 X 1 − r 0 2 + l 2 1 + ω 1 − ω , X 2 2 − X 1 2 − 2 m 2 = 8 m X 1 − 2 m Q 2 X 1 − 4 m 2 r 0 2 + l 2 Q 2 X 1 − 2 m − 4 m 2 r 0 2 + l 2 , X 2 2 − X 1 2 − 2 m 2 1 8 m r 0 2 + l 2 3 2 − Q 2 X 1 32 m 3 r 0 2 + l 2 2 = 1 λ 2 1 − ω 2 .

By exactly following the same procedure as in Section 3.3, one can write down the expression of deflection angle upto O ( l 2 ) , and the correction term is given by ξ m , a , λ , Q = ∫ 0 u 2 d u u 2 G + u + − u + G − u − − u + G Q + u u + − u + G Q − u u − − u 1 − Q 2 1 − ω 2 u − u 1 u − u 2 u − u 3 u − u 4 = G + g u + Δ F + u 1 − u 4 Δ Π α 2 + u 4 Δ F − u + 2 1 u + − u 1 u 1 − u 4 u + − u 4 Δ Π α + 3 2 + u + − u 1 u + − u 4 Δ F + G − g u − Δ F + u 1 − u 4 Δ Π α 2 ) + u 4 Δ F − u − 2 1 u − − u 1 u 1 − u 4 u − − u 4 Δ Π α − 3 2 + u − − u 1 u − − u 4 Δ F + G Q + g u 4 2 − 2 u + u 4 + u + 2 − u + 3 u + − u 4 Δ F + u + u 1 − u 4 − u + 3 u 1 − u 4 u + − u 1 u + − u 4 Δ Π α + 3 2 − 2 u 4 u 1 − u 4 χ 1 α 2 − u 1 − u 4 2 χ 2 α 2 + G Q − g u 4 2 − 2 u − u 4 + u − 2 − u − 3 u − − u 4 Δ F + u − u 1 − u 4 − u − 3 u 1 − u 4 u − − u 1 u − − u 4 Δ Π α − 3 2 − 2 u 4 u 1 − u 4 χ 1 α 2 − u 1 − u 4 2 χ 2 α 2 , where Δ F = F ψ 0 , k − F π 2 , k , Δ Π α 2 = Π ψ 0 , α 2 , k − Π π 2 , α 2 , k , Δ Π α ± 3 2 = Π ψ 0 , α ± 3 2 , k − Π π 2 , α ± 3 2 , k and χ 1 and χ 2 are given by χ 1 α 2 = Π π 2 , α 2 , k − Π ψ 0 , α 2 , k , χ 2 α 2 = 1 2 α 2 − 1 k 2 − α 2 α 2 E π 2 , k − E ψ 0 , k + k 2 − α 2 F π 2 , k − F ψ 0 , k , 2 α 2 k 2 + 2 α 2 − α 4 − 3 k 2 Π π 2 , α 2 , k − Π ψ 0 , α 2 , k − α 4 sin ψ 0 1 − sin 2 ψ 0 1 − k 2 sin 2 ψ 0 1 − α 2 sin 2 ψ 0 , and g = 2 Q 2 1 − ω 2 u 4 − u 2 u 3 − u 1 , α 2 = u 1 − u 2 u 4 − u 2 < 0 , α ± 3 2 = α 2 u ± − u 4 u ± − u 1 , k 2 = u 4 − u 3 u 2 − u 1 u 4 − u 2 u 3 − u 1 , ψ 0 = arcsin X 2 + 2 m − X 1 4 m 2 r 0 2 + l 2 − Q 2 X 1 + 2 m X 2 + 6 m − X 1 4 m 2 r 0 2 + l 2 − Q 2 X 1 .

Π(ψ 0, α 2, k) and Π ( π 2 , α 2 , k ) are the incomplete and complete elliptic integrals of the third kind, respectively. On the other hand, F(ψ 0, k) and F ( π 2 , k ) are the incomplete and complete first-kind elliptic integrals, respectively. In addition, we have E(ψ 0, k) and E ( π 2 , k ) in V 1,2, which are incomplete and complete elliptic integrals of the second kind, respectively. The computation details are given in [120]. So, similar to the Kerr–Newman case, we provide an analytical expression of the equatorial deflection angle for the Kerr–Newman black-bounce metric upto O ( l 2 ) in this section.

 4.3 Strong deflection analysis

Following the same analysis given in Section 3.4, we can find out the deflection angle of equatorial light rays in the strong deflection limit. We will summarize the steps and results as follows:

• First, we write down the deflection angle as follows:

ϕ r 0 = ∫ 0 1 R ̄ z 2 , r 0 F ̄ z 2 , r 0 d z 2 , with R ̄ z 2 , r 0 = 2 1 − y 0 A ′ B ̄ A ̄ 0 D ̄ + 2 L A ̄ 4 A ̄ C 2 + C ̄ D 2 , F ̄ z 2 , r 0 = 1 1 C ̄ C ̄ A 0 − A ̄ C 0 + L A ̄ D 0 − D ̄ A 0 = 1 H ̄ , H ̄ = 1 C ̄ C ̄ A 0 − A ̄ C 0 + L A ̄ D 0 − D ̄ A 0 .

The integral is potentially divergent at r 0 = r c .

• Second, we separate the convergent and the divergent integral as

ϕ r 0 = ϕ R ̄ r 0 + ϕ F ̄ r 0 , where ϕ F ̄ r 0 = ∫ 0 1 R ̄ z 2 = 0 , r c F ̄ 0 z 2 , r 0 d z 2 . →  the divergent part  ϕ R ̄ r 0 = ∫ 0 1 R ̄ z 2 , r 0 F ̄ z 2 , r 0 d z 2 − ∫ 0 1 R ̄ z 2 = 0 , r c F ̄ 0 z 2 , r 0 d z 2 . →  the finite part  .

• Next, we find out the critical turning point r c by solving the following equation:

 A ̄ 0 C ̄ 0 ′ − A ̄ 0 ′ C ̄ 0 − L A ̄ 0 D ̄ 0 ′ − A ̄ 0 ′ D ̄ 0 r 0 = r c = 0 , which gives 2 a 2 Q 2 − m l 2 + r c 2 − 2 a Q 2 − m l 2 + r c 2 a 2 − 2 m l 2 + r c 2 + l 2 + Q 2 + r c 2 + l 4 + l 2 − 5 m l 2 + r c 2 + 6 m 2 + 3 Q 2 + 2 r c 2 − 7 m Q 2 l 2 + r c 2 − 5 m r c 2 l 2 + r c 2 + 6 m 2 r c 2 + 2 Q 4 + 3 Q 2 r c 2 + r c 4 = 0 We reproduce the results of [120] for the dependence of the critical turning point on different spacetime parameters (a, Q, l) in Figure 1.

• At the end, we calculate the integrals in (4.23) at the limit r 0 → r c , and we will find the logarithmic nature of the deflection angle as shown in Figure 2. Again, we reproduce the result of [120] here.

Variation of radii of the photon sphere for charged, rotating Kerr–Newman black-bounce metrics for various values of (a, Q, l). In the leftmost figure, we vary l keeping a and Q fixed. In the middle figure, we vary a keeping l and Q, and in the rightmost figure, we vary Q keeping a and l [120].

Variation of the deflection angle α ̂ with respect to the impact parameter λ for fixed values of Q, l, a [120].

 4.4 Observational signature in the strong deflection limit

Now, we will discuss some observational consequences. The first step for this is to relate the deflection angle and the angular radius of the Einstein ring. This is performed by using a lens equation. In this paper, we will use the following lens equation [148]: β = θ − D L S D O S Δ α n , where D O S = D O L + D L S .

Also, β is the angular separation between the source and the lens. D LS , D OS , D OL are the distances between the lens to the source, the observer to the source, and the observer to the lens, respectively. Finally, Δ α n = α ̂ θ − 2 n π .

The angular separation between the lens and the n th image can be written as follows: θ n = θ n 0 + Δ θ n , where θ n 0 is the angular separation between the n image and the lens when the extra deflection angle (Δθ n ) over 2nπ is negligible.

For the perfect alignment, that is, when β = 0 and assuming Δ θ n ≪ θ n 0 , the angular separation (angular radius) can be written as [79] θ n Einstein = λ c D O L 1 + exp b ̄ a ̄ − 2 n π a ̄ .

n = 1 corresponds to the outermost Einstein ring. Following [120], we plot some of these rings for different values of Q, a and l in Figure 3.

Polar plots showing the angular radius θ 1 for different Q, l, a values [120].

From the leftmost plot in Figure 3, we can conclude that as the charge of the black hole Q increases (for fixed l and a), the radii of the ring decrease. On the other hand, the effect of changing l (for fixed a and Q) on the ring radius is negligible. This is evident from the rightmost plot in Figure 3.

To make this observation more concrete, we carry out a detailed study, as shown in Tables 1, 2. Some of the values provided in (1) and (2) are reproduced from [120]. We have listed some values regarding the representative percentage change in the angular radius of the outermost Einstein rings with respect to Q (for fixed a and l) and l (for fixed a and Q) in Tables 1, 2, respectively. These values corroborate perfectly the conclusion drawn previously.

Percentage change in the angular radius of the first Einstein ring for different values of charge Q for fixed a and l. Some of the numerical values presented here are reproduced from [120].

Angular separation Percentage change Values of charge Fixed parameter
θ 1 ( 1 ) ( Q 2 ) θ 1 ( 2 ) ( Q 1 ) Δ θ = ( θ 1 ( 2 ) − θ 1 ( 1 ) ) θ 1 ( 1 ) × 100 (Q 1, Q 2) l a
0.1686 0.1970 16.9013 (0, 0.6) 0.7 0.8
0.18548 0.21590 16.399 (0.6,0.8) 0.4 0.5
0.22691 0.24236 6.8088 (0.45,0.65) 0.4 0.5

Percentage change in the angular radius of the first Einstein ring for different values of the regularization parameter l for a = 0.8 and Q = 0.6. Some of the numerical values presented here are reproduced from [120].

Angular separation Percentage change Values of the regularization parameter
θ 1 ( 1 ) ( l 1 ) θ 1 ( 2 ) ( l 2 ) Δ θ = ( θ 1 ( 2 ) − θ 1 ( 1 ) ) θ 1 ( 1 ) × 100 (l 1, l 2)
0.16856 0.1703 1.07059 (0.2,0.4)
0.16756 0.16760 0.023872 (0.25,0.75)
0.166678 0.166680 0.001199 (0.15,0.25)

Before closing this section, a few comments regarding possible avenues to constraints on the spacetime parameters should be in order, utilizing observational data. One of the ways is to look into the ratio of mass to distance, the mass being the mass of the central object (e.g., Sagittarius A*), which is around 4.4 × 106 M ⊙ and its distance being 8.5 kpc, the ratio turns out to be around 2.4734 × 10−11. One can use these data to provide the angular position of the relativistic images and the angular separation between the two Einstein rings. On the other hand, we can compute the angular separation of two successive Einstein rings from (4.27) for different values of n. Then, one can utilize it to build a parameter space for the spacetime parameters. Interested readers are referred to, for example, [107] for a more comprehensive discussion. In future detections, if one can better resolve the angular separation of various Einstein rings, we will get better constraints on the charge of the underlying black-bounce metric.

 4.5 Results for non-equatorial lensing: Caustic points

In Section 3.5, we discussed about the non-equatorial lensing for the Kerr–Newman black hole. For this case also, the analysis would be the same. The only difference is that the scaling factor ω(r) will be the function of the regularization parameter l apart from Q and a. ω r = λ ̄ a 2 + r 2 + l 2 r 2 + l 2 − 2 r 2 + l 2 2 a + λ r 2 + l 2 − 2 + Q 2 λ − a , with λ ̄ = λ 2 − a 2 .

Next, we show how the angular radius of the first Einstein ring depends on Q and l in Figure 4 following [120]. It is evident from Figure 4 that the dependence of the angular radius on the charge parameter (Q) of the black-bounce metric is significantly more than that on the regularization parameter (l) similar to the case of equatorial lensing.

Polar plots showing the angular radius θ 1 for different Q, l, a values [120].

In Table 3, we investigate the variation of the second caustic point with respect to Q for fixed l and a.

Angular position of the second caustic point for different values of Q. The first three entities are for direct photons, and the last three are for retrograde photons. Some of the numerical values presented here are reproduced from [120].

(l, a) Q σ ̄ 2
(0.4, 0.5) 0.55 −6.28
(0.4, 0.5) 0.65 −4.98
(0.4, 0.5) 0.75 −0.137
(0.4, −0.5) 0.5 −5.69
(0.4, −0.5) 0.6 −5.84
(0.4, −0.5) 0.7 −5.90

Before closing this section, we further investigate the variation of the second caustic point with respect to l for fixed Q and a. It is demonstrated in Table 4. Comparing the values in Tables 3, 4, we can conclude that the change in caustic points for different Q and different l is not that robust [120].

Angular position of the second caustic point for different values of l. The first three entities are for direct photons, and the last three are for retrograde photons. Some of the numerical values presented here are reproduced from [120].

(Q, a) l σ ̄ 2
(0.5, 0.5) 0.45 −6.646
(0.5, 0.5) 0.55 −6.702
(0.5, 0.5) 0.65 −6.79
(0.5, −0.5) 0.7 −5.92
(0.5, −0.5) 0.8 −5.98
(0.5, −0.5) 1.0 −6.20
5 Conclusion

As mentioned earlier, gravitational lensing studies provide an excellent tool to provide insight into the structure of spacetime itself. In light of these advantages, this review provides a brief tour of the analytic methods used to calculate observables which can be measured. First, we review some facts about lensing in Schwarzschild geometry and observe the following points:

• The geodesics in this geometry are studied. As expected, owing to the symmetry of the spacetime itself, we have two constants of motion. Looking closely into the geodesic structures, we can find the location of the photon rings around them. Not only that, one can see that these photon rings located at 3M are unstable, and one can seek out some non-trivial physics once you deviate infinitesimally from this range.

• After we have understood how geodesics behave in such a geometry, we can calculate quantities using the equations at hand. A comprehensive and self-explanatory calculation is provided which gives an estimate of the deflection angle in terms of the central massive object responsible for this deflection. The contribution of this central massive object in the formula is through the mass of the object. There is also a contribution from an r min term in the denominator which indirectly also depends on the metric structure around the central massive object. We have also listed down some salient features related to the calculation in the bullet points in the following (Eq. 2.15).

After giving a brief overview of the time delay suffered by these geodesics and also going on to calculate the diameter of the Einstein ring in this setup, we move on to a more general spacetime having extra rotation parameters. As expected, all the aforementioned observables will have a non-trivial rotation parameter-dependent term. The calculations are all in the strong deflection limit and analytical. We also consider a more general class of proposed solutions called black-bounce spacetimes, where there is an extra parameter involved as a deviation of the already known solutions in GR. We list the salient features of our findings as follows:

• We presented a method of calculating the deflection angle analytically by performing a perturbation in l. The results are given in terms of elliptic integrals of various kinds. Also, we have restricted ourselves to the equatorial plane while performing this analysis. We observe that for non-zero l, the value of the deflection angle for a fixed impact parameter decreases. Our calculation provides a general methodology to compute deflection angles analytically in a perturbative series. In future, it will be interesting to go beyond this small l expansion. This will require a thorough numerical analysis.

• Next, we study the strong deflection limit of the equatorial deflection angle. This has direct observational implications because it provides information about the Einstein rings. We can conclude that the effect of the charge (Q) on the size of the Einstein ring is much more pronounced than that of the regularization parameter (l) for a fixed value of spin parameter a. We discovered that decreasing the charge (Q) considerably increases the ring’s size. This observation remained the same even when we computed the ring’s radius for a small polar inclination.

• Furthermore, we extend our analysis for non-equatorial lensing. This enables us to compute the location of the caustic points. We again observed that the effect of the charge (Q) on the position of caustic points is more pronounced than the regularization parameter (l). Again, one can apply this method to different black hole spacetimes to obtain an exact analytical expression. The study of non-equatorial lensing presented in this review assumes a small inclination angle. Again, it will be interesting to go beyond this regime.

Another interesting direction we could not include in this review is the analysis of the structure of the shadow. Interested readers are referred to this review [149] for more details about this topic. Analysis of the shadow structure complements the analysis of the Einstein ring, which is presented here. It provides further constraints on the different black hole parameters and various theories of gravity. Interested readers are again referred to [149] for relevant references for this.

The analysis of the deflection angle presented in this review can straightforwardly be repeated to other black-bounce spacetimes, for example, [150]. Finally, there are several other avenues which have been pursued recently in the context of strong deflection of light rays. One such thing is the study of multilevel images. This helps one to predict how much resolution is required to distinguish between different Einstein rings. One important aspect is also to study the two-point correlation function of intensity fluctuations on the photon ring, which result from the photon traveling through several orbits around the central object. This plays a significant role from the perspective of image analysis. This black-bounce metric might be subject to these studies along the lines of [151]. These investigations will support our efforts to communicate with plausible astrophysical settings.

 Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

 Funding

Research of AC was supported by the Prime Minister’s Research Fellowship (PMRF-192002-1174) of Government of India. AB was supported by the Mathematical Research Impact Centric Support Grant (MTR/2021/000490) by the Department of Science and Technology Science and Engineering Research Board (India) and the Relevant Research Project grant (202011BRE03RP06633-BRNS) by the Board of Research in Nuclear Sciences (BRNS), Department of Atomic Energy (DAE), India.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

This list is by no means exhaustive. Interested readers are referred to this review [24] and citations there for more details.

Again, this list is by no means exhaustive. Interested readers are referred to the references and citations of these papers.

For deflection of the massive particle, one needs to consider time-like trajectories given by g μ ν x ̇ μ x ̇ ν = − 1 .

 References Dyson FW Eddington AS Davidson C . A determination of the deflection of light by the sun’s gravitational field, from observations made at the total eclipse of may 29, 1919. Phil Trans Roy Soc Lond A (1920) 220:291–333. Carroll SM . Lecture notes on general relativity (1997). gr-qc/9712019. Jackiw R Pi S-Y . Chern-Simons modification of general relativity. Phys Rev D (2003) 68:104012. 10.1103/physrevd.68.104012 Kanti P Mavromatos NE Rizos J Tamvakis K Winstanley E . Dilatonic black holes in higher curvature string gravity. Phys Rev D (1996) 54:5049–58. [hep-th/9511071]. 10.1103/physrevd.54.5049 Horndeski GW . Second-order scalar-tensor field equations in a four-dimensional space. Int J Theor Phys (1974) 10:363–84. 10.1007/bf01807638 Eling C Jacobson T Mattingly D . Einstein-Aether theory. In: Deserfest: A celebration of the life and works of stanley deser (2004). p. 163–79. 10, 2004.gr-qc/0410001. Moffat JW . Black holes in modified gravity (MOG). Eur Phys J C (2015) 75(4):175. [1412.5424]. 10.1140/epjc/s10052-015-3405-x Gao X . Unifying framework for scalar-tensor theories of gravity. Phys Rev D (2014) 90:081501. [1406.0822]. 10.1103/physrevd.90.081501 Gao X Yamaguchi M Yoshida D . Higher derivative scalar-tensor theory through a non-dynamical scalar field. JCAP (2019) 03:006. [1810.07434]. 10.1088/1475-7516/2019/03/006 Gross DJ Witten E . Superstring modifications of Einstein’s equations. Nucl Phys B (1986) 277:1–10. 10.1016/0550-3213(86)90429-3 Gross DJ Sloan JH . The quartic effective action for the heterotic string. Nucl Phys B (1987) 291:41–89. 10.1016/0550-3213(87)90465-2 de Roo M Suelmann H Wiedemann A . The Supersymmetric effective action of the heterotic string in ten-dimensions. Nucl Phys B (1993) 405:326–66. [hep-th/9210099]. 10.1016/0550-3213(93)90550-9 Lovelock D . The Einstein tensor and its generalizations. J Math Phys (1971) 12:498–501. 10.1063/1.1665613 Lovelock D . The four-dimensionality of space and the einstein tensor. J Math Phys (1972) 13:874–6. 10.1063/1.1666069 Lodge OJ . Gravitation and light. Nature (1919) 104:354. 10.1038/104354a0 Chwolson O . Traité de Physique, 5 vols. Paris: Hermann (1906-1919). Chwolson O Chwolson O . Astr Nach (1924) 221:329–30. Hewitt JN Turner EL Schneider DP Burke BF Langston GI Lawrence CR . Unusual radio source MG1131+0456: A possible einstein ring. Nature (1988) 333:537–40. 10.1038/333537a0 Hewitt JN Turner EL Lawrence CR Schneider DP Brody JP . A gravitational lens candidate with an unusually red optical counterpart. Astron J (1992) 104:968. 10.1086/116290 Hewitt JN Chen GH Messier MD . Variability in the einstein ring gravitational lens MG 1131+0456. Astron J (1995) 109:1956. 10.1086/117421 Hammer F Le Fevre O Angonin MC Meylan G Smette A Surdej J 1131+0456: Discovery of the optical einstein ring with the NTT, Astron Astrophys 250 (1991) L5. James Webb Space Telescope Goddard Space Flight Center. James Webb space telescope goddard space flight center (2022). Available from: https://webb.nasa.gov/ . Schneider P Ehlers J Falco EE . Gravitational lenses (1992). Perlick V . Gravitational lensing from a spacetime perspective (2010). 1010.3416. Wambsganss J . Gravitational lensing in astronomy. Living Rev Relativity (1998) 1:12. [astro-ph/9812021]. 10.12942/lrr-1998-12 Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Asada K Azulay R Baczko A-K First M87 event horizon telescope results. I. The shadow of the supermassive black hole. Astrophys J Lett (2019) 875:L1. [1906.11238]. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Asada K Azulay R Baczko A-K First M87 event horizon telescope results. II. Array and instrumentation. Astrophys J Lett (2019) 875(1):L2. [1906.11239]. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Asada K Azulay R Baczko A-K First M87 event horizon telescope results. III. Data processing and calibration. Astrophys J Lett (2019) 875(1):L3. [1906.11240]. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Asada K Azulay R Baczko A-K First M87 event horizon telescope results. IV. Imaging the central supermassive black hole. Astrophys J Lett (2019) 875(1):L4. [1906.11241]. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Asada K Azulay R Baczko A-K First M87 event horizon telescope results. V. Physical origin of the asymmetric ring. Astrophys J Lett (2019) 875(1):L5. [1906.11242]. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Asada K Azulay R Baczko A-K First M87 event horizon telescope results. VI. The shadow and mass of the central black hole. Astrophys J Lett (2019) 875(1):L6. [1906.11243]. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Algaba JC Anantua R Asada K First Sagittarius A* event horizon telescope results. I. The shadow of the supermassive black hole in the center of the milky way. Astrophys J Lett (2022) 930(2):L12. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Algaba JC Anantua R Asada K First Sagittarius A* event horizon telescope results. II. EHT and multiwavelength observations, data processing, and calibration. Astrophys J Lett (2022) 930(2):L13. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Algaba JC Anantua R Asada K First Sagittarius A* event horizon telescope results. III. Imaging of the galactic center supermassive black hole. Astrophys J Lett (2022) 930(2):L14. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Algaba JC Anantua R Asada K First Sagittarius A* event horizon telescope results. IV. Variability, morphology, and black hole mass. Astrophys J Lett (2022) 930(2):L15. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Algaba JC Anantua R Asada K First Sagittarius A* event horizon telescope results. V. Testing astrophysical models of the galactic center black hole. Astrophys J Lett (2022) 930(2):L16. Event Horizon Telescope Collaboration, Akiyama K Alberdi A Alef W Algaba JC Anantua R Asada K First Sagittarius A* event horizon telescope results. VI. Testing the black hole metric. Astrophys J Lett (2022) 930(2):L17. Event Horizon Telescope Collaboration, Kocherlakota P Rezzolla L Falcke H Fromm MC Kramer M Mizuno Y Constraints on black-hole charges with the 2017 EHT observations of M87*. Phys Rev D (2021) 103(10):104047. [2105.09343]. LIGO Scientific Virgo Collaboration, Abbott BP Abbott R Abbott TD Abernathy MR Acernese F Ackley K Observation of gravitational waves from a binary black hole merger. Phys Rev Lett (2016) 116(6):061102. [1602.03837]. 10.1103/PhysRevLett.116.061102 LIGO Scientific Virgo Collaboration, Abbott BP Abbott R Abbott TD Abernathy MR Acernese F Ackley K Properties of the binary black hole merger GW150914. Phys Rev Lett (2016) 116(24):241102. [1602.03840]. 10.1103/PhysRevLett.116.241102 LIGO Scientific Virgo Collaboration, Abbott BP Abbott R Abbott TD Abernathy MR Acernese F Ackley K GW151226: Observation of gravitational waves from a 22-solar-mass binary black hole coalescence. Phys Rev Lett (2016) 116(24):241103. [1606.04855]. 10.1103/PhysRevLett.116.241103 LIGO Scientific VIRGO Collaboration, Abbott BP Abbott R Abbott TD Acernese F Ackley K Adams C GW170104: Observation of a 50-solar-mass binary black hole coalescence at redshift 0.2. Phys Rev Lett 118 (2017), no. 22, 221101 [1706.01812], [Erratum: Phys.Rev.Lett. 121, 129901 (2018)]. LIGO Scientific KAGRAVIRGO Collaboration, Abbott R Abbott TD Abraham S Acernese F Ackley K Adams C Observation of gravitational waves from two neutron star–black hole coalescences. Astrophys J Lett (2021) 915(1):L5. [2106.15163]. KAGRA Collaboration, Kokeyama K . Observing the universe from underground gravitational wave telescope KAGRA. In: 3rd world summit on exploring the dark side of the universe (2020). p. 41–8. LIGO Scientific Virgo Collaboration, Abbott BP Abbott R Abbott TD Acernese F Ackley K Adams C Tests of general relativity with the binary black hole signals from the LIGO-virgo catalog GWTC-1. Phys Rev D (2019) 100(10):104036. [1903.04467]. LIGO Scientific VIRGOKAGRA Collaboration, Abbott R Abe H Acernese F Ackley K Adhikari N Adhikari RX Tests of general relativity with GWTC-3 (2021). 2112.06861. Psaltis D Abbott R Abe H Acernese F Ackley K Adhikari N Testing general relativity with the event horizon telescope. Gen Rel Grav (2019) 51(10):137. [1806.09740]. 10.1007/s10714-019-2611-5 Himwich E Johnson MD Lupsasca A Strominger A . Universal polarimetric signatures of the black hole photon ring. Phys Rev D (2020) 101(8):084020. [2001.08750]. 10.1103/physrevd.101.084020 Gralla SE . Can the EHT M87 results be used to test general relativity? Phys Rev D (2021) 103(2):024023. [2010.08557]. 10.1103/physrevd.103.024023 Vagnozzi S Roy R Tsai YD Visinelli L . Horizon-scale tests of gravity theories and fundamental physics from the Event Horizon Telescope image of Sagittarius A* (2022). 2205.07787. Völkel SH Barausse E Franchini N Broderick AE . EHT tests of the strong-field regime of general relativity. Class Quant Grav (2021) 38(21):21LT01. [2011.06812]. 10.1088/1361-6382/ac27ed Bambi C Freese K Vagnozzi S Visinelli L . Testing the rotational nature of the supermassive object M87* from the circularity and size of its first image. Phys Rev D (2019) 100(4):044057. [1904.12983]. 10.1103/physrevd.100.044057 Vagnozzi S Visinelli L . Hunting for extra dimensions in the shadow of M87*. Phys Rev D (2019) 100(2):024020. [1905.12421]. 10.1103/physrevd.100.024020 Allahyari A Khodadi M Vagnozzi S Mota DF . Magnetically charged black holes from non-linear electrodynamics and the Event Horizon Telescope. JCAP (2020) 02:003. [1912.08231]. 10.1088/1475-7516/2020/02/003 Khodadi M Allahyari A Vagnozzi S Mota DF . Black holes with scalar hair in light of the Event Horizon Telescope. JCAP (2020) 09:026. [2005.05992]. 10.1088/1475-7516/2020/09/026 Cunha PVP Herdeiro CAR Radu E Runarsson HF . Shadows of Kerr black holes with scalar hair. Phys Rev Lett (2015) 115(21):211102. [1509.00021]. 10.1103/physrevlett.115.211102 Wang M Chen S Jing J . Shadow casted by a Konoplya-Zhidenko rotating non-Kerr black hole. JCAP (2017) 10:051. [1707.09451]. 10.1088/1475-7516/2017/10/051 Younsi Z Zhidenko A Rezzolla L Konoplya R Mizuno Y . New method for shadow calculations: Application to parametrized axisymmetric black holes. Phys Rev D (2016) 94(8):084025. [1607.05767]. 10.1103/physrevd.94.084025 Bisnovatyi-Kogan GS Tsupko OY . Shadow of a black hole at cosmological distances. Phys Rev D (2018) 98(8):084020. [1805.03311]. 10.1103/physrevd.98.084020 Banerjee I Chakraborty S SenGupta S . Silhouette of M87*: A new window to peek into the world of hidden dimensions. Phys Rev D (2020) 101(4):041301. [1909.09385]. 10.1103/physrevd.101.041301 Tsupko OY Bisnovatyi-Kogan GS . First analytical calculation of black hole shadow in McVittie metric. Int J Mod Phys D (2020) 29(09):2050062. [1912.07495]. 10.1142/s0218271820500625 Tsupko OY Fan Z Bisnovatyi-Kogan GS . Black hole shadow as a standard ruler in cosmology. Class Quant Grav (2020) 37(6):065016. [1905.10509]. 10.1088/1361-6382/ab6f7d Mishra AK Chakraborty S Sarkar S . Understanding photon sphere and black hole shadow in dynamically evolving spacetimes. Phys Rev D (2019) 99(10):104080. [1903.06376]. 10.1103/physrevd.99.104080 Vagnozzi S Bambi C Visinelli L . Concerns regarding the use of black hole shadows as standard rulers. Class Quant Grav (2020) 37(8):087001. [2001.02986]. 10.1088/1361-6382/ab7965 Li P-C Guo M Chen B . Shadow of a spinning black hole in an expanding universe. Phys Rev D (2020) 101(8):084041. [2001.04231]. 10.1103/physrevd.101.084041 Perlick V Tsupko OY . Light propagation in a plasma on Kerr spacetime: Separation of the Hamilton-Jacobi equation and calculation of the shadow. Phys Rev D (2017) 95(10):104003. [1702.08768]. 10.1103/physrevd.95.104003 Wei SW Liu YX Mann RB . Intrinsic curvature and topology of shadows in Kerr spacetime. Phys Rev D (2019) 99(4):041303. [1811.00047]. 10.1103/physrevd.99.041303 Chowdhuri A Bhattacharyya A . Shadow analysis for rotating black holes in the presence of plasma for an expanding universe. Phys Rev D (2021) 104(6):064039. [2012.12914]. 10.1103/physrevd.104.064039 Papnoi U Atamurotov F Ghosh SG Ahmedov B . Shadow of five-dimensional rotating Myers-Perry black hole. Phys Rev D (2014) 90(2):024073. [1407.0834]. 10.1103/physrevd.90.024073 Lin F-L Patel A Pu HY . Black hole shadow with soft hair (2022). 2202.13559. Chandrasekhar S . The mathematical theory of black holes. Oxford classic texts in the physical sciences. Oxford: Oxford Univ. Press (2002). Adler SL Virbhadra KS . Cosmological constant corrections to the photon sphere and black hole shadow radii (2022). 2205.04628. Virbhadra KS . Compactness of supermassive dark objects at galactic centers (2022). 2204.01792. Roy R Vagnozzi S Visinelli L . Superradiance evolution of black hole shadows revisited. Phys Rev D (2022) 105(8):083002. [2112.06932]. 10.1103/physrevd.105.083002 Chen Y Roy R Vagnozzi S Visinelli L . Superradiant evolution of the shadow and photon ring of Sgr A* (2022). 2205.06238. Bozza V . Gravitational lensing in the strong field limit. Phys Rev D (2002) 66:103001. [gr-qc/0208075]. 10.1103/physrevd.66.103001 Virbhadra KS Ellis GFR . Schwarzschild black hole lensing. Phys Rev D (2000) 62:084003. [astro-ph/9904193]. 10.1103/physrevd.62.084003 Frittelli S Kling TP Newman ET . Space-time perspective of Schwarzschild lensing. Phys Rev D (2000) 61:064021. [gr-qc/0001037]. 10.1103/physrevd.61.064021 Bozza V . Quasiequatorial gravitational lensing by spinning black holes in the strong field limit. Phys Rev D (2003) 67:103006. [gr-qc/0210109]. 10.1103/physrevd.67.103006 Bozza V . Comparison of approximate gravitational lens equations and a proposal for an improved new one. Phys Rev D (2008) 78:103005. [0807.3872]. 10.1103/physrevd.78.103005 Einstein A . Lens-like action of a star by the deviation of light in the gravitational field. Science (1936) 84:506–7. 10.1126/science.84.2188.506 Eiroa EF Sendra CM . Gravitational lensing by a regular black hole. Class Quant Grav (2011) 28:085008. [1011.2455]. 10.1088/0264-9381/28/8/085008 Bin-Nun AY . Strong gravitational lensing by Sgr A*. Quant Grav (2011) 28:114003. [1011.5848]. 10.1088/0264-9381/28/11/114003 Amarilla L Eiroa EF Giribet G . Null geodesics and shadow of a rotating black hole in extended Chern-Simons modified gravity. Phys Rev D (2010) 81:124045. [1005.0607]. 10.1103/physrevd.81.124045 Stefanov IZ Yazadjiev SS Gyulchev GG . Connection between black-hole quasinormal modes and lensing in the strong deflection limit. Phys Rev Lett (2010) 104:251103. [1003.1609]. 10.1103/physrevlett.104.251103 Wei S-W Liu YX Fu CE Yang K . Strong field limit analysis of gravitational lensing in Kerr-Taub-NUT spacetime. JCAP (2012) 10:053. [1104.0776]. 10.1088/1475-7516/2012/10/053 Chen S Liu Y Jing J . Strong gravitational lensing in a squashed Kaluza-Klein Gödel black hole. Phys Rev D (2011) 83:124019. [1102.0086]. 10.1103/physrevd.83.124019 Gyulchev GN Stefanov IZ . Gravitational lensing by phantom black holes. Phys Rev D (2013) 87(6):063005. [1211.3458]. 10.1103/physrevd.87.063005 Tsukamoto N Harada T Yajima K . Can we distinguish between black holes and wormholes by their Einstein ring systems? Phys Rev D (2012) 86:104062. [1207.0047]. 10.1103/physrevd.86.104062 Sahu S Patil M Narasimha D Joshi PS . Can strong gravitational lensing distinguish naked singularities from black holes? Phys Rev D (2012) 86:063010. [1206.3077]. 10.1103/physrevd.86.063010 Chen S Jing J . Strong gravitational lensing by a rotating non-Kerr compact object. Phys Rev D (2012) 85:124029. [1204.2468]. 10.1103/physrevd.85.124029 Wei SW Liu YX . Observing the shadow of Einstein-Maxwell-Dilaton-Axion black hole. JCAP (2013) 11:063. [1311.4251]. 10.1088/1475-7516/2013/11/063 Atamurotov F Abdujabbarov A Ahmedov B . Shadow of rotating non-Kerr black hole. Phys Rev D (2013) 88(6):064004. 10.1103/physrevd.88.064004 Eiroa EF Sendra CM . Regular phantom black hole gravitational lensing. Phys Rev D (2013) 88(10):103007. [1308.5959]. 10.1103/physrevd.88.103007 Wei SW Yang K Liu YX . Black hole solution and strong gravitational lensing in Eddington-inspired Born–Infeld gravity. Eur Phys J C (2015) 75:253. [1405.2178], [Erratum: Eur.Phys.J.C 75, 331 (2015)]. 10.1140/epjc/s10052-015-3469-7 Tsukamoto N Kitamura T Nakajima K Asada H . Gravitational lensing in Tangherlini spacetime in the weak gravitational field and the strong gravitational field. Phys Rev D (2014) 90(6):064043. [1402.6823]. 10.1103/physrevd.90.064043 Liu X Jia J Yang N . Gravitational lensing of massive particles in Schwarzschild gravity. Class Quant Grav (2016) 33(17):175014. [1512.04037]. 10.1088/0264-9381/33/17/175014 Wei SW Liu YX Fu CE . Null geodesics and gravitational lensing in a nonsingular spacetime. Adv High Energ Phys. (2015) 2015:1–11. [1510.02560]. 10.1155/2015/454217 Sotani H Miyamoto U . Strong gravitational lensing by an electrically charged black hole in Eddington-inspired Born-Infeld gravity. Phys Rev D (2015) 92(4):044052. [1508.03119]. 10.1103/physrevd.92.044052 Bisnovatyi-Kogan GS Tsupko OY . Gravitational lensing in plasmic medium. Plasma Phys Rep (2015) 41:562–81. [1507.08545]. 10.1134/s1063780x15070016 Sharif M Iftikhar S . Strong gravitational lensing in non-commutative wormholes. Astrophys Space Sci (2015) 357(1):85. 10.1007/s10509-015-2231-9 Younas A Hussain S Jamil M Bahamonde S . Strong gravitational lensing by kiselev black hole. Phys Rev D (2015) 92(8):084042. [1502.01676]. 10.1103/physrevd.92.084042 Chen S Jing J . Strong gravitational lensing for the photons coupled to Weyl tensor in a Schwarzschild black hole spacetime. JCAP (2015) 10:002. [1502.01088]. 10.1088/1475-7516/2015/10/002 Schee J Stuchlík Z Ahmedov B Abdujabbarov A Toshmatov B . Gravitational lensing by regular black holes surrounded by plasma. Int J Mod Phys D (2017) 26(5):1741011. 10.1142/s0218271817410115 Man J Cheng H . Analytical discussion on strong gravitational lensing for a massive source with af(R)global monopole. Phys Rev D (2015) 92(2):024004. [1205.4857]. 10.1103/physrevd.92.024004 Tsukamoto N . Deflection angle in the strong deflection limit in a general asymptotically flat, static, spherically symmetric spacetime. Phys Rev D (2017) 95(6):064035. [1612.08251]. 10.1103/physrevd.95.064035 Wang S Chen S Jing J . Strong gravitational lensing by a Konoplya-Zhidenko rotating non-Kerr compact object. JCAP (2016) 11:020. [1609.00802]. 10.1088/1475-7516/2016/11/020 Tsukamoto N . Strong deflection limit analysis and gravitational lensing of an Ellis wormhole. Phys Rev D (2016) 94(12):124001. [1607.07022]. 10.1103/physrevd.94.124001 Zhao F Tang J He F . Gravitational lensing effects of a Reissner–Nordstrom–de Sitter black hole. Phys Rev D (2016) 93(12):123017. 10.1103/physrevd.93.123017 Aldi GF Bozza V . Relativistic iron lines in accretion disks: The contribution of higher order images in the strong deflection limit. JCAP (2017) 02:033. [1607.05365]. 10.1088/1475-7516/2017/02/033 Lu X Yang FW Xie Y . Strong gravitational field time delay for photons coupled to Weyl tensor in a Schwarzschild black hole. Eur Phys J C (2016) 76(7):357. [1606.02932]. 10.1140/epjc/s10052-016-4218-2 Cavalcanti RT da Silva AG da Rocha R . Strong deflection limit lensing effects in the minimal geometric deformation and Casadio–Fabbri–Mazzacurati solutions. Class Quant Grav (2016) 33(21):215007. [1605.01271]. 10.1088/0264-9381/33/21/215007 Zhao S-S Xie Y . Strong field gravitational lensing by a charged Galileon black hole. JCAP (2016) 07:007. [1603.00637]. 10.1088/1475-7516/2016/07/007 Shaikh R Kar S . Gravitational lensing by scalar-tensor wormholes and the energy conditions. Phys Rev D (2017) 96(4):044037. [1705.11008]. 10.1103/physrevd.96.044037 Rahman M Sen AA . Astrophysical signatures of black holes in generalized proca theories. Phys Rev D (2019) 99(2):024052. [1810.09200]. 10.1103/physrevd.99.024052 Virbhadra KS . Distortions of images of Schwarzschild lensing (2022). 2204.01879. Hsieh T Lee DS Lin CY . Strong gravitational lensing by Kerr and Kerr-Newman black holes. Phys Rev D (2021) 103(10):104063. [2101.09008]. 10.1103/physrevd.103.104063 Hsieh T Lee DS Lin CY . Gravitational time delay effects by Kerr and Kerr-Newman black holes in strong field limits. Phys Rev D (2021) 104(10):104013. [2108.05006]. 10.1103/physrevd.104.104013 Edery A Godin J . Second order Kerr deflection. Gen Rel Grav (2006) 38:1715–22. 10.1007/s10714-006-0347-5 Ghosh S Bhattacharyya A . Analytical study of gravitational lensing in Kerr-Newman black-bounce spacetime. JCAP (2022) 11:006. [2206.09954]. 10.1088/1475-7516/2022/11/006 Balazs NL . Effect of a gravitational field, due to a rotating body, on the plane of polarization of an electromagnetic wave. Phys Rev (1958) 110:236–9. 10.1103/physrev.110.236 Plebanski J . Electromagnetic waves in gravitational fields. Phys Rev (1960) 118:1396–408. 10.1103/physrev.118.1396 Atkinson RD . Light delays near a relativistic star. Astron J (1969) 74:320. 10.1086/110813 de Felice F . On the gravitational field acting as an optical medium. Gen Relativity Gravitation (1971) 2:347–57. 10.1007/bf00758153 Mashhoon B . Scattering of electromagnetic radiation from a black hole. Phys Rev D (1973) 7:2807–14. 10.1103/physrevd.7.2807 Mashhoon B . Influence of gravitation on the propagation of electromagnetic radiation. Phys Rev D (1975) 11:2679–84. 10.1103/physrevd.11.2679 Fischbach E Freeman BS . Second-order contribution to the gravitational deflection of light. Phys Rev D (1980) 22:2950–2. 10.1103/physrevd.22.2950 Sereno M . Gravitational lensing in metric theories of gravity. Phys Rev D (2003) 67:064007. 10.1103/physrevd.67.064007 Ye X-H Lin Q . A Simple optical analysis of gravitational lensing. J Mod Opt (2008) 55:1119–26. [0704.3485]. 10.1080/09500340701618395 Sen AK . A more exact expression for the gravitational deflection of light, derived using material medium approach. Astrophysics (2010) 53:560–9. 10.1007/s10511-010-9148-3 Roy S Sen AK . Trajectory of a light ray in Kerr field: A material medium approach. Astrophys Space Sci (2015) 360(1):23. [1408.3212]. 10.1007/s10509-015-2538-6 Chakraborty S Sen AK . Deflection of light due to rotating mass – A comparison among the results of different approaches. J Phys Conf Ser (2014) 481:012008. 10.1088/1742-6596/481/1/012008 Chakraborty S Sen AK . Light deflection in Kerr field for off-equatorial source (2015). 1504.03124. Simpson A Visser M . Black-bounce to traversable wormhole. JCAP (2019) 02:042. [1812.07114]. 10.1088/1475-7516/2019/02/042 Lobo FSN Rodrigues ME Silva MVDS Simpson A Visser M . Novel black-bounce spacetimes: Wormholes, regularity, energy conditions, and causal structure. Phys Rev D (2021) 103(8):084052. [2009.12057]. 10.1103/physrevd.103.084052 Mazza J Franzin E Liberati S . A novel family of rotating black hole mimickers. JCAP (2021) 04:082. [2102.01105]. 10.1088/1475-7516/2021/04/082 Franzin E Liberati S Mazza J Simpson A Visser M . Charged black-bounce spacetimes. JCAP (2021) 07:036. [2104.11376]. 10.1088/1475-7516/2021/07/036 Tsukamoto N . Gravitational lensing in the Simpson-Visser black-bounce spacetime in a strong deflection limit. Phys Rev D (2021) 103(2):024033. [2011.03932]. 10.1103/physrevd.103.024033 Nascimento JR Petrov AY Porfirio PJ Soares AR . Gravitational lensing in black-bounce spacetimes. Phys Rev D (2020) 102(4):044021. [2005.13096]. 10.1103/physrevd.102.044021 Guerrero M Olmo GJ Rubiera-Garcia D Gómez DSC . Shadows and optical appearance of black bounces illuminated by a thin accretion disk. JCAP (2021) 08:036. [2105.15073]. 10.1088/1475-7516/2021/08/036 Islam SU Kumar J Ghosh SG . Strong gravitational lensing by rotating Simpson-Visser black holes. JCAP (2021) 10:013. [2104.00696]. 10.1088/1475-7516/2021/10/013 Will CM . Gravity: Newtonian, post-Newtonian, and general relativistic (2016). p. 9–72. Iyer SV Hansen EC . Light’s bending angle in the equatorial plane of a Kerr black hole. Phys Rev D (2009) 80:124023. [0907.5352]. 10.1103/physrevd.80.124023 Hsiao YW Lee DS Lin CY . Equatorial light bending around Kerr-Newman black holes. Phys Rev D (2020) 101(6):064070. [1910.04372]. 10.1103/physrevd.101.064070 Claudel CM Virbhadra KS Ellis GFR . The Geometry of photon surfaces. J Math Phys (2001) 42:818–38. [gr-qc/0005050]. 10.1063/1.1308507 Müller T . Einstein rings as a tool for estimating distances and the mass of a Schwarzschild black hole. Phys Rev D (2008) 77:124042. 10.1103/physrevd.77.124042 Virbhadra KS . Relativistic images of Schwarzschild black hole lensing. Phys Rev D (2009) 79:083004. [0810.2109]. 10.1103/physrevd.79.083004 Bozza V Capozziello S Iovane G Scarpetta G . Strong field limit of black hole gravitational lensing. Gen Rel Grav (2001) 33:1535–48. [gr-qc/0102068]. 10.1023/a:1012292927358 Perlick V Tsupko OY . Calculating black hole shadows: Review of analytical studies. Phys Rept (2022) 947:1–39. [2105.07101]. 10.1016/j.physrep.2021.10.004 Barrientos J Cisterna A Mora N Viganò A . AdS-Taub-NUT spacetimes and exact black bounces with scalar hair. Phys Rev D (2022) 106(2):024038. [2202.06706]. 10.1103/physrevd.106.024038 Hadar S Johnson MD Lupsasca A Wong GN . Photon ring autocorrelations. Phys Rev D (2021) 103(10):104038. [2010.03683]. 10.1103/physrevd.103.104038 Appendix A: Details of the computation of the deflection angle for Kerr–Newman spacetime

Here, we give details of computing the integral mentioned in (3.16). First, we rewrite B(u) in the factorized form. B u = − Q 2 1 − w 2 u − u 1 u − u 2 u − u 3 u − u 4 , where u 1 = X 1 − 2 m − X 2 4 m r 0 , u 2 = 1 r 0 , u 3 = X 1 − 2 m + X 2 4 m r 0 , u 4 = 2 m Q 2 − X 1 2 m r 0 .

Then, by choosing the constants X 1 and X 2, we can write down the roots in the following order: u 1 < u 2 < u 3 < u 4. Here, the positive roots are u 2, u 3, u 4, while the negative root is u 1. Then, we substitute Eqs A2–A5 into (A1) and compare it with the coefficients of u 0, u 2, u 3, u 4 in (3.14). This way, we can eventually extract the roots. That gives [144] Q 2 X 2 2 − X 1 − 2 m X 1 + 6 m + 4 X 1 2 = 16 m 2 r 0 X 1 − r 0 1 + ω 1 − ω , X 2 2 − X 1 2 − 2 m 2 = 8 m X 1 − 2 m Q 2 X 1 − 4 m 2 r 0 Q 2 X 1 − 2 m − 4 m 2 r 0 , X 2 2 − X 1 2 − 2 m 2 1 8 m r 0 3 − Q 2 X 1 32 m 3 r 0 4 = 1 λ 2 1 − ω 2 .

Combining the first and second equation of (A6), we obtain the equation for X 1, which is given by Q 2 2 m X 1 3 − Q 2 + 4 m r 0 X 1 2 + 4 m 2 r 0 + 2 m Q 2 + 8 m 3 r 0 2 Q 2 + 2 m r 0 2 1 + ω 1 − ω X 1 = 4 m 2 Q 2 + 4 m 2 r 0 2 1 + ω 1 − ω + 8 m 2 r 0 3 1 + ω Q 2 1 − ω .

Eq. A7 can be exactly solved. The positive real root X 1(m, ω, Q, r 0) turns out to be the one mentioned in (3.19) of Section 3.3.

From (3.20), one can easily check that when Q becomes zero, δ becomes π. In that limit, (3.19) can be written as follows: X 1 m , ω , Q = 0 , r 0 = r 0 1 + ω 1 − ω .

It then exactly reproduces the result for the Kerr black hole [143]. Furthermore, we can reproduce the result for the Reissner–Nordstrom black hole in the limit, a = 0, and the root of the equation (A7) then reduces to X 1 m , Q , ω = 1 , r c = 2 m 2 m r 0 Q 2 − 1 r c , where r c is defined in (3.11), and the result (A9) matches with the result of [144].

Now, we try to factorize the remaining part of the integrand of (3.16). We obtain 1 − 2 m u 1 − ω + Q 2 u 2 1 − ω 1 − 2 m u + a 2 + Q 2 u 2 = G + u + − u + G − u − − u + G Q + u u + − u + G Q − u u − − u , where u ± = m ± m 2 − a 2 + Q 2 a 2 + Q 2 .

G −, G +, G Q+, G Q− are defined in (3.18) of Section 3.3. Also, we can easily see from (A11) u ± are positive.