Integrable Singularities in Black Hole Interiors

• Integrable singularities are regions in black hole interiors where curvature invariants diverge mildly yet yield finite volume integrals and allow a C⁰ metric extension.
• They provide finite tidal forces along radial paths, enabling geodesic completeness in contrast to destructive Schwarzschild-type singularities.
• Key challenges include stability under non-radial perturbations, matter accumulation at the core, and reconciling these features with cosmic censorship.

An integrable singularity in a black hole interior is defined as a locus where curvature invariants such as the Ricci scalar or Kretschmann scalar diverge as the singular point is approached, but the divergence is sufficiently mild that their volume integrals remain finite and the metric potentials and their first derivatives are regular at the core. This structure interpolates between classical Schwarzschild-type singularities, where geodesics terminate and destructive tidal forces diverge, and fully regular 'core' geometries with de Sitter-like interiors shielded by an inner (Cauchy) horizon. The principal features, physical mechanisms, and theoretical challenges of integrable singularities are reviewed and critically examined, with particular emphasis on their role in black-hole interiors, the (non-)traversability of these loci, and their stability under perturbations (Arrechea et al., 24 Apr 2025).

1. Definition and Metric Characterization

The canonical setting for integrable singularities involves static, spherically symmetric line elements with a Misner–Sharp mass function $m(r)$ parameterizing the interior geometry: $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$ An integrable singularity is present when the central energy density $\epsilon(r)$, entering $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$, is not necessarily bounded but fulfills $\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$ and $m(r) = m_1\, r + \mathcal{O}(r^2)$ as $r \to 0$ with $m_1 > 1/2$. The function $f(r)$ thus stays finite at $r=0$ with $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$0.

Curvature invariants diverge in a controlled manner: $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$1 while $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$2 is locally integrable. This is the minimal departure from regularity that allows radial geodesic extension through the core without infinite tidal distortion (Arrechea et al., 24 Apr 2025, Estrada et al., 2023, Estrada, 3 Feb 2026, Lukash et al., 2012).

2. Physical Properties and Tidal Structure

Radially infalling observers experience Riemann components in their local frame given by: $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$3 Since $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$4 and $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$5 are finite, all orthonormal tidal tensor components remain finite at $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$6. The geodesic deviation equation admits no divergence in the radial separation vector, and thus Tipler and Ori's criteria for weak (deformational) singularities are satisfied along radial trajectories (Arrechea et al., 24 Apr 2025, Estrada et al., 2023, Lukash et al., 2011). For these trajectories, both the proper-time and affine-parameter extensions through $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$7 are regular.

However, non-radial geodesics cannot, in general, traverse $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$8: they terminate at the singular locus since the effective potential diverges unless angular momentum vanishes (Arrechea et al., 24 Apr 2025).

3. Backreaction, Matter Accumulation, and Stability

The stability of integrable singularities under linear perturbations and realistic physical processes is problematic. Scalar or higher-spin test fields on these backgrounds satisfy generalized Regge–Wheeler–Zerilli equations, with the effective potential near the core behaving as $$ds^{2} = -f(r) dt^{2} + f(r)^{-1} dr^{2} + r^{2}(d\theta^{2} + \sin^{2}\theta d\phi^{2}),\quad f(r) = 1 - \frac{2m(r)}{r}.$$9 for multipole $\epsilon(r)$0. For $\epsilon(r)$1, the field energy density is integrable, but for $\epsilon(r)$2 the stress tensor diverges as $\epsilon(r)$3 with $\epsilon(r)$4, making the integrated energy density at the core non-integrable. Nonradial perturbations therefore destabilize the integrable structure (Arrechea et al., 24 Apr 2025).

Matter with angular momentum, or any non-radial component, accumulates at $\epsilon(r)$5 due to focusing by diverging expansion scalars $\epsilon(r)$6. Without anisotropizing or dissipating this angular momentum, the integrability condition fails as $\epsilon(r)$7 deviates from the required $\epsilon(r)$8 scaling (Arrechea et al., 24 Apr 2025).

Extended bodies, governed by the Mathisson–Papapetrou–Dixon equations, are generically deflected onto nonradial trajectories by spin–curvature coupling and develop internal differential stresses that become destructive. Even with finite tidal forces along a central trajectory, deformations and energy accumulations preclude the survival of realistic objects through the singularity (Arrechea et al., 24 Apr 2025).

4. Examples, Model Realizations, and Dimensional Extensions

Constructed models span classical general relativity, Lovelock gravity in higher dimensions, fluids with screened energy profiles, and phenomenological "astrogenic" universes:

