Representations of Quantum Gravity
Part IV · Transition Amplitudesgr-qc

Covariant Quantum Geometry as a State Sum: Spin Foams, Vertex Amplitudes, and the Refinement Problem of Boundary Amplitudes

Abstract

A spin foam is a two-complex whose faces and edges carry representation labels and whose state sum defines a transition amplitude between spin-network boundary states. We treat this construction as the covariant, sum-over-histories layer of a reconstruction program in which spacetime, causality, and measurement are read off invariant representation structures rather than quantized on a fixed background. The object of interest is not a metric but a triple: a boundary Hilbert space, an amplitude map defined by a colored two-complex, and the behavior of that map under refinement of the complex. We give a self-contained account of colored two-complexes dual to triangulated four-manifolds, the face/edge/vertex decomposition of the state sum Z(F)={jf},{ιe}fAfeAevAv\Zc(\Ff)=\sum_{\{j_f\},\{\iota_e\}}\prod_f A_f\,\prod_e A_e\,\prod_v A_v, and the constraint route from topological BFBF theory to gravity: strong simplicity in the Barrett–Crane model and weak (linear) simplicity through the Engle–Pereira–Rovelli–Livine Yγ\Ygam map in the model that carries the Barbero–Immirzi parameter. Our results are stated as labeled propositions with each claim tagged by the project's S/H/P epistemic calculus. That the truncated state sum is a well-defined finite amplitude map, that the topological BFBF sum is a triangulation invariant after quantum-group regularization, and that the Engle–Pereira–Rovelli–Livine four-simplex vertex has a large-spin asymptotic proportional to the cosine of the Regge action are standard mathematics ([S]). That the resulting vertex is the correct gravitational dynamics, and above all that the family of amplitudes over refinements has a continuum limit reproducing general relativity, are heuristic ([H]) and, in the case of the continuum limit, openly unresolved. We are explicit about the flatness problem, the accidental curvature constraint that obstructs a naive large-spin recovery of curved geometry, and about the effective-spin-foam reframing that partially addresses it. A worked example computes the vertex structure dual to a single four-simplex and states the asymptotic formula; a companion Haskell development exhibits the refinement invariant exactly in the solvable topological case and supplies checkable recoupling bookkeeping. We close by placing spin foams among the other regimes of the library: spin networks supply the boundary data, classical Lorentzian geometry is the target the missing continuum limit must reproduce, and the weakest-link status calculus bounds every such composite by its bare-[H]{} spin-foam link.

1 Introduction

1.1 The reconstruction thesis for covariant dynamics

The canonical picture of loop quantum gravity gives a kinematics: a state of quantum three-geometry is a spin network, a graph whose edges carry SU(2)\mathrm{SU}(2) representations and whose vertices carry intertwiners, and geometry appears as the discrete spectrum of area and volume operators built from those labels. That picture is the subject of a companion paper. What it does not give, on its own, is dynamics. A theory of gravity must say how one three-geometry evolves into another, and general covariance forbids doing so by a preferred time coordinate. The covariant answer is a sum over histories: instead of evolving a spin network by a Hamiltonian, one sums over all “quantum four-geometries” interpolating between a fixed incoming and outgoing spin network. A spin foam is precisely such a history, and the sum over them is the spin foam amplitude.

This paper works in that covariant layer, under the organizing principle shared across a modular library of twelve quantum-gravity regimes:

Quantum gravity is not primarily the quantization of material objects in spacetime. It is the reconstruction of spacetime, matter, causality, and measurement from invariant representation structures.

In the present regime the invariant structure is a colored two-complex. A history is a two-complex F\mathcal{F} whose faces are labeled by irreducible representations and whose edges are labeled by intertwiners; the amplitude of that history is a product of local recoupling invariants, one per face, edge, and vertex; and the physical content is the boundary amplitude obtained by summing over the bulk labels. Nowhere in this data is there a spacetime metric. A smooth four-geometry, with its Einstein–Hilbert or Regge action, is not put in by hand: it is recovered, if at all, only in an asymptotic (large-spin) limit of the representation-theoretic sum. Whether that recovery survives refinement of the two-complex into a genuine continuum theory is the central open problem of the field, and we treat it as open throughout.

1.2 What is settled and what is not

Honesty about status is the point of this series, so we state it at the outset. The following are theorems or standard constructions in a well-defined mathematical setting, and we tag them [S]:

  • For a finite two-complex with a cutoff on the representation labels, the state sum is a finite sum of products of finite-dimensional recoupling invariants, hence a well-defined amplitude map between boundary Hilbert spaces (5).

  • The topological BFBF state sum built from the {15j}\{15j\} symbol is invariant under Pachner moves once regularized by a quantum group, so it is a piecewise-linear invariant of the triangulated four-manifold (7).

  • Amplitudes glue along shared boundary spin networks, giving the state sum the composition law of a functor on a cobordism-like category (6).

  • The Engle–Pereira–Rovelli–Livine four-simplex vertex amplitude, evaluated on boundary data peaked on a nondegenerate Regge geometry, has a large-spin stationary-phase asymptotic whose leading term is proportional to cos(SRegge)\cos(S_{\mathrm{Regge}}) (12).

The following are physical hypotheses, well motivated but not proved, and we tag them [H]:

  • That the Engle–Pereira–Rovelli–Livine or Freidel–Krasnov vertex is the correct quantum dynamics of gravity. The vertex is fixed by imposing a discretized simplicity constraint, but the choice of how to impose it (strongly, weakly, on which variables), the face amplitude, and the measure are physical inputs, not consequences of an axiom.

  • That summing over refinements of the two-complex, or over two-complexes themselves, defines a continuum theory that reproduces general relativity. This is open. The known obstruction, the flatness problem (7.2), is that a naive large-spin limit at fixed complex is dominated by flat geometries.

The dictionary of the library assigns the entry spin foam two-complex with face/edge labels the bare status [H], the only bare-[H] entry among the core geometric regimes, precisely because the dynamics and its continuum limit remain a heuristic proposal rather than a settled structure. We honor that label: the abstract and this introduction state the continuum-limit status as unresolved, and every claim in the body carries its warrant explicitly.

1.3 Contributions and outline

We do not claim new physics. The contribution is a rigorous, status-labeled synthesis of the covariant-dynamics regime organized around one invariant — the boundary amplitude and its refinement behavior — together with a companion computation that exhibits that invariant exactly in the solvable topological case.

2 fixes the reconstruction pipeline and the status calculus. 3 defines colored two-complexes dual to triangulations, the face/edge/vertex data, and the state sum, and states the well-definedness and gluing results. 4 gives the constraint route from BFBF theory to gravity: the topological BFBF state sum and its invariance, the simplicity constraints, the Barrett–Crane model, and the YγY_{\gamma} map defining the Engle–Pereira–Rovelli–Livine vertex. 5 states the semiclassical asymptotics. 6 works the four-simplex example explicitly. 7 treats refinement, cylindrical consistency, the flatness problem, and the effective-spin-foam reframing honestly. 8 locates the regime among the others through the weakest-link calculus. 9 describes the Haskell model, whose core module is reproduced in 13. 10 lists the open problems, and [sec:discussion,sec:conclusion] conclude.

2 Mathematical framework: reconstruction from representation data

2.1 The realization pipeline

The library models every physical representation as a translation from a mathematical object MM to a physical target, factored through a realization: M    Φ(M)    Realα(Φ(M))    physical representation    Obs.M \;\longmapsto\; \Phi(M) \;\longmapsto\; \mathrm{Real}_\alpha(\Phi(M)) \;\longmapsto\; \text{physical representation} \;\longmapsto\; \mathrm{Obs}. For spin foams the mathematical object MM is a colored two-complex; the intermediate Φ(M)\Phi(M) is its state sum, a number attached to boundary data; the realization Realα\mathrm{Real}_\alpha chooses a model (which group, which vertex amplitude, which face amplitude and measure); and the observable is the boundary transition amplitude Z(F):HΓinHΓout\mathcal{Z}(\mathcal{F}):\mathcal{H}_{\Gamma_{\mathrm{in}}}\to\mathcal{H}_{\Gamma_{\mathrm{out}}}. The index α\alpha matters: Barrett–Crane, Engle–Pereira–Rovelli–Livine, and Freidel–Krasnov are different realizations of the same abstract state-sum object, and they differ in physically consequential ways. Keeping α\alpha explicit is what prevents the honest statement “the topological state sum is a triangulation invariant” from silently upgrading into the false statement “the gravitational state sum is a continuum theory.”

2.2 The status calculus

We use three ordered warrants S<H<P\mathsf{S} < \mathsf{H} < \mathsf{P}: S\mathsf{S} for a standard mathematical or mathematical-physics correspondence, H\mathsf{H} for a strong heuristic physical representation, and P\mathsf{P} for a speculative ontological extension. A status is a tuple rendered [S], [H], [P], or a composite such as [S/H]; composition takes the worst component.

