Next Article in Journal
Observational Constraints on Dynamical Dark Energy with Pivoting Redshift
Next Article in Special Issue
EHT Constraint on the Ultralight Scalar Hair of the M87 Supermassive Black Hole
Previous Article in Journal
Continuous Gravitational Waves from Neutron Stars: Current Status and Prospects
Previous Article in Special Issue
Magnification Cross Sections for the Elliptic Umbilic Caustic Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Effects of Finite Distance on the Gravitational Deflection Angle of Light

Graduate School of Science and Technology, Hirosaki University, Aomori 036-8561, Japan
*
Author to whom correspondence should be addressed.
Universe 2019, 5(11), 218; https://doi.org/10.3390/universe5110218
Submission received: 19 July 2019 / Revised: 18 October 2019 / Accepted: 21 October 2019 / Published: 1 November 2019
(This article belongs to the Special Issue Gravitational Lensing and Optical Geometry: A Centennial Perspective)

Abstract

:
In order to clarify the effects of the finite distance from a lens object to a light source and a receiver, the gravitational deflection of light has been recently reexamined by using the Gauss–Bonnet (GB) theorem in differential geometry (Ishihara et al. 2016). The purpose of the present paper is to give a short review of a series of works initiated by the above paper. First, we provide the definition of the gravitational deflection angle of light for the finite-distance source and receiver in a static, spherically symmetric and asymptotically flat spacetime. We discuss the geometrical invariance of the definition by using the GB theorem. The present definition is used to discuss finite-distance effects on the light deflection in Schwarzschild spacetime for both the cases of weak deflection and strong deflection. Next, we extend the definition to stationary and axisymmetric spacetimes. We compute finite-distance effects on the deflection angle of light for Kerr black holes and rotating Teo wormholes. Our results are consistent with the previous works if we take the infinite-distance limit. We briefly mention also the finite-distance effects on the light deflection by Sagittarius A*.

1. Introduction

In 1919, the experimental confirmation of the theory of general relativity [1] succeeded [2]. It is the measurement of the gravitational deflection angle of light. Since then, the gravitational deflection angle of light has attracted a lot of attention. Many authors have studied the gravitational deflection of light by black holes [3,4,5,6,7,8,9,10,11,12,13,14,15,16]. The gravitational lens by other objects such as wormholes and gravitational monopoles also has attracted a lot of interest [17,18,19,20,21,22,23,24,25,26,27,28,29,30]. Very recently, the Event Horizon Telescope (EHT) team has reported a direct image of the inner edge of the hot matter around the black hole candidate at the center of M87 galaxy [31,32,33,34,35,36]. The direct imaging of black hole shadows must again and steeply raise the importance of the gravitational deflection of light.
Most of those calculations are based on the coordinate angle. The angle respects the rotational symmetry of the spacetime. Gibbons and Werner (2008) made an attempt at defining, in a more geometrical manner, the deflection angle of light [37]. In their paper, the source and receiver are needed to be located at an asymptotic Minkowskian region. The Gauss–Bonnet theorem was applied to a spatial domain by introducing the optical metric, for which a light ray is expressed as a spatial geodesic curve. Ishihara et al. have successfully extended Gibbons and Werner’s idea such that the source and receiver can be at a finite distance from the lens object [38]. They extend the earlier work to the case of the strong deflection limit, in which the winding number of the photon orbits may be larger than unity [39]. In particular, the asymptotic receiver and source are not needed. Arakida [40] made an attempt to apply the Gauss–Bonnet theorem to quadrilaterals that are not extending to infinity and proposed a new definition of the deflection angle of light, though a comparison between two different manifolds that he proposed is an open issue. Proposing an alternative definition of the deflection angle of light, Crisnejo et al. [41] has recently made a comparison between the alternative definitions in References [38,39,40] and showed by explicit calculations that the definition by Arakida in Reference [40] is different from that by Ishihara et al. [38,39]. Their definition has been applied to study gravitational lensing with a plasma medium [41].
The earlier works [38,39] are restricted within the spherical symmetry. Ono et al. have extended the Gauss–Bonnet method with the optical metric to axisymmetric spacetimes [42]. This extension includes mathematical quantities and calculations, with which most physicists are not very familiar. Therefore, the purpose of this paper provides a review of the series of papers on the gravitational deflection of light for finite-distance sources and receivers. In particular, we hope that the detailed calculations in this paper will be helpful for readers to compute the gravitational deflection of light by the new powerful method. For instance, this new technique has been used to study the gravitational lensing in rotating Teo wormholes [43] and also in Damour–Solodukhin wormholes [44]. This formulation has been successfully used to clarify the deflection of light in a rotating global monopole spacetime with a deficit angle [45].
This paper is organized as follows. Section 2 discusses the definition of the gravitational deflection angle of light in static and spherically symmetric spacetimes. Section 3 considers the weak deflection of light in Schwarzschild spacetime. Section 4 discusses the weak deflection of light in the Kottler spacetime and the Weyl conformal gravity model. The strong deflection of light is examined in Section 5. Sagittarius A* (Sgr A*) is also discussed as an example for possible candidates. In Section 6, we discuss the strong deflection of light with finite-distance corrections in Schwarzschild spacetime. Section 7 proposes the definition of the gravitational deflection angle of light in stationary and axisymmetric spacetimes. Sgr A* is also discussed. The weak deflection of light is discussed for Kerr spacetime in Section 8 and for rotating Teo wormholes in Section 9. Section 10 is a summary of this paper. Appendix A provides the detailed calculations for the Kerr spacetime. Throughout this paper, we use the unit of G = c = 1 and the observer may be called the receiver in order to avoid confusion between r O and r 0 by using r R .

2. Definition of the Gravitational Deflection Angle of Light: Static and Spherically Symmetric Spacetimes

Notation

Following Ishihara et al. [38], this section begins by considering a static and spherically symmetric (SSS) spacetime. The metric of this spacetime can be written as:
d s 2 = g μ ν d x μ d x ν = A ( r ) d t 2 + B ( r ) d r 2 + r 2 d Ω 2 ,
where d Ω 2 d θ 2 + sin 2 θ d ϕ 2 and t, θ , and ϕ are associated with the symmetries of the SSS spacetime. For a metric of the form in Equation (1), we always have to restrict to the domain where A ( r ) and B ( r ) are positive such that a static emitter and a static receiver can exist. The spacetime has a spherical symmetry. Therefore, the photon orbital plane is chosen without loss of generality as the equatorial plane ( θ = π / 2 ). We follow the usual definition of the impact parameter of the light ray as:
b L E = r 2 A ( r ) d ϕ d t .
From d s 2 = 0 for the light ray, the orbit Equation is derived as:
d r d ϕ 2 + r 2 B ( r ) = r 4 b 2 A ( r ) B ( r ) .
Light rays are described by the null condition d s 2 = 0 , which is solved for d t 2 as:
d t 2 = γ I J d x I d x J = B ( r ) A ( r ) d r 2 + r 2 A ( r ) d ϕ 2 ,
where I and J denote 1 and 2 and we used Equation (1). We refer to γ I J as the optical metric. The optical metric can be used to describe a two-dimensional Riemannian space. This Riemannian space is denoted as M opt . The light ray is a spatial geodetic curve in M opt .
In the optical metric space M opt , let Ψ denote the angle between the light propagation direction and the radial direction. A straightforward calculation gives:
cos Ψ = b A ( r ) B ( r ) r 2 d r d ϕ .
This is rewritten as:
sin Ψ = b A ( r ) r ,
where we used Equation (3).
We denote Ψ R and Ψ S as the directional angles of light propagation. Ψ R and Ψ S are measured at the receiver position (R) and the source position (S), respectively. We denote ϕ R S ϕ R ϕ S as the coordinate separation angle between the receiver and source. By using angles Ψ R , Ψ S , and ϕ R S , we define the following:
α Ψ R Ψ S + ϕ R S .
This is a basic tool that was invented in Reference [38]. In the following, we shall prove that the definition by Equation (7) is geometrically invariant [38,39].
Here, we briefly mention the Gauss–Bonnet theorem. T is a two-dimensional orientable surface. Differentiable curves T a ( a = 1 , 2 , , N ) are its boundaries. Please see Figure 1 for the orientable surface. We denote the jump angles between the curves as θ a ( a = 1 , 2 , , N ). The Gauss–Bonnet theorem is as follows [46]:
T K d S + a = 1 N T a κ g d + a = 1 N θ a = 2 π ,
where means the line element of the boundary curve, d S denotes the area element of the surface, K means the Gaussian curvature of the surface T, and κ g is the geodesic curvature of T a . The sign of is chosen to be consistent with the surface orientation.
Suppose a quadrilateral R S . Please see Figure 2 for this. This is made of four lines: (1) the spatial curve for the light ray, (2, 3) two outgoing radial lines from R and from S, and (4) a circular arc segment C r that is centered at the lens with the coordinate radius r C ( r C ) and intersects the radial lines at the receiver or the source. We restrict ourselves within the asymptotically flat spacetime. Then, κ g 1 / r C and d r C d ϕ as r C (See, e.g., Reference [37]). By using them, we find C r κ g d ϕ R S . Applying this result to the Gauss–Bonnet theorem for R S , we obtain:
α = Ψ R Ψ S + ϕ R S = R S K d S .
Therefore, α is shown to be invariant for transformations of the spatial coordinates. In addition, α is well defined even when L is a singular point. This is because the point L does not appear in the surface integral nor in the line integral. Furthermore, α vanishes in Euclidean space. This means α is a measure of the deviation from the flat space.
Here, we explain that α defined in Equation (7) is observable in principle. For simplicity, let us imagine the following ideal situation. The positions of a source and a receiver are known. For instance, we assume that the lens object is the Sun, the receiver is located at the Earth, and the source is a pulsar which radiates radio signals with a constant period in an anisotropic manner. In particular, we assume that the source is one of the known pulsars of which the spin period and pulse signal behaviors such as pulse profiles are well understood. By very accurate radio observations such as Very Long Baseline Interferometry (VLBI), the relative positions of the Earth, Sun, and the pulsar can be determined from the ephemeris. (1) From this, we can know ϕ R S in principle. (2) We can directly measure the angle Ψ R at the Earth between the solar direction and the pulsar direction. (3) More importantly, the direction of radiating pulses that reaches the receiver can be also determined in principle because the viewing angle of the pulsar seen by the receiver is known from the pulse profiles. The viewing angle is changing with time because of the Earth’s motion around the Sun. By using the pulsar position and the pulse radiation direction, we can determine Ψ S . Please see Figure 3 for this situation. We explain in more detail how Ψ S at S can be measured by the observer at R. We consider a pulsar of which the spin axis is known from some astronomical observations. A point is that the spin axis of an isolated pulsar is constant with time. The pulse shape and profile depend on the viewing angle with respect to the spin axis of the pulsar. The Earth moves around the Sun, and hence, the observer sees the same pulsar with different viewing angles with time. Accordingly, the observed pulse shape changes. By observing such a change in the pulse shape, we can in principle determine the intrinsic direction of the radio emission, namely the angle between the spin axis and the direction of the emitted light to the observer. In addition, we can know the intrinsic position (including the radial direction from the lens) of such a known pulsar from the ephemeris. By using the intrinsic position (its radial direction) and emission direction at S, Ψ S can be determined in principle, though it is very difficult with current technology. As a result, we can determine in principle Ψ R Ψ S + ϕ R S from astronomical observations. Namely, α in Equation (7) is observable. Note that this procedure does not need to assume a different spacetime, while such a fiducial spacetime was assumed by Arakida (2018) [40], though the receiver in our universe cannot observe the fiducial different spacetime but can assume (or make theoretical calculations of) some quantities on the different spacetime.
One can easily see that, in the far limit of the source and the receiver, Equation (9) agrees with the deflection angle of light as:
α = 2 0 u 0 d u F ( u ) π .
Here, we define u and u 0 as as the inverse of r and the inverse of the closest approach (often denoted as r 0 ), respectively. F ( u ) is defined as:
F ( u ) d u d ϕ 2 .
F ( u ) can be computed by using Equation (3).
The present paper wishes to avoid the far limit for the following reason. Every observed star and galaxy is never located at infinite distance from us. For instance, we observe finite-redshift galaxies in cosmology. We cannot see objects at infinite redshift (exactly at the horizon). Except for a few rare cases in astronomy, the distance to the light source is much larger than the size of the lens. Therefore, we find a strong motivation for studying a situation in which the distance from the source to the receiver is finite. We define u R and u S as the inverse of r R and r S , respectively, where r R and r S are finite. Equation (7) is rewritten in an explicit form as [38,39]:
α = u R u 0 d u F ( u ) + u S u 0 d u F ( u ) + Ψ R Ψ S .
Here, we assume light rays that have only one local minimum of the radius coordinate between r S and r R . This is valid for normal situations in astronomy. However, we should note that multiple local minima are possible, e.g., if the emitter or the receiver (or both) are between the horizon and the light sphere in the Schwarzschild spacetime or if the emitter and receiver are at different sides of the throat of a wormhole spacetime. For such a case of multiple local minima, Equation (12) has to be modified because it assumes only the local minimum at u = u 0 .

3. Weak Deflection of Light in Schwarzschild Spacetime

In this section, we consider the weak deflection of light in Schwarzschild spacetime, for which the line element becomes:
d s 2 = 1 r g r d t 2 + d r 2 1 r g r + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ,
where r g = 2 M in the geometrical unit. Then, F ( u ) is:
F ( u ) = 1 b 2 u 2 + r g u 3 .
By using Equation (6), Ψ R Ψ S in the Schwarzschild spacetime is expanded as:
Ψ R Sch Ψ S Sch [ arcsin ( b u R ) + arcsin ( b u S ) π ] 1 2 b r g u R 2 1 b 2 u R 2 + u S 2 1 b 2 u S 2 + O ( b r g 2 u S 3 , b r g 2 u R 3 ) .
Note that Ψ R Ψ S π in the Schwarzschild spacetime as u S 0 and u R 0 .

