Mixing “Magnetic” and “Electric” Ehlers–Harrison transformations: The Electromagnetic Swirling Spacetime and Novel Type I Backgrounds

: In this paper, we obtain a complete list of stationary and axisymmetric space-times, generated from a Minkowski spacetime using the Ernst technique. We do so by operating on the associated seed potentials with a composition of Ehlers and Harrison transformations. In particular, assigning an additional “electric” or “magnetic” tag to the transformations, we investigate the new spacetimes obtained either via a composition of magnetic Ehlers and Harrison transformations (first part) or via a magnetic-electric combination (second part). In the first part, the resulting type D spacetime, dubbed electromagnetic swirling universe, features key properties, separately found in swirling and (Bonnor–)Melvin spacetimes, the latter recovered in appropriate limits. A detailed analysis of the geometry is included, and subtle issues are addressed. A detailed proof that the spacetime belongs to the Kundt family, is included, and a notable relation to the planar-Reissner-Nordstr¨om-NUT

1 Introduction Since the advent of Einstein's equations, the quest for exact solutions to this set of coupled nonlinear partial differential equations has played an important part in modern physics [1].Exact solutions have significantly aided us in understanding many classical and semiclassical properties of gravity, and this is exactly what makes them integral to our comprehension of diverse phenomena occurring at astrophysical and cosmological scales.Furthermore, these solutions have set the stage for the theoretical exploration of many groundbreaking concepts, including, but not limited to, black hole thermodynamics [2][3][4][5], the information paradox [6], and holography [7].
Dealing with the field equations of General Relativity (GR) poses a nontrivial challenge, which reduces to a tractable problem only when a high amount of symmetry is imposed.In particular, exploiting the Lie point symmetries inherent in Einstein's equations, serves as a robust method for generating exact solutions-solutions that would be practically impossible to integrate with brute force.Two especially interesting Lie point symmetries of the Einstein-Maxwell system, Ehlers [8,9] and Harrison [10] transformations, suffice for the construction of novel stationary and axially symmetric spacetimes.These two transformations are part of a larger set of Lie point symmetries, which exist as hidden symmetries of the Einstein-Maxwell system of field equations, and which are revealed only when one formulates the theory in terms of the complex so-called Ernst potentials [11,12]; they are indeed potential-space symmetries, the parameters of which comprise an eightparameter isometry group, a nonlinear representation of which was originally given in [13].The linear representation of this group and the apparent isomorphism with SU(2,1), was only a few years later delivered by Kinnersley [14], who extended the previous results of Geroch [15] to include electromagnetism.
Recently, stationary and axisymmetric spacetimes have received a considerable amount of attention within the framework of the Ernst description [11,12] of gravity. 1 In particular, the significant effect that Ehlers and Harrison transformations have on accelerating spacetimes, has been explored in detail [17][18][19][20].Using Ehlers or Harrison transformations, or a combination of them for that matter, it has lately been demonstrated that certain algebraically special accelerating spacetimes can be mapped to novel algebraically general solutions.For example, some of us have successfully constructed a complete hierarchy of type I spacetimes [20] obtained via this generating technique, generalizing the well-known Plebański-Demiański family, i.e., the most general family of type D solutions in Einstein-Maxwell theory.
To better understand these developments, we shall briefly review the effects of these transformations.For seed spacetimes cast into the electric form of the Weyl-Lewis-Papapetrou (WLP) metric,2 the standard lore is that Ehlers transformations introduce an additional real parameter which in certain cases can be associated with a NUT parameter in the target metric.On the other hand, Harrison transformations introduce an additional complex parameter, whose real and imaginary parts, can in certain cases be associated with monopolic electromagnetic charges in the target configuration.When the seed is a type D accelerating spacetime, it has been observed that, on top of the above effects, the transformations also bring about a change in the algebraic character of the generated solution, namely the target spacetimes are of type I.An explanation for this peculiar effect seemingly lies in the way that the transformation parameters enter the new metrics.
To get a good grasp on this, let us for a moment consider the Schwarzschild black hole as a prototypical seed example.Casting it into the electric form of the WLP metric and operating on the seed potentials with an Ehlers transformation, it is known that the resulting spacetime is the Taub-NUT black hole, modulo coordinate transformations, and parameter redefinitions.In other words, the transformation-without changing the Petrov type [17,18]-has introduced a new parameter, which now, together with the mass, determines the location of the black hole horizon.In sharp contrast to this, if we let our seed be the C-metric [21][22][23], an Ehlers transformation will not only affect the usual event horizon in the above sense, but also the Rindler one [19,20].As a byproduct of this, the resulting spacetime turns out to be an accelerating Schwarzschild black hole with a NUTlike parameter, which however is of type I, and as such, it cannot be found within the Plebański-Demiański type D hierarchy.Similarly, a Harrison transformation has its two real parameters entering both, black hole and Rindler, horizons of the C-metric, thereby leading to a charged accelerating black hole of type I and not in a Reissner-Nordström-Cmetric of type D, as one may have perhaps expected.
Important insight into the modified Rindler horizons in these type I accelerating black holes can be obtained by viewing the new solutions as a particular limit of black hole binaries [19].Recall that the near-horizon geometry of a Schwarzschild black hole is described by the Rindler metric, characteristic of an accelerating observer.A mathematically equivalent way to "zoom in" on the event horizon is to take the infinite mass limit of the solution.An accelerating Schwarzschild black hole can be conceptualized as a binary system of two Schwarzschild black holes, effectively described by the Bach-Weyl solution [24], where one of the two grows infinitely large (becomes infinitely massive) while retaining a finite distance from the other.The event horizon of the "big" black hole appears then as an accelerating horizon to its small sibling.Consequently, the Bach-Weyl solution ends up appearing as a C-metric in this limit.Analogously, these type I accelerating black holes, featuring Rindler horizons which depend on the transformation parameters among others, can be thought about as a limit of the NUTty and/or charged extension of the Bach-Weyl spacetime.A complete hierarchy of these novel type I spacetimes, including also a seed angular momentum parameter, can be found in [20]. 3n light of these recent developments, this study aims to shed light on the remaining spacetimes one can generate by composing Ehlers and Harrison transformations.To further elaborate on our agenda, it is best if we use the following terminology, which will be formally introduced in Sec.2.1.An electric (resp.magnetic) Ehlers transformation is an Ehlers transformation of the Ernst potentials associated with a seed metric cast into the electric (resp.magnetic) form of the WLP metric, and ditto for Harrison.To date, and to the best of our knowledge, only the combination of electric transformations has been investigated.In this work, we wish to fill the gaps, and we do so by first combining a magnetic Ehlers transformation with a magnetic Harrison one (first part), and then by taking all possible combinations of an electric and a magnetic transformation (second part).It is a firmly established fact that operating with magnetic Harrison and magnetic Ehlers transformations on the seed potentials of Minkowski spacetime, one obtains two interesting asymptotically nontrivial backgrounds, commonly known as Bonnor-Melvin [25,26] and swirling [27,28] spacetimes, respectively.We will briefly review them later on.In the first part of this work, we present a more general spacetime, dubbed electromagnetic swirling universe, from which, as the name suggests, the aforementioned solutions follow in appropriate limits.We study this geometry in detail, also addressing a rather subtle issue concerning the uniqueness of timelike Killing vectors, relevant also in the swirling case, which has been unfortunately neglected so far.Moreover, we prove that the new metric is Kundt via an intricate chain of coordinate reparametrizations and parameter redefinitions.We also explicitly prove the existence of an intriguing relation between this new background and a planar Reissner-Nordström-NUT spacetime.It is this particular relation which we exploit to also analytically derive the cosmological extension of the electromagnetic swirling universe, i.e., the form of the solution in the presence of a cosmological constant.As a finale to the first part, we embed a Schwarzschild black hole into the new background, giving some emphasis on the dragging effect and the deformation of the horizon surface.
In the second part, we give attention to electric-magnetic mixtures, registering the complete list of novel spacetimes one can generate from Minkowski spacetime by combining an electric with a magnetic transformation.Four type I families are obtainable in this way, and we discuss the conditions under which they can be legitimately called backgrounds.Due to the number of solutions, the analysis will not be exhaustive.We show that all spacetimes in the second part of this work, feature closed timelike/null curves (CTCs/CNCs).These appear inside regions, the boundaries of which are also surfaces where the framedragging angular velocity becomes singular.We then argue that one may perhaps ascribe the occurrence of nonchronal regions to the intensity of rotation building up very close to these singular surfaces; inertial frames get dragged so strongly that the light cones end up being tilted in the direction of the circumference.We remark that two of the backgrounds carry everywhere finite electric and magnetic fields which decay in all directions at infinity.Such behavior is in stark contrast to what happens in the Bonnor-Melvin solution, where the fields are uniform in the vicinity of the symmetry axis.Whether CTC-free parts of these two backgrounds could perhaps be as suitable as the former spacetime, for example in describing astrophysical black holes surrounded by strong magnetic fields, remains an open question.The presence of closed timelike/null curves requires much deeper scrutiny.Although they render the entirety of each solution unrealistic for modeling (astro)physical phenomena, the fact that these pathologies beset solutions to the Einstein-Maxwell theory, perhaps makes these backgrounds interesting from a totally different perspective.
This paper is structured as follows.In Sec.2.1 we communicate the very basics of the Ernst formalism and briefly discuss the symmetry transformations.We use this section to also derive the nomenclature, and we further lay down the steps we follow to generate new solutions, giving an algorithmic description of the generating technique we use in this work.In Sec.2.2, we review, in double-quick time, a highly convenient method for the Petrov classification, which we will use throughout.In Sec. 3, we combine magnetic Ehlers and Harrison transformations to obtain the electromagnetic swirling universe, analyzing its geometric properties and making its relation to a planar Reissner-Nordström-NUT spacetime manifest.We use the latter link to derive the electromagnetic-swirling-(A)dS solution.Finally, a Schwarzschild black hole is embedded into the new background, with its most interesting features discussed in some depth.In Sec. 4, we direct our efforts towards generating new spacetimes by combining electric and magnetic Ehlers and Harrison transformations.We start with a Minkowski seed and work our way up to the four different spacetimes one can obtain by combining a magnetic transformation with an electric one.We investigate whether (and under which conditions) these asymptotically nonflat geome-tries can be characterized as backgrounds.Finally, in Sec. 5 we summarize our findings and conclude, also suggesting new possible avenues for further research.