Definition 1 (Weakest-link composition). For warrants a,ba,b with S<H<P\mathsf{S}<\mathsf{H}<\mathsf{P}, define ab=max(a,b)a\ast b=\max(a,b). A reasoning chain that composes representations of statuses a1,,ana_1,\dots,a_n carries status a1ana_1\ast\cdots\ast a_n. Composition is monotone and S\mathsf{S} is a unit: Sa=a\mathsf{S}\ast a=a.

The rule is not decoration. It is the precise sense in which the library is modular rather than unified: chaining spin-network kinematics (status [S/H]) to spin-foam dynamics (status [H]) yields a composite bounded above by [H], no matter how many rigorous theorems live inside the kinematical factor. One cannot buy credibility for the dynamics with the rigor of the kinematics.

2.3 The relevant dictionary entries

Three entries of the twenty-six-row dictionary bear on this paper.

Spin network: graph labeled by SU(2)\mathrm{SU}(2) representations == [S/H]. The boundary data. The kinematical Hilbert space and the area/volume spectra are [S] theorems; the physical identification with the geometry of our universe is [H].

Spin foam two-complex with face/edge labels == [H]. The subject of this paper. No assumptions field is populated; the limitation clause reads “continuum/ refinement limit and phenomenology are active research problems.”

Cobordism category == [S/H]. The categorical grammar in which a spin foam is a morphism from an incoming to an outgoing spin network, and gluing is composition.

The composite of the first two, under 1, is [H]: a complete covariant quantum-gravity dynamics built from spin-network kinematics plus spin-foam amplitudes inherits the bare-[H] label of the dynamics. This is exactly why the roadmap’s task for this module demands “record entries with explicit continuum-limit limitations” rather than asserting the package as settled.

2.4 Gauge and boundary as retained structure

A recurring structural fact across the library is that physical content is a groupoid, not a coarse quotient: gluing of local data is by equivalence, so automorphism and stabilizer data (gauge redundancy, symmetry) is retained rather than erased. In the spin-foam regime this appears twice. First, the edge intertwiners are exactly the gauge-invariant data at a node, the SU(2)\mathrm{SU}(2) (or SL(2,C)\mathrm{SL}(2,\mathbb{C})) invariant subspace of a tensor product of representations; the state sum integrates over the group at each edge, which is the discrete avatar of imposing the Gauss constraint. Second, the boundary of a spin foam is a spin network taken modulo spatial diffeomorphisms, and the amplitude is a pairing on that quotient. The reconstruction never works with a metric; it works with representation labels and invariants, and the geometry is a derived quantity.

3 Colored two-complexes and the state sum

3.1 Two-complexes dual to triangulations

Definition 2 (Two-complex). A (combinatorial, oriented) two-complex F\mathcal{F} is a set of vertices VV, edges EE, and faces FF, together with incidence data: each edge has a source and target vertex, and each face is bounded by a cyclic sequence of edges. We write fef\ni e when the edge ee lies on the boundary of the face ff, and eve\ni v when the vertex vv is an endpoint of ee. The boundary F\partial\mathcal{F} is the graph formed by the boundary edges and vertices (those lying on a single face on one side).

The spin foams relevant to four-dimensional gravity arise as the two-skeleton of the cell complex dual to a triangulation Δ\Delta of a four-manifold M\mathcal{M}. Under the duality, 4-simplices of Δ    vertices v,tetrahedra    edges e,triangles    faces f.\text{4-simplices of }\Delta \;\leftrightarrow\; \text{vertices }v,\quad \text{tetrahedra} \;\leftrightarrow\; \text{edges }e,\quad \text{triangles} \;\leftrightarrow\; \text{faces }f . A four-simplex has five tetrahedra and ten triangles on its boundary; dually, each vertex vv is met by five edges and ten faces, and the local combinatorics at vv is the complete graph K5K_5 whose five nodes are the foam edges and whose ten links are the foam faces. This K5K_5 pattern is what fixes the vertex amplitude to be a {15j}\{15j\}-type recoupling invariant. When M\mathcal M has boundary, F\partial\mathcal{F} is a pair of graphs Γin,Γout\Gamma_{\mathrm{in}},\Gamma_{\mathrm{out}}, the spin networks between which the foam interpolates.

The local structure at a spin-foam vertex vv dual to a four-simplex: five edges e1,,e5e_1,\dots,e_5 (dual to the five tetrahedra) meeting at vv, with the ten faces (dual to the ten triangles) running as the ten links of the complete graph K5K_5. The vertex amplitude contracts the five edge intertwiners along this graph, giving a {15j}\{15j\}-type invariant.

3.2 Face, edge, and vertex labels

Definition 3 (Coloring). A coloring of F\mathcal{F} assigns to each face ff an irreducible unitary representation and to each edge ee an intertwiner. In the Riemannian models the face label is an irreducible ρf\rho_f of SU(2)\mathrm{SU}(2) (equivalently a spin jf12Z0j_f\in\tfrac12\mathbb Z_{\ge0}, carrier VjfV_{j_f} of dimension 2jf+12j_f+1) or of Spin(4)=SU(2)×SU(2)\mathrm{Spin}(4)=\mathrm{SU}(2)\times\mathrm{SU}(2); in the Lorentzian models it is a unitary principal-series representation of SL(2,C)\mathrm{SL}(2,\mathbb{C}), labeled by a pair (ρf,kf)(\rho_f,k_f) with ρfR\rho_f\in\mathbb R and kf12Zk_f\in\tfrac12\mathbb Z. The edge label ιe\iota_e is an element of the invariant space Inv(feVjf)  =  G(C,  feVjf),G{SU(2),SL(2,C)},\mathrm{Inv}\Big(\bigotimes_{f\ni e} V_{j_f}\Big) \;=\;\mathop{\mathrm{Hom}}_{G}\Big(\mathbb C,\;\bigotimes_{f\ni e} V_{j_f}\Big), \qquad G\in\{\mathrm{SU}(2),\mathrm{SL}(2,\mathbb{C})\}, that is, an intertwiner among the representations on the faces meeting the edge.

The invariant space is finite-dimensional for SU(2)\mathrm{SU}(2); for SL(2,C)\mathrm{SL}(2,\mathbb{C}) it is infinite-dimensional, and the models handle this by an embedding (the YγY_{\gamma} map of 4.3) that lands the SU(2)\mathrm{SU}(2) intertwiner inside an SL(2,C)\mathrm{SL}(2,\mathbb{C}) one. On the boundary, the face and edge labels restrict to the edge and vertex labels of the boundary spin networks Γin,Γout\Gamma_{\mathrm{in}},\Gamma_{\mathrm{out}}: a boundary face is dual to a boundary spin-network edge, and a boundary edge is dual to a boundary node.

3.3 The state sum

Definition 4 (State-sum amplitude). Fix a model, i.e. a choice of face amplitude AfA_f, edge amplitude AeA_e, and vertex amplitude AvA_v. For a two-complex F\mathcal{F} with fixed boundary coloring b=(j,ι)b=(j_{\partial},\iota_{\partial}), the state-sum amplitude is Z(F)[b]  =  {jf},{ιe}  fFAf(jf)  eEAe(jfe,ιe)  vVAv(jfv,ιev),\begin{equation} \mathcal{Z}(\mathcal{F})[b] \;=\; \sum_{\{j_f\},\{\iota_e\}}\; \prod_{f\in F} A_f(j_f)\; \prod_{e\in E} A_e(j_{f\ni e},\iota_e)\; \prod_{v\in V} A_v\big(j_{f\ni v},\iota_{e\ni v}\big), \end{equation} where the sum runs over all colorings of the bulk (interior) faces and edges compatible with the fixed boundary bb, and the products run over faces, edges, and vertices. As bb varies, [eq:statesum] is a kernel that defines a linear map Z(F)  :  HΓinHΓout,\mathcal{Z}(\mathcal{F})\;:\;\mathcal{H}_{\Gamma_{\mathrm{in}}}\longrightarrow\mathcal{H}_{\Gamma_{\mathrm{out}}}, between the spin-network boundary Hilbert spaces.

The face amplitude is most often Af(jf)=Vjf=2jf+1A_f(j_f)=\mathop{\mathrm{dim}}V_{j_f}=2j_f+1 (the BFBF choice), though the correct power of the dimension is itself a model choice that affects the divergence structure and the semiclassical limit. The vertex amplitude AvA_v is the crucial dynamical input; the entire content of a spin-foam model is the choice of AvA_v, which is fixed by imposing a discretized simplicity constraint on a topological BFBF theory, as we now describe.

Theorem 5 (Well-definedness of the truncated amplitude, [S]). Let F\mathcal{F} be a finite two-complex and fix a cutoff Λ\Lambda on the representation labels (for SU(2)\mathrm{SU}(2): jfΛj_f\le\Lambda; for SL(2,C)\mathrm{SL}(2,\mathbb{C}): a compact set of principal-series labels). Then for any fixed boundary coloring bb, the truncated state sum ZΛ(F)[b]\mathcal{Z}_\Lambda(\mathcal{F})[b] obtained by restricting the sum in [eq:statesum] to labels below Λ\Lambda is a finite sum of finite products of finite-dimensional recoupling invariants, hence an absolutely convergent, well-defined complex number. The map ZΛ(F):HΓinHΓout\mathcal{Z}_\Lambda(\mathcal{F}):\mathcal{H}_{\Gamma_{\mathrm{in}}}\to\mathcal{H}_{\Gamma_{\mathrm{out}}} is linear and bounded on the finite-dimensional boundary label sectors.

Proof. With F\mathcal{F} finite there are finitely many faces, edges, and vertices. With the cutoff, each bulk face ranges over finitely many labels and each bulk edge over a finite-dimensional intertwiner space (finite-dimensional after the YγY_{\gamma} embedding in the SL(2,C)\mathrm{SL}(2,\mathbb{C}) case), so the index set of the sum is finite. Each factor Af,Ae,AvA_f,A_e,A_v is a contraction of finite-dimensional intertwiners — a Clebsch–Gordan or {nj}\{nj\} recoupling coefficient — and hence a finite complex number. A finite sum of finite products of finite numbers is finite. Linearity in the boundary data is immediate from [eq:statesum], which is multilinear in the boundary intertwiners; boundedness on a fixed finite-dimensional boundary sector is automatic. ◻

5 is deliberately modest. It says the amplitude is a well-defined number once a complex and a cutoff are fixed. It says nothing about removing the cutoff (the sum over bulk spins may diverge; see 7) or about refining the complex. Those are exactly the open questions, and separating the settled part from them is the discipline the status calculus enforces.

3.4 Gluing and the composition law

Proposition 6 (Gluing, [S/H]). Suppose F=F1ΓF2\mathcal{F}=\mathcal{F}_1\cup_\Gamma\mathcal{F}_2 is obtained by gluing two two-complexes along a common boundary spin-network graph Γ\Gamma, so that F1=ΓinΓ\partial\mathcal{F}_1=\Gamma_{\mathrm{in}}\sqcup\Gamma and F2=ΓΓout\partial\mathcal{F}_2=\Gamma\sqcup\Gamma_{\mathrm{out}}. Then, with the natural face and edge amplitudes on the shared boundary distributed to avoid double counting, Z(F)=c(Γ)Z(F1)[,c(Γ)]  Z(F2)[c(Γ),],\begin{equation} \mathcal{Z}(\mathcal{F}) =\sum_{c(\Gamma)} \mathcal{Z}(\mathcal{F}_1)\big[\,\cdot\,,c(\Gamma)\big]\; \mathcal{Z}(\mathcal{F}_2)\big[c(\Gamma),\,\cdot\,\big], \end{equation} where the sum runs over colorings c(Γ)c(\Gamma) of the shared boundary. Equivalently, Z(F)=Z(F2)Z(F1)\mathcal{Z}(\mathcal{F})=\mathcal{Z}(\mathcal{F}_2)\circ\mathcal{Z}(\mathcal{F}_1) as maps of boundary Hilbert spaces. The assignment ΓHΓ\Gamma\mapsto\mathcal{H}_{\Gamma}, FZ(F)\mathcal{F}\mapsto\mathcal{Z}(\mathcal{F}) is a functor from a cobordism-like category of spin networks and two-complexes to (finite-dimensional) vector spaces.

Proof sketch. The state sum [eq:statesum] factorizes over the two pieces: every bulk face, edge, and vertex belongs to exactly one of F1,F2\mathcal{F}_1,\mathcal{F}_2, while the shared faces and edges lie on Γ\Gamma. Fixing the shared coloring c(Γ)c(\Gamma) makes each piece an independent state sum with Γ\Gamma as (part of its) boundary, and summing over c(Γ)c(\Gamma) reconstitutes the sum over the interior labels of F\mathcal{F} that lie on the gluing surface. Reindexing gives [eq:glue]. The tag is [S/H] rather than [S]: for the topological BFBF models this is a genuine gluing axiom of an extended topological field theory ([S]), but for the gravitational models the boundary Hilbert spaces, the correct distribution of face amplitudes across the seam, and the convergence of the sum over c(Γ)c(\Gamma) are model-dependent, which is the heuristic part. ◻

6 is the structural reason spin foams sit in the cobordism/TQFT column of the dictionary. A spin foam is a morphism ΓinΓout\Gamma_{\mathrm{in}}\to \Gamma_{\mathrm{out}}; composition is gluing; the identity is a trivial cylinder. The qualifier “like” matters: a genuine TQFT would assign a finite-dimensional space to each boundary and be diffeomorphism (indeed triangulation) invariant, and the gravitational models are not triangulation invariant — if they were, they would be topological and contain no local degrees of freedom, hence no gravitons. The tension between wanting a triangulation invariant (for a background-independent continuum limit) and wanting local curvature degrees of freedom is the conceptual core of the refinement problem (7).

4 From BFBF theory to gravity: simplicity constraints

4.1 BFBF theory and its state sum

The starting point is a topological field theory. Let GG be a Lie group with Lie algebra g\mathfrak g. BFBF theory has fields a connection AA (with curvature F(A)F(A)) and a g\mathfrak g-valued two-form BB, and action SBF[A,B]=M(BF(A)).\begin{equation} S_{BF}[A,B]=\int_{\mathcal M}\mathop{\mathrm{tr}}\big(B\wedge F(A)\big). \end{equation} The equations of motion are F(A)=0F(A)=0 (flat connection) and dAB=0\mathrm{d}_A B=0; the theory has no local degrees of freedom. Discretizing on the two-complex F\mathcal{F}, one assigns a group element (holonomy) geGg_e\in G to each edge and integrates the constraint F(A)=0F(A)=0 face by face. Expanding the group delta function on the holonomy gf=efgeg_f=\prod_{e\subset\partial f} g_e around each face by the Peter–Weyl theorem, δG(gf)=jf(Vjf)χjf(gf),\begin{equation} \delta_G(g_f)=\sum_{j_f}\mathop{\mathrm{dim}}(V_{j_f})\,\chi_{j_f}(g_f), \end{equation} and performing the edge integrals, one obtains a state sum in which each face carries an irreducible jfj_f, each edge carries an intertwiner produced by the Haar projection, and each vertex carries the recoupling invariant contracting the intertwiners on its edges. For G=SU(2)G=\mathrm{SU}(2) on a four-dimensional complex the vertex invariant is the {15j}\{15j\} symbol, and the partition function is the Ooguri model, ZBF(F)={jf}  f(Vjf)  v{15j}(jfv).\begin{equation} \mathcal{Z}_{BF}(\mathcal{F})=\sum_{\{j_f\}}\;\prod_{f}\mathop{\mathrm{dim}}(V_{j_f})\;\prod_{v}\{15j\}(j_{f\ni v}). \end{equation}

Theorem 7 (Triangulation invariance of the topological state sum, [S]). The three-dimensional SU(2)\mathrm{SU}(2) state sum [eq:peterweyl][eq:ooguri] (Ponzano–Regge in dimension three, Ooguri in dimension four) is formally invariant under the Pachner moves relating any two triangulations of a fixed piecewise-linear manifold. After regularization by replacing SU(2)\mathrm{SU}(2) with the quantum group Uq(su2)U_q(\mathfrak{su}_2) at a root of unity (Turaev–Viro in dimension three, Crane–Yetter in dimension four), the resulting sum is finite and is a genuine piecewise-linear invariant of the manifold.

Proof sketch. Pachner-move invariance reduces, via the Biedenharn–Elliott (pentagon) identity and the orthogonality of 6j6j (respectively 15j15j) symbols, to algebraic identities among recoupling coefficients; these are the defining coherence relations of the representation category. The classical sums diverge because the trivial representation appears with infinite multiplicity in an unbounded sum over spins; the quantum-group truncation at level kk makes the set of admissible labels finite, and Turaev–Viro (dimension three) and Crane–Yetter (dimension four) proved the truncated sum converges and is move-invariant, hence a topological invariant. ◻

7 is the sharpest rigorous statement in the vicinity: a state sum of representation labels reconstructs a topological invariant of the manifold, exactly. It is also a warning. A theory that is fully triangulation invariant is topological — it has no local degrees of freedom and cannot be gravity. Gravity must break the topological invariance in a controlled way, keeping local curvature while retaining enough structure to have a continuum limit. That controlled breaking is the job of the simplicity constraint.

4.2 Simplicity constraints

General relativity in the Plebanski (or Palatini–Cartan) formulation is BFBF theory for G=SO(4)G=\mathrm{SO}(4) (Riemannian) or SO(3,1)\mathrm{SO}(3,1) (Lorentzian) supplemented by a constraint forcing the two-form BB to come from a tetrad: BIJ=±(ee)IJB^{IJ}=\pm(e\wedge e)^{IJ} or its dual. This simplicity constraint is what turns the topological BFBF theory into gravity with local degrees of freedom. In terms of the bivector BfIJB^{IJ}_f attached to a face (a discretized triangle), simplicity requires BfB_f to be a simple bivector, one of the form uwu\wedge w for vectors u,wu,w.

For Spin(4)=SU(2)×SU(2)\mathrm{Spin}(4)=\mathrm{SU}(2)\times\mathrm{SU}(2), the universal (double) cover of SO(4)\mathrm{SO}(4), a bivector splits into self-dual and anti-self-dual parts, carrying spins (jf+,jf)(j^+_f,j^-_f); the two SU(2)\mathrm{SU}(2) factors are exactly the self-dual and anti-self-dual sectors of the covered rotation group. Discrete simplicity, in its strong (Barrett–Crane) form, sets jf+=jf,\begin{equation} j^+_f=j^-_f, \end{equation} the balanced or simple representations. Geometrically, [eq:bcsimplicity] says the self-dual and anti-self-dual areas agree, which is the statement that the bivector is simple. The Barrett–Crane model imposes [eq:bcsimplicity] strongly, as an operator equation annihilating states, and builds the vertex from the SO(4)\mathrm{SO}(4)-invariant contraction of balanced representations.

Remark 8 (Why strong imposition is too strong, [H]). Imposing all simplicity constraints strongly, as operator equations, over-constrains the intertwiner degrees of freedom: the Barrett–Crane vertex fixes the edge intertwiner uniquely (the Barrett–Crane intertwiner), removing the intertwiner labels that should survive as independent quantum numbers. The consequence, recognized around 2007, is that the boundary states of the Barrett–Crane model do not match the spin-network states of canonical loop quantum gravity, whose nodes carry a full intertwiner space. This mismatch motivated the weak imposition of the Engle–Pereira–Rovelli–Livine and Freidel–Krasnov models. The status here is [H]: the mismatch is a physical-adequacy judgment about boundary data, not a theorem that Barrett–Crane is mathematically inconsistent.

4.3 The Engle–Pereira–Rovelli–Livine vertex and the YγY_{\gamma} map

The Engle–Pereira–Rovelli–Livine (EPRL) model imposes simplicity weakly, in the sense of the expectation value of a linear constraint, and in doing so incorporates the Barbero–Immirzi parameter γ\gamma. The construction rests on a map embedding SU(2)\mathrm{SU}(2) representations into SL(2,C)\mathrm{SL}(2,\mathbb{C}) (Lorentzian) or Spin(4)\mathrm{Spin}(4) (Riemannian) ones.

In the Lorentzian case the relevant SL(2,C)\mathrm{SL}(2,\mathbb{C}) representations are the unitary principal series H(ρ,k)\mathcal H_{(\rho,k)}, labeled by ρR\rho\in\mathbb R and k12Zk\in\tfrac12\mathbb Z. Each such representation decomposes under the SU(2)\mathrm{SU}(2) subgroup as H(ρ,k)=jkVj\mathcal H_{(\rho,k)}=\bigoplus_{j\ge k} V_j. The linear simplicity constraint, imposed weakly, selects the lowest SU(2)\mathrm{SU}(2) component and ties the labels through γ\gamma: ρf=γjf,kf=jf.\begin{equation} \rho_f=\gamma\, j_f,\qquad k_f=j_f . \end{equation}

Definition 9 (The YγY_{\gamma} map). The Engle–Pereira–Rovelli–Livine map is the isometric embedding Yγ:  Vj    H(γj,j),Yγj,m=(γj,j);j,m,\begin{equation} Y_{\gamma}:\;V_j\;\hookrightarrow\;\mathcal H_{(\gamma j,\,j)},\qquad Y_{\gamma}\,\lvert j,m \rangle=\lvert (\gamma j, j);\,j,m \rangle, \end{equation} sending the spin-jj SU(2)\mathrm{SU}(2) representation to the lowest (j=k=jj'=k=j) SU(2)\mathrm{SU}(2)-isotypic component of the principal series H(γj,j)\mathcal H_{(\gamma j,j)}. In the Riemannian model YγY_{\gamma} sends VjV_j into Vj+VjV_{j^+}\otimes V_{j^-} with j±=121±γjj^\pm=\tfrac{1}{2}|1\pm\gamma|\,j (assuming γ<1\gamma<1; the general case carries signs).

Definition 10 (EPRL vertex amplitude). For a vertex vv dual to a four-simplex, with the ten faces carrying spins jfj_{f} and the five edges a=1,,5a=1,\dots,5 carrying SU(2)\mathrm{SU}(2) intertwiners ιa\iota_a, the EPRL vertex amplitude is AvEPRL(jf,ιa)=SL(2,C)4a=14dga  TrK5 ⁣(a=15gaYγιa),\begin{equation} A_v^{\mathrm{EPRL}}\big(j_f,\iota_a\big) =\int_{\mathrm{SL}(2,\mathbb{C})^{4}} \prod_{a=1}^{4}\mathrm{d}g_a\; \operatorname{Tr}_{K_5}\!\bigg(\bigotimes_{a=1}^{5}\, g_a\triangleright Y_{\gamma}\iota_a\bigg), \end{equation} where nominally one SL(2,C)\mathrm{SL}(2,\mathbb{C}) element is associated to each of the five edges, but one group element (say g5=1g_5=\mathbf{1}) is held fixed to make the noncompact SL(2,C)\mathrm{SL}(2,\mathbb{C}) integral finite, leaving the four integrations over SL(2,C)4\mathrm{SL}(2,\mathbb{C})^{4} shown. Here gaYγιag_a\triangleright Y_{\gamma}\iota_a denotes the group element gag_a acting, in the representations (γjf,jf)(\gamma j_f,j_f) carried by the faces faf\ni a, on the YγY_{\gamma}-embedded edge intertwiner, and TrK5\operatorname{Tr}_{K_5} contracts the five resulting invariant tensors along the ten faces of the vertex — the ten edges of the complete graph K5K_5 of 1. The contraction is holistic and does not factorize into per-face matrix elements: an intertwiner ιaInv(faVjf)\iota_a\in\mathrm{Inv}\big(\bigotimes_{f\ni a}V_{j_f}\big) is an entangled invariant tensor, so a single face carries no gauge-invariant matrix element of its own; only the full K5K_5 trace is well defined. In the Riemannian model SL(2,C)\mathrm{SL}(2,\mathbb{C}) is replaced by Spin(4)\mathrm{Spin}(4) and the integrals are over a compact group, so no element need be fixed.

The vertex [eq:eprlvertex] is the recoupling invariant of the K5K_5 pattern of 1: five edges, ten faces, contracted along the complete graph. In the Riemannian case it reduces to a γ\gamma-weighted {15j}\{15j\} symbol. The Freidel–Krasnov model reaches an amplitude of the same form through a coherent-state (rather than projector) route; for γ<1\gamma<1 the two agree, and we treat them together as the weak-simplicity vertex.

Remark 11 (What is chosen, [H]). Three things in [eq:eprlvertex] are physical choices, not derivations: the weak (rather than strong) imposition of simplicity; the identification [eq:eprllabels] of the continuous SL(2,C)\mathrm{SL}(2,\mathbb{C}) label ρ\rho with γ\gamma times the discrete spin; and the face amplitude AfA_f, taken here as a power of Vjf\mathop{\mathrm{dim}}V_{j_f}. Different choices give different models with different divergence and continuum behavior. The literature motivates [eq:eprllabels] by matching the loop-gravity area spectrum AfγjfA_f\propto \gamma\, j_f, which is itself a semiclassical input. We tag the identification of [eq:eprlvertex] as “the” gravitational vertex with status [H].

5 Semiclassical asymptotics

The strongest evidence that the EPRL vertex is a gravity amplitude is its large-spin limit. Uniformly rescale all ten spins jfλjfj_f\mapsto\lambda j_f and let λ\lambda\to\infty.

Proposition 12 (Regge asymptotics of the EPRL four-simplex vertex, [S]). Let the boundary data of a single four-simplex vertex be a coherent (Livine–Speziale) spin-network state peaked on the geometry of a nondegenerate Euclidean or Lorentzian four-simplex with triangle areas Af=γjfA_f=\gamma\,j_f (in Planck units). Then as λ\lambda\to\infty the rescaled vertex amplitude has the stationary-phase asymptotic AvEPRL(λjf)    1λ12[N+eiλSRegge(jf)+NeiλSRegge(jf)]  +  (degenerate terms),\begin{equation} A_v^{\mathrm{EPRL}}(\lambda j_f)\;\sim\; \frac{1}{\lambda^{12}}\Big[\,N_+\,e^{\,i\lambda\,S_{\mathrm{Regge}}(j_f)} +N_-\,e^{-i\lambda\,S_{\mathrm{Regge}}(j_f)}\,\Big]\;+\;(\text{degenerate terms}), \end{equation} where SRegge(jf)=fvAfΘfS_{\mathrm{Regge}}(j_f)=\sum_{f\ni v} A_f\,\Theta_f is the (single-simplex) Regge action of the four-simplex, Θf\Theta_f the interior dihedral angle at the triangle ff (a single simplex has no deficit angle), and N±N_\pm are λ\lambda-independent prefactors. The leading term is proportional to cos(SRegge+const)\cos\big(S_{\mathrm{Regge}}+\text{const}\big).

Status and provenance. This is the theorem of Barrett, Dowdall, Fairbairn, Gomes, and Hellmann (Riemannian EPRL) and Barrett, Dowdall, Fairbairn, Hellmann, and Pereira (Lorentzian EPRL). The proof is a stationary-phase analysis: the vertex [eq:eprlvertex] is an integral over SL(2,C)\mathrm{SL}(2,\mathbb{C}) (or Spin(4)\mathrm{Spin}(4)) of a product of matrix elements written, on coherent boundary data, as eS(g)e^{S(g)} with SS an “action” on group variables; the critical points of SS are in bijection with reconstructions of a geometric four-simplex from the boundary areas and normals, and the on-shell value of SS at a nondegenerate critical point is ±iSRegge\pm i\,S_{\mathrm{Regge}}. Standard stationary-phase then gives [eq:asymp]: the integration runs over four copies of the six-dimensional group SL(2,C)\mathrm{SL}(2,\mathbb{C}) (or Spin(4)\mathrm{Spin}(4)), so the Hessian at a nondegenerate critical point is 2424-dimensional, and the stationary-phase prefactor is λ24/2=λ12\lambda^{-24/2}=\lambda^{-12}. The result is [S]: a proved asymptotic in a well-defined integral-asymptotics setting. ◻

Two honest caveats accompany 12, both [H]. First, the two terms of the cosine are the two parity-related geometric reconstructions of the four-simplex from the boundary areas and normals: in the Lorentzian model these are the two time orientations, so the amplitude does not distinguish a geometry from its time-reverse, and extracting a causal, single-exponential propagator requires additional input. The Riemannian model is more delicate still: its stationary-phase locus contains, besides the two geometric (Regge) critical points, additional degenerate “vector geometry” configurations that contribute non-Regge terms unless the boundary coherent state is chosen to suppress them, so [eq:asymp] is the leading behavior only on suitably peaked geometric boundary data. Second, and more consequentially, [eq:asymp] is the asymptotics of a single vertex on fixed boundary data. It does not by itself say what happens when one sums over the bulk spins of a many-vertex complex, and it is exactly there that the flatness problem appears (7.2).

6 Worked example: the four-simplex amplitude

We make the smallest nontrivial spin foam explicit: the two-complex dual to a single four-simplex, which has one bulk vertex, five edges (dual to the five boundary tetrahedra), and ten faces (dual to the ten boundary triangles). Its boundary is the spin network on the pentagram graph K5K_5: five nodes (tetrahedra), ten links (triangles).

6.1 Boundary data and the amplitude

Label the ten faces by spins jab=jbaj_{ab}=j_{ba} for 1a<b51\le a<b\le5 (the triangle shared by tetrahedra aa and bb), and the five edges by intertwiners ιaInv(baVjab)\iota_a\in\mathrm{Inv}\big(\bigotimes_{b\ne a}V_{j_{ab}}\big), each a four-valent intertwiner. The vertex amplitude [eq:eprlvertex] contracts the five intertwiners along K5K_5. In the Riemannian, γ0\gamma\to 0 or topological limit it degenerates to the {15j}\{15j\} symbol Av(jab,ιa)={15j}({jab};{ιa}),\begin{equation} A_v(j_{ab},\iota_a)=\{15j\}\big(\{j_{ab}\};\{\iota_a\}\big), \end{equation} the unique (up to normalization) SU(2)\mathrm{SU}(2)-invariant contraction of five four-valent intertwiners whose external legs are the ten jabj_{ab}. The full state sum for this one-vertex complex with fixed boundary is just the single vertex weighted by its face and edge amplitudes: Z(FΔ4)[b]=a<bAf(jab)  aAe(ιa)  Av(jab,ιa).\begin{equation} \mathcal{Z}(\mathcal{F}_{\Delta_4})[b] =\prod_{a<b}A_f(j_{ab})\;\prod_{a}A_e(\iota_a)\;A_v(j_{ab},\iota_a). \end{equation} There is no bulk sum for a single four-simplex with fully specified boundary: every face and edge is a boundary object. The sum over bulk labels first appears when two or more four-simplices are glued and share an internal triangle or tetrahedron.

6.2 Equilateral boundary and the asymptotic angle

Take the symmetric boundary jab=jj_{ab}=j for all ten faces, the quantum analog of an equilateral four-simplex. By 12, as jj\to\infty, Av(λj)    1λ122Ncos ⁣(λSRegge(j)+φ),SRegge(j)=10γjΘ,\begin{equation} A_v(\lambda j)\;\sim\;\frac{1}{\lambda^{12}}\,2\,|N|\, \cos\!\Big(\lambda\,S_{\mathrm{Regge}}(j)+\varphi\Big),\qquad S_{\mathrm{Regge}}(j)=10\,\gamma\, j\,\Theta, \end{equation} where Θ=arccos(1/4)\Theta=\arccos(1/4) is the interior dihedral angle of a regular Euclidean four-simplex (the angle between two of its tetrahedral faces), and the total is ten triangles times the area γj\gamma j times that dihedral angle. The amplitude oscillates in jj with a frequency set by the Regge action; the appearance of the regular four-simplex dihedral angle is the sense in which “the representation labels know the geometry.”

6.3 A solvable topological shadow of the refinement invariant

The four-simplex amplitude [eq:15j] is not something we can sum in closed form over refinements; that is the open problem. To exhibit the refinement invariant exactly, we drop to the solvable case: a finite-group BFBF (Dijkgraaf–Witten) state sum, the abelian, discrete-group avatar of the topological theory of 4.1. Let GG be a finite abelian group. On a closed oriented surface Σ\Sigma triangulated with VV vertices, EE edges, and FF triangles, assign a group element heGh_e\in G to each edge and define ZG(Σ)=1GV{he}  f[ef±he=0],\begin{equation} \mathcal{Z}_G(\Sigma)=\frac{1}{|G|^{V}} \sum_{\{h_e\}}\;\prod_{f}\big[\textstyle\sum_{e\subset\partial f}\pm h_e=0\big], \end{equation} where [][\cdot] is 11 if the oriented holonomy around the face is trivial and 00 otherwise, and the normalization GV|G|^{-V} divides out the gauge redundancy hehe+(ϕt(e)ϕs(e))h_e\mapsto h_e+(\phi_{t(e)}-\phi_{s(e)}). The following is a clean, fully [S] statement.

Proposition 13 (Exact refinement invariance in the solvable case, [S]). For finite abelian GG, the state sum [eq:dw] depends only on the topology of Σ\Sigma: ZG(Σ)=G1χ(Σ),\begin{equation} \mathcal{Z}_G(\Sigma)=|G|^{\,1-\chi(\Sigma)}, \end{equation} with χ(Σ)=VE+F\chi(\Sigma)=V-E+F the Euler characteristic. In particular it is invariant under any subdivision (refinement) of the triangulation: ZG(S2)=G1\mathcal{Z}_G(S^2)=|G|^{-1} for every triangulation of the sphere, and ZG(T2)=G\mathcal{Z}_G(T^2)=|G| for every triangulation of the torus.

Proof. The face constraints say the coloring is a GG-valued discrete flat connection, i.e. a cocycle in Z1(Σ;G)Z^1(\Sigma;G); the number of them is Z1|Z^1|. The gauge redundancy is B0B^0-worth of coboundaries, and for a connected Σ\Sigma the map C0Z1C^0\to Z^1 has kernel the constants, so the count of gauge orbits is Z1/GV1|Z^1|/|G|^{V-1}. Hence ZG=Z1/GV\mathcal{Z}_G=|Z^1|/|G|^{V}. For abelian GG, Z1=H1B1=G2gGV1|Z^1|=|H^1|\cdot|B^1|=|G|^{\,2g}\cdot |G|^{\,V-1} on a genus-gg surface (using dimB1=V1\dim B^1=V-1 for connected Σ\Sigma and H1=2g|H^1|=2g-dimensional). The count H1=G2g|H^1|=|G|^{2g} is the topological input: a flat GG-connection up to gauge is a homomorphism π1(Σg)G\pi_1(\Sigma_g)\to G, and since π1(Σg)\pi_1(\Sigma_g) has 2g2g generators subject to the single relation [a1,b1][ag,bg]=1[a_1,b_1]\cdots[a_g,b_g]=1, which is automatically satisfied once GG is abelian, one has H1=(π1(Σg),G)=G2g|H^1|=|\mathop{\mathrm{Hom}}(\pi_1(\Sigma_g),G)|=|G|^{2g}. Therefore ZG=G2g+V1/GV=G2g1=G1χ\mathcal{Z}_G=|G|^{2g+V-1}/|G|^{V}=|G|^{2g-1}=|G|^{1-\chi}, since χ=22g\chi=2-2g. Because χ\chi is a topological invariant, the value is unchanged under refinement. ◻

13 is the exactly-solvable shadow of what the gravitational state sum is being asked to do: produce an amplitude that stabilizes under refinement. In the topological case it stabilizes perfectly, because the theory has no local degrees of freedom. Gravity must stabilize while carrying local curvature, and whether it can is open. The companion Haskell code (9, 13) computes [eq:dw] by brute force for G=ZnG=\mathbb Z_n on several triangulations and checks [eq:dwvalue] and refinement invariance, alongside the SU(2)\mathrm{SU}(2) recoupling bookkeeping that supports [eq:15j].

7 Refinement, the continuum limit, and the flatness problem

7.1 Refinement and cylindrical consistency

A single two-complex is a truncation of the degrees of freedom, like a single lattice in lattice gauge theory. Because gravity is background independent, there is no fixed lattice spacing to send to zero; refinement means subdividing the complex, and the continuum theory (if it exists) is a limit over the directed set of refinements.

Definition 14 (Refinement and cylindrical consistency). A refinement FF\mathcal{F}\to\mathcal{F}' is a subdivision: each cell of F\mathcal{F} is partitioned into cells of F\mathcal{F}', inducing an embedding of boundary graphs ΓΓ\Gamma\hookrightarrow\Gamma' and hence an isometry ιΓΓ:HΓHΓ\iota_{\Gamma\Gamma'}:\mathcal{H}_{\Gamma}\hookrightarrow\mathcal{H}_{\Gamma'} of boundary Hilbert spaces. The family of amplitudes {Z(F)}\{\mathcal{Z}(\mathcal{F})\} is cylindrically consistent if for every refinement, Z(F)=ιΓΓZ(F)ιΓΓ,\begin{equation} \mathcal{Z}(\mathcal{F})=\iota_{\Gamma\Gamma'}^{\dagger}\,\mathcal{Z}(\mathcal{F}')\,\iota_{\Gamma\Gamma'}, \end{equation} i.e. coarser amplitudes are restrictions of finer ones. When [eq:cylconsist] holds, the amplitudes define a single map on the projective limit HΓ=HΓ\mathcal{H}_{\Gamma_\infty}=\varinjlim\mathcal{H}_{\Gamma}, a genuine continuum amplitude.

Remark 15 (The continuum limit is open, [H]). Generic spin-foam models are not cylindrically consistent as written: refining the complex changes the amplitude. The proposed remedies are (i) to construct improved, consistent amplitudes as fixed points of a coarse-graining (renormalization group) flow on the space of amplitudes, following Dittrich and collaborators, or (ii) to sum over all two-complexes, e.g. via a group field theory whose Feynman expansion generates spin foams. Neither program has produced a proof that a continuum limit exists and reproduces general relativity. We state this as the primary open problem of the regime and tag every claim about the existence of the limit [H].

7.2 The flatness problem

The most concrete obstruction to a naive continuum limit is the flatness (or “cosine,” or “accidental curvature”) problem. Consider a complex with internal faces and take the large-spin limit while summing over the bulk spins. On coherent boundary data each vertex contributes e±iS(v)e^{\pm i S^{(v)}} by 12, and the bulk amplitude becomes a stationary-phase integral over the group elements geg_e on the internal edges (together with the bulk spin sums). The critical-point equations of these group integrals are the geometric content of the flatness problem: at each internal face ff they force the product of edge holonomies around ff to be trivial, Hf=efge=1,\begin{equation} H_f=\prod_{e\subset\partial f} g_e=\mathbf{1}, \end{equation} and a trivial holonomy around every internal hinge is exactly the statement that the reconstructed geometry has vanishing deficit angle there, εf=2πvfθf(v)=0\varepsilon_f=2\pi-\sum_{v\supset f}\theta_f^{(v)}=0, with θf(v)\theta_f^{(v)} the interior dihedral angle contributed by the simplex at vertex vv. That is, the dominant configurations in the summed large-spin limit are flat. It is the group-integral saddle [eq:flatness], not the variation of the phase with respect to the areas Af=γjfA_f=\gamma j_f, that carries the obstruction: for a generic Barbero–Immirzi parameter the naive area-variation condition holds only modulo 2π2\pi and does not by itself force εf=0\varepsilon_f=0, whereas the holonomy condition does. Curved geometries, which carry nonzero deficit angles and are exactly the physical content of general relativity, are therefore suppressed. This was flagged by Bonzom and given a partition-function-level (wavefront-set) derivation by Hellmann and Kaminski, and confirmed numerically in EPRL and BFBF cases by Donà and collaborators.

Remark 16 (Status of the flatness problem, [H]). The precise scope of [eq:flatness] is contested: whether the flatness holds for all refinements or is an artifact of particular limits and boundary conditions, and whether it survives at finite spin, remain under active investigation. What is not contested is that it is a genuine obstruction to the simplest hoped-for scenario in which refining a spin foam and taking large spins directly yields curved classical gravity. The status is [H]: it is an asymptotic/numerical property of the models, sharply posed, and unresolved.

7.3 Effective spin foams

A recent reframing (Asante, Dittrich, Haggard) confronts the flatness problem by building an effective spin foam directly from the discrete gravity (Area–Regge) action, taking the quantum area spectrum as fundamental and imposing the gluing (shape-matching) constraints as strongly as the uncertainty relations of the area variables allow. In these models curved Regge geometries are not suppressed; the flat dominance of [eq:flatness] is traced to over-imposing constraints that the area uncertainties forbid imposing sharply. This is real progress on the flatness problem, and it clarifies which choices in 11 drive it. It is not a resolution of the continuum-limit problem: effective spin foams are still defined on a fixed complex, and the limit over refinements remains open. We tag the effective-spin-foam resolution of flatness [H] and the continuum limit itself open.

7.4 Divergences and the sum over complexes

Two further honest points. First, at fixed complex the sum over bulk spins can diverge: internal “bubbles” (closed surfaces in the complex) produce factors resembling j(Vj)p\sum_j (\mathop{\mathrm{dim}}V_j)^{p} that diverge for the natural face amplitudes, the analog of loop divergences in lattice gauge theory. The quantum-group truncation of 7 regulates the topological case; for gravity the divergence structure depends on the face amplitude and is part of what a renormalization treatment must control. Second, summing over two-complexes (rather than refining a fixed one) is the group-field-theory strategy; its convergence and its independence of the chosen family of complexes are not established. We tag the sum over complexes [H/P]: heuristic where a concrete group field theory is specified, speculative as a claimed definition of the full nonperturbative theory.

8 Relation to the other regimes and the S/H/P calculus

The library treats twelve quantum-gravity regimes as a modular site, composed by the weakest-link calculus of 1 rather than unified. We locate spin foams among them.

8.1 Spin networks as boundary data (QG-III)

The most direct relation is internal to module QG-III. The boundary Hilbert spaces HΓ\mathcal{H}_{\Gamma} on which spin-foam amplitudes act are the spin-network states of the companion regime: graphs with SU(2)\mathrm{SU}(2) edge labels and intertwiner vertex labels. The kinematics is settled (the kinematical Hilbert space and area/volume spectra are [S]; the representation is unique by the Lewandowski–Okołów–Sahlmann–Thiemann theorem). The dynamics — the spin-foam amplitude between these states — is the bare-[H] link. The composite “spin-network kinematics \ast spin-foam dynamics” is therefore [S/H]    [H]  =  [H],\textsf{[S/H]}\;\ast\;\textsf{[H]}\;=\;\textsf{[H]}, the precise statement that a complete covariant quantum-gravity dynamics inherits the heuristic status of its dynamical half. This is not pessimism; it is bookkeeping. The kinematics really is settled, and the calculus records exactly that the settled part does not settle the whole.

8.2 Classical geometry as the target (QG-I)

Module QG-I — Lorentzian manifolds and the moduli of Einstein solutions — is the target that the missing continuum limit must reproduce. The one rigorous rung of the bridge is the semiclassical asymptotics of 12: a single vertex, on Regge boundary data, reproduces e±iSReggee^{\pm iS_{\mathrm{Regge}}}, the discrete Einstein–Hilbert action. Extending this from one vertex on fixed data to a refined complex with summed bulk, and thence to smooth solutions, is exactly the continuum-limit problem, obstructed by flatness. The composite (classical geometry, [S/H])    (spin foam, [H])=[H],\text{(classical geometry, }\textsf{[S/H]})\;\ast\;\text{(spin foam, }\textsf{[H]})=\textsf{[H]}, and it is optimistic even at [H]: no theorem currently establishes that QG-I is the continuum limit of the state sum.

8.3 Asymptotic safety and renormalization (QG-VII)

The coarse-graining program for spin foams (15) is a renormalization group flow on amplitudes, and it is natural to ask whether a fixed point of that flow could coincide with the non-Gaussian fixed point of asymptotic safety (module QG-VII). Both are bare-[H] entries, so any composite is bounded by [H]; but here the weakest-link bound is an upper bound only. No rigorous bridge exists between a combinatorial refinement fixed point and a continuum theory-space fixed point, so the [H] of the composite should be read as “neither endpoint exceeds heuristic, hence no chain built from them can,” not as “the bridge is established.” The recent tensor-network and Monte Carlo methods make the refinement fixed-point question computationally accessible, which is where a genuine comparison might eventually be made.

8.4 Gauge, BFBF, and cobordism (QG-II, TQFT)

The route from BFBF theory to gravity (4) ties this regime to the Palatini–Cartan and Plebanski formulations of classical gravity (module QG-I) and, through the constraint structure, to the gauge/BV–BRST regime (module QG-II): simplicity is a second-class constraint whose proper treatment is a homological question. The gluing law (6) places spin foams in the cobordism/TQFT column: a spin foam is a morphism between spin networks. The composite with the ([S/H]) cobordism entry stays [S/H] for the topological BFBF models and degrades to [H] for the gravitational ones, whose triangulation dependence breaks the strict TQFT axioms.

8.5 Explicit non-identifications

To keep the modularity honest we record what spin foams are not. A spin foam is not a spacetime: it is a combinatorial object whose relation to a smooth four-manifold is, beyond the topological/asymptotic level, exactly the open problem. The state sum is not a proven path integral for quantum gravity: it is a proposal whose vertex is chosen, not derived. And the sum over refinements is not known to converge: cylindrical consistency fails generically, and the continuum theory is a program, not a theorem. The weakest-link calculus exists to prevent these non-identifications from being quietly elided.

9 A Haskell model of the state sum

To make the reconstruction concrete and checkable we accompany the paper with a small Haskell development, split into base-only modules (association lists, no external containers) so that it compiles with a bare ghc. It models the two layers of the paper that are exactly computable: the SU(2)\mathrm{SU}(2) recoupling bookkeeping underlying the vertex amplitude, and the finite-group BFBF (Dijkgraaf–Witten) state sum that exhibits the refinement invariant exactly.

The recoupling layer (Core) implements SU(2)\mathrm{SU}(2) spins as half-integers, irreducible dimensions 2j+12j+1, the triangle admissibility condition, Clebsch–Gordan fusion, and the dimension of an intertwiner space Inv(iVji)\mathrm{Inv}(\bigotimes_i V_{j_i}) as the multiplicity of the trivial representation, by two independent routes (iterated fusion and, for four legs, the recoupling count). A checkable property asserts the two routes agree, and that a four-simplex boundary labeled by admissible spins has the expected nonzero intertwiner support — the data on which the {15j}\{15j\} vertex [eq:15j] is defined.

The state-sum layer (StateSum, Complexes) implements the amplitude [eq:dw] for G=ZnG=\mathbb Z_n by brute-force enumeration of edge colorings with the face-holonomy constraints, on three explicit triangulations: the boundary of a tetrahedron and an octahedron (two triangulations of S2S^2) and a minimal two-triangle torus. A checkable property asserts 13: that both sphere triangulations give Z=G1\mathcal{Z}=|G|^{-1} (refinement invariance) and that the torus gives Z=G\mathcal{Z}=|G|, matching G1χ|G|^{1-\chi}. This is the exact, solvable version of the paper’s core invariant — the refinement behavior of a state-sum amplitude — in the one case where it is fully under control.

The Main module runs the demonstrations and the property suite and exits with a nonzero status if any property fails. The core module is reproduced in 13; the full sources (Core, StateSum, Complexes, Properties, Main) accompany the submission as ancillary files.

10 Limitations and open problems

We collect the honest status.

10.0.0.1 Continuum/refinement limit (open, [H]).

The central open problem: whether the family {Z(F)}\{\mathcal{Z}(\mathcal{F})\} over refinements is cylindrically consistent, or admits a renormalization-group fixed point, defining a continuum amplitude that reproduces general relativity. Coarse-graining/tensor-network and group-field-theory programs exist but have not closed this. No theorem asserts that module QG-I is the continuum limit of the state sum.

10.0.0.2 The flatness problem (open, [H]).

The naive summed large-spin limit is dominated by flat configurations [eq:flatness]. Effective spin foams reframe the amplitude so that curved configurations survive, which is progress on this specific obstruction, but not a resolution of the continuum limit.

10.0.0.3 Choice of vertex and measure ([H]).

The EPRL/FK vertex is fixed by a chosen (weak) imposition of simplicity, a chosen label identification [eq:eprllabels], and a chosen face amplitude. Different choices give different models; there is no derivation from first principles selecting one.

10.0.0.4 Divergences ([H]).

At fixed complex the bulk sum can diverge through internal bubbles; the divergence structure depends on the face amplitude and is part of the renormalization problem.

10.0.0.5 Lorentzian signature and causality ([H]).

The cosine in [eq:asymp] contains both time orientations; extracting a causal, single-exponential propagator, and controlling the noncompact SL(2,C)\mathrm{SL}(2,\mathbb{C}) integrals, are technical and conceptual issues beyond the scope of the settled asymptotics.

10.0.0.6 Sum over topologies ([H/P]).

Whether and how to sum over the topology of the underlying manifold, and whether the group-field-theory generating function defines the nonperturbative theory, is speculative.

11 Discussion

The reconstruction thesis is sharp in this regime, and so is its limit. What is reconstructed, rigorously, is a family of amplitudes: a colored two-complex, through its state sum, defines a well-defined map between spin-network boundary Hilbert spaces (5), the map composes by gluing (6), the topological version is an exact manifold invariant (7), and a single gravitational vertex reproduces the discrete Einstein–Hilbert action in its large-spin asymptotics (12). At this level the slogan “a spacetime history is a pattern of representation labels on a two-complex, and its amplitude is a product of recoupling invariants” is not a metaphor: there is no metric in the data.

What the thesis has not earned, in this regime, is the continuum. A smooth curved spacetime with its Einstein dynamics is supposed to emerge as the two-complex is refined, and it is precisely there that the construction is unfinished: cylindrical consistency fails generically, the naive large-spin limit is flat, and the renormalization and group-field-theory programs that would fix this are open. The value of the status calculus is that it holds these two facts apart. The settled part is genuinely settled, and the calculus forbids it from lending its credibility to the unsettled continuum limit: a chain is as strong as its weakest warrant, and the weakest warrant here is bare [H].

Seen this way, spin foams are the covariant-dynamics counterpart of the spin-network kinematics: the kinematics gives the boundary states, the foams give the amplitudes between them, and together they would be a background-independent path integral for gravity — if the continuum limit existed. The neighboring regimes supply the rest of the picture: classical Lorentzian geometry is the target, asymptotic safety offers a continuum fixed point to compare against, and the BFBF/cobordism structure supplies the categorical grammar in which a history is a morphism.

12 Conclusion

We have presented covariant quantum geometry as a state sum on a colored two-complex and organized the regime around one invariant: the boundary amplitude and its refinement behavior. The proved statements — well-definedness of the truncated amplitude (5), the topological state-sum invariant (7), the gluing/composition law (6), the exact refinement invariance of the solvable finite-group model (13), and the Regge asymptotics of the EPRL four-simplex vertex (12) — are [S]: standard mathematics in well-defined settings. The physical readings — that the EPRL/FK vertex is the correct gravitational dynamics, and above all that the family of amplitudes has a continuum limit reproducing general relativity — are [H], and the continuum limit is openly unresolved, obstructed at the simplest level by the flatness problem.

The regime composes modularly with its neighbors through the weakest-link calculus, which bounds every composite by its bare-[H] spin-foam link, and the accompanying Haskell model turns the exactly-solvable part of the invariant — refinement invariance of a topological state sum — into a machine-checkable property. The honest headline is narrow and firm: a spin foam is a well-defined amplitude between representation-theoretic boundary states, and its single-vertex semiclassical limit is the Regge action; both are theorems. The broad headline — that this is a continuum theory of quantum gravity — remains a hypothesis, and we have been careful to say so.

13 The core Haskell module

For reproducibility we reproduce the recoupling core, which models SU(2)\mathrm{SU}(2) irreducibles, Clebsch–Gordan fusion, and intertwiner-space dimensions by two independent routes. It uses only base (association lists), compiles with a bare ghc, and underlies the checkable properties of 9. The remaining modules (StateSum, Complexes, Properties, Main) accompany the submission as ancillary files.

{-# LANGUAGE ScopedTypeVariables #-}
-- Module Core: SU(2) recoupling bookkeeping for spin-foam vertices.
-- All content here is the [S] (standard mathematical) layer.
module Core
  ( Spin, spin, halfInt, dimIrrep, casimir, triangle
  , fuse, fuseMany, intertwinerDim, intertwinerDim4
  ) where

-- A spin j is a nonnegative half-integer stored as twice its value (an Int),
-- so j = 1/2 is HalfInt 1, j = 1 is HalfInt 2. Base-only, no Data.Ratio needed.
newtype Spin = Spin Int deriving (Eq, Ord)   -- Spin k represents j = k/2

instance Show Spin where
  show (Spin k)
    | even k    = show (k `div` 2)
    | otherwise = show k ++ "/2"

-- Smart constructor from twice-the-spin (rejects negatives).
halfInt :: Int -> Spin
halfInt k | k < 0     = error ("halfInt: negative " ++ show k)
          | otherwise = Spin k

-- Convenience: integer spin.
spin :: Int -> Spin
spin j = halfInt (2 * j)

-- Dimension of the carrier V_j: 2j + 1 = k + 1.
dimIrrep :: Spin -> Int
dimIrrep (Spin k) = k + 1

-- Casimir j(j+1), returned as twice-numerator over 4 to stay in Int-space:
-- j(j+1) = (k/2)(k/2 + 1) = k(k+2)/4.
casimir :: Spin -> Rational'
casimir (Spin k) = R (k * (k + 2)) 4

-- A tiny rational to avoid Data.Ratio (containers/ratio kept base-simple).
data Rational' = R Int Int
instance Show Rational' where show (R a b) = show a ++ "/" ++ show b

-- Triangle admissibility for j1 (x) j2 -> j3: |j1-j2| <= j3 <= j1+j2 and
-- j1 + j2 + j3 integral (i.e. k1 + k2 + k3 even).
triangle :: Spin -> Spin -> Spin -> Bool
triangle (Spin a) (Spin b) (Spin c) =
  abs (a - b) <= c && c <= a + b && even (a + b + c)

-- Clebsch-Gordan fusion V_a (x) V_b = sum_{k=|a-b|}^{a+b} V_k, step 1 in j
-- (step 2 in the twice-spin integer). Returned as (spin, multiplicity) alist.
fuse :: Spin -> Spin -> [(Spin, Int)]
fuse (Spin a) (Spin b) =
  [ (Spin k, 1) | k <- [abs (a - b), abs (a - b) + 2 .. a + b] ]

-- Iterated fusion of a list of irreducibles, accumulating multiplicities in
-- a base-only association list.
fuseMany :: [Spin] -> [(Spin, Int)]
fuseMany []       = [(spin 0, 1)]
fuseMany (x : xs) = foldl step [(x, 1)] xs
  where
    step acc j = collect
      [ (k, m * mult) | (a, m) <- acc, (k, mult) <- fuse a j ]

-- Sum multiplicities of equal keys in an association list.
collect :: [(Spin, Int)] -> [(Spin, Int)]
collect = foldr ins []
  where ins (k, m) [] = [(k, m)]
        ins (k, m) ((k', m') : rest)
          | k == k'   = (k', m + m') : rest
          | otherwise = (k', m') : ins (k, m) rest

-- Intertwiner-space dimension Inv(V_{j1} (x) ... (x) V_{jn}) = multiplicity
-- of the trivial rep (spin 0) in the full tensor product.
intertwinerDim :: [Spin] -> Int
intertwinerDim js =
  maybe 0 id (lookup (spin 0) (fuseMany js))

-- Independent four-valent count via the recoupling channel:
-- #{ k : max(|j1-j2|,|j3-j4|) <= k <= min(j1+j2,j3+j4), parities match }.
intertwinerDim4 :: Spin -> Spin -> Spin -> Spin -> Int
intertwinerDim4 (Spin a) (Spin b) (Spin c) (Spin d) =
  length
    [ () | k <- [lo, lo + 2 .. hi]
         , even (a + b + k), even (c + d + k) ]
  where lo = max (abs (a - b)) (abs (c - d))
        hi = min (a + b) (c + d)

99

G. Ponzano and T. Regge, “Semiclassical limit of Racah coefficients,” in Spectroscopic and Group Theoretical Methods in Physics (F. Bloch et al., eds.), North-Holland, Amsterdam, 1968, pp. 1–58.

H. Ooguri, “Topological lattice models in four dimensions,” Mod. Phys. Lett. A 7 (1992) 2799–2810; arXiv:hep-th/9205090.

V. G. Turaev and O. Y. Viro, “State sum invariants of 3-manifolds and quantum 6j6j-symbols,” Topology 31 (1992) 865–902.

L. Crane and D. N. Yetter, “A categorical construction of 4D topological quantum field theories,” arXiv:hep-th/9301062; in Quantum Topology, World Scientific, 1993.

M. P. Reisenberger and C. Rovelli, “‘Sum over surfaces’ form of loop quantum gravity,” Phys. Rev. D 56 (1997) 3490–3508; arXiv:gr-qc/9612035.

J. C. Baez, “Spin foam models,” Class. Quantum Grav. 15 (1998) 1827–1858; arXiv:gr-qc/9709052.

J. C. Baez, “An introduction to spin foam models of BFBF theory and quantum gravity,” Lect. Notes Phys. 543 (2000) 25–93; arXiv:gr-qc/9905087.

J. W. Barrett and L. Crane, “Relativistic spin networks and quantum gravity,” J. Math. Phys. 39 (1998) 3296–3302; arXiv:gr-qc/9709028.

J. W. Barrett and R. M. Williams, “The asymptotics of an amplitude for the 4-simplex,” Adv. Theor. Math. Phys. 3 (1999) 209–215; arXiv:gr-qc/9809032.

J. Engle, R. Pereira, and C. Rovelli, “The loop-quantum-gravity vertex amplitude,” Phys. Rev. Lett. 99 (2007) 161301; arXiv:0705.2388.

J. Engle, E. Livine, R. Pereira, and C. Rovelli, “LQG vertex with finite Immirzi parameter,” Nucl. Phys. B 799 (2008) 136–149; arXiv:0711.0146.

L. Freidel and K. Krasnov, “A new spin foam model for 4d gravity,” Class. Quantum Grav. 25 (2008) 125018; arXiv:0708.1595.

J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, H. Gomes, and F. Hellmann, “Asymptotic analysis of the EPRL four-simplex amplitude,” J. Math. Phys. 50 (2009) 112504; arXiv:0902.1170.

J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, F. Hellmann, and R. Pereira, “Lorentzian spin foam amplitudes: graphical calculus and asymptotics,” Class. Quantum Grav. 27 (2010) 165009; arXiv:0907.2440.

V. Bonzom, “Spin foam models for quantum gravity from lattice path integrals,” Phys. Rev. D 80 (2009) 064028; arXiv:0905.1501.

F. Hellmann and W. Kaminski, “Holonomy spin foam models: asymptotic geometry of the partition function,” JHEP 10 (2013) 165; arXiv:1307.1679.

B. Dittrich, “The continuum limit of loop quantum gravity — a framework for solving the theory,” arXiv:1409.1450; in Loop Quantum Gravity: The First 30 Years (A. Ashtekar and J. Pullin, eds.), World Scientific, 2017.

A. Perez, “The spin-foam approach to quantum gravity,” Living Rev. Relativity 16 (2013) 3; arXiv:1205.2019.

C. Rovelli and F. Vidotto, Covariant Loop Quantum Gravity: An Elementary Introduction to Quantum Gravity and Spinfoam Theory, Cambridge University Press, Cambridge, 2014.

P. Donà, F. Gozzini, and G. Sarno, “Numerical methods for EPRL spin foam transition amplitudes and Lorentzian recoupling theory,” Gen. Rel. Grav. 52 (2020) 63; arXiv:1807.03066.

P. Donà, F. Gozzini, and G. Sarno, “Numerical analysis of spin foam dynamics and the flatness problem,” Phys. Rev. D 102 (2020) 106003; arXiv:2004.12911.

F. Gozzini, “A high-performance code for EPRL spin foam amplitudes,” Class. Quantum Grav. 38 (2021) 225010; arXiv:2107.13952.

S. Steinhaus, “Coarse graining spin foam quantum gravity — a review,” Front. Phys. 8 (2020) 295; arXiv:2007.01315.

S. K. Asante, B. Dittrich, and H. M. Haggard, “Effective spin foam models for four-dimensional quantum gravity,” Phys. Rev. Lett. 125 (2020) 231301; arXiv:2004.07013.

S. K. Asante, B. Dittrich, and H. M. Haggard, “Discrete gravity dynamics from effective spin foams,” Class. Quantum Grav. 38 (2021) 145023; arXiv:2011.14468.

S. Steinhaus, “A Monte Carlo algorithm for spin foam intertwiners,” Phys. Rev. D 110 (2024) 026022; arXiv:2403.04836.

C. Rovelli and L. Smolin, “Discreteness of area and volume in quantum gravity,” Nucl. Phys. B 442 (1995) 593–619; arXiv:gr-qc/9411005.

J. C. Baez, “Spin network states in gauge theory,” Adv. Math. 117 (1996) 253–272; arXiv:gr-qc/9411007.

J. Lewandowski, A. Okołów, H. Sahlmann, and T. Thiemann, “Uniqueness of diffeomorphism invariant states on holonomy-flux algebras,” Commun. Math. Phys. 267 (2006) 703–733; arXiv:gr-qc/0504147.

A. Ashtekar and J. Lewandowski, “Background independent quantum gravity: a status report,” Class. Quantum Grav. 21 (2004) R53–R152; arXiv:gr-qc/0404018.