4. Other Examples

This section discusses two examples for a non-asymptotically flat spacetime. One is the Kottler solution to the Einstein Equation. The other is an exact solution in the Weyl conformal gravity. The aim of this study is to give us a suggestion or a speculation. We note that the present formulation is limited within the asymptotic flatness, rigorously speaking. As mentioned in the Introduction, Arakida [40] made an attempt to apply the Gauss–Bonnet theorem to quadrilaterals that are not extending to infinity, though a comparison between two different manifolds that he proposed is an open issue. A more careful study that gives a justification for this speculation or perhaps disproves it will be left for the future.
In this section, we do not assume the source at the past null infinity ( r S ) or the receiver at the future null infinity ( r R ) because A ( r ) diverges or does not exist as r . We keep in mind that the source and receiver are located at finite distances from the lens object. Therefore, we use Equation (12). As mentioned already, Equation (6) is more useful in calculating Ψ R and Ψ S than Equation (5) because Equation (6) requires only the local quantities but not any differentiation. By straightforward calculations, we obtain the following results for the above two models.

4.1. Kottler Solution

We consider the Kottler solution [47]. This solution is written as:
d s 2 = 1 2 M r Λ 3 r 2 d t 2 + d r 2 1 2 M r Λ 3 r 2 + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ,
where the cosmological constant is denoted by Λ .
We use Equation (6) such that Ψ R Ψ S can be expanded in terms of 2 M and Λ as:
Ψ R Ψ S = Ψ R Sch Ψ S Sch b Λ 6 u R 1 b 2 u R 2 b Λ 6 u S 1 b 2 u S 2 + b u R ( 1 + 2 b 2 u R 2 ) 8 ( 1 b 2 u R 2 ) 3 / 2 4 M 2 u R 2 + 4 M Λ 3 u R + Λ 2 9 u R 4 + b u S ( 1 + 2 b 2 u S 2 ) 8 ( 1 b 2 u S 2 ) 3 / 2 4 M 2 u S 2 + 4 M Λ 3 u S + Λ 2 9 u S 4 + O ( M 3 , M 2 Λ , M Λ 2 , Λ 3 ) .
Here, Ψ R Sch Ψ S Sch is a pair of the terms that appear also in a case of the Schwarzschild spacetime. The above expansion of Ψ R Ψ S has a divergent term in the limit as u S 0 and u R 0 . The reason for this divergent behavior is that the spacetime is not asymptotically flat and, therefore, the limit of u S 0 and u R 0 is no longer allowed. Hence, the power series in Equation (17) is mathematically valid only within a convergence radius.
For the Kottler spacetime, F ( u ) becomes:
F ( u ) = 1 b 2 u 2 + r g u 3 + Λ 3 .
We obtain:
ϕ R S = π arcsin ( b u R ) arcsin ( b u S ) + r g b 1 1 b 2 u R 2 1 1 2 b 2 u R 2 + 1 1 b 2 u S 2 1 1 2 b 2 u S 2 + Λ b 3 6 u R 1 b 2 u R 2 + u S 1 b 2 u S 2 + r g Λ b 12 2 3 b 2 u R 2 ( 1 b 2 u R 2 ) 3 2 + 2 3 b 2 u S 2 ( 1 b 2 u S 2 ) 3 2 + O ( r g 2 , Λ 2 ) .
By using Equations (17) and (19), α is obtained as:
α = r g b 1 b 2 u R 2 + 1 b 2 u S 2 Λ b 6 1 b 2 u R 2 u R + 1 b 2 u S 2 u S + r g Λ b 12 1 1 b 2 u R 2 + 1 1 b 2 u S 2 + O ( r g 2 , Λ 2 ) .
This Equation has several divergent terms as b u R 0 and b u S 0 . The apparent divergent is problematic only in the case that the source or receiver is located at the horizon. In other words, all the terms in Equation (20) are finite and, thus, harmless for astronomical situations.

4.2. Weyl Conformal Gravity Case

Next, we consider the Weyl conformal gravity model. This theory was originally suggested by Bach [48]. The SSS solution in this model is expressed by introducing three new parameters that are often denoted as β , γ , and k. For this generalized solution in conformal gravity, Birkhoff’s theorem still holds [49]. The SSS solution in the Weyl gravity model is [50]:
d s 2 = A ( r ) d t 2 + 1 A ( r ) d r 2 + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) , A ( r ) = 1 3 m γ 2 m r + γ r k r 2 .
Here, m β ( 2 3 β γ ) / 2 . k r 2 in the metric plays the same role as the cosmological constant in the Kottler spacetime that has been studied above. Therefore, we omit the r 2 term for simplifying our analysis.
By using Equation (6), we expand Ψ R Ψ S in β and γ . The result is:
Ψ R Ψ S Ψ R Sch Ψ S Sch + b γ 2 u R 1 b 2 u R 2 + u S 1 b 2 u S 2 m γ 2 b u R ( 2 b 2 u R 2 ) ( 1 b 2 u R 2 ) 3 / 2 + b u S ( 2 b 2 u S 2 ) ( 1 b 2 u S 2 ) 3 / 2 + O ( m 2 , γ 2 ) .
We should note that this expansion of Ψ R Ψ S is divergent as u S 0 and u R 0 . This divergent behavior is not so problematic because the limit of u S 0 and u R 0 is not allowed in this spacetime. Hence, we note that, rigorously speaking, Equation (22) is mathematically valid only within a convergence radius.
For the present case omitting k, we obtain:
F ( u ) = 1 b 2 u 2 + 2 m u 3 + γ u 2 γ u .
ϕ R S is computed as:
ϕ R S = [ π arcsin ( b u R ) arcsin ( b u S ) ] + m b 2 b 2 u R 2 1 b 2 u R 2 + 2 b 2 u S 2 1 b 2 u S 2 γ 2 b 1 b 2 u R 2 + b 1 b 2 u R 2 + m γ 2 b 3 u R 3 ( 1 b 2 u R 2 ) 3 / 2 + b 3 u S 3 ( 1 b 2 u S 2 ) 3 / 2 + O ( m 2 , γ 2 ) .
Consequently, we obtain α as:
α = 2 m b 1 b 2 u R 2 + 1 b 2 u S 2 m γ b u R 1 b 2 u R 2 + b u S 1 b 2 u S 2 + O ( m 2 , γ 2 ) .
The linear terms in γ cancel out each other, and they do not appear in the final expression for the deflection angle of light. This result may suggest a correction to the results in previous papers [51,52,53] that reported nonzero contributions from γ .

4.3. Far Source and Receiver

Next, we investigate a situation of a distant source and receiver from the lens object: b u S 1 and b u R 1 . Divergent terms in the deflection angle appear in the limit as b u S 0 . Therefore, we carefully investigate the leading part in a series expansion where the infinite limit is not taken. As a result, approximate expressions for the deflection of light are obtained as follows.
(1) Kottler model:
The expression for ϕ R S in this approximation is the same as the seventh and eighth terms of Equation (5) in Reference [54], the third and fifth terms of Equation (15) in Reference [55], and the second term of Equation (14) in Reference [56]. On the other hand, they [54,55,56] did not take account Ψ R Ψ S . In the far approximation, Equation (20) becomes:
α 4 M b 1 6 Λ b 1 u R + 1 u S + 1 3 M Λ b .
This expression suggestions a correction to the earlier works [54,55,56]. For instance, only the term of ϕ R S was considered in Sereno (2009).
(2) Weyl conformal gravity model:
Next, we consider the Weyl conformal gravity model. The deflection angle of light in the far approximation is computed as:
α 4 m b + O ( m 2 , γ 2 ) ,
where m γ parts from Ψ R Ψ S and from ψ R S cancel out each other. Please see also Equations (22) and (24). For instance, Reference [57] gives the exact expression of the deflection angle for the asymptotic receiver and source in the Kottler and Weyl conformal gravity spacetimes.

5. Extension to the Strong Deflection of Light

In the previous sections, we considered the weak deflection of light: A light ray from the source to the receiver is expressed by a spatial curve. The curve is simply connected. In the strong deflection limit, on the other hand, it is possible that the spatial curve has a winding number with intersection points. We thus divide the whole curve into segments. It is easier to investigate each simple segment.

Loops in the Photon Orbit

We begin with one loop case of the light ray curve. This case is shown by Figure 4.
First, we consider the two quadrilaterals (1) and (2) in Figure 5. They can be constructed by introducing an auxiliary point (P) and, next, by adding auxiliary outgoing radial lines (solid line in this figure) from the point P in quadrilaterals (1) and (2). The point P does not need to be the periastron. The directions of the two auxiliary lines in (1) and (2) are opposite to each other. The two auxiliary lines thus cancel out to make no contributions to α . Here, θ 1 and θ 2 denote the inner angle at the point P in the quadrilateral (1) and that in the quadrilateral (2), respectively. We can see that θ 1 + θ 2 = π . This is because the line from the source to the receiver is a geodesic and the point P is located in this line.
For a quadrilateral in Figure 5, the method in Section 2 is still applicable. By the same way of obtaining Equation (9), we obtain:
α ( 1 ) = ( π θ 1 ) Ψ S + ϕ R S ( 1 ) , α ( 2 ) = Ψ R θ 2 + ϕ R S ( 2 ) .
Here, ϕ R S is divided into two parts: One is ϕ R S ( 1 ) for one quadrilateral, and the other is ϕ R S ( 2 ) for the other quadrilateral.
If r S = r R , quadrilaterals (1) and (2) are symmetric for reflection and ϕ R S ( 1 ) = ϕ R S ( 2 ) = ϕ R S / 2 . If not, ϕ R S ( 1 ) is not the same as ϕ R S ( 2 ) . In any case, however, ϕ R S ( 1 ) + ϕ R S ( 2 ) = ϕ R S . Ψ S and ( π Ψ R ) are the inner angles at S and R, respectively. Therefore:
α = α ( 1 ) + α ( 2 ) = Ψ R Ψ S + ϕ R S ,
where we use θ 1 + θ 2 = π and ϕ R S ( 1 ) + ϕ R S ( 2 ) = ϕ R S . This result is the same as Equation (7), though the validity domain is different.
Next, we investigate a case of two loops shown by Figure 6. For this case, we add lines in order to divide the shape into four quadrilaterals as shown by Figure 7. We immediately find:
α ( 1 ) = ( π θ 1 ) Ψ S + ϕ R S ( 1 ) , α ( 2 ) = ( π θ 3 ) θ 2 + ϕ R S ( 2 ) , α ( 3 ) = ( π θ 5 ) θ 4 + ϕ R S ( 3 ) , α ( 4 ) = Ψ R θ 6 + ϕ R S ( 4 ) ,
where ϕ R S ( 1 ) + ϕ R S ( 2 ) + ϕ R S ( 3 ) + ϕ R S ( 4 ) = ϕ R S . Hence, we obtain:
α = α ( 1 ) + α ( 2 ) + α ( 3 ) + α ( 4 ) = Ψ R Ψ S + ϕ R S ,
where we use θ 1 + θ 2 = θ 3 + θ 4 = θ 5 + θ 6 = π . Equation (31) is obtained for the two-loop case in the same form as Equation (7). A loop does make contributions to α only through the terms of ϕ R S ( 2 ) + ϕ R S ( 3 ) .
Finally, we shall complete the proof. We consider the arbitrary winding number, say W. For this case, we prepare 2 W quadrilaterals. We denote the inner angles at finite distance from L as θ 0 , , θ 2 W in order from S to R as shown by Figure 8. Here, θ 0 = Ψ S and θ 2 W = π Ψ R . Neighboring quadrilaterals (N) and (N+1) make the contribution to α only through ϕ R S ( N ) + ϕ R S ( N + 1 ) . We can understand this by noting that θ 2 N 1 + θ 2 N = θ 2 N + 1 + θ 2 N + 2 = π , and the auxiliary lines cancel out. By induction, therefore, we complete the proof; Equation (7) holds for any winding number.
Equation (7) is equivalent to Equation (12). This is shown by using the orbit Equation. This expression is rearranged as:
α = Ψ R Ψ S + ϕ R S = Ψ R Ψ S + u R 0 d u F ( u ) + u S 0 d u F ( u ) + 2 0 u 0 d u F ( u ) .
We define the difference between the asymptotic deflection angle and the deflection angle for the finite distance case as δ α :
δ α α α .
The meaning of this is the finite-distance correction to the deflection angle of light. By substituting Equations (10) and (32) into Equation (33), we get:
δ α = ( Ψ R Ψ S + π ) + u R 0 d u F ( u ) + u S 0 d u F ( u ) .
This expression implies two origins of the finite-distance corrections. One origin is Ψ R and Ψ S . They are angles that are defined in a curved space. The other origin is the two path integrals. They contain information on the curved space. If we consider a receiver and source in the weak gravitational field (as common in astronomy), the finite-distance correction reflects only the weak field region, even if the light ray passes through a strong field region.

6. Strong Deflection of Light in Schwarzschild Spacetime