The complex-potential formalism
In this section, we review the formulation of the Einstein-Maxwell field equations in terms of two complex potentials, as presented by Ernst in his seminal works [11,12].Making a stationary and axisymmetric ansatz, one can introduce two complex potentials, E, and Φ, and write the Einstein-Maxwell field equations as a pair of complex equations.Doing so, a set of Lie point symmetries is revealed, the realization of which eludes one in the usual tensorial formalism.This set of symmetry transformations can then be exploited to generate new solutions from old ones.Let us briefly review the scheme.
The first step is to make a stationary and axisymmetric metric ansatz, the Weyl-Lewis-Papapetrou (WLP) ansatz, where f, ω, and γ are functions of ρ and z.We further take our gauge field to have the symmetry-compatible form where the scalar potentials are also functions of the Weyl coordinates ρ and z.Defining one may show, after some cumbersome algebra, that the Einstein-Maxwell field equations for stationary axisymmetric fields are equivalent to the complex Ernst equations together with a pair of equations determining γ via integration by quadrature. 4 The Laplacian and the gradient are understood as operators in three-dimensional Euclidean space in cylindrical coordinates.The potentials χ and A ϕ are twist potentials satisfying the real equations respectively, where φ is the unit normal in the azimuthal direction.Interestingly, if we define our complex potentials as and considering the metric ansatz instead of (2.1), then for a gauge field as in eq. ( 2.2), we also arrive at eqs. (2.4), with the twist potentials now satisfying the equations (2.8b) As previously mentioned, eqs.(2.4) are invariant under a bunch of symmetry transformations in potential space, whose domain and target (as maps) are potentials associated with stationary and axisymmetric Einstein-Maxwell fields.Their finite forms read [13] G where a, c are real parameters and α, β, ϵ complex.Transformations G 1 and G 2 are "gauge" transformations which transform the potentials, but leave the metric and the gauge field invariant, D is a duality-rescaling transformation, E denotes the Ehlers transformation [8], and H stands for the Harrison [10] one.Note that a composition of G 1 , D and E gives the inversion map which thus is also a (discrete) symmetry of the Ernst equations.In particular, and we can easily verify that The transformations (2.9) are associated with eight Killing vectors (KVs) locally generating an isometry group of the potential space whose linear representation is a representation of SU(2,1) [14].As discussed in [20]  Let us also agree on the terminology to be used in this work.The metric (2.1) will be called the electric form of the WLP, whereas (2.7) will be the magnetic form.These names [29], which unfortunately are misleading, are based on the simple observation that a Harrison transformation (with real parameter) of the potentials (2.3) associated with a vacuum spacetime cast into (2.1),introduces an electric charge, whereas the same transformation acting on the potentials (2.6) associated with the very same vacuum spacetime cast into (2.7),gives new potentials associated with a magnetized version of the seed.Nevertheless, this way of naming things is convenient for the task at hand, and we adhere to it henceforth.Therefore, the potentials (2.3) will be called electric, as they are associated with a seed metric cast into the electric WLP, while the potentials (2.6) will be dubbed magnetic by the same reasoning.Finally, a symmetry transformation of electric potentials will be addressed as electric, and ditto for the magnetic case.
Since this is a solution-generating technique, the process is purely algorithmic.To the aid of the interested reader, we list the steps below for the so-called electric case.The same steps ought to be followed in the magnetic case using the relevant equations.
2. These should be fed to eq. (2.5b), together with the components of the seed gauge field, which can then be integrated for the twist potential A ϕ .Results should in turn be fed to eq. (2.5a) to obtain χ.
3. Substitute everything into the definitions (2.3) to get exact expressions for the seed potentials.Next, do your favorite transformations and read off the target functions f, A t , Ãϕ , and χ again from (2.3).
4. These should now be fed to eq. (2.5a), which can be integrated for ω.Having ω, plug everything into eq.(2.5b) to get the azimuthal component of the target gauge field.
5. Substitute the target f, ω, A t , A ϕ , and the seed γ into the metric (2.1) and the gauge field (2.2), and voilà.
Although this procedure looks simple per se, the computational complexity involved may become quite intractable.

Petrov classification
A fundamental way to distinguish gravitational fields, independent of the coordinate system, is to classify them according to their Petrov type [30].To do so, we need to study the algebraic structure of the Weyl tensor, which in four dimensions has ten independent components, encoded in the five complex Weyl-NP scalars Ψ 0 , . . ., Ψ 4 , within the framework of the Newman-Penrose formalism.The Petrov classification is a useful tool in our case, for one can directly prove the possible nonequivalence between certain solutions appearing in this manuscript, by simply looking at their Petrov type.This being the case, it is worth including a double-quick review of the classification algorithm used in this work.First, we consider an arbitrary complex null tetrad (CNT) e = {k, l, m, m} with k, l real null vectors and m, m complex conjugate null vectors, such that k a l a = −1, m a m a = 1, and all other products zero.Latin indices are lowered/raised with the use of the metric g ab = g µν e µ a e ν b and its inverse.Having a complex null tetrad, we next consider the definition of the five complex Weyl-NP scalar as given in [1].The problem of finding the Petrov type of a given spacetime will be attacked using the d'Inverno and Russell-Clark method [31].
Starting with an arbitrary null basis, we wish to find the specific Lorentz transformation leading to a new basis, in which the number of vanishing Weyl scalars is maximal.For Ψ 4 ̸ = 0, it is known that this is equivalent to the problem of finding the roots of the complex quartic equation Based on the number and multiplicity of these roots, we can then determine the Petrov type.To make things simpler, we will choose our CNT such that Ψ 1 = 0 = Ψ 3 and Ψ 4 ̸ = 0.In App.A, we show that such a choice is always possible for the general WLP metric, and we explicitly suggest the way to construct it.Given this CNT, eq.(2.14) becomes a quadratic for z = ϵ 2 , and dividing it by Ψ 4 , we get with discriminant proportional to 9Ψ 2 2 − Ψ 0 Ψ 4 . (2.16) The two roots of the quadratic are (2.17) Thus, if 9Ψ

The electromagnetic swirling universe
Since this section is dedicated to a sort of composite background, we shall start it with a quick discussion about the separate building blocks of the latter, the electromagnetic universe and the swirling spacetime.The magnetic universe, also known as the Bonnor-Melvin solution, was first found by Bonnor [25], and it was only later rediscovered by Melvin [26].It describes a static and cylindrically symmetric magnetic field immersed in its own gravitational field.In other words, it can be seen as describing a magnetic flux tube held together by its own gravitational pull.The magnetic field lines are parallel to the axis of symmetry, and the field can be treated as a uniform one only near the vicinity of the axis.Since it contributes to the stress tensor, and since stress-energy acts as a gravitating mass, its intensity must be falling off far away from the symmetry axis to prevent a collapse under its own gravity; and this is the case indeed.
Here, we shall present the solution with both, magnetic and electric, external fields.We will refer to it as the electromagnetic universe.In cylindrical coordinates, the metric describing it reads where V (ρ) := 1 + XXρ 2 and X is the complex conjugate of X := (E + iB)/2, with E (resp.B) controlling the intensity of the electric (resp.magnetic) field.This metric is accompanied by the gauge field and it belongs to the Kundt class of Petrov type D electrovac spacetimes.Asymptotically, it approached the Levi-Civita spacetime5 To see this, a rescaling of the noncompact coordinates is necessary.Its motion group is locally generated by the four KVs These KVs do not commute, thus forming a nonabelian Lie algebra g with nonvanishing brackets The center of this Lie algebra is the one-dimensional subspace span{T 3 } which is obviously isomorphic to the real line.This is a solvable algebra with its derived subalgebra being two-dimensional abelian.One can actually observe that g is a trivial extension of e(1, 1), the latter being the Lie algebra of the pseudo-Euclidean group E(1, 1) of rigid motions in Minkowski 2-space.
Another interesting feature of this solution is that, much like what happens in antide Sitter space (AdS), timelike geodesics are forbidden to escape to radial infinity due to the strong attraction towards the axis of symmetry.Yet, here the source of this extreme gravitational pull is not a negative cosmological constant, but rather the electromagnetic field itself (see [26] for the study of geodesic motion in the magnetic universe).Finally, let us remark that the electromagnetic universe can be obtained from a Minkowski seed via a magnetic Harrison transformation with parameter iX.
Definitely, less has been said about the swirling spacetime.This is a stationary vacuum solution of Einstein's field equations, with the above expression first reported in [28], to the best of our knowledge, as an analytic continuation of the Bianchi II cosmological metric of Taub [33].This metric belongs to the Kundt family of Petrov type D vacuum solutions, 6 and its associated isometry group is locally generated by the four KVs Besides them, there is an additional irreducible rank-2 Killing tensor obtained from the Killing-Yano 2-form Since the KVs do not commute, their linear span ought to be a nonabelian Lie algebra.To identify this algebra, it is best if we choose another set of basis vectors, {T + , T − , T, T 3 }, with T ± = (T 1 ± T 4 )/ √ 2 and T = 4jT 2 .Concerning the new basis, the nonvanishing Lie brackets read and it is now easy to see that the derived subalgebra, spanned by {T ± , T }, is the threedimensional Heisenberg algebra.It turns out that the full Lie algebra, with center span{T }, is solvable and nondecomposable.If it bears a special name, then this name unfortunately eludes us.It features as A 4,8 in the classification of four-dimensional Lie algebras by Patera and Winternitz [35].The limit of (3.6) to the Levi-Civita metric is discussed in [27].In the same work, a numeric treatment suggests that geodesic motion is vortex-like. 7The swirling spacetime is free of curvature singularities, a Misner string, and nonchronal regions.Finally, it is interesting to remark that the metric function ω grows infinitely large as |z| → ∞.Being linear in z, it is constant on fixed-z planes and zero on the equatorial plane z = 0, where it changes sign.Do also note that this solution can be obtained from a Minkowski seed via a magnetic Ehlers transformation with parameter j.