• In classical GR, mass profiles such as $\epsilon(r)$9 generate a single-horizon, nondegenerate black hole solution with Ricci scalar $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$0 as $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$1, whose four-volume integral remains finite. The near-core geometry is warped AdS (Ovalle, 2023).
• Lovelock theories (Einstein-Gauss-Bonnet in $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$2, cubic Lovelock in $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$3) admit vacuum solutions with integrable singularities, characterized by $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$4 and $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$5 for $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$6; the absence of inner (Cauchy) horizons eliminates mass inflation instabilities (Estrada, 2024, Estrada et al., 2023).
• Fluids of strings with screened energy densities, $m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$7, yield interior solutions matched to Schwarzschild (or Reissner–Nordström) exteriors. These provide integrable singularities free of internal horizons and satisfy both strong and dominant energy conditions in certain parameter regimes (Estrada, 3 Feb 2026).
• In dynamical black-hole simulations, singularities arise in the interior metric when mass-profile parameters approach a critical value; if the divergence is logarithmic (rather than power-law), the singularity is integrable, otherwise it is not (Ovalle, 31 Aug 2025, Casadio et al., 3 May 2026).

A comparative summary is provided below:

Model Class	Key Feature	Inner Horizon	Curvature Divergence	Integrability
Schwarzschild	Destructive singularity	Absent	$m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$8	Not integrable
Regular (Bardeen/Hayward)	Cauchy horizon, regular	Present	Finite	Not applicable
Integrable (classical GR)	C$m(r) = 4\pi \int_{0}^{r} \epsilon(x)x^{2}dx$9 extension, weak	Absent	$\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$0 ($\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$1)	Volume integrals finite; tidal forces finite
Lovelock (EGB, cubic)	Vacuum IS in $\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$2	Absent	$\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$3 ($\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$4)	Integrable for $\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$5
String-fluid	Screened energy interior	Absent	$\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$6	Total energy finite

5. Causal and Global Structure

The causal structure of an integrable singularity differs crucially from both Schwarzschild and regular black holes. With only a single (event) horizon and no inner Cauchy horizon, geodesics extend from asymptotic infinity through the horizon, across $\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$7, and into a time-reversed (white hole or cosmological) region. This global extension is continuous in the metric, though not in curvature invariants (Lukash et al., 2012, Lukash et al., 2013, Lukash et al., 2011). Penrose diagrams for these solutions replace the classical spacelike singularity with a lightlike, geodesically complete boundary, potentially permitting cosmological flows or "baby universes" generated from black hole interiors ("astrogenic universes").

Null singularities of the weak (integrable) type, realized in strong cosmic censorship scenarios for Reissner–Nordström and Kerr, also exhibit this behavior: while the metric extends continuously, the Christoffel symbols are not $\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$8-integrable, preventing further classical extension but not geodesic completeness (Luk, 2013).

6. Theoretical Challenges and Open Problems

While integrable singularities provide a C$\int_{0}^{r_h} \epsilon(r)\, r^{2}\,dr < \infty$9 extension beyond traditional black-hole interiors and circumvent inner horizon (mass inflation) and core instabilities, several substantial difficulties remain:

• Generic perturbations: Linear test fields with nonzero angular momentum (or higher-spin) generically lead to non-integrable stress tensors at $m(r) = m_1\, r + \mathcal{O}(r^2)$0 (Arrechea et al., 24 Apr 2025).
• Matter with angular momentum: Accumulation at the focal point ($m(r) = m_1\, r + \mathcal{O}(r^2)$1) leads to violations of the integrability requirement for the mass function (Arrechea et al., 24 Apr 2025).
• Extended objects: Spin–curvature and finite-size effects destroy extended bodies at or before crossing $m(r) = m_1\, r + \mathcal{O}(r^2)$2 even if tidal tensors are finite for idealized radial-pointlike probes (Arrechea et al., 24 Apr 2025).
• Cosmic censorship: Inclusion of $m(r) = m_1\, r + \mathcal{O}(r^2)$3 as a manifold point would violate strong cosmic censorship since the spacetime is $m(r) = m_1\, r + \mathcal{O}(r^2)$4 but not $m(r) = m_1\, r + \mathcal{O}(r^2)$5-extendible through the singularity.
• Quantum gravity completion: Classical integrable singularities may be stabilized or regularized by quantum gravitational effects (e.g., through smearing, probabilistic regularization, or dimensional reduction), but no fully consistent semiclassical or quantum implementation exists to date (Stoica, 2014, Casadio et al., 3 May 2026).

7. Broader Implications and Physical Significance

Integrable singularities offer an alternative endpoint for gravitational collapse in classical and semi-classical gravity, characterized by finite tidal forces along privileged directions and the avoidance of mass-inflation instabilities associated with internal horizons. These solutions have consequences beyond singularity avoidance—enabling geodesic completeness, potentially allowing black-hole spacetimes to seed new cosmological regions, and providing a controlled context for exploring the interface between classical gravity and quantum gravity phenomenology, including effects such as metric/measure reduction and possible dimensional reduction at Planckian curvatures (Lukash et al., 2012, Stoica, 2014). Nonetheless, their physical relevance is contingent on stabilizing the interior structure against perturbations, understanding matter and field evolution near and through the singular locus, and integrating quantum effects in a consistent formulation (Arrechea et al., 24 Apr 2025).