In this section, we consider the Schwarzschild black hole. By using F ( u ) given by Equation (14), we solve Equation (32) in an analytic manner. The exact expressions involve incomplete elliptic integrals of the first kind. When the distances from the lens to the source and the receiver are much larger than the impact parameter of light ( r S b , r R b ) but the light ray passes near the photon sphere ( r 0 3 M ), Equation (32) becomes approximately:
α = 2 M b 1 b 2 u R 2 + 1 b 2 u S 2 2 + 2 log 12 ( 2 3 ) r 0 r 0 3 M π + O M 2 r R 2 , M 2 r S 2 , 1 3 M r 0 ,
where we used a logarithmic term [8] in the last term of Equation (32). Here, the dominant terms in Ψ R and Ψ S cancel the terms in the integrals. As a consequence, Ψ R and Ψ S do not appear in the approximate expression of Equation (35).
As mentioned above, it follows that the logarithmic term by the strong gravity is free from finite-distance corrections such as 1 ( b u S ) 2 . By chance, δ α in the strong deflection limit (See Equation (32)) is apparently the same as that for the weak deflection case (See, e.g., Equation (29) in Reference [39]). Therefore, the finite-distance correction in the strong deflection limit is again:
δ α O M b r S 2 + M b r R 2 .
This is the same expression as that for the weak field case (e.g., Reference [38]). Namely, the correction is linear in the impact parameter. The finite-distance correction in the weak deflection case (large b) is thus larger than that in the strong deflection one (small b), if the other parameters remain the same.

Sagittarius A*

Next, we briefly mention an astronomical implication of the strong deflection. One of the most feasible candidates for the strong deflection is Sagittarius * (Sgr A*) that is located at our galactic center. In this case, the receiver distance is much larger than the impact parameter of light and a source star may live in the bulge of our Galaxy.
The apparent size of Sgr A* is expected to be nearly the same as that of the central massive object of M87. However, the finite-distance correction to Sgr A* becomes much larger than that to the M87 case because Sgr A* is much closer to us than M87.
For Sgr A*, Equation (36) is evaluated as:
δ α M b r S 2 10 5 arcsec . × M 4 × 10 6 M b 3 M 0.1 p c r S 2 ,
where the central black hole mass is assumed as M 4 × 10 6 M and we take the limit of strong deflection b 3 M . Rather interestingly, this correction as 10 5 arcsec . will be reachable by the Event Horizon Telescope [31,32,33,34,35,36] and near-future astronomy.
See Figure 9 for numerical estimations of the finite-distance correction by the source distance. This figure and Equation (37) suggest that δ α is ∼ ten (or more) micro arcseconds if a source star is sufficiently close to Sgr A* for instance within a tenth of one parsec from Sgr A*. For such a case, the infinite-distance limit does not hold even though the source is still in the weak field. We should take account of finite-distance corrections that are discussed in this paper.
In the strong deflection case, each orbit around the black hole will have a slightly different r 0 , thereby producing a number of “ghost” images (often called relativistic images). In this paper, detailed calculations about it for the finite-distance source and receiver are not done. It is left for the future.

7. Defining the Gravitational Deflection Angle of Light for a Stationary and Axially Symmetric Spacetime

7.1. Optical Metric for the Stationary, Axisymmetric Spacetime

In this section, a stationary and axisymmetric spacetime is considered, for which we shall discuss how to define the gravitational deflection angle of light especially by using the Gauss–Bonnet theorem [42]. The line element in this spacetime is [58,59,60]:
d s 2 = g μ ν d x μ d x ν = A ( y 1 , y 2 ) d t 2 2 H ( y 1 , y 2 ) d t d ϕ + F ( y 1 , y 2 ) γ p q d y p d y q + D ( y 1 , y 2 ) d ϕ 2 .
Here, p and q mean 1 and 2, γ p q is a two-dimensional symmetric tensor, μ , ν are from 0 to 3, and the t and ϕ coordinates respect the Killing vectors. We rewrite this metric into a form such that γ p q can be diagonal. We prefer to use the polar coordinates rather than the cylindrical ones because the Kerr metric and the rotating Teo wormhole one are usually expressed in the polar coordinates. In this paper, we thus use the polar coordinates. In the cylindrical coordinates, the line element is known as the Weyl-Lewis-Papapetrou form [58,59,60]. Equation (38) is rewritten as
d s 2 = A ( r , θ ) d t 2 2 H ( r , θ ) d t d ϕ + B ( r , θ ) d r 2 + C ( r , θ ) d θ 2 + D ( r , θ ) d ϕ 2 ,
where a local reflection symmetry is assumed with respect to the equatorial plane θ = π 2 .
This assumption is expressed as:
g μ ν θ θ = π 2 = 0 .
The functions are A ( r , θ ) > 0 , B ( r , θ ) > 0 , C ( r , θ ) > 0 , D ( r , θ ) > 0 , and H ( r , θ ) > 0 . This assumption by Equation (40) is needed for the existence of a photon orbit on the equatorial plane. Note that we do not assume the global reflection symmetry with respect to the equatorial plane.
The null condition d s 2 = 0 is solved for d t as [61,62]:
d t = γ i j d x i d x j + β i d x i ,
where i and j denote from 1 to 3 and γ i j and β i are defined as:
γ i j d x i d x j B ( r , θ ) A ( r , θ ) d r 2 + C ( r , θ ) A ( r , θ ) d θ 2 + A ( r , θ ) D ( r , θ ) + H 2 ( r , θ ) A 2 ( r , θ ) d ϕ 2 ,
β i d x i H ( r , θ ) A ( r , θ ) d ϕ .
This spatial metric γ i j ( g i j ) is used in order to define the arc length () along the photon orbit as:
d 2 γ i j d x i d x j ,
for which we define γ i j by γ i j γ j k = δ k i . γ i j defines a 3-dimensional Riemannian space ( 3 ) M , where the photon orbit is a spatial curve. In the appendix of Reference [62], they show that is an affine parameter of a light ray.
If the spacetime is static, spherically symmetric, and asymptotically flat, β i is zero and γ i j is nothing but the optical metric. The photon orbit follows a geodesic in a 3-dimensional Riemannian space. In this section and after, we refer to γ i j as the generalized optical metric. Note that the metric γ i j has been called the Fermat metric and the one-form β i is called the Fermat one-form by some authors.
We apply the Gauss–Bonnet theorem to a surface (See Figure 1). The Gauss-Bonnet theorem is expressed as:
R R S S K d S + R S κ g d + S R κ ¯ g d + [ Ψ R + ( π Ψ S ) + π ] = 2 π ,
where we note that the geodesic curvatures of the path from S to S and the path from R to R are both 0 because these paths are geodesic. κ g is the geodesic curvature of the photon orbit, and κ ¯ g is the geodesic curvature of the circular arc segment with an infinite radius.

7.2. Gaussian Curvature

In this subsection, we examine whether the rotational part ( β i ) of the spacetime makes a contribution to the Gaussian curvature. The Gaussian curvature on the equatorial plane is expressed by using the 2-dimensional Riemann tensor ( 2 ) R r ϕ r ϕ :
K = ( 2 ) R r ϕ r ϕ det γ i j ( 2 ) = 1 det γ i j ( 2 ) ϕ det γ i j ( 2 ) γ r r ( 2 ) ( 2 ) Γ ϕ r r r det γ i j ( 2 ) γ r r ( 2 ) ( 2 ) Γ ϕ r ϕ ,
where ( 2 ) R r ϕ r ϕ and ( 2 ) Γ ß j k are defined by using the generalized optical metric γ i j on the equatorial plane. det γ i j ( 2 ) is the determinant of the generalized optical metric in the equatorial plane.
d S in Equation (45) becomes:
d S = det γ ( 2 ) d r d ϕ .
The surface integration of the Gaussian curvature in Equation (45) is rewritten explicitly as:
R R S S K d S = ϕ S ϕ R r O E K det γ ( 2 ) d r d ϕ ,
where r O E means the solution of the orbit Equation.

7.3. Geodesic Curvature

Let us imagine a parameterized curve in a surface. Roughly speaking, the geodesic curvature of the parameterized curve is a measure of how different the curve is from the geodesic. The geodesic curvature of the parameterized curve is defined as the surface-tangential component of the acceleration (namely the geodesic curvature) of the curve. The normal curvature is defined as the surface-normal component of the acceleration. The normal curvature does not appear in the present paper because we consider only the curves on the equatorial plane.
The geodesic curvature in the vector form is defined as (see, e.g., Reference [63,64]):
κ g T · T × N ,
where, for a parameterized curve, T denotes the unit tangent vector for the curve by reparameterizing the curve using its arc length, T means its derivative with respect to the parameter, and N indicates the unit normal vector for the surface. The geodesic curvature of a curve vanishes if the curve follows the geodesic. This zero is because the acceleration vector T vanishes.

7.4. Photon Orbit with the Generalized Optical Metric

In this subsection, we discuss geometrical aspects of a photon orbit in terms of the generalized optical metric. The unit vector tangent to the spatial curve is generally expressed as:
e i d x i d ,
where a parameter is defined by Equation (44).
The flight time T of a light from the source to the receiver is obtained by performing the integral of Equation (41):
T = t S t R d t = S R γ i j d e i d e j + β i d e i d .
The light ray follows the Fermat’s principle, namely δ T = 0 [65]. The Lagrangian for a photon can be expressed as:
L = γ i j e i e j + β i e i .
From this, we obtain:
d d L e k = γ i k e i , l e l + γ i k , l e i e l + β k , i e i ,
L x k = 1 2 γ i j , k e i e j + β i , k e i ,
where we used γ i j e i e j = 1 and the comma (,) defines the partial derivative. The Euler–Lagrange Equation is calculated as:
e j , l e l + γ k j γ i k , l e i e l 1 2 γ i l , k e i e l = γ k j ( β l , k β k , l ) e l .
This leads to the Equation for the light ray [62]:
d e i d = γ i l ( γ l j , k 1 2 γ j k , l ) e j e k + γ i j ( β k , j β j , k ) e k .
Therefore, the geodesic Equation is equivalent to:
e i | j e j = d e i d + ( 3 ) Γ i j k e j e k = d e i d + γ i l ( γ l j , k 1 2 γ j k , l ) e j e k = γ i j ( β k , j β j , k ) e k ,
where we define | as the covariant derivative with respect to γ i j . ( 3 ) Γ j k i means the Christoffel symbol by γ i j .
The acceleration vector a i is defined by:
a i e i | j e j = γ i j ( β k | j β j | k ) e k = γ i j ( β k , j β j , k ) e k .
By using the Levi–Civita symbol ε i j k , we express the cross (outer) product: Of A and B in the covariant manner:
γ ε i j k A j B k = ( A × B ) i .
The Levi–Civita tensor ϵ i j k is defined by ϵ i j k γ ε i j k , where and ε i j k is the Levi–Civita symbol ( ε 123 = 1 ).
The Levi–Civita tensor ϵ i j k in a three-dimensional satisfies:
ϵ s j k ϵ s l m = γ ε s j k 1 γ ε s l m = δ j l δ k m δ j m δ k l ,
ϵ s j k ϵ l m s = γ j l γ k m γ j m γ k l .
By using Equations (58)–(60), Equation (57) is rewritten as:
a i = γ i j e k ϵ s j k ( × β ) s .
Vector a i is the spatial vector representing the acceleration due to β i . In particular, a i is caused in gravitomagnetism [66]. To be more precise, the gravitomagnetic vector has an analogy to the Lorentz force in electromagnetism v × ( × A m ), in which A m denotes the vector potential. The vector potential is defined as B = × A m , E = ϕ A m t , where E and B are the electric and magnetic fields, respectively, and the electric potential is ϕ .
γ i j is not an induced metric but the generalized optical metric. If β i is nonvanishing, the photon orbit may be different from a geodesic in ( 3 ) M with γ i j , even though the light ray in the four-dimensional spacetime follows the null geodesic.
In a stationary and axisymmetric spacetime, it is always possible to find out coordinates such that g 0 i can vanish and a i = 0 . In this case, the photon orbit is considered a spatial geodesic curve in ( 3 ) M .
We study axisymmetric cases which allow g 0 i 0 . Therefore, geodesic curvature κ g does not always vanish in the photon orbit in the Gauss–Bonnet theorem because the geodesic curvature κ g for a photon orbit is due to the gravitomagnetic effect. This nonvanishing κ g for the photon orbit leads to a crucial difference from the SSS case [38,39].

7.5. Geodesic Curvature of a Photon Orbit

Equation (49) is rearranged to be in the tensor form:
κ g = ϵ i j k N i a j e k ,
where T and T correspond to e k and a j , respectively.
In this paper, the acceleration vector of the photon orbit depends on β i . Hence, the geodesic curvature for the photon orbit also depends on it. A nonvanishing integral of the geodesic curvature along the light ray appears in the Gauss–Bonnet theorem in Equation (8).
Substituting Equation (57) into a i in Equation (62) leads to:
κ g = ϵ i j k N i γ j l ( β n | l β l | n ) e n e k = γ j a N i e k e b ϵ i j k ϵ s a b ϵ s m l β l | m = N i e k e b ( δ i s δ k b δ i b δ k s ) ϵ s m l β l | m = ϵ i j k N i β j | k ,
where we used γ i j e i e j = 1 and γ i j e i N j = 0 . The unit vector normal to the equatorial plane is:
N p = 1 γ θ θ δ p θ ,
where the upward direction is chosen without loss of generality.
For the equatorial plane, we obtain:
ϵ θ p q β q | p = 1 γ β ϕ , r ,
where we use ϵ θ r ϕ = 1 / γ and β r , ϕ = 0 because of the axisymmetry.
By using Equations (64) and (65), κ g in Equation (63) becomes:
κ g = 1 γ γ θ θ β ϕ , r .
By using Equation (44), the line element in the path integral is obtained as:
d = γ r r d r d ϕ 2 + γ ϕ ϕ d ϕ ,
where θ = π / 2 .