The geometry
Let us then construct a new spacetime which features both, an external electromagnetic field and swirling rotation. 8We will create this from a Minkowski seed via a composition of magnetic Ehlers and Harrison transformations, in particular (3.10) where the complex X was defined directly below eq.(3.1).Since this is the first solution we present in this work, we will execute the steps listed in Sec.2.1 one by one.Our seed metric is Minkowski in cylindrical coordinates, This can be cast into the metric (2.7) with nonvanishing seed functions Eqs. (2.8a) and (2.8b) then yield vanishing twist potentials up to the choice of integration constants.The seed potentials are then the simplest possible, E 0 = −ρ 2 and Φ 0 = 0. We act upon them with the transformation (3.10) to obtain the new potentials from which we may read off where V was defined directly below (3.1), and (A ϕ ) (E,B)→(−B,E) denotes the value of A ϕ with B exchanged with E and E exchanged with −B.We use |w| = √ ww for the modulus of a complex variable w.Plugging the above into eq.(2.8a), we obtain a pair of differential equations, first-order in derivatives of ω, which we can integrate for ω = −4jz. (3.15) With ω available, everything shall be fed to eq. (2.8b) which now yields a pair of differential equations, first-order in derivatives of A t , the solution of which reads It follows that the metric describing the electromagnetic swirling universe (EMS) is accompanied by a gauge field It is straightforward to see that when the Harrison parameter vanishes, i.e., X = 0 or equivalently, E = 0 = B, the gauge field vanishes and, taking into account that V = 1 in such a case, we recover the swirling metric (3.6).On the other hand, when the Ehlers parameter vanishes, it is also easy to verify that the resulting spacetime is the electromagnetic universe with metric (3.1) and gauge field (3.2).The metric (3.17) admits a nonabelian group of motions G 4 , locally generated by the KVs in the swirling case, eq.(3.7).Of course, equality at the level of the algebras does not in general imply a group isomorphism (consider covering groups for example).In addition, we have a different Killing-Yano 2-form, which reduces to the Killing-Yano 2-form (3.8) when we switch off the Harrison parameter (after a harmless overall division by j).On the other hand, if we make j vanish, we get a Killing tensor which is just a trivial Killing tensor in the case of the electromagnetic universe.
Let us now have a closer look at the metric (3.17).First of all, observe that the ρ coordinate is not the so-called reduced circumference.The latter reads which goes to zero both as ρ → 0 and ρ → ∞.In fact, its maximum is at a Hence, this would make a very restricted coordinate, and this is why we will stick to the use of the initial ρ coordinate.There are only two metric functions that change sign, g tt and g tϕ .For the former, the surface where the change of sign happens, that is the surface on which g tt = 0, is given by the equation S = 0, where This actually defines two surfaces S ± = 0 (the + for positive z and the − for negative), with on which ∂ t is null, and whose unit normals are spacelike, i.e., n µ ± n ν ± g µν S ± = 1.This means that the surfaces are timelike.Observe that S(ρ, −z) = S(ρ, z), or, equivalently, that S ± (ρ, −z) = S ∓ (ρ, z), and that S ± (ρ, 0) ̸ = 0.This implies that the equator, which is also the timelike surface where g tϕ = 0, acts as a plane of reflection, with S − = 0 in the z < 0 half-space being the mirror image of S + = 0 in the positive z half-space.Since it is clear that for z > 0, S − > 0, and that if then g tt > 0, viz.S + < 0 gives a region in which ∂ t is spacelike.Similarly, for z < 0, S + > 0, and if this provides another region, this time S − < 0, in which ∂ t is again spacelike.These two regions pretty much fulfill the criteria to be called ergoregions, with S ± = 0 giving the ergosurfaces.However, there might be a caveat with this interpretation, which requires further investigation and a deeper understanding.
In the familiar Kerr geometry, ∂ t can be selected as the unique timelike and normalized KV at infinity.Asymptotic flatness of a metric guarantees that the ergoregion (if it exists) is confined.Here, the metric exhibits a peculiar asymptotic behavior (it is basically asymptotically swirling as we will soon see).In fact, it is not hard to see that regions, where ∂ t is spacelike, extend to infinity.Indeed, take z to grow faster than ρ 3 , and notice that the second term in (3.23) becomes dominant, yielding On the other hand, if z grows slower than ρ 3 , it is the first term in S ± that prevails, yielding Therefore, there are regions at infinity where ∂ t is not timelike.Actually, there is simply no such KV (or a linear combination of them) in our case.To see this, consider the most general linear combination ξ = 4 i=1 C i T i , where the T i 's are given in eq.(3.7), with the C i 's being constant coefficients.Being a linear combination of KVs with constant coefficients, this is obviously another KV.Fix a t, do (ρ, z) → (ρϵ −1 , zϵ −4 ), and Taylor expand about ϵ = 0.This ensures that we are probing a case where z grows faster than ρ 3 , as fast as ρ 4 in particular.Doing so, one confronts the following situation: it is impossible to choose the constants {C i } in a way such that the leading-order term in the expansion is negative!In other words, there is a region at asymptotic infinity where no KV can be timelike.Therefore, the concept of t as a time of distant observers or a "time at infinity" seems to be somehow problematic, to say the least.
Well, even in Kerr spacetime, it is true that the interpretation of t as a "time", universal in the entirety of "space" (as a time of distant observers), is meaningful only outside the ergosphere [38].Here, it just happens that the ergoregion unfortunately extends to infinity.
Of course, it is always possible to find a KV which is timelike at ρ infinity provided that z grows slower than ρ 3 .Truly, ∂ t is such a KV, but so are other combinations ξ, e.g.
Clearly, a normalization condition at infinity cannot be used here to single out a unique combination, for there is no ξ satisfying this at all; indeed, ξ One may however demand that ξ • ξ ∼ −1 as ρ → 0 (after all this is a physical region).This forces C 2 1 = C 2 4 + 1, but still leaves C 2 completely arbitrary.Moreover, the condition that ξ • ξ is time independent further fixes C 4 = 0, so we are left with ξ = ∂ t + C 2 ∂ ϕ .There is no other condition, based on limits, which we can use to somehow fix C 2 .Note that one confronts the same situation also in the swirling spacetime.
It is certainly tempting to consider the discrete symmetries of the metric as a means to fix C 2 .Recall that in a Kerr geometry, reflection of time t → −t is not a symmetry unless it is accompanied by a change in the direction of rotation, namely ϕ → −ϕ, and vice versa.This is true also for the EMS metric (3.17), except now we have additional ways to restore time-reversal symmetry.As a matter of fact, time reversal here, if accompanied by a transformation z → −z (which maps one semi-axis to the other), is another (simultaneous) discrete symmetry transformation; it leaves the metric invariant.Do also notice that a transformation ϕ → −ϕ, again accompanied by z → −z, is yet another discrete symmetry.These additional discrete symmetries are a key characteristic of the EMS and swirling metrics.They do not exist for example in the Kerr case.It then appears natural to ponder whether demanding that the norm of the KV candidate is invariant under all the aforementioned discrete symmetries of the metric, namely (t, ϕ) → −(t, ϕ), (t, z) → −(t, z), and (z, ϕ) → −(z, ϕ), is a good condition (on top of the previous ones), one that does the job.This indeed yields C 2 = 0, resulting in ξ = ∂ t .However, this does not prove uniqueness, although it naively appears to do so.
To see this, consider the harmless coordinate transformation (t ′ , z ′ , ϕ ′ ) = (t, z − α, ϕ + 4jαt) where α is some real constant.The transformation is (metric-)form preserving, and therefore, the metric has the same discrete symmetries.Since the inner product ξ • ξ is a coordinate scalar, we may directly write it in terms of the prime coordinates.Now, ∂ t • ∂ t was invariant under (t, z) → −(t, z), but it will not be invariant under (t ′ , z ′ ) → −(t ′ , z ′ ).However, the norm of another KV, namely ∂ t ′ = ∂ t − 4jα∂ ϕ , will.In fact, ∂ t ′ also satisfies all the preceding criteria by default, and we see that discrete symmetries cannot help us single out a candidate KV after all.Our naively "unique" ∂ t is as good as any other member of the family ∂ t − 4jα∂ ϕ .Therefore, in the absence of a robust selection mechanism, we argue that one should indeed use the whole family ξ = ∂ t + C∂ ϕ as the timelike KV, where C is an arbitrary real constant.It then follows that ergosurfaces would be understood as (the timelike) surfaces where ξ is null, ergoregions as regions where ξ is spacelike, and the frame-dragging angular velocity would be given by The latter deserves a few more comments.In particular, consider for a moment the general WLP metric in the form (2.7), for which we have the convenient equation It becomes evident, that using Ω to measure the value of the angular speed at each point is not really a meaningful practice, for C is completely arbitrary.Instead, the meaningful quantity to look at is ∆Ω(ρ, z) := Ω(ρ, z) − Ω 0 = ω(ρ, z) − ω 0 .Remarkably, it then seems that for such spacetimes, in which the timelike KV is the one-parameter family ξ, rotation can only be understood in a relative manner.For example, in the EMS case we currently study, Ω = −C −4jz.It is clear that, taking Ω 0 to be the value of Ω on an arbitrary z slice, this slice will be the surface where the difference ∆Ω changes sign. 9This ultimately implies the presence of counter-rotating regions, though the exact localization of these regions is obviously observer-dependent.After this elucidating aside, let us now, once and for all, choose coordinates adapted to a C = 0 observer (our timelike KV being ξ = ∂ t ), and let us proceed with discussing further features of the solution.
Going to a rectangular coordinate system, one can easily show that the metric has no coordinate singularities.There is no event horizon, and the absence of a conical singularity can be verified by the fact that The absence of a Misner string is also evident since the norm of the azimuthal KV vanishes as ∼ ρ 2 near the symmetry axis.The electromagnetic swirling universe is also free of Closed Timelike Curves (CTCs), for ∂ ϕ is everywhere spacelike.Probing the spacetime for curvature singularities, we shall have a look at R ab cd , i.e., (the components of) the Riemann tensor in the orthonormal basis {e a }, defined in App. A. If this tensor is regular near a coordinate singularity, then the singularity is just due to a poor choice of chart.Indeed, all curvature invariants up to arbitrary polynomial order can be constructed using this particular tensor.Thus, if the tensor itself is regular, the regularity of the invariants follows.On the other hand, if R ab cd is singular near a locus of interest, this does not automatically imply the existence of a curvature singularity, for the poles appearing in the tensor components could, in theory, not appear when taking traces to form curvature invariants.We find that R ab cd depends solely on ρ, and that it is regular everywhere since the denominator of all components is just (V 2 + j 2 ρ 4 ) 3 . 10We also verify that it falls off quite fast as ρ → ∞, which reassures us that tidal forces are diminishing as one moves far away from the axis.To give an example, we mention that the Kretschmann scalar goes as 64(5|X| 4 − 3j 2 ) + O(ρ) in the neighborhood of the symmetry axis, whereas it falls off as ∼ ρ −12 when ρ → ∞.
The new metric is asymptotic to provided that z grows slower than ρ 3 .A coordinate rescaling {t, ρ, z} = |j| brings the above to a form asymptotic to a swirling metric with parameter However, the gauge field strength 2-form does not vanish as ρ ′ → ∞.In fact, Therefore, the complete solution is not asymptotic to the swirling spacetime because the latter is a vacuum solution.Note that * (F ∧ * F ) → 0 as ρ ′ → ∞.
Concerning the Petrov type of the EMS spacetime, it is straightforward to conclude that it is D, for we find that 9Ψ 2 2 = Ψ 0 Ψ 4 (see the reasoning and other details in Sec.2.2).On the contrary, it is not trivial at all to prove that (3.17) actually belongs to Kundt's class.Solutions in the Kundt family admit a shearfree, nonexpanding, and nontwisting null geodesic congruence, with the general metric being [1,32,39,40] where P, H are real functions, and W is complex.Now, consider the specific functions where and where z is supposed to be given in terms of Re ζ via √ 2P 2 d z = dζ + dζ.Let us then perform the coordinate redefinitions in order to express (3.37) in the coordinate system {τ, z, q, ψ}.We get which we readily recognize as a member of the general family of nonexpanding type D solutions (see (16.27) in [32]).At this stage, yet another coordinate transformation with finally brings us to the metric (3.17), and that is all.Therefore, we conclude that the EMS universe is also Kundt, as a combination of two Kundt spacetimes, the swirling one (k = 0 via X = 0) and the EM universe (γ = 0 via j = 0).At the same time, there is a gauge field that we completely neglected so far.To find its form in coordinates {u, v, ζ, ζ}, we start the other way around.Let where k = − E 2 − B 2 , and perform the coordinate transformations to arrive at which is the form of the gauge field (3.18) in the coordinate system {τ, z, q, ψ} with the spacetime metric being (3.41).At this stage, we shall consider which bring us to where and where z is implicitly given in terms of Re ζ via the second equation in the right column of (3.46).From here, one can check that using the backwards coordinate transformations (3.40) and (3.42), together with (3.39) and one indeed reaches eq.(3.18).
Finally, let us have a look at the electric and magnetic fields in the EMS universe.Following the procedure presented in App.B, we find that Both depend only on ρ with field lines parallel to the axis of symmetry, in the vicinity of which they acquire a constant profile, −E ẑ and B ẑ, respectively.Moreover, they fall off as ∼ ρ −4 when ρ grows large, with the extrema of their magnitude given by the positive real roots of a hexic polynomial in ρ.

