Euclidean Double Kerr–NUT Solution

• The Euclidean Double Kerr–NUT solution is a composite vacuum Einstein geometry that unifies two Kerr–NUT structures through a recursive, nonlinear superposition within the Ernst framework.
• The methodology employs a variation-of-parameters substitution in a stationary axisymmetric setup to transform flat euclidon seeds into curved Kerr–NUT configurations.
• Global regularity is attained by interpreting the solution as an ALF gravitational instanton characterized by a precise rod structure, bridging local geometric constructions with complete manifold descriptions.

The Euclidean Double Kerr–NUT solution denotes a closely related set of exact vacuum Einstein geometries rather than a single universally normalized metric. In the recent euclidon literature, it arises from a nonlinear superposition scheme in which a flat stationary “euclidon” is used as a building block for Kerr–NUT and multi-center Kerr–NUT-type fields; in a separate Euclidean Ricci-flat setting, it appears as an ALF gravitational instanton interpreted as two touching Kerr–NUTs (Shaideman et al., 5 Feb 2025, Shaideman et al., 9 Mar 2026, Chen et al., 2015). Across these formulations, the common theme is the generation of Kerr–NUT structure from composite stationary or Euclidean data, typically within the stationary axisymmetric Ernst framework.

1. Stationary axisymmetric setting

The standard local framework is the Papapetrou form

$$ds^{2}=f^{-1}\!\left[e^{2\gamma}(d\rho^{2}+dz^{2})+\rho^{2}d\varphi^{2}\right]-f(dt-\omega\,d\varphi)^{2},$$

with unknown functions $f(\rho,z)$, $\gamma(\rho,z)$, and $\omega(\rho,z)$. Introducing the twist potential $\Phi$ through

$$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$

the vacuum equations become

$$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$

and can be written in terms of the complex Ernst potential

$\varepsilon=f+i\Phi$

satisfying

$$(\varepsilon+\varepsilon^*)\Delta \varepsilon=2(\nabla\varepsilon)^2 .$$

This formulation is the common base of the euclidon construction, the Kerr–NUT specialization, and the two-NUT superposition picture (Shaideman et al., 9 Mar 2026, 0901.3168).

Within this setting, “Euclidean” enters in two different ways. In the euclidon method, the elementary seed is flat and functions as a stationary non-inertial frame in Minkowski space, although the generated solutions are stationary vacuum metrics. In the Ricci-flat instanton literature, by contrast, the geometry is genuinely Euclidean-signature and is analyzed through rod structure, ALF/ALE asymptotics, and instanton regularity (Shaideman et al., 5 Feb 2025, Chen et al., 2015).

2. The stationary euclidon and the variation-of-parameters method

The basic one-center stationary euclidon is given by

$$f=\frac{(z-z_i)+r_i\tanh U_0}{C_1}, \qquad \Phi=\frac{r_i}{C_1\cosh U_0}+C_2,$$

$f(\rho,z)$0

with constants $f(\rho,z)$1. The defining property emphasized in the euclidon papers is that all Riemann tensor components vanish, so the metric is flat. It is therefore interpreted not as a gravitating source but as a stationary representation of a relativistic noninertial frame of reference in flat space-time (Shaideman et al., 9 Mar 2026, Shaideman et al., 5 Feb 2025).

The constructive step is a variation-of-parameters substitution. One replaces

$f(\rho,z)$2

where $f(\rho,z)$3 is a stationary vacuum seed. This yields

$f(\rho,z)$4

$f(\rho,z)$5

The function $f(\rho,z)$6 satisfies a coupled first-order system; its integrability condition is automatically satisfied if the seed solves the Einstein equations. A useful linearization is

$f(\rho,z)$7

In this sense, the euclidon method converts a nonlinear composition problem into an integrable parameter-variation problem organized by the auxiliary variables $f(\rho,z)$8 and $f(\rho,z)$9 (Shaideman et al., 9 Mar 2026).

This construction is central to the double Kerr–NUT topic because the euclidon is not introduced as an isolated curiosity. It is explicitly treated as a building block from which Kerr–NUT and higher multi-center stationary configurations are generated by recursion (Shaideman et al., 5 Feb 2025).

3. From two euclidons to Kerr–NUT and double Kerr–NUT

The two-stationary-euclidon solution is the first nontrivial application. Its Ernst potential is written as

$\gamma(\rho,z)$0

with

$\gamma(\rho,z)$1

