Application of the Gauss-Bonnet theorem to lensing in the NUT metric

We show with the help of Fermat's principle that every lightlike geodesic in the NUT metric projects to a geodesic of a two-dimensional Riemannian metric which we call the optical metric. The optical metric is defined on a (coordinate) cone whose opening angle is determined by the impact parameter of the lightlike geodesic. We show that, surprisingly, the optical metrics on cones with different opening angles are locally (but not globally) isometric. With the help of the Gauss-Bonnet theorem we demonstrate that the deflection angle of a lightlike geodesic is determined by an area integral over the Gaussian curvature of the optical metric. A similar result is known to be true for static and spherically symmetric spacetimes. The generalisation to the NUT spacetime, which is neither static nor spherically symmetric (at least not in the usual sense), is rather non-trivial.

and l in the following. m detemines the mass of the central object, whereas l determines its gravitomagnetic charge. In an analogy to electromagnetism, m corresponds to the electric charge whereas l, which is also known as the NUT parameter, corresponds to a magnetic (monopole) charge.
The NUT metric has so many unusual properties that Misner [2] called it a "counter-example to almost anything". This, however, does not necessarily mean that the NUT metric is unphysical. On the contrary, history has taught us to take solutions to Einstein's vacuum equation seriously, even if they have apparently "exotic" properties. So we believe that the existence in Nature of sources with a non-zero NUT parameter is a possibility that should not be rejected off-hand.
Following this line of thought, we have to ask ourselves what kind of observable features a hypothetical NUT source would have. In essence, there are two types of such features. Firstly, we could consider the influence of such a source on the motion of massive particles that come close to it, secondly we could consider its influence on light rays. For the discussion of observable effects of a NUT source on massive particles we refer in particular to papers by Hackmann and Lämmerzahl [3] and Jefremov and Perlick [4]. We will not discuss this subject here. As to the influence on light rays, i.e., to the lensing features of a NUT source, there is a paper by Nouri-Zonoz and Lynden Bell [5]. It is the purpose of the present paper to develop their analysis somewhat further.
Our main goal is to find out whether the lensing features in the NUT spacetime can be characterised in terms of two-dimensional Riemannian geometry. It is well known that such a characterisation is possible in spacetimes that are static and spherically symmetric: In this case it suffices to consider lightlike geodesics in the equatorial plane ϑ = π/2 because all the other ones are then determined by the spherical symmetry. Here (t, r, ϑ, ϕ) are the spherical polar coordinates in which a static and spherically symmetric spacetime is usually given. Then one finds that the spatial paths, in the two-dimensional manifold with coordinates (r, ϕ), of lightlike geodesics are the geodesics of a Riemannian metric which is known as the Fermat metric or the optical metric. Gibbons and Werner [6] have used the Gauss-Bonnet theorem to show how the optical metric determines lensing features. In particular, they have demonstrated that the deflection angle is determined by the Gaussian curvature of the optical metric. In this paper we want to investigate if and how this result carries over to the NUT spacetime. As the NUT spacetime is not static but only stationary, and as in the NUT spacetime geodesics are not contained in a plane, this makes several non-trivial modifications necessary.
The paper is organised as follows. In Sec. 2 and Sec. 3 we summarise some basic features of the NUT spacetime and its geodesics. Throughout we do this for a version of the NUT metric that is more general than the one introduced in the original NUT paper [1] because it depends not only on the above-mentioned parameters m and l but also on a third parameter, C, that was introduced in 2005 by Manko and Ruiz [7]. In Sec. 4 and Sec. 5 we make use of the fact that in the NUT spacetime every lightlike geodesic is contained in a cone whose opening angle is determined by the impact parameter of the lightlike geodesic and we calculate the deflection angle. For the original NUT metric with C = −1 the fact that every geodesic is contained in a cone was derived already by Zimmerman and Shahir [8]; this fact plays a crucial role also in the above-mentioned paper by Nouri-Zonoz and Lynden-Bell [5] who consider the NUT metric with C = 0. In Sec. 6 we use Fermat's principle to define a Riemannian metric, which we call the optical metric, on each of these cones and we demonstrate that a lightlike geodesic with the corresponding impact parameter projects, indeed, to a geodesic of the optical metric. We illustrate the optical metric by way of an embedding diagram, we demonstrate that it has negative Gaussian curvature and we show that the optical metrics on cones with different opening angles are locally isometric. Finally, in Sec. 7 we use the Gauss-Bonnet theorem to rewrite the deflection angle as an integral over the Gaussian curvature of the optical metric. As the latter is negative, this result demonstrates that all light rays are deflected towards the centre and that the deflection angle increases with decreasing impact parameter.
We use the metric signature (-,+,+,+) and we set the speed of light c = 1. We use Einstein's summation convention with greek indices for the four spacetime coordinates and with latin indices for the two coordinates of the two-dimensional manifold on which the optical metric lives.

