On the properties of a deformed extension of the NUT space-time

We consider a class of space-times given by a stationary extension of the Zipoy-Voorhees metric that was found by Halilsoy. We show that the solutions do not describe rotating sources but must be interpreted, similarly to the NUT case, as deformed sources endowed with a gravitomagnetic charge. We show that the Halilsoy family is directly linked to the NUT space-time, which can be obtained in the limit of vanishing deformations. We investigate the motion of test particles and photons in this class of space-times, in particular the innermost stable circular orbits and photon capture radius. Finally we show that this class of solutions possesses a sub-manifold where closed time-like curves are allowed.


I. INTRODUCTION
Vacuum, exact solutions describing stationary spacetimes are of great importance in General Relativity since they can describe the field outside a rotating compact source. The most famous of such solutions is the well known Kerr metric which describes the field of a rotating black hole [1].
However other solutions do exist and investigating their physical properties is important in order to establish whether astrophysical black hole candidates are well described by mathematical black hole solutions [2]. An important class of exact solutions of Einstein's equations that describes deviations from black hole space-times is the so-called Weyl class of static axially symmetric vacuum solutions [3]. Given the one-to-one correspondence between metrics belonging to Weyl's class and solutions of Laplace equation in flat two-dimensional space, all static axially symmetric solutions are in principle known [4]. For example in recent times some attention has been given to the Zipoy-Voorhees (ZV) metric which is a static generalization of the Schwarzschild solution to include higher multipole moments and describes the field outside prolate or oblate spheroids [5]. The properties of the motion of test particles in the ZV space-time and the possibility of testing the geometry from astrophysical observations has been discussed in several articles [6]. However, astrophysical compact objects typically rotate and therefore it would be more interesting to study the properties of stationary solutions. This has led many authors in the past to consider 'rotating' generalization of static axially symmetric solutions (see for example [7]). In particular, stationary generalizations of the ZV metric have been studied in [8]. However, as it turns out, not * nbakhtiyor18@fudan.edu.cn † daniele.malafarina@nu.edu.kz ‡ ahmadjon@astrin.uz § bambi@fudan.edu.cn all such generalizations describe rotating objects. There exist a number of solution generating techniques that can be adopted in order to obtain a stationary solution from a known static one [9]. Such techniques have been employed in the past to derive a large number of 'rotating' generalizations of known solutions. However, the physical properties of these new solutions have rarely been investigated. One notable exception, i.e. a stationary solution that has been thoroughly investigated, is the socalled Newman-Unti-Tamburino (NUT) space-time [10], which, due to its peculiar properties, has been referred to as the 'counterexample to everything' [11]. Similarly to Kerr, the NUT space-time is fully characterized by two parameters, one related to the mass of the source, and the other, the so-called NUT parameter, related to the gravitomagnetic charge of the source. However, and differently from the Kerr case, the NUT parameter does not characterise the rotation of the source.
In the present article we consider another stationary space-time that was obtained by Halilsoy in [12] and investigate its properties. The stationary solution was constructed using the ZV space-time as the 'seed' metric. We show that, contrary to what is claimed in the title of Halilsoy's paper, the line element does not describe a spinning massive source. On the contrary, the 'spin' parameter in the solution behaves similarly to the NUT parameter and we show explicitly that the metric in question reduces to the NUT space-time in the case of vanishing deformations.
The paper is organised as follows: In section II we outline the basic features of the stationary extension of the ZV space-time, show its relation with the NUT spacetime and discuss its properties as the exterior of a gravitating compact source. Section III is devoted to the study of the motion of test particles in the stationary ZV geometry, with particular emphasis on its relation to the Schwarzschild, Kerr and NUT space-times. Finally in section IV we briefly discuss the results and their implications for astrophysical black holes and exotic compact objects.
Throughout the paper we make use of geometrized units setting G = c = 1.
Specific solutions of this class can be obtained from known static axially-symmetric solutions through a variety of procedures [9]. However, it is well known that not all stationary solutions of the above class describe the exterior gravitational field of a rotating object. The most famous of such examples is the Newmann-Unti-Tamburino (NUT) metric [10].
The NUT space-time belongs to the above class and its metric functions are given by where M and l are positive parameters and In the case l = 0 the metric becomes static and it reduces to the Schwarzschild geometry. From the asymptotic expansion of the metric the parameter M is immediately interpreted as the gravitational mass of the source. However, the interpretation of the NUT solution and of the NUT parameter l is less trivial and it has attracted lots of attention over the years. In [13] the solution was interpreted as representing the exterior field of a mass located at the origin together with a semi-infinite massless source of angular momentum. On the other hand in [14] the metric was interpreted as "describing the exterior field of two counterrotating semiinfinite sources possessing negative masses and infinite angular momenta". In [15] the parameter was interpreted as a gravitomagnetic charge bestowed upon the central mass. Conversely in [16] the parameter was interpreted as a property of the surrounding space-time. Another famous solution of the above class is the ZV space-time describing the gravitational field outside a static deformed object [5]. The metric functions are given by where m and γ are positive parameters and The space-time reduces to Schwarzschild for γ = 1 and the parameter γ is easily interpreted as a deformation parameter. The total gravitational mass of the source is M = mγ and from the computation of the quadrupole moment Q = γm 3 (1 − γ 2 )/3 one can see that values of γ < 1 (γ > 1) correspond to prolate (oblate) deformations. The properties of the ZV space-time were studied in [6], while interior solutions for the ZV metric were obtained in [17]. A stationary generalization of the ZV metric was obtained by Halilsoy in [12]. If we perform a change of coordinates from cylindrical to prolate spheroidal coordinates {t, x, y, φ} given by with κ = m/γ, the line element (1) takes the form and the metric functions can be given as The line element depends on three parameters m, γ and q with p given by p 2 + q 2 = 1. For q = 0 one retrieves the ZV space-time so that for q = 0 (i.e. p = 1) and γ = 1 one obtains the Schwarzschild metric. Similarly to the space-times discussed above the parameter m is related to the gravitational mass of the source, γ is related to the deformation of the source while the parameter q describes the departure from the static solution. However, one needs to be careful when interpreting q as a rotation parameter similar to the angular momentum in Kerr's solution, since a stationary space-time may be obtained also by means of introducing a NUT-like parameter. Indeed, this turns out to be the case here. In order to relate the above line element to the Schwarzschild and NUT space-times it is useful to make use of Erez-Rosen coordinates {t, r, θ, φ} using the transformations or equivalently so that the line element (1) takes the form and the metric functions become with It can be easily checked that setting q = 0 one recovers the ZV metric in the usual form given by Erez-Rosen coordinates. However, as mentioned before, for q = 0 and γ = 1 one does not retrieve the Kerr solution. In fact the line element for the case γ = 1 becomes It is not difficult to see that for p = 1 (i.e. q = 0) we retrieve the Schwarzschild solution in Schwarzschild coordinates. However, to retrieve the NUT solution in the usual form when γ = 1 and q = 0 one needs to perform the change of coordinates and identify the parameters from Then the line element (27) takes the familiar form It is easy to see that the metric (20) reduces to Minkowski space-time for any value of γ when m = 0. This also shows that Minkowski can be obtained as the zero mass limit for the NUT metric, if we set m = 0, instead of the usually discussed limit with H(r) = (r 2 − l 2 )/(r 2 + l 2 ), which is obtained by setting M = 0 while keeping l = 0. Also it is worth noticing that the coordinate transformation (29) maps the horizons of the NUT line elementr h = M ± √ M 2 + l 2 into r i = 0 and r o = 2m.
In the case of small deviation from spherical symmetry (δ = γ −1 << 1) and small deviations from staticity (q 0), the components of the metric (20) can be expanded as If we look at the asymptotic properties of the metric (20) we find the following expressions From equation (38) one can see that the gravitational mass of the source is given by M = γmp. However, from the expansion of g tφ we see how the parameter q is not related to the angular momentum of the central object, as, otherwise, it should appear as q/r at the leading order in g tφ . The above statement can be made more precise by evaluating the Komar integrals in the NUT limit, i.e. in the case γ = 1. The Komar angular momentum J K is given by (see [18]): where φ α is a space-like Killing vector, while n α and r α are the unit vectors normal to the t = const and r = const hypersurfaces, respectively. In the Kerr spacetime the Komar angular momentum is J = ma. On the other hand, in the NUT space-time the Komar angular momentum is equal to zero. The fact that the integral above for the metric (27) vanishes confirms that the parameter q does not describe rotation but is better interpreted as some NUT-like charge that we shall call for simplicity 'quasi-NUT' parameter, which is directly related to the NUT parameter in the case γ = 1 from equation (30). Similarly one can calculate the Komar mass M K of the source for the metric (27). This is given by As expected we see that the result coincides with the NUT mass parameter. Similarly to the ZV metric the Halilsoy's metric exhibits a curvature singularity for r = 2m for all values of γ = 1, as can be seen from evaluation of the Kretschmann scalar K = R αβγσ R αβγσ , which characterizes the space-time curvature. On the other hand, one can check that for γ = 1, similarly to Schwarzschild, g tt changes sign at r = 2m while g rr diverges, as can be seen from Fig. 1. Evaluation of the Kretschmann scalar in this case shows that r = 2m is regular and coincides with the horizon of the NUT spacetime. It is easy to check that for γ = 1 the Kretschmann scalar K does not diverge for any value of q = 0, thus showing that the singularity given by the divergence of g rr at r = 0 is also a coordinate singularity and the spacetime can be extended to the range r ∈ (−∞, +∞). In this case K is given by with h(r) = 8m 6 (p − 1) 3 p + 24m 5 (p − 1) 3 r + 2p 2 − 1 r 6 + +60m 4 (p − 1) 2 r 2 − 40m 3 (p − 1) 2 r 3 + +30m 2 (p − 1)pr 4 + 6m −2p 2 + p + 1 r 5 .
Also, it is easy to notice that in the coordinates used in (27) for p = 1 we have K(0) = −192m 6 p(p−1) 2 < 0, thus suggesting that repulsive effects appear in the vicinity of the center of the space-time. Further we can investigate the departure from spherical symmetry of the source of the metric (27) by calculating the area of the surfaces with constant r surrounding the source. Therefore we must use the three-dimensional line element given by [19] where Then the surface becomes It is clear that in the static case (i.e. p = 1) we recover the revolution surfaces in the Schwarzschild geometry, S Sch = 4πr 2 , while for p = 1 we obtain S < S Sch . Also for p = 1 there is a limiting radius given by at which S(r 0 ) = 0.
Finally, when thinking about the physical validity of the line element (27) it is worth checking whether pathologies such as closed time-like curves do appear. One simple way to do so is to study the properties of geodesics with constant values of the coordinates t, r and θ. The interval reduces to and it is space-like for g φφ > 0. It is easy to see that for γ = 1 and q = 0 there exist regions where g φφ turns from space-like to time-like thus showing that closed time-like curves can appear in those cases. For any fixed value of r one may find the regions where closed time-like curve are allowed by finding the zeroes of g φφ as a function of q and θ (see Fig. 2). On the other hand, for q = 0, i.e. the ZV metric, and for γ = 1, i.e. the NUT metric, we see that g φφ remains positive everywhere for positive values of r.

III. PARTICLE MOTION
In order to investigate the viability of the line element as a possible exterior gravitational field of a physical source it is worth looking at the motion of test particles. Therefore we shall now outline the formalism to describe test particle motion in the metric given by equation (20). Since the metric does not depend on the coordinates t and φ we have the two usual time-like and space-like Killing vector fields associated with time translations and rotations. Therefore the action S for the particle can be written in the usual form where ε and L are the constants of motion associated with the two Killing vectors and representing the energy and angular momentum per unit mass of the test particle. Using Hamilton-Jacobi equation we find the effective potential V ef f (r) for test particles in the equatorial plane (θ = π/2). The radial component of the four velocity of the particle readṡ and the effective potential takes the form where e ψ is evaluated at θ = π/2. The radial dependence of the effective potential for various values of γ and various values of the 'quasi-NUT' parameter p is shown in FIG. 3. In the case γ = 1 the effective potential takes the simple form and particle motion in the NUT space-time was studied in [20]. The innermost stable circular orbit (ISCO) of test particles in the equatorial plane can be calculated using the following standard conditions obtained from equation (50) [21]: Numerical values for the ISCO radius depending on γ and p are presented in TAB.I. From the table one can see that in the cases when γ = 1 and γ = 1.5 the ISCO has a minimum for a value of p ∈ (0, 1) (see right panel in Fig. 4). One can also compare the stationary ZV metric with Kerr and NUT in the usual coordinates (see Fig. 4). It is interesting to notice how the actual parametrization affects the ISCO, since it is a coordinate dependent quantity. In fact the case γ = 1 and the NUT metric represent the same space-time in different coordinates and the two plots can be made to coincide making use of the change of coordinates (29).  It is worth remembering that the NUT space-time posses circular time-like geodesics outside the equatorial plane [22]. This can be seen also in the newly introduced coordinates for the NUT space-time given in equation (27). In fact the components of the four velocity of a particle moving in the space-time described by the line element (27) arė where K is a Carter's constant, b(r) = r 2 − 2m(1 − p) + 2m 2 (1−p), and a(r) = −b(r) 2 /(r 2 −2mr). From equation (59) one can show that there is off equatorial motion on a plane parallel to the equatorial one when Since we know that in the absence of parameter q the metric must reduce to Schwarzschild we are led to choose Carter's constant as K = L 2 . Then the above expression reduces to where the minus (−) sign gives us the motion on equatorial plane while the plus (+) sign gives the orbit confined to a cone with the opening angle θ given by Let us now investigate the photon motion in the spacetime characterized by the line element (20). Using Hamilton-Jacobi equation of motion for massless particle one may derive the effective potential for photons as where and L again define the conserved energy and angular momentum of photons. The radial dependence of this effective potential is plotted in Fig. 5. One may check that for the effective potential of photons the sign of V ef f (r) is negative at points near the photon sphere showing that, as expected, photon orbits are unstable. From the symmetry of the space-time one may calculate the radius of circular photon orbits on the equatorial plane (θ = π/2) from either of the following geodesic equations x µ + Γ µ αβẋ αẋβ = 0 . Settingr =ṙ =θ = 0 we obtain following set of equations g 00 + 2g 03 dφ dt + g 33 dφ dt Solving the above equations analytically in the general case is complicated, however, some insights can be obtained from the plot of the photon capture radius in the equatorial plane as a function of γ and q, as shown in Fig. 6. It can be seen that the location of the photon capture radius increases with γ similarly to what happens in the static case. Also the dependence of the photon sphere on the 'quasi-NUT' parameter shows how the space-time differs from Kerr. However, one needs to be careful with the interpretation of the radial coordinate as a radial distance, since a simple change of coordinate like the one in equation (29) and a redefinition of the parameters like the ones given in equations (30) and (31) may completely change the behaviour of the photon sphere. Finally one may consider the appearance of the source for far away observers by evaluating the shape of the shadow of the line element (20). The shadow of the Kerr-Taub-NUT black hole was considered in [23] and the shadow of the static ZV metric was studied in [24]. The shadow of the stationary ZV metric may be related to the shadow of the static ZV metric and that of the Kerr-NUT metric once the 'quasi-NUT' parameter and the deformation parameter vanish. The shadow of the stationary ZV metric has been obtained using the ray tracing code that was developed in [25][26][27][28] and it is shown in in Fig. 7. It is clearly seen how the deformation parameter γ and the 'quasi-NUT' parameter affect the appearance of the source, thus suggesting that the nature geometry could in principle be tested via observations.

IV. CONCLUSION
We considered a stationary extension of the Zipoy-Voorhees space-time that was originally found by Halilsoy in [12] in order to investigate if it could be taken as a suitable candidate to describe the field in the exterior of an astrophysical compact object or a black hole mimicker.
We found that the line element does not describe a rotating source but rather a deformed NUT space-time that reduces to the NUT metric in the limit of vanishing deformations, i.e. γ = 1. Therefore the presence of the 'quasi-NUT' parameter q does not affect the motion of test particles in a way comparable to the Kerr parame-ter. In fact, in terms of the appearance of the source for far away observers the behaviour remains qualitatively similar to the ZV case, for all values of the 'quasi-NUT' parameter. Furthermore, for γ = 1, the metric exhibits the presence of closed time-like curves in a finite region of the space-time thus suggesting that the field may not represent the exterior of a viable source.
The effective potentials describing the motion of massive test particles and photons in the equatorial plane have been used to determine the location of the innermost stable circular orbit and the photon sphere. The dependence of these two radii on the deformation parameter and the 'quasi-NUT' parameter shows that it would be in principle possible to observationally distinguish the Halilsoy geometry from Kerr and Schwarzschild. It can be seen that in the case γ = 1 the introduction of the 'quasi-NUT' parameter does not alter the shape of the shadow, thus showing again that q must not be interpreted as a rotation parameter but is instead related to the NUT charge. On the other hand, departures from the NUT case, i.e. γ = 1 significantly alter the shape of the shadow, as expected.