7.6. Geodesic Curvature of a Circular arc Segment

In a flat space, the geodesic curvature κ of the circular arc segment of radius R is obtained as:
κ = 1 R .
The geodesic curvature κ ¯ g of a circular arc segment of radius R c = R is obtained as:
κ ¯ g = 1 R c ,
where the radius R c is sufficiently larger than r R and r S and the circular arc segment is in the asymptotically flat region.
Equation (44) becomes d 2 = d r 2 + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) because we assume an asymptotically flat spacetime. Hence, the line element in the path integral of κ ¯ g is obtained as:
d = R c d ϕ ,
where we choose θ = π / 2 and r = R c for the circular arc segment.
Therefore, the path integral of κ ¯ g in Equation (45) is rewritten as:
S R κ ¯ g d = ϕ S ϕ R d ϕ = ϕ R ϕ S = ϕ R S ,
where we denote the angular coordinate values of the receiver and the source as ϕ R and ϕ S , respectively.

7.7. Impact Parameter and Light Rays

By using Equation (39), we study the orbit Equation on the equatorial plane. The Lagrangian for a photon in the equatorial plane is obtained as:
L ^ = A ( r ) t ˙ 2 2 H ( r ) t ˙ ϕ ˙ + B ( r ) r ˙ 2 + D ( r ) ϕ ˙ 2 ,
where the dot denotes the derivative with respect to the affine parameter and the functions A ( r ) , B ( r ) , D ( r ) , and H ( r ) mean, to be rigorous, A ( r , π / 2 ) , B ( r , π / 2 ) , D ( r , π / 2 ) , and H ( r , π / 2 ) respectively.
The metric (or the Lagrangian L ^ in the 4-dimensional spacetime) is independent from t and ϕ . Therefore:
d d L ^ t ˙ = 0 , d d L ^ ϕ ˙ = 0 .
Then, associated with the two Killing vectors ξ μ = ( 1 , 0 , 0 , 0 ) and ξ ¯ μ = ( 0 , 0 , 0 , 1 ) , respectively:
L ^ t ˙ = g μ ν ξ μ k ν , L ^ ϕ ˙ = g μ ν ξ ¯ μ k ν ,
where k μ = d x μ d is the vector tangent to the light ray in the four-dimensional spacetime. There are two constants of motion:
E = A ( r ) t ˙ + H ( r ) ϕ ˙ ,
L = D ( r ) ϕ ˙ H ( r ) t ˙ ,
where E denotes the energy of the photon and L means the angular momentum of the photon. The impact parameter of the photon is defined as:
b L E = H ( r ) t ˙ + D ( r ) ϕ ˙ A ( r ) t ˙ + H ( r ) ϕ ˙ = H ( r ) + D ( r ) d ϕ d t A ( r ) + H ( r ) d ϕ d t .
In terms of the impact parameter b, L ^ = 0 can be considered as the orbit Equation:
d r d ϕ 2 = A ( r ) D ( r ) + H 2 ( r ) B ( r ) D ( r ) 2 H ( r ) b A ( r ) b 2 H ( r ) + A ( r ) b 2 ,
where we used Equation (39). By introducing u 1 / r , we rewrite the orbit Equation as:
d u d ϕ 2 = F ( u ) ,
where F ( u ) is:
F ( u ) = u 4 ( A D + H 2 ) ( D 2 H b A b 2 ) B ( H + A b ) 2 .
We examine the angles ( Ψ R and Ψ S in Figure 10) at the receiver position and the source one. The unit vector tangent to the photon orbit in ( 3 ) M is e i . Its components on the equatorial plane are expressed as:
e i = 1 ξ d r d ϕ , 0 , 1 ,
where ξ satisfies:
1 ξ = A ( r ) [ H ( r ) + A ( r ) b ] A ( r ) D ( r ) + H 2 ( r ) .
This can be derived also from γ i j e i e j = 1 by using Equation (77).
In the equatorial plane, the unit radial vector is:
R i = 1 γ r r , 0 , 0 ,
where the outgoing direction is chosen for a sign convention.
By using the inner product between e i and R i , we therefore define the angle as:
cos Ψ γ i j e i R j γ r r A ( r ) [ H ( r ) + A ( r ) b ] A ( r ) D ( r ) + H 2 ( r ) d r d ϕ ,
where Equations (80)–(82) are used. This is rewritten as:
sin Ψ = H ( r ) + A ( r ) b A ( r ) D ( r ) + H 2 ( r ) ,
where Equation (77) is used. We should note that sin Ψ in Equation (84) is more useful in practical calculations, because it needs only the local quantities. On the other hand, cos Ψ by Equation (83) needs the derivative d r / d ϕ . In addition, the domain of this Ψ is 0 Ψ π and, hence, sin Ψ is always positive.
By substituting r R and r S into r of Equation (84), we obtain sin Ψ R and sin Ψ S , respectively. We note that the range of the principal value of y = arcsin x is π 2 y π 2 as usual. However, the range of Ψ R ( Ψ S ) is 0 Ψ R ( Ψ S ) π . By using the usual principal value, Equation (84) for ( Ψ R ) and ( Ψ S ) becomes:
sin Ψ R = H ( r R ) + A ( r R ) b A ( r R ) D ( r R ) + H 2 ( r R ) ,
sin ( π Ψ S ) = H ( r S ) + A ( r S ) b A ( r S ) D ( r S ) + H 2 ( r S ) ,
respectively because Ψ R is an acute angle and Ψ S is an obtuse angle as shown by Figure 10.

7.8. Gravitational Deflection Light in the Axisymmetric Case

We define:
α Ψ R Ψ S + ϕ R S
for the equatorial plane in the axisymmetric spacetime. This definition apparently depends on the angular coordinate ϕ . By using the Gauss–Bonnet theorem in Equation (8), this Equation is rearranged as:
α = R S K d S R S κ g d .
Here, d is positive when the photon is in the prograde motion, whereas it is negative for the retrograde case. Equation (88) means that α is coordinate-invariant for the axisymmetric case. Up until now, we did not use any Equations for gravitational fields. Therefore, the above discussion and results still stand not only in the theory of general relativity but also in a general class of metric theories of gravity only if the light ray in the four-dimensional spacetime is a null geodesic.

8. Weak Deflection of Light in Kerr Spacetime

8.1. Kerr Spacetime and γ i j

In this section, we focus on the weak deflection of light in the Kerr spacetime as an axisymmetric example. Kerr metric in the Boyer–Lindquist form is expressed as:
d s 2 = 1 2 M r Σ d t 2 4 a M r sin 2 θ Σ d t d ϕ + Σ Δ d r 2 + Σ d θ 2 + r 2 + a 2 + 2 a 2 M r sin 2 θ Σ sin 2 θ d ϕ 2 ,
where Σ and Δ are defined as:
Σ r 2 + a 2 cos 2 θ ,
Δ r 2 2 M r + a 2 .
Using the Gauss–Bonnet theorem, the deflection angle of light in the Kerr spacetime was calculated for the asymptotic source and receiver by Werner [67]. However, his method based on the osculating metric is limited within the asymptotic case. Later, Ono et al. developed a different approach using the Gauss–Bonnet theorem that enables the calculation of the deflection angle for the finite distance case in the Kerr spacetime [42].
By using Equations (42) and (43), the generalized optical metric γ i j and the gravitomagnetic term β i for the Kerr metric are obtained as:
γ i j d x i d x j = Σ 2 Δ ( Σ 2 M r ) d r 2 + Σ 2 ( Σ 2 M r ) d θ 2    + r 2 + a 2 + 2 a 2 M r sin 2 θ ( Σ 2 M r ) Σ sin 2 θ ( Σ 2 M r ) d ϕ 2 ,
β i d x i = 2 a M r sin 2 θ ( Σ 2 M r ) d ϕ .
Note that γ i j has no linear terms in the Kerr spin parameter a because only g 0 i in g μ ν has a linear term in a and g 0 i H contributes to γ i j through a quadratic term g 0 i g 0 j H 2 , as shown by Equation (42).
In order to calculate the Gaussian curvature K of the equatorial plane, the geodesic curvature κ g of the light ray and the geodesic curvature κ g ¯ of the circular arc of an infinite radius and of the angles Ψ R and Ψ S , we use two approximations for the weak field and slow rotation, where M and a play roles as book-keeping parameters though they are dimensional quantities.
By using Equation (77), we obtain the orbit Equation:
d r d ϕ 2 = b 2 a 2 b 2 + r b ( r b 2 M b ) 2 a 2 b 2 ( 2 M b + r b ) 4 a M b 2 + 2 M b r b + r 3 b 3 r b { 2 a M b 2 + r b 2 M b } 2 = r 4 b 2 r 2 + 2 M r 4 r 3 b 3 a M + O ( a 2 ) ,
where the weak-field and slow-rotation approximations are used in the last line. There are no M-squared terms in the last line. The orbit Equation becomes:
d u d ϕ 2 = F ( u ) = 1 b 2 u 2 + 2 M u 3 4 u b 3 a M + O ( a 2 u 4 ) .
We solve iteratively Equation (95). In order to find the zeroth-order solution, we solve the truncated Equation (95):
d u d ϕ 2 = 1 b 2 u 2 + O ( M u 3 , a M u 4 , a 2 u 4 ) .
The zeroth-order solution for this Equation is:
u = sin ϕ b ,
where we use d u d ϕ ϕ = π / 2 = 0 as the boundary condition. This condition means that the closest approach of the photon orbit is expressed as r = r 0 = 1 / u 0 , ϕ = π / 2 . We assume that the linear-order solution with M is u = sin ϕ b + u 1 ( ϕ ) M . In order to obtain u 1 ( ϕ ) , we substitute this expression of u into the Equation (95) with terms linear in M:
d u d ϕ 2 = 1 b 2 u 2 + 2 M u 3 + O ( a M u 4 , a 2 u 4 ) .
u 1 ( ϕ ) is thus obtained as:
u 1 ( ϕ ) = 1 b 2 ( 1 + cos 2 ϕ ) ,
where we used the boundary condition mentioned above. The solution with a is in a form of u = sin ϕ b + M b 2 ( 1 + cos 2 ϕ ) + u 2 ( ϕ ) a . Since Equation (95) does not include any linear term in a, we find u 2 ( ϕ ) = 0 . The solution with a M is u = sin ϕ b + M b 2 ( 1 + cos 2 ϕ ) + u 3 ( ϕ ) a M . We substitute this solution into Equation (95):
a M b b 3 d u 3 ( ϕ ) d ϕ cos ϕ + b 3 u 3 ( ϕ ) sin ϕ + 2 sin ϕ + O ( a 2 u 4 ) = 0 .
Hence, u 3 ( ϕ ) is obtained as:
u 3 ( ϕ ) = 2 b 3 .
Bringing the above results together, the iterative solution of Equation (95) is expressed as:
u = sin ϕ b + M b 2 ( 1 + cos 2 ϕ ) 2 a M b 3 + O M 2 b 3 , a 2 b 3 .
Next, we solve Equation (102) for ϕ . We obtain ϕ as:
ϕ = arcsin ( b u ) + 2 + b 2 u 2 b 1 b 2 u 2 M + 2 a M b 2 1 b 2 u 2 + O M 2 b 3 , a 2 b 3 ( | ϕ | < π 2 ) π arcsin ( b u ) 2 + b 2 u 2 b 1 b 2 u 2 M 2 a M b 2 1 b 2 u 2 + O M 2 b 3 , a 2 b 3 ( π 2 < | ϕ | ) ,
where we can choose the domain of ϕ to be π ϕ < π without loss of generality. In the following, the range of the angular coordinate value ϕ S at the source point is π 2 ϕ S < π 2 and the range of the angular coordinate value ϕ R at the receiver point is | ϕ R | > π 2 . We find | b u | < 1 because the square root in Equation (103) must be real and nonzero, and the values of b and u are positive. Therefore, b u satisfies 0 < b u < 1 in our calculation.

8.2. Gaussian Curvature on the Equatorial Plane

Let us explain how to compute the Gaussian curvature by using Equation (46). In the Kerr case, it becomes:
K = M 6 r a 2 + M 2 + 6 a 2 M + 7 M r 2 2 r 3 r 5 ( r 2 M ) = 2 M r 3 + O M 2 r 4 , a 2 M r 5 ,
where the weak-field and slow-rotation approximations are used in the last line.
Next, we discuss the area element on the equatorial plane by using Equation (47). In the Kerr case, the area element of the equatorial plane is expressed as:
d S = [ r + 3 M + O ( M 2 / r ) ] d r d ϕ .
By using Equations (104) and (105), the surface integral of the Gaussian curvature in Equation (88) is performed as:
R R S S K d S = ϕ S ϕ R r O E ( 2 M r 3 r ) d r d ϕ + O M 2 b 2 , a M 2 b 3 , a 2 M b 3 = 2 M ϕ S ϕ R 0 1 b sin ϕ + M b 2 ( 1 + cos 2 ϕ ) 2 a M b 3 d u d ϕ + O M 2 b 2 , a M 2 b 3 , a 2 M b 3 = 2 M ϕ S ϕ R 1 b sin ϕ d ϕ + O M 2 b 2 , a M 2 b 3 , a 2 M b 3 = 2 M b cos ϕ S cos ϕ R + O M 2 b 2 , a M 2 b 3 , a 2 M b 3 = 2 M b 1 b 2 u S 2 + 1 b 2 u R 2 + O M 2 b 2 , a M 2 b 3 , a 2 M b 3 ,
where r O E in the first line is the solution of Equation (94), we transform the integral variable as r = 1 / u in the second line, and we used cos ϕ S = 1 b 2 u S 2 + O ( M / b ) and cos ϕ R = 1 b 2 u R 2 + O ( M / b ) from Equation (103) in the last line.

