Self-Dual Black Hole: LQG and Geometric Duality

• Self-dual black hole is a quantum-corrected solution inspired by loop quantum gravity that replaces the classical singularity with a regular core and connects two asymptotically flat regions via a radial inversion symmetry.
• Key methodologies involve polymerization of the Kantowski–Sachs sector and trigonometric substitution for canonical variables to derive effective metrics that reveal modified horizon structures and thermodynamic properties.
• Observational and perturbative analyses constrain the polymeric parameter P through shifts in orbital dynamics, quasinormal mode behavior, and gravitational lensing, highlighting its potential relevance in dark matter phenomenology.

A self-dual black hole is, in the loop-quantum-gravity-inspired literature, a static, spherically symmetric, quantum-corrected black-hole geometry whose metric is invariant in functional form under the radial inversion $r\to a_0/r$, with $a_0$ set by the minimum area gap of loop quantum gravity; the same expression is also used, in a distinct geometric sense, for Ricci-flat Euclidean or Kleinian metrics with self-dual curvature, notably self-dual Taub–NUT and related instantons (Yan et al., 2022, Adamo et al., 8 Jan 2026). In the LQG-inspired case, the classical singularity is replaced by a regular core and a second asymptotically flat region, while phenomenology is largely controlled by a dimensionless polymeric function $P$, with $a_0$ typically Planck-suppressed at astrophysical radii (Yan et al., 2022).

1. LQG-inspired self-dual geometry

In the effective mini-superspace construction based on polymerization of the Kantowski–Sachs sector, the canonical variables are replaced by trigonometric functions,

$$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$

and the resulting effective Hamiltonian yields a static quantum-corrected metric (Yan et al., 2022). In one common parameterization, the line element is

$$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$

with

$$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$

$h(r)=r^2+\frac{a_0^2}{r^2},$

and

$$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$

The same geometry is also written with $r_+=2m$, $a_0$0, $a_0$1, and $a_0$2 (Brown et al., 2010).

The polymeric correction is encoded by

$a_0$3

where $a_0$4 is the Barbero–Immirzi parameter and $a_0$5 is the polymeric parameter (Yan et al., 2022). The area scale is

$a_0$6

so $a_0$7 represents the minimum area gap of LQG (Yan et al., 2022).

Self-duality is realized by the radial inversion

$a_0$8

together with the corresponding map on the horizon parameters; after a suitable time rescaling, the metric retains the same functional form (Brown et al., 2010). The angular sector

$a_0$9

is explicitly invariant under $P$0, and this enforces a minimum areal radius, eliminating the classical $P$1 singularity and rendering the spacetime geodesically complete (Yan et al., 2022). In this sense, the LQG self-dual black hole is a regular, two-horizon geometry connecting the usual exterior to another asymptotically flat region (Brown et al., 2010).

2. Global structure, temperature, and entropy

The LQG self-dual black hole has an outer horizon at $P$2, an inner Cauchy horizon at $P$3, and a minimal two-sphere at $P$4 (Silva, 2012). The areal radius

$P$5

attains its minimum at $P$6, with $P$7, so the interior is regular and effectively wormhole-like (Santos et al., 2021). The event-horizon area is

$P$8

which differs from Schwarzschild by the $P$9-dependent term (Anacleto et al., 2020).

Using the Hamilton–Jacobi tunneling method, the Hawking temperature of the self-dual black hole is

$a_0$0

which reduces to the Schwarzschild value in the limit $a_0$1 and $a_0$2 (Silva, 2012). A distinctive feature is the small-mass behavior: for sub-Planckian masses, the temperature tends to zero as $a_0$3, rather than diverging as in Schwarzschild (Hossenfelder et al., 2012).

Entropy calculations in the tunneling formalism with generalized uncertainty principles recover the Bekenstein–Hawking term and add logarithmic and further subleading corrections (Anacleto et al., 2015). For linear and quadratic GUPs, the corrected entropy contains logarithmic terms and additional contributions of several other types, reflecting the interplay between the horizon scale and the self-dual area scale $a_0$4 (Anacleto et al., 2015).

These thermodynamic properties feed directly into evaporation models. For ultralight self-dual black holes, the relevant mass range is below the Planck mass and the temperature is small compared to the mass; in that regime, the effective emitting surface is a Planck-size “pinhole,” and the low-energy approximation for greybody factors is valid (Hossenfelder et al., 2012). This has led to the proposal that ultralight self-dual black holes may behave as long-lived dark matter candidates (Hossenfelder et al., 2012).