Double Wick rotation and the planar Reissner-Nordström-NUT spacetime
It is a well-established result that the Bonnor-Melvin solution can be mapped to a planar Reissner-Nordström (RN) spacetime [41].It is also known that the swirling solution can be mapped to a planar Taub-NUT spacetime [27].These mappings are in general achieved by employing a double Wick rotation, coordinate transformations, and parameter redefinitions.Therefore, one may reasonably expect to be able, in a similar fashion, to map the electromagnetic swirling universe, which utterly is a combination of the above, to a planar RN-NUT spacetime.Indeed, we prove that this is the case.First, let us show how to "planarize" the standard RN-NUT metric via a limiting process.Here, the word "standard" refers to the metric and the gauge field where with e and g being the electric and magnetic charge parameters, respectively, l denoting the NUT parameter, and m standing for the mass parameter.We can "planarize" this solution, i.e., "flatten" the S 2 into R 2 , by performing the rescalings and sending λ → 0 with m, e, g, l → ∞ while m, e, g, and l are kept fixed.Here, y = {τ, x, φ}.This procedure results in the metric and in a gauge field equal to (3.52).We have dropped the use of tilde accents for convenience.Notice that the gauge field is invariant under this process.Of course, the result is not necessarily guaranteed to be a solution; this is something that has to be checked.Nevertheless, one can verify that the full solution, comprised of the metric (3.55) and the gauge field (3.52), does indeed satisfy the Einstein-Maxwell field equations.Note that if m > 0 and r > 0, there is a Killing horizon at r = q 2 /(2m) (same as in the case of the planar RN solution), separating the inner region 0 < r < q 2 /(2m) from the outer one, where t is timelike and r is spacelike in the former and the other way around in the latter.
On the other hand, if m < 0 and r > 0, and in contrast to the situation in the planar RN spacetime, the planar RN-NUT does not suffer from a naked singularity exactly due to the presence of the NUT parameter; the curvature scalars are everywhere regular.Of course, the planar RN-NUT solution is also plagued with a Misner string, for the symmetry axis cannot be well-behaved both at x = 1 and x = −1.Now, let us perform a double Wick rotation (t, ϕ) = i( ϕ, t ) of (3.17) and (3.18), also doing E = i E and B = i B such that X := E + i B = iX.After the coordinate transformations and the parameter redefinitions we find out that the resulting metric is exactly the planar RNN metric (3.55).However, to bring the resulting gauge field into the form (3.52), an additional gauge transformation is necessary, the purpose of which is to shift the temporal component by the specific constant 2m(2glm − eq 2 )(4m 2 l 2 + q 4 ) −1 .Only then, the further parameter redefinitions which satisfies the right equation in (3.58), leading to the desired result, namely the mapping of the gauge field in the EMS universe to (3.52).Consequently, we conclude that the full solution can be consistently mapped to a planar RN-NUT spacetime via the above sequence of operations.The various limits are then clear.Killing l is tantamount to switching off j, and vice versa; this provides the (bijective) mapping of the electromagnetic universe to the planar RN spacetime [41].Killing e, g is tantamount to switching off X, and the other way around; this gives the mapping of the swirling solution to the planar Taub-NUT spacetime [27].

Adding a cosmological constant
The previous result motivates one to use the extension of the RN-NUT spacetime for a nonvanishing cosmological constant Λ, to derive a generalization of the EMS solution which includes a Λ ̸ = 0.In general, when a cosmological constant is included, the system of field equations is no longer integrable.Equations that were homogeneous in the absence of Λ, become inhomogeneous in its presence.In particular, the WLP metric itself is not suitable for stationary axisymmetric fields in the presence of a cosmological constant, 11 ergo the machinery used so far is not applicable.This is why the task of extending stationary and axisymmetric solutions of the Ernst equations to account for the presence of Λ, is a highly nontrivial one.
Having said that, let us now attack the problem of generalizing the EMS solution.It can be straightforwardly checked that the form of the planar RN-NUT spacetime in the presence of Λ, is given by (3.55) and (3.52), with Considering the inverse forms of the coordinate transformations (3.57), together with the parameter redefinitions and doing a double Wick rotation ( t, ϕ) = i(ϕ, t) of the resulting spacetime, also setting X = iX, we obtain the generalized metric where Up to gauge transformations, the resulting gauge field is (3.18) as expected.Of course, there is no guarantee that the metric (3.62) and the gauge field (3.18) solve the Einstein-Maxwell-Λ field equations, but we verify that this is the case indeed.Therefore, we have successfully constructed the cosmological extension of the EMS spacetime, which also is of Petrov type D, as well as a member of the Kundt class. 12bserve that for ∂ ϕ to remain spacelike at infinity, we must consider Λ < 0. We also see that there is a spinning string at ρ = 0, for lim ρ→0 Let us then do (t, ϕ) = ( t + a ϕ, ϕ), with a being a constant, together with a reidentification of the new coordinates such that ϕ is 2π periodic, to see whether we can obtain a new spacetime free of it.It is not difficult to check that g ϕ ϕ has a z dependence, and that the axis can be made regular only at a single z (by fixing a), meaning that the spinning string will persist.Therefore, we deduce that it is not artificial (it cannot be removed everywhere along the z axis); it rather corresponds to an actual Misner string!Fortunately, we can remove it by tuning the Ehlers and Harrison parameters as before any regluing.Doing so, we then observe that the induced metric with t = cte = z, after the optional rescaling ρ = ρ 4|X| 2 − Λ/(2|X|), assumes the form close to the symmetry axis. 13The above expression suggests the presence of an infinite strut with negative mass per unit length where 2πδ is the excess angle of the line source.Thankfully, this can be made to vanish via the rescaling if we reidentify ϕ as our new azimuthal coordinate with period 2π.Consequently, it is always possible to obtain a regular spacetime with a negative cosmological constant, free of a Misner string, conical singularities, and CTCs, with the caveat of having the particular relation (3.65) between j and X.
After imposing the tuning (3.65) by replacing j, and with ϕ being our new azimuthal coordinate, the R function in the metric (3.62) acquires the form This is a positive function, with the reduced circumference being (dropping the tilde accent) (3.70) In contrast to R in the case Λ = 0, here the reduced circumference, being a monotonically increasing function of ρ, has the same range as the ρ coordinate.Moreover, the presence of a negative cosmological constant has a significant impact on the full extent of the ergoregions as we see in the instructive Fig. 1.Indeed, if we let both ρ and |z| approach infinity, we find that where j, X are understood as being related via eq.(3.65).In fact, as ρ → ∞, it suffices that |z| ⪆ 1.05/ √ −Λ for ∂ t to be spacelike.Close to the axis, that is as ρ → 0, we have where we let z approach infinity at least as fast as ρ −2 .As in the EMS case, the ergoregions do not touch the z axis.Finally, for completeness, let us mention that the metric (3.62) with Λ > 0, after taking care of the Misner string by imposing eq.(3.65), features a cosmological horizon located at the largest positive real root of a polynomial hexic in ρ.This can be utterly written as a polynomial cubic in ρ 2 which, using Descartes' rule of signs, appears to have a single positive real root if and only if X 2 > Λ/4.This is then a double root of the hexic equation.Alas, the whole situation (with positive Λ) gets more complicated, if we notice that as we approach this root, call it ρ + , the reduced circumference vanishes, meaning that ρ = ρ + behaves as a sort of "axis" besides ρ = 0.This observation is rather expected, simply because the root of g ρρ is necessarily a root of g ϕϕ , as can be seen from the line element (3.62).
Let us also briefly go through the various limits.When X = 0, the gauge field vanishes, and the cosmological EMS metric (3.62) reduces to the cosmological extension of the swirling solution presented in [27].The latter is free of curvature singularities and free of any horizons if Λ < 0. Notwithstanding these good features, a fact completely neglected in [27], is that the cosmological swirling solution actually features a cosmic spinning string, evident from lim which proves to be irremovable through coordinate transformations and regluing.As such, it shall again be understood as a Misner string.Since we no longer have the freedom to tune parameters in order to remove it (as we did previously), one comes to the unfortunate conclusion that the swirling-Λ solution of [27] is in general plagued with a Misner string and the CTCs accompanying it.When j = 0, the metric (3.62) acquires the static form and the gauge field becomes (3.2).We remark that this is not the metric presented in the perhaps pertinent cases [43,44].The above spacetime seems to have a spinning string, for lim ρ→0 g ϕϕ = −Λ|X| −4 /12. (3.75) If Λ > 0, it is impossible to remove this string by any means.If on the other hand Λ < 0, we can do and reglue the spacetime to get rid of it.The new spacetime is then also free of conical singularities.An expansion of the induced metric with t = cte = z near the symmetry axis attests to that.Moreover, ∂ ϕ is clearly spacelike everywhere if we adhere to the use of a negative cosmological constant, while we also remark that there are no curvature singularities, nor are any ergoregions present.