The NUT metric
The NUT metric [1] is a solution to Einstein's vacuum field equation. It depends on two parameters both of which have the dimension of a length, a mass parameter m and a gravitomagnetic charge, also known as the NUT parameter, l. We assume m > 0 and −∞ < l < ∞ throughout. For l = 0 the NUT metric reduces to the Schwarzschild metric. In 2005 Manko and Ruiz [7] brought forward a generalised version of the NUT metric which involves a third parameter, C, that is dimensionless and may take any value −∞ < C < ∞. With the Manko-Ruiz parameter included, the metric reads Here, (t, r, ϑ, ϕ) are Boyer-Lindquist-type coordinates. We assume that t ranges over all of R, ϑ and ϕ are standard coordinates on the two-sphere, and the radius coordinate is restricted to the domain m + √ m 2 + l 2 < r < ∞. At r = m + √ m 2 + l 2 the metric features a black-hole horizon. If analytically extended beyond the horizon, in the domain m− √ m 2 + l 2 < r < m+ √ m 2 + l 2 the metric is isometric to a cosmological solution that was found already in 1951 by Taub [9]. The entire spacetime is therefore known as the Taub-NUT spacetime. In this paper, however, we will not consider the Taub region.
By a coordinate transformation t = t − 2lCϕ, r = r, ϕ = ϕ, ϑ = ϑ the metric (1) is transformed to the corresponding metric with C = 0. This, however, is not a globally well-defined transformation: As we assume that ϑ and ϕ are standard coordinates on the sphere, ϕ is 2π-periodic; therefore, the transformation is well-defined only locally, on any neighbourhood that does not include an entire ϕ-line, unless we make the time coordinate periodic with period 4π|lC|. As the latter would lead to the most drastic causality violation, giving a closed timelike curve through each point of the considered domain of the spacetime, we will not do this. Therefore, NUT spacetimes with different values of C are locally but not globally isometric.
The original NUT metric [1] is the metric (1) with C = −1. It features a conic singularity on the half-axis ϑ = π which was interpreted by Bonnor [10] as a spinning rod. For C = 1 the singularity is on the other half-axis ϑ = 0, and for all other values of C it is on both half-axes, symmetrically distributed for C = 0 and asymmetrically for C = 0. For more details on the physical and geometrical meaning of the parameter C we refer to Manko and Ruiz [7] and also to Jefremov and Perlick [4].
For any values of m, l and C the NUT metric admits a four-dimensional Lie algebra of Killing vector fields, i.e., a NUT metric with l = 0 has as many symmetries as the Schwarzschild metric. For C = −1 the Killing vector fields have been given already in the original NUT article [1] and the resulting symmetry properties have been discussed in detail by Misner [11] who advocated the idea of making the t coordinate periodic. Here we give the four linearly independent Killing vector fields for the NUT metric with an arbitrary value of C: The Lie brackets are This demonstrates that ξ 1 , ξ 2 and ξ 3 generate a three-dimensional group of isometries which is isomorphic to the rotation group SO(3, R), whereas ξ 0 generates a one-dimensional group of isometries that expresses stationarity. Note that for l = 0 the orbits of the rotations are not two-dimensional spacelike spheres; they are rather three-dimensional submanifolds with topology R × S 2 and the metric has the signature (− + +) on these submanifolds. For this reason many authors find it inappropriate to call the NUT metric "spherically symmetric". It is, however, safe to call it "rotationally symmetric about any radial direction" or to say that "it admits an SO(3, R) symmetry".

