Kerr-Taub-NUT Spacetime Analysis

• Kerr-Taub-NUT spacetime is a general-relativistic solution characterized by mass, rotation (Kerr parameter), and a NUT parameter that introduces novel gravitomagnetic effects.
• The metric exhibits distinctive horizon and ergosphere structures, with shifts in event horizon radius and ergosphere volume modulated by the NUT charge.
• Astrophysical implications include modified geodesic motion, altered accretion flow dynamics, non-monotonic frame-dragging, and varied energy extraction via the Penrose process.

Searching arXiv for recent and foundational papers on Kerr-Taub-NUT spacetime to support the article. Kerr-Taub-NUT spacetime is a stationary, axisymmetric vacuum solution to Einstein’s equations characterized by the mass $M$, the Kerr parameter $a=J/M$, and a NUT parameter $n$ or $l$, often described in the literature as a “magnetic mass,” “dual mass,” or gravitomagnetic charge. It reduces to Kerr for $n=0$, to Taub-NUT for $a=0$, and to Schwarzschild for $a=n=0$. The NUT parameter introduces additional gravitomagnetic effects with no Newtonian analog, and the physical interpretation of this parameter remains subtle and debated; the spacetime is also associated with asymptotic non-flatness and Misner-string structure in standard presentations (Chakraborty, 2014, Abdujabbarov et al., 2011, Liu et al., 2010).

1. Metric structure, parameters, and horizons

A standard Boyer-Lindquist-type form of the Kerr-Taub-NUT metric is

$$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$

with

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$

$$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$

An equivalent notation used in other analyses writes

$a=J/M$0

so that the metric coefficients differ only by symbol choice for the NUT parameter (Chakraborty, 2014, Abdujabbarov et al., 2011).

The horizon structure is determined by the roots of $a=J/M$1. In the $a=J/M$2-notation,

$a=J/M$3

and the static limit or ergosurface is

$a=J/M$4

In the $a=J/M$5-notation, the same relations are written as

$a=J/M$6

One study emphasizes that increasing $a=J/M$7 increases the radius of the event horizon but decreases the relative volume of the ergosphere, whereas a later repetitive-Penrose analysis states that the event horizon and ergosphere radii increase under the influence of the gravitomagnetic charge $a=J/M$8 (Abdujabbarov et al., 2011, Alloqulov et al., 26 Mar 2026). A special case discussed in the Lense-Thirring analysis is $a=J/M$9, equivalently $n$0, for which

$n$1

so the inner horizon moves to $n$2 while the event horizon remains at $n$3 (Chakraborty, 2014).

Global regularity is nontrivial. A non-singular extension of the broader Kerr-NUT-(anti) de Sitter family was constructed by refining Misner’s method through a $n$4-principal bundle structure over the non-singular space of non-null orbits, with two distinguished Killing vector fields selected by regularity at both poles (Lewandowski et al., 2021). Within that framework, the only possible residual singularity is Kerr-like and occurs when $n$5 (Lewandowski et al., 2021).

2. Geodesics, effective potentials, and accretion structure

For quasi-equatorial timelike motion, one form of the radial equation is

$n$6

with effective potential

$n$7

where, on the equatorial plane, $n$8 and $n$9. The NUT parameter appears explicitly in all terms. The corresponding analysis reports that increasing $l$0 slightly lowers the effective potential, which can facilitate access to negative-energy orbits inside the ergosphere (Abdujabbarov et al., 2011).

The equatorial circular-geodesic problem has a particularly sharp result in the extremal case. For extremal Kerr-Taub-NUT, defined by $l$1, both null and timelike ISCOs occur at

$l$2

coincident with the event horizon and independent of the NUT charge $l$3. In non-extremal cases the ISCO lies outside the horizon and depends on both $l$4 and $l$5. In the Taub-NUT limit $l$6, the radius of the direct ISCO increases with NUT charge, while in massless Taub-NUT there is no stable circular orbit for timelike geodesics (Chakraborty, 2013).

Accretion-flow analyses recast these geodesic features into hydrodynamic structure. For relativistic, steady, axisymmetric, low angular momentum, inviscid, advective, geometrically thin accretion in Kerr-Taub-NUT spacetime, the equatorial effective potential is written as

$l$7