Embedding a Schwarzschild black hole
Having analyzed the background, it is time to discuss the EMS black hole which is just the Schwarzschild black hole embedded in the previously discussed spacetime.For this reason, we will refer to it also as Schwarzschild-EMS.This embedding shall be understood as a composition of magnetic Ehlers and Harrison transformations acting on the potentials associated with a Schwarzschild spacetime.We will work with the spherical-like coordinates {t, r, x = cos θ, ϕ},14 since these prove to be the most convenient for integration.It is quite straightforward to carry out the first steps which result in identifying the Schwarzschild metric with the magnetic WLP metric (2.7), the nonvanishing metric functions being with Weyl's coordinates ρ, z given in terms of r, x via Since this is a static vacuum solution, it is evident that Φ 0 = 0 and E 0 = −f 0 .Observe that the mass does not appear in the potentials, but rather in the function γ and the coordinate transformations (3.79).
Having retrieved the seed data, we now act with the transformation (3.13) on the seed potentials to obtain-after integrating the twist equations-a new spacetime, comprised of the metric with functions together with the gauge field where the function V is now defined as This spacetime describes the exterior of a black hole with a Killing horizon, located at r + = 2M , dressing a singularity at the radial origin.Indeed, g rr (r + , x) = 0, and the Kretschmann scalar blows up as ∼ M 2 r −6 only when r → 0. The surface area of the horizon is that of a 2-sphere of radius r + , namely 4πr 2 + .However, its circumference at the equator x = 0 is not 2πr + in the presence of the transformation parameters, i.e., 2π 0 This suggests that the horizon surface, its embedding in particular, can be visualized as a 2-sphere of area 4πr 2 + , which is deformed when we switch on j and/or X.In particular, the deformation caused by these parameters can be visualized as squeezing the 2-sphere on the equatorial plane, progressively obtaining an egg-shaped surface which is deformed into a peanut-like one as we squeeze stronger.Moreover, the solution is stationary; the black hole gets dragged due to the swirling property of the background.It also got "electromagnetized", in the sense that it is no longer a vacuum solution, but rather an electrovac one, a fact ascribed to the presence of the external electric and magnetic fields in the background geometry.Note that (3.80) enjoys the good features of the EMS background, i.e., it is free of topological singularities and nonchronal regions, 15 while it also shares the same asymptotic behavior with the latter.
Concerning ω, we display various isolines over a heatmap in Fig. 2 (see caption for details).We do so in coordinates {u x , u y , v}, introduced previously.In fact, we first switch to coordinates (3.79) via the inverse transformations In this chart, the horizon surface is understood as a closed line segment of length 2M on the z axis (with center at z = 0).We can formally approach it for |z| ≤ M by taking ρ → 0.Then, we once again employ the convenient coordinates in which |z| infinity sits at |v| = 1 and ρ infinity at u 2 x + u 2 y = 1.The horizon is again understood as a line segment, now of length (4/π) atan M , on the v axis (with center at v = 0).In Fig. 2 we (indirectly) see that |ω| grows infinitely large as |z| → ∞.Since it describes the frame-dragging angular velocity, the fact that it is linear in z, implies that the two z half-spaces (positive and negative) counterrotate.We also observe that it gets smaller and smaller as we approach the equator, where it vanishes.It also decreases as we approach the horizon where it also vanishes in the respective limit.Note that, for j > 0, ω < 0 in the half-space defined by z > 0 (or v > 0), positive otherwise; the exact opposite holds true when j < 0. Regarding the ergoregions, these prove to be more or less insensitive to the introduction of a mass.Actually, their behavior at asymptotic infinity is exactly the same as that in the case of the EMS solution we previously discussed, and there is nothing, really worth reporting, going on close to the horizon.
Attacking the Petrov classification next, and following the reasoning presented in Sec.(2.2), we can deduce that the general Petrov type is I, for we have that 9Ψ 2  2 ̸ = Ψ 0 Ψ 4 .In particular, This becomes zero when M = 0, or j = 0 = X, or at the poles x = 1 and x = −1, or at the horizon r = 2M .In the limit of vanishing mass, we recover the EMS universe (discussed in the previous sections), albeit in the spherical-like coordinate system {t, r, x, ϕ}.When j = 0 = X, we obtain the Schwarzschild seed which is type D everywhere.At the poles, we observe that Ψ 4 ̸ = 0, meaning that the Petrov type is D on the axis (at least in the exterior).On the horizon surface, the spacetime is algebraically special.Since r = 2M turns out to be a pole of Ψ 4 , it would be erroneous to make a definite claim that the Petrov type is D; we can only argue (as we just did) that the solution is algebraically special, because 9Ψ 2 2 − Ψ 0 Ψ 4 , 16 which vanishes at r = 2M , utterly corresponds to a combination of curvature invariants and thus, it is independent of the reference frame.Note that our choice (A.10) of the orthonormal tetrad, from which we constructed the CNT (see App.A for details), corresponds to a zero angular momentum observer in circular motion; it is not suitable for studying the horizon limit.To determine the actual Petrov type at r = 2M , one must rather consider a falling observer [45].Finally, when X = 0, we obtain the swirling black hole, the Petrov type of which is also I, whereas for j = 0 we recover the static type I Schwarzschild-Melvin spacetime [46], which describes a Schwarzschild black hole embedded into an (electro)magnetic universe.The latter solution is of Petrov type II on the horizon surface (again located at r = 2M ), and it also features another locus of interest given by where all Weyl scalars vanish and the Petrov type is, therefore, O [45].

Mixing electric and magnetic transformations
In this section, we wish to scrutinize another possible route.Instead of composing transformations of the same "kind", i.e., either electric or magnetic, we shall explore their mixture.
Once again, we will consider only Ehlers and Harrison transformations.In general, there are eight possible mixed compositions, with β ̸ = δ, where α, . . ., δ = 1, 2 and However, if the seed is Minkowski, the number of available compositions reduces to four, namely This happens because where ∼ denotes a rough equivalence relation here, in the sense that electric transformations of the seed potentials associated with Minkowski space, result in a metric which is also Minkowski modulo coordinate transformations.
Our first course of action is to study the novel spacetimes obtained via these mixed compositions in the case of a Minkowski seed.For starters, we wish to see if these geometries, presented here for the first time, turn out to be backgrounds, i.e., that they are first and foremost free of curvature singularities, topological defects, and other sorts of pathologies.Since we know that where j, c are real parameters, X was previously introduced in the case of the electromagnetic universe, namely X = (E + iB)/2, and we also define a new complex parameter Q := (q e + iq m )/2.This would be tantamount to casting the EMS metric in the form (2.1) and acting on the associated electric potentials with E e [c] • H e [Q].We could then consider all possible limits, in which we remain with two layers of transformations, one magnetic and one electric.Nevertheless, computationally speaking, this would not be a wise strategy to pursue.For this reason, we will generate each case separately.We will probe the new geometries only for curvature singularities, spinning strings, conical singularities, and nonchronal regions (regions with CTCs).We will call them backgrounds if they are free of curvature and topological singularities, proper backgrounds if they are also free of CTCs/CNCs.

Electromagnetic universe and electric Harrison transformations
Here, we operate on the magnetic potentials E 0 = −ρ 2 and Φ 0 = 0, associated with Minkowski spacetime, with 17 Of course, we do not have to truly do the composition; we may directly act on the electric potentials, associated with the electromagnetic universe, with H e [Q].Hence, our seed potentials are where we recall that V = 1 + |X| 2 ρ 2 .We thus follow the prescription presented at the end of Sec.2.1, skipping the integration details to get the target metric (2.1), with where we further defined V(ρ, z) := 1 − |Q| 2 E 0 for convenience.Soon, we will also display the gauge field; before that, let us study this new metric which has Petrov type I.
The first thing that one observes, is that on the equatorial plane z = 0, f has a double (and single) pole at ρ = ρ * , where ρ * denotes the positive real root of V = 0, or |Q| 2 V 2 = 1 equivalently.We find that which is clearly real if-f |Q| ≤ 1. Saturating the bound, this coordinate singularity is put exactly at ρ = 0 = z.For Q < 1, the locus is a ring.For q e B ̸ = −q m E, this is the only pole of f .In any case, if this is a true curvature singularity (be it a point or a ring), and since f ̸ = 0 for all ρ, z, indicating the absence of a horizon, it must be that it is a naked one, with the general consensus being that such configurations are unphysical.Note that the metric does not have any other singular limits besides the one we just reported, which unfortunately happens to be a curvature singularity, for we find that the Kretschmann scalar at the equator blows up as ∼ (ρ − ρ * ) −8 in the limit ρ → ρ * .Parameter tuning cannot be a remedy to this; there are simply not enough parameters to tune in order to get rid of all poles up to octic order in the pertinent Taylor expansion.The situation is worse when |Q| = 1, in which case R µν ρσ R ρσ µν ∼ ρ −16 when ρ → 0. Consequently, it is mandatory to restrict |Q| > 1 in order to expel ρ * from the physical range of ρ.Still, regularity (of curvature invariants) is guaranteed only if the components of the Riemann tensor in the orthonormal basis of App.A are everywhere regular.Fortunately, we find that their denominators are proportional to f −1 , with the proportionality factors being nonvanishing functions of ρ.Since f −1 = 0 has no positive real solutions for |Q| > 1, it follows that the spacetime is indeed everywhere regular.
Having successfully tackled this important issue, it is time to address another subtlety.Observe that at fixed z, where the limit of f as ρ → 0 is a nonvanishing expression that involves the parameters and z, appearing as a denominator of the leading term in the above.Once again, we are confronted with a spinning string, which we can remove by fixing the integration constant C ω as Of course, keeping the integration constant free in this case, only to fix it now, was after all a proactive action.Had we set it to zero, we would again have the string removed via a coordinate transformation, followed by a regluing of spacetime.Nevertheless, the final expression becomes and we now also have well-defined limits X → 0 (Minkowski), or Q → 0 (electromagnetic universe). 18Additionally, we remark that after taking care of the string, the induced metric with t = cte = z becomes ∼ C(dρ 2 + ρ 2 dϕ 2 ) in the vicinity of the symmetry axis, where C denotes a proportionality factor depending on the parameters and the chosen value of z.This tells us that there are no conical singularities to be bothered with.Henceforth, we consider only the family with |Q| > 1 and C ω as in eq.(4.9).For |Q| > 1, the solution is everywhere stationary (not just static) provided that q e B ̸ = −q m E. If q e B + q m E > 0, we have that ω < 0. If the former is negative, the latter is positive.It is also clear that as ρ → 0, it holds that ω → 0. However, the metric In panel [a] we plot a heatmap of (2/π) atan | ω| on the u xv -plane for q e = 1.05, q m = 2.06, E = 1, and B = −1.We do so in coordinates (3.86).The lighter the hue, the greater | ω| is, with white denoting the singular surface.Regions enclosed by the white curves are regions of negative ω.In panel [b] we display a u y = 0 cross-section of the g ϕϕ = 0 surface for q e = 2.07, q m = 2.33, E = −0.31,B = 0.06 (solid), q e = 1.05, q m = 2.06, E = 1, B = −1 (dashed), q e = 0.01, q m = −3.14,E = −0.07,B = 0.92 (dotted).Regions bounded by the curves are nonchronal regions (filled with CTCs).Keep in mind that the full picture is obtained via a complete revolution of these profiles about the v axis, whereby infinity is depicted as a cylinder of unit height and radius.
function ω is no longer the object of interest here, because the frame-dragging velocity (in coordinates adapted to ∂ t ) is actually given by This function is too lengthy to write it down explicitly.However, it is obvious that its zeroes are the zeroes of ω (given that f has no poles), its poles the zeroes of i.e., the surface where the reduced circumference shrinks to zero.Since ω has a fixed sign depending on the choice of q e B+q m E, and since f > 0 everywhere, it follows that the sign of ω may change if and only if g ϕϕ changes sign, namely if there exist CTCs in this spacetime.The boundary of these nonchronal regions, that is the surface on which the norm of ∂ ϕ vanishes, will then be a surface of infinite | ω|.These observations are manifest in Fig. 3. Indeed, observe that the white curves in panel [a] are exactly the dashed ones in panel [b].Note that appearances can be deceiving here, for it looks like the nonchronal regions extend to infinity.This is not true; by doing ρ → ρ/λ, z → z/λ α with α real positive, and expanding about λ = 0, one may check that the leading term is always positive, whichever the value of α is.Interestingly, the fact that the chronology horizons [47] (the surfaces where the nonchronal regions meet the chronal ones) coincide with the surfaces of singular framedragging angular velocity, perhaps admits a physical interpretation.Rotation becomes very rapid very close to these singular surfaces, thereby dragging inertial frames so strongly that the light cones are completely tilted in the direction of the circumference!Such a situation is not unfamiliar generally speaking.A similar interpretation, roughly speaking, appears, for example, in the case of the Van Stockum solution [48,49].Now, observe that if we tune our transformation parameters such that q m E = −q e B, ω vanishes and thus we obtain a static metric, which also is of Petrov type I.The new metric is now free of CTCs, because g ϕϕ = ρ 2 /f with f > 0 everywhere; there is no rotation taking place any longer to tilt the light cones.However, our claim that the spacetime is regular is not valid under the particular tuning.Indeed, the assumptions for that were q m E ̸ = −q e B and |Q| > 1.Unfortunately, when q m E = −q e B, there are two spatial surfaces, one for z < 0 and the other for positive z, on which the curvature invariants become singular.The explicit surface equations are found by requesting the vanishing of the denominator of f .Therefore, although we got rid of rotation and the CTCs, we ended up with a far worse situation, namely a naked singularity with a weird disconnected geometry.It goes without saying that there is no reason to further discuss this scenario, and one should stick to the previous assumptions which at least guarantee regularity.
Let us finally have a look at the gauge field.Its components read modulo gauge transformations.A visual of the field lines (and more) can be found in Fig. 4. We mention that the electric and magnetic fields vanish at asymptotic infinity in all directions.Contrasting this with the behavior of the fields in the electromagnetic universe, in which they are uniform close to the axis for all z, we can argue that the electric and magnetic fields in this spacetime are better-behaved, at least in terms of asymptotic behavior.