8.3. Path Integral of κ g

Substituting Equation (93) into β i in Equation (66) leads to:
κ g = 2 a M r 2 ( r 2 M ) 1 2 M r + a 2 r 2 1 + a 2 r 2 + 2 a 2 M r 3 1 / 2 = 2 a M r 3 + O a M 2 r 4 ,
where the weak-field and slow-rotation approximations are used in the last line. We stress that the terms of a n M ( n 2 ) do not exist in this expression.
The line element for the path integral by Equation (67) becomes:
d = b sin 2 ϕ + O ( M ) d ϕ ,
where Equation (102) was used for a relation between r and ϕ .
By using Equations (107) and (108), the path integral of κ g in Equation (88) is performed as:
R S κ g d = S R 2 a M r 3 d + O a M 2 r 4 = 2 a M b 2 ϕ S ϕ R sin ϕ d ϕ + O a M 2 r 4 = 2 a M b 2 [ 1 b 2 u R 2 + 1 b 2 u S 2 ] + O a M 2 b 3 .
Here, we assumed d > 0 , such that the orbital angular momentum can be parallel with the spin of the black hole, and we used a linear approximation of the photon orbit as 1 / r = u = sin ϕ / b + O ( M / b 2 , a M / b 3 ) from Equation (102). In the retrograde case, d becomes negative and the magnitude of the above path integral thus remains the same but the sign of the integral is opposite.

8.4. ϕ R S Part

The displacement of the angular coordinate ϕ in Equation (87) is computed as:
ϕ R S = S R d ϕ = 2 0 u 0 1 F ( u ) d u + u S 0 1 F ( u ) d u + u R 0 1 F ( u ) d u ,
where the orbit equation of Equation (78) was used. We substitute Equation (95) into F ( u ) in Equation (110) to obtain:
ϕ R S = u S u 0 1 u 0 2 u 2 + M u 0 3 u 3 ( u 0 2 u 2 ) 3 / 2 2 a M u 0 3 ( u 0 u ) ( u 0 2 u 2 ) 3 / 2 d u + u R u 0 1 u 0 2 u 2 + M u 0 3 u 3 ( u 0 2 u 2 ) 3 / 2 2 a M u 0 3 ( u 0 u ) ( u 0 2 u 2 ) 3 / 2 d u + O ( M 2 u 0 2 , a 2 u 0 2 ) = π 2 arcsin u S u 0 + M ( 2 u 0 + u S ) u 0 2 u S 2 u 0 + u S 2 a M u 0 3 u 0 2 u S 2 u 0 2 + u 0 u S + π 2 arcsin u R u 0 + M ( 2 u 0 + u R ) u 0 2 u R 2 u 0 + u R 2 a M u 0 3 u 0 2 u R 2 u 0 2 + u 0 u R + O M 2 u 0 2 , a 2 u 0 2 ,
where the prograde case is assumed. In the retrograde motion, the sign of the linear term in a is opposite. In Equation (111), the impact parameter b is rewritten in terms of the closest approach u 0 for the integration from u S (or u R ) to u 0 . Namely, Equation (95) tells us the relation between the impact parameter b and the inverse of the closest approach u 0 as b = u 0 1 + M 2 a M u 0 + O ( M 2 u 0 , a 2 u 0 ) in the weak-field and slow-rotation approximations. By making use of this relation, Equation (111) is rearranged as:
ϕ R S = π arcsin ( b u S ) arcsin ( b u R ) + M ( 2 b 2 u S 2 ) b 1 b 2 u S 2 + M ( 2 b 2 u R 2 ) b 1 b 2 u R 2 2 a M b 2 1 1 b 2 u S 2 + 1 1 b 2 u R 2 + O M 2 / b 2 , a 2 / b 2 .
The first line of this equation recovers Equation (32) of Reference [38].

8.5. Ψ Parts

In the Kerr spacetime by Equation (89), Equation (85) is:
sin Ψ R = b r R × 1 2 M r R + 2 a M b r R 1 2 M r R + a 2 r R 2 , = b r R 1 M r R + 2 a M b r R + O M 2 r R 2 , a 2 r R 2 , a M 2 r R 3 = b u R 1 M u R + 2 a M u R b + O M 2 u R 2 , a 2 u R 2 , a M 2 u R 3 ,
and Equation (86) is calculated as:
sin ( π Ψ S ) = b u S 1 M u S + 2 a M u S b + O M 2 u S 2 , a 2 u S 2 , a M 2 u S 3 ,
where r R = 1 / u R , r S = 1 / u S and we used the weak-field and slow-rotation approximations. By combining Equations (113) and (114), we obtain Ψ R and Ψ S as:
Ψ R = arcsin b u R 1 M u R + 2 a M u R b + O M 2 u R 2 , a 2 u R 2 , a M 2 u R 3 = arcsin ( b u R ) M b u R 2 1 b 2 u R 2 + 2 a M u R 2 1 b 2 u R 2 + O M 2 u R 2 , a 2 u R 2 , a M 2 u R 3 , π Ψ S = arcsin ( b u S ) M b u S 2 1 b 2 u S 2 + 2 a M u S 2 1 b 2 u S 2 + O M 2 u S 2 , a 2 u S 2 , a M 2 u S 3 .
By combining these relations, we obtain the Ψ part in Equation (87) as:
Ψ R Ψ S = arcsin ( b u R ) + arcsin ( b u S ) π M b u R 2 1 b 2 u R 2 M b u S 2 1 b 2 u S 2 + 2 a M u R 2 1 b 2 u R 2 + 2 a M u S 2 1 b 2 u S 2 + O M 2 u R 2 , M 2 u S 2 , a 2 u R 2 , a 2 u S 2 , a M 2 u R 3 , a M 2 u S 3 .

8.6. Deflection of Light in Kerr Spacetime

On the equatorial plane in the Kerr spacetime, the deflection angle of light is described by Equations (87) and (88). Let us examine whether the two results agree with each other.
First, we substitute Equations (112) and (116) into Equation (87). We obtain the deflection angle of light as:
α p r o g = arcsin ( b u R ) + arcsin ( b u S ) π M b u R 2 1 b 2 u R 2 M b u S 2 1 b 2 u S 2 + 2 a M u R 2 1 b 2 u R 2 + 2 a M u S 2 1 b 2 u S 2 + π arcsin ( b u S ) arcsin ( b u R ) + M ( 2 b 2 u S 2 ) b 1 b 2 u S 2 + M ( 2 b 2 u R 2 ) b 1 b 2 u R 2 2 a M b 2 1 1 b 2 u S 2 + 1 1 b 2 u R 2 + O M 2 b 2 = 2 M b 1 b 2 u R 2 + 1 b 2 u S 2 2 a M b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O M 2 b 2 ,
where the prograde orbit of light is assumed. For the retrograde motion, we obtain:
α r e t r o = 2 M b 1 b 2 u R 2 + 1 b 2 u S 2 + 2 a M b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O M 2 b 2 .
Next, we substitute Equations (106) and (109) into Equation (88). Then, we obtain the deflection angle of light in the prograde motion as:
α p r o g = 2 M b 1 b 2 u R 2 + 1 b 2 u S 2 2 a M b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O M 2 b 2 ,
and the deflection angle for the retrograde case as:
α r e t r o = 2 M b 1 b 2 u R 2 + 1 b 2 u S 2 + 2 a M b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O M 2 b 2 .
Note that the a 2 terms in the deflection angle in Equation (87) cancel out thanks to Equation (88).
Here, we consider the limit as u R 0 and u S 0 . In this limit, we get:
α p r o g 4 M b 4 a M b 2 + O M 2 b 2 ,
α r e t r o 4 M b + 4 a M b 2 + O M 2 b 2 .
This shows that Equations (117) and (118) agree with the asymptotic deflection angles that are known in earlier works [4,68,69,70]. Precise analytic treatments of the deflection angle of light were done in a conventional approach, on the equatorial plane of a Kerr black hole [70] and for generic photon orbits in terms of the generalized hypergeometric functions of Appell and Lauricella [71]. They assume that both the source and the receiver are located at the null infinity.
If we wish to consider the deflection angle of light in a case where the receiver point is closer to the source point than the closest approach point, Equations (117) and (118) become:
α p r o g = 2 M b 1 b 2 u S 2 1 b 2 u R 2 2 a M b 2 1 b 2 u S 2 1 b 2 u R 2 + O M 2 b 2 , α r e t r o = 2 M b 1 b 2 u S 2 1 b 2 u R 2 + 2 a M b 2 1 b 2 u S 2 1 b 2 u R 2 + O M 2 b 2 .
If we wish to consider the deflection angle of light in such a case that the source point is closer to the receiver than the closest approach point, Equations (117) and (118) become:
α p r o g = 2 M b 1 b 2 u R 2 1 b 2 u S 2 2 a M b 2 1 b 2 u R 2 1 b 2 u S 2 + O M 2 b 2 , α r e t r o = 2 M b 1 b 2 u R 2 1 b 2 u S 2 + 2 a M b 2 1 b 2 u R 2 1 b 2 u S 2 + O M 2 b 2 .

8.7. Finite-Distance Corrections

In the previous subsections so far, we discussed an effect of the spin of the lens object to the deflection of light. In particular, we do not require the receiver and the source to be located at infinity. The finite-distance correction to the deflection angle of light is defined as δ α . This is the difference between the asymptotic deflection angle α and the deflection angle for the finite distance case. Namely:
δ α α α .
Equations (117) and (118) tell us the magnitude of the finite-distance correction to the gravitomagnetic bending angle due to the spin. The result is:
δ α G M O a M r S 2 + a M r R 2 O J r S 2 + J r R 2 ,
where b u R , b u S < 1 is assumed, J a M denotes the spin angular momentum of the lens, and the subscript G M means the gravitomagnetic part. We introduce the dimensionless spin parameter as: s a / M . Hence, Equation (124) is rearranged as:
δ α G M O s M r S 2 + s M r R 2 .
This implies that δ α G M is of the same order as the second post-Newtonian effect (with the dimensionless spin parameter).
The second-order Schwarzschild contribution to α is 15 π M 2 / 4 b 2 . This contribution can be obtained also by using the present method, especially by using a relation between b and r 0 in M 2 in calculating ϕ R S . Appendix A provides detailed calculations at the second order of M and a. We explain detailed calculations for the integrals of K and κ g in the present formulation. Note that δ α G M in the above approximations is free from the impact parameter b. We can see this fact from Figure 11, Figure 12 and Figure 13 below.

8.8. Possible Astronomical Applications

What are possible astronomical applications? As a first example, we consider the Sun, in which its higher multipole moments are ignored for simplicity. Its spin angular momentum denoted as J is ∼ 2 × 10 41 m 2 kg s 1 [72,73]. This means G J c 2 5 × 10 5 m 2 , for which the dimensionless spin parameter becomes s 10 1 .
Here, our assumption is that a receiver on the Earth observes the light deflected by the Sun, while the distant source is safely in the asymptotic region. For the light ray passing near the Sun, Equation (125) allows us to make an order-of-magnitude estimation of the finite-distance correction. The result is:
δ α G M O J r R 2 10 12 arcsec . × J J 1 AU r R 2 ,
where 4 M / R 1.75 arcsec . 10 5 r a d . , where M o d o t means the solar mass and R denotes the solar radius. This correction is nearly a pico-arcsecond. Therefore, the correction is beyond the reach of present and near-future technology [74,75].
Figure 11 shows the finite-distance correction to the light deflection. Our numerical calculations are consistent with the above order-of-magnitude estimation. This figure shows also the very weak dependence of δ α on b.
See Figure 12 and Figure 13 for the deflection angle with finite-distance corrections for the prograde motion and retrograde one, respectively, where we choose r S 1.5 × 10 8 km and r R . The finite-distance correction reduces the deflection angle of light. As the impact parameter b increases, the finite-distance correction also increases.
As a second example, we discuss Sgr A* that is located at our galactic center. This object is a good candidate for measuring the strong gravitational deflection of light. The distance to the receiver is much larger than the impact parameter of light. On the other hand, some of the source stars may live in our galactic center.
For Sgr A*, Equation (125) becomes:
δ α G M s M r S 2 10 7 arcsec . × s 0.1 M 4 × 10 6 M 2 0.1 p c r S 2 ,
where we assume that the mass of the central black hole is M 4 × 10 6 M . This correction is nearly at a sub-microarcsecond level. Therefore, it is beyond the capability of present technology (e.g., [31,32,33,34,35,36]).
See Figure 9 for the finite-distance correction due to the source location. The result in this figure is in agreement with the above order-of-magnitude estimation. This figure suggests the very weak dependence on the impact parameter b.

9. Rotating Teo Wormhole: Another Example

9.1. Rotating Teo Wormhole and Optical Metric