The parameter region admitting multiple critical points in the $l$8 plane is modified in opposite ways by $l$9 and $n=0$0: increasing $n=0$1 shrinks and shifts the region to higher $n=0$2 and lower $n=0$3, whereas increasing $n=0$4 shifts it to lower $n=0$5 and higher $n=0$6. The study concludes that the NUT parameter effectively mitigates the effect of black-hole rotation in deciding accretion-flow structure, and that the maximum disc luminosity $n=0$7 decreases with increasing $n=0$8 (Dihingia et al., 2020). This indicates that the NUT charge is not dynamically equivalent to ordinary rotation, even though both enter the metric through off-diagonal structure.

3. Frame dragging, Lense-Thirring anomalies, and axial dynamics

The exact Lense-Thirring precession in Kerr-Taub-NUT spacetime departs from the weak-field inverse-cube law in the strong-gravity region. Along the pole, the precession rate is reported to be maximal just outside the horizon, then to fall sharply to zero at a certain radius, then to rise again to a local maximum before eventually decaying according to the usual $n=0$9 law. This anomalous non-monotonic behavior is maximal near the pole, disappears beyond a critical angle toward the equator, and is absent on the equator (Chakraborty, 2014).

Along $a=0$0, the precession rate reduces to

$a=0$1

and vanishes at

$a=0$2

which lies outside the horizon for physically reasonable values (Chakraborty, 2014). In the Taub-NUT limit $a=0$3, the exact frame-dragging effect does not vanish for spinning gyroscopes: $a=0$4 and for $a=0$5,

$a=0$6

This directly shows that the NUT charge can source gyroscope precession even in the absence of conventional Kerr rotation (Chakraborty et al., 2013).

A complementary dynamical signature appears in the motion of unbound high-energy particles along the rotation axis. In Kerr-Taub-NUT spacetime, only particles satisfying

$a=0$7

can escape along the symmetry axis $a=0$8; in Kerr, the condition is $a=0$9. The proper axial acceleration on the axis is written as

$a=n=0$0

with a critical Carter constant $a=n=0$1 determining whether the acceleration is positive. The paper finds that, for small $a=n=0$2, the NUT charge expands the range of repelled trajectories on the axis, while away from the axis its effect may strengthen or weaken repulsion depending on the Carter constant, position, and velocity. Numerically, increasing $a=n=0$3 moves the jet trajectory farther from the axis, increases observable axial velocity $a=n=0$4, and modifies both proper and observable acceleration in a non-uniform way (Zhang et al., 2017).

4. Penrose process, particle collisions, and superradiant scattering

The Penrose process in Kerr-Taub-NUT spacetime has been analyzed from several angles. In the quasi-equatorial single-splitting analysis, the lower limit for the change in black-hole mass at the horizon is written as

$a=n=0$5

with the corresponding energy-extraction condition

$a=n=0$6

The same study reports that the maximum extracted energy increases with the dimensionless NUT parameter $a=n=0$7, even though the relative volume of the ergosphere decreases and the effective potential is shifted downward (Abdujabbarov et al., 2011).

A later repetitive-Penrose analysis adopts the irreducible-mass viewpoint. There the horizon and ergosphere radii are

$a=n=0$8

and the irreducible mass is written as

$a=n=0$9

so that

$$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$0

That study finds that extractable energy decreases as the gravitomagnetic charge $$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$1 increases, that $$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$2 increases monotonically across iterations, and that repetitive extraction terminates once the spin falls below the threshold required by the iterative constraints (Alloqulov et al., 26 Mar 2026). This suggests that the single-event maximum-extraction calculation and the iterative irreducible-mass accounting are probing different operational notions of extractability.

Near-horizon collisions provide a different energetic probe. For particles moving in the equatorial plane, the center-of-mass energy satisfies

$$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$3

In extremal Kerr-Taub-NUT spacetime, unlimited center-of-mass energy can be approached if $$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$4 lies in the range

$$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$5

with $$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$6, provided one particle reaches the horizon with critical angular momentum. This differs from Kerr and Kerr-Newman and shows explicitly that the center-of-mass energy depends on both the rotation parameter and the NUT charge (Liu et al., 2010).

Wave scattering extends the energetic picture to field theory. For massless fields of arbitrary spin $$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$7, a generalized Teukolsky master equation separates into radial and angular parts in Kerr-Taub-NUT spacetime, and the angular sector contains an effective spin

$$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$8

together with a shifted azimuthal number

$$ds^2 = -\frac{\Delta}{p^2}(dt - A d\phi)^2 + \frac{p^2}{\Delta} dr^2 + p^2 d\theta^2 + \frac{1}{p^2} \sin^2\theta (a dt - B d\phi)^2$$9