Electromagnetic universe and electric Ehlers transformations
Next, let us discuss an alternative possibility that generates a new type I axisymmetric stationary electrovac field, starting again with the electromagnetic universe as our seed.This scenario involves acting with E[c], where c is a real parameter, upon the seed potentials (4.5).Skipping the integration details (we just follow the algorithmic process outlined at the end of Sec.2.1), we obtain the target metric (2.1), with functions together with a gauge field, whose components read It is easy to see that if c ̸ = 0, the Petrov type is I.Of course, for c = 0 the solution reduces to the electromagnetic universe which has Petrov type D. The metric does not exhibit any coordinate singularities, and the axis is not plagued with a spinning string.In particular, the induced metric with t = cte = z behaves as near the symmetry axis, which also proves that there is no conical singularity there.The components R ab cd of the Riemann tensor in the pertinent orthonormal basis of App.A, have a denominator of the general form V n f −k , where n and k are irrelevant positive integers.Since this denominator cannot be made to vanish, we conclude that curvature invariants up to arbitrary polynomial order will be everywhere regular.Moreover, lim ρ→∞ R ab cd = 0 (ditto for z → ∞) further ensures that tidal forces vanish as we move far away from the z axis and/or the equator.We may then call this geometry a background, although not a proper one.
Note that f is strictly positive continuous, meaning that there are no ergoregions in this stationary spacetime.However, there are surfaces where the reduced circumference shrinks to zero.These are, once again, exactly the surfaces where the frame-dragging angular velocity (4.11) blows up.Actually, there are two separate g ϕϕ = 0 surfaces, one in each z half-space, which moreover behave as chronology horizons, in the sense that they separate CTC-free regions from CTC-full ones.Therefore, the interpretation we gave in the previous section applies also here. 19Their structure is more or less similar to the one in the previous case, although now the equator functions as a plane of reflection, simply because g ϕϕ is invariant under z → −z.For this reason, we do not bother plotting them.Now, regarding the electric and magnetic fields, the electric Ehlers transformation has-besides modifying the z component-generated a nontrivial component in the ρ direction which, as expected, is proportional to c. Interestingly, there are surfaces where the radial component of the electric field (with respect to the unit basis) vanishes.These are given by provided that the right hand side is positive.On these surfaces, which is exactly the form of the electric (or magnetic) field in the electromagnetic universe!Note that both, electric and magnetic, fields vanish as ρ → ∞.This was also the case in the seed spacetime.Remarkably, however, here they also vanish far away in the ẑ direction, which was not at all the case in the Bonnor-Melvin solution, where the fields did not depend on z.In particular, we have that