In this section, we consider a rotating Teo wormhole [76] in order to examine how our method can be applied to a wormhole spacetime. The spacetime metric for this wormhole is:
d s 2 = N 2 d t 2 + d r 2 1 b 0 r + r 2 H 2 d θ 2 + sin 2 θ ( d ϕ ω d t ) 2 ,
where we denote:
N = H = 1 + d ( 4 a ¯ cos θ ) 2 r ,
ω = 2 a ¯ r 3 .
Here, b 0 means the throat radius of this wormhole, a ¯ is corresponding to the spin angular momentum, and d is a positive constant.
For the rotating Teo wormhole of Equation (128), the components of the generalized optical metric are [43]:
γ i j d x i d x j = r 7 ( r b 0 ) r 4 4 a ¯ 2 sin 2 θ 16 d a ¯ 2 cos 2 θ + r 2 d r 2 + r 6 r 4 4 a ¯ 2 sin 2 θ d θ 2 + r 10 sin 2 θ r 4 4 a ¯ 2 sin 2 θ 2 d ϕ 2 .
Here, γ i j is not the induced metric in the Arnowitt-Deser-Misner(ADM) formulation. The components of β i are obtained as:
β i d x i = 2 a ¯ r 3 sin 2 θ r 4 4 a ¯ 2 sin 2 θ d ϕ .
In this section, we restrict ourselves within the equatorial plane, namely θ = π / 2 . On the equatorial plane, the constant d in the metric always vanish because d is always associated with cos θ .
We employ the same way for the Kerr case; we first derive the orbit Equation on the equatorial plane from Equation (77) as:
d r d ϕ 2 = r 5 ( b 0 r ) 4 a ¯ 2 b 2 4 a ¯ b r 3 b 2 r 4 + r 6 4 a ¯ 2 b + 2 a ¯ r 3 + b r 4 2 = r 4 b 2 r 2 b 0 r 3 b 2 + b 0 r 4 a ¯ r 3 b 3 + 4 a ¯ b 0 r 2 b 3 + O ( a ¯ 2 / b 2 ) ,
where b denotes the impact parameter of the light ray and we use the weak field and slow rotation approximations in the last line. There are no b 0 squared terms in the last line. The orbit Equation thus becomes:
d u d ϕ 2 = 1 b 2 u 2 b 0 u b 2 + b 0 u 3 4 a ¯ u b 3 4 a ¯ b 0 u 2 b 3 + O ( a ¯ 2 / b 6 ) .
This Equation is iteratively solved as:
u = sin ϕ b + cos 2 ϕ 2 b 2 b 0 2 b 3 a ¯ + O b 0 2 b 3 , a ¯ b 0 b 4 .
Solving Equation (135) for ϕ S and ϕ R , we obtain ϕ S and ϕ R as:
ϕ S = arcsin ( b u S ) b 0 1 b 2 u S 2 2 b + 2 a ¯ b 2 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 ,
ϕ R = π arcsin ( b u R ) + b 0 1 b 2 u R 2 2 b 2 a ¯ b 2 1 b 2 u R 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .

9.2. Gaussian Curvature

In the weak-field approximation, the Gaussian curvature of the equatorial plane is:
K = b 0 2 r 3 56 a ¯ 2 r 6 + O a ¯ 2 b 0 r 7 , a ¯ 4 r 10 ,
where a ¯ and b 0 play roles as book-keeping parameters in the weak-field approximation. It is not surprising that this Gaussian curvature deviates from Equation (26) in Jusufi and Övgün [77], because their Gaussian curvature describes a different surface that is defined by using the Randers–Finsler metric. The Randers–Finsler metric is quite different from our generalized optical metric γ i j .
When we perform the surface integral of the Gaussian curvature in Equation (88), we use Equation (135) for a boundary of the integration domain. The surface integral of the Gaussian curvature in Equation (88) is thus calculated as:
R S K d S = ϕ S ϕ R r ( ϕ ) b 0 2 r 2 d r d ϕ + O b 0 2 b 2 , a ¯ b 0 b 3 = b 0 2 ϕ S ϕ R 0 sin ϕ b + cos 2 ϕ 2 b 2 b 0 2 b 3 a ¯ d u d ϕ + O b 0 2 b 2 , a ¯ b 0 b 3 = b 0 2 ϕ S ϕ R sin ϕ b d ϕ + O b 0 2 b 2 , a ¯ b 0 b 3 = b 0 2 cos ϕ b ϕ = ϕ S ϕ R + O b 0 2 b 2 , a ¯ b 0 b 3 = b 0 2 b 1 b 2 u R 2 + 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 ,
where we use sin ϕ R = b u R + O ( a ¯ b 2 , b 0 b 1 ) and sin ϕ S = b u S + O ( a ¯ b 2 , b 0 b 1 ) by Equations (137) and (136) in the last line.

9.3. Geodesic Curvature of Photon Orbit

We study the geodesic curvature of the photon orbit on the equatorial plane in the stationary and axisymmetric spacetime by using the generalized optical metric. It generally becomes [42]:
κ g = 1 γ γ θ θ β ϕ , r .
In the Teo wormhole, this expression is rearranged as:
κ g = 2 a ¯ r 3 + a ¯ b 0 r 4 + a ¯ b 0 2 4 r 5 + a ¯ b 0 3 8 r 6 + O a ¯ 3 r 7 , a ¯ 3 b 0 r 8 .
We compute the path integral of the geodesic curvature of the photon orbit. The detailed calculations and result are:
S R κ g d = R S 2 a ¯ r 3 d + O b 0 2 b 2 , a ¯ b 0 b 3 = π / 2 ϕ R π / 2 ϕ S 2 a ¯ cos ϑ b 2 d ϑ + O b 0 2 b 2 , a ¯ b 0 b 3 = 2 a ¯ b 2 sin π 2 ϕ S sin π 2 ϕ R + O b 0 2 b 2 , a ¯ b 0 b 3 = 2 a ¯ b 2 1 b 2 u S 2 + 1 b 2 u R 2 + O b 0 2 b 2 , a ¯ b 0 b 3 ,
for the retrograde orbit of the photon. In the last line, we used sin ϕ R = b u R + O ( a ¯ b 2 , b 0 b 1 ) and sin ϕ S = b u S + O ( a ¯ b 2 , b 0 b 1 ) from Equation (135). The above result becomes 4 a ¯ / b 2 , as r R and r S . The sign of the right-hand side in Equation (142) is opposite if the photon is in prograde motion.

9.4. ϕ R S Part

The rotating Teo wormhole is an asymptotically flat spacetime, as seen from Equation (128). Therefore, the integral of the geodesic curvature of the circular arc segment with an infinite radius can be expressed simply as ϕ R S . By using Equations (136) and (137), ϕ R S is obtained as:
ϕ R S = ϕ R ϕ S = π arcsin ( b u R ) arcsin ( b u S ) + b 0 1 b 2 u R 2 2 b + b 0 1 b 2 u S 2 2 b 2 a ¯ b 2 1 b 2 u R 2 2 a ¯ b 2 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .

9.5. Ψ Parts

For the rotating Teo wormhole in Equation (128), Equation (85) is computed as:
sin Ψ R = b u R + 2 a ¯ u R 2 4 a ¯ 2 b u R 5 ,
and Equation (86) becomes:
sin ( π Ψ S ) = b u S + 2 a ¯ u S 2 4 a ¯ 2 b u S 5 ,
where the slow-rotation approximation is not needed.
Therefore, we obtain Ψ R and Ψ S as:
Ψ R = arcsin ( b u R ) + 2 a ¯ u R 2 1 b 2 u R 2 + 2 a ¯ 2 b u R 5 2 b 2 u R 2 1 b 2 u R 2 1 3 / 2 + O ( a ¯ 3 / b 6 ) ,
π Ψ S = arcsin ( b u S ) + 2 a ¯ u S 2 1 b 2 u S 2 + 2 a ¯ 2 b u S 5 2 b 2 u S 2 1 b 2 u S 2 1 3 / 2 + O ( a ¯ 3 / b 6 ) ,
where we used the slow-rotation approximation.

9.6. Deflection Angle of Light

We combine Equations (139) and (142) to obtain the deflection angle of light in the prograde orbit as:
α prog = b 0 2 b 1 b 2 u R 2 + 1 b 2 u S 2 2 a ¯ b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .
The deflection angle of the retrograde light is:
α retro = b 0 2 b 1 b 2 u R 2 + 1 b 2 u S 2 + 2 a ¯ b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .
Next, by using Equations (143), (146), and (147), we obtain the deflection angle of the prograde light as:
α prog = π arcsin ( b u R ) arcsin ( b u S ) + b 0 1 b 2 u R 2 2 b + b 0 1 b 2 u S 2 2 b 2 a ¯ b 2 1 b 2 u R 2 2 a ¯ b 2 1 b 2 u S 2 + arcsin ( b u R ) + 2 a ¯ u R 2 1 b 2 u R 2 π + arcsin ( b u S ) + 2 a ¯ u S 2 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 = b 0 2 b 1 b 2 u R 2 + 1 b 2 u S 2 2 a ¯ b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .
The deflection angle of light in the retrograde orbit is:
α retro = b 0 2 b 1 b 2 u R 2 + 1 b 2 u S 2 + 2 a ¯ b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .
The deflection of light in the prograde (retrograde) orbit is weaker (stronger) when increasing the angular momentum of the Teo wormhole. The reason is as follows. The local inertial frame in which the light travels at the light speed c in general relativity moves faster (slower). Hence, the time-of-flight of light becomes shorter (longer). On light propagation, a similar explanation is done by using the dragging of the inertial frame also by Laguna and Wolsczan [78]. They discussed the Shapiro time delay. The expression of the deflection angle of light by a rotating Teo wormhole is similar to that by Kerr black hole. This implies that it is hard to distinguish a Kerr black hole from a rotating Teo wormhole by the gravitational lens observations.
In Equations (150) and (151), the source and receiver can be located at finite distances from the wormhole. In the limit as r R and r S , Equations (148) and (149) become:
α prog b 0 b 4 a ¯ b 2 + O b 0 2 b 2 , a ¯ b 0 b 3 , α retro b 0 b + 4 a ¯ b 2 + O b 0 2 b 2 , a ¯ b 0 b 3 .
They are in complete agreement with Equations (39) and (56) in Jusufi and Övgün [77], where they restrict themselves within the asymptotic source and receiver ( r R and r S ).

9.7. Finite-Distance Corrections in the Teo Wormhole Spacetime

To be precise, we define the finite-distance correction to the deflection angle of light as the difference between the asymptotic deflection angle α and the deflection angle for the finite distance case. It is denoted as δ α .
We consider the following situation. An observer on the Earth sees the light deflected by the solar mass. The source of light is located in a practically asymptotic region. In other words, we choose b 0 = M , a ¯ = J , r R 1.5 × 10 8 km, r S . See Figure 14 for the finite-distance correction due to the impact parameter b. In Figure 14, the green curve means the difference between the asymptotic bending angle and the deflection angle with finite-distance corrections, the blue curve denotes the asymptotic deflection angle, and the orange curve is the deflection angle with finite-distance corrections by the rotating Teo wormhole. The deflection angle is decreased by the finite-distance correction. If the impact parameter b increases, the finite-distance correction also increases.
See also Figure 15 for numerical calculations of the finite-distance correction due to the impact parameter b. In Figure 15, the blue curve is the deflection angle with finite-distance correction by a Kerr black hole and the red curve is the deflection angle with finite-correction by a rotating Teo wormhole. The deflection of light is stronger in a Kerr black hole case for the chosen values.

10. Summary

In this paper, we provided a brief review of a series of works on the deflection angle of light for a light source and receiver in a non-asymptotic region. [38,39,42,43]. The validity and usefulness of the new formulation come from the GB theorem in differential geometry. First, we discussed how to define the gravitational deflection angle of light in a static, spherically symmetric, and asymptotically flat spacetime, for which we assume the finite-distance source and receiver. We examined whether our definition is invariant geometrically by using the GB theorem. By using our definition, we carefully computed finite-distance corrections to the light deflection in Schwarzschild spacetime. We considered both the cases of weak deflection and the strong one. Next, we extended the definition to stationary and axisymmetric spacetimes. This extension allows us to compute finite-distance corrections for Kerr black holes and rotating Teo wormholes. We verified that these results are consistent with previous works in the infinite-distance limit. We mentioned also the finite-distance corrections to the light deflection by Sagittarius A*. It is left as future work to apply the present formulation to other interesting spacetime models and also to extend it to a more general spacetime structure.

Funding

This work was supported in part by JSPS research fellowship for young researchers (T.O.); in part by Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research, No. 18J14865 (T.O.) and No. 17K05431 (H.A.); and in part by by the Ministry of Education, Culture, Sports, Science, and Technology, No. 17H06359 (H.A.).

Acknowledgments

We are grateful to Marcus Werner for the stimulating and very fruitful discussions. We thank Takao Kitamura, Asahi Ishihara, and Yusuke Suzuki for useful conversations. We would like to thank Yuuiti Sendouda, Ryuichi Takahashi, Yuya Nakamura, and Naoki Tsukamoto for useful conversations.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Detailed Calculations at O(M2/b2) and O(a2/b2) in Kerr Spacetime