3. Geodesics and observational constraints

For equatorial timelike geodesics, the conserved quantities are

$a_0$5

and the radial motion obeys

$a_0$6

so the polymeric correction $a_0$7 shifts orbital frequencies and periapsis precession relative to Schwarzschild (Yan et al., 2022). The leading periapsis advance is

$a_0$8

which makes stellar dynamics a direct probe of the polymeric sector (Yan et al., 2022).

A direct comparison with three decades of astrometric positions, radial velocities, and orbital precession data for the S0-2 star orbiting Sgr A$a_0$9 found no significant evidence for the self-dual spacetime and yielded

$$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$0

at $$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$1 confidence level (Yan et al., 2022). With $$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$2, these correspond to

$$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$3

respectively (Yan et al., 2022).

A rotating self-dual black hole metric, obtained by the Newman–Janis algorithm, has also been constrained through high-frequency quasi-periodic oscillations in X-ray binaries. Within the relativistic precession model, the strongest bound came from GRO J1655-40,

$$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$4

which translates to $$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$5 for $$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$6 (Liu et al., 2023). This is a strong-field, rotating-system bound, complementary to solar-system and Galactic-center analyses.

A later study combining timelike and null geodesics, Mercury’s perihelion shift, and the orbit of the S2 star around Sgr A$$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$7 reported

$$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$8

while also finding that the photon critical impact parameter decreases with increasing $$b \rightarrow \frac{\sin(\delta_b\,b)}{\delta_b},\qquad c \rightarrow \frac{\sin(\delta_c\,c)}{\delta_c},$$9, implying a smaller shadow (Xamidov et al., 13 Sep 2025). Across these analyses, the data consistently favor small polymeric corrections.

4. Stability, quasinormal modes, and wave scattering

The existence of an inner horizon makes perturbative stability a central issue. For the original LQG self-dual metric, linear analysis of a massless scalar field shows that the Cauchy horizon is stable only if

$$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$0

equivalently

$$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$1

where $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$2 are the outer and inner surface gravities (Brown et al., 2010). A distinct “symmetric” self-dual metric, form-invariant under $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$3, satisfies

$$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$4

and in that case the Cauchy horizon is protected for any value of $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$5 (Brown et al., 2010).

In the perturbation spectrum, sector dependence is pronounced. For massless scalar perturbations, sixth-order WKB and time-domain calculations showed that as $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$6 increases, the real part of the quasinormal frequency initially increases and then decreases, while $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$7 decreases, implying slower damping (Santos et al., 2015). For a massive scalar field nonminimally coupled to gravity, the spectrum depends strongly on the field mass $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$8, the coupling $$ds^2=-\,f(r)\,dt^2+\frac{dr^2}{g(r)}+h(r)\left(d\theta^2+\sin^2\theta\,d\phi^2\right),$$9, and the LQG parameters; the self-dual black hole is stable under scalar perturbations for small parameters, and the nonzero Ricci scalar $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$0 breaks the degeneracy between $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$1 and $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$2 even in the massless case (Santos et al., 2021). By contrast, axial gravitational perturbations display increasing $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$3 and $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$4 as $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$5 grows, and no unstable axial modes were found in the parameter range studied (Yang et al., 2023).

Low-energy scattering and absorption also carry distinctive signatures. For a massless scalar field, the small-angle differential cross section becomes

$$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$6

and the low-frequency absorption cross section is approximately

$$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$7

with $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$8 the self-dual horizon area (Anacleto et al., 2020). Unlike Schwarzschild, the differential scattering and absorption cross sections remain nonzero in the zero-mass limit because the throat geometry set by $$f(r)=\frac{(r-r_+)(r-r_-)(r+r_*)^2}{r^4+a_0^2},\qquad g(r)=\frac{(r-r_+)(r-r_-)r^4}{(r+r_*)^2(r^4+a_0^2)},$$9 persists (Anacleto et al., 2020).

5. Particle dynamics, ISCO shifts, and accretion observables

When the self-dual black hole is immersed in an external asymptotically uniform magnetic field, both electrically charged particles and magnetic dipoles exhibit modified circular motion. In the neutral geodesic sector with $h(r)=r^2+\frac{a_0^2}{r^2},$0, increasing $h(r)=r^2+\frac{a_0^2}{r^2},$1 shrinks the ISCO radius and raises the orbital frequency; for example,

$h(r)=r^2+\frac{a_0^2}{r^2},$2