The superradiant threshold remains of the Kerr form,

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$0

but the horizon area and thermodynamic factors are NUT-deformed. Analytical low-energy observables were obtained for scalar, neutrino, electromagnetic, Rarita-Schwinger, and gravitational fields, with superradiance present for integer-spin bosonic fields and absent for half-integer-spin fermionic fields (Lee et al., 2023).

5. Separability, hidden symmetries, and wave equations

The Dirac equation in Kerr-Taub-NUT spacetime is separable in Boyer-Lindquist coordinates. With

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$1

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$2

the spinor ansatz

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$3

leads to separated radial and angular ordinary differential equations with separation constant $$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$4. Exact angular solutions are obtained in several special cases in terms of associated Legendre, hypergeometric, or confluent hypergeometric functions (Cebeci et al., 2012).

The radial system can be transformed into Schrödinger-type equations,

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$5

with

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$6

At large $$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$7,

$$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$8

and the NUT parameter first enters at order $$\Delta = r^2 - 2Mr + a^2 - n^2, \qquad p^2 = r^2 + (n + a \cos\theta)^2,$$9, described in that analysis as a dipole-type contribution. The corresponding potential plots show that increasing $$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$0 flattens the barrier structure exterior to the event horizon (Cebeci et al., 2012).

In the Kaluza-Klein generalization, obtained by uplifting Kerr-Taub-NUT to five dimensions, boosting in the fifth direction, and reducing back to four dimensions, the Hamilton-Jacobi equation for geodesics separates completely only for massless particles. The resulting hidden symmetry is therefore generated by a conformal Killing tensor rather than an ordinary Killing tensor. Using a conformally related “effective” metric $$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$1, the explicit conformal Killing tensor is constructed, showing that the extra dimension appears through companion electromagnetic and dilaton fields and modifies the separability structure in a precise way (Esmer, 2013).

6. Stability, extensions, and broader theoretical context

The overspinning problem provides a stringent test of black-hole stability. For Kerr-Taub-NUT spacetime, the horizon condition is

$$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$2

A first-order test-field analysis found that overspinning is possible only for a non-generic, finely tuned sub-collection of nearly extremal backgrounds and only for frequencies in a narrow interval above the superradiant threshold. A second-order analysis using the Sorce-Wald condition,

$$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$3

then yields

$$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$4

showing that the event horizon is preserved and that the final state is driven away from extremality (Düztaş, 14 Jul 2025).

Several extensions place Kerr-Taub-NUT geometry inside larger solution families. In Einstein-Maxwell theory, a magnetized Kerr-Taub-NUT solution can be generated by the Ernst magnetization or Harrison transformation, producing a Melvin-type background with external magnetic field parameter $$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$5. For the extremal magnetized geometry, the near-horizon region has $$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$6 isometry, and the Kerr/CFT computation reproduces the extremal entropy through the Cardy formula (Siahaan, 2021). In the low-energy limit of heterotic string theory, an accelerating Kerr-Taub-NUT spacetime generated through the Hassan-Sen transformation carries Maxwell, dilaton, and axion fields, has black-hole and acceleration horizons, unavoidable conical singularities, and closed timelike curves associated with the NUT parameter (Siahaan, 2024). A related Kerr-Sen-Taub-NUT solution similarly combines rotation, electric charge, and gravitomagnetic monopole moment; in that spacetime, the NUT parameter forbids equatorial circular geodesics and light rings (Siahaan, 2019).

Modern mathematical reformulations further enlarge the significance of the Kerr-Taub-NUT family. In amplitude language, the Kerr-Taub-NUT solution corresponds to a gravitational electric-magnetic duality rotation acting on Kerr, and at linearized order its Newman-Penrose scalar is

$$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$7

The associated three-point amplitude is

$$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$8

which packages duality and spin as elementary operations (Emond et al., 2020). In a distinct but related self-dual setting, the self-dual analog of Kerr has been shown via the Cartan-Karlhede algorithm to be locally equivalent to the self-dual Taub-NUT instanton in $$A = a \sin^2\theta - 2n\cos\theta, \qquad B = r^2 + a^2 + n^2.$$9 Kleinian signature (Desai et al., 2024). These developments place Kerr-Taub-NUT spacetime simultaneously within classical exact-solution theory, strong-gravity phenomenology, perturbation theory, holography, and amplitude-based descriptions.