Geodesics in the NUT metric
With the Killing vector fields known, it is easy to solve the geodesic equation in the NUT spacetime. For any geodesic x µ (s), parametrised by an affine parameter s, each of the four Killing vector fields ξ A gives us a constant of motion g µν ξ µ Aẋ ν : Hence and The latter equation, in which the central dot means the standard Euclidean scalar product in R 3 , demonstrates that a geodesic with constants of motion E, J 1 , J 2 and J 3 is contained in a cone whose vertex is at the coordinate origin, whose symmetry axis is spanned by the vector with Cartesian components J 1 , J 2 and J 3 and whose opening angle α is given by Note that this statement refers to a three-dimensional Euclidean space that is defined by the chosen coordinates. We will later discuss the intrinsic geometries of these coordinate cones. For the NUT metric with C = −1, the fact that each geodesic is contained in a cone was discovered already by Zimmerman and Shahir [8]. (Note, however, that in their paper the g tt component of the metric is misprinted.) Eq. (14) demonstrates that the opening angle of this cone is independent of C if the geodesic is labelled by the constant of motion J/E. Of course, for l = 0 (14) gives α = π/2, i.e., it reproduces the well-known fact that in the Schwarzschild spacetime each geodesic is contained in a plane through the coordinate origin.
In addition to the constants of motion that arise from the Killing symmetries, also the Lagrangian is a constant of motion. With the help of the four constants of motion E, J 2 , J 3 and L the geodesic equation can be written in first-order form. To that end we have to solve, in this order, (11), (8), (12) and (15) forφ,ṫ,θ 2 andṙ 2 which results inφ = Here J 2 − 4 l 2 E 2 is the Carter constant which arises as a separation constant if the geodesic equation is written in Hamilton-Jacobi form. Note that the spatial paths of the geodesics, which are determined by eqs. (16), (18) and (19), are not affected by the Manko-Ruiz parameter C.

Lightlike geodesics contained in a cone with ϑ = constans
In the Schwarzschild spacetime it suffices to consider geodesics with J 1 = J 2 = 0, which are exactly the geodesics in the equatorial plane ϑ = π/2. If one has determined all such geodesics, one gets all the other ones by applying all possible rotations, i.e., the action of the group generated by the Killing vector fields ξ 1 , ξ 2 and ξ 3 . An analogous statement holds for the NUT metric: Also in this case it suffices to consider geodesics with J 1 = J 2 = 0 which now are the geodesics contained in cones of the form ϑ = constans; the opening angle α of such a cone equals the coordinate angle ϑ. Again, one gets all the other geodesics by applying to such geodesics all possible rotations. If J 1 = J 2 = 0 and ϑ = constans, (9) and (10) requirė and (13) yields The latter equation demonstrates that, in a NUT spacetime with l = 0, the opening angle ϑ of the cone is determined by the impact parameter J 3 /E and vice versa. In the Schwarzschild case l = 0 eqs. (20) and (21) determine neitheṙ ϕ nor J 3 /E because with l also cos ϑ goes to zero. For lightlike geodesics (L = 0), inserting (21) into (19) results iṅ Differentiating this equation with respect to the affine parameter and dividing by 2ṙ givesr Although we have divided byṙ, this equation is valid by continuity also iḟ r = 0. By solving simultaneously the equationsṙ = 0 andr = 0 we see that there is a circular lightlike geodesic at the intersection of the sphere r = r ph with the cone ϑ = ϑ ph where and It is easy to see that the cubic equation (24) has three real solutions exactly one of which lies in the considered domain, m + √ m 2 + l 2 < r ph . Applying all possible rotations demonstrates that the sphere at radius coordinate r ph is filled with circular lightlike geodesics. We refer to it as to the photon sphere. The existence of the photon sphere is crucial for determining the socalled shadow of a NUT black hole, see Grenzebach et al. [20] where this is discussed for a class of spacetimes that contains the NUT spacetime as a special case.
From (20) and (22) we get the orbit equation that determines the shape of a lightlike geodesic in the cone ϑ = constans, From this equation we read that, for l = 0, there are no lightlike geodesics contained in the equatorial plane ϑ = π/2 because the expression under the square-root must be non-negative.