in the $h(r)=r^2+\frac{a_0^2}{r^2},$3 case reported for the neutral/dipole sector (Uktamov et al., 2024). For magnetic dipoles, the ISCO radius is greater than that of electrically charged particles because of the magnetic-field interaction (Uktamov et al., 2024).

The same study found that the quantum correction parameter shifts the electromagnetic flux and disk temperature profiles toward the central object, leading to a slight increase in these quantities (Uktamov et al., 2024). In the Novikov–Thorne framework, the inner edge is set by the ISCO, so the inward shift produced by increasing $h(r)=r^2+\frac{a_0^2}{r^2},$4 modifies the radiative efficiency and differential luminosity (Uktamov et al., 2024).

Thin-disk imaging calculations likewise show that increasing $h(r)=r^2+\frac{a_0^2}{r^2},$5 shrinks direct and secondary images and reduces the critical impact parameter. One representative sequence is

$h(r)=r^2+\frac{a_0^2}{r^2},$6

so the shadow angular radius decreases as $h(r)=r^2+\frac{a_0^2}{r^2},$7 increases (Xamidov et al., 13 Sep 2025). The same analysis found that increasing $h(r)=r^2+\frac{a_0^2}{r^2},$8 enhances the peak flux and shifts it inward, while the redshift pattern remains qualitatively Schwarzschild-like but somewhat weaker at larger $h(r)=r^2+\frac{a_0^2}{r^2},$9 (Xamidov et al., 13 Sep 2025). This suggests that disk spectroscopy, direct imaging, and timing probe different combinations of the geodesic and radiative structure.

6. Euclidean, Kleinian, and lower-dimensional uses of the term

The expression “self-dual black hole” is also used for a distinct class of vacuum metrics whose curvature is self-dual. In four dimensions, this means that one chiral Weyl spinor vanishes,

$$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$0

or equivalently $$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$1 in suitable conventions (Adamo et al., 24 Jul 2025). In this usage, the self-dual Taub–NUT metric in Gibbons–Hawking form,

$$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$2

is a canonical example (Adamo et al., 24 Jul 2025). Recent work has shown that such backgrounds admit exact quasi-momentum eigenstates and background-exact MHV graviton amplitudes, with the holomorphic collinear splitting functions equal to those of flat space (Adamo et al., 24 Jul 2025).

A dual-twistor reformulation goes further: asymptotically flat self-dual black holes, including self-dual Taub–NUT, Eguchi–Hanson, and the self-dual Plebański–Demiański metric, can be encoded by holomorphic quadrics in flat dual twistor space, from which their hyperkähler structure, Kerr–Schild form, Gibbons–Hawking form, Killing vectors, and Killing tensors are read off directly (Adamo et al., 8 Jan 2026). In a related Kleinian program, the self-dual Kerr–Taub–NUT solution with $$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$3 is self-similar to its near-horizon region, hidden conformal symmetries become exact, and the separated wave equation maps exactly to the hydrogen atom, with elementary polynomial wavefunctions and a hyperfine-splitting structure away from self-duality (Guevara et al., 2023).

Celestial-holographic studies of linearized rotating self-dual Kerr–Taub–NUT in $$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$4 signature construct coherent states on the celestial torus carrying an infinite tower of $$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$5 charges, interpreted as soft hair (Crawley et al., 2023). Near Kleinian horizons of the self-dual Schwarzschild–Taub–NUT solution, the geometry admits a local infinite-dimensional symmetry generated by supertranslations and superrotations, with integrable Noether charges (Giribet et al., 16 May 2025). A separate three-dimensional usage appears in topological massive gravity, where self-dual warped AdS$$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$6 black holes possess a chiral asymptotic symmetry algebra with one Virasoro copy and

$$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$7

and their entropy is reproduced by the Cardy formula (Chen et al., 2010, Li et al., 2010).

Taken together, these usages define a broad but coherent research area. In LQG-inspired black-hole phenomenology, the self-dual black hole is a regular polymer-corrected deformation of Schwarzschild or Kerr, controlled primarily by $$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$8 and $$r_+=\frac{2GM}{(1+P)^2},\qquad r_-=\frac{2GM\,P^2}{(1+P)^2},\qquad r_*=\frac{2GM\,P}{(1+P)^2}.$$9. In Euclidean and Kleinian gravity, it is a self-dual curvature background with hidden integrability, twistor structure, and exact amplitude technology. The shared terminology reflects a common emphasis on duality, regularity, and enhanced geometric structure, but the underlying constructions are not the same.