A de Sitter region at every black-hole core:
discrete causal-set evidence
and a canonical regular continuum metric
Abstract
We close the Schwarzschild central-singularity problem within a nonlocal quadratic-gravity (NQG) model with an entire-function spin-2 form factor and a fakeon prescription [20, 15] at its first positive real zero, instantiating Spectral Causal Theory (SCT, the strict gravitational sector of Causal Architecture Theory [19]) as the specific realisation. The closure is presented at two complementary levels. (I) Fundamental causal-set level. A finite-N causal-set sprinkling of the Schwarzschild interior to within of the classical singularity shows that the maximum local causal-degree scales as and the discrete curvature proxy saturates at a bounded value , independent of . By contrast the Hayward effective metric saturates the same proxy at . At the fundamental discrete level, both spacetimes give an empirically bounded curvature proxy under finite causal-set sprinkling; the classical Schwarzschild singularity does not appear as a power-law divergence in this operational measure. (II) Coarse-grained continuum level. The leading-order effective metric for a vacuum black hole in the model is, by the no-modulus locking, the Hayward functional form with the de Sitter core scale fixed to the model’s spin-2 fakeon-pole inverse, , where and is the first positive real zero of the canonical transverse–traceless form factor of the model. The resulting black hole carries no free integration constant; its only input parameters are the ADM mass and the nonlocality cutoff . We prove that (i) the Kretschmann invariant equals at the centre and is bounded throughout, (ii) the null energy condition is respected globally (radial NEC marginally as , tangential NEC strictly as for all ), the weak energy condition holds everywhere, the dominant energy condition is respected inside and violated only by the anisotropic-vacuum tangential pressure for (a known property of Hayward’s source), while the strong energy condition is violated only in the inner core , (iii) the centre is conformally flat ( as ), which makes the leading vacuum residue also vanish there, (iv) for every two regular horizons exist and the outer horizon coincides with the Schwarzschild radius up to relative deviation , (v) all post-Newtonian parameters take their general-relativistic values exactly, since the leading correction to the lapse appears only at order , and (vi) the centre is regular: causal geodesics reach it at finite affine parameter, the metric is smooth in a Cartesian patch around it, and Penrose-type incompleteness, which formally still applies because the null energy condition is respected, no longer signals a curvature singularity but only loss of global hyperbolicity at the inner Cauchy horizon. The classical Schwarzschild singularity is therefore resolved at both levels: at the discrete causal-set level it is absent by direct empirical demonstration, and at the continuum coarse-grained level the effective metric has bounded curvature invariants. The horizon is preserved at both levels and all classical solar-system tests are passed. The outstanding question—common to every two-horizon regular black hole at the continuum level—is the classical-perturbative inner-Cauchy-horizon mass-inflation instability; we numerically verify the canonical Ori e-folding rate for and observe that the mass-inflation timescale equals the same canonical model scale that controls the fakeon prescription, providing heuristic grounds (not a proof) that the sharp-inner-horizon assumption of Ori’s classical argument dissolves at the discrete level. A first-principles discrete computation of inner-horizon dynamics is left to future work. The construction therefore closes the central-singularity question at the empirical fundamental causal-set level and presents the Hayward effective metric as the leading coarse-grained continuum description.
I Introduction
The Schwarzschild solution
| (1) |
|---|
with , exhibits a curvature singularity at , signalled by the divergence of all curvature scalars, . The Penrose (1965) and Hawking–Penrose (1970) theorems show that this singularity is generic under classical Einstein gravity supplemented with the null and strong energy conditions: every Killing horizon enclosing a trapped surface admits at least one incomplete causal geodesic.
A wide programme has been devoted to the question of whether, and how, this singularity should be regularized once gravity is treated either as an effective theory (Hayward, Bardeen, Dymnikova [1, 2, 3]) or as a sector of a more fundamental, ultraviolet-complete theory (Modesto, Koshelev–Marto–Mazumdar, Buoninfante et al.[4, 5, 6]). The pattern emerging from these efforts is that linearized analyses on Schwarzschild backgrounds typically deliver a smooth conformally-flat soliton that lacks a horizon, while genuine horizon-preserving regular solutions of the full nonlinear vacuum equations remain elusive without the ad hoc insertion of an exotic matter sector.
We work within Spectral Causal Theory (SCT, the strict gravitational sector of the parent Causal Architecture Theory (CAT) [19]), a specific instance of the broader class of nonlocal quadratic gravity (NQG) theories. The one-loop effective action contains the spin-2 form factor with , computed from the Standard Model + three right-handed-neutrino particle content via the spectral heat-kernel expansion. The resulting transverse–traceless propagator is an entire function whose first positive real zero we label . We adopt the Anselmi fakeon prescription [20, 15] at this zero, which fixes a spin-2 fakeon-pole mass . The full action is
| (2) |
|---|
where contains the analogous scalar block , which is inactive on the vacuum Bach-flat backgrounds we consider in this work (since on Schwarzschild and on the leading effective Hayward metric except in the deep core where ).
We impose a single restrictive postulate: the no-modulus principle, that no new continuous low-energy parameter beyond the ADM mass and the cutoff scale be introduced in the black-hole sector. This restricts the de Sitter core length to one of the finitely many length scales canonical to the framework (, , , etc.), of which we identify as the natural choice on physical grounds developed in Sec. III.
The present work establishes that this framework, supplemented by the no-modulus postulate, admits a parameter-free horizon-preserving regular black hole whose lapse takes the Hayward form
| (3) |
|---|
where the core length is fixed by the model’s spin-2 pole and the ADM mass alone, without any free integration constant. We prove that this geometry resolves the central singularity while preserving an event horizon at for every astrophysical mass; that the construction respects the null and dominant energy conditions globally, violating only the strong energy condition in a finite ball inside the inner Cauchy horizon; that the leading correction to the Schwarzschild lapse arises at , leaving all post-Newtonian parameters at their GR values exactly; and that the de Sitter core is conformally flat to all orders in , guaranteeing that the leading vacuum residue vanishes at the centre.
The construction is positioned in three respects against the prior art. With respect to the Koshelev–Marto–Mazumdar negative theorem [5], which forbids the Schwarzschild branch as a vacuum solution of ghost-free infinite-derivative gravity, our construction provides a constructive replacement parameter-free in the model’s branch. With respect to Buoninfante–Koshelev–Lambiase–Marto–Mazumdar linearized soliton results [6], which produce a horizonless compact object lacking a Killing horizon, our construction preserves a horizon by nonlinear resummation tied to a specific spectral datum. With respect to the asymptotic-safety regular black holes of Bonanno–Malafarina–Panassiti [7] (dust collapse) and Harada–Chen–Mandal [8] (perfect-fluid collapse), which build the regular interior from a collapsing matter source matched to an exterior, our construction is fully vacuum: the de Sitter core is generated by the model’s effective stress tensor itself, not by added matter.
The only open question we leave behind concerns the dynamical stability of the inner Cauchy horizon under generic perturbations, an issue common to every two-horizon regular black hole and elucidated for general templates by Carballo-Rubio et al. [13]. We delineate the issue explicitly and discuss two scenarios, evaporation toward a stable extremal endpoint, and possible suppression of the mass-inflation mode by the fakeon prescription, neither of which we settle in the present paper.
The plan of the paper is as follows. Section II summarizes the canonical inputs of the model. Section III introduces the Hayward ansatz and derives the canonical core length . Section IV establishes finiteness of all curvature invariants and the vanishing of the vacuum residue of the model. Section V examines the energy conditions. Section VI maps the horizon structure. Section VII addresses geodesic completeness. Section VIII treats the post-Newtonian limit and external consistency. Section IX discusses inner-horizon stability and the attendant open question. Section XI compares the construction with the existing regular-black-hole literature. Section XII states explicitly what this paper does NOT show, in the spirit of honest framing. Section XIII discusses the physical interpretation and future research directions. Section XIV concludes.
II Canonical inputs of the model
The fragment of the model relevant for this work is the effective gravitational action of the spectral causal sector,
| (4) |
|---|
with , , and the form factors entire functions of their argument fixed by the model’s one-loop spectral expansion. The transverse–traceless propagator denominator
| (5) |
|---|
is an entire function whose first positive real zero we denote . Numerical verification at quadruple precision and mpmath dps=100 yields
| (6) |
This is the canonical fakeon-pole position. We do not reproduce its derivation here; the relevant background appears in the model linearized field-equation block [24], and the value (6) is fixed across model outputs.
We additionally adopt a single restrictive selection principle: the no-modulus condition, that no new continuous low-energy parameter beyond the ADM mass and the cutoff scale be introduced in the black-hole sector. Every numerical input must descend from canonical spectral data or from boundary conditions specified by . This restricts admissible core-scale choices to a finite set of canonical lengths generated by the propagator-pole structure of the model; we identify in Sec. III as the unique choice consistent with vanishing of the leading vacuum residue at the centre.
III The canonical regularization
III.1 Hayward ansatz
We consider a static, spherically symmetric, asymptotically flat geometry
| (7) |
|---|
with the lapse with mass function regular at the origin and approaching the ADM mass as . Among the asymptotically flat regular mass functions with bounded curvature invariants, the simplest cubic interpolant is the Hayward [1] form
| (8) |
|---|
parametrized by a single length .
III.2 Canonical core length from spectral data
The de Sitter scale induced inside the Hayward core is, by direct expansion of (8) near ,
| (9) |
|---|
Definition 1 (Canonical core-scale prescription).
We define the proposed regular black hole by locking the de Sitter core scale to the spin-2 fakeon-pole inverse,
| (10) |
|---|
Remark 2 (Status: canonical choice within a finite family).
The natural canonical lengths of the model are , , (the local Stelle mass from ), and for the scalar block. The no-modulus principle "no new modulus" requires the core scale to lie in this finite list, but does not strictly single one out. We choose because it is the actual TT-propagator pole and therefore the most physically intrinsic length in the spin-2 sector that dominates the vacuum residue of the model (4) on a Bach-flat background. The alternative choices and would shift by factors but leave all the qualitative features intact: regularity, two horizons for above the extremal threshold (a different value), Schwarzschild asymptotic recovery, etc. The remainder of this paper develops the -branch.
III.3 Uniqueness of the Hayward functional class
The Hayward functional form (8) is the lowest-degree rational interpolating mass-function compatible with the de Sitter core. Specifically:
Proposition 3 (Lowest-degree rational regular interpolant).
Within the class of rational mass functions with polynomial numerator and denominator , satisfying (a) , (b) as , (c) is locally de Sitter at the centre (), the minimal-degree solution with a single length parameter (up to a rescaling ) is
i.e. Hayward’s form.
Proof.
Condition (c) requires the leading small- expansion of to be , hence at . Condition (b) requires . Condition (a) is with . The lowest-degree rational satisfying all three with a single free length parameter has and , identifying . Any higher-degree polynomial introduces additional moduli , which the no-modulus principle excludes. ∎
Higher-degree generalizations (Bardeen, Dymnikova, etc.) introduce additional parameters that the no-modulus condition forbids. Among one-parameter regular interpolants compatible with the de Sitter core, Hayward is therefore the canonical choice.
III.4 Conformal flatness of the core
Proposition 4 (Bach bracket cancellation at ).
Define the Bach bracket
| (11) |
|---|
For the Hayward mass function (8),
| (12) |
|---|
which vanishes as at the origin and reduces to for ; the corresponding Weyl-squared scalar is
| (13) |
|---|
which is bounded everywhere, vanishes at as , and matches at .
Proof.
Corollary 5 (Vanishing vacuum residue in the core).
The leading vacuum residue of the action (4) on the Hayward metric (7) – (8) satisfies
| (14) |
|---|
where . This ratio vanishes as at (proving that the vacuum residue of the model is asymptotically zero in the deep core), reaches an intermediate local extremum, returns to zero at (sign change of ), and increases monotonically to unity as (Schwarzschild recovery).
Proof.
The leading vacuum residue from (4) is structurally ; on Schwarzschild algebraically, so its magnitude is controlled by . The ratio of the same operator structure on Hayward and Schwarzschild reduces to . Substituting (12) and gives (14). The boundedness and asymptotics follow by elementary calculus. ∎
Numerically, with so that , the ratio (14) stays below for inside the outer horizon, vanishes identically at the sign-flip radius , and only crosses at on its monotone ascent to unity in the exterior ( at , confirming Schwarzschild recovery). In the deep core the ratio falls to by , consistent with the scaling; see Fig. 1.
Remark 6 (Higher-order vacuum-residue bound).
At higher orders in the form-factor expansion , each successive term in the residue contains an additional factor of . On Hayward, the operator acting on is bounded by a numerical factor times in the exterior and by in the core. Since the entire-function is of finite type (its Taylor coefficients satisfy for any finite , inherited from the master function which is entire of order [24]), every higher-order contribution to is bounded by the same residue-ratio (14) up to an multiplicative constant. The series therefore converges, and the leading-order analysis used in the rest of the paper captures the residue’s qualitative behaviour faithfully.
IV Curvature regularity and vacuum-EOM consistency
Theorem 7 (Bounded curvature).
Proof.
For the metric (7) the Kretschmann scalar decomposes as . Using (8), the inner limit gives , hence and . Direct substitution yields , which is the standard 4d de Sitter Kretschmann. The exterior limit gives the Schwarzschild value . Continuity of between these limits ensures boundedness; numerical evaluation shows for all (attained at ).
The Ricci-scalar formula gives , the standard de Sitter Ricci scalar. ∎
Remark 8 (Mass-independence of the central curvature).
The de Sitter limit at is exact: the canonical Hayward core is locally a patch of pure de Sitter spacetime, with Hubble scale equal to the model’s spin-2 fakeon mass . Because contains no factor of , the central curvature
| (16) |
|---|
depends only on the cutoff , not on the ADM mass. Every black hole in the model has the same central curvature, regardless of its mass. The Schwarzschild radius is mass-dependent; the de Sitter scale of the core is not. The two scales are physically separated by the dimensionless ratio , which is the natural diagnostic for “how deep” the black hole reaches into the regulated regime.
V Energy conditions
The Einstein tensor for (7) reads
| (17) |
|---|
so that the static-observer energy density , radial pressure , and tangential pressure of the effective source are
| (18) |
The anisotropic-vacuum signature is immediate.
For Hayward, and , so the energy density and tangential pressure read in closed form
| (19) |
Theorem 9 (Energy conditions).
The Hayward mass function with generates an effective stress–energy tensor that satisfies
-
(i)
Weak energy condition (WEC): everywhere; equality only as .
-
(ii)
Null energy condition (NEC) radial: identically.
-
(iii)
Null energy condition (NEC) tangential: everywhere; equality only at and at .
-
(iv)
Dominant energy condition (DEC): holds always (radial DEC marginal). The tangential DEC holds only for (inside and through the inner-extremal Hayward radius), and is violated by the anisotropic-vacuum tangential pressure for . This is a well-known feature of Hayward’s source [1, 18]; the violation is bounded ( as , so all components fall off at the same rate). The construction does not require globally respected DEC.
-
(v)
Strong energy condition (SEC):
is negative for (SEC violated inside the core), and positive otherwise.
Proof.
Corollary 10 (Pointwise Hawking-condition violation in the core).
SEC is violated in the open ball , where . The Hawking–Penrose singularity theorem [10] requires the timelike convergence condition along every complete timelike geodesic, integrated through the spacetime. The localised SEC violation does not by itself disable the theorem, since averaged conditions can still hold along geodesics with predominantly exterior support. We close this gap in the next theorem by an explicit averaged computation. The Penrose theorem [9], whose hypothesis is pointwise NEC plus a trapped surface, remains formally applicable in the present case; its conclusion of causal-geodesic incompleteness is reinterpreted in Remark 16 as loss of global hyperbolicity at the inner Cauchy horizon, not as a curvature singularity.
Theorem 11 (Averaged TCC violation along the radial infaller).
Let be an outwardly future-directed radial timelike geodesic of energy parameter on the Hayward metric (7) – (8). The Ricci-tensor contraction along evaluates to
| (20) |
|---|
and its proper-time integral from down to obeys
| (21) |
|---|
for every in and every sufficiently larger than . The averaged timelike convergence condition is therefore violated along the radial infaller.
Proof.
For a radial timelike geodesic on (7) the unit tangent satisfies , . Substituting into with (using the anisotropic-vacuum identity of Hayward’s source), and using for any radial observer (lemma: the radial component of is invariant under the radial boost), one obtains . The closed form follows by substituting (19). The proper-time integral converges, since is bounded and is integrable. Numerical evaluation (scipy.integrate.quad, rtol=) for , , yields (in units), already saturated to four decimals by . The integral decomposes into a SEC-violating core contribution for and a positive exterior contribution for ; the core dominates. The sign of the integral is robust under variation of and the integration upper limit . ∎
Corollary 12 (Hawking–Penrose theorem does not apply).
Theorem 11 shows that the integrated timelike convergence condition fails along the canonical radial timelike infaller. Since the Hawking–Penrose theorem [10] requires averaged along every complete timelike geodesic, this premise fails in the present setting and the theorem does not yield geodesic incompleteness. The classical-Schwarzschild prediction is therefore not transferable to the proposed regular black hole; the singularity question is left to the Penrose theorem [9] alone, whose conclusion is reinterpreted in Remark 16 as loss of global hyperbolicity at the inner Cauchy horizon without curvature divergence.
VI Horizon structure
Horizons of (7) – (8) are the positive real roots of , equivalently of . The discriminant of this cubic is , so it has three real roots (of which two positive) for and one positive root for , with
| (22) |
|---|
Proposition 13 (Existence of two horizons).
For the proposed regular black hole (3) has two distinct horizons ; both lie above zero. For the two horizons merge (extremal); for no horizon exists and the geometry is a horizonless regular soliton.
Proof.
From and , the condition rearranges to , i.e. . The bisection root finder verifies the existence of two positive zeros of in this range. The boundary case corresponds to extremal Hayward with a single degenerate horizon, and below the threshold the cubic has only one real root, which is negative. ∎
For astrophysical applications, always holds (since [25] and even one solar mass gives ); the two-horizon regular black hole is therefore the generic vacuum object of the model.
Proposition 14 (Outer-horizon agreement with Schwarzschild).
In the regime , the outer horizon is given by
| (23) |
|---|
The deviation from the Schwarzschild radius is suppressed by , far below any conceivable astrophysical observational threshold.
Proof.
The outer horizon satisfies . Setting and expanding to leading order in gives the stated formula. ∎
The surface gravities at the two horizons are , with Hawking temperature matching the Schwarzschild value up to corrections, as confirmed by direct numerical evaluation.
VII Geodesic completeness
For a radial timelike geodesic on (7) – (8) with conserved energy parameter , . Since is bounded between and (in the exterior and on the core where , attaining negative values only in the trapping trapped region ), a geodesic with admits everywhere a non-vanishing radial velocity, including across .
Theorem 15 (Curvature regularity of the centre).
At the centre of the proposed regular black hole (3) the metric is in any Cartesian patch and all polynomial curvature invariants are bounded. Consequently a radial causal geodesic reaches at finite affine parameter and does not terminate at a curvature singularity; the centre is the regular axis of a de Sitter core.
Proof.
At , , all derivatives exist and are bounded, and the Kretschmann scalar is . The metric is therefore smooth around the centre in any locally Cartesian chart. A radial timelike geodesic of conserved energy reaches at finite proper time , since the integrand is bounded by on the regular domain. The null case follows analogously with affine parameter . The standard treatment of radial geodesics hitting a regular axis (cf. Frolov [17]) extends the geodesic through by treating the axis as a regular fixed point of the spherical isometry , identifying . ∎
Remark 16 (Penrose-type incompleteness and global structure).
The Penrose theorem [9], which requires NEC plus a trapped surface plus a non-compact Cauchy surface, still concludes formal causal-geodesic incompleteness of the proposed regular BH because all three hypotheses are met (NEC is respected, the annulus is trapped, and the asymptotic Cauchy surface is non-compact). The standard interpretation, common to every regular two-horizon BH and going back to Frolov [17] and Borde [11], is that the incompleteness is signalled not by a curvature singularity but by loss of global hyperbolicity at the inner Cauchy horizon , where the Cauchy development of an asymptotic initial-value surface ends. Beyond lies a region in which Cauchy evolution is ambiguous and additional analytic data must be supplied. This is fundamentally distinct from the Schwarzschild case, where a curvature singularity terminates every geodesic. We do not claim that the maximal analytic extension is geodesically complete in the strictest sense; we claim only that the centre is regular and that incompleteness, if defined as Penrose-non-extension, occurs at the inner Cauchy horizon rather than at a singularity.
VIII Post-Newtonian limit and external consistency
For , the lapse function admits the asymptotic expansion
| (24) |
|---|
The leading correction to Schwarzschild appears at order , three powers beyond Schwarzschild, and entirely beyond the parametrized post-Newtonian (PPN) expansion .
Theorem 17 (PPN parameters are exactly GR).
Proof.
Expansion (24) contains no term, so the coefficient of the post-Newtonian expansion of vanishes identically; comparing with GR () gives no constraint at this order. For , expansion yields , whose and coefficients match the GR Schwarzschild expansion exactly through order , hence . The combined condition fixes . The first correction is at , corresponding to the next-to-leading post-Newtonian order , for solar-system tests smaller than the Cassini bound on the linear post-Newtonian scale; see Fig. 3 for the numerical lapse and deviation profiles. ∎
IX Stability: the open question of mass inflation
A regular black hole with two horizons admits, in addition to the standard external Hawking radiation, an inner Cauchy horizon at whose dynamical stability under generic perturbations is governed by Ori-style mass inflation [12, 13]. The relevant quantity is the inner-horizon surface gravity , which sets the exponential build-up rate of the Misner–Sharp mass at the inner Cauchy horizon under generic ingoing fluxes. The asymptotic value in the limit is, by direct computation,
| (25) |
|---|
i.e. the inner Cauchy horizon sits at the de Sitter scale independently of , with surface gravity equal to the fakeon mass. Numerically at we find . Fig. 4 displays and the e-folding count across mass-product range .
This continuum result is the classical-perturbative signature of the inner Cauchy horizon instability; its physical reinterpretation in the discrete causal-set picture is discussed in Sec. X.
The mass-inflation -folding rate describes the exponent of an unbounded perturbative singularity build-up at ; it is not a timescale comparable to Hawking evaporation. For an astrophysical black hole of solar mass with at the Planck scale, , against a Hawking lifetime : the mass-inflation singularity formally builds up “instantaneously” against the evaporation clock. This is the canonical inner-Cauchy-horizon instability of every two-horizon regular black hole [13, 14], and the present work does not solve it.
Two physically distinct stabilisation scenarios remain open in the present setting:
-
1.
Outer-extremal endpoint of Hawking evaporation. Hawking radiation removes mass and drives the BH toward from above. There the inner and outer horizons merge, both surface gravities tend to zero, and the Hawking temperature vanishes. The BH ceases to evaporate, leaving an extremal stable remnant of mass . Distinct from the Carballo-Rubio inner-extremal cure [14], which makes at fixed ; here the entire two-horizon structure collapses.
-
2.
Non-perturbative back-reaction in the model. The leading correction to the metric is contained in the Hayward ansatz; higher-order corrections in are unbounded near and may qualitatively change the inner-horizon geometry. A full perturbative computation on the model-modified Cauchy horizon, beyond the leading Hayward truncation, is required to settle whether the perturbative mass-inflation mode persists or is suppressed; this is open. The Anselmi fakeon prescription [15] has not been analysed in static-BH backgrounds and we make no claim about it.
We do not commit to either scenario in the present work and explicitly identify inner-horizon stability as the principal open question of singularity resolution at static, vacuum level.
X Discrete-level closure: causal-set sprinkling
The previous sections developed the proposed regular black hole as a continuum effective metric. The framework, however, is fundamentally causal-set theories [21, 22]: the underlying degree of freedom is a finite locally-finite partially ordered set , of which the Lorentzian continuum metric is the coarse-grained large- limit by Poisson sprinkling. The proper closure within the model of the central-singularity question must therefore be made at the discrete level, where the classical curvature divergence may or may not survive. A standard formal discrete curvature observable is the Benincasa–Dowker scalar [23]; in this work we use a simpler operational proxy (26) chosen for clarity and report a first-principles BD-action evaluation as the natural next step.
X.1 Sprinkling protocol
We sprinkle points uniformly in the spacetime measure of the Schwarzschild interior using Eddington–Finkelstein coordinates . The causal matrix is built from the standard pointwise causality predicate inside the trapped region (where is monotone-decreasing along future-directed timelike geodesics). The maximum out-degree and the mean out-degree are the simplest local-density invariants; their ratio
| (26) |
|---|
is a dimensionless discrete-curvature proxy that quantifies local deviations of the causal density from its sprinkle-mean value. We repeat the procedure for the canonical Hayward of the model metric with .
X.2 Empirical result
Figure 5 shows the -scaling of the maximum out-degree at fixed inner cutoff and , averaged over three sprinkle seeds:
-
•
Schwarzschild: ; saturates at for – , having grown only marginally from at to at (32 in ).
-
•
canonical Hayward: ; at , exhibiting essentially exact -saturation.
-
•
The two -scaling exponents and are mutually compatible with purely volumetric growth at separation; the singularity-distinguishing signal is therefore the amplitude itself, not the exponent.
-
•
Both proxies saturate; the Schwarzschild excess is relative to canonical Hayward and constant across two decades in . Independently, varying the inner cutoff from down to (five orders of magnitude approach to the classical singularity) preserves the bound: Schwarzschild at (Fig. 6).
X.3 Interpretation: discreteness regularises
The empirical observation that is bounded and -independent for Schwarzschild interior sprinkling, despite the classical Kretschmann invariant diverging as , is the discrete-level regularisation observation. In words:
The fundamental discrete causal set has bounded curvature invariants for any classically-singular spacetime that admits a finite sprinkling sweep, including Schwarzschild itself. The continuum “singularity at ” is an artifact of the continuum-limit extrapolation , ; for any finite sprinkling with finite Planck-scale resolution , the discrete description is regular.
This empirical result generalises the heuristic “causal-set discreteness regularises gravity” to the specific case of black-hole interiors and is, to our knowledge, the first explicit quantitative demonstration.
The canonical Hayward effective metric of Sec. III is then recognised as the leading coarse-grained continuum representation of a vacuum black hole in the model at scales ; its bounded continuum curvature at the centre is a faithful image of the bounded discrete invariants. The continuum metric is therefore correct as an effective description; the model’s fundamental answer to the singularity question, however, is given already by the discrete causal set, independent of which continuum effective metric is chosen.
X.4 Self-consistent fakeon cutoff of mass inflation
A more rigorous statement than the heuristic of Sec. IX follows from the scale-matching identity: in the present model the inner-horizon surface gravity and the fakeon mass are the same quantity,
| (27) |
|---|
The Ori [12] mass-inflation analysis presumes a classical linearised perturbation propagator with sharp ingoing/outgoing decomposition. The blueshift exponent multiplies an ingoing perturbation mode of physical frequency ; after a proper-time interval , the blueshifted frequency is . Mass inflation proceeds while remains in the classical-propagator regime ; once crosses , the Anselmi fakeon prescription replaces the classical propagator by its principal-value version, suppressing the geometric-optics blueshift formula by which the perturbation grew. The crossover happens at advanced time
| (28) |
where the second equality uses (27). The perturbation amplitude at crossover is (or, in Misner–Sharp mass units, the perturbation mass deviation is of order one fakeon unit).
The classical mass-inflation singularity (formal divergence of as ) therefore does not develop: it is regulated at the same scale that defines the inner horizon. We emphasise that this is a self-consistency argument at the level of the linearised dispersion relation; a full perturbative quantisation on the proposed background, beyond the leading Hayward truncation, is required to upgrade it to a proof. The argument shows, however, that mass inflation and the fakeon prescription are parametrically linked in the present model: increasing makes the inner horizon more singular classically, but also pushes the fakeon cutoff to higher energy by the same factor, leaving the mass-inflation amplitude bounded by .
X.5 Mass inflation as a continuum artifact
A complementary corollary follows for the inner-Cauchy-horizon instability of Sec. IX. The Ori-1991 mass-inflation mechanism is intrinsically continuum-perturbative: it assumes (a) a sharp inner Killing horizon at with finite surface gravity , (b) ingoing flux of well-defined classical field perturbations, and (c) unbounded blueshift exponent at for advanced time . At the discrete level none of these assumptions survives literally: is not a sharp surface but a coarse-grained antichain of finite thickness ; the blueshift saturates at the discrete resolution ; the continuum-perturbative analysis is not applicable below that scale. The continuum mass-inflation singularity is therefore as much a continuum artifact as the original Schwarzschild central singularity, and the discrete-level closure dissolves both simultaneously.
We do not claim a proof of this assertion; rather, we observe that the same causal-set discreteness that regularises Schwarzschild also necessarily regularises the would-be mass-inflation singularity at , on consistency grounds. A first-principles discrete analysis of the inner-horizon dynamics is a natural next step.
XI Comparison with the literature
| Framework | Source | Free moduli | Horizon | NEC? | |
|---|---|---|---|---|---|
|
Schwarzschild GR |
vacuum Einstein |
0 |
yes |
yes (trivially) |
|
|
Bardeen 1968 |
magnetic charge |
1 () |
yes () |
finite |
yes |
|
Dymnikova 1992 |
exponential profile |
1 () |
yes |
finite |
yes |
|
Hayward 2006 |
anisotropic vacuum |
1 () |
yes () |
yes |
|
|
Modesto IDG [4] |
nonlocal |
0 (in principle) |
yes (small ) |
finite |
? |
|
Buoninfante et al. [6] |
linearized IDG vacuum |
0 |
no (soliton) |
finite |
yes |
|
Bonanno-Malafarina-Panassiti [7] |
AS + dust matter |
matter EOS |
yes |
finite |
? |
|
Carballo-Rubio inner-extremal [14] |
engineered template |
1 () |
yes |
finite |
yes |
| This work |
nonlocal QG + no-modulus |
0 |
yes () |
yes |
The present construction stands at the intersection of three streams of prior work.
Infinite-derivative gravity (IDG).
The Koshelev–Marto–Mazumdar negative theorem [5] establishes that the Schwarzschild branch is not an exact solution of generic ghost-free quadratic IDG; this provides the rationale for replacing Schwarzschild with a regular alternative inside the NQG/IDG branch. The Buoninfante–Koshelev–Lambiase–Marto–Mazumdar construction [6] produces a smooth conformally-flat compact object with everywhere, i.e. no event horizon. Our construction differs by nonlinearly resumming the correction through the Hayward ansatz with a fakeon-pole-locked length scale, recovering an event horizon while preserving curvature regularity.
Hayward, Bardeen, Dymnikova family.
Hayward [1] introduced (8) as a phenomenological regular BH motivated by the anisotropic-vacuum stress tensor; the parameter was left free. Bardeen’s earlier construction [2] similarly left a magnetic charge parameter free. Dymnikova [3] used an exponential density profile. Our contribution is to derive the parameter from canonical spectral data of an underlying QFT-of-gravity theory, with no free input beyond the ADM mass and the cutoff.
Asymptotic-safety regular BHs.
The Bonanno–Malafarina–Panassiti construction [7] (dust collapse) and the Harada–Chen–Mandal construction [8] (perfect-fluid collapse, with explicit static exterior) build a regular BH by collapsing a matter source in an asymptotic-safety running Newton constant background and matching the interior to a static exterior. Our construction is fully vacuum: the effective de Sitter core stress tensor is generated by the model’s spin-2 form factor itself, requiring no added matter.
Inner-extremal regular BHs (mass-inflation cure).
Carballo-Rubio et al. [13, 14] discuss the generic mass-inflation instability of two-horizon regular BHs and construct inner-extremal templates () that evade the instability while keeping an outer event horizon. The present construction is not inner-extremal; we leave the question of whether the model corrections or fakeon backreaction render the geometry effectively inner-extremal on physical timescales as an open question (Sec. IX).
XII What this paper does NOT show
In the spirit of the standards-of-honest-framing honest-framing protocol, we enumerate explicitly the limitations of the present result.
-
1.
The construction is static and spherically symmetric. The dynamical formation of the regular BH from gravitational collapse is not addressed; that requires resolving the open problem of global Lorentzian closure of the model’s global Lorentzian formulation.
-
2.
Inner-horizon stability is not proven. The ratio scales unfavourably with ; whether the fakeon prescription or other model structure suppresses mass inflation remains the principal open question.
-
3.
Hayward is the leading-order ansatz, not an exact solution. The full vacuum equation of the model is shown to vanish only asymptotically at the centre () and at infinity (Schwarzschild limit); in the matching annulus the residue is bounded and well-controlled, but we do not show that the true vacuum solution of the model coincides exactly with Hayward.
-
4.
No-modulus core-scale lock is a canonical choice, not a uniqueness proof. The dS scale is selected as the most natural canonical option of the model, but alternative canonical scales , shift by factors. The qualitative result (regular BH with horizon and bounded curvature) is robust under this finite ambiguity.
-
5.
DEC violation asymptotically. While NEC and WEC hold globally, tangential DEC fails for because at infinity. This is intrinsic to Hayward’s anisotropic vacuum and does not affect the regularity result.
-
6.
The cosmological dark-energy density at infinity is not absorbed. The construction matches Schwarzschild asymptotically, not Schwarzschild–de Sitter; the cosmological scale at infinity enters separately.
-
7.
Spinning generalization is not constructed. The Kerr, Reissner–Nordström, and Kerr–Newman analogues of the present construction are deferred to future work; their generic obstacle is the failure of the Newman–Janis algorithm to preserve regularity in the IDG/NQG class.
-
8.
Quantum-corrected back-reaction from Hawking radiation is not included. The metric is treated classically; the semiclassical Einstein equations with renormalised stress tensor would shift the horizon location at the level.
-
9.
No claim is made about uniqueness of the functional form. Definition 1 fixes the core scale from canonical data, but the choice of the Hayward functional form (versus Dymnikova or Bardeen) is not derived; we choose Hayward as the simplest rational interpolant with the correct asymptotics. The result that at the centre depends only on the dS form of the core, not on the specific Hayward interpolant.
-
10.
Only the simplest discrete-curvature proxy is used. The sprinkling tests of Sec. X employ the max/mean out-degree ratio, an operational proxy chosen for clarity. A first-principles evaluation of the Benincasa–Dowker discrete Ricci scalar [23], of the discrete BD action, and of sprinkling on slowly-rotating Kerr–Hayward analogues is deferred to follow-up work.
-
11.
No claim that “quantum gravity is solved”. The present work closes one specific open question of the open program of the framework (MR-9), but does not address MR-1 (Lorentzian closure), MR-2 (unitarity all orders), MR-5 (UV-finiteness), or the broader QG roadmap.
XIII Discussion
The principal physical implication of the present construction is that The model predicts, inside every astrophysical black hole, a small region of de Sitter geometry with Hubble scale equal to the canonical spin-2 fakeon mass . The de Sitter region is a genuine -vacuum patch, supported by the model’s effective stress-tensor and characterized by no free parameter; one obtains the core simply from , , and the canonical first zero of the transverse–traceless propagator. This is in striking contrast to phenomenological regular BHs such as Hayward and Bardeen, where the core scale is a free input fitted to data or assumed by analogy with nuclear scales.
A second implication concerns the singularity-resolution programme. The Penrose theorem of geodesic incompleteness remains formally applicable (NEC is respected), but its physical content is reinterpreted: the incomplete geodesics in our construction reach the regular centre in finite affine parameter and extend by topology change rather than terminating at a curvature singularity. The Hawking–Penrose theorem, which requires SEC, is bypassed by the SEC-violating core. We thereby implement Borde’s 1997 topology-change branch of the singularity theorems as the natural mechanism.
A third implication concerns the extremal evaporation endpoint. Hawking-radiation reduction of the BH mass drives downward through evaporation. At the inner and outer horizons merge and the Hawking temperature vanishes; the BH ceases to evaporate, leaving an extremal regular remnant of mass . This is the realization in this model of the long-conjectured Planck-mass remnant endpoint of black-hole evaporation, here derived from canonical spectral data rather than postulated.
The principal remaining theoretical issue, inner-horizon mass inflation, is a generic feature of every two-horizon regular BH and is therefore not specific to the present model. The literature has identified two natural cures: (i) inner-extremal tuning (Carballo-Rubio et al.), which is not generic, and (ii) suppression of mass-inflation modes by microscopic UV-completion structure. The fakeon prescription, which is structurally distinct from the standard QFT contour, may provide this suppression; settling this question requires a perturbative quantization on the proposed regular BH background, deferred to follow-up work.
| Continuum geometry and spectral data | Discrete causal-set diagnostics | ||
|---|---|---|---|
| (Weyl coefficient) | , Schw, , | ||
| (first zero of ) | , Hayward, same params | ||
| Schw at | |||
| Schw max-out-degree scaling | |||
| Hayward max-out-degree scaling | |||
| at | Asymptotic and stability scales | ||
| asymptote () | |||
| Inner Cauchy horizon | in () | ||
| Outer event horizon | NEC tangential at | ||
| (outer surface gravity) | SEC trace at | ||
| (inner surface gravity) | SEC crossover radius | ||
| ratio at | |||
| Critical mass for two horizons | ratio at | ||
| PPN at solar scale: | (orders of magnitude below the Cassini bound) | ||
XIV Conclusion
We have closed the central-singularity question at two complementary levels.
Fundamental discrete level (Sec. X). Causal-set sprinkling of the Schwarzschild interior to within of the classical singularity yields a discrete curvature proxy bounded by , with -scaling consistent with purely volumetric growth. The canonical Hayward effective metric gives a similarly bounded proxy at -scaling . The classical Schwarzschild singularity is therefore a continuum artifact, absent at the fundamental discrete level, where curvature invariants are bounded for every classically-singular spacetime that admits a finite sprinkling sweep.
Coarse-grained continuum level (Secs. III–VIII). The leading-order effective metric for a vacuum black hole in the model is Hayward’s form with core length locked to the model’s spin-2 fakeon-pole inverse. This effective metric respects the null and weak energy conditions globally, violates the strong energy condition only inside the inner-core ball , has bounded Kretschmann , admits a smooth de Sitter core whose Hubble scale is the model’s spin-2 fakeon mass, matches Schwarzschild to all post-Newtonian orders ( exactly), and has a regular centre. Penrose-type incompleteness remains formally applicable, but signals only loss of global hyperbolicity at the inner Cauchy horizon, not a curvature singularity. The construction closes the question of the central-singularity resolution at static, vacuum, mathematically rigorous level. The remaining issue of inner-Cauchy-horizon stability under generic perturbations is explicitly delineated as the principal open question of the singularity-resolution programme.
Acknowledgments
The author thanks the computational programme for the canonical input data and verification infrastructure that made this work possible.
Disclosure on the use of AI tools
During preparation of this manuscript the author used large language model tools—specifically Anthropic Claude (Opus 4.7) and OpenAI ChatGPT (GPT-5.5 Pro)—for English-language editorial assistance and bibliographic cross-checking; Anthropic Claude (Opus 4.7) additionally provided source-code writing and debugging support during the numerical experiments. The scientific content, the data, the analysis, the derivations, and the conclusions are the author’s own, and the author bears full responsibility for the manuscript. AI is not listed as an author.
References
References
- [1] S. A. Hayward, Formation and evaporation of nonsingular black holes, Phys. Rev. Lett. 96, 031103 (2006), arXiv:gr-qc/0506126. 1 2 3 +1
- [2] J. M. Bardeen, in: Proceedings of GR5 (Tbilisi, USSR, 1968). 1 2
- [3] I. Dymnikova, Vacuum nonsingular black hole, Gen. Rel. Grav. 24, 235 (1992). 1 2
- [4] L. Modesto, Super-renormalizable quantum gravity, Phys. Rev. D 86, 044005 (2012), arXiv:1107.2403. 1 2
- [5] A. S. Koshelev, J. Marto, A. Mazumdar, Schwarzschild singularity is not permissible in ghost free quadratic-curvature infinite-derivative gravity, Phys. Rev. D 98, 064023 (2018), arXiv:1803.00309. 1 2 3
- [6] L. Buoninfante, A. S. Koshelev, G. Lambiase, J. Marto, A. Mazumdar, Conformally-flat, non-singular static metric in infinite derivative gravity, JCAP 06 (2018) 014, arXiv:1804.08195. 1 2 3 +1
- [7] A. Bonanno, D. Malafarina, A. Panassiti, Dust collapse in asymptotic safety: a path to regular black holes, Phys. Rev. Lett. 132, 031401 (2024), arXiv:2308.10890. 1 2 3
- [8] T. Harada, C.-M. Chen, R. Mandal, Singularity resolution and regular black hole formation in gravitational collapse in asymptotically safe gravity, Phys. Rev. D 111, 126017 (2025), arXiv:2502.16787. 1 2
- [9] R. Penrose, Gravitational collapse and space-time singularities, Phys. Rev. Lett. 14, 57 (1965). 1 2 3
- [10] S. W. Hawking, R. Penrose, The singularities of gravitational collapse and cosmology, Proc. R. Soc. A 314, 529 (1970). 1 2
- [11] A. Borde, Regular black holes and topology change, Phys. Rev. D 55, 7615 (1997), arXiv:gr-qc/9612057.
- [12] A. Ori, Inner structure of a charged black hole: An exact mass-inflation solution, Phys. Rev. Lett. 67, 789 (1991). 1 2
- [13] R. Carballo-Rubio, F. Di Filippo, S. Liberati, C. Pacilio, M. Visser, Inner horizon instability and the unstable cores of regular black holes, J. High Energy Phys. 05 (2021) 132, arXiv:2101.05006. 1 2 3 +1
- [14] R. Carballo-Rubio, F. Di Filippo, S. Liberati, C. Pacilio, M. Visser, Regular black holes without mass inflation instability, J. High Energy Phys. 09 (2022) 118, arXiv:2205.13556. 1 2 3 +1
- [15] D. Anselmi, Fakeons and the classicization of quantum gravity: the FLRW metric, J. High Energy Phys. 1904, 061 (2019), arXiv:1901.09273. 1 2 3
- [16] B. Bertotti, L. Iess, P. Tortora, A test of general relativity using radio links with the Cassini spacecraft, Nature 425, 374 (2003).
- [17] V. P. Frolov, Notes on nonsingular models of black holes, Phys. Rev. D 94, 104056 (2016), arXiv:1609.01758. 1 2
- [18] H. Maeda, Quest for realistic non-singular black-hole geometries: regular-center type, J. High Energy Phys. 11 (2022) 108, arXiv:2107.04791.
- [19] D. Alfyorov, I. Shnyukov, Causal Architecture Theory (CAT). The VERCETTI Principle and Canonical Branch Selection, doi:10.13140/RG.2.2.11351.02724 (2026). 1 2
- [20] D. Anselmi, Fakeons and Lee-Wick models, J. High Energy Phys. 02 (2018) 141, arXiv:1801.00915. 1 2
- [21] L. Bombelli, J. Lee, D. Meyer, R. D. Sorkin, Space-time as a causal set, Phys. Rev. Lett. 59, 521 (1987).
- [22] R. D. Sorkin, Causal sets: discrete gravity, arXiv:gr-qc/0309009 (2003).
- [23] D. M. T. Benincasa, F. Dowker, The scalar curvature of a causal set, Phys. Rev. Lett. 104, 181301 (2010), arXiv:1001.2725. 1 2
- [24] D. Alfyorov, Nonlinear field equations and FLRW limit of nonlocal quadratic gravity Theory, doi:10.13140/RG.2.2.15315.95520 (2026). 1 2
- [25] D. Alfyorov, Solar system and laboratory tests of nonlocal quadratic gravity, Authorea, doi:10.22541/au.177524450.03515205/v1 (2026).