The deflection angle
We want to consider, in the cone ϑ = constans, a lightlike geodesic that comes in from infinity, goes through a minimum radius value r = r m and then escapes back to infinity. Note that such a light ray may make arbitrarily many turns around the centre. At r = r m , the right-hand side of (22) must be zero, and the right-hand side of (23) must be negative. Therefore, lightlike geodesics that go through a minimal radius value r m exist for all r m outside of the photon sphere, r m > r ph , and the corresponding values of ϑ converge towards ϑ ph for r m → r ph . If we define, for fixed ϑ, the function then r m is determined by the equation F (r m ) = 0, where the largest positive solution to this fourth-order equation is the relevant one. In Fig. 1 the function F (r) is plotted for a fixed value l = 0 and various values of ϑ. This figure confirms our earlier observation that for l = 0 there are no light rays in the equatorial plane ϑ = π/2. Fig. 2 shows a plot of r m as a function of l for various values of ϑ.
. (29) Integration of the orbit equation (29) over the light ray from its point of closest approach to infinity yields .
For defining the deflection angle of a light ray we introduce a new azimuthal coordinateφ which is the angle defined around the surface of the cone. On each circle r = constans on the cone, ϕ runs from 0 to 2π whereasφ runs from 0 to 2π sin ϑ, i.e., over a smaller interval. This deficit angle can be visualised by cutting the cone open and flattening it. The deflection angle, or bending angle, is defined as the angle under which the asymptotes to the light ray intersect in the cut and flattened cone, see Fig. 3, where ∆φ = sin ϑ ∆ϕ .
Here we have used (27). Our definition of the deflection angle is in agreement with that of Nouri-Zonoz and Lynden-Bell [5]. Note, however, that in their paper there is a notational inconsistency. In the beginning they define the NUT metric (with C = 0) in the same coordinates as we do: Their r is our r and their φ is our ϕ. But later, from their eqs. (4) and (5) onwards, their r is our √ r 2 + l 2 and their φ is ourφ. For light rays that make many turns around the centre, the deflection angle becomes arbitrarily large, i.e. δ → ∞ for r m → r ph . If, however, r m is big, δ is small and we get a valid approximation of δ if we do a Taylor expansion of a low order with respect to the dimensionless parameters m/r m and l/r m . We have to go at least up to the second order for having a non-vanishing influence of l. Then (33) can be written as where O(3) stands for terms of third or higher order in m/r m and l/r m . As a consequence, the bending angle to within this order is In linear order we have the same result as in the Schwarzschild spacetime which is clear from the fact that l enters quadratically in all relevant equations.