First, we investigate the Gaussian curvature K of the equatorial plane in the Kerr spacetime. Here, we assume the weak-field and slow-rotation approximations. Up to the second order, K is expanded as:
K = R r ϕ r ϕ γ = 2 M r 3 + 3 M 2 r 4 + O a 2 M r 5 ,
where γ denotes det ( γ i j ) . There are no a 2 terms in K. More interestingly, only the a 2 M term at the third-order level exists in K. By noting that K begins with O ( M ) , what we need for the second-order calculations is only the linear-order term in the area element on the equatorial plane. This is obtained as:
d S γ d r d ϕ r + 3 M + O M 2 r d r d ϕ ,
where terms at O ( a ) and at O ( a 2 ) do not exist in d S . This is because all terms including the spin parameter cancel out in γ for θ = π / 2 and γ thus depends only on M, as shown by direct calculations.
By using Equations (A1) and (A2), the surface integration of the Gaussian curvature is performed as:
K d S = r O E d r ϕ S ϕ R d ϕ 2 M r 3 + 3 M 2 r 4 ( r + 3 M ) + O M 3 b 3 , a M 2 b 3 , a 2 M b 3 = 0 1 b sin ϕ + M b 2 ( 1 + cos 2 ϕ ) d u ϕ S ϕ R d ϕ ( 2 M + 3 u M 2 ) + O M 3 b 3 , a M 2 b 3 , a 2 M b 3 = ϕ S ϕ R 2 M b sin ϕ + M 2 2 b 2 ( 7 + cos 2 ϕ ) d ϕ + O M 3 b 3 , a M 2 b 3 , a 2 M b 3 = 2 M b cos ϕ ϕ R ϕ S + M 2 2 b 2 30 ϕ + sin ( 2 ϕ ) 4 ϕ S ϕ R + O M 3 b 3 , a M 2 b 3 , a 2 M b 3 = 2 M b 1 b 2 u S 2 + 1 b 2 u R 2 + 15 M 2 4 b 2 [ π arcsin ( b u S ) arcsin ( b u R ) ] + M 2 4 b 2 [ b u S ( 15 7 b 2 u S 2 ) 1 b 2 u S 2 + b u R ( 15 7 b 2 u R 2 ) 1 b 2 u R 2 ] + O M 3 b 3 , a M 2 b 3 , a 2 M b 3 ,
where we use, in the second line, an iterative solution for the orbit equation of Equation (77) in the Kerr spacetime.
Next, we study the geodesic curvature. On the equatorial plane, we find:
κ g = 1 Σ 2 Δ ( Σ 2 M r ) r 2 + a 2 + 2 a 2 M r sin 2 θ Σ Σ sin 2 θ ( Σ 2 M r ) β ϕ , r = 2 a M r 3 + O a M 2 r 3 .
Note that a 2 terms do not exist. Therefore, we obtain:
c p κ g d = S R d 2 a M r 2 + O a M 2 r 3 = 2 a M b 2 ϕ S ϕ R cos ϑ d ϑ + O a M 2 b 3 = 2 a M b 2 [ 1 b 2 u R 2 + 1 b 2 u S 2 ] + O a M 2 b 3 ,
where we use sin ϕ S = r S 2 b 2 / r S + O ( M / r S ) and sin ϕ R = r R 2 b 2 / r R + O ( M / r R ) .
By combining Equations (A3) and (A5), we obtain:
α R S K d S R S κ g d = 2 M b 1 b 2 u S 2 + 1 b 2 u R 2 + 15 M 2 4 b 2 π arcsin ( b u S ) arcsin ( b u R ) + M 2 4 b 2 b u S ( 15 7 b 2 u S 2 ) 1 b 2 u S 2 + b u R ( 15 7 b 2 u R 2 ) 1 b 2 u R 2 2 a M b 2 1 b 2 u R 2 + 1 b 2 u S 2 + O M 3 b 3 , a M 2 b 3 , a 2 M b 3 .
Note that a 2 terms and a 3 ones do not appear in α for the finite distance situation as well as in the infinite distance limit. If we assume the infinite distance limit u R , u S 0 , Equation (A6) becomes:
α 4 M b + 15 π M 2 4 b 2 4 a M b 2 .
This agrees with the known results, especially on the numerical coefficients at the order of M 2 and a M .

References

  1. Einstein, A. Die Grundlage der allgemeinen Relativitatstheorie. Ann. Phys. (Berlin) 1916, 49, 769–822. [Google Scholar] [CrossRef]
  2. Dyson, F.W.; Eddington, A.S.; 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. Philos. Trans. R. Soc. A 1920, 220, 291. [Google Scholar] [CrossRef]
  3. Hagihara, Y. Theory of the relativistic trajectories in a gravitational field of Schwarzschild. Jpn. J. Astron. Geophys. 1931, 8, 67. [Google Scholar]
  4. Chandrasekhar, S. The Mathematical Theory of Black Holes; Oxford University Press: New York, NY, USA, 1998. [Google Scholar]
  5. Misner, C.W.; Thorne, K.S.; Wheeler, J.A. Gravitation; Freeman: New York, NY, USA, 1973. [Google Scholar]
  6. Darwin, C. The gravity field of a particle. Proc. R. Soc. A 1959, 249, 180. [Google Scholar]
  7. Bozza, V. Gravitational lensing in the strong field limit. Phys. Rev. D 2002, 66, 103001. [Google Scholar] [CrossRef] [Green Version]
  8. Iyer, S.V.; Petters, A.O. Light’s bending angle due to black holes: From the photon sphere to infinity. Gen. Relativ. Gravit. 2007, 39, 1563. [Google Scholar] [CrossRef]
  9. Bozza, V.; Scarpetta, G. Strong deflection limit of black hole gravitational lensing with arbitrary source distances. Phys. Rev. D 2007, 76, 083008. [Google Scholar] [CrossRef] [Green Version]
  10. Frittelli, S.; Kling, T.P.; Newman, E.T. Spacetime perspective of Schwarzschild lensing. Phys. Rev. D 2000, 61, 064021. [Google Scholar] [CrossRef] [Green Version]
  11. Virbhadra, K.S.; Ellis, G.F.R. Schwarzschild black hole lensing. Phys. Rev. D 2000, 62, 084003. [Google Scholar] [CrossRef] [Green Version]
  12. Virbhadra, K.S. Relativistic images of Schwarzschild black hole lensing. Phys. Rev. D 2009, 79, 083004. [Google Scholar] [CrossRef] [Green Version]
  13. Virbhadra, K.S.; Narasimha, D.; Chitre, S.M. Role of the scalar field in gravitational lensing. Astron. Astrophys. 1998, 337, 1. [Google Scholar]
  14. Virbhadra, K.S.; Ellis, G.F.R. Gravitational lensing by naked singularities. Phys. Rev. D 2002, 65, 103004. [Google Scholar] [CrossRef]
  15. Virbhadra, K.S.; Keeton, C.R. Time delay and magnification centroid due to gravitational lensing by black holes and naked singularities. Phys. Rev. D 2008, 77, 124014. [Google Scholar] [CrossRef] [Green Version]
  16. Zschocke, S. A generalized lens equation for light deflection in weak gravitational fields. Class. Quantum Gravity 2011, 28, 125016. [Google Scholar] [CrossRef]
  17. Eiroa, E.F.; Romero, G.E.; Torres, D.F. Reissner-Nordstrom black hole lensing. Phys. Rev. D 2002, 66, 024010. [Google Scholar] [CrossRef]
  18. Perlick, V. Exact gravitational lens equation in spherically symmetric and static spacetimes. Phys. Rev. D 2004, 69, 064017. [Google Scholar] [CrossRef] [Green Version]
  19. Abe, F. Gravitational Microlensing by the Ellis Wormhole. Astrophys. J. 2010, 725, 787. [Google Scholar] [CrossRef]
  20. Toki, Y.; Kitamura, T.; Asada, H.; Abe, F. Astrometric Image Centroid Displacements due to Gravitational Microlensing by the Ellis Wormhole. Astrophys. J. 2011, 740, 121. [Google Scholar] [CrossRef]
  21. Nakajima, K.; Asada, H. Deflection angle of light in an Ellis wormhole geometry. Phys. Rev. D 2012, 85, 107501. [Google Scholar] [CrossRef]
  22. Gibbons, G.W.; Vyska, M. The Application of Weierstrass elliptic functions to Schwarzschild Null Geodesics. Class. Quant. Grav. 2012, 29, 065016. [Google Scholar] [CrossRef]
  23. DeAndrea, J.P.; Alexander, K.M. Negative time delay in strongly naked singularity lensing. Phys. Rev. D 2014, 89, 123012. [Google Scholar] [CrossRef]
  24. Kitamura, T.; Nakajima, K.; Asada, H. Demagnifying gravitational lenses toward hunting a clue of exotic matter and energy. Phys. Rev. D 2013, 87, 027501. [Google Scholar] [CrossRef]
  25. Tsukamoto, N.; Harada, T. Signed magnification sums for general spherical lenses. Phys. Rev. D 2013, 87, 024024. [Google Scholar] [CrossRef]
  26. Izumi, K.; Hagiwara, C.; Nakajima, K.; Kitamura, T.; Asada, H. Gravitational lensing shear by an exotic lens object with negative convergence or negative mass. Phys. Rev. D 2013, 88, 024049. [Google Scholar] [CrossRef]
  27. Kitamura, T.; Izumi, K.; Nakajima, K.; Hagiwara, C.; Asada, H. Microlensed image centroid motions by an exotic lens object with negative convergence or negative mass. Phys. Rev. D 2014, 89, 084020. [Google Scholar] [CrossRef]
  28. Nakajima, K.; Izumi, K.; Asada, H. Negative time delay of light by a gravitational concave lens. Phys. Rev. D 2014, 90, 084026. [Google Scholar] [CrossRef] [Green Version]
  29. 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, 064043. [Google Scholar] [CrossRef] [Green Version]
  30. Azreg-Ainou, M. Confined-exotic-matter wormholes with no gluing effects – Imaging supermassive wormholes and black holes. J. Cosmol. Astropart. Phys. 2015, 07, 037. [Google Scholar] [CrossRef]
  31. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball1, D.; Baloković, M.; Barrett, J.; Bintley, D.; et al. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole. Astrophys. J. 2019, 875, L1. [Google Scholar]
  32. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball1, D.; Baloković, M.; Barrett, J.; Bintley, D.; et al. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. II. Array and Instrumentation. Astrophys. J. 2019, 875, L2. [Google Scholar]
  33. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball1, D.; Baloković, M.; Barrett, J.; Bintley, D.; et al. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. III. Data Processing and Calibration. Astrophys. J. 2019, 875, L3. [Google Scholar]
  34. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball1, D.; Baloković, M.; Barrett, J.; Bintley, D.; et al. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole. Astrophys. J. 2019, 875, L4. [Google Scholar]
  35. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball1, D.; Baloković, M.; Barrett, J.; Bintley, D.; et al. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring. Astrophys. J. 2019, 875, L5. [Google Scholar]
  36. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball1, D.; Baloković, M.; Barrett, J.; Bintley, D.; et al. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole. Astrophys. J. 2019, 875, L6. [Google Scholar]
  37. Gibbons, G.W.; Werner, M.C. Applications of the Gauss-Bonnet theorem to gravitational lensing. Class. Quant. Grav. 2008, 25, 235009. [Google Scholar] [CrossRef]
  38. Ishihara, A.; Suzuki, Y.; Ono, T.; Kitamura, T.; Asada, H. Gravitational bending angle of light for finite distance and the Gauss-Bonnet theorem. Phys. Rev. D 2016, 94, 084015. [Google Scholar] [CrossRef] [Green Version]
  39. Ishihara, A.; Suzuki, Y.; Ono, T.; Asada, H. Finite-distance corrections to the gravitational bending angle of light in the strong deflection limit. Phys. Rev. D 2017, 95, 044017. [Google Scholar] [CrossRef] [Green Version]
  40. Arakida, H. Light deflection and Gauss-Bonnet theorem: Definition of total deflection angle and its applications. Gen. Rel. Grav. 2018, 50, 48. [Google Scholar] [CrossRef]
  41. Crisnejo, G.; Gallo, E.; Rogers, A. Finite distance corrections to the light deflection in a gravitational field with a plasma medium. Phys. Rev. D 2019, 99, 124001. [Google Scholar] [CrossRef] [Green Version]
  42. Ono, T.; Ishihara, A.; Asada, H. Gravitomagnetic bending angle of light with finite-distance corrections in stationary axisymmetric spacetimes. Phys. Rev. D 2017, 96, 104037. [Google Scholar] [CrossRef] [Green Version]
  43. Ono, T.; Ishihara, A.; Asada, H. Deflection angle of light for an observer and source at finite distance from a rotating wormhole. Phys. Rev. D 2018, 98, 044047. [Google Scholar] [CrossRef] [Green Version]
  44. Ovgun, A. Light deflection by Damour-Solodukhin wormholes and Gauss-Bonnet theorem. Phys. Rev. D 2018, 98, 044033. [Google Scholar] [CrossRef] [Green Version]
  45. Ono, T.; Ishihara, A.; Asada, H. Deflection angle of light for an observer and source at finite distance from a rotating global monopole. Phys. Rev. D 2019, 99, 124030. [Google Scholar] [CrossRef] [Green Version]
  46. Carmo, M.P.D. Differential Geometry of Curves and Surfaces; Prentice-Hall: Upper Saddle River, NJ, USA, 1976; pp. 268–269. [Google Scholar]
  47. Kottler, F. Uber die physikalischen Grundlagen der Einsteinschen Gravitationstheorie. Annalen. Phys. 1918, 361, 401–462. [Google Scholar] [CrossRef]
  48. Bach, R. Zur Weylschen Relativitatstheorie und der Weylschen Erweiterung des Krummungstensorbegriffs. Math. Zeit. 1921, 9, 110–135. [Google Scholar] [CrossRef]
  49. Riegert, R.J. Birkhoff’s Theorem in Conformal Gravity. Phys. Rev. Lett. 1984, 53, 315. [Google Scholar] [CrossRef]
  50. Mannheim, P.D.; Kazanas, D. Exact vacuum solution to conformal Weyl gravity and galactic rotation curves. Astrophys. J. 1989, 342, 635–638. [Google Scholar] [CrossRef]
  51. Edery, A.; Paranjape, M.B. Classical tests for Weyl gravity: Deflection of light and time delay. Phys. Rev. D 1998, 58, 024011. [Google Scholar] [CrossRef] [Green Version]
  52. Sultana, J.; Kazanas, D. Bending of light in conformal Weyl gravity. Phys. Rev. D 2010, 81, 127502. [Google Scholar] [CrossRef]
  53. Cattani, C.; Scalia, M.; Laserra, E.; Bochicchio, I.; Nandi, K.K. Correct light deflection in Weyl conformal gravity. Phys. Rev. D 2013, 87, 047503. [Google Scholar] [CrossRef]
  54. Sereno, M. Role of Λ in the Cosmological Lens Equation. Phys. Rev. Lett. 2009, 102, 021301. [Google Scholar] [CrossRef] [PubMed]
  55. Bhadra, A.; Biswas, S.; Sarkar, K. Gravitational deflection of light in the Schwarzschild–de Sitter space-time. Phys. Rev. D 2010, 82, 063003. [Google Scholar] [CrossRef]
  56. Arakida, H.; Kasai, M. Effect of the cosmological constant on the bending of light and the cosmological lens equation. Phys. Rev. D 2012, 85, 023006. [Google Scholar] [CrossRef]
  57. Lim, Y.; Wang, Q. Exact gravitational lensing in conformal gravity and Schwarzschild–de Sitter spacetime. Phys. Rev. D 2017, 95, 024004. [Google Scholar] [CrossRef]
  58. Lewis, T. Some Special Solutions of the Equations of Axially Symmetric Gravitational Fields. Proc. Roy. Soc. A 1932, 136, 176. [Google Scholar] [CrossRef]
  59. Levy, H.; Robinson, W.J. The rotating body problem. Proc. Camb. Philos. Soc. 1964, 60, 279. [Google Scholar] [CrossRef]
  60. Papapetrou, A. Champs gravitationnels stationnaires a symetrie axiale. Ann. Inst. H. Poincare A 1966, 4, 83–105. [Google Scholar]
  61. Levi-Civita, T. Absolute Differential Calculus; Blackie and Son: Glasgow, UK, 1927. [Google Scholar]
  62. Asada, H.; Kasai, M. Can We See a Rotating Gravitational Lens? Prog. Theor. Phys. 2000, 104, 95–102. [Google Scholar] [CrossRef] [Green Version]
  63. Belton, A.C. Geometry of Curves and Surfaces. 2015, p. 38. Available online: www.maths.lancs.ac.uk/~belton/www/notes/geom_notes.pdf (accessed on 10 June 2019).
  64. Oprea, J. Differential Geometry and Its Applications, 2nd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2003; p. 210. [Google Scholar]
  65. Perlick, V. Ray Optics, Fermat’s Principle, and Applications to General Relativity; Springer: Berlin, Germany, 2000. [Google Scholar]
  66. Kopeikin, S.; Mashhoon, B. Gravitomagnetic effects in the propagation of electromagnetic waves in variable gravitational fields of arbitrary-moving and spinning bodies. Phys. Rev. D 2002, 65, 064025. [Google Scholar] [CrossRef] [Green Version]
  67. Werner, M.C. Gravitational lensing in the Kerr-Randers optical geometry. Gen. Rel. Grav. 2012, 44, 3047–3057. [Google Scholar] [CrossRef] [Green Version]
  68. Epstein, R.; Shapiro, I.I. Post-post-Newtonian deflection of light by the Sun. Phys. Rev. D 1980, 22, 2947. [Google Scholar] [CrossRef]
  69. Ibanez, J. Gravitational lenses with angular momentum. Astron. Astrophys. 1983, 124, 175–180. [Google Scholar]
  70. Iyer, S.V.; Hansen, E.C. Light’s bending angle in the equatorial plane of a Kerr black hole. Phys. Rev. D 2009, 80, 124023. [Google Scholar] [CrossRef]
  71. Kraniotis, G.V. Precise analytic treatment of Kerr and Kerr-(anti) de Sitter black holes as gravitational lenses. Class. Quant. Grav. 2011, 28, 085021. [Google Scholar] [CrossRef]
  72. Pijpers, F.P. Helioseismic determination of the solar gravitational quadrupole moment. Mon. Not. R. Astron. Soc. 1998, 297, L76. [Google Scholar] [CrossRef]
  73. Bi, S.L.; Li, T.D.; Li, L.H.; Yang, W.M. Solar Models with Revised Abundance. Astrophys. J. Lett. 2011, 731, L42. [Google Scholar] [CrossRef]
  74. Gaia. Available online: http://sci.esa.int/gaia/ (accessed on 10 June 2019).
  75. JASMINE. Available online: http://www.jasmine-galaxy.org/index-en.html (accessed on 10 June 2019).
  76. Teo, E. Rotating traversable wormholes. Phys. Rev. D 1998, 58, 024014. [Google Scholar] [CrossRef] [Green Version]
  77. Jusufi, K.; Övgün, A. Exact traversable wormhole solution in bumblebee gravity. Phys. Rev. D 2018, 97, 024042. [Google Scholar] [CrossRef]
  78. Laguna, P.; Wolszczan, A. Pulse Arrival Times from Binary Pulsars with Rotating Black Hole Companions. Astrophys. J. 1997, 486, L27. [Google Scholar] [CrossRef]