Swirling universe and electric Ehlers transformations
Having explored the scenarios with an electromagnetic universe as our seed, we shall now discuss our options when considering a swirling seed.Let us initiate this discussion with the following composition.In theory, we start with Minkowski space and operate with E e [c] • E m [j] on the associated seed potentials.In practice, we will just act with E[c] upon the seed potentials associated with the swirling spacetime cast into the electric WLP form.Consequently, it is necessary to identify the seed quantities anew; the swirling spacetime can be described by a metric (2.1) with functions ) where we have defined S(ρ) := 1 + j 2 ρ 4 .Given the fact that the seed is stationary, the solution to (2.5a) is a nontrivial seed potential With the previous quantities at hand, the only nonvanishing seed Ernst potential is found to be After acting with E[c], we arrive at the metric (2.1) with functions The former is directly read off from the target potential E since it corresponds to Re E in the absence of Φ.The latter requires integrating eq.(2.5a).Concerning the Petrov type, here we can explicitly write down the relevant expression because it is fairly short, in particular.Therefore, the Petrov type-based on the analysis we did in Sec.2.2of the solution is I.The swirling universe is obtained in the limit c → 0. In this case, expression (4.26) vanishes, but Ψ 4 ̸ = 0. On the other hand, in the limit j → 0, all five complex Weyl-NP scalars become zero, indicating a Petrov type O. Indeed, the solution reduces to flat spacetime.This totally agrees with the fact that electric transformations, when applied to the potentials of Minkowski spacetime, give rise to a target spacetime which is again Minkowski modulo coordinate rescalings.For the case at hand, these rescalings read The new metric is free of a spinning string, for lim ρ→0 g ϕϕ = 0 at fixed arbitrary z.Moreover, it is also free of conical singularities; the induced metric with t = cte = z behaves as ∼ C(ρ 2 + ρ 2 dϕ 2 ) near the symmetry axis, where C is a constant depending on the parameters and the fixed value of z.Note that there are no coordinate singularities evident.However, the components of the Riemann tensor in the orthonormal basis, i.e., R ab cd , come with a denominator of the general form [S 2 − (4jρz , where the exact values of the integers m ≥ 0 < n are utterly unimportant at this stage.When m = 0 (and for some components it is), the denominator can never vanish in the admissible coordinate range.On the other hand, when m ̸ = 0 (true for some components), there are potential poles on the surfaces S ± = 0 (the ergosurfaces as we will soon see), with S ± given in eq.(3.23) for X = 0. Therefore, one cannot argue that the spacetime is everywhere regular, at least not in the fashion we previously did; one needs to compute curvature invariants explicitly.Such behavior is solely due to the swirling nature of the target spacetime, for S ± is independent of c.Nevertheless, calculating the Kretschmann scalar, one finds that the denominator of the latter has m = 0 and n = 6, meaning that it is everywhere regular.Cubic polynomials also have m = 0 (n=9).Consequently, at least up to third-order curvature polynomials (which are coordinate scalars), the absence of singularities is verified.The appearance of poles in higher order polynomials is highly unlikely then, although regularity is not guaranteed in the robust sense of having an everywhere regular R ab cd .Therefore, we may consider this as a background geometry, although we will immediately see that it cannot be proper.
Note that the denominator of f is everywhere positive, meaning that the sign of g tt = −f only depends on the numerator.Its vanishing happens on loci satisfying the surface equation which gives the ergosurfaces.Consequently, the ergosurfaces in this spacetime are exactly the same as the ones in the swirling universe.Finally, let us once again probe for CTCs.By now, the narrative should be clear.For the metric (2.1), the frame-dragging angular velocity is given by (4.11).Rotation becomes infinitely rapid on the surfaces where the denominator vanishes, provided that the zeroes of the denominator are not zeroes of the numerator, or if they are, that the denominator grows faster than ωf 2 close to the surface.Then, these surfaces are necessarily zeroes of g ϕϕ , and if g ϕϕ changes sign there, these act as chronology horizons, separating the chronal from the nonchronal parts of spacetime.In the previous examples of this section, the denominator of g ϕϕ , namely the function f , was strictly positive.Thus, the sign of g ϕϕ solely depended on the sign of the numerator ρ 2 − (ωf ) 2 .
Here, because there are ergoregions, both the numerator and the denominator are allowed to change sign.In fact, we expect that the ergoregions and the nonchronal regions containing CTCs, partially overlap, namely that chronology horizons and ergosurfaces cross each other.Indeed, this can be seen in Fig. 5. Quite interestingly, it turns out that in the spacetime under study, there can actually be up to four disconnected regions filled with CTCs for certain parameter ranges, two toroidal regions with finite volume, and two other regions that extend to infinity, as can be seen from where we let z grow exactly as fast as ρ 2 .Do note that although it seems that the nonchronal "tori" closer to the equator comes into contact with the ergosurfaces, this is not the case.There are only two rings where the two surfaces intersect.Looking at the left panel of Fig. 5, these would be at |v| ≈ 0.47 with radius |u x | ≈ 0.66.

Swirling universe and electric Harrison transformations
The remaining spacetime to consider involves the action of the composition on the seed potentials of Minkowski spacetime.In practice, we just act with H[Q] on the seed potentials (4.24).The integration details are "left as an exercise"; the whole process is already described in the introduction of this work.We obtain the target metric (2.1) with functions By doing this transformation, we have further excited a gauge field with nonvanishing components One can immediately check that 9Ψ 2 2 ̸ = Ψ 0 Ψ 4 with Ψ 4 ̸ = 0 for j, Q ̸ = 0. Therefore, the Petrov type of the solution is I.The swirling universe solution is recovered in the limit Q → 0 (recall that this implies q e , q m → 0).Indeed, also 9Ψ 2 2 = Ψ 0 Ψ 4 in this case, with Ψ 4 ̸ = 0, which gives the Petrov type we expect, that is D. Minkowski spacetime (up to irrelevant coordinate rescalings) is obtained in the limit j → 0. In this limit, all Weyl-NP scalars vanish and the type is O as expected.This spacetime is also free of a spinning string.The fact that g ϕϕ vanishes as ∼ ρ 2 near the symmetry axis, proves the claim.It is also free of conical singularities, for the induced metric with t = cte = z behaves as ∼ C(dρ 2 + ρ 2 dϕ 2 ) in the vicinity of the axis, with C being a constant depending on the parameters and the fixed value of z.What about coordinate singularities?
Let us focus on f and probe it for poles.The function is expected to blow up at solutions to the equation (1 − |Q| 2 f 0 ) 2 + |Q| 4 χ 2 0 = 0.This is an equation quadratic in z 2 , which does not have any real solution unless z = 0. On the equatorial plane, the form is assumed, where ρ := ρ 1/4 .It is clear that the above admits only a single positive real solution if-f |Q| = 1, the solution being ρ = 0.This localizes the coordinate singularity to a single point ρ = 0 = z (the origin).We then need to check whether if this is an honest curvature singularity, or just due to a poor choice of coordinates.Once again, we turn to R ab cd in order to draw conclusions.In the case under study, the components have a denominator of the general form [S 2 − (4jρz 20 Clearly, regularity cannot be directly deduced from R ab cd , because this guy simply blows up on surfaces (4.28) and at the origin ρ, z = 0.This is not problematic per se, as long as the invariants are regular.Upon examining the behavior of the Kretschmann scalar for example, it becomes evident that the thing can be singular only at the origin, which it approaches as ∼ ρ −12 , and only when |Q| = 1.Consequently, if |Q| = 1, the singularity sitting at ρ = 0 = z constitutes a true curvature singularity.However, there is absolutely no null hypersurface to act as a sort of event horizon here, and one thus fails to comply with the censorship hypothesis; the singularity is naked.Moreover, the zeroes of the function f , i.e., the surfaces (4.28), which represent the ergosurfaces in this spacetime, are completely independent of Q (they coincide with the ergosurfaces in the swirling universe), meaning that they will also exist for |Q| = 1.A naked singularity with ergoregions sounds indeed like things really took the wrong turn, and they most probably did.However, such a situation is not that unfamiliar, with the case of the hyperextreme Kerr black hole first coming to mind [32].As previously mentioned, naked singularities are generally regarded as unphysical.Fortunately, all we have to do to avoid this worrying issue, is to simply exclude the value |Q| = 1 from the pool of allowed parameter values.Doing so, we can then argue that spacetime is regular, of course in the less robust sense of the foregoing section, namely that curvature invariants at least up to some low order are regular everywhere.Then, since we have a geometry free of all sorts of "physical" singularities, we can again call this a background.
It turns out that this background cannot get the "proper" attribute, because, once again, we are confronted with CTCs.The nonchronal regions in this spacetime are more or less similar to the ones in the previous case (see Fig. 5); there can be up to four disjoint regions, two in the z > 0 half-space and their mirror images in the negative half-space, with the equator acting as the plane of reflection.As before, closer to the equator we have compact toroidal regions full of CTCs which exist if |Q| > 1 regardless of the value of j.As |z| grows larger, one enters the other nonchronal regions which again extend to infinity, as can be seen from where we take z to grow exactly as fast as ρ 2 .The two regions in each half-space merge exactly when |Q| = 1, which was previously excluded to avoid the curvature singularity at the origin.Concerning the electric and magnetic fields in this electrovac solution, their expressions are quite lengthy, and we prefer to plot them instead.We do so in Fig. 6.There, we see that the fields decay at infinity and that their magnitude is everywhere finite besides on the ergosurfaces, the latter indeed appearing as poles in the expressions of the norms.We mention here that the graphic inspection is certainly not sufficient, and that the behavior of the electric and magnetic fields in this spacetime deserves further investigation which we will unfortunately not include in this work.

Conclusions
This paper focused on obtaining the complete list of stationary and axisymmetric spacetimes, generated from Minkowski spacetime by operating on the associated seed potentials with a composition of Ehlers and Harrison transformations.Adopting the terminology laid down in Sec.2.1, it is clear that the metric functions in the electric and magnetic forms of the WLP metric are completely different functions, though ultimately related via function redefinitions.This implies that the Ernst potentials in the two cases are also going to differ, and thus, any symmetry transformation in potential space is going to provide us also with different target potentials, and ultimately different spacetimes.This is especially clarified in the case of a Minkowski seed, where Ehlers or Harrison transformations of the electric potentials eventually lead us again to seed spacetime (in cylindrical coordinates), modulo rescalings of the noncompact coordinates.On the contrary, an Ehlers or Harrison transformation of the magnetic potentials gives rise to the swirling or electromagnetic universe, respectively.Hence, the form of the WLP metric, electric or magnetic, with which we identify our seed metric, crucially determines the output of the generating process, a deep-seated fact already.
Consequently, taking into account that there are two ways to identify a seed metric and two kinds of transformations we are interested in, we should a priori expect sixteen different spacetimes.However, our seed is Minkowski and thus, this number quickly reduces to ten because a composition of electric transformations does nothing, and while also excluding cases where the resulting spacetime is either the swirling or the electromagnetic universe, we should be expecting at most five distinct target spacetimes.Indeed, one of them was extensively studied in Sec.3.1, and the other four in Sec.4.1.
In Sec. 3, we started by reviewing the electromagnetic universe (the extension of the Bonnor-Melvin solution including an additional external electric field) and the swirling solution.These can be obtained from a Minkowski spacetime by operating on its associated magnetic seed potentials with a Harrison and an Ehlers transformation, respectively.Since our first goal was to combine these two, it was reasonable to expect that the resulting spacetime correctly reduces to its building blocks in the appropriate limits, and that it also inherits properties from both.Indeed, applying a magnetic "Ehlers of a Harrison" map, we were led to the electromagnetic swirling universe, which met the above expectations.We first probed for surfaces where the metric functions change sign.This revealed the presence of a timelike surface, on which ∂ t is null.Had we been able to single out ∂ t as a unique timelike Killing vector, we would call this the ergosurface.However, we pointed out that the usual selection mechanism one uses in Kerr geometry, for example, is not applicable here; there is simply no KV, neither timelike nor normalized, everywhere at infinity.Again, this was not necessarily problematic, since also in the case of Kerr, the definition of t as a "time" is meaningful outside the ergosphere; it just happens that here the ergoregions eventually extend to infinity in particular directions.
We then pondered whether-excluding these regions at infinity where no KV can be timelike-there are other sensible criteria, specific to our case, which can be used to advocate for the existence of a single special timelike Killing vector, to which we should adapt our coordinates.Without being able to rigorously prove that there is no such set of conditions, we nevertheless were not able to find one.The best we could do is to restrict the list of candidates to the one-parameter family ∂ t + C∂ ϕ , with C being an arbitrary real parameter.It was then natural to argue that notions such as ergosurfaces, ergoregions, and frame-dragging angular velocity, should be understood with respect to the complete one-parameter family and not just ∂ t ; the latter is as special as any other member of the former.Moreover, taking the above into consideration, we emphasized that rotation in these spacetimes can only be perceived in a relative manner.The framedragging angular velocity Ω, defined with respect to the timelike KV ∂ t + C∂ ϕ in eq.(3.30), cannot provide an absolute measure, for there is a free parameter C roaming around.However, we remarked that differences (in the value of Ω on two different (ρ, z) surfaces) are independent of C, and indeed, one may use these differences to (i) prove the existence of rotating regions in the first place, and (ii) to further probe for counterrotating regions, with their exact localization being observer-dependent.Notwithstanding this interesting observation, relevant also in the swirling case (though not addressed so far), we eventually decided to adapt our coordinates to a C = 0 observer for clarity and simplicity.However, the particular issue definitely requires further investigation, a task we plan to undertake in future work.
Next, we showed that the EMS spacetime is free of curvature singularities, a Misner string, and conical singularities.We also verified, by checking the components of the Riemann tensor in the suitable orthonormal basis of App.A, that tidal forces diminish at infinity in all directions.Actually, the EMS metric is asymptotic to a swirling metric, given that the growth rate of the ratio ρ/z satisfies a certain inequality.However, the complete solution is not, for the gauge field strength does not in general vanish.We applied the method described in Sec.2.2 to conclude that the Petrov type of the solution is D, and then we gave a detailed sequence of coordinate transformations and parameter redefinitions, which ultimately proves that the EMS spacetime belongs to the Kundt family of solutions admitting a shearfree, nonexpanding, and nontwisting null geodesic congruence.The electric and magnetic fields in this spacetime were found to decay far away from the symmetry axis, but, exactly as in the case of the electromagnetic universe, they were uniform in the proximity of the latter, meaning that they do not vanish far away from the equator near the axis.
Motivated by the well-established relations of the Bonnor-Melvin solution to a planar Reissner-Nordström spacetime, and of the swirling solution to a planar Taub-NUT spacetime, we managed to successfully demonstrate the relation of our EMS solution to a planar Reissner-Nordström-NUT spacetime; we did so by performing a double Wick rotation, and by appropriately redefining our coordinates and parameters.This relation, interesting in its own right, proves to be even more intriguing, if one notices that it can be exploited to directly derive the cosmological extension of the solution.The presence of Λ renders the field equations into a set of inhomogeneous equations, and the system is no longer integrable.The generating technique no longer applies, for the potential-space symmetries are lost.Thus, the only course of action practically is direct integration, which can be a daunting task.Remarkably, starting from a planar Reissner-Nordström-NUT spacetime with a cosmological constant, and applying the previous transformations and redefinitions in reverse order, we arrive at a new spacetime, which for Λ = 0 reduces to the EMS spacetime, and which also is a solution to the field equations of the theory.We study this spacetime in some detail, highlighting the presence of a Misner string which can fortunately be exorcised if we properly tune our parameters.Note that the string exists also in the swirling-Λ solution of [27], despite having been overlooked therein.The crucial difference is that in the latter solution, which can be recovered from the EMS-Λ solution we presented here in the limit X → 0, there is absolutely no way to "banish" it.
Next, we immersed a Schwarzschild black hole in the EMS background.The presence of a mass source in the seed brought in a singularity and an event horizon also in the target spacetime, their locations not modified.We made it clear that whilst the surface area of the horizon does not change, the surface itself is deformed, with the deformation controlled by the transformation parameters.The structure of the ergoregions was qualitatively the same, while no topological singularities or closed timelike curves were present.The asymptotic behavior of the black hole was that of the background.The solution was found to have Petrov type I almost everywhere besides the axis, where the type was D, and the horizon surface, where we proved that it is algebraically special.We remarked that in order to determine the exact Petrov type on the horizon, a frame attached to a freely falling observer is needed.
In the second part of this manuscript, we extracted all the spacetimes one can obtain obtain from Minkowski spacetime, by operating on its associated seed potentials with one electric and one magnetic transformation, Ehlers or Harrison.We showed that the order of operations matters.We excluded cases with U m • U e for reasons previously explained, though we remark that one must consider these when the seed spacetime is other than Minkowski.For example, given a Schwarzschild seed, a combination U m • U e would lead to a completely different spacetime.We registered four novel type I asymptotically nontrivial spacetimes: I.A four-parameter family containing the real parameters E, B, found in the electromagnetic universe, and two additional real parameters q e , q m introduced via the electric Harrison transformation.
II.A three-parameter family containing E, B and one additional real parameter c entering via the electric Ehlers transformation.
III.A two-parameter family containing j and one additional parameter c introduced via the electric Ehlers transformation.
IV.A three-parameter family containing j, found in the swirling spacetime, and two additional parameters q e , q m brought in by the electric Harrison transformation.
As expected, the new spacetimes are complicated modifications of either the electromagnetic universe or the swirling solution.However, these are neither NUTty nor charged, extensions of the latter.To be precise, there is at least no direct evidence to support the association of c with a NUT parameter, or the association of q e , q m with monopolic charges, and the "physical" meaning of these parameters deserves further scrutiny.
Starting with the electrovac spacetime I, here the interaction of q e , q m with the parameters E, B controlling the magnitude of the external electric and magnetic fields in the electromagnetic universe, remarkably produced a stationary spacetime without ergoregions.We argued that for q 2 e + q 2 m ≤ 4, the spacetime describes a naked singularity (a ring singularity on the equatorial plane which can be contracted to a point at ρ = 0 = z when the bound is saturated), otherwise it is regular.We then showcased the presence of a spinning string on the axis which was removable, and once we took care of that, we showed that there are no further topological singularities to be bothered with.Limits of the new geometry to Minkowski spacetime and the electromagnetic universe were checked.By looking at the frame-dragging angular velocity (in coordinates adapted to ∂ t ), we deduced the presence of counterrotating regions and of surfaces where the it cannot even be defined.These last surfaces were then identified with the so-called chronology horizons which bound the nonchronal regions appearing in this spacetime (regions filled with CTCs).Indeed, we claimed that such regions are expected in the configurations we study if the angular velocity is singular on some surface, and a physical interpretation was given, namely that very close to these surfaces, light cones are completely tipped in the direction of angular motion due to the extreme intensity of rotation building up there.The electromagnetic fields in the target spacetime, finite everywhere, were found to decay in all directions at infinity, a behavior contrasting the one in the electromagnetic universe, where the fields appear uniformly close to the axis.
Electrovac spacetime II can be thought of as another modification of the electromagnetic universe.We showed that this is yet another stationary background without ergoregions, but with nonchronal regions which exist for arbitrary values of the involved parameters.The electric Ehlers parameter c has nontrivially modified the gauge field of the electromagnetic universe, with the target electric and magnetic fields acquiring an additional component in the ρ direction.They again were everywhere nonsingular, and they vanished in all directions at infinity.We remarked that there exists a (ρ, z) surface, on which the fields behave exactly like in the electromagnetic universe.Next, we extracted vacuum spacetime III, practically a modification of the swirling universe.We showed that this stationary spacetime inherited its ergoregion structure from the swirling solution and that it was regular, at least in the sense of regular curvature invariants up to some low order.The absence of topological singularities was verified.As in the preceding backgrounds, this one also featured nonchronal regions, with their geometry however being quite different.In fact, we showed that there can be up to four disconnected regions with CTCs.We proved that there are always two of them, which extend to infinity in certain directions, and that for certain values of the parameters, two additional toroidal regions can appear closer to z = 0.The fact that the ergoregions partially overlap with the nonchronal regions, made this particular solution even more puzzling.It was shown that there are exactly two rings where chronology horizons and ergosurfaces meet.Finally, electrovac spacetime IV was registered last.For q 2 e + q 2 m = 4, we found that the new solution, with ergoregions as in the swirling case, describes a naked singularity sitting at ρ = 0 = z, otherwise, it is regular.Topological defects were again absent, but CTCs were present, with their structure being more or less the same as in the case of spacetime III.The expressions of the electromagnetic fields were again too lengthy to explicitly write them down, and thus we plotted them instead.For the corresponding figure, we were able to tell that the fields fall off asymptotically, but that they are not regular everywhere.In particular, we observed that their norm is singular on the ergosurfaces.The specific behavior of the field lines definitely requires a deeper study in order to conclude whether this electromagnetic setup can be eventually of some use.For example, and very roughly speaking, one could say that the plots are somehow reminiscent of two current loops at a given distance, perhaps in the presence of other external electric and magnetic fields, but a convincing explanation of the particular singular surfaces where the fields cannot even be defined definitely eludes us.
We saw that a common factor in all of the new type I solutions in Sec.4.1, was the presence of nonchronal parts.The "physical" interpretation given for the emergence of CTCs in the case of spacetime I, also applies to the other cases; their occurrence seems to be, causally enough, related to the fact that the frame-dragging angular velocity (which would also be the velocity of a zero angular momentum observer) blows up, not in some asymptotic region, but actually in finite regions of these spacetimes.As mentioned in the introduction, the problem of CTCs/CNCs requires a much deeper investigation, i.e., whether they are geodesics or, in any case, whether they are actually traversable under sensible conditions.The whole topic is delicate anyway.Protection mechanisms [47,50] forbidding causality violation have yet to acquire the form of formal theorems, and it also has been conjectured that causality violating curves exist only in the classical theory (and would disappear in the quantum version).Nevertheless, at the classical level these pathologies exist, and the very fact that they beset honest solutions to the Einstein-Maxwell field equations (or even to pure general relativity), namely that they do not arise due to some artificial distribution of matter or some weird gluing of spacetimes, actually makes them interesting in our opinion, though from a very different perspective.
The questions opened are perhaps more than the ones answered, and this implies that there is still room and need for further exploration.It would be interesting to see if and how, for example, the consideration of different seed spacetimes would affect fundamental properties of the target spacetimes, obtained via this combination of electric and magnetic transformations.One simple example is the Schwarzschild black hole embedded into spacetime I, which can be thought of as the counterpart of the Reissner-Nordström black hole immersed in the electromagnetic universe [28], in the sense that the latter can be obtained from a Schwarzschild seed via the combination H m • H e .In the former case, it is expected that the solution, representing a black hole with a deformed horizon, will not be free of CTCs in the exterior, while it will not feature ergoregions at all.On the other hand, nonchronal regions are absent in [41], while ergoregions are present, developing towards asymptotic infinity in some directions.These are two spacetimes with vastly different properties, and with parameters whose physical meaning most probably differs greatly.It is then compelling to deeper understand how and why a simple switch in the order of operations, can bring about so drastically different results, and to shed light on the physical meaning of the parameters in the spacetimes listed here, if any.The existence of a one-parameter family of equally good timelike Killing vectors, instead of a unique one, is also something that requires further research.Moreover, when dealing with purely electric transformations, it is by now understood that any combination of more than two transformations will not yield something new.This does not seem to be the case when mixing transformations, and it is certainly interesting to see if there is a particular number of mixed transformations which gives the same spacetime, regardless of the ordering.Finally, based on the fact that the emergence of CTCs/CNCs proved to be systematic, it is worth investigating-besides scrutinizing the corresponding regions in the spacetimes presented here-whether a set of conditions for the seed spacetime can be mathematically formulated, under which such pathologies can be avoided when applying these transformations.
, Ehlers transformations form a one-dimensional subgroup, i.e., E[a] • E[b] = E[a + b], and they commute with Harrison transformations because [G 1 , G 2 ] = 0. On the other hand, two Harrison transformations do not in general yield a Harrison, simply because

Figure 3 .
Figure 3.In panel[a]  we plot a heatmap of (2/π) atan | ω| on the u xv -plane for q e = 1.05, q m = 2.06, E = 1, and B = −1.We do so in coordinates (3.86).The lighter the hue, the greater | ω| is, with white denoting the singular surface.Regions enclosed by the white curves are regions of negative ω.In panel[b]  we display a u y = 0 cross-section of the g ϕϕ = 0 surface for q e = 2.07, q m = 2.33, E = −0.31,B = 0.06 (solid), q e = 1.05, q m = 2.06, E = 1, B = −1 (dashed), q e = 0.01, q m = −3.14,E = −0.07,B = 0.92 (dotted).Regions bounded by the curves are nonchronal regions (filled with CTCs).Keep in mind that the full picture is obtained via a complete revolution of these profiles about the v axis, whereby infinity is depicted as a cylinder of unit height and radius.

Figure 4 .
Figure 4. Field lines plotted over the heatmap of (a function of) the norm for the electric field (left panel) and the magnetic field (right panel).The size of the arrows is proportional to the magnitude of the vector at each sampled point.The heatmaps show the values of atan |E| (left) and atan |B| (right).The lighter the color, the larger the number.The norms are everywhere finite, and the atan function is used for display purposes.Plots are for q e = −2, q m = 131/151, E = −2/17, and B = −131/2567.These values satisfy eq.(??).

Figure 6 .
Figure 6.Field lines plotted over the heatmap of (a function of) the norm for the electric field (left panel) and the magnetic field (right panel).The heatmaps show the values of atan |E| (left) and atan |B| (right).The lighter the color, the larger the number.Pure white color denotes an infinite vector magnitude.Plots are for q e = −3, q m = −1/2, and j = −1/2.
) so, unless αβ = βα, a "Harrison of Harrison" gives an "Ehlers of Harrison" with the Ehlers parameter fixed in terms of the Harrison parameters α and β.If we recall that the Harrison map generates electrovac solutions from vacuum ones, then eq.(2.13) implies that two Harrison transformations map static vacuum solutions to stationary electrovac ones.