*Article* **An Invariant Characterization of the Levi-Civita Spacetimes**

**Cooper K. Watson 1,2,\*, William Julius1,2, Matthew Gorban 1,2, David D. McNutt 3, Eric W. Davis <sup>1</sup> and Gerald B. Cleaver 1,2**


**Abstract:** In the years 1917–1919 Tullio Levi-Civita published a number of papers presenting new solutions to Einstein's equations. This work, while partially translated, remains largely inaccessible to English speaking researchers. In this paper we review these solutions, and present them in a modern readable manner. We will also compute both Cartan–Karlhede and Carminati–Mclenaghan invariants such that these solutions are invariantly characterized by two distinct methods. These methods will allow for these solutions to be totally and invariantly characterized. Because of the variety of solutions considered here, this paper will also be a useful reference for those seeking to learn to apply the Cartan–Karlhede algorithm in practice.

**Keywords:** Levi-Civita metric; general relativity; curvature invariant

### **Citation:** Watson, C.K.; Julius, W.; Gorban, M.; McNutt, D.D.; Davis, E.W.; Cleaver, G.B. An Invariant Characterization of the Levi-Civita Spacetimes. *Symmetry* **2021**, *13*, 1469. https://doi.org/10.3390/sym13081469

Academic Editor: Roberto Passante

Received: 23 July 2021 Accepted: 4 August 2021 Published: 11 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

#### **1. Introduction**

In the years 1917–1919 Tullio Levi-Civita (LC) published nearly a dozen papers introducing and analyzing a variety of new solutions to Einstein's field equations (collected works in Italian available in Volume IV at [1]). Recently, several key papers have been republished in English, including two of Levi-Civita's original papers [2,3], and [4] containing an overview of several solutions not included in any of the other translations. In [2], a homogeneous Einstein–Maxwell spacetime is derived; in [3], a spacetime with a potential analogous to the logarithmic Newtonian gravitational potential is derived; and [4] discusses derivation of several degenerate vacuum spacetimes. There are additional spacetime solutions in literature not translated into English [5,6] which are similar, but distinct from the other degenerate vacuum solutions.

The age and structure of these papers has resulted in more contemporary works citing these papers in confusing or incorrect ways. Here, we clarify the structure of works on exact solutions that are of interest. First, a homogeneous solution was published in 1917 as a standalone paper [7]. Then, in the years from 1917 to 1919, a series of nine notes were published starting with [8] and ending with [9]. It is not uncommon to find the different papers in this series cited by the general heading of the entire series or by referencing to only the first article in the series.

Additionally, we will provide an invariant (local) characterization of these solutions via two different methods. First, we will utilize the Cartan–Karlhede (CK) algorithm [10] to generate an invariant coframe and the corresponding scalar quantities which uniquely characterize these spacetimes. The variety of solutions considered in this work will result in the CK algorithm running in several markedly different ways. We will present an overview of the algorithm itself, as well as a comprehensive guide that fully outlines and computes each step of the CK algorithm for the different spacetimes. Thus, this paper should serve as a useful resource for those attempting to learn to apply the CK algorithm in practice.

Second, we calculate the Carminati–McLenaghan (CM) [11] scalar invariants to construct a coframe independent classification of each solution. The set of CM invariants is advantageous as they are of the lowest possible degree and are generally the minimal independent set for any valid Petrov and Segre type spacetime. All spacetimes considered herein are of a Segre type, such that only CM invariants are needed, i.e., we do not need the extended set of invariants given in [12]. In fact, in several of the cases considered the spacetime is sufficiently special, such that only a subset of the CM invariants are needed [13]. We do note explicitly that the CM invariants will only uniquely characterize these solutions to zeroth order (in derivatives), but such invariants are useful for distinguishing LC solutions. In several cases, we will also present "I" invariants [14], as these invariants are distinct from the CM invariants and may contain information regarding algebraically special surfaces [14,15]. For completeness, we note that all spacetimes considered here are I non-degenerate, as the only case considered with constant scalar invariants is homogeneous [16]. We also note that the CK algorithm will always generate a complete classification of the spacetime, thus there is no possibility this may fail for the specific cases considered here.

We will also consider several generalizations of these solutions in cases where our methods of characterization extend directly, and in an instructive manner. In particular, we are interested in the work given by [17], which generalizes the solution in [2,7] to one which is not conformally flat, and in [18] which generalizes the work in [3,9] to a solution which is not generally static.

Throughout, the (−1, 1, 1, 1) signature convention will be used. Greek indices will be taken to run from 0 to 3, and follow the Einstein summation convention. Parentheses (brackets) around indices will denote the usual (anti)-symmetrized indices. Partial derivatives in the form *<sup>∂</sup> <sup>∂</sup><sup>x</sup>* will always be taken to be covariant, and null tetrads will always be listed in covariant form. The null vectors {*l*, *k*, *m*, *m*} will be taken to have normalization, such that *<sup>l</sup>μk<sup>μ</sup>* <sup>=</sup> <sup>−</sup>1 and *<sup>m</sup>νm<sup>ν</sup>* <sup>=</sup> 1.

#### **2. Overview of the CK Algorithm**

Here, we provide an overview of the practical CK algorithm used throughout this paper. For a review of the theoretic underpinning of the general Cartan process, see [19]. For a review of this process's application to general relativity, refer to [10,20,21]. In this algorithm we will use *q* to denote the order of differentiation, which tracks the current iteration of the algorithm. In the steps given here, we will also depart from the "standard" description of the algorithm by treating the *q* = 0 order step as a distinct "initialization" step and all steps with *q* ≥ 1 as the repeated part of the algorithm. We do this as the zeroth order step is the only step in which we will require full knowledge of the algebraic type of the tensors considered as it will usually be the step at which the parameters of the isotropy group are fixed the most.

The CK algorithm will run as follows:


or Weyl tensor into an invariant form, as it is not generally possible to find a frame which fixes both tensors into their canonical form;


#### **3. The Homogeneous Levi-Civita Solution (1917)**

In 1917, Levi-Civita presented a solution to Einstein's field equations which described a space permeated by a homogeneous, non-null Maxwell field [2,7]. This solution was later rediscovered independently (and nearly simultaneously) in [17,25], and, as such, is often called the Bertotti–Robinson metric in literature. It was [17] that presented a slight generalization (discussed below) which has a non-vanishing cosmological constant. It was shown in [26] that this metric was generally singularity free. This spacetime can also be shown to be a limiting case of the more general Petrov solution [27].

Later considerations of more general Maxwell spacetimes have also revealed several interesting properties regarding this solution. First, this is the only homogeneous nonnull Maxwell solution [28]. Generalizations where the Maxwell field does not share the homogeneous symmetry also give interesting solutions (not discussed here) which are algebraically more general [29,30].

#### *3.1. Forms of the Metric and Nature of the Coordinates*

In present literature, there are at least six equivalent forms of the homogeneous 1917 Levi-Civita line element. Here, we will introduce these different forms and discuss the relations between them when possible and mention any relevant coordinate artifacts that might be present. We will also present several new forms of this metric and discuss several cases in which no coordinate transforms exist in the literature, nor can it be derived without the use of complex transforms.

The original form of the line element [2] was given in cylindrical coordinates as:

$$ds^2 = -\left(c\_1 e^{z/a} + c\_2 e^{-z/a}\right)^2 dt^2 + d\rho^2 + a^2 \sin(\rho/a)^2 d\phi^2 + dz^2,\tag{1}$$

where *c*<sup>1</sup> and *c*<sup>2</sup> are integration constants originally derived by Levi-Civita and *a* is a constant usually associated with an electric or magnetic field strength [2,7]. If *c*<sup>1</sup> and *c*<sup>2</sup> are both non-zero and of opposite sign, this metric will be singular for *z* = *<sup>a</sup>* <sup>2</sup> ln | *c*2 *c*1 |. Additionally, this metric is singular for *<sup>ρ</sup>* <sup>=</sup> *<sup>n</sup>πa*, where *<sup>n</sup>* <sup>∈</sup> <sup>Z</sup>.