A closed expression for the intermediate function $\gamma(\rho,z)$2 is also given, and the construction becomes especially transparent after imposing

$\gamma(\rho,z)$3

and introducing prolate spheroidal coordinates

$\gamma(\rho,z)$4

With the parameter identifications written in equation (5.6) of the 2026 paper, the resulting metric becomes the Kerr–NUT family; the limit $\gamma(\rho,z)$5 with $\gamma(\rho,z)$6 yields Kerr (Shaideman et al., 9 Mar 2026).

The 2025 paper presents the same point in slightly different language: the stationary two-euclidon solution coincides with the Kerr–NUT solution in form, and the Kerr solution is recovered by setting

$\gamma(\rho,z)$7

That statement makes the intended interpretation explicit: Kerr–NUT is not treated as an isolated exact solution but as the two-body member of a hierarchy generated from the stationary euclidon seed (Shaideman et al., 5 Feb 2025).

The same recursive machinery then produces higher composites. In the $\gamma(\rho,z)$8-center construction, one starts from a Zipoy-like static seed and iterates the euclidon addition law. The recursive form is

$\gamma(\rho,z)$9

$\omega(\rho,z)$0

$\omega(\rho,z)$1

The authors state that, when $\omega(\rho,z)$2, $\omega(\rho,z)$3, $\omega(\rho,z)$4, and $\omega(\rho,z)$5, the construction becomes an $\omega(\rho,z)$6-center solution describing $\omega(\rho,z)$7 rotating axially symmetric masses. In the absence of rotation it reduces to an $\omega(\rho,z)$8-center Zipoy-type configuration, and without distortion it becomes a collection of $\omega(\rho,z)$9 Kerr–NUT-type sources (Shaideman et al., 9 Mar 2026).

Within this hierarchy, the explicit “double Kerr–NUT” case is the two-center specialization

$\Phi$0

with

$\Phi$1

For

$\Phi$2

the resulting metric functions become the two-Kerr–NUT solution

$\Phi$3

and if

$\Phi$4

they reduce to the two-soliton Kerr–NUT form (Shaideman et al., 9 Mar 2026). In this literature, “double Kerr–NUT” therefore means a two-center rotating, NUT-charged composite obtained as a four-euclidon specialization of the recursive construction.

4. Euclidean Ricci-flat instanton as two touching Kerr–NUTs

A distinct but directly relevant Euclidean realization is the five-parameter Ricci-flat solution with Euclidean signature given in C-metric-like coordinates $\Phi$5. It is asymptotically locally flat, has a finite asymptotic NUT charge, and becomes asymptotically locally Euclidean when that charge is sent to infinity. The authors interpret it as a system consisting of two touching Kerr–NUTs: the south pole of one Kerr–NUT touches the north pole of the other (Chen et al., 2015).

The metric possesses four rods and three turning points. In Weyl–Papapetrou coordinates,

$\Phi$6

the turning points are located at

$\Phi$7

with corresponding Weyl $\Phi$8-coordinates $\Phi$9. In natural Killing coordinates $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$0, the rod directions are

$$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$1

This is identified as the rod pattern of two Kerr–NUTs placed one above the other, each with its own horizon rod, with the outer rods carrying opposite asymptotic NUT charges $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$2 (Chen et al., 2015).

The geometric interpretation is precise. The first Kerr–NUT occupies rods $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$3, the second occupies rods $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$4, and the common middle structure means that the south pole of one touches the north pole of the other. The authors relate this to the inverse-scattering construction: the seed is the double-Schwarzschild geometry, and the BZ transformation eliminates the inner axis by joining it to a horizon. The result is therefore not two separated Euclidean Kerr–NUTs with an intervening axis segment, but a touching configuration (Chen et al., 2015).

The asymptotic NUT charge $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$5 is a global invariant of this ALF instanton, and the local turning-point charges satisfy

$$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$6

The principal limiting cases are also explicit: $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$7 gives an AF solution, $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$8 gives an ALE limit, $$\frac{\partial \omega}{\partial \rho}=\frac{\rho}{f^2}\frac{\partial \Phi}{\partial z}, \qquad \frac{\partial \omega}{\partial z}=-\frac{\rho}{f^2}\frac{\partial \Phi}{\partial \rho},$$9 yields the Ricci-flat Plebański–Demiański solution, and lower-turning-point degenerations include Kerr–NUT and Taub–NUT limits (Chen et al., 2015).

