Geometric Horizon Inﬂation - A Precise Prediction for Next-Generation Gravitational Wave Detectors

Geometric Horizon Inﬂation: An Eﬀective

Field Theory for Binary Black Hole Mergers

in an Isometric Tensor Network

Raghu Kulkarni

SSMTheory Group, IDrive Inc., Calabasas, CA 91302, USA

raghu@idrive.com

March 8, 2026

Abstract

The Hawking Area Theorem states that a black hole horizon area cannot de-

crease. Current gravitational wave observations conﬁrm this constraint to within

∼ 10%. Discrete quantum gravity models suggest the continuum approximation of

General Relativity fails at the horizon. In this paper, we construct an Eﬀective Field

Theory (EFT) using the Selection-Stitch Model (SSM). We model the vacuum as an

Isometric Tensor Network (isoTNS) on a saturated Face-Centered Cubic (K = 12)

lattice. We write the explicit lattice strain action S

isoT NS

and perform the complete

functional variation δS/δg

µν

to derive the stress-energy tensor T

µν

at the horizon.

We establish consistency with the Israel-Darmois thin-shell junction conditions by

computing the extrinsic curvature discontinuity [K

] across the horizon shell. The

area inﬂation is derived from the Raychaudhuri equation with the modiﬁed source

term, yielding a scale-invariant band of ∆A ≈ 6.86 − 7.13%. We compute the

quasinormal mode (QNM) frequency shift for the dominant (2,2,0) mode using per-

turbation theory on the Teukolsky equation: δω/ω = −ϵ/2 = −3.4%. For a 100M

⊙

remnant, this produces a ringdown frequency decrease of 4 − 6 Hz and a damping

time increase of ∼ 3.5%. These shifts match GW250114 parameter estimations and

sit within the projected sensitivity of upcoming XG observatories.

1 Introduction

Testing the Kerr hypothesis requires precision measurements of the post-merger ringdown

spectrum. Observations show broad consistency with the No-Hair theorem of General

Relativity (GR) [1]. Uniﬁcation frameworks suggest the continuum approximation breaks

down at extreme horizon curvature. Discrete vacuum models [2, 3] deﬁne the vacuum

using a lattice-like microstructure that recovers Bekenstein-Hawking entropy (S ∝ A) in

the continuum limit [4] but predicts speciﬁc structural corrections at ﬁnite lattice spacing.

Multiple frameworks exist to parametrize potential deviations from GR in gravitational

wave observations. The parametrized post-Einsteinian (ppE) framework [10] introduces

generic amplitude and phase deviations to the inspiral waveform. Parametrized tests

of the ringdown spectrum [11, 12] constrain fractional deviations in QNM frequencies

and damping times. Using O3 data, the LIGO-Virgo-KAGRA Collaboration constrained

fractional QNM frequency deviations to |δf

220

| < 0.1 at 90% credibility for the

loudest events [1]. The high-SNR detection of GW250114 has tightened these bounds.

Grimaldi et al. report an amplitude deviation of δA

= 0.06

+0.13

−0.11

[7].

We derive a quantitative prediction from the Selection-Stitch Model (SSM) [5], where

the vacuum operates as an Isometric Tensor Network (isoTNS) based on a saturated FCC

lattice (K = 12). The geometry of close-packed spheres creates a kinematic exclusion

limit—a rigid metric wall at L/

√

3 [6]—that modiﬁes the horizon boundary conditions

during a merger. Standard perturbative quantum gravity corrections scale as (l

)

and vanish for macroscopic black holes. This geometric correction functions as a macro-

scopic equation of state and persists at all mass scales.

2 The Lattice Stiﬀness from the Bekenstein-Hawking

Entropy Bound

The ﬁrst ingredient of the eﬀective action is the elastic stiﬀness κ of the lattice at the

horizon. We derive this from the Bekenstein-Hawking entropy bound, which constrains

the number of microscopic degrees of freedom per unit area.

2.1 Degrees of Freedom per Unit Area

The Bekenstein-Hawking entropy of a black hole with horizon area A is [4]:

(1)

where l

ℏG/c

is the Planck length. In natural units (k

= 1), the entropy counts

the number of independent microscopic degrees of freedom N

dof

accessible to the horizon:

dof

= S

(2)

This is the holographic bound: one binary degree of freedom per 4l

of horizon area.

In the isoTNS framework, each such degree of freedom corresponds to one lattice bond

crossing the horizon surface. The areal density of bonds is:

bond

(3)

2.2 Energy per Degree of Freedom

Each lattice bond at the horizon acts as a quantum mechanical degree of freedom in its

ground state. By the holographic principle, a horizon-crossing bond that is maximally

compressed (i.e., at the saturation limit of the lattice) carries an energy on the order of

the Planck energy E

= ℏc/l

. More precisely, the energy per bond is the energy scale

at which the lattice spacing equals the Planck length:

bond

ℏc

= E

(4)

This identiﬁcation follows from dimensional analysis constrained by the requirement that

the total energy stored in the horizon must reproduce the ADM mass M = Ac

/(16πG)

of the black hole in the continuum limit. The lattice energy density evaluates to n

bond

= E

/(4l

) = ℏc/(4l

), which must equal the stiﬀness κ.

2.3 Derivation of the Stiﬀness κ

The elastic stiﬀness is deﬁned as the energy density per unit strain. For a lattice with

bond density n

bond

and bond energy E

bond

, the energy stored per unit volume when the

lattice is strained by ϵ is:

u =

κϵ

(5)

At maximal strain (ϵ = 1, where the lattice reaches its topological limit), the stored

energy per unit area must equal the Bekenstein-Hawking surface energy density. This

requirement ﬁxes:

κ = n

bond

× E

bond

ℏc

(6)

We verify dimensional consistency: [κ] = [Energy/Length

] = [Pressure], as required for

an elastic stiﬀness. In SI units, κ = ℏc/(4l

) ≈ 4.6 × 10

113

Pa. This corresponds to

the Planck pressure divided by 4, consistent with the maximum stress-energy density

permitted by the holographic bound.

3 Eﬀective Action and Explicit Functional Variation

3.1 The Total Action

We add a localized lattice strain term to the standard Einstein-Hilbert action:

eff

√

−g



16πG

+ L

matter

+ L

isoT NS



(7)

The lattice strain Lagrangian density is:

isoT NS

κ(ξ − 1)

Θ(r

min

− r)δ(r − R

) (8)

where:

• κ = ℏc/(4l

) is the lattice stiﬀness (Eq. 6).

• ξ = 13/12 is the topological strain trace at saturation.

• Θ(r

min

−r) is the Heaviside step function that activates only when radial compression

reaches the metric wall at r

min

= L/

√

• δ(r − R

) localizes the geometric correction explicitly to the horizon boundary.

The Dirac delta function localizes the strain energy to a thin shell at the horizon. This

satisﬁes the standard thin-shell formalism utilized in the Israel junction conditions. In

the bulk spacetime (r = R

), L

isoT NS

= 0 identically, and GR is recovered exactly.

3.2 Functional Variation δS

isoT N S

/δg

µν

We perform the explicit variation of S

isoT NS

with respect to the inverse metric g

µν

. The

isoTNS action reads:

isoT NS

√

−g



κ(ξ − 1)



Φ(x) (9)

where Φ(x) = Θ(r

min

−r)δ(r −R

) is the activation function. We separate the variation

into three distinct contributions:

Term 1: Variation of

√

−g. The standard identity yields:

√

−g = −

√

−gg

µν

δg

µν

(10)

Term 2: Variation of Φ(x). The Heaviside-delta product depends on the metric

through the radial coordinate r. In Schwarzschild-like coordinates, r is a scalar function of

the spacetime point and does not depend on g

µν

at a ﬁxed coordinate position. However,

the horizon radius R

itself depends on the metric. We handle this by treating R

determined self-consistently by the modiﬁed Einstein equations. At leading order in the

strain (ξ − 1) ≪ 1, the variation of Φ(x) contributes at O((ξ − 1)

) and evaluates as

negligible. We therefore treat Φ(x) as ﬁxed under the variation.

Term 3: Variation of the strain scalar. The strain trace ξ depends on the metric

through the localized compression of the lattice. For a radially infalling conﬁguration, the

strain is related to the radial metric component:

ξ =

(0)

1/2

(11)

where g

(0)

is the unstrained reference metric. At the saturation limit, ξ reaches its topo-

logical maximum 13/12 and is mechanically clamped: δξ/δg

µν

= 0 at saturation. This

deﬁnes the metric wall—the strain cannot increase beyond the topological limit regardless

of further metric compression.

Combining these terms, the variation yields:

δS

isoT NS

δg

µν

= −

√

−g



κ(ξ − 1)

µν



Φ(x) (12)

This generates a stress-energy tensor of the form T

isoT NS

µν

= −(2/

√

−g)δS

isoT NS

/δg

µν

isoT NS

µν

κ(ξ − 1)

Φ(x)g

µν

(13)

This manifests as an isotropic pressure conﬁned strictly to the horizon shell. To decompose

it into normal and tangential components, we project using the unit normal n

to the

horizon and the induced metric γ

µν

= g

µν

+ n

isoT NS

µν

= σδ(r − R

)[ηn

+ (1 − η)γ

µν

] (14)

where σ = κ(ξ − 1) is the surface energy density and η = σ

is the projection of the

maximum compression direction onto the boundary normal. The modiﬁed Einstein ﬁeld

equation reads:

µν

= 8πG(T

matter

µν

+ T

isoT NS

µν

) (15)

4 Consistency with the Israel-Darmois Junction Con-

ditions

Any action containing δ(r − R

) structurally describes a thin shell at the horizon. The

standard GR formalism for thin shells evaluates the Israel-Darmois junction conditions

[19, 21]. We demonstrate that the stress-energy tensor (Eq. 14) is consistent with this

formalism.

4.1 Setup: The Horizon as a Thin Shell

Let Σ denote the horizon hypersurface at r = R

. The spacetime is divided into an

exterior region (r > R

) and an interior region (r < R

). The induced metric on Σ is:

= g

µν

(16)

where e

= ∂x

/∂y

are the tangent vectors to Σ and y

(i = 0, 1, 2) are the intrinsic

coordinates on the shell. The extrinsic curvature on each respective side is:

= −n

µ;ν

(17)

where n

is the unit normal pointing from − to +.

4.2 The Israel Junction Condition

The Israel junction condition structurally relates the discontinuity in extrinsic curvature

across the shell directly to the surface stress-energy tensor S

] − [K]h

= −8πGS

(18)

where [K

] = K

− K

−

is the jump in extrinsic curvature and [K] = h

] represents

the trace of the jump.

4.3 Evaluation for the isoTNS Shell

The surface stress-energy tensor S

of the isoTNS shell is obtained by integrating T

isoT NS

µν

across the shell thickness:

isoT NS

µν

dn = σ(1 − η)h

(19)

where the normal-normal component σηn

does not project onto Σ. The surface stress-

energy evaluates purely tangential: S

= ρ

with surface energy density ρ

= σ(1−η) =

κ(ξ − 1)(1 − η). Substituting this into the Israel condition (Eq. 18):

] − [K]h

= −8πGκ(ξ −1)(1 − η)h

(20)

Taking the trace (h

= 3 for a 3D hypersurface):

− 3[K] = −24πGκ(ξ − 1)(1 − η) (21)

Because [K

= [K], this simpliﬁes to −2[K] = −24πGκ(ξ −1)(1 − η). Therefore, the

extrinsic curvature jump evaluates to:

[K] = 12πGκ(ξ − 1)(1 − η) (22)

This manifests as a well-deﬁned, ﬁnite discontinuity strictly proportional to the topological

strain (ξ −1) = 1/12. The Israel junction conditions are explicitly satisﬁed by the isoTNS

stress-energy tensor. The shell exhibits a positive surface energy density (ξ > 1 and

η < 1), corresponding geometrically to a stretched horizon that inﬂates the area.

4.4 Veriﬁcation: No-Hair Theorem Compatibility

The No-Hair theorem mandates that the asymptotic spacetime is characterized solely by

M, J, and Q. Because T

isoT NS

µν

= 0 for r = R

, the exterior spacetime remains an exact

Kerr metric. The ADM mass M

and dimensionless spin χ

measured at spatial inﬁnity

are unmodiﬁed. The isoTNS shell modiﬁes only the speciﬁc relationship between M

and the local irreducible mass M

irr

at the horizon, inﬂating the latter by the scale factor

√

1 + ϵ.

5 Geometric Derivation of the Area Inﬂation

5.1 The Kinematic Exclusion Limit

The body diagonal of the FCC unit cell establishes a strict minimum radial compression

limit. For close-packed hard spheres of diameter L, the deepest (111) void possesses a

depth of L/

√

3 ≈ 0.577L. No adjacent pair of nodes can compress below this physical

limit [6]. This constitutes a kinematic exclusion boundary analogous to hard-core nuclear

repulsion in neutron star equations of state. The maximum internal compression vector

aligns precisely with the (111) diagonal: σ = (1, 1, 1)/

√

5.2 The Topological Strain Trace ξ = 13/12

A saturated FCC unit cell maintains K = 12 nearest-neighbor nodes surrounding 1 cen-

tral interstitial void. Under extreme gravitational compression, the lattice structure ap-

proaches the theoretical limit of ﬁlling this 13th coordinate position. The maximum

volumetric strain is ξ = V

final

initial

= 13/12, representing the exact excess volume ratio

when the interstitial void is saturated. This topological invariant is independent of the

absolute lattice spacing L.

5.3 Derivation from the Raychaudhuri Equation

We derive the precise area inﬂation from the Raychaudhuri equation, which dictates the

evolution of the expansion scalar θ of a congruence of null generators mapped onto the

horizon [20]. The equation evaluates to:

dθ

dλ

= −

− σ

µν

− R

µν

(23)

where θ = ∇

deﬁnes the expansion, σ

µν

is the shear tensor, and R

µν

is the Ricci tensor.

Applying the modiﬁed Einstein equation (Eq. 15), the Ricci tensor becomes:

µν

= 8πG



µν

−

µν



+ Λg

µν

(24)

The speciﬁc contribution from the isoTNS shell to the geometric focusing term evaluates

to:

isoT NS

µν

= 8πG



isoT NS

µν

−

µν

isoT NS



(25)

For a strict null vector, g

µν

= 0, nullifying the trace term. Substituting the localized

decomposition (Eq. 14):

isoT NS

µν

= σδ(r − R

)[η(n

)

+ (1 − η)γ

µν

] (26)

Integrating the Raychaudhuri equation across the horizon shell evaluates the discrete

jump in expansion:

∆θ = −8πGσ[η(n · k)

+ (1 − η)(γ · k · k)] (27)

The expansion parameter θ = (1/A)(dA/dλ) maps the integrated ∆θ directly to the

fractional area change:

∆A

= ∆θ ×∆λ (28)

5.4 Orthogonal Strain Partition

Because the rigid metric wall halts further radial compression, the stress structurally

partitions into two orthogonal geometric modes. The normal mode maps the fraction

η = σ ·n of the compression vector projected onto the boundary normal. This component

carries the linear volumetric strain ξ = 13/12 and is strictly blocked by the metric wall.

The transverse mode evaluates the remaining fraction (1 −η) shunted into the 2D horizon

surface. A 2D area accommodating a 3D volumetric strain scales precisely as A ∝ V

2/3

The transverse area correction evaluates to ξ

2/3

= (13/12)

2/3

The geometric area inﬂation factor operates as the linear superposition of these or-

thogonal contributions:

∆A(η)

= η





+ (1 − η)





2/3

− 1 (29)

For a Cartesian boundary patch with η = 1/

√

3 ≈ 0.577, ∆A = 7.13%. For an isotropic

spherical average (η =

cos θ sin θdθ/

sin θdθ = 1/2), ∆A = 6.86%. This establishes

the scale-invariant prediction band:

SSM

= A

Kerr

, χ

) × [1.0686 to 1.0713] (30)

6 Quasinormal Mode Frequency Shift from Teukol-

sky Perturbation Theory

6.1 The Teukolsky Equation and Boundary Conditions

The QNM spectrum of a Kerr black hole is strictly governed by the Teukolsky equation [13]

for spin-weighted perturbations. For the dominant (2,2,0) mode, the radial equation

formulates as:

∆

+ V (r)R = 0 (31)

where ∆ = r

− 2Mr + a

and V (r) deﬁnes the eﬀective potential. The QNM boundary

conditions require purely ingoing waves at the horizon. In the tortoise coordinate r

∗

, this

condition evaluates as R ∼ e

−iωr

∗

as r

∗

→ −∞ (r → r

6.2 Modiﬁed Horizon Position

The area inﬂation ϵ = ∆A/A mechanically modiﬁes the eﬀective irreducible mass:

eff

irr

= M

irr

√

1 + ϵ (32)

The horizon area deﬁnes the relationship between irreducible mass and the horizon radius.

For a Schwarzschild conﬁguration, A = 16πM

and r

= 2M. The inﬂated horizon radius

becomes:

eff

= 2M

√

1 + ϵ ≈ r



1 +

+ O(ϵ

)



(33)

The tortoise coordinate mapping near the modiﬁed horizon adjusts to:

eff

∗

= r + 2M

√

1 + ϵ ln



√

1 + ϵ

− 1



(34)

6.3 Perturbative QNM Frequency Shift

QNM frequencies manifest as precise eigenvalues of the Teukolsky equation. For the

Schwarzschild case, the (2,2,0) mode scales fundamentally as [14]:

220

+ iω

(35)

where ω

and ω

are dimensionless constants. The frequency maps inversely proportional

to the horizon radius. Substituting r

eff

√

1 + ϵ

≈ ω



1 −



(36)

The precise fractional frequency shift evaluates to:

δω

= −

+ O(ϵ

) (37)

This derivation holds symmetrically for both the real (oscillation) and imaginary (damp-

ing) components. At the midpoint of the prediction band (ϵ = 0.07), the shift yields

δω/ω ≈ −3.4%. The damping time τ = 1/|ω

| increases concurrently by +3.5%. The

identical scaling (ω ∝ 1/M) holds rigidly at a ﬁxed dimensionless spin χ = a/M [14, 15].

The isoTNS area inﬂation modiﬁes the eﬀective mass at ﬁxed χ, preserving the exact

fractional shift:

δω



Kerr

= −

+ O(ϵ

, ϵχ

) (38)

6.4 Numerical Evaluation

7 Observational Consistency with GW250114

The GW250114 detection provides the highest-SNR test of post-merger dynamics to date

[1]. The SSM analytic predictions map consistently to four independent observational

bounds derived from the event:

Parameter Schwarzschild Kerr (χ = 0.7) Units

M (GR) 0.3737 0.5349 1/M

M (GR) -0.0890 -0.1768 1/M

δω

/ω

-3.4% -3.3% –

δω

/ω

-3.4% -3.3% –

220

= 100M

⊙

) 121 → 117 173 → 167 Hz

220

= 100M

⊙

) 5.52 → 5.71 2.83 → 2.93 ms

Table 1: QNM frequency shifts from geometric area inﬂation for the dominant (2,2,0)

mode. Values are derived using standard unperturbed QNM frequencies [14,15] scaled by

δω/ω = −ϵ/2.

• Amplitude shift: Grimaldi et al. [7] report δA

= 0.06

+0.13

−0.11

. The SSM area

inﬂation geometrically expands the boundary source, dictating a direct amplitude

scaling of

√

1 + ϵ − 1 ≈ +3.4%. This shift aligns symmetrically with the +6%

median of the posterior distribution.

• Inferred mass shift: The LVK collaboration bounds ∆M

= 0.02 ± 0.04 [1].

Standard waveform templates assume A ∝ M

. Parameter estimation executed

against a ∼ 7% inﬂated area mathematically forces an apparent mass inﬂation of

√

1.07 − 1 = 3.4%, remaining well within the 1σ margin.

• Pre-merger anomaly: Lu et al. [8] isolate a localized residual strain peak at t ≤

−7M

. In the SSM geometric framework, this precise timing marks the threshold

where local compression forces contact with the L/

√

3 metric wall, shunting latent

kinetic energy into the 2D horizon prior to peak emission.

• Absence of echoes: Wu et al. [9] locate no statistical evidence for late-time echoes.

The metric wall dissipates extreme stress exclusively via transverse, in-plane lattice

stretch (the area inﬂation) rather than reﬂective bulk bounces, mechanically sup-

pressing echo generation.

8 Scale Invariance and the Continuum Limit

The SSM area inﬂation diverges structurally from standard perturbative quantum gravity

corrections. Standard loop corrections scale strictly as (l

)

and vanish asymptomat-

ically for macroscopic astrophysical black holes. The SSM area inﬂation maps as a rigid

geometric constant of 6.86 − 7.13%, entirely independent of remnant mass.

The metric wall operates mathematically as a macroscopic equation of state. The

macroscopic horizon area is the integral sum of discrete fundamental cells: A =

cell

The fractional inﬂation evaluates to ∆A/A = ∆a

cell

. The absolute lattice spacing

L strictly cancels out of this ratio, yielding a pure topological invariant mapped by the

strain limit ξ = 13/12. The macroscopic Lorentz invariance of the bulk is strictly preserved

because the underlying 3D network operates as an isometric tensor network generated by

a continuous 2D boundary [6]. In the unstrained bulk exterior, the isoTNS stress trace

evaluates exactly to zero (T

isoT NS

µν

= 0), and classical GR is recovered ﬂawlessly.

9 Conclusion and Falsiﬁability

This paper formulates a rigorous Eﬀective Field Theory (EFT) for the K = 12 isoTNS

vacuum precisely at the horizon of a binary black hole merger. The lattice elastic stiﬀness

κ = ℏc/(4l

) was derived directly from the Bekenstein-Hawking entropy bound without

introducing free parameters. The explicit functional variation of the localized lattice

strain action yields a horizon stress-energy tensor that strictly satisﬁes the Israel-Darmois

junction conditions for a dynamic thin shell.

Integrating the Raychaudhuri equation against this modiﬁed source dictates a rigid,

mass-independent area inﬂation of ∆A = 6.86 − 7.13%. Teukolsky perturbation theory

translates this inﬂation into a fractional Quasinormal Mode shift of δω/ω = −3.4%. For

a 100M

⊙

remnant at χ = 0.7, the frequency structurally decreases by 5.8 Hz, while the

damping time increases by 0.10 ms. The recent GW250114 event data accommodates

both the strict GR null hypothesis (0%) and the precise SSM prediction. Next-generation

observatories (Cosmic Explorer, Einstein Telescope) are projected to bound fractional

QNM deviations to < 2% [16]. A deﬁnitive, high-SNR measurement of |δf

220

| < 2%

lacking systematic oﬀsets will explicitly falsify this geometric inﬂation framework.

A Self-Contained SSM Summary

A.1. K = 12 Lattice Saturation. The FCC lattice deﬁnes the mathematically unique

solution to the Kepler conjecture [17]. The absolute densest packing of identical uniform

spheres in three dimensions maps a coordination of K = 12. The vacuum tensor network

naturally saturates at this theoretical maximum, yielding exactly 12 nearest-neighbor

bonds per node of length L/

√

A.2. The Metric Wall at L/

√

3. The FCC unit cell houses its deepest interstitial

void positioned along the internal (111) body diagonal. For rigid spheres of geometric

diameter L, the absolute minimum center-to-center distance permitted along this diagonal

evaluates to L/

√

3. This threshold deﬁnes a strict kinematic exclusion limit [6].

A.3. Isometric Tensor Networks and Lorentz Invariance. The 3D bulk lattice

is generated as a quasilocal isometric projection of a purely continuous 2D boundary,

mirroring the Ryu-Takayanagi formal prescription [18]. The exact isometry functionally

maps boundary entanglement entropy to bulk geodesic areas. The macroscopic poly-

crystalline grain structure of the bulk statistically averages across all lattice orientations,

extinguishing preferred spatial directions at scales exceeding L ∼ 1 fm.

References

[1] A. G. Abac et al. (LVK Collaboration), “Black Hole Spectroscopy and Tests of Gen-

eral Relativity with GW250114,” Phys. Rev. Lett. 136, 041403 (2026).

[2] G. ’t Hooft, The Cellular Automaton Interpretation of Quantum Mechanics, Springer

(2016).

[3] G. E. Volovik, The Universe in a Helium Droplet, Oxford University Press (2003).

[4] J. D. Bekenstein, “Black holes and entropy,” Phys. Rev. D 7, 2333 (1973).

[5] R. Kulkarni, “Late-Universe Dynamics from Vacuum Geometry,” Zenodo:

10.5281/zenodo.18753874 (In Review) (2026).

[6] R. Kulkarni, “Constructive Veriﬁcation of K=12 Lattice Saturation,” Zenodo:

10.5281/zenodo.18294925 (In Review) (2026).

[7] L. Grimaldi et al., “Plunge-Merger-Ringdown Tests of GR with GW250114,”

arXiv:2601.13173 (2026).

[8] N. Lu et al., “GW250114 reveals black hole horizon signatures,” arXiv:2510.01001

(2025).

[9] D. Wu et al., “Model-independent search of gravitational wave echoes,”

arXiv:2512.24730 (2025).

[10] N. Yunes and F. Pretorius, “Fundamental theoretical bias in gravitational wave as-

trophysics and the ppE framework,” Phys. Rev. D 80, 122003 (2009).

[11] A. Ghosh et al., “Testing general relativity using golden black-hole binaries,” Phys.

Rev. D 94, 021101(R) (2016).

[12] M. Isi et al., “Testing the black-hole area law with GW150914,” Phys. Rev. Lett.

127, 011103 (2021).

[13] S. A. Teukolsky, “Perturbations of a rotating black hole. I,” Astrophys. J. 185, 635

(1973).

[14] E. Berti, V. Cardoso, and A. O. Starinets, “Quasinormal modes of black holes and

black branes,” Class. Quant. Grav. 26, 163001 (2009).

[15] E. W. Leaver, “An analytic representation for the quasi-normal modes of Kerr black

holes,” Proc. R. Soc. Lond. A 402, 285 (1985).

[16] M. Punturo et al., “The Einstein Telescope: a third-generation gravitational wave

observatory,” Class. Quant. Grav. 27, 194002 (2010).

[17] T. C. Hales, “A proof of the Kepler conjecture,” Annals Math. 162, 1065 (2005).

[18] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from

AdS/CFT,” Phys. Rev. Lett. 96, 181602 (2006).

[19] W. Israel, “Singular hypersurfaces and thin shells in general relativity,” Nuovo Cim.

B 44, 1 (1966); Erratum: Nuovo Cim. B 48, 463 (1967).

[20] R. M. Wald, General Relativity, University of Chicago Press (1984).

[21] E. Poisson, A Relativist’s Toolkit, Cambridge University Press (2004).

Geometric Horizon Inﬂation: A Universal Prediction for Binary Black Hole Mergers