Optical metric and Gaussian curvature
It is known that in an arbitrary spacetime lightlike geodesics satisfy a variational principle that can be viewed as a general-relativistic Fermat principle, see Temple [12] for a version restricted to a local normal neighbourhood, Kovner [13] for a formulation of the general principle and Perlick [14] for a complete proof. For static or stationary spacetimes, there are simpler and older versions of this variational principle. In the static case it was shown already in 1917 by Weyl [16] that, as a consequence of Fermat's principle, the spatial paths of lightlike geodesics are the geodesics of a Riemannian metric which is now known as the Fermat metric or the optical metric. Levi-Civita [17] generalised Weyl's result to the case of a stationary spacetime. Then the spatial paths of lightlike geodesics are determined by the combined action of a Riemannian metric and a one-form which is sometimes called the Fermat one-form; for examples we refer to Perlick [15]. Because of the influence of the Fermat one-form, the lightlike geodesics do not project to geodesics of a Riemannian metric. They rather project to geodesics of a Finsler metric. Although, in hindsight, this relation to Finsler geometry is fairly obvious, to the best of our knowledge it was realised only recently, namely by Caponio, Javaloyes and Masiello [18] in a paper that was put on the arXiv in 2007. It is Levi-Civita's version of Fermat's principle that we will apply to the NUT metric now. According to the foregoing results we may restrict our consideration to a cone ϑ = constans. From (1) we read that along every lightlike curve (not necessarily a geodesic) the differential of the coordinate time along its spatial projection equals where i, j ∈ {r, ϕ}. Here β i andḡ ij are given by If one fixes two points in the cone and integrates dt over lightlike curves that project to curves in the cone between these two fixed points, then the integral over the one-form β i dx i gives the same value for all these curves, By Levi-Civita's version of Fermat's principle, for an actual lightlike geodesic the variation of T vanishes, i.e.
This demonstrates that the spatial path of the lightlike geodesic is a geodesic of the Riemannian metricḡ ij defined in (38). We call it the "optical metric" henceforth. Here we have to keep in mind that we have to choose the opening angle ϑ of the cone as determined by the impact parameter of the lightlike geodesic according to (21). As the NUT metric is not static, but only stationary, it was not to be expected that the lightlike geodesics project to geodesics of a Riemannian metric.
We have already emphasised that calling the surface ϑ = constans a "cone" refers to the chosen coordinates. The optical metric, which has a coordinateindependent meaning, makes this surface into a two-dimensional Riemannian manifold. In order to visualise its intrinsic geometry we try to isometrically embed it into Euclidean 3-space. In cylindrical polar coordinates (Z, R, ϕ), we have to satisfy the condition with embedding functions Z(r) and R(r). After inserting (38) and comparing coefficients of dϕ 2 and dr 2 we find and dZ(r) dr (43) An isometric embedding into Euclidean 3-space is possible in that part of the domain m + √ m 2 + l 2 < r < ∞ where dZ(r)/dr 2 is non-negative. The boundary value, r = r b , of the embeddable part is the greatest zero of the sixth-order polynomial in the numerator on the right-hand of (43). In the Schwarzschild limit l = 0, which requires ϑ = π/2, we recover the well-known result r b = 9m/4. From Fig. 4 we read that the intrinsic geometry of the coordinate cone ϑ = constans becomes that of a Euclidean cone only asymptotically for r → ∞. It develops a "neck" at the intersection with the photon sphere r = r ph , then opens out again before the boundary of the embeddable part is reached at r = r b . We also read from the picture that at each point the principal curvatures have opposite signs, i.e., that the Gaussian curvature is negative which means that geodesics of the optical metric must locally diverge. This is true not only for the embeddable part but everywhere on the domain m+ √ m 2 + l 2 < r < ∞ on which the optical metric is defined. To demonstrate this we calculate the Gaussian curvature of the optical metric analytically. By definition, the Gaussian curvature K of a two-dimensional Riemannian metricḡ ij satisfies the equation whereR ijkl are the covariant components of the Riemannian curvature tensor. For a metric of the form (38), the Gaussian curvature can then be calculated as which gives, with the special metric coefficients from (38), For a plot of K as a function of r, for different NUT-parameters, see Fig. 5. Eq. (46) confirms that, for any value of l, the Gaussian curvature of the optical metric is indeed negative on the entire domain m + √ m 2 + l 2 < r < ∞. To prove this, substitute r = m + √ m 2 + l 2 + x; then (l 2 + r 2 ) 4 K becomes a fifth-order polynomial in x whose coefficients are manifestly negative. Also, (46) implies that for r → m + √ m 2 + l 2 the Gaussian curvature K approaches a strictly negative value and its derivative dK/dr approaches zero. Here it is important that we consider a NUT spacetime with m > 0. The observation that K is negative is also true for m = 0 as long as l = 0. For m = 0 and l = 0, we have Minkowski spacetime and the cones are Euclidean cones, i.e, they are locally flat with K = 0.
Remarkably, K is independent of ϑ. That is to say, on all cones with their different opening angles the Gaussian curvature depends on r in exactly the same way. A careful look at (38) shows that this is, actually, not so surprising: If one changes from the coordinates (r, ϕ) to (r = r,φ = ϕ sin ϑ), the metric coefficients become independent of ϑ. Hence it is clear that the optical geometries of any two cones with different opening angles are locally isometric. They are, however, not globally isometric because the range of the coordinatẽ ϕ depends on ϑ. By the same token, their embedding diagrams depend on ϑ because in the ambient Euclidean space the azimuthal coordinate is assumed to be 2π-periodic. If, on the other hand, we use the representation in a plane where (r,φ) are the polar coordinates, as in Fig. 3, then the optical metrics are represented by a metric on this plane that is independent of ϑ. We can thus show all the geodesics of the optical metrics, i.e., all lightlike geodesics, in one and the same plane; we just have to keep in mind that the deficit angle, which is marked in grey in Fig. 3, is different for geodesics with different impact parameters.

Gravitational lensing and Gauss-Bonnet theorem
For the Schwarzschild spacetime, and any other spherically symmetric and static spacetime, the Gauss-Bonnet theorem can be used for characterising some lensing features. This was demonstrated in a pioneering paper by Gibbons and Werner [6] who in particular derived from the Gauss-Bonnet theorem a formula for the bending angle in such spacetimes. Their analysis was based on the well-known facts that in a spherically symmetric and static spacetime it suffices to consider lightlike geodesics in the equatorial plane ϑ = π/2 and that the spatial projections of lightlike geodesics in this plane are the geodesics of a two-dimensional Riemannian metric. Werner [19] demonstrated that a similar analysis is possible for lightlike geodesics in the equatorial plane of the Kerr metric. There are several differences in comparison to the spherically symmetric and static case. Firstly, in the Kerr metric the lightlike geodesics do not project to geodesics of a Riemannian metric; therefore Werner had to use the idea of an "osculating Riemannian metric" for achieving his goal. Secondly, in the Kerr metric confining oneself to the equatorial plane is, of course, a strong restriction of generality. We will now demonstrate that in the NUT spacetime the Gauss-Bonnet theorem can be applied to all lightlike geodesics, without any loss of generality, and that there is no need to introduce an osculating metric or any other fiducial quantities. We follow the original idea of Gibbons and Werner as closely as possible. For background material on the Gauss-Bonnet theorem the reader may consult any text-book on Riemannian geometry, e.g. Klingenberg [21]. In its standard text-book version, the Gauss-Bonnet theorem is valid for a compact domain with boundary in a two-dimensional Riemannian manifold. For applying this theorem to lightlike geodesics in the NUT spacetime, we use for the two-dimensional Riemannian manifold the cone ϑ = constans with the optical metric (38). For the compact domain we use a quadrangle D R whose boundary consists of four smooth parts, see Fig. 6. As in Fig. 3, also in this picture the cone is shown cut open and flattened, i.e., what is shown is a diagram using (r,φ) as polar coordinates.
The four smooth parts of the boundary of D R are constructed in the following way. The first part is a finite section of the projection of a lightlike geodesic, denoted γ in Fig. 6; we assume that γ has no self-intersections and that it comes from infinity, goes through a minimal radius value r = r m and then escapes to infinity. Along the section of γ that is part of the boundary of D R the angleφ runs over an interval of lengthψ. Another part of the boundary of D R is a section of a circle r = R, denoted γ R in the picture. The two remaining parts are radial. Now the Gauss-Bonnet theorem says that Here K is the Gaussian curvature and dS is the area element of the optical metric. ∂D R is the boundary of D R and κ is the geodesic curvature of this boundary, to be integrated over the four smooth parts of ∂D R where s denotes arclength with respect to the optical metric. The α A are the four jump angles marked in Fig. 6 and χ(D R ) is the Euler characteristic of D R . As D R is connected and has no holes, χ(D R ) = 1. Moreover, obviously α 1 = α 2 = π/2. The integral over the boundary of D R reduces to an integral over the circular arc because the other three parts are geodesic, hence κ = 0 there. We parametrise the circular arc by arclength, i.e., we write it in the form γ R (s) with g γ R ,γ R ) = 1. The geodesic curvature of γ R is then to be calculated as where ∇ is the Levi-Civita connection of the optical metric g and n is a unit vector perpendicular toγ R . This puts (47) into the following form:
(50) Therefore, the limit version of the Gauss-Bonnet theorem reads Comparison with (32) yields As K is negative, we read from this formula that the deflection angle in the NUT spacetime is always positive, i.e., that every light ray is deflected towards the centre. This is not obvious if the deflection angle is represented with the help of (33) because one does not know whether the integral on the right-hand side of this equation is bigger or smaller than π/2 before one has actually calculated it. Moreover, as K is independent of ϑ, it is also evident from (52) that a second light ray has a smaller deflection angle than the first one if on its entire path it stays farther away from the centre. Again, this is not obvious from the representation with (33) because the minimal radius r m occurs not only as the lower limit of the integral but also in the integrand. When deriving (52) we have assumed that the light ray has no self-intersection, i.e., we excluded the case that the light ray makes a full turn, or several full turns, around the centre. If one wants to include such cases, one has to consider a region D ∞ that is bounded by the parts of the light ray between infinity and the outermost self-intersection point, one has to add 2nπ sin ϑ on the right-hand side of (52) where n is the number of full turns performed by the light ray, and one has to take the jump angle at the self-intersection point into account. In particular the last point makes the use of the Gauss-Bonnet theorem to light rays with self-intersections rather awkward.
To make (52) more explicit, we write the area element of the optical metric as The bending angle now is where r = z(φ) is the polar representation of the deflected light ray, with the angleφ as the parameter. Werner [19] has calculated the bending angle in the equatorial plane of the Kerr metric to within linear order in the mass parameter and the spin parameter. In the NUT metric we have to go at least up to second order if we want to have a non-zero contribution of the NUT parameter. We have calculated the bending angle already up to this order in (35). As a cross-check, we want to reproduce this result with the Gauss-Bonnet theorem. To that end we need the integrand in (54) up to second order, As this expression has a vanishing zeroth-order term, we need ∆φ and z(φ) only up to first order. ∆φ was calculated in (34), For calculating the orbit r = z(φ), we integrate on the right-hand side of (34) not from r m but from z(φ) to infinity, where we have assumed that the light ray passes atφ = 0 through the point of closest approach. Solving for z(φ) results in z(φ) = r m cosφ 1 − 2 + cosφ sin 2φ 1 + cosφ cosφ Inserting (55), (56) and (58) into (54) reproduces, indeed, (35). We emphasise that the line integral in (33) is easier to evaluate than the area integral in (52). Therefore, if the only goal is to actually calculate the deflection angle, then there is no point in using the Gaus-Bonnet theorem. The merit of the latter is in the fact that it immediately gives some qualitative lensing features in terms of geometric quantities.

Discussion and conclusions
When Gibbons and Werner [6] applied the Gauss-Bonnet theorem to the spatial paths of lightlike geodesics in static and spherically symmetric spacetimes, many readers found this idea attractive because it related the lensing features in such spacetimes to geometric quantities such as the Gaussian curvature of a two-dimensional Riemannian metric. Naturally, the question arises if a similar result can be found in spacetimes that are not static and spherically symmetric. Werner [19] considered the Kerr metric and he demonstrated that, at least for light rays in the equatorial plane, the Gauss-Bonnet theorem is applicable. However, in contrast to the static and spherically symmetric case, the light rays do not project to geodesics of a Riemannian metric but rather to geodesics of a Finsler metric of Randers type; this made it necessary to introduce a socalled "osculating metric" to which the standard Gauss-Bonnet theorem could be applied.
Here we have shown that the original method by Gibbons and Werner, without a restriction to special lightlike geodesics and without the need of introducing an osculating metric, can be applied to the NUT metric. The main difference to the cases considered earlier is in the fact that now the spatial projection of each lightlike geodesic is in a cone, rather than in a plane. Throughout we have allowed for an arbitrary value of the Manko-Ruiz parameter C, and we have verified that C has no influence on the spatial paths of lightlike geodesics. Our analysis revealed two facts that could not have been anticipated before the calculation was done: Firstly, we showed that the projection of every lightlike geodesic is indeed a geodesic of an optical metric that is Riemannian; so there is no need of considering Finsler metrics. Secondly, on cones with different opening angles the optical metrics turned out to be locally isometric; in particular, the Gaussian curvature of the optical metric, as a function of the radius coordinate, turned out to be independent of the opening angle of the cone. Together with the observation that the Gaussian curvature of the optical metric is negative, this allowed us to use the Gauss-Bonnet theorem for representing qualitative lensing features of the NUT spacetime in a geometric way that is as significant as in the Schwarzschild spacetime.