Figure 1. Gauss–Bonnet theorem: We consider a closed curve in a surface.
Figure 1. Gauss–Bonnet theorem: We consider a closed curve in a surface.
Universe 05 00218 g001
Figure 2. R S is a quadrilateral embedded in a curved space.
Figure 2. R S is a quadrilateral embedded in a curved space.
Universe 05 00218 g002
Figure 3. Observable α in Equation (7): In this schematic figure, the lens, receiver, and source are respectively the Sun, the Earth, and a pulsar that periodically radiates radio signals in a specific anistropic manner. From the pulse profile, we can determine the radiation direction at the source. By using the ephemeris, we know the relative positions of the Sun, Earth, and the pulsar. Hence, we can determine ϕ R S and Ψ S . By observing the pulsar, we can measure Ψ R . In principle, therefore, we can determine Ψ R Ψ S + ϕ R S from these astronomical observations.
Figure 3. Observable α in Equation (7): In this schematic figure, the lens, receiver, and source are respectively the Sun, the Earth, and a pulsar that periodically radiates radio signals in a specific anistropic manner. From the pulse profile, we can determine the radiation direction at the source. By using the ephemeris, we know the relative positions of the Sun, Earth, and the pulsar. Hence, we can determine ϕ R S and Ψ S . By observing the pulsar, we can measure Ψ R . In principle, therefore, we can determine Ψ R Ψ S + ϕ R S from these astronomical observations.
Universe 05 00218 g003
Figure 4. A one-loop case for the photon orbit in M opt .
Figure 4. A one-loop case for the photon orbit in M opt .
Universe 05 00218 g004
Figure 5. Quadrilaterals: They are made from the photon orbit in a non-Euclidean space. See Figure 4.
Figure 5. Quadrilaterals: They are made from the photon orbit in a non-Euclidean space. See Figure 4.
Universe 05 00218 g005
Figure 6. Two loops for the light ray in M opt .
Figure 6. Two loops for the light ray in M opt .
Universe 05 00218 g006
Figure 7. Quadrilaterals (1)–(4): They are in a non-Euclidean plane M opt . See also Figure 6.
Figure 7. Quadrilaterals (1)–(4): They are in a non-Euclidean plane M opt . See also Figure 6.
Universe 05 00218 g007
Figure 8. A quadrilateral in any loop number: This case is discussed when we prove by induction that Equation (7) holds in any loop number.
Figure 8. A quadrilateral in any loop number: This case is discussed when we prove by induction that Equation (7) holds in any loop number.
Universe 05 00218 g008
Figure 9. The finite-distance correction for Sgr A* as δ α G M : The horizontal axis denotes the source distance r S . The vertical one means the finite-distance correction to the light deflection. The solid line (blue in color) and dashed one (red in color) mean b = 10 2 M and b = 10 4 M , respectively. The dotted curve (yellow in color) denotes the leading term of δ α G M given by Equation (33). These three lines are substantially overlapped with each other. This implies that δ α G M is weakly dependent on the impact parameter b.
Figure 9. The finite-distance correction for Sgr A* as δ α G M : The horizontal axis denotes the source distance r S . The vertical one means the finite-distance correction to the light deflection. The solid line (blue in color) and dashed one (red in color) mean b = 10 2 M and b = 10 4 M , respectively. The dotted curve (yellow in color) denotes the leading term of δ α G M given by Equation (33). These three lines are substantially overlapped with each other. This implies that δ α G M is weakly dependent on the impact parameter b.
Universe 05 00218 g009
Figure 10. Ψ R and Ψ S : Ψ R is the angle between the radial direction and the light ray at the receiver position. Ψ S is that at the source position.
Figure 10. Ψ R and Ψ S : Ψ R is the angle between the radial direction and the light ray at the receiver position. Ψ S is that at the source position.
Universe 05 00218 g010
Figure 11. δ α G M for the Sun: The horizontal axis is the distance of the receiver distance r R . The vertical axis means the finite-distance correction due to the gravitomagnetic deflection angle of light. The solid curve (blue in color) and dashed one (red in color) denote b = R and b = 10 5 R , respectively. The dotted line (black in color) corresponds to the leading term in δ α G M given by Equation (124). These three curves are overlapped. This implies the very weak dependence of δ α G M on b.
Figure 11. δ α G M for the Sun: The horizontal axis is the distance of the receiver distance r R . The vertical axis means the finite-distance correction due to the gravitomagnetic deflection angle of light. The solid curve (blue in color) and dashed one (red in color) denote b = R and b = 10 5 R , respectively. The dotted line (black in color) corresponds to the leading term in δ α G M given by Equation (124). These three curves are overlapped. This implies the very weak dependence of δ α G M on b.
Universe 05 00218 g011
Figure 12. α in the prograde motion: The horizontal axis is the impact parameter for a photon orbit. The vertical axis means the deflection angle of light. The blue curve is the asymptotic deflection angle by a Kerr black hole. The orange curve means the deflection angle with finite-corrections by a Kerr black hole. The green curve shows the difference between the asymptotic bending angle and the deflection angle with finite-corrections by a Kerr black hole.
Figure 12. α in the prograde motion: The horizontal axis is the impact parameter for a photon orbit. The vertical axis means the deflection angle of light. The blue curve is the asymptotic deflection angle by a Kerr black hole. The orange curve means the deflection angle with finite-corrections by a Kerr black hole. The green curve shows the difference between the asymptotic bending angle and the deflection angle with finite-corrections by a Kerr black hole.
Universe 05 00218 g012
Figure 13. α for light of retrograde motion: The horizontal axis denotes the impact parameter for a photon orbit, and the vertical axis denotes the deflection angle of light. The blue curve is the asymptotic deflection angle by the Kerr black hole. The orange curve is the deflection angle with finite-correction by the Kerr black hole. The green curve shows the difference between the asymptotic bending angle and the deflection angle with finite-correction by the Kerr black hole.
Figure 13. α for light of retrograde motion: The horizontal axis denotes the impact parameter for a photon orbit, and the vertical axis denotes the deflection angle of light. The blue curve is the asymptotic deflection angle by the Kerr black hole. The orange curve is the deflection angle with finite-correction by the Kerr black hole. The green curve shows the difference between the asymptotic bending angle and the deflection angle with finite-correction by the Kerr black hole.
Universe 05 00218 g013
Figure 14. α in the Teo wormhole: The blue curve is the asymptotic deflection angle by the rotating Teo wormhole. The orange curve is the deflection angle with finite-distance corrections by the rotating Teo wormhole. The blue curve shows the difference between the asymptotic deflection angle and the deflection angle with finite-distance corrections by the rotating Teo wormhole.
Figure 14. α in the Teo wormhole: The blue curve is the asymptotic deflection angle by the rotating Teo wormhole. The orange curve is the deflection angle with finite-distance corrections by the rotating Teo wormhole. The blue curve shows the difference between the asymptotic deflection angle and the deflection angle with finite-distance corrections by the rotating Teo wormhole.
Universe 05 00218 g014
Figure 15. α for prograde motion of light: The horizontal axis is the impact parameter of photon orbit. The vertical axis means the deflection angle of light. The blue curve means the deflection angle with finite-distance corrections by the Kerr black hole. The red curve corresponds to that by the rotating Teo wormhole. For the purpose of this comparison, the mass of a Kerr black hole M and the throat radius of a rotating Teo wormhole b 0 are chosen as M = b 0 = M . The spin angular momentum of a Kerr black hole and that of a rotating Teo wormhole are chosen as the same as that the Sun for simplicity.
Figure 15. α for prograde motion of light: The horizontal axis is the impact parameter of photon orbit. The vertical axis means the deflection angle of light. The blue curve means the deflection angle with finite-distance corrections by the Kerr black hole. The red curve corresponds to that by the rotating Teo wormhole. For the purpose of this comparison, the mass of a Kerr black hole M and the throat radius of a rotating Teo wormhole b 0 are chosen as M = b 0 = M . The spin angular momentum of a Kerr black hole and that of a rotating Teo wormhole are chosen as the same as that the Sun for simplicity.
Universe 05 00218 g015

Share and Cite

MDPI and ACS Style

Ono, T.; Asada, H. The Effects of Finite Distance on the Gravitational Deflection Angle of Light. Universe 2019, 5, 218. https://doi.org/10.3390/universe5110218

AMA Style

Ono T, Asada H. The Effects of Finite Distance on the Gravitational Deflection Angle of Light. Universe. 2019; 5(11):218. https://doi.org/10.3390/universe5110218

Chicago/Turabian Style

Ono, Toshiaki, and Hideki Asada. 2019. "The Effects of Finite Distance on the Gravitational Deflection Angle of Light" Universe 5, no. 11: 218. https://doi.org/10.3390/universe5110218

APA Style

Ono, T., & Asada, H. (2019). The Effects of Finite Distance on the Gravitational Deflection Angle of Light. Universe, 5(11), 218. https://doi.org/10.3390/universe5110218

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop