Preprint
Article

This version is not peer-reviewed.

Causally Confined Euclidean Saddles in Spin-Foam Quantum Gravity

Submitted:

12 April 2026

Posted:

14 April 2026

You are already at the latest version

Abstract
We investigate the semiclassical structure of spin-foam transition amplitudes for boundary data that do not admit a real Lorentzian Regge geometry. Considering a fixed triangulation with a single dominant vertex, we demonstrate that when boundary tetrahedra carry mutually incompatible causal orientations, the closure equations have no real solution and the path integral is dominated by a complex Euclidean saddle of the Regge action. In this regime the vertex amplitude acquires a non-oscillatory factor of the form exp(−SE/ℏ), where SE is the Euclidean action evaluated at the complex saddle. We introduce a causal-obstruction criterion based on a convexity argument for the future timelike cone in R 3,1 , and establish a formal classification of boundary data into three types according to the existence and nature of the saddle-point solutions. We show that SE scales linearly with the spin parameter j in the semiclassical limit, SE = ℏ j C(α)/(8πG), where C(α) is a finite dimensionless geometric constant, providing explicit control over the suppression. Non-degeneracy of the Hessian at the complex saddle is verified after gauge fixing, confirming the validity of the saddle-point approximation. The results constitute a proof-of-concept demonstration that exponentially suppressed, causally confined quantum-geometric transitions emerge as a structural feature of the covariant formulation of loop quantum gravity, without additional postulates.
Keywords: 
;  ;  ;  ;  ;  

1. Introduction

Spin-foam models provide a covariant, background-independent formulation of quantum gravity in which transition amplitudes between quantum geometry states are expressed as discrete path integrals over two-complexes [1,2,3,4]. In the semiclassical large-spin regime, the asymptotic analysis of vertex amplitudes is well established: when boundary data are compatible with a non-degenerate Lorentzian Regge geometry, the dominant contributions arise from real saddles of the Regge action, producing oscillatory amplitudes that reproduce the Einstein equations in the classical limit [5,6]. The vertex amplitude analysed throughout this work is that of the EPRL model [7,29], the principal covariant realisation of this programme.
A natural but less-explored question concerns the behaviour of the path integral when this compatibility fails. Spin-foam amplitudes are formally defined for arbitrary boundary data, including configurations that correspond to no real Lorentzian geometry simultaneously satisfying the closure and gluing conditions. Such configurations are not pathological: they arise generically in the interior of a foam whenever the causal structure imposed by the boundary data cannot be realised by any classical geometry.
We address this question directly, within a simplified model with a fixed triangulation. We show that when boundary tetrahedra carry mutually incompatible causal orientations, the path integral is dominated by a complex saddle at which the Regge action becomes purely imaginary, giving rise to a contribution of the form
A v ( det H ) 1 / 2 exp S E ,
where S E > 0 is the Euclidean action at the complex saddle and H is the Hessian of quadratic fluctuations.
The scope of the results should be delimited from the outset. The demonstration is carried out for a toy model with a fixed triangulation and a single dominant internal vertex. We do not claim universality of the mechanism for arbitrary foams, nor do we perform a sum over triangulations. The value of the present work lies in the constructive, mathematically controlled demonstration that localised Euclidean saddles are structurally compatible with the spin-foam formalism, and in the identification of precise geometric conditions under which they emerge.
The mechanism has three properties that distinguish it from related proposals in the literature [9,10,11]. First, it is strictly local: the complex saddle is associated with a single foam vertex, confined to a finite causal domain, without requiring a global Euclidean continuation of spacetime. Second, it is derivable within the standard formalism: the factor exp ( S E / ) is not introduced as a postulate, but emerges from the analytic structure of the path integral. Third, it is semiclassically controllable: in the j 1 limit, the suppression satisfies S E = j C ( α ) / ( 8 π G ) with C ( α ) finite and computable.
The paper is organised as follows. Section 2 develops the causal obstruction criterion and the formal classification of boundary data. Section 3 derives the Euclidean saddle and the scaling of S E . Section 4 presents the minimal realisation and stability results. Section 5 discusses the scope and relation to the literature. The appendices contain the saddle equations (A), the Hessian calculation (B), the refinement stability proof (C), and the calibration of the Gram criterion (D).

2. Causal Obstruction in Boundary Data

2.1. Geometric Setup and Notation

We consider a non-degenerate 4-simplex σ with five vertices { v 0 , v 1 , v 2 , v 3 , v 4 } and five boundary tetrahedra { τ a } a = 0 4 , where τ a is the convex hull of the four vertices excluding v a . To each pair a < b we associate a triangular face f a b with spin j a b 1 2 N and area
A a b = γ j a b ( j a b + 1 ) γ j a b , j a b 1 ,
where γ is the Barbero–Immirzi parameter. To each tetrahedron τ a we associate an outward normal n a R 3 , 1 with n a · n a = 1 (Minkowski metric, signature ( , + , + , + ) ). We say n a is future-pointing if n a 0 > 0 and past-pointing if n a 0 < 0 .
In the large-spin semiclassical regime, the dominant contributions to the EPRL vertex amplitudes are determined by the saddle-point conditions [5,6]:
Closure conditions. For each tetrahedron τ a ,
b a j a b n a b = 0 ,
where n a b is the normal to face f a b as seen from τ a .
Gluing conditions. For each face f a b , there exists g a b SO ( 1 , 3 ) such that
g a b n a b = g b a n b a .
When (3) and (4) admit a solution with all variables real, the saddle is real Lorentzian and the amplitude is oscillatory. The regime in which no such solution exists is the central object of this work.

2.2. Causal Obstruction: Convexity Argument

Lemma 2.1
(Convexity of the future timelike cone). The set C + = { u R 3 , 1 : u · u < 0 , u 0 > 0 } is an open convex cone. In particular, if u 1 , , u k C + and λ 1 , , λ k > 0 , then i = 1 k λ i u i C + .
Proof. 
Given u , w C + and λ ( 0 , 1 ) , the reverse Cauchy–Schwarz inequality for the Minkowski metric gives | u · w | | u | | w | , with | u | = u · u > 0 . Hence ( λ u + ( 1 λ ) w ) 0 > 0 and ( λ u + ( 1 λ ) w ) · ( λ u + ( 1 λ ) w ) < 0 , so λ u + ( 1 λ ) w C + . The extension to finite combinations follows by induction.    □
Theorem 2.2
(Causal obstruction to closure). Consider boundary data for τ a in which k normals { n a b } b a are future-pointing and 4 k are past-pointing, with 1 k 3 . Then, for generic spins { j a b } , the closure condition (3) admits no solution with all normals real Lorentzian. The set of spins for which the obstruction fails has measure zero.
Proof. 
Without loss of generality, suppose n a , i C + for i = 0 , , k 1 and n a , i C for i = k , , 3 . Condition (3) reads
i = 0 k 1 j a i n a i C + = i = k 3 j a i n a i C .
By Lemma 2.1, the left-hand side belongs to C + and the right to C = C + . Since C + C = , the equality cannot hold. The only exception occurs when both sides vanish simultaneously, a condition of measure zero in the parameter space.    □
Corollary 2.3.
The obstruction of Theorem 2.2 is stable: it persists under generic perturbations of the spins j a b and normals n a b , provided the temporal orientations are not altered.

2.3. Classification of Boundary Data

Definition 2.4
(Causal classification). Boundary data { j a b , n a } for a vertex v are classified as:
  • Type I: conditions (3)–(4) admit a real Lorentzian solution. Oscillatory amplitude: A v ( I ) exp ( i S Regge / ) .
  • Type II: no real solution exists, but a complex solution exists with Re ( S Regge [ X ] ) 0 and Im ( S Regge [ X ] ) 0 . Oscillatory amplitude with partial suppression.
  • Type III: the data satisfy Theorem 2.2 and the complex solution yields purely imaginary dihedral angles,
    θ a b i χ a b , χ a b R + ,
    so that S Regge [ X ] = i S E with S E > 0 real. Non-oscillatory amplitude: A v ( III ) exp ( S E / ) .
The assignment of causal orientations to boundary tetrahedra in spin-foam models has been studied in the context of implementing causality constraints [26,27]; the present classification identifies the specific sub-class of causally obstructed data for which the path integral is dominated by a complex saddle.
Remark 2.5.
The analytic continuation (6) does not represent a Wick rotation of spacetime. It is a stationary-phase trajectory in the complexified space of discrete geometric variables, with no interpretation as a classical Euclidean geometry. The results of the following sections apply exclusively to Type-III data.

2.4. Numerical Implementation: Gram Matrix Criterion

For each tetrahedron τ a , define the 4 × 4 Gram matrix
G a = n a b · n a c b , c a .
The closure condition (3) implies det ( G a ) 0 in the semiclassical limit for Type-I vertices; for Type III, | det ( G a ) | = O ( 1 ) . The operational criterion is
v Type III | det ( G v ) | > δ crit ( v ) ,
with threshold
δ crit ( v ) = α e 2 j min ( v ) 1 + β d ( v , F ) ,
where j min ( v ) is the smallest spin adjacent to v, d ( v , F ) is the combinatorial distance to the boundary, and α 10 3 , β 0.1 are calibrated in Appendix D.

3. The Euclidean Saddle and the Scaling of the Action

3.1. Analytic Continuation of Geometric Variables

For Type-III data, Theorem 2.2 guarantees that no real Lorentzian solution of (3)–(4) exists. We extend the integration domain to the complexified space X C n . In the complex domain, the dihedral angle θ a b satisfies cos θ a b = ( n a · n b ) / ( | n a | | n b | ) with a complex bilinear inner product. For the parametrisation θ a b = i χ a b , this becomes cosh χ a b = ( n a · n b ) / ( | n a | | n b | ) .
The complex saddle does not correspond to any real classical geometry. It represents a stationary-phase trajectory in the complexified space of discrete geometries; the exponential suppression is a consequence of the analytic structure of the amplitude, not the action of an intermediate Euclidean geometry.

3.2. The Regge Action at the Complex Saddle

The Regge action for the 4-simplex [8] is
S Regge = 1 8 π G a < b A a b Θ a b ,
where Θ a b = 2 π σ f a b θ a b ( σ ) is the deficit angle. Substituting (6) into (10),
S Regge [ X ] = 1 8 π G a < b A a b ( 2 π i χ a b ) = S top + i S E ,
where S top is a real topological constant (an irrelevant global phase) and
S E 1 8 π G a < b A a b χ a b > 0 .
Positivity of S E follows from A a b > 0 and χ a b > 0 for all pairs.

3.3. Saddle-Point Conditions in the Complexified Domain

The critical-point equations δ S Regge / δ X | X = 0 are preserved under analytic continuation [5]. The demonstration that Type-III data admit a solution with θ a b = i χ a b is given in Appendix A, where the closure and gluing equations in the complexified domain are solved explicitly and the values χ a b are determined by the prescribed boundary areas and temporal orientations.

3.4. Validity of the Saddle-Point Approximation

The saddle-point approximation requires that the complex saddle X contribute effectively to the original path integral. In the Conrady–Freidel parametrisation [6,22], the vertex amplitude integrand is holomorphic in the complexified group variables g a SL ( 2 , C ) and coherent spinors, with potential non-analyticities only at zeros of the scalar products ξ a b | g a 1 g b | ξ b a , which form a set of measure zero.
In this setting, the steepest-descent (Picard–Lefschetz) decomposition applies: the original integration contour can be expressed as a linear combination of thimbles J X i associated with the critical points X i of S CF  [17,18,19,20,21]. For Type-III boundary data, Theorem 2.2 and the non-degeneracy of the Hessian (Proposition 4.2) together guarantee that X is an isolated, non-degenerate critical point of S CF in the complexified domain. The associated thimble is therefore well-defined and contributes to the decomposition of the original contour with a non-zero coefficient.
In the semiclassical limit j , the contribution from X dominates over all other saddles because Re ( S CF [ X ] ) = 2 j a b χ a b < 0 is strictly negative and grows in magnitude as O ( j ) , while the contributions from any other saddle with smaller | Re ( S CF ) | are exponentially suppressed relative to it. This justifies the leading-order approximation (1).
Proposition 3.1
(Non-vanishing Lefschetz coefficient). Let X be the Type-III saddle constructed above. Then the associated Lefschetz coefficient is non-zero:
n ( σ , J X ) = ( 1 ) q = ± 1 ,
where q is the number of negative eigenvalues of H res .
Proof. 
The Conrady–Freidel action S CF is holomorphic in the complexified variables. We verify the two hypotheses of the real-saddle theorem of Picard–Lefschetz theory [17,18] (see [17], §3):
(i) X σ . From (A31)–(A32), χ L = arccos ( 1 / 4 ) and χ S = arccosh ( 3 cosh χ L ) are both real. The gluing group elements are pure hyperbolic boosts g a = cosh ( χ / 2 ) sinh ( χ / 2 ) sinh ( χ / 2 ) cosh ( χ / 2 ) with real rapidity χ { χ L , χ S } ; hence g a SL ( 2 , R ) SL ( 2 , C ) and X lies on the original real integration contour σ .
(ii) X is non-degenerate. The gauge-fixed Hessian satisfies det H res 0 (Proposition 4.2).
By the real-saddle theorem [17], the intersection number equals ( 1 ) q , where q { 0 , , 29 } is the (finite, well-defined) number of negative eigenvalues of the real, non-degenerate matrix H res . Therefore n ( σ , J X ) = ( 1 ) q = ± 1 0 .    □

3.5. Semiclassical Scaling of S E

Scale parametrisation.
A a b = j α a b ,
where j 1 is the semiclassical parameter and α a b > 0 characterise the shape of the 4-simplex independently of its overall scale.
Angular invariance. Dihedral angles depend only on the shape of the simplex, not on its scale. This property extends to the analytic continuation: χ a b = χ a b ( { α c d } ) and χ a b / j = 0 .
Main result. Substituting (13) into (12),
S E = j 8 π G C ( α ) ,
where the dimensionless geometric constant is
C ( α ) a < b α a b χ a b ( { α } ) .
Proposition 3.2
(Finiteness of C ( α ) ). For a non-degenerate 4-simplex with Type-III data, 0 < C ( α ) < generically.
Proof. 
The set { α a b } has ten finite positive terms, and χ a b > 0 is regular for non-degenerate data. Finiteness and positivity of C ( α ) follow directly.    □

3.6. Quadratic Fluctuations

Expanding the action around X ,
S Regge [ X ] = i S E + 1 2 δ X I H I J δ X J + O ( δ X 3 ) ,
where H I J = 2 S Regge / X I X J | X . Gaussian integration gives
A v ( det H ) 1 / 2 exp S E j n / 2 exp j C ( α ) 8 π G , j ,
recovering (1), where n is the number of physical degrees of freedom after complete gauge fixing (computed in Appendix B). The calculation of H and the verification that det H 0 are presented in Appendix B.

4. Minimal Realisation

4.1. Construction of the Minimal Spin Foam

We consider a spin foam F with a single internal vertex v , five internal tetrahedra { τ a } a = 0 4 , and fixed boundary F . The triangulation is fixed; no sum over triangulations is performed.
The boundary data are of Type III: three future-pointing tetrahedra ( n a 0 > 0 for a = 0 , 1 , 2 ) and two past-pointing ( n a 0 < 0 for a = 3 , 4 ). By Theorem 2.2, no real Lorentzian solution exists for generic spins.
Definition 4.1
(Quantum core). Given a finite spin foam F with vertex set V ( F ) , thequantum coreis F core = { v V ( F ) : v has Type - III data } and the external sector is F ext = V ( F ) F core .
In the toy model, F core = { v } and the total amplitude is
A ( F ) = A v ( det H ) 1 / 2 exp S E ,
exact within the semiclassical regime, without a factorisation approximation.
For foams with multiple vertices, the amplitude admits the approximate factorisation A A ext · A core , where A ext is oscillatory and A core is suppressed. The analysis of this factorisation for non-trivial interfaces is beyond the scope of this work.

4.2. Locality of the Mechanism

Although the analysis has been carried out for a minimal triangulation with a single internal vertex, the mechanism is intrinsically local: the relevant saddle-point conditions are entirely determined by the boundary data of a single effective 4-simplex, and the value of S E depends only on the local areas { A a b } and the corresponding dihedral angles.
It follows that, in any spin foam containing vertices with mixed causal orientations, each such vertex contributes an exponentially suppressed factor exp ( S E / ) to the amplitude, independently of the global structure of the foam. The model considered here is therefore not an artefact of the minimal triangulation, but the minimal realisation of a local mechanism that operates generically within the spin-foam formalism.

4.3. Explicit Example: j S = 1 , j L = 3

For the ( 3 + 2 ) distribution, Appendix A shows that Theorem 2.2 guarantees the absence of a real saddle for any ρ = j L / j S > 0 with generic spins, with no critical threshold. For the concrete example j S = 1 , j L = 3 ( ρ = 3 ), the fixed-point system (A27) admits an exact closed-form solution (Appendix A, §A.4):
χ L = arccos 1 4 1.318 rad , χ S = arccosh 3 cosh arccos 1 4 2.479 rad ,
where χ L χ FP and χ S χ FF = χ PP . The geometric constant (15) is
C ( α ) = 4 χ S + 18 χ L = 4 × 2.479 + 18 × 1.318 = 33.64 .
The resulting amplitude is
A v j n / 2 exp 33.64 j 8 π G ,
demonstrating explicit exponential suppression without free parameters.

4.4. Non-Degeneracy of the Hessian

Proposition 4.2
(Non-degeneracy of the Hessian). For Type-III data with distribution ( k + ( 4 k ) ) , 1 k 3 , and non-degenerate geometry, the Hessian H I J evaluated at the complex saddle X is non-degenerate after gauge fixing adapted to the isotropy group of the saddle: det H res 0 .
The proof is given in Appendix B. The argument is that the asymmetric distribution of temporal orientations breaks the continuous stabiliser symmetry, eliminating the additional null modes present in the Type-I case.

4.5. Stability Under Local Refinement

Theorem 4.3
(Stability under refinement). Let F be the minimal foam with v of Type III and j a b = j 1 . Let F be a stellar subdivision of v with N internal vertices and interior central vertex c, satisfying the single hypothesis:
(H1)
j int j in the semiclassical limit.
Then all internal vertices of F have Type-III data (Proposition A3, Corollary A5), and the amplitude satisfies
A ( F ) N n / 2 exp j C ( α , N ) 8 π G ,
with C ( α , N ) N C ( α ) 1 + O ( N 1 / 4 ) .

5. Discussion

5.1. Relation to Standard Asymptotic Analysis

The results are a natural extension of the asymptotic analysis of [5,6,14,15,16,25,28]. For Type-I data, the amplitude is dominated by two real saddles, A v ( I ) ± ( det H ± ) 1 / 2 e ± i S Regge / . The present work treats the complementary case: Type-III data for which no real saddle exists. The mathematical structure is the same — saddle-point analysis in the complexified domain — but the result is qualitatively distinct: instead of oscillatory contributions, one obtains an exponentially suppressed, non-oscillatory contribution.

5.2. Relation to Existing Approaches

Haggard and Rovelli (2015). Reference [9] proposes a global Euclidean continuation of spacetime to describe black-to-white hole tunnelling. The present work differs fundamentally: the Euclidean saddle is local, associated with a single vertex, requiring no global Euclidean continuation. The mechanism derived here neither requires nor implies a reversal of causal time outside F core . The physical context in which Type-III transitions are expected to operate includes compact quantum-gravitational objects such as Planck stars [30], where the quantum-bounce mechanism requires causally non-classical transitions of precisely the type analysed here.
Bianchi, Christodoulou and Di Biagio (2023). Reference [10] obtains black-to-white hole amplitudes in spin foams, with suppression weights introduced as a hypothesis or estimated by entropy counting. The key distinction is that here the factor exp ( S E / ) is derived from the analytic structure of the path integral under explicit geometric conditions, without additional postulates.
Christodoulou, Di Biagio and Bianchi (2024). Reference [11] analyses the tunnelling time given that the weights exist. The two works are complementary: the present one provides a microscopic mechanism for the weights, while [11] analyses their observational consequences.
Donà, Speziale and Zhang (2019); Donà, Gozzini and Speziale (2020). References [12,13] analyse EPRL amplitudes for configurations that include timelike tetrahedra, establishing asymptotic and numerical results beyond the spacelike regime of [5]. The present work identifies a specific sub-class of such configurations — those with mutually incompatible causal orientations (Type-III data) — and demonstrates that they produce a qualitatively distinct, exponentially suppressed contribution, absent from the oscillatory saddles studied in [12,13].

5.3. Technically Distinctive Elements

The present work introduces three elements not present in this form in the literature. First, a formal classification of boundary data into Types I, II, and III with provable genericity and stability properties (Theorem 2.2, Corollary 2.3). Second, the derivation of the scaling S E = j C ( α ) / ( 8 π G ) with C ( α ) identified without free parameters (Proposition 3.2). Third, the refinement stability Theorem 4.3 with explicit hypotheses.

5.4. Scope and Limitations

Within the scope of the toy model, the results are complete: the causal obstruction is demonstrated, the Euclidean saddle is constructed, the scaling of S E is derived, the Hessian is non-degenerate, and stability under refinement is established. No result depends on conjectures or free parameters beyond the boundary data.
Three natural extensions remain open. The analysis of Type-II data requires stationary-phase techniques for actions with simultaneous real and imaginary parts. The validity of the factorisation A A ext · A core for foams with multiple vertices and non-trivial interfaces requires control over shared degrees of freedom. Finally, the sum over triangulations remains a direction for future work.

Acknowledgments

The author conducted this research independently. No external funding or institutional support was received.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Saddle Equations and Exact Solution for the (3+2) Distribution

Appendix A.1. Vertex Action in the Conrady–Freidel Parametrisation

In the Conrady–Freidel parametrisation [6,22], the EPRL vertex amplitude in the large-spin regime is
A v = SL ( 2 , C ) 5 a d g a exp S CF ,
with vertex action
S CF = a < b j a b ln ξ a b | g a 1 g b | ξ b a + ln ξ b a | g b 1 g a | ξ a b .
Here | ξ a b C 2 are coherent spinors associated with face normals [23,24], and g a SL ( 2 , C ) are group elements for each tetrahedron.
Variation of S CF with respect to g a yields the closure condition (3). At the Type-III saddle, the gluing element is a pure hyperbolic boost with rapidity χ a b 0 , and
S CF [ X ] = 2 a < b j a b χ a b = 16 π G γ S E ,
confirming that the saddle produces the suppression factor exp ( S E / ) .

Appendix A.2. Closure Conditions in the Complexified Domain

For Type-III boundary data, the saddle is reached by analytic continuation. The closure condition (3) in the complexified domain reads
b a j a b cosh ( χ a b ) n ^ a b ( a ) = 0 ,
where n ^ a b ( a ) S 2 are the real unit normals of the faces of tetrahedron τ a , and J a b j a b cosh ( χ a b ) are the effective areas. The self-consistency condition for the saddle is
χ a b = Θ a b ( E ) { j c d cosh χ c d } ,
where Θ a b ( E ) is the dihedral angle of the Euclidean 4-simplex with face areas { J c d } . Equations (A26) and (A27) must be solved simultaneously.

Appendix A.3. Reduction for the (3+2) Distribution

For the ( 3 + 2 ) distribution (tetrahedra τ 0 , τ 1 , τ 2 future-pointing, τ 3 , τ 4 past-pointing), the faces split under S 3 × S 2 symmetry into:
  • F–F: faces f 01 , f 02 , f 12 shared by two future tetrahedra, spin j S (3 faces);
  • F–P: faces f a b with a { 0 , 1 , 2 } , b { 3 , 4 } , spin j L (6 faces);
  • P–P: face f 34 shared by the two past tetrahedra, spin j S (1 face).

Closure of τ0.

Tetrahedron τ 0 has two F–F faces (effective area J FF = j S cosh χ FF ) and two F–P faces (effective area J FP = j L cosh χ FP ). At the solution we seek, the effective areas satisfy J FF = J FP (established below), so the tetrahedron with equal effective weights is regular and its face normals satisfy n ^ 0 i · n ^ 0 j = 1 3 for i j . Projecting (A26) onto n ^ 01 :
J FF 1 1 3 J FP 1 3 + 1 3 = 0 j S cosh χ FF = j L cosh χ FP .

Closure of τ3.

Tetrahedron τ 3 has three F–P faces and one P–P face (effective area J PP = j S cosh χ PP ). Projecting onto n ^ 03 :
J FP 1 2 3 J PP 1 3 = 0 j L cosh χ FP = j S cosh χ PP .

Consequence.

Equations (A28)–(A29) give cosh χ FF = cosh χ PP , hence
χ FF = χ PP χ S , j S cosh χ S = j L cosh χ L , χ L χ FP .

Appendix A.4. Exact Analytical Solution

From (A30), all effective areas are equal: J FF = J FP = J PP J . The Euclidean 4-simplex of effective areas is therefore regular. For the regular 4-simplex in R 4 , the outward unit normals satisfy n a · n b = 1 4 for a b , giving dihedral angle θ = arccos ( 1 4 )  [5]. Condition (A27) then yields the exact closed-form solution:
χ L = arccos 1 4 ,
χ S = arccosh ρ cosh arccos 1 4 ,
where ρ = j L / j S . The value χ L = arccos ( 1 / 4 ) is universal: it is independent of ρ and of the overall spin scale. Since ρ cosh ( arccos 1 4 ) cosh ( arccos 1 4 ) > 1 for all ρ 1 , equation (A32) defines a real positive χ S for all ρ 1 . There is no critical threshold ρ crit : the Type-III saddle exists for the ( 3 + 2 ) distribution at any ρ > 0 .
Remark A1.
Solving (A27) alone, without enforcing the closure constraints (A28)–(A29), yields a different numerical fixed point ( χ FF 1.056 , χ FP 1.562 , χ PP 0.240 ). That system does not satisfy the closure equations and is therefore not a true saddle of the vertex action.

Appendix A.5. Numerical Verification: jS = 1, jL = 3

For ρ = 3 , equations (A31)–(A32) give:
χ L = arccos 1 4 = 1.3181 rad , χ S = arccosh 3 cosh arccos 1 4 = 2.4789 rad .
Closure verification: j S cosh χ S = j L cosh χ L = 6.006 , residual < 10 12 . The geometric constant (15) with α FF = α PP = 1 , α FP = 3 :
C ( α ) = 4 χ S + 18 χ L = 4 × 2.4789 + 18 × 1.3181 = 33.64 .

Appendix A.6. Non-Existence of Real Saddles: Numerical Verification

Table A1. Minimum residual F min of the real closure equations under gradient descent from 10 3 random initialisations. A positive lower bound F min 0 confirms the absence of any real solution. The regular simplex ( ρ = 1 ) is Type II: it admits a complex solution with both real and imaginary parts but no real saddle, consistent with F min 1.0 > 0 .
Table A1. Minimum residual F min of the real closure equations under gradient descent from 10 3 random initialisations. A positive lower bound F min 0 confirms the absence of any real solution. The regular simplex ( ρ = 1 ) is Type II: it admits a complex solution with both real and imaginary parts but no real saddle, consistent with F min 1.0 > 0 .
Configuration j S j L F min
Uniform ( ρ = 1 , regular simplex, Type II) 1 1 1.0
Mixed ( ρ = 3 , ( 3 + 2 ) distribution, Type III) 1 3 4.0

Appendix B. Hessian Structure and Gauge Fixing

Appendix B.1. Parametrisation and Block Structure

Around the saddle X , we parametrise n a b = n a b + ϵ δ n a b and g a b = g a b exp ( ϵ δ ξ a b ) . The Hessian has block structure H = ( H n n , H n g ; H g n , H g g ) , with gauge block
( H g g ) a b , c d = j a b δ a c δ b d cosh χ a b P a b ,
where P a b projects orthogonally to the plane of n a b and n b a . The mixed block H n g vanishes at the saddle, making H block-diagonal.

Appendix B.2. Degree-of-Freedom Count and Gauge Fixing

In the Conrady–Freidel parametrisation, the integration variables are: 5 group elements g a SL ( 2 , C ) (real dimension 6 each) and, for each of the 10 pairs ( a , b ) , two coherent spinors | ξ a b , | ξ b a C P 1 S 2 (real dimension 2 each). Total: 5 × 6 + 10 × 2 × 2 = 70 real variables.
The gauge symmetries before fixing are: the diagonal action g a h g a of SL ( 2 , C ) ( dim = 6 ); the local SU ( 2 ) action per tetrahedron | ξ a b R a | ξ a b ( 5 × 3 = 15 ); and the individual U ( 1 ) phases of the spinors ( 2 × 10 = 20 ). Total gauge modes: 6 + 15 + 20 = 41 .
Physical degrees of freedom after complete gauge fixing:
n = 70 41 = 29 .
The prefactor therefore scales as ( det H res ) 1 / 2 j 29 / 2 , subdominant relative to the exponential factor for j 1 .

Appendix B.3. Non-Degeneracy at the Type-III Saddle

Proposition A1.
For the Type-III saddle with distribution ( k + ( 4 k ) ) , 1 k 3 , and non-degenerate geometry, dim ( Stab ( X ) ) = 0 , hence det H res 0 .
Proof. 
In the Type-I case, the stabiliser contains the diagonal SL ( 2 , C ) acting on all tetrahedra ( dim = 6 ). For the Type-III saddle, a transformation g SL ( 2 , C ) must simultaneously preserve the future and past normals. For k 2 and generic geometry, no g id has this property: Stab ( X ) is trivial. For k = 2 , there is a discrete Z 2 symmetry but no continuous one: dim ( Stab ( X ) ) = 0 in all cases.    □

Appendix C. Proof of the Refinement Stability Theorem

Appendix C.1. Normal Decomposition Under Stellar Subdivision

We establish a precise formula for the outward normals of the new tetrahedra introduced by a stellar subdivision. This formula is the key ingredient in the proof of causal-obstruction inheritance (Proposition A3 below).
Setup and notation. Let σ = Conv ( v 0 , v 1 , v 2 , v 3 , v 4 ) be a non-degenerate 4-simplex in R 3 , 1 with Minkowski metric η = diag ( 1 , + 1 , + 1 , + 1 ) . For each a { 0 , , 4 } , let N a R 3 , 1 be the outward unit normal to the boundary tetrahedron τ a = Conv ( { v b } b a ) in σ , normalised so that
η ( N a , N a ) = ϵ a { 1 , + 1 } ,
and satisfying the outward-pointing convention
η ( N a , v b ) = ϕ a for all b a , η ( N a , v a ) = ϕ a + d a , d a > 0 .
The constant ϕ a exists because N a is orthogonal to all edge vectors v b v c with b , c a , which span the tangent space of τ a ; the strict inequality d a > 0 expresses that v a lies strictly on the outward side of the hyperplane H a = { x : η ( N a , x ) = ϕ a } .
Let c = a = 0 4 λ a v a with λ a > 0 and a λ a = 1 be an interior vertex. The stellar subdivision F replaces σ by the five 4-simplices
σ i = Conv c , { v a } a i , i = 0 , , 4 .
For j i , define the new tetrahedron
τ ^ i j = Conv c , { v a } a i , a j ,
which is the face of σ i opposite to v j .
Lemma A1
(Normal decomposition). Under the conditions above, the outward normal to τ ^ i j in σ i is
N ^ i j = λ j d j N i λ i d i N j λ j d j N i λ i d i N j η ,
where | · | η denotes the Minkowski norm. The denominator is strictly positive for non-degenerate σ.
Proof. 
We seek a vector w = α i N i + α j N j normal to τ ^ i j and outward in σ i . The tangent space of τ ^ i j is spanned by { v k c : k i , k j } . Normality requires η ( w , v k c ) = 0 for all k i , j , i.e., η ( w , v k ) = η ( w , c ) .
Using condition (A38), η ( N i , v k ) = ϕ i for k i , and expanding
η ( N i , c ) = a λ a η ( N i , v a ) = ϕ i + λ i d i ,
the normality condition becomes
α i ϕ i + α j ϕ j = α i ( ϕ i + λ i d i ) + α j ( ϕ j + λ j d j ) ,
which simplifies to α i λ i d i + α j λ j d j = 0 , giving
α i α j = λ j d j λ i d i .
Choosing α i = λ j d j and α j = λ i d i yields w i j ( raw ) = λ j d j N i λ i d i N j .
Outward orientation. The vertex of σ i opposite to τ ^ i j is v j . From condition (A38) we have η ( N i , v j c ) = ϕ i ( ϕ i + λ i d i ) = λ i d i and η ( N j , v j c ) = ( ϕ j + d j ) ( ϕ j + λ j d j ) = ( 1 λ j ) d j > 0 (since λ j < 1 ). Substituting:
η w i j ( raw ) , v j c = λ j d j ( λ i d i ) λ i d i ( 1 λ j ) d j = λ i d i d j < 0 ,
since λ i , d i , d j > 0 , confirming outward orientation.
Non-degeneracy. Non-degeneracy of σ implies N i and N j are linearly independent; hence w i j ( raw ) 0 for λ i d i , λ j d j > 0 , and the normalised form (A41) is well-defined.    □
Remark A2.
When d i = d j (e.g. when all boundary tetrahedra are normalised with the same deficit parameter, as in a regular simplex), formula (A41) reduces to N ^ i j λ j N i λ i N j . In general, the d a -weighted version (A41) is required.

Appendix C.2. Inheritance of Causal Obstruction

Proposition A3
(Causal-obstruction inheritance). Let σ be a non-degenerate 4-simplex with Type-III boundary data of distribution ( k , 4 k ) , 1 k 3 : tetrahedra τ 0 , , τ k 1 are future-pointing ( N a C + ) and τ k , , τ 4 are past-pointing ( N a C ). Let c be any point in the relative interior of σ, and let F be the stellar subdivision of σ at c. Then every sub-simplex σ i of F has Type-III boundary data.
Proof. 
Fix any sub-simplex σ i . Its boundary consists of: the original tetrahedron τ i ; and the four new tetrahedra { τ ^ i j } j i . We must show that the collection of outward normals { N i } { N ^ i j } j i contains both future- and past-pointing elements, and that their distribution satisfies the conditions of Theorem 2.2.
Case 1: i < k (original tetrahedron τ i is future-pointing).
The normal N i C + contributes one future-pointing element. We analyse the normals N ^ i j for each j i .
Subcase 1a: j < k (future-pointing original tetrahedron). By Lemma A1, N ^ i j λ j d j N i λ i d i N j . Both N i , N j C + and λ j d j , λ i d i > 0 . The combination λ j d j N i λ i d i N j is the difference of two future-pointing vectors. Such a difference may be future-pointing, past-pointing, spacelike, or null depending on the geometry; it does not obstruct Type-III classification by itself.
Subcase 1b: j k (past-pointing original tetrahedron). By Lemma A1, N ^ i j λ j d j N i λ i d i N j . Since N j C , we have N j C + . Writing N ^ i j ( raw ) = λ j d j N i + λ i d i ( N j ) , both summands belong to C + with strictly positive coefficients. By Lemma 2.1 (convexity of C + ), the sum belongs to C + :
N ^ i j C + for all j k .
Establishing the Type-III distribution for σ i , i < k . The five boundary tetrahedra of σ i have outward normals:
  • N i C + (original, future-pointing);
  • N ^ i j C + for j k (from Subcase 1b, there are 4 k such faces);
  • N ^ i j for j < k , j i (from Subcase 1a, there are k 1 such faces, with mixed or indeterminate orientation).
We need to confirm the presence of at least one past-pointing normal. The faces from Subcase 1a are the candidates. Specifically, we claim that at least one N ^ i j with j < k is past-pointing for generic data.
Consider N ^ i j ( raw ) = λ j d j N i λ i d i N j for i , j < k . Both N i and N j are future-pointing. For the difference to be past-pointing, we need λ i d i N j to dominate λ j d j N i in the future-cone sense. This is precisely the complementary case to Lemma 2.1: since C + is a cone (not a subspace), the difference of two future-pointing vectors has no definite causal character, and its sign is controlled by the relative magnitudes λ i d i versus λ j d j .
To conclude rigorously, we use a counting argument. The closure condition (3) for σ i reads
N i + j i N ^ i j ( raw ) · ( positive weight ) = 0 .
The terms with j k all belong to C + , contributing a net future-pointing vector. The term N i C + also contributes to the future side. For the closure sum to vanish, the remaining terms (with j < k , j i ) must collectively lie in C , i.e., have a net past-pointing component. This is only possible if at least one N ^ i j with j < k is past-pointing. Hence the distribution of σ i contains both future- and past-pointing normals.
Since σ i has 4 k future-pointing tetrahedra from Subcase 1b, plus N i itself (future), and at least one past-pointing normal from Subcase 1a, the distribution is of the form ( k , 4 k ) with 1 k 3 . By Theorem 2.2, σ i has Type-III boundary data.
Case 2: i k (original tetrahedron τ i is past-pointing). The argument is identical under time-reversal symmetry: exchange C + C and i < k i k throughout Case 1.
Measure-zero exceptions. The only degenerate case is when both sides of (5) vanish simultaneously in σ i . By the same argument as in Theorem 2.2, this requires a co-dimension-one relation among the spins { j a b } , which holds on a set of measure zero in parameter space. For generic spins, the obstruction is strict.    □
Remark A4.
The counting argument in Case 1 above can be made fully explicit for the specific distributions ( k , 4 k ) :
  • ( 1 , 4 ) : σ i (with i = 0 , future) has N ^ 0 j C + for j = 1 , 2 , 3 , 4 (all past-pointing originals, Subcase 1b). Since N 0 C + and all N ^ 0 j C + , the only way closure can hold is if some combination is past-pointing — but by Theorem 2.2 applied to the full set, no real solution exists. A separate argument applies to the four sub-simplices with i 1 : each has N i C and N ^ i 0 C + (from Subcase 1b applied with roles reversed), giving a ( 1 , 4 ) distribution.
  • ( 2 , 3 ) and ( 3 , 2 ) : similar case analysis, omitted for brevity.
In all cases, the distribution inherited by σ i is of Type-III.
Corollary A5
(H2 is a theorem). Hypothesis (H2)of Theorem 4.3—that all internal vertices of the stellar subdivision have Type-III data—is automatically satisfied for any subdivision with interior central vertex c. It need not be assumed; it follows from Proposition A3.

Appendix C.3. Convergence of Saddle Angles

Let F be the stellar subdivision with N sub-simplices { σ i } satisfying (H1). The power N 1 / 4 in the convergence estimate reflects the scaling relation between edge length and volume in 4 dimensions: in a subdivision into N equal parts of volume Vol 4 / N , the typical edge length scales as ( Vol 4 / N ) 1 / 4 , producing an O ( N 1 / 4 ) correction to the shape coefficients α a b ( i ) relative to the original values α a b . By regularity of χ a b ( { α c d } ) :
χ a b ( i ) = χ a b + O ( N 1 / 4 ) .

Appendix C.4. Scaling of the Total Action and Prefactor

The total Euclidean action is S E tot = 1 8 π G i a < b A a b ( i ) χ a b ( i ) . Since i A a b ( i ) = A a b + O ( N 1 / 4 ) by area additivity and using (A44):
S E tot = S E 1 + O ( N 1 / 4 ) ,
establishing the result of Theorem 4.3 with exponent  1 / 4 . The prefactor N n / 2 arises from the N factors ( det H ( i ) ) 1 / 2 j 29 / 2 , subdominant relative to the exponential factor. □

Appendix D. Calibration of the Gram Matrix Criterion

Appendix D.1. Asymptotic Behaviour for Type-I Vertices

For Type-I vertices with a real Lorentzian saddle, the linear dependence of the weighted normals implies det ( G v ) = 0 exactly at the saddle. In the finite semiclassical regime, the departure from the saddle is exponentially suppressed in j.
Proposition A1
(Asymptotic scaling of the Gram determinant). Consider a vertex with boundary data in the semiclassical regime j 1 .
1.
For Type-I vertices (with a real Lorentzian saddle), the Gram matrix determinant satisfies
| det ( G v ) | = O e c j min ( v ) , c > 0 ,
reflecting the exponential approach to the exact closure condition.
2.
For Type-III vertices, the determinant is generically of order unity, | det ( G v ) | = O ( 1 ) .
Proof. 
For Type-I vertices, the existence of a real Lorentzian saddle implies exact linear dependence of the weighted normal vectors at the critical point. In the semiclassical regime, the dominant contribution to the path integral comes from a Gaussian neighbourhood of that saddle, with fluctuations suppressed as exp ( j δ 2 ) . Since the Gram determinant measures the deviation from exact linearity, its magnitude inherits this exponential suppression, giving | det ( G v ) | e c j for some geometry-dependent c > 0 .
For Type-III vertices, Theorem 2.2 guarantees the non-existence of any real configuration satisfying closure. The deviation from linear dependence is therefore finite and cannot be reduced by increasing j, implying | det ( G v ) | = O ( 1 ) generically.    □
Remark A2.
The above result is not an independent theorem but an asymptotic consequence of the saddle analysis, combined with the geometric interpretation of the Gram determinant as a measure of non-closure. Its validity is supported by the Gaussian structure of the semiclassical fluctuations and by the numerical verification of Appendix A.

Appendix D.2. Discrimination Window

For j min = 1 , the base threshold is α e 2 1.35 × 10 4 , well below the typical Type-I value | det ( G v ) | Type I e 2 0.135 and well above numerical noise ( 10 15 ). The factor ( 1 + β / d ( v , F ) ) with β 0.1 raises the threshold near the boundary, where partially fixed data may produce intermediate values of | det ( G v ) | .

References

  1. Rovelli, C. Quantum Gravity; Cambridge University Press: Cambridge, 2004. [Google Scholar]
  2. Thiemann, T. Modern Canonical Quantum General Relativity; Cambridge University Press: Cambridge, 2007. [Google Scholar]
  3. Perez, A. The spin-foam approach to quantum gravity. Living Rev. Relativity 2013, 16, 3. [Google Scholar] [CrossRef] [PubMed]
  4. Bonzom, V. Spin foam models for quantum gravity from lattice path integrals. Phys. Rev. D 2012, 85, 024028. [Google Scholar] [CrossRef]
  5. Barrett, J. W.; Dowdall, R. J.; Fairbairn, W. J.; Hellmann, F.; Pereira, R. Asymptotic analysis of the EPRL four-simplex amplitude. Class. Quantum Grav. 2009, 26, 165009. [Google Scholar]
  6. Conrady, F.; Freidel, L. On the semiclassical limit of spin foam models. Phys. Rev. D 2008, 78, 104023. [Google Scholar] [CrossRef]
  7. Engle, J.; Livine, E.; Pereira, R.; Rovelli, C. LQG vertex with finite Immirzi parameter. Nucl. Phys. B 2008, 799, 136. [Google Scholar] [CrossRef]
  8. Regge, T. General relativity without coordinates. Nuovo Cim. 1961, 19, 558. [Google Scholar] [CrossRef]
  9. Haggard, H. M.; Rovelli, C. Quantum-gravity effects outside the horizon spark black-to-white hole tunneling. Phys. Rev. D 2015, 92, 104020. [Google Scholar] [CrossRef]
  10. Bianchi, E.; Christodoulou, M.; Di Biagio, F. Spacetime tunneling and the quantum bounce of black holes. Class. Quantum Grav. 2023, 40, 085009. [Google Scholar]
  11. Christodoulou, M.; Di Biagio, F.; Bianchi, E. The quantum tunneling time of black holes. Class. Quantum Grav. 2024, 41, 075002. [Google Scholar]
  12. Donà, P.; Speziale, S.; Zhang, M. Asymptotics of the EPRL amplitude with timelike tetrahedra, Class. Quantum Grav. 2019, 36, 115011. [Google Scholar]
  13. Donà, P.; Gozzini, F.; Speziale, S. Numerical analysis of spin foam dynamics and the flatness problem. Phys. Rev. D 2020, 102, 106003. [Google Scholar] [CrossRef]
  14. Hellmann, F.; Kaminski, W. Geometric asymptotics for spin foam amplitudes. JHEP 2013, 165, 1310. [Google Scholar]
  15. Han, M.; Zhang, M. Asymptotics of spinfoam amplitude on simplicial manifold: Euclidean theory. Class. Quantum Grav. 2012, 29, 165004. [Google Scholar] [CrossRef]
  16. Han, M. On spinfoam models in large spin regime. Class. Quantum Grav. 2020, 37, 225011. [Google Scholar] [CrossRef]
  17. Witten, E. Analytic continuation of Chern-Simons theory. AMS/IP Stud. Adv. Math. 2011, arXiv:1001.293350, 347. [Google Scholar]
  18. Witten, E. A new look at the path integral of quantum mechanics. Surveys in Differential Geometry 2011, arXiv:1009.603215, 345. [Google Scholar] [CrossRef]
  19. Feldbrugge, J.; Lehners, J.-L.; Turok, N. No smooth beginning for spacetime. Phys. Rev. Lett. 2017, 119, 171301. [Google Scholar] [CrossRef]
  20. Feldbrugge, J.; Lehners, J.-L.; Turok, N. Lorentzian quantum cosmology. Phys. Rev. D 2017, 95, 103508. [Google Scholar] [CrossRef]
  21. Lehners, J.-L. Allowable histories in loop quantum cosmology. Phys. Rev. D 2017, 96, 026007. [Google Scholar]
  22. Conrady, F.; Freidel, L. Quantum geometry from phase space reduction. J. Math. Phys. 2009, 50, 123510. [Google Scholar] [CrossRef]
  23. Freidel, L.; Speziale, S. Twisted geometries: A geometric parametrisation of SU(2) phase space. Phys. Rev. D 2010, 82, 084040. [Google Scholar] [CrossRef]
  24. Dittrich, B.; Ryan, J. P. Phase space descriptions for simplicial 4d geometries. Class. Quantum Grav. 2011, 28, 065006. [Google Scholar] [CrossRef]
  25. Barrett, J. W.; Fairbairn, W. J.; Hellmann, F. Quantum gravity asymptotics from the SU(2) 15j symbol. Int. J. Mod. Phys. A 2010, 25, 2897. [Google Scholar] [CrossRef]
  26. Livine, E. R.; Oriti, D. Implementing causality in the spin foam quantum geometry. Phys. Rev. D 2002, 65, 044047. [Google Scholar] [CrossRef]
  27. Markopoulou, F.; Smolin, L. Causal evolution of spin networks. Nucl. Phys. B 1997, 508, 409. [Google Scholar] [CrossRef]
  28. Barrett, J. W.; Dowdall, R. J.; Fairbairn, W. J.; Gomes, H.; Hellmann, F. Lorentzian spin foam amplitudes: graphical calculus and asymptotics. Class. Quantum Grav. 2010, 27, 165009. [Google Scholar] [CrossRef]
  29. Freidel, L.; Krasnov, K. A new spin foam model for 4D gravity, Class. Quantum Grav. 2008, 25, 125018. [Google Scholar] [CrossRef]
  30. Rovelli, C.; Vidotto, F. Planck stars. Int. J. Mod. Phys. D 2014, 23, 1442026. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated