1. Introduction
In 1972, nine years after Kerr obtained the metric associated to a spinning mass, Tomimatsu and Sato were able to derive a new solution to the Einstein field equations representing a stationary axisymmetric and asymptotically flat metric describing the geometry around a deformed spinning mass with deformation parameter
[
1,
2]. Such a solution was gradually generalized to the case of an arbitrary positive integer distortion parameter in the series of papers [
2,
3,
4,
5,
6]. In the present work, we are interested in the case
for which [
7] showed that the corresponding space–time is characterized by four distinct principal null directions and therefore, it is of type I according to Petrov’s classification. We recall that the Kerr metric is used instead of type D and as such, it admits a Carter constant, which plays a fundamental role in the separation process and reduction of order for the geodesic equation while the TS metric, being of general Petrov type, will not allow for Carter’s integral of motion [
8]. Another striking difference between the aforementioned manifolds is that, in contrast to the Kerr metric where an event horizon shields the ring-like curvature singularity, the TS space–time has a naked spinning singularity [
7] in the vicinity of which closed time-like curves appear [
7]. If from one side this feature is at odds with the Cosmic Censorship Conjecture (CCC) [
9] and we may be induced to dismiss the TS solution as unphysical, on the other side, we should be aware that CCC is still a conjecture and up to now, several counterexamples are known [
10,
11,
12,
13,
14,
15,
16,
17].
There has been a long debate in the literature on whether the TS metric may admit an event horizon. Both Kerr and TS have an ergoregion, but TS has a ring-like curvature singularity at the intersection of the inner surface of the ergoregion with the equatorial plane. Everybody uses prolate spheroidal coordinates while we use Boyer-Lindquist coordinates adapted to the present problem since the Kerr metric is also usually written in terms of BL-like coordinates making it easier to highlight similarities and differences between these two metrics.
One of the reasons leading us to study the TS metric is that the literature on this topic is characterized by several mistakes, misprints happening from the derivation of the TS metric and extending until Bambi where the shadow of a TS gravitational object with
is studied and instead of using the BL coordinates corresponding to that case which can be found in Yamazaki Bambi took the BL coordinates for
, i.e., the Kerr metric. This led us to go over the derivation of the TS metric in
Section 2 where we pointed out some mistakes and typos occurring in the literature. In
Section 3, we focused our attention on the study of null surfaces. We found that in addition to the ring-shaped naked singularity, there is an event horizon placed in the inner region inside the aforementioned curvature singularity. In In
Section 4, we draw our conclusions and discuss future research directions related to the TS metric and its possible extension in the presence of a positive cosmological constant.
2. Derivation of the Tomimatsu–Sato Metric with Deformation Parameter
The general form of the Lewis-Papapetrou line element describing cylindrical solutions to the vacuum Einstein field equations in cylindrical coordinates
is [
18,
19]
where the four unknown functions appearing in (
1) depend on the variables
and
z only. In the reminder of this section, we will assume that the metric coefficients satisfy the condition
. Note that if we want the line element (
1) to go over to the Minkowski metric in cylindrical coordinates as
, i.e.,
we need to require in such limit that
,
and
. If we impose that the Ricci tensor
vanishes, we find by means of Maple that the only components that do not identically vanish are
,
,
,
,
and
. More precisely, we end up with the following overdetermined coupled system of PDEs
with
Note that Equations (
3)–(
5) agree with the corresponding Equations (2.2a)–(2.2c) in [
20] where the corresponding equations arising from the Ricci components
,
and
have not been given. As already observed by [
18], the combination
gives rise to the Laplace equation
. This is a remarkable property satisfied by the function
D, which leads to a further simplification of the line element (
1). More precisely,
D is the real part of some complex-valued function
where the functions
D and
E are intertwined through the Cauchy-Riemann equations
and
. If we introduce the coordinate transformation
together with the invertibility condition
, we can rewrite the line element (
1) as
where the functions
f,
and
ℓ depend on the new variables
and
. At this point, we can exploit the freedom of choosing
D so that it satisfies the Laplace equation. The simplest choice is
, which implies that
while the Cauchy-Riemann equations lead to
modulo, an integration constant that can be chosen to be zero. Hence, we have
and as a consequence, the transformed metric (
12) coincides again with (
1) but the unknown functions
f,
and
ℓ are not independent anymore because they are connected through the relation
If we implement (
13) into (
3)–(
8), the resulting Einstein field equations read
As a consistency check, it is gratifying to observe that Equations (
14)–(
16) coincide with Equations (A)–(C) in [
19]. Moreover, (
17)–(
19) agree with (2.10a) and (2.10b) in [
20], respectively. Finally, if we consider the combinations (
17) ± (
19), we obtain Equations (2.2) and (2.2a) in [
19], namely
Note that Equations (
14)–(
16) control the functions
k,
f and
ℓ and once they are found, we can compute
by quadratures from (
18) and (
21). In order to derive the Tomimatsu–Sato metric with deformation parameter
, it is convenient to recast (
1) in the so-called Weyl-Lewis-Papapetrou form. This is achieved with the help of (
13) and by introducing the new function
As a result, the line element (
1) becomes
The corresponding Einstein field equations can be readily obtained from (
14), (
15), (
18) and (
21). They are
At this point a couple of remarks are necessary. First of all, it is not difficult to check that the Equations (
26) and (
27) are consistent, i.e.,
. Moreover,
defined through (
26) and (
27) satisfies the Equations (
17) and (
19). This can be verified by replacing
and
therein. Finally, Equation (
16) after the aforementioned substitution coincides with (
24). Hence, the new Einstein field equations are represented by (
24)–(
27). It is worth mentioning that [
19] solved the above system of equations under the additional assumption that the r.h.s. of (
25) vanishes. Such an approach leads to a solution describing a gravitational field where either the mass or the angular momentum can be different from zero. Last but not least, instead of working with the functions
f and
w, [
21,
22] introduced a new function
and derived a fourth order quasi-linear PDE for
. A less restrictive procedure than the one adopted in [
19] relies on the observation that the divergence of the vector field
defined as follows
coincides with the LHS of (
25). This is easily accomplished by showing that the equation
coincides with (
25) if
. Hence, it is possible to construct a function
such that
By means of (
29), it is not difficult to derive from (
24)–(
27) the following PDEs governing the functions
f and
u, namely
where
is the Laplace operator in cylindrical coordinates. We draw the attention of the reader to the fact that there is a typo in Equation (3.3a) in [
20] where the term
should appear squared. The Ernst equation [
23] emerges from (
30) and (
31) by introducing some complex-valued function
and realizing that the aforementioned equations are the real and imaginary parts of the complex PDE
Hence, (
34) is an equivalent form of the Einstein field Equations (
30) and (
31). In order to derive the Tomimatsu–Sato metric, it is convenient to introduce the Ansatz
with
a complex-valued function yet to be determined. The Ernst equation becomes
where the star denotes complex conjugation. Note that our notation departs from that in [
1,
2,
20] where the lowercase Greek letter
is used instead of
. The reason behind our choice is that the letter
will denote one of the prolate spheroidal coordinates. Finally, the Tomimatsu–Sato line element can be constructed by searching for solutions to (
36) in prolate spheroidal coordinates
,
with
and
such that
The corresponding inversion formulae are
We alert the reader that our notation for the prolate spheroidal coordinates
differs from that employed in [
1,
2,
20] where the lowercase Latin letters
x and
y stand for
and
, respectively. If we apply the coordinate transformation (
37) to (
36), we end up with the Ernst equation in prolate spheroidal coordinates, i.e.,
Moreover, the Equations (
29), (
32) and (
33) become
It is well-known that the vacuum stationary axisymmetric solutions to Einstein’s field equations can be generated from the Ernst Equation (
34), (
36) or equivalently from (
39). For instance, the Kerr metric emerges from the case
where
M is the mass of the gravitational source and
J its total angular momentum. A further solution to (
39) was found by [
1,
2] in the form
with
It is a straightforward exercise with Maple to verify that (
45) indeed satisfies (
39). Since
for
, in this limit, the Tomimatsu–Sato solution coincides with the class of Weyl’s metrics generated by
when
. In the case
, the function
leads to the Schwarzschild metric [
24]. Therefore,
and
q can be viewed as positive parameters measuring the deviation from spherical symmetry. If we insert (
45) into (
35) and recall that
, we can easily compute the real and imaginary parts of
with the help of Maple. More precisely, we find that
with
It is gratifying to observe that the expressions for the coefficients
A and
B represented by (
49) and (
50) agree with the corresponding ones offered in [
1,
2,
20]. Moreover,
f exhibits the desired behaviour at space-like infinity, that is
as
. Concerning the function
w, it can be derived from the Equations (
40) and (
41). In particular, we try, as in [
1,
2,
20], the Ansatz
where
M is the total mass of the gravitational source and
C is some polynomial in
and
of maximum degree seven. This requirement together with the fact that
A is a polynomial of degree eight in
is necessary in order to ensure that
as
. If we substitute (
52) into (
40) and (
41), we obtain the following system
in agreement with Equations (7) and (8) in [
3]. Setting
and substituting (
52) into (
53) leads to a linear system for the unknown coefficients
, which is easily handled by Maple. The final result is
As a consistency check, we used Maple to verify that (
55) is indeed a solution to both PDEs (
53) and (
54). It is gratifying to see that (
55) agrees with the corresponding expressions given in [
2,
3]. However, a comment is in order. [
2] also derived a system of first order PDEs for
C represented by Equations (3.3) and (3.4) therein. Even though the l.h.s.’s of (
53) and (3.3) in [
2] coincide, the same cannot be said for the r.h.s.’s. Despite the fact that the expression for
C given in [
2] agrees with the corresponding expressions presented here and in [
3], it does not satisfy Equation (3.3) in [
2] as it can be easily verified by means of Maple. The only conclusion is that such a discrepancy is due to a typo in the aforementioned equation. Finally, [
20], instead of deriving a formula for
C, gives the same result as in [
2] but with a typo in the last term. It should also be mentioned that [
20] makes an opposite choice for the signs entering in front of the l.h.s. of the equations represented by (
29) in the present work. It is not difficult to check that such a choice would alter the sign in front of the l.h.s.’s in (
53) and (
54) and send
C into
. A formula for the metric coefficient
was derived in [
2] while [
20] gives an expression for
which is linked to the aforementioned metric coefficient in [
2] by the relation
. We checked with Maple that
is indeed a solution of the system (
42) and (
43). Finally, the line element (
23) can be written in prolate spheroidal coordinates as follows
where
A,
B and
C are given by (
49), (
50) and (
55), respectively. If, instead, we use the standard definition of prolate spheroidal coordinates
with
and
where
and
, (
57) can be cast into the form
where
3. Analysis of the Metric
In this section, we improve the singularity and ergoregion analysis provided in [
7,
25]. First of all, a direct inspection of (
49), (
50), (
55) and (
56) shows that all metric coefficients are invariant under the transformation
. The fact that
B is always nonnegative implies that the sign of the metric coefficient
f is uniquely controlled by the sign of the polynomial function
A. Moreover, it can be immediately seen that the functions (
49), (
50) and (
55) have common roots at
and
. If we Taylor expand
B around
and
as follows
with
and we recall that
, we conclude that the metric coefficient
has a regular behaviour on the line segment described by
and
. The points
making
f singular have the interpretation of quasi-regular singularities [
25]. At this step, it is not obvious whether the polynomial functions
A and
B have some other common zeroes in addition to those mentioned above. Inspired by the parametric surface representation of the ergosphere of a Kerr black hole, we introduce the ansatz [
7]
with
. This allows us to cast
A into the form
First of all, we observe that a quartic polynomial similar to
has been also given in [
7]. However, there is a typo there and the factor 4 multiplying the term
should be replaced by 2. The polynomial
does not admit any real root for
and therefore, if there are some other real roots of
A distinct from those already identified above, they must come from the quartic
. This observation is important because it signalizes that
may not be everywhere positive definite for
and
. In other words, for decreasing values of
the Killing vector field
may already loose its property of being time-like before
reaches the value one and an ergoregion may arise as pointed out by [
7,
25]. Concerning the zeroes of the polynomial
, there are two sign changes and Descartes’ rule of signs predicts that the number of positive roots is either equal to two or zero. Moreover, after the transformation
, the corresponding polynomial has no sign changes and therefore,
does not have negative real roots. Finally, by applying (1.163a) and (1.163b) in [
26] the two positive real roots are
to which the following values of
correspond, determined with the help of (
66), i.e.,
where the subscript
stands for the ergoregion
Note that on the equator, i.e.,
, the inner and outer boundaries of the ergoregion are located at
In order to investigate whether or not the metric coefficient
f becomes singular, it is necessary to study the zeroes of
B. If we restrict our attention to the equatorial plane as in [
7], we end up with the problem of finding the roots of the polynomial equation
Since
exhibits only one sign change, there is only one positive real root, here denoted as
, and it must be a zero of order two for
B. Moreover, the polynomial obtained by means of the transformation
has three sign changes meaning that there can be three or one negative real roots. As it can be seen from
Table 1, it turns out that
, in agreement with [
7], where, however, no numerical/analytic evidence was given for such a result. Since
is a simple zero for
A, while
is a root of order two for
B, we can conclude that the metric coefficient
f exhibits a singularity at
on the equatorial plane. We checked numerically that
A and
B have no common real roots away from the equatorial plane (see
Appendix A for typical values of
when
). More precisely, all zeroes of
B are complex whenever
. This indicates that the metric coefficient
f can only have a singularity at the point
in addition to the singular points at
. If we recall that prolate spheroidal and cylindrical coordinates are connected by (
37), it can be easily checked that
corresponds to a ring singularity with radius
This is where the inner boundary of the ergoregion intersects the equatorial plane. Hence, such a singularity resides inside the ergoregion. The ring singularity was proved to be a curvature singularity by analysing the behaviour of the Weyl curvature invariants there (see [
25] for a detailed discussion). Ref. [
25] also showed that such a ring singularity has zero Komar mass: an unexpected result if we think that an axisymmetric space–time such as the one described by a Kerr black hole has instead a positive Komar mass [
27] despite the fact it also exhibits a ring singularity. This aspect of the TS metric and its physical implications seem not to have been addressed in the related existing literature and they are definitely worth being studied in more detail elsewhere. Finally, the ring singularity is naked because it is not hidden by a Cauchy or event horizon as it has been already discussed in [
25] (see Figure 10 therein).
It is interesting to understand how the aforementioned singularities are mapped when we switch from prolate to Boyer-Lindquist (BL) coordinates. When we reviewed the existing literature, we observed that there has been some confusion regarding the transformations from cylindrical to BL coordinates in the presence of an arbitrary deformation parameter. For instance, Ref. [
1] gives without proof the formulae
where, strangely enough, the deformation parameter completely disappeared from the expressions written in terms of the BL coordinates and the function
is the same as the one entering in the Kerr metric. Moreover, Ref. [
1] refers to [
24] for more details, but [
24] provides only the transformations from prolate spheroidal to cylindrical coordinates. Finally, the relations
and
appear in the work of [
23] where Ernst showed how the metric derived therein is equivalent to the Kerr solution if the coordinate transformation
and
is introduced. Hence, those parts in the expressions (
74) and (
75) where the
coordinates enter are correct only in the special case
. At this point, we warn the reader about the fact that [
28], when analyzing the shadow of a TS manifold with
, made use of the coordinate transformations (
74) and (
75) and therefore, the numerical results presented there should be taken with some caution. The correct transformation for arbitrary
has been given by [
4], namely
Note that the condition
requires that
Then, it is straightforward to verify that the cylindrical coordinates can be expressed as
Note that
corresponds to
and
, that is to an infinitesimally thin rod of length
. In the present work, we are interested in the case
. Hence, (
78) becomes
It is not difficult to verify that the prolate spheroidal coordinates can be expressed in terms of the BL coordinates as
and as a consistency check, it can be easily shown that under the condition
, Formulae (
81) and (82) give
in agreement with (
76). Finally, the TS line element in BL coordinates is
with
and
A,
B and
C, given by (
49), (
50) and (
55) where
and
, must be replaced according to (
83). To find the radii of the ergoregion on the equatorial plane, we can use (
69) with
and (
81) with
together with the constraint
where
is defined according to (
77) to obtain the formulae
with
given as in (
71). Here, the subscript
e stands for equator. For typical numerical values of these radii we refer to
Table 2. In
Figure 1 we represented the ergoregion for the case
where cusp singularities emerge at the top and bottom of the inner surface of the ergoregion.
At this point, some comments are in order. First of all,
remains a curvature singularity after we performed the coordinate transformation on the TS metric. Moreover, according to the discussion here below, it turns out that switching to Boyer-Lindquist coordinates has the effect of introducing an event horizon located in the inner region of the ring singularity. This behaviour is not surprising if we keep in mind the example of the Schwarzschild metric. In Kruskal-Szekeres coordinates there is no horizon, but if we rewrite the line element in spherical coordinates, an horizon appears at
. In order to proceed further, we recall that an event horizon can be characterized in terms of a null surface
, that is
[
20]. In the case of the TS metric in Boyer-Lindquist coordinates, such an equation reduces to
Let us consider the surface
, i.e.,
. Call these surfaces
and
. It can be easily verified that these surfaces indeed satisfy (
87) so that they are null surfaces. However, due to the condition
only the surface
should be taken into account. Since
when
(see
Table 1), it follows from (
86) that
. Hence, in addition to the ring-shaped naked singularity, there is an event horizon placed in the inner region inside the aforementioned curvature singularity. However, such a null surface corresponds to the segment
. Finally, we observe that the term in the first bracket of (
87) vanishes at
with
and
coinciding with
for
and
, respectively, and otherwise,
. Since
becomes purely imaginary on the interval
, we can neglect this case.