5. Global regularity, Euclideanization, and complex-geometric structure

The Euclidean Double Kerr–NUT topic is not exhausted by local line elements. A separate issue is how such geometries are globally completed and how Euclidean continuation is encoded.

For Kerr–NUT–(A)dS, a non-singular extension is obtained by replacing naive Misner compactification with a principal $$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$0-bundle construction. The metric admits a two-dimensional algebra of commuting Killing fields generated by $$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$1 and $$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$2, and for $$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$3 there are exactly two admissible bundle generators,

$$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$4

with

$$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$5

Using either branch, the construction glues two charts with a nontrivial transition function and produces a globally defined manifold diffeomorphic to

$$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$6

That work does not explicitly construct a Euclidean Double Kerr–NUT instanton, but it provides the relevant global bundle mechanism and clarifies that non-singular completion is not a naive single identification along a Misner fiber (Lewandowski et al., 2021).

A complementary perspective comes from complex geometry. In the Euclidean continuation of the vacuum Plebański–Demiański family, the geometry carries two commuting complex structures $$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$7 and $$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$8 of opposite orientation together with two commuting Killing vector fields. The resulting linear-algebraic constraints force the metric into the ambitoric ansatz

$$f\Delta f=(\nabla f)^2-(\nabla \Phi)^2, \qquad \nabla\!\cdot\!(f^{-2}\nabla\Phi)=0,$$9

which is then reduced by the Einstein equations to the Euclidean Plebański–Demiański family. Kerr is recovered as a two-parameter subfamily by setting

$\varepsilon=f+i\Phi$0

The paper explicitly remarks that this ansatz is “half-way” to the Euclidean double Kerr–NUT / Plebański–Demiański form (Krasnov et al., 2024).

Taken together, these works show that the Euclidean Double Kerr–NUT problem has at least three layers: local exact metric form, nonlinear composition law, and global/topological completion.

6. Terminology, related constructions, and common confusions

The expression “Euclidean Double Kerr–NUT” is used across adjacent but non-identical constructions. The following summary organizes the main usages already present in the literature cited here.

Framework	Exact object	Relation to the term
Euclidon recursion	Two-euclidon Kerr–NUT; four-euclidon/two-center two-Kerr–NUT	Constructive stationary-vacuum usage (Shaideman et al., 5 Feb 2025, Shaideman et al., 9 Mar 2026)
Euclidean ALF instanton	Five-parameter Ricci-flat solution interpreted as two touching Kerr–NUTs	Direct Euclidean instanton realization (Chen et al., 2015)
Symmetric two-NUT superposition	Two identical counter-rotating NUT objects yielding Kerr for $\varepsilon=f+i\Phi$1	Conceptual analogue, not explicitly Euclidean (0901.3168)
Non-singular Kerr–NUT–(A)dS extension	Principal $\varepsilon=f+i\Phi$2-bundle completion with two admissible generators	Global mechanism, not a direct double instanton (Lewandowski et al., 2021)

A frequent source of confusion is the assumption that every double Kerr–NUT construction is Euclidean and globally regular. The 2009 two-source solution is explicitly a Lorentzian stationary vacuum metric built from axis data via Sibgatullin’s method. Its axis data are

$\varepsilon=f+i\Phi$3

describing two identical counter-rotating NUT objects placed symmetrically at $\varepsilon=f+i\Phi$4, and the special choice

$\varepsilon=f+i\Phi$5

makes the metric collapse exactly to Kerr (0901.3168). The paper itself states that it does not explicitly derive the construction as a Euclidean double Kerr–NUT solution and does not use a Euclidean continuation.

Another common misunderstanding is to treat the euclidon as a gravitating source. In the euclidon method, the seed is explicitly flat: its apparent gravitational structure comes from its stationary, non-inertial coordinate form, not from curvature. The genuinely curved solutions arise only after the nonlinear composition with a seed vacuum field (Shaideman et al., 5 Feb 2025, Shaideman et al., 9 Mar 2026).

This suggests that the term “Euclidean Double Kerr–NUT solution” is best read as a family resemblance label. In one branch it names a Euclidean Ricci-flat instanton interpretable by rod structure as two touching Kerr–NUTs; in another it denotes the two-center Kerr–NUT composite generated by the euclidon algebra; and in a third it functions as a conceptual bridge between symmetric double-NUT superpositions and the Kerr limit.