One may define an angular parameter *θ* = *ρ*/*a*, such that (1) can be rewritten in the form given in [31] as:

$$ds^2 = -\left(c\_1 e^{z/a} + c\_2 e^{-z/a}\right)^2 dt^2 + dz^2 + a^2 d\Omega^2,\tag{2}$$

where *<sup>d</sup>*Ω<sup>2</sup> is the line element for the 2-sphere. We explicitly note that *<sup>θ</sup>* <sup>∈</sup> <sup>R</sup>. This metric can also be be put into a secondary set of cylindrical coordinates [31] by taking *ρ* /*a* = sin *θ* giving:

$$ds^2 = -\left(c\_1 e^{z/a} + c\_2 e^{-z/a}\right)^2 dt^2 + \frac{d\rho'^2}{\left(1 - \left(\rho'/a\right)^2\right)} + \rho'^2 d\phi^2 + dz^2,\tag{3}$$

We note here that *ρ* ∈ [0, *a*]. In these coordinates the countably infinite singularities in the original radial coordinate have been reduced to only *ρ* = 0 and *ρ* = *a*.

We may once again rewrite (1), this time in the associated Cartesian coordinates, previously exploited in [32], as:

$$ds^2 = -\left(c\_1 \varepsilon^{z/a} + c\_2 \varepsilon^{-z/a}\right)^2 dt^2 + dx'^2 + dy'^2 + dz^2 + \frac{\left(\mathbf{x'} d\mathbf{x'} + \mathbf{y'} dy'\right)^2}{a^2 - \left(\mathbf{x'}^2 + \mathbf{y'}^2\right)'},\tag{4}$$

where we explicitly note that the primed coordinates *x* and *y* are restricted to take values subject to *<sup>x</sup>* <sup>2</sup> <sup>+</sup> *<sup>y</sup>* <sup>2</sup> <sup>≤</sup> *<sup>a</sup>*2.

It is also possible to eliminate *c*<sup>1</sup> and *c*<sup>2</sup> via a transform on both *t* and *z*. By scaling *t* and translating *z*, the original metric can be rewritten into one of two equivalent forms:

$$ds^2 = -\sinh(z)^2dt^2 + a^2dz^2 + a^2d\Omega^2\tag{5}$$

or

$$ds^2 = -\cosh(z)^2dt^2 + a^2dz^2 + a^2d\Omega^2,\tag{6}$$

where, for real coordinate transforms, one can get either a *sinh* or *cosh* solution depending on the relative signs of *c*<sup>1</sup> and *c*2. Interestingly, it appears that the only method of connecting these two *equivalent* solutions is via complex transformations on both *t* and *z*, although it is not presently understood why this method works. We do note that this property of certain complex transformations, resulting in the same solution in different forms, is remarkably similar to the application of the Newman–Janis trick when applied to Minkowski space [33].

Here we also present a new form of this metric which is not related to previous solutions by any known transform, real or complex:

$$ds^2 = -e^{2z}dt^2 + a^2dz^2 + a^2d\Omega^2. \tag{7}$$

Using this form and making the coordinate transform *r* = *e*−*<sup>z</sup>* (along with a rescaling in *t*) we get the form seen in [20], given as:

$$ds^2 = \frac{a^2}{r^2} \left(dr^2 - dt^2\right) + a^2 d\Omega^2. \tag{8}$$

This reference also presents the form

$$ds^2 = -\left(1 + \frac{z^2}{a^2}\right)dt^2 + \left(1 - \frac{y^2}{a^2}\right)dx^2 + \left(1 - \frac{y^2}{a^2}\right)^{-1}dy^2 + \left(1 + \frac{z^2}{a^2}\right)^{-1}dz^2,\tag{9}$$

which we see is a special case of the more general metric considered in (17).

Throughout the rest of this paper we will work with the metric and coordinates given by (8). A convenient choice of covariant null tetrad is given by:

$$\begin{split} l &= \frac{a}{\sqrt{2}r} \left( \frac{\partial}{\partial r} - \frac{\partial}{\partial t} \right), \quad k = -\frac{a}{\sqrt{2}r} \left( \frac{\partial}{\partial r} + \frac{\partial}{\partial t} \right), \\ m &= \frac{a}{\sqrt{2}} \left( \frac{\partial}{\partial \theta} + i \sin \left( \theta \right) \frac{\partial}{\partial \phi} \right), \quad \overline{m} = \frac{a}{\sqrt{2}} \left( \frac{\partial}{\partial \theta} - i \sin \left( \theta \right) \frac{\partial}{\partial \phi} \right). \end{split} \tag{10}$$

#### *3.2. Curvature Invariants and CK Classification*

Using the null tetrad defined in (10), one may calculate the only non-vanishing curvature component

$$
\Phi\_{00} = \frac{1}{2a^2}.\tag{11}
$$

This constant solution is therefore both conformally flat and Ricci flat. Since this curvature component is constant and is an invariantly defined frame, (11) is the only non-vanishing CK invariant [24].

Additionally, this result fixes all CM invariants to be either zero or of the form *ca*2*<sup>n</sup>* for *<sup>a</sup>*, *<sup>c</sup>* <sup>∈</sup> <sup>R</sup> and *<sup>n</sup>* <sup>∈</sup> <sup>Z</sup>. For example, the only two non-vanishing CM invariants are

*r*<sup>1</sup> = 2 <sup>√</sup>*r*<sup>3</sup> <sup>=</sup> <sup>1</sup> *<sup>a</sup>*<sup>4</sup> , (12)

where only *r*<sup>1</sup> is independent since the spacetime is of warped product type *B*<sup>2</sup> [13].

#### *3.3. Regarding Electromagnetic "Wormholes"*

It has been suggested that this solution may in some sense constitute a wormhole supported by electromagnetic stress [34,35]. It was shown in [31,36,37] that, despite there being a coordinate singularity at *ρ* = *a* in (3), this does not correspond to a wormhole throat, as an appropriate choice of coordinate transform can be made, such that the spatial part of the metric becomes that of a hypercylinder. Here, we will highlight two different methods of characterizing this surface, which provide a secondary method of determining that this solution is not a wormhole.

Working with the coordinates and metric given by (3), we choose the following null tetrad to analyze the surface *ρ* = *a*, where we have dropped the prime out of convenience,

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( - \left| c\_1 e^{z/a} + c\_2 e^{-z/a} \right| \frac{\partial}{\partial t} + \frac{a}{\sqrt{a^2 - \rho^2}} \frac{\partial}{\partial \rho} \right), \\ k &= \frac{1}{\sqrt{2}} \left( - \left| c\_1 e^{z/a} + c\_2 e^{-z/a} \right| \frac{\partial}{\partial t} - \frac{a}{\sqrt{a^2 - \rho^2}} \frac{\partial}{\partial \rho} \right), \\ m &= \frac{1}{\sqrt{2}} \left( \rho \frac{\partial}{\partial \phi} + i \frac{\partial}{\partial z} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \rho \frac{\partial}{\partial \phi} - i \frac{\partial}{\partial z} \right). \end{split} \tag{13}$$

The null expansions along *l* and *k* are, respectively,

$$
\theta\_{(l)} = \dot{q}^{\mu\nu} \nabla\_{\mu} l\_{\nu} \quad \text{and} \quad \theta\_{(k)} = \dot{q}^{\mu\nu} \nabla\_{\mu} k\_{\nu} \tag{14}
$$

where *qμν* = *gμν* + 2*l*(*μnν*) is a local, induced two-metric and ∇*<sup>μ</sup>* is the standard covariant derivative. In terms of the given null frame, these are explicitly:

$$
\theta\_{(l)} = -\theta\_{(n)} = -\frac{\sqrt{a^2 - \rho^2}}{\sqrt{2}a\rho}.\tag{15}
$$

It can be seen that the surface *ρ* = *a* is indeed a surface on which the expansion of these null directions vanish. For this surface to correspond to a wormhole throat, the derivative of the expansion projected along the null direction must be positive on the surface [38]. Computing these terms at the surface we can show that

$$k^{\mu}\nabla\_{\mu}\theta\_{(l)} = k^{\mu}\nabla\_{\mu}\theta\_{(k)} = -\frac{1}{2a^2}.\tag{16}$$

Thus, this surface cannot correspond to a wormhole throat. It is in fact only a maximal surface relating to the given range of coordinates [39,40]. This surface is also not detectable by any scalar invariants (since all invariants are constant) and, thus, is not a geometric surface, as would be expected for a wormhole throat [15].

#### *3.4. Bertotti Generalization*

In [17], a generalization of the 1917 solution is used, with the line element taking the form of

$$ds^2 = -\left(1 + \frac{y^2}{r\_+^2}\right)dt^2 + \left(1 - \frac{z^2}{r\_-^2}\right)dx^2 + \left(1 + \frac{y^2}{r\_+^2}\right)^{-1}dy^2 + \left(1 - \frac{z^2}{r\_-^2}\right)^{-1}dz^2\,,\tag{17}$$

where *x*, *y*, *z* are Cartesian coordinates. If *r*<sup>+</sup> = *r*−, this solution is no longer conformally flat or Ricci flat. In the case where *r*<sup>+</sup> = *r*<sup>−</sup> = *a*, the spacetime becomes conformally flat and reduces to (9).

Using the null tetrad:

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\sqrt{1 + \frac{\mathbf{x}^{2}}{r\_{+}^{2}}} \frac{\partial}{\partial t} + \frac{1}{\sqrt{1 + \frac{\mathbf{x}^{2}}{r\_{+}^{2}}}} \frac{\partial}{\partial x} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\sqrt{1 + \frac{\mathbf{x}^{2}}{r\_{+}^{2}}} \frac{\partial}{\partial t} - \frac{1}{\sqrt{1 + \frac{\mathbf{x}^{2}}{r\_{+}^{2}}}} \frac{\partial}{\partial x} \right), \\\ m &= \frac{1}{\sqrt{2}} \left( \sqrt{1 - \frac{\mathbf{x}^{2}}{r\_{-}^{2}}} \frac{\partial}{\partial y} + \frac{i}{\sqrt{1 - \frac{\mathbf{x}^{2}}{r\_{-}^{2}}}} \frac{\partial}{\partial z} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \sqrt{1 - \frac{\mathbf{x}^{2}}{r\_{-}^{2}}} \frac{\partial}{\partial y} - \frac{i}{\sqrt{1 - \frac{\mathbf{x}^{2}}{r\_{-}^{2}}}} \frac{\partial}{\partial z} \right), \end{split} \tag{18}$$

the only non-vanishing curvature components are found to be

$$\begin{aligned} R &= 2\left(\frac{1}{r\_{-}^{2}} - \frac{1}{r\_{+}^{2}}\right), \\ \Phi\_{11} &= \frac{1}{4}\left(\frac{1}{r\_{+}^{2}} + \frac{1}{r\_{-}^{2}}\right), \quad \Psi\_{2} = \frac{1}{6}\left(\frac{1}{r\_{+}^{2}} - \frac{1}{r\_{-}^{2}}\right), \end{aligned} \tag{19}$$

which are again all constants. Since both terms in (19) are constant, and this is already an invariantly defined frame to zeroth order, these are the only three non-vanishing Cartan invariants.

Since this spacetime is a warped product spacetime of type *B*2, the four invariants given below are the complete set [13]. These are:

$$R = 2\left(\frac{1}{r\_{-}^{2}} - \frac{1}{r\_{+}^{2}}\right), \quad r\_{1} = \frac{1}{4}\left(\frac{1}{r\_{-}^{2}} + \frac{1}{r\_{+}^{2}}\right)^{2}, \quad r\_{2} = 0, \quad w\_{2} = \frac{1}{36}\left(\frac{1}{r\_{-}^{2}} - \frac{1}{r\_{+}^{2}}\right)^{3}.\tag{20}$$

#### **4. The Cylindrical Levi-Civita Solution (1919)**

In [3,9], a solution to Einstein's field equations is presented that serves as the analog to the Newtonian logarithmic potential. In [41], the CK invariants have already been computed, and in [42] these invariants were used to show that this solution can be found as a limiting subcase of the *γ* (or "Zipoy-Voorhees") [43,44] solution. Here, for completeness, we will independently compute the CK invariants.

This spacetime solution has also been of interest as an exterior vacuum solution to various physical sources, see [41,45]. There has also been interest in global (topological) properties of this solution as certain parameters (which do no affect local properties) are related to cosmic strings, see [45–47] for discussion of these properties as they pertain to this solution, and see [48] for a general review of cosmic strings.

#### *4.1. Forms of the Metric and Nature of the Coordinates*

The line element for this spacetime is given by

$$ds^2 = -e^{2\upsilon}dt^2 + e^{-2\upsilon} \left( e^{2\lambda} \left( d\rho^2 + dz^2 \right) + r^2 d\phi^2 \right),\tag{21}$$

in standard cylindrical coordinates, where *ν*(*ρ*, *z*) is a solution to

$$\frac{1}{\rho} \frac{\partial}{\partial \rho} \left( \rho \frac{\partial \nu}{\partial \rho} \right) + \frac{\partial^2 \nu}{\partial z^2} = 0,\tag{22}$$

and *λ* is defined, up to a constant of integration, by the differential relation:

$$d\lambda = \rho \left(\frac{\partial \nu}{\partial \rho}^2 - \frac{\partial \nu}{\partial z}^2\right) dr + 2\rho \nu\_1 \nu\_2 dz. \tag{23}$$

Once solved, particular solutions for *ν* and *λ* are given by

$$\varepsilon^{-\nu} = \left(\frac{\rho\_0}{\rho}\right)^h, \quad \varepsilon^{\lambda-\nu} = \left(\frac{\rho}{\rho\_0}\right)^{h^2-h},\tag{24}$$

where *h* is an arbitrary real constant. Note that, with the form of line element given by (21), it is possible to choose coordinates, such that *ρ*<sup>0</sup> may be eliminated from the metric, although this will result in *φ* not being parameterized from (0, 2*π*). This will correspond to a global angular defect discussed in [47,49]. Since the metric is independent of the angular coordinate, such a reparameterization will not necessarily affect the local classification given by either CM or CK invariants. Explicitly, we will take this metric to be

$$ds^2 = -\rho^{2h}dt^2 + \rho^{2\left(h^2 - h\right)} \left(d\rho^2 + dz^2\right) + \rho^{2\left(1 - h\right)} dx\_{3\prime}^2 \tag{25}$$

and use it for all calculations going forward. We note, once again for clarity, that *x*<sup>3</sup> ∈ [0, *a*] and *a* is determined by the specific value of *r*<sup>0</sup> taken above.

Going forward we will take a null frame given by:

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\rho^{h} \frac{\partial}{\partial t} + \rho^{h^{2}-h} \frac{\partial}{\partial \rho} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\rho^{h} \frac{\partial}{\partial t} - \rho^{h^{2}-h} \frac{\partial}{\partial \rho} \right), \\ m &= \frac{1}{\sqrt{2}} \left( \rho^{h^{2}-h} \frac{\partial}{\partial z} + i\rho^{h-1} \frac{\partial}{\partial \mathbf{x}\_{3}} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \rho^{h^{2}-h} \frac{\partial}{\partial z} - i\rho^{h-1} \frac{\partial}{\partial \mathbf{x}\_{3}} \right). \end{split} \tag{26}$$

### *4.2. Curvature Invariants and CK Classification*

In the above frame, the only non-vanishing components of curvature are

$$\Psi\_0 = \Psi\_4 = \frac{1}{2}h(h^2 - 1)\rho^{-2\left(h^2 - h + 1\right)} \quad \text{and} \quad \Psi\_2 = -\frac{1}{2}h(h - 1)^2 \rho^{-2\left(h^2 - h + 1\right)}.\tag{27}$$

This frame is invariantly defined and thus these are the zeroth order CK invariants. Additionally, this frame fixes out all isotropy. Taking the covariant derivatives, the only non-vanishing terms are

$$\begin{split} \mathbb{C}\_{\text{kn}\mathbb{N};\eta;k} &= -\mathbb{C}\_{\text{ll}\mathbb{M}\mathbb{N};l} = \frac{(h^{5}-h)\rho^{-3}(1+h^{2}-h)}{\sqrt{2}},\\ \mathbb{C}\_{\text{kn}\mathbb{N};l} = \mathbb{C}\_{\text{kn}\mathbb{M};m} = -\mathbb{C}\_{\text{kl}\mathbb{M};\mathbb{N}} &= -\mathbb{C}\_{\text{ll}\mathbb{M}\mathbb{M};\mathbb{N}} = -\frac{(h^{2}+h)(h-1)^{3}\rho^{-3}(1+h^{2}-h)}{\sqrt{2}},\\ \mathbb{C}\_{\text{kn}\mathbb{M};\mathbb{N}} = \mathbb{C}\_{\text{kn}\mathbb{M};k} = -\mathbb{C}\_{\text{kn}\mathbb{M};l} = -\mathbb{C}\_{\text{kl}\mathbb{M};m} = -\frac{(h-1)^{2}\left(h^{3}-h^{2}+h\right)\rho^{-3}(1+h^{2}-h)}{\sqrt{2}}. \end{split} \tag{28}$$

none of which are functionally independent from the zeroth order invariants. Thus, the algorithm terminates here.

We do note that for certain special cases this spacetime is of type D *<sup>h</sup>* <sup>=</sup> <sup>−</sup>1, <sup>1</sup> <sup>2</sup> , 2 [20,42]. In these cases, the algorithm must be approached separately as the frame given is only invariantly defined when *h* = −1. Additionally, since the Weyl tensor is now more algebraically special, there will be boost and spin isotropy remaining at zeroth order. This means that, in general, one will need to compute higher order derivatives to fully classify these special cases. Here, we will work each of these three cases explicitly.

When *h* = −1, the frame given by (25) becomes

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\rho^{-1} \frac{\partial}{\partial t} + \rho^{2} \frac{\partial}{\partial \rho} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\rho^{-1} \frac{\partial}{\partial t} - \rho^{2} \frac{\partial}{\partial \rho} \right), \\ m &= \frac{1}{\sqrt{2}} \left( \rho^{2} \frac{\partial}{\partial z} + i\rho^{2} \frac{\partial}{\partial x\_{3}} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \rho^{2} \frac{\partial}{\partial z} - i\rho^{2} \frac{\partial}{\partial x\_{3}} \right). \end{split} \tag{29}$$

In this frame, the only non-vanishing curvature component is

$$
\Psi\_2 = 2\rho^{-6},
\tag{30}
$$

where we have one functionally independent term, and remaining isotropy

$$
\begin{pmatrix} \alpha & 0 \\ 0 & \alpha^{-1} \end{pmatrix}, \text{with } \alpha \in \mathbb{C} \text{ and } \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix}. \tag{31}
$$

At first order, the non-vanishing derivatives are

$$\mathcal{C}\_{kmkl;\mathfrak{M}} = \mathcal{C}\_{km\mathfrak{M};k} = -\mathcal{C}\_{km\mathfrak{M};l} = -\mathcal{C}\_{kl\mathfrak{M};m} = 6\sqrt{2}\rho^{-9},\tag{32}$$

which are not functionally independent of the zeroth order components. At first order, the isotropy group is reduced to just

$$
\begin{pmatrix} e^{i\theta} & 0 \\ 0 & e^{-i\theta} \end{pmatrix} \text{, with } \theta \in \mathbb{R}. \tag{33}
$$

The non-vanishing second derivative components are

$$\begin{aligned} \mathbf{C}\_{kmlm;\overline{m}\overline{m}} &= \mathbf{C}\_{kmkl;km} = \mathbf{C}\_{kmkl;\overline{m}k} = \mathbf{C}\_{km\overline{m}l;kl} = \mathbf{C}\_{km\overline{m}l;lm} = \mathbf{C}\_{kl\overline{m}l;lm} = \mathbf{C}\_{kl\overline{m}l;ml} \\ &= \mathbf{C}\_{\mathfrak{m}l\mathfrak{M};mm} = -\mathbf{C}\_{kl\mathfrak{M};\mathfrak{M}m} = -\mathbf{C}\_{km\mathfrak{M};l\mathfrak{M}} = -\mathbf{C}\_{km\mathfrak{M};\mathfrak{M}m} = -\mathbf{C}\_{kl\mathfrak{M};km} \\ &= -\frac{4}{5}\mathbf{C}\_{kmkl;\overline{m}l} = -\frac{4}{5}\mathbf{C}\_{km\overline{m}l;kl} = -\frac{4}{5}\mathbf{C}\_{km\overline{m}l;lk} = -\frac{4}{5}\mathbf{C}\_{kl\overline{m}l;mk} = \frac{48}{\rho^{12}} \end{aligned} \tag{34}$$

which do not reduce the above isotropy. Thus, the algorithm stops at second order.

In the following two cases, the functional independence follows exactly as above. At zeroth order the isotropy is identical, but in both cases, at first order, it is reduced to

$$
\begin{pmatrix} \beta & 0 \\ 0 & \beta^{-1} \end{pmatrix}, \text{with } \beta \in \mathbb{R}. \tag{35}
$$

For *h* = <sup>1</sup> <sup>2</sup> , using the frame:

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\sqrt{\rho} \frac{\partial}{\partial t} + \sqrt{\rho} \frac{\partial}{\partial \mathbf{x}\_3} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\sqrt{\rho} \frac{\partial}{\partial t} - \sqrt{\rho} \frac{\partial}{\partial \mathbf{x}\_3} \right), \\ m &= \frac{1}{\sqrt{2}} \left( \rho^{-1/4} \frac{\partial}{\partial \rho} + i\rho^{-1/4} \frac{\partial}{\partial z} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \rho^{-1/4} \frac{\partial}{\partial \rho} - i\rho^{-1/4} \frac{\partial}{\partial z} \right), \end{split} \tag{36}$$

we have the CK invariants:

$$\begin{aligned} \Psi\_{2} &= \frac{1}{8} \rho^{-3/2}, \\ \mathbf{C}\_{kmkl;l} = \mathbf{C}\_{km\mathfrak{M};m} = \mathbf{C}\_{km\mathfrak{M};\mathfrak{M}} = \mathbf{C}\_{kl\mathfrak{M};k} = -\frac{3}{16\sqrt{2}} \rho^{-9/4}, \\ \mathbf{C}\_{km\mathfrak{M};l\mathfrak{l}} = \mathbf{C}\_{km\mathfrak{l};l\mathfrak{m}} = \mathbf{C}\_{km\mathfrak{l};m\mathfrak{l}} = \mathbf{C}\_{km\mathfrak{l};\mathfrak{M}} = \mathbf{C}\_{km\mathfrak{M};\mathfrak{l}k} = \mathbf{C}\_{km\mathfrak{M};\mathfrak{l}k} \\ = \mathbf{C}\_{km\mathfrak{M};\mathfrak{m}\mathfrak{l}} = \mathbf{C}\_{km\mathfrak{M};\mathfrak{l}\mathfrak{M}} = \mathbf{C}\_{kl\mathfrak{M};\mathfrak{l}\mathfrak{M}} = \mathbf{C}\_{kl\mathfrak{M};\mathfrak{l}\mathfrak{M}} = \mathbf{C}\_{kl\mathfrak{M};\mathfrak{l}\mathfrak{M}} = \mathbf{C}\_{\mathfrak{M};\mathfrak{l}\mathfrak{M};\mathfrak{m}} = \mathbf{C}\_{\mathfrak{M};\mathfrak{l}\mathfrak{M};\mathfrak{m}} \\ = \frac{4}{5} \mathbf{C}\_{km\mathfrak{M};\mathfrak{l}\mathfrak{M}} = \frac{4}{5} \mathbf{C}\_{km\mathfrak{M};m\mathfrak{l}} = \frac{4}{5} \mathbf{C}\_{km\mathfrak{M};\mathfrak{l}\mathfrak{m}} = \frac{4}{5} \mathbf{C}\_{kl\mathfrak{M};\mathfrak{m}} = \frac{3}{16} \rho^{-3}. \end{aligned} \tag{37}$$

For *h* = 2, using the frame:

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\rho^2 \frac{\partial}{\partial t} + \rho^2 \frac{\partial}{\partial z} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\rho^2 \frac{\partial}{\partial t} - \rho^2 \frac{\partial}{\partial z} \right), \\ m &= \frac{1}{\sqrt{2}} \left( \rho^2 \frac{\partial}{\partial \rho} + \frac{i}{\rho} \frac{\partial}{\partial \mathbf{x}\_3} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \rho^2 \frac{\partial}{\partial \rho} - \frac{i}{\rho} \frac{\partial}{\partial \mathbf{x}\_3} \right), \end{split} \tag{38}$$

we have the CK invariants:

$$\begin{aligned} \Psi\_2 &= 2\rho^{-6}, \\ \mathbb{C}\_{km\mathbb{M};l} = \mathbb{C}\_{km\overline{m}l;m} = \mathbb{C}\_{km\mathbb{M};\overline{m}} = \mathbb{C}\_{kl\overline{m}l;k} = -6\sqrt{2}\rho^{-9}, \\ \mathbb{C}\_{km\mathbb{M};l} = \mathbb{C}\_{km\mathbb{M};m} = \mathbb{C}\_{km\mathbb{M};\mathbb{M}} = \mathbb{C}\_{km\mathbb{M};\mathbb{M}} = \mathbb{C}\_{km\mathbb{M};l\mathbb{k}} = \mathbb{C}\_{km\mathbb{M};\mathbb{l}k} \\ = \mathbb{C}\_{km\mathbb{M};\mathbb{m}m} = \mathbb{C}\_{km\mathbb{M};\mathbb{l}\mathbb{M}} = \mathbb{C}\_{kl\mathbb{M};\mathbb{l}\mathbb{M}} = \mathbb{C}\_{kl\mathbb{M};\mathbb{l}\mathbb{M}} = \mathbb{C}\_{kl\mathbb{M};\mathbb{l}k} = \mathbb{C}\_{\mathbb{M};l\mathbb{M};\mathbb{k}k} \\ = \frac{4}{5}\mathbb{C}\_{km\mathbb{M};\mathbb{l}\mathbb{M}} = \frac{4}{5}\mathbb{C}\_{km\mathbb{M};\mathbb{m}m} = \frac{4}{5}\mathbb{C}\_{km\mathbb{M};\mathbb{l}\mathbb{M}} = \frac{4}{5}\mathbb{C}\_{kl\mathbb{M};\mathbb{k}l} = 48\rho^{-12}. \end{aligned} \tag{39}$$

Here, we can see that the *h* = 1/2 and *h* = 2 solutions are in fact identical as their CK invariants can be found to be compatible.

In the case where *h* = 0, 1, the spacetime is (locally) flat, and, thus, all invariants will vanish, although globally there are topological properties mentioned above not captured by this approach.

The only two non-vanishing CM invariants are

$$w\_1 = 2\left(h^2 - h\right)^2 \left(h^2 - h + 1\right) \rho^{-4\left(h^2 - h + 1\right)}, \quad w\_2 = -3\left(h^2 - h\right)^4 \rho^{-6\left(h^2 - h + 1\right)}, \tag{40}$$

which reduces as expected for the special values of *h*.

#### *4.3. Kasner Generalization*

In [18], a generalization of the 1919 solution was given as

$$ds^2 = -r^{2D}dt^2 + t^{2A}dr^2 + r^{2E}t^{2B}dz^2 + \alpha^2r^{2F}t^{2C}d\phi^2,\tag{41}$$

which is in general stationary (rather than static). The constants in this solution are related via the Kasner constraints

$$A + B + C = A^2 + B^2 + C^2 = 1, \quad D + E + F = D^2 + E^2 + F^2 = 1,\tag{42}$$

and take specific values given by:

$$A = \frac{2s + H}{S + H}, \quad B = \frac{(2s - 1)(2s + \epsilon(2s - 1))}{S + H}, \quad C = \frac{(1 - 2s)(1 + \epsilon(2s - 1))}{S + H},$$

$$D = \frac{2s}{S}, \quad E = \frac{2s(2s - 1)}{S}, \quad F = \frac{1 - 2s}{S},$$

where *s* is a real constant, = ±1, and

$$H = \varepsilon \left( 4s^2 - 1 \right) + \left( 1 - 2s \right)^2, \quad S = 4s^2 - 2s + 1. \tag{44}$$

Computing the curvature with respect to the null tetrad

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -r^D \frac{\partial}{\partial t} + t^A \frac{\partial}{\partial r} \right), \quad k = \frac{1}{\sqrt{2}} \left( -r^D \frac{\partial}{\partial t} - t^A \frac{\partial}{\partial r} \right), \\ m &= \frac{1}{\sqrt{2}} \left( r^F t^B \frac{\partial}{\partial z} + ia^2 r^F t^C \frac{\partial}{\partial \phi} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( r^F t^B \frac{\partial}{\partial z} - ia^2 r^F t^C \frac{\partial}{\partial \phi} \right), \end{split} \tag{45}$$

we have the non-vanishing Ricci components:

$$\begin{aligned} \Phi\_{00} &= -\frac{1}{2} (B(D-E) + \mathcal{C}(D-F) + A(E+F))t^{-2+E+F}r^{-2+B+C}, \\ \Phi\_{22} &= \frac{1}{2} (B+\mathcal{C}+E-3BE-2\mathcal{C}E+F-2BF-3CF)r^{-2+E+F}t^{-2+B+C}, \end{aligned} \tag{46}$$

and

$$\begin{split} \Psi\_{0} &= \frac{1}{2r^{2}t^{2}} \Big( \mathbb{C}(-2+B+2C)r^{2(E+F)} \\ &+ (\mathbb{C}-E+B(3E-1)+F-3\mathbb{C}F)r^{E+F}t^{B+C} + F(-2+E+2F)t^{2(B+C)} \Big), \\ \Psi\_{2} &= \frac{1}{2} \Big( -\frac{BCr^{-2D}}{t^{2}} + \frac{EFt^{-2A}}{r^{2}} \Big), \\ \Psi\_{4} &= \frac{1}{2r^{2}t^{2}} \Big( \mathbb{C}(-2+B+2C)r^{2(E+F)} \\ &+ (B+E-3BE-F+\mathbb{C}(3F-1))r^{E+F}t^{B+C} + F(-2+E+2F)t^{2(B+C)} \Big). \end{split} \tag{47}$$

This frame is not canonical, as this spacetime is Petrov type I, and Ψ<sup>0</sup> = Ψ4, but the boost

$$
\begin{pmatrix}
(\Psi\_0/\Psi\_4)^{1/8} & 0\\0 & (\Psi\_4/\Psi\_0)^{1/8}
\end{pmatrix},
\tag{48}
$$

will bring it to the canonical form

$$\begin{aligned} \Psi\_0' &= \Psi\_4' = \sqrt{\Psi\_0 \Psi\_4} \\ \Psi\_2' &= \Psi\_2. \end{aligned} \tag{49}$$

For compactness we will write curvature components of the new canonical frame in terms of the old frame, where the new frame will be primed and the old frame unprimed. At zeroth order we have no remaining isotropy and have only two functionally independent terms. The first order derivatives, shown below, cannot reduce the isotropy further and contain no new functionally independent terms, and so the algorithm terminates.

$$\begin{split} \mathbb{C}\_{\text{knw};k} &= \frac{1}{2\sqrt{2}r^{3}t^{3}} \Big( 2\mathbb{C} \Big( 4 - 3\mathbb{B} - 5\mathbb{C} + 2\mathbb{B}\mathbb{C} + \mathbb{C}^{2} \Big) r^{3(E+F)} + (4(E-F) \\ + \mathbb{B}(1 - 6E + \mathbb{C}(7 - 4E - 13F) + 3F) + \mathbb{C}(-15 + 14E + \mathbb{C}(14 - 17E - 17F)) \\ + 23F)) r^{2(E+F)} t^{B+\mathbb{C}} &+ (-E + 15F - 7F(E + 2F) + \mathbb{C}(4 + E(-3 + 13F) \\ + F(-23 + 17F)) + \mathbb{B}(-4 + E(6 + 4F) + F(-14 + 17F))) r^{E+F} t^{2(B+\mathbb{C})} \\ &- 2F \Big( 4 - 3E - 5F + 2EF + F^{2} \Big) t^{3(B+\mathbb{C})} \Big) \Big( \frac{\Psi\_{4}}{\Psi\_{0}} \Big)^{3/4}, \end{split}$$

$$\begin{split} \mathbb{C}\_{\text{knuku};l} &= -\frac{1}{2\sqrt{2}r^3t^3} \Big( 2\mathbb{C}((-1+\mathbb{C})\mathbb{C}+B(-1+2\mathbb{C}))r^{3(E+F)} + (\mathbb{C}(-1+2E) \\ &+ \mathbb{C}(2-3E-3F) + F) + B(-1+\mathbb{C}+2E+F-3\mathbb{C}F))r^{2(E+F)}t^{B+C} \\ &+ (E(-1+2B+\mathbb{C}+F-3\mathbb{C}F)+F(-1+2B+\mathbb{C}+2F-3(B+\mathbb{C})F))r^{E+F}t^{2(B+C)} \\ &+ 2F((-1+F)F+E(-1+2F))t^{3(B+C)} \Big) \left(\frac{\mathbb{V}\_4}{\mathbb{V}\_0}\right)^{1/4}, \end{split}$$

$$\begin{split} \mathbb{C}\_{\text{kml};m} &= \frac{1}{4\sqrt{2}r^{3}t^{3}} \Big( -4\mathbb{C}((-1+\mathbb{C})\mathbb{C} + B(-1+2\mathbb{C}))r^{3(E+F)} \\ &+ (-B(-1+\mathbb{C}+2E+4\mathbb{C}E+F-5\mathbb{C}F) + \mathbb{C}(1+F+\mathbb{C}(-2+E+F))) \\ &\times r^{2(E+F)}t^{B+\mathbb{C}} + (F(1+\mathbb{C}+(-2+B+\mathbb{C})F)-E(-1+\mathbb{C}+F-5\mathbb{C}F) \\ &+ B(2+4F)))r^{E+F}t^{2(B+\mathbb{C})} - 4F((-1+F)F+E(-1+2F))t^{3(B+\mathbb{C})} \Big) \Big( \frac{\mathbb{Y}\_{4}}{\mathbb{Y}\_{0}} \Big)^{1/4} \Big) \end{split}$$

$$\begin{split} \mathbb{C}\_{\text{knl};\mathbb{H}} &= \frac{1}{4\sqrt{2}}r^{-3+E+F}t^{-3+B+C}((B+\mathbb{C}-2BE+(-2+\mathbb{C})\mathbb{C}E) \\ &- BF-\mathbb{C}(1+\mathbb{C})F+BC(-3+2E+F))r^{E+F}+4B\mathbb{C}r^{2(E+F)}t^{-B-\mathbb{C}} \\ &+ (E-\mathbb{C}E+2BE(-1+F)+F+(-3+\mathbb{C})EF+B(-2+F)F-\mathbb{C}F(1+F))) \\ &\times t^{B+\mathbb{C}}+4EFr^{-E-F}t^{2(B+\mathbb{C})}\Big{)}\left(\frac{\mathbb{V}\_{4}}{\mathbb{V}\_{0}}\right)^{1/4} \end{split}$$

$$\mathcal{C}\_{km\overline{m}l;k} = \frac{1}{\sqrt{2}} r^{-3(1+D)} t^{-3(1+A)} \left( EFr^{3D}t^3 + BCr^3t^{3A} \right.$$

$$+ (-1 + B + C) EFr^{1+2D}t^{2+A} + B\mathcal{C}(-1+E+F)r^{2+D}t^{1+2A} \left( \left( \frac{\Psi\_4}{\Psi\_0} \right)^{1/4} \right)$$

$$\begin{split} \mathbb{C}\_{\text{kun}\mathbb{W};\mathcal{I}} &= \frac{r^{-3(1+D)}t^{-3(1+A)}}{\sqrt{2}\left(\frac{\mathbb{V}\_{4}}{\mathbb{V}\_{0}}\right)^{1/4}} \Big( -EFr^{3D}t^{3} + BCr^{3}t^{3A} \\ &+ (-1+B+C)EFr^{1+2D}t^{2+A} - BC(-1+E+F)r^{2+D}t^{1+2A} \Big) . \end{split}$$

,

$$\begin{split} \mathbb{C}\_{kl\mathfrak{M},m} &= \frac{r^{-3+E+F\_{\mathfrak{f}}-3+B+C}}{4\sqrt{2}\left(\frac{\mathbb{V}\_{4}}{\Psi\_{0}}\right)^{1/4}}(-(B+\mathcal{C}-2BE+(-2+\mathcal{C})\mathcal{CE}-BF-\mathcal{C}(1+\mathcal{C})F) \\ &+ BC(-3+2E+F))r^{E+F}+4BCr^{2(E+F)}t^{-B-\mathcal{C}}+(E-\mathcal{C}E+2BE(-1+F)) \\ &+ F+(-3+\mathcal{C})EF+B(-2+F)F-\mathcal{C}F(1+F))t^{B+\mathcal{C}}-4EFr^{-E-F}t^{2(B+\mathcal{C})}), \end{split}$$

$$\begin{split} \mathbb{C}\_{kl\|\mathbb{M};\mathbb{R}} &= \frac{1}{4\sqrt{2}r^{3}t^{3}\binom{\frac{\mathbf{v}\_{4}}{\mathbf{v}\_{0}}}} \Big( -4\mathbb{C}((-1+\mathbb{C})\mathbb{C} + B(-1+2\mathbb{C}))r^{3(E+F)} \\ &+ (B(-1+\mathbb{C}+2E+4\mathbb{C}E+F-5\mathbb{C}F) - \mathbb{C}(1+F+\mathbb{C}(-2+E+F)))r^{2(E+F)}t^{B+C} \\ &+ (F(1+\mathbb{C}+(-2+B+\mathbb{C})F) - E(-1+\mathbb{C}+F-5\mathbb{C}F+\mathcal{B}(2+4F)))r^{E+F}t^{2(B+C)} \\ &+ 4F((-1+F)F+E(-1+2F))t^{3(B+C)} \Big), \end{split}$$

$$\mathbb{C}\_{\mathsf{Milml};k} = \frac{1}{2\sqrt{2}r^3t^3\binom{\frac{\mathsf{F}\_4}{\mathsf{F}\_0}}}(-2\mathbb{C}((-1+\mathsf{C})\mathsf{C}+B(-1+2\mathsf{C}))r^{3(E+F)}$$

$$+ (\mathbb{C}(-1+2E+\mathsf{C}(2-3E-3F)+F)+B(-1+\mathsf{C}+2E+F-3\mathsf{C}F))r^{2(E+F)}t^{B+C}$$

$$+ (-E(-1+2B+\mathsf{C}+F-3\mathsf{C}F)+F(1-2B-\mathsf{C}-2F+3(B+\mathsf{C})F))r^{E+F}t^{2(B+\mathsf{C})}$$

$$+ 2F((-1+F)F+E(-1+2F))t^{3(B+\mathsf{C})}\Big)\_{\mathsf{F}}$$

$$\begin{split} \mathbf{C}\_{\text{RIM};1} &= \frac{1}{2\sqrt{2}r^{3}t^{3}\left(\frac{\mathbf{V}\_{4}}{\Psi\_{0}}\right)^{3/4}} \Bigl( 2\mathbf{C}\Big(4 - 3\mathbf{B} - 5\mathbf{C} + 2\mathbf{B}\mathbf{C} + \mathbf{C}^{2} \Big) r^{3(E+F)} \\ &+ (-4E + 4F + B(-1 + 6E - 3F + \mathbf{C}(-7 + 4E + 13F)) + \mathbf{C}(15 - 14E - 23F) \\ &+ \mathbf{C}(-14 + 17E + 17F)) (r^{2(E+F)}t^{B+C} + (-E + 15F - 7F(E + 2F) \\ &+ \mathbf{C}(4 + E(-3 + 13F) + F(-23 + 17F)) + B(-4 + E(6 + 4F) \\ &+ F(-14 + 17F))) r^{E+F}t^{2(B+C)} + 2F\Big(4 - 3E - 5F + 2EF + F^{2}\Big)t^{3(B+C)}). \end{split} \tag{50}$$

Here, the non-vanishing CM invariants are:

$$\begin{aligned} r\_1 &= \sqrt{2r\_3} = -\frac{1}{2}(B + \mathbb{C} + E - 3BE - 2CE + F - 2BF - 3CF) \\ &\times (B(D - E) + \mathbb{C}(D - F) + A(E + F))r^{2(E + F - 2)}t^{2(B + \mathbb{C} - 2)}, \end{aligned}$$

$$w\_{1} = \frac{1}{2r^{4}t^{4}} \left( 3\left(B\mathcal{C}r^{2(E+F)} - EFt^{2(B+C)}\right)^{2} + \left(\mathbb{C}(B+2\mathcal{C}-2)r^{2(E+F)}\right) \right)$$

$$+ (\mathbb{C} - E + B(3E-1) + F - 3\mathcal{C}F)r^{E+F}t^{B+C} + F(E+2F-2)t^{2(B+C)}$$

$$+ \left(\mathbb{C}(B+2\mathcal{C}-2)r^{2(E+F)} + (B+E-3BE-F+\mathcal{C}(3F-1))r^{E+F}t^{B+C} + 1\right)$$

$$+ F(E+2F-2)t^{2(B+C)}\Big)\Big),$$

$$w\_{2} = \frac{3}{4r^{6}t^{6}} \left( \left( BCr^{2(E+F)} - EFt^{2(B+C)} \right)^{3} - \left( BCr^{2(E+F)} - EFt^{2(B+C)} \right) \right)$$

$$\times \left( \mathbb{C} (-2+B+2C)r^{2(E+F)} + (\mathbb{C} - E + B(-1+3E)) \right.$$

$$+F-3CF)r^{E+F}t^{B+C} + F(-2+E+2F)t^{2(B+C)}$$

$$+\times \left( \mathbb{C} (-2+B+2C)r^{2(E+F)} + (B+E-3BE-F+\mathbb{C}(-1+3F))r^{E+F}t^{B+C} \right)$$

$$+F(-2+E+2F)t^{2(B+C)} \Big) \Big).$$

$$m\_1 = \frac{1}{4}(B(E-D) + \mathcal{C}(F-D) - A(E+F))$$

$$\times(-E-F+B(-1+3E+2F)+\mathcal{C}(-1+2E+3F))r^{4(-2+E+F)}t^{4(-2+B+C)}$$

$$\times\left(-EFr^{2D}t^2 + B\mathcal{C}r^2t^{2A}\right),$$

$$\begin{split} m\_{2} &= \frac{1}{8} (B + \mathbb{C} + E - 3BE - 2\mathbb{C}E + F - 2BF - 3\mathbb{C}F) \\ &\times (B(D - E) + \mathbb{C}(D - F) + A(E + F))r^{2(-4 + E + F)}t^{2(-4 + B + \mathbb{C})} \\ &\times \left( -\left( B\mathbb{C}r^{2(\mathbb{E} + F)} - EFt^{2(\mathbb{B} + \mathbb{C})} \right)^{2} - \left( \mathbb{C}(-2 + B + 2\mathbb{C})r^{2(\mathbb{E} + F)} \right. \\ &+ (\mathbb{C} - E + B(-1 + 3E) + F - 3CF)r^{\mathbb{E} + F}t^{B + \mathbb{C}} + F(-2 + E + 2F)t^{2(B + \mathbb{C})} \right) \\ &\times \left( \mathbb{C}(-2 + B + 2\mathbb{C})r^{2(\mathbb{E} + F)} + (B + E - 3BE - F + \mathbb{C}(-1 + 3F))r^{\mathbb{E} + F}t^{B + \mathbb{C}} \\ &\quad + F(-2 + E + 2F)t^{2(B + \mathbb{C})} \right) \Big), \end{split}$$

$$\begin{split} m\_{3} &= \frac{1}{16}r^{2(E+F-4)}t^{2(B+C-4)}(-2(B(E-D)+C(F-D)-A(E+F)) \\ \times (-E-F+B(-1+3E+2F)+C(-1+2E+3F))r^{-4D}t^{-4A}\left(EFr^{2D}t^{2}-BCr^{2}t^{2A}\right)^{2} \\ &+ (B+C+E-3BE-2CE+F-2BF-3CF)^{2} \\ \times \left(2C^{2}r^{2(E+F)}+Cr^{E+F}\left((B-2)r^{E+F}+(1-3F)t^{B+C}\right) \\ &+ t^{B+C}\left(B(3E-1)r^{E+F}+F\left(r^{E+F}+2(F-1)t^{B+C}\right)+E\left(-r^{E+F}+Ft^{B+C}\right)\right)\right)^{2} \\ &+ (B(D-E)+C(D-F)+A(E+F))^{2} \\ &\times\left(2C^{2}r^{2(E+F)}+Cr^{E+F}\left((B-2)r^{E+F}+(3F-1)t^{B+C}\right) \\ &+ t^{B+C}\left(-B(3E-1)r^{E+F}-Fr^{E+F}+2(F-1)t^{B+C}+E\left(r^{E+F}+t^{B+C}\right)\right)\right)^{2}\right)^{2} \\ \end{split}$$

*<sup>m</sup>*<sup>5</sup> <sup>=</sup> <sup>1</sup> 8 *r* 2(*E*+*F*−4) *t* 2(*B*+*C*−4) 1 2 (*B*(*E* − *D*) + *C*(*F* − *D*) − *A*(*E* + *F*)) <sup>×</sup>(−*<sup>E</sup>* <sup>−</sup> *<sup>F</sup>* <sup>+</sup> *<sup>B</sup>*(3*<sup>E</sup>* <sup>+</sup> <sup>2</sup>*<sup>F</sup>* <sup>−</sup> <sup>1</sup>) <sup>+</sup> *<sup>C</sup>*(2*<sup>E</sup>* <sup>+</sup> <sup>3</sup>*<sup>F</sup>* <sup>−</sup> <sup>1</sup>))*r*−2−<sup>6</sup>*Dt* −2−6*A BCr*2*t* <sup>2</sup>*<sup>A</sup>* <sup>−</sup> *EFr*<sup>2</sup>*Dt* 2 3 + 1 2 (*B* + *C* + *E* − 3*BE* − 2*CE* + *F* − 2*BF* − 3*CF*) 2 *r* 2(*E*+*F*−2) *t* 2(*B*+*C*−2) × *EFr*<sup>2</sup>*Dt* <sup>2</sup> <sup>−</sup> *BCr*2*<sup>t</sup>* 2*A* 2*C*2*<sup>r</sup>* <sup>2</sup>(*E*+*F*) + *CrE*+*<sup>F</sup>* (*<sup>B</sup>* <sup>−</sup> <sup>2</sup>)*rE*+*<sup>F</sup>* <sup>+</sup> (<sup>1</sup> <sup>−</sup> <sup>3</sup>*F*)*<sup>t</sup> B*+*C* +*t B*+*C <sup>B</sup>*(3*<sup>E</sup>* <sup>−</sup> <sup>1</sup>)*rE*+*<sup>F</sup>* <sup>+</sup> *<sup>F</sup> <sup>r</sup>E*+*<sup>F</sup>* <sup>+</sup> <sup>2</sup>(*<sup>F</sup>* <sup>−</sup> <sup>1</sup>)*<sup>t</sup> B*+*C* + *E* <sup>−</sup>*rE*+*<sup>F</sup>* <sup>+</sup> *FtB*+*<sup>C</sup>* <sup>2</sup> + 1 2 (*B*(*E* − *D*) + *C*(*F* − *D*) − *A*(*E* + *F*))(−*E* − *F* + *B*(3*E* + 2*F* − 1) <sup>+</sup>*C*(2*<sup>E</sup>* <sup>+</sup> <sup>3</sup>*<sup>F</sup>* <sup>−</sup> <sup>1</sup>))*r*−2(1+*D*) *t* −2(1+*A*) <sup>−</sup>*EFr*<sup>2</sup>*Dt* <sup>2</sup> + *BCr*2*t* 2*A* 2*C*2*<sup>r</sup>* 2(*E*+*F*) +*CrE*+*<sup>F</sup>* (*<sup>B</sup>* <sup>−</sup> <sup>2</sup>)*rE*+*<sup>F</sup>* <sup>+</sup> (<sup>1</sup> <sup>−</sup> <sup>3</sup>*F*)*<sup>t</sup> B*+*C* + *t B*+*C <sup>B</sup>*(3*<sup>E</sup>* <sup>−</sup> <sup>1</sup>)*rE*+*<sup>F</sup>* +*F <sup>r</sup>E*+*<sup>F</sup>* <sup>+</sup> <sup>2</sup>(−<sup>1</sup> <sup>+</sup> *<sup>F</sup>*)*<sup>t</sup> B*+*C* + *E* <sup>−</sup>*rE*+*<sup>F</sup>* <sup>+</sup> *FtB*+*<sup>C</sup>* × 2*C*2*r* <sup>2</sup>(*E*+*F*) + *CrE*+*<sup>F</sup>* (−<sup>2</sup> <sup>+</sup> *<sup>B</sup>*)*rE*+*<sup>F</sup>* <sup>+</sup> (3*<sup>F</sup>* <sup>−</sup> <sup>1</sup>)*<sup>t</sup> B*+*C* + *t B*+*C* <sup>−</sup>*B*(3*<sup>E</sup>* <sup>−</sup> <sup>1</sup>)*rE*+*<sup>F</sup>* <sup>−</sup>*FrE*+*<sup>F</sup>* <sup>+</sup> <sup>2</sup>(*<sup>F</sup>* <sup>−</sup> <sup>1</sup>)*FtB*+*<sup>C</sup>* <sup>+</sup> *<sup>E</sup> rE*+*<sup>F</sup>* + *FtB*+*<sup>C</sup>* −1 2 (*B*(*<sup>D</sup>* <sup>−</sup> *<sup>E</sup>*) <sup>+</sup> *<sup>C</sup>*(*<sup>D</sup>* <sup>−</sup> *<sup>F</sup>*) <sup>+</sup> *<sup>A</sup>*(*<sup>E</sup>* <sup>+</sup> *<sup>F</sup>*))<sup>2</sup> *r* 2(−2+*E*+*F*) *t* 2(−2+*B*+*C*) × <sup>−</sup>*EFr*<sup>2</sup>*Dt* <sup>2</sup> + *BCr*2*t* 2*A* 2*C*2*<sup>r</sup>* <sup>2</sup>(*E*+*F*) + *CrE*+*<sup>F</sup>* (−<sup>2</sup> <sup>+</sup> *<sup>B</sup>*)*rE*+*<sup>F</sup>* <sup>+</sup> (−<sup>1</sup> <sup>+</sup> <sup>3</sup>*F*)*<sup>t</sup> B*+*C* +*t B*+*C* <sup>−</sup>*B*(3*<sup>E</sup>* <sup>−</sup> <sup>1</sup>)*rE*+*<sup>F</sup>* <sup>−</sup> *FrE*+*<sup>F</sup>* <sup>+</sup> <sup>2</sup>(*<sup>F</sup>* <sup>−</sup> <sup>1</sup>)*FtB*+*<sup>C</sup>* <sup>+</sup> *<sup>E</sup> rE*+*<sup>F</sup>* + *FtB*+*<sup>C</sup>* <sup>2</sup> , (51)

where we point out that the vanishing of the Ricci invariant *r*<sup>2</sup> and the mixed invariant *m*<sup>4</sup> are useful for identifying this spacetime.

#### **5. The Longitudinal, Quadrantal, and Oblique Levi-Civita Solutions (1918)**

There are three classes of solutions that are less commonly discussed in the literature than the proceeding two we presented. These solutions are derived in [5,6,50] and all three are reviewed and summarized in [6]. These three papers are not published in English and the latter two remain nearly uncited directly, thus we will carefully review these solutions in more depth. Throughout this section we will also present the non-vanishing I invariants (see [14]) as these scalars contain information not present in the CM invariants, due to being constructed from derivatives of the Weyl tensor.

#### *5.1. The Longitudinal Solutions*

Here we characterize the solutions described in [50], which Levi-Civita refers to as longitudinal solutions, described by the metric

$$ds^2 = -(\mu - \epsilon \eta)dt^2 + \frac{d\sigma^2}{\eta^2} + \frac{d\eta^2}{K\_0 \eta^4 (\mu - \epsilon \eta)},\tag{52}$$

where *dσ*<sup>2</sup> is the line element for a two-dimensional space with constant Gaussian curvature and *μ*, *η* are coordinates which run over intervals such that the spacetime is consistent with the chosen signature convention, *μ* is a real constant, *K*<sup>0</sup> is a positive constant, and = ±1.

This solution can be split into three distinct subcases depending on the sign of *μ*. Here, we will move to work with the coordinates given by [4], as they are closer to those typically employed in modern references and explicitly reduce the number of free constants to one. The three forms given are

$$(\mu > 0): \quad ds^2 = -\left(1 - \frac{2m}{r}\right)dt^2 + \left(1 - \frac{2m}{r}\right)^{-1}dr^2 + r^2\left(d\theta^2 + \sin(\theta)^2d\phi^2\right),\tag{53}$$

where *m* > 0 and 2*m* < *r* < ∞, or *m* < 0 and 0 < *r* < ∞,

$$\left(\mu < 0\right): \quad ds^2 = -\left(\frac{2m}{z} - 1\right)dt^2 + \left(\frac{2m}{z} - 1\right)^{-1}dz^2 + z^2\left(dr^2 + \sinh(r)^2d\phi^2\right), \tag{54}$$

where *m* > 0 and 0 < *z* < 2*m*,

$$\left(\mu = 0\right): \quad ds^2 = -\frac{dt^2}{z} + zdz^2 + z^2 \left(dr^2 + r^2d\phi^2\right),\tag{55}$$

where *z* > 0.

For completeness, we note that (55) is not the only solution for *μ* = 0 (we also note that the degenerate static vacuum fields are also listed in Table 2–3.1 in [51], in which Equation (53) is classified as "A1", (54) is classified as "A2", and (55) is classified as "A3", where *b* = 2*m*). In particular, this solution takes Gaussian flat two-spaces to be a twodimensional plane. An equally valid choice would be to take these two spaces to be cylinders,

$$ds^2 = -\frac{dt^2}{z} + zdz^2 + z^2 \left(d\rho^2 + a^2 d\phi^2\right),\tag{56}$$

but these solutions are only different globally, and, thus, both the CK algorithm and CM invariants will not be able to detect this difference.

The solution given by (53) is the Schwarzschild solution which has been invariantly characterized via the CK algorithm in [10], via CM invariants in [52], and via I invariants as a subcase in [14]. The solutions given by (54) and (55) are distinct and, as such, we will explicitly state CK and scalar invariants.

For (54), the null frame:

$$l = \frac{1}{\sqrt{2}} \left( -\sqrt{\frac{2m}{z} - 1} \frac{\partial}{\partial t} + \frac{1}{\sqrt{\frac{2m}{z} - 1}} \frac{\partial}{\partial z} \right),$$

$$k = \frac{1}{\sqrt{2}} \left( -\sqrt{\frac{2m}{z} - 1} \frac{\partial}{\partial t} - \frac{1}{\sqrt{\frac{2m}{z} - 1}} \frac{\partial}{\partial z} \right), \tag{57}$$

$$m = \frac{1}{\sqrt{2}} \left( z \frac{\partial}{\partial r} + iz \sinh(r) \frac{\partial}{\partial \phi} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( z \frac{\partial}{\partial r} - iz \sinh(r) \frac{\partial}{\partial \phi} \right).$$

gives the only non-vanishing curvature component to be

$$
\Psi\_2 = mz^{-3}.\tag{58}
$$

This frame is, therefore, invariant to zeroth order, and the remaining isotropy is given by (31). The independent, first order, non-zero derivatives (components) are

$$\mathbb{C}\_{km\mathbb{I};\mathfrak{M}} = \mathbb{C}\_{km\mathfrak{M};\mathbb{k}} = -\mathbb{C}\_{km\mathfrak{M};\mathbb{l}} = -\mathbb{C}\_{kl\mathfrak{M};\mathbb{m}} = \frac{3m}{\sqrt{2}z^4} \sqrt{\frac{2m}{z} - 1},\tag{59}$$

which are not functionally independent of the zeroth order components. At first order, the isotropy group is reduced to just (33). Since the isotropy group has been reduced at first order, the algorithm proceeds to second order, where the independent non-vanishing terms are

$$\begin{aligned} \mathbf{C}\_{kmkl;lm} &= \mathbf{C}\_{km\mathfrak{M};m\mathfrak{R}} = \mathbf{C}\_{km\mathfrak{M};lm\mathfrak{R}} = \mathbf{C}\_{kl\mathfrak{M};lm} = -\mathbf{C}\_{km\mathfrak{M};\mathfrak{R}\mathfrak{R}} = -\mathbf{C}\_{km\mathfrak{M};\mathfrak{R}\mathfrak{R}} = -\mathbf{C}\_{km\mathfrak{M};\mathfrak{R}\mathfrak{R}} \\ &= -\mathbf{C}\_{km\mathfrak{M};jk} = -\mathbf{C}\_{km\mathfrak{M};jl} = -\mathbf{C}\_{kl\mathfrak{M};lm} = -\mathbf{C}\_{kl\mathfrak{M};ml} = -\mathbf{C}\_{\mathfrak{M}l\mathfrak{M};mm} = \frac{3m(2z-4m)}{z^6}, \\ \mathbf{C}\_{km\mathfrak{M};\tilde{m}\mathfrak{I}} = \mathbf{C}\_{km\tilde{m};lk} = \mathbf{C}\_{km\tilde{m};lk} = \frac{3m(2z-5m)}{z^6}, \end{aligned} \tag{60}$$

which possess the same isotropy and produces no new functionally independent terms. Thus the zeroth, first, and second order terms are the CK invariants needed to fully characterize this spacetime.

In this case, the non-vanishing CM invariants are:

$$w\_1 = -\frac{1}{6\sqrt{w\_2}} = -\frac{z^3}{m'} \tag{61}$$

and the non-vanishing I invariants are

$$I\_1 = 48 \frac{m^2}{z^6}, \quad I\_3 = 720 \frac{m^2(2m-z)}{z^9}, \quad I\_5 = 82944 \frac{m^4(2m-z)}{z^{15}}.\tag{62}$$

We note in this case that the *I*<sup>3</sup> and *I*<sup>5</sup> invariants both vanish for the extremal value *z* = 2*m*, which indicates this hypersurface is invariantly defined in much the same way that the Schwarzschild solution's event horizon is.

For the spacetime given by (55), the algorithm proceeds identically to the preceding case with the null frame given by

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\frac{1}{\sqrt{z}} \frac{\partial}{\partial t} + \sqrt{z} \frac{\partial}{\partial z} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\frac{1}{\sqrt{z}} \frac{\partial}{\partial t} - \sqrt{z} \frac{\partial}{\partial z} \right), \\ m &= \frac{1}{\sqrt{2}} \left( z \frac{\partial}{\partial r} + izr \frac{\partial}{\partial \phi} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( z \frac{\partial}{\partial r} - izr \frac{\partial}{\partial \phi} \right), \end{split} \tag{63}$$

which gives CK invariants at zeroth, first, and second order as:

$$\Psi\_2 = \frac{1}{2}z^{-3},$$

$$\mathbf{C}\_{kmkl;\mathfrak{m}} = \mathbf{C}\_{km\mathfrak{m}l;k} = -\mathbf{C}\_{km\mathfrak{m}l;l} = -\mathbf{C}\_{kl\mathfrak{m}l;m} = \frac{3}{2\sqrt{2}}z^{-9/2},$$

$$\mathbf{C}\_{kmkm;\mathfrak{m}\mathfrak{R}} = \mathbf{C}\_{kmkl;k\mathfrak{m}} = \mathbf{C}\_{kmkl;\mathfrak{m}k} = \mathbf{C}\_{km\mathfrak{R}l;kl} = \mathbf{C}\_{km\mathfrak{R}l;lm} = \mathbf{C}\_{kl\mathfrak{R}l;lm} = \mathbf{C}\_{kl\mathfrak{R}l;ml} \tag{64}$$

$$= \mathbf{C}\_{\mathfrak{R}l\mathfrak{R}l;nm} = -\mathbf{C}\_{kmkl;l\mathfrak{R}} = -\mathbf{C}\_{km\mathfrak{R}l;ml} = -\mathbf{C}\_{km\mathfrak{R}l;ml} = -\mathbf{C}\_{kl\mathfrak{R}l;km}$$

$$= -\frac{4}{5}\mathbf{C}\_{kmkl;\mathfrak{M}l} = -\frac{4}{5}\mathbf{C}\_{km\mathfrak{R}l;kl} = -\frac{4}{5}\mathbf{C}\_{km\mathfrak{R}l;kl} = -\frac{4}{5}\mathbf{C}\_{kl\mathfrak{R}l;mk} = \frac{3}{z^6}.$$

A particularly interesting remark is that this special case coincides (at least locally) with (21) where *h* = −1.

We also have the following non-vanishing CM invariants:

$$w\_1 = \frac{3}{2z^6}, \quad w\_2 = -\frac{3}{4z^9} \tag{65}$$

and I invariants

$$I\_1 = \frac{12}{z^6}, \quad I\_3 = \frac{180}{z^9}, \quad I\_5 = \frac{5184}{z^{15}}.\tag{66}$$

#### *5.2. Quadrantal Solutions*

In [5], Levi-Civita gave a second set of vacuum solutions. These are similar to, but distinct from, the longitudinal solutions (which he called quadrantal solutions). He gave the following line element

$$ds^2 = -\frac{e^{2\zeta(\psi)}}{\xi^2}dt^2 + \frac{1}{K\_0 \xi^2} \left(\frac{d\xi^2}{\Xi(\xi)} + \Xi(\xi)d\phi^2 + d\psi^2\right),\tag{67}$$

where <sup>Ξ</sup>(*ξ*) <sup>=</sup> *μξ*<sup>3</sup> <sup>+</sup> *ξ*2, *<sup>K</sup>*<sup>0</sup> is a positive constant, *<sup>μ</sup>* is a real constant, and <sup>=</sup> <sup>±</sup>1. Note that the coordinate range of *ξ* is restricted to the subset of R<sup>+</sup> such that Ξ(*ξ*) > 0. The undetermined function *ζ*(*ψ*) is defined by the differential equation

$$\frac{\partial^2}{\partial \boldsymbol{\psi}^2} \boldsymbol{e}^{\mathbb{C}(\boldsymbol{\psi})} + \mu \boldsymbol{e}^{\mathbb{C}(\boldsymbol{\psi})} = \mathbf{0},\tag{68}$$

which will have solutions of the form

$$\begin{aligned} \mu > 0 &: \quad \varepsilon^{\zeta}(\psi) = \cos(\sqrt{\mu}\psi), \\ \mu < 0 &: \quad \varepsilon^{\zeta}(\psi) = \cosh(\sqrt{|\mu|}\psi), \\ \mu = 0 &: \quad \varepsilon^{\zeta}(\psi) = \psi. \end{aligned} \tag{69}$$

Working in the null frame

$$\begin{split} l &= \frac{1}{\sqrt{2}} \left( -\frac{e^{\overline{\zeta}(\boldsymbol{\varphi})}}{\overline{\zeta}} \frac{\partial}{\partial t} + \frac{1}{\sqrt{K\_{0}\overline{\zeta}}} \frac{\partial}{\partial \boldsymbol{\psi}} \right), \quad k = \frac{1}{\sqrt{2}} \left( -\frac{e^{\overline{\zeta}(\boldsymbol{\varphi})}}{\overline{\zeta}} \frac{\partial}{\partial t} - \frac{1}{\sqrt{K\_{0}\overline{\zeta}}} \frac{\partial}{\partial \boldsymbol{\psi}} \right), \\ m &= \frac{1}{\sqrt{2}} \left( \frac{1}{\overline{\zeta}\sqrt{K\_{0}\Xi(\boldsymbol{\xi})}} \frac{\partial}{\partial \boldsymbol{\xi}} + \frac{i}{\overline{\zeta}} \sqrt{\frac{\Xi(\boldsymbol{\xi})}{K\_{0}}} \frac{\partial}{\partial \boldsymbol{\psi}} \right), \quad \overline{m} = \frac{1}{\sqrt{2}} \left( \frac{1}{\overline{\zeta}\sqrt{K\_{0}\Xi(\boldsymbol{\xi})}} \frac{\partial}{\partial \boldsymbol{\uppsi}} - \frac{i}{\overline{\zeta}} \sqrt{\frac{\Xi(\boldsymbol{\xi})}{K\_{0}}} \frac{\partial}{\partial \boldsymbol{\psi}} \right), \end{split} \tag{70}$$

the only non-vanishing curvature component is

$$\Psi\_2 = \frac{1}{2} \epsilon K\_0 \mathfrak{J}^3 \, \_\prime \tag{71}$$

which is in an invariant form to first order. The remaining isotropy is given by (31). The independent first derivative components are

$$\mathbf{C}\_{km\mathbf{l}\mathbf{l}\mathbf{l}\mathbf{l}} = \mathbf{C}\_{km\mathbf{l}\mathbf{l}\mathbf{l}\mathbf{m}} = \mathbf{C}\_{km\mathbf{l}\mathbf{l}\mathbf{l}\mathbf{l}} = \mathbf{C}\_{\mathbf{k}\mathbf{l}\mathbf{l}\mathbf{l}\mathbf{l}\mathbf{k}} = \frac{3\epsilon K\_0^{\frac{3}{2}}\xi^4}{2\sqrt{2}}\sqrt{\mu + \epsilon \xi\_\mathbf{v}^2} \tag{72}$$

which are not functionally independent of Ψ<sup>2</sup> and the isotropy is reduced to (35). The second derivative components are

$$\begin{aligned} \mathbf{C}\_{kn\mathbf{m};ll} &= \mathbf{C}\_{kn\mathbf{m}!;lm} = \mathbf{C}\_{kn\mathbf{l}!;m\mathbf{l}} = \mathbf{C}\_{kn\mathbf{m}!;\mathbf{l}\mathbf{l}} = \mathbf{C}\_{kn\mathbf{m}\mathbf{l};jk} = \mathbf{C}\_{kn\mathbf{m}\mathbf{l};jk} = \mathbf{C}\_{km\mathbf{m}!;mm} \\ &= \mathbf{C}\_{kn\mathbf{m}\mathbf{l}!;\mathbf{m}\mathbf{\tilde{n}}} = \mathbf{C}\_{kl\mathbf{\tilde{n}}!;k\mathbf{\tilde{n}}} = \mathbf{C}\_{kl\mathbf{\tilde{n}}!;mk} = \mathbf{C}\_{kl\mathbf{\tilde{n}}!;\mathbf{\tilde{n}}k} = \mathbf{C}\_{\mathbf{\tilde{m}}!l\mathbf{\tilde{n}};jk} = 3\epsilon K\_0^2 \zeta^5 (\mu + \epsilon \xi), \tag{73} \\ \mathbf{C}\_{kn\mathbf{\tilde{n}}!;\mathbf{\tilde{m}}} = \mathbf{C}\_{kn\mathbf{\tilde{m}}!;m\mathbf{\tilde{n}}} = \mathbf{C}\_{km\mathbf{\tilde{m}}!;m\mathbf{\tilde{n}}} = \mathbf{C}\_{kl\mathbf{\tilde{m}}!;km} = \frac{3}{4}\epsilon K\_0^2 \zeta^5 (4\mu + 5\epsilon \xi), \end{aligned}$$

which have the same isotropy and introduce no new functionally independent terms. Thus the algorithm terminates here. In the special case where *μ* = 0 and = 1, this solution will coincide with (21) where *h* = 1/2 or *h* = 2.

In this case, the non-vanishing CM invariants are given by

$$w\_1 = \frac{3}{2} \mathcal{K}\_0^2 \mathfrak{f}^6, \quad w\_2 = -\frac{3}{4} \epsilon \mathcal{K}\_0^3 \mathfrak{f}^9,\tag{74}$$

and the non-vanishing I invariants are given by

$$I\_1 = 12K\_0^2 \xi^6, \quad I\_3 = 180K\_0^3 \xi^8 (\mu + \epsilon \xi^\imath), \quad I\_5 = 5184K\_0^5 \xi^{14} (\mu + \epsilon \xi^\imath). \tag{75}$$

Once again the derivative I invariants will vanish at the extremal value of the coordinate *ξ* (if that coordinate is restricted due to the signs of *μ* and ).

We note that, despite this solution appearing very similar to the longitudinal solutions, the fact that the first order isotropy groups are distinct guarantees that these are inequivalent solutions.

### *5.3. Oblique Solutions*

In [6], a third vacuum solution was given which Levi-Civita described as oblique. This solution is characterized by the line element:

$$ds^2 = -\frac{H(\eta)}{\left(\xi + \eta\right)^2}dt^2 + \frac{1}{K\_0 \xi^2} \left(\frac{d\xi^2}{\Xi(\xi)} + \frac{d\eta^2}{H(\eta)} + \Xi(\xi)d\phi^2\right),\tag{76}$$

where *K*<sup>0</sup> is a positive constant and *H* and Ξ are given by

$$\Xi(\mathfrak{f}) = 4\mathfrak{f}^3 - \mathfrak{g}\_2\mathfrak{f} - \mathfrak{g}\_3, \quad H(\eta) = 4\eta^3 - \mathfrak{g}\_2\eta + \mathfrak{g}\_3. \tag{77}$$

with *g*<sup>2</sup> and *g*<sup>3</sup> being real numbers. Working in the tetrad

$$\begin{split} l &= \frac{1}{\sqrt{2}} \Biggl( -\frac{\sqrt{g\_3 - g\_2}\eta + 4\eta^3}{\eta + \xi} \frac{\partial}{\partial t} + \frac{1}{\left(\eta + \xi\right)\sqrt{K\_0(g\_3 - g\_2\eta + 4\eta^3)}} \frac{\partial}{\partial \eta} \Biggr), \\ k &= \frac{1}{\sqrt{2}} \Biggl( -\frac{\sqrt{g\_3 - g\_2}\eta + 4\eta^3}{\eta + \xi} \frac{\partial}{\partial t} - \frac{1}{\left(\eta + \xi\right)\sqrt{K\_0(g\_3 - g\_2\eta + 4\eta^3)}} \frac{\partial}{\partial \eta} \Biggr), \\ m &= \frac{1}{\sqrt{2}K\_0} \Biggl( \frac{1}{\left(\eta + \xi\right)\sqrt{\left(-g\_3 - g\_2\xi + 4\xi^3\right)}} \frac{\partial}{\partial \xi} + i\frac{\sqrt{-g\_3 - g\_2\xi + 4\xi^3}}{\left(\eta + \xi\right)} \frac{\partial}{\partial \phi} \Biggr), \\ \tilde{m} &= \frac{1}{\sqrt{2}K\_0} \Bigl( \frac{1}{\left(\eta + \xi\right)\sqrt{\left(-g\_3 - g\_2\xi + 4\xi^3\right)}} \frac{\partial}{\partial \xi} - i\frac{\sqrt{-g\_3 - g\_2\xi^3 + 4\xi^3}}{\left(\eta + \xi\right)} \frac{\partial}{\partial \phi} \Bigr). \end{split} \tag{78}$$

the only non-vanishing curvature component is

$$\Psi\_2 = 2K\_0(\eta + \xi)^3. \tag{79}$$

This frame is canonical and has remaining isotropy given by (31) as before. The first order derivative components are

$$\begin{aligned} \mathbf{C}\_{kmkl;l} &= \mathbf{C}\_{km\overline{m}l;m} = \mathbf{C}\_{km\overline{m}l;\overline{m}} = \mathbf{C}\_{kl\overline{m}l;k} = 3\sqrt{2}\sqrt{K\_0^3(\eta+\overline{\xi})^6(-\mathbf{g}\_3 - \mathbf{g}\_2\overline{\xi} + 4\xi^3)},\\ \mathbf{C}\_{km\overline{m}l;l} = \mathbf{C}\_{kl\overline{m}l;m} = -\mathbf{C}\_{km\overline{m}l;\overline{m}} = -\mathbf{C}\_{km\overline{m}l;k} = 3\sqrt{2}\sqrt{K\_0^3(\eta+\overline{\xi})^6(\mathbf{g}\_3 - \mathbf{g}\_2\eta + 4\eta^3)}.\end{aligned} \tag{80}$$

These terms introduce one new functionally independent term. Additionally, since the derivatives are mixed between those seen in (59) and (72), we see that all isotropy will be fixed out at first order. The second order derivatives are

*Ckmkm*;*ll* = *Ckmkl*;*lm* = *Ckmkl*;*ml* = *Ckmml*;*mm* = *Ckmml*;*mm* = *Cklml*;*km* = *Cklml*;*mk* = *Cmlml*;*kk* = 12*K*<sup>2</sup> <sup>0</sup>(*η* + *ξ*) 3 <sup>−</sup>*g*<sup>3</sup> <sup>−</sup> *<sup>g</sup>*2*<sup>ξ</sup>* <sup>+</sup> <sup>4</sup>*ξ*<sup>3</sup> , *Ckmkm*;*mm* = *Ckmkl*;*km* = *Ckmkl*;*mk* = *Ckmml*;*kk* = *Ckmml*;*ll* = *Cklml*;*lm* = *Cklml*;*ml* = *Cmlml*;*mm* = 12*K*<sup>2</sup> <sup>0</sup>(*η* + *ξ*) 3 *<sup>g</sup>*<sup>3</sup> <sup>−</sup> *<sup>g</sup>*2*<sup>η</sup>* <sup>+</sup> <sup>4</sup>*η*<sup>3</sup> , *Cklml*;*mm* = *Cklml*;*mm* = *Cmlml*;*km* = *Cmlml*;*mk* = *Cklml*;*lk* = *Cklml*;*kl* = *Ckmml*;*ml* = *Ckmml*;*lm* = −*Ckmkm*;*lm* = −*Ckmkm*;*mm* = −*Ckmkl*;*kl* = −*Ckmkl*;*lk* = −*Ckmml*;*km* = −*Ckmkl*;*mm* = −*Ckmml*;*mk* = −*Ckmkl*;*mm*

$$\begin{split} &= -C\_{\text{kmfl};lm} = -C\_{\text{kmfl};lm} = -C\_{\text{kmfl};ml} = -C\_{\text{kmfl};ml} \\ &= \frac{3}{2}C\_{\text{kmfl};l} = \frac{3}{2}C\_{\text{kmfl};lm} = \frac{3}{2}C\_{\text{kmfl};ml} = \frac{3}{2}C\_{\text{kfl};mm} \\ &= -\frac{3}{2}C\_{\text{kmfl};m\mathbb{R}} = -\frac{3}{2}C\_{\text{kmfl};lm} = -\frac{3}{2}C\_{\text{kmfl};lm} = -\frac{3}{2}C\_{\text{kfl};lm} \\ &= 12K\_0^2(\eta + \xi)^3 \sqrt{(g\_3 - g\_2\eta + 4\eta^3)(-g\_3 - g\_2\xi + 4\xi^3)}, \\ \mathcal{L}\_{\text{kml};lm} = \mathcal{C}\_{\text{kml};m\mathbb{R}} = \mathcal{C}\_{\text{kml};lm} = \mathcal{C}\_{\text{kfl};lm} \\ &= 3K\_0^2(\eta + \xi)^3 \Big(-6g\_3 + 3g\_2\eta - 3g\_2\xi - 16\eta^3 + 12\eta\xi^2 + 20\xi^3\Big), \\ \mathcal{L}\_{\text{kml};lm} = \mathcal{C}\_{\text{kml};lm} = \mathcal{C}\_{\text{kml};lm} = \mathcal{C}\_{\text{kml};lm} \\ &= 3K\_0^2(\eta + \xi)^3 \Big(-6g\_3 + 3g\_2\eta - 3g\_2\xi + 16\xi^3 - 12\xi\eta^2 - 20\eta^3\Big). \end{split}$$

which introduce no new functionally independent terms, and, thus, the algorithm stops here.

The non-vanishing CM invariants are

$$w\_1 = -\frac{1}{\sqrt{6}} w\_3^{2/3} = 24\mathcal{K}\_0^2 (\eta + \xi)^6,\tag{82}$$

,

and the non-vanishing I invariants are

$$\begin{split} I\_1 &= 192K\_0(\eta + \underline{\varsigma}), \quad I\_3 = -2880K\_0^3(\eta + \underline{\varsigma})^7 \left( \mathfrak{g}\_2 - 4 \left( \eta^2 - \eta \underline{\varsigma} + \underline{\varsigma}^2 \right) \right), \\ I\_5 &= -1327104K\_0^5(\eta + \underline{\varsigma})^{13} \left( \mathfrak{g}\_2 - 4 \left( \eta^2 - \eta \underline{\varsigma} + \underline{\varsigma}^2 \right) \right). \end{split} \tag{83}$$

Note that there are cases where the I invariants will vanish and detect some special surface, although the physical meaning of these surfaces has yet to be understood.

#### **6. Conclusions**

<sup>0</sup>(*η* + *ξ*)

In this paper, we discussed a number of historical solutions originally presented by Tullio Levi-Civita, outlining these solutions in a modern fashion. We also presented two different invariant characterizations of these solutions, using both the CK algorithm and CM invariants. As these methods are fully coordinate invariant they present a unique method of identifying these solutions for future work.

The CK algorithm is used here to generate invariantly defined frames, along with curvature components projected along those frames. These invariant quantities are used to classify these spacetimes and provide insight into certain special cases which can now be seen to be identical. The CK invariants are of particular interest, as this method will extend to higher dimensional generalizations of these solutions in a natural way [19,53]. This paper should also serve as a useful reference for those who want to learn to use this algorithm in their own work, as the numerous cases presented here allow for one to see how the algorithm is applied under a number of different circumstances.

The CM invariants, and in some cases I invariants, for these solutions are also constructed. These invariants provide a coframe independent classification of these solutions (to zeroth order), and also offer a useful characterization in cases where the CK invariants become complicated. In particular, the vanishing of certain invariants can be used to distinguish solutions that might otherwise appear deceivingly similar. The I invariants are also of interest, as there are cases in which these invariants vanish for certain coordinate values. This indicates that there may be interesting, and invariantly defined, surfaces present in these solutions which have not yet been analyzed.

As gravitational wave astronomy moves into the limelight of modern physics, solutions to Einstein's General Theory of Relativity have become increasingly important as a field of study even as an example for proof-by-contradiction. The high non-linearity of Einstein's field equations impose difficulties on finding solutions and an even greater difficulty in the interpretation of them. It is possible that invariantly analyzing these solutions could prove useful in understanding gravitational wave signals from novel sources. Though the 1917 Levi-Civita solution was not a wormhole, cylindrical symmetry seems to be a way to avoid topological censorship and consequently give hope to obtaining phantom-free wormholes, asymptotically flat in the radial direction [45]. Of interest is the possibility that Levi-Civita's metric could provide the pathway towards the engineering of artificial gravitational fields for human spaceflight.

**Author Contributions:** The individual contributions of the authors are: conceptualization, C.K.W. and W.J.; methodology, C.K.W. and W.J.; software, C.K.W. and W.J.; validation, D.D.M.; formal analysis, C.K.W. and W.J.; investigation, C.K.W. and W.J.; resources, E.W.D. and G.B.C.; data curation, C.K.W., W.J.; writing—original draft preparation, C.K.W. and W.J.; writing—review and editing, M.G., D.D.M., E.W.D., and G.B.C; supervision, G.B.C.; project administration, E.W.D. and G.B.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** The authors would like to thank Anzhong Wang for beneficial discussions. The authors would like to thank the reviewers at MDPI symmetry.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

