Journal of Gravity

Volume 2013 (2013), Article ID 549824, 28 pages

http://dx.doi.org/10.1155/2013/549824

## Spin Foam Models with Finite Groups

^{1}DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, UK^{2}MPI for Gravitational Physics, Albert Einstein Institute, Am Mühlenberg 1, 14476 Potsdam, Germany

Received 30 March 2013; Accepted 3 June 2013

Academic Editor: Kazuharu Bamba

Copyright © 2013 Benjamin Bahr et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Spin foam models, loop quantum gravity, and group field theory are discussed as quantum gravity candidate theories and usually involve a continuous Lie group. We advocate here to consider quantum gravity-inspired models with finite groups, firstly as a test bed for the full theory and secondly as a class of new lattice theories possibly featuring an analogue diffeomorphism symmetry. To make these notes accessible to readers outside the quantum gravity community, we provide an introduction to some essential concepts in the loop quantum gravity, spin foam, and group field theory approach and point out the many connections to the lattice field theory and the condensed-matter systems.

#### 1. Introduction

Spin foam models [1–8] arose as one possible quantization method for gravity. These models can be seen in the tradition of lattice approaches to quantum gravity [9], in which gravity is formulated as a statistical system on either a fixed or fluctuating lattice. In a rather independent fashion, lattice methods are also ubiquitous in condensed-matter and (Yang-Mills) gauge theories. It is interesting to note that the lattice structures appearing in the strong coupling expansion of QCD or in the high-temperature expansion for the Ising gauge [10] systems are of spin foam type.

In any lattice approach to gravity, one has to discuss the significance of the (choice of) lattice for physical predictions. Although this point arises also for other lattice field theories, it carries extreme importance for general relativity. When it comes to matter fields on a lattice, we can interpret the lattice itself as providing the background geometry and therefore background space time. This could not contrast more with the situation in general relativity, where geometry, and therefore space time itself, is a dynamical variable and has to be an outcome rather than an ingredient of the theory. This is one aspect of background independence for quantum gravity [11, 12].

One possibility, to alleviate the dependence of physical results on the choice of lattice, is to turn the lattice itself into a dynamical variable and to sum over (a certain class of) lattices. This is one motivation for the dynamical triangulation approach [13–15], causal dynamical triangulations [16, 17], and tensor model/tensor group field theory [18–26]. See also [27, 28] for another possibility to introduce a dynamical lattice. Once one decides to sum over lattice structure, one must provide a prescription to do so. The class of triangulations to be summed over can be restricted, for instance, in order to implement causality [16, 17, 29] or to symmetry—reduce models [30–32].^{1}

Another possibility is to ask for models which are *per se* lattice or discretization independent. Such models might arise as fixed points of a Wilsonian renormalization flow and represent the so-called perfect discretizations [33–38]. For gravity, this concept is intertwined with the appearance of discrete representation of diffeomorphism symmetry in the lattice models [36–41]. Such a discrete notion of diffeomorphism symmetry can arise in the form of vertex translations [42–46]. Discrete geometries can be represented by a triangulation carrying geometric data, for instance, in the Regge calculus [47, 48], the lengths of the edges in the triangulation. Vertex translations are then symmetries that act (locally) on the vertices of the triangulations by changing the geometric data of the adjacent building blocks. If the symmetry is fully implemented into the model, then this action should not change the weights in the partition function. Such vertex translation symmetries are realized in 3D gravity, where general relativity is a topological theory [49]. This means that the theory has no local degrees of freedom, only a topology-dependent finite number of global ones. The first-order formalism of 3D gravity coincides with a 3D version of the *BF*-theory [50]. The *BF*-theory is a gauge theory and can be formulated with a gauge group given by a Lie group. The 3D gravity is obtained by choosing . However, one can also choose finite groups [51–53]. In this case, one obtains models used in quantum computing [54] as well as models describing topological phases of condensed-matter systems, so called string-net models [55]. Moreover the *BF*-systems can be seen as Yang-Mills-like systems for a special choice of the coupling parameter (corresponding to zero temperature). Hence, in 3D, we have Yang-Mills-like theories with the usual (lattice) gauge symmetry as well as the topological *BF*-theories with the additional translation symmetries.

The *BF*-theory (in any dimension) is a topological theory because it has a large gauge group of symmetries, which reduces the physical degrees of freedom to a finite number. It is characterized not only by the usual (lattice) gauge symmetry determined by the gauge group but also by the so-called translational symmetries, based on the -cells (i.e., cubes) of the lattice and parametrized by the Lie algebra elements (if we work with the Lie groups). For 3D gravity, these symmetries can be interpreted as being associated to the vertices of the dual lattice (for instance, given by a triangulation) and are hence termed vertex translations. However, for 4D, the symmetries are still associated to the -cells of the lattice and hence to the edges of the dual lattice or triangulation. As in general there are more edges than vertices in a triangulation, we obtain much more symmetries than the vertex translations one might look for in 4D gravity. To obtain vertex translations of the dual lattice in 4D, one would rather need a symmetry based on the -cells of the lattice.

In 4D, gravity is not a topological theory. There is however a formulation due to Plebański [56] that starts with the *BF*-theory and imposes the so-called simplicity constraints on the variables. These break the symmetries of the *BF*-theory down to a subgroup that can be interpreted as diffeomorphism symmetry (in the continuum). This Plebanski formulation is also the one used in many spin foam models. Here, one of the main problems is to implement the simplicity constraints [57–64] into the *BF*-partition function. The status of the symmetries is quite unclear in these models; however, there are indications [40] that in the discrete the reduction of the *BF*-translation symmetries does not leave a sufficiently large group to be interpretable as diffeomorphisms. Nevertheless, if these models are related to gravity, then there is a possibility that diffeomorphism symmetry as the fundamental symmetry of general relativity arises as a symmetry for large scales. In the abovementioned method of regaining symmetries by coarse graining, one can then expect that this diffeomorphism symmetry arises as vertex translation symmetry (on the vertices of the dual lattice).

Hence, we want to emphasize that, in 4D, spin foam models are candidates for a new class of lattice models in addition to Yang-Mills-like systems and the *BF*-theory. This new class would be characterized by a symmetry group in between those for the Yang-Mills theory (usual lattice gauge symmetry) and *BF*-symmetry (usual lattice gauge symmetry and translation symmetry associated to the -cells). Here, we would look for a symmetry given by the usual lattice gauge symmetry and translation symmetries associated to the -cells.

As we will show in this work, the 4D spin foam models, which as gravity models are based on , , or , can be easily generalized to finite groups (or more generally tensor categories [65, 66]). Although the immediate interpretation as gravity models is lost, one can nevertheless ask whether translation symmetries are realized, and if not, how these could be implemented. This question is much easier to answer for finite groups than in the full gravity case. One can therefore see these models as a test bed for the full theory. This also applies for renormalization and coarse graining techniques which need to be developed to access the large-scale limit of spin foams. With finite group models, it might be in particular possible to access the many-particle (that is many simplices or building blocks in the triangulation) and small-spin (corresponding to small geometrical size of the building blocks) regime.^{2} This is in contrast to the few-particle and large-spin (semiclassical) regime [67–69] which is accessible so far.

In the emergent gravity approaches [70–72], one attempts to construct models, which do not necessarily start from gravitational or even geometrical variables but nevertheless show features typical of gravity. The models proposed here can also be seen as candidates for an emergent gravity scenario. A related proposal is the truncated Regge models in [73–76], in which the continuous length variables of the Regge calculus are replaced by values in , with for “the Ising quantum gravity”.

Apart from providing a wealth of test models for quantum gravity researchers, another intention of this work is to give an introduction to some of the spin foam concepts and ideas to researchers outside quantum gravity. Indeed, spin foams arise as graphical tools in high-temperature expansions of lattice theories, or more generally in the construction of dual models. We will review these ideas in Section 2 where we will discuss the Ising-like systems and introduce spin nets as a graphical tool for the high-temperature expansion. Next, in Section 3, we will discuss lattice gauge theories with the “Abelian” finite groups. Here, spin foams arise in the high-temperature expansion. The zero-temperature limit of these theories gives topological *BF*-theories which we will discuss in more detail in Section 4. In Section 5, we detail spin foam models with the non-Abelian finite groups. Following a strategy from quantum gravity, we will also discuss the so-called constrained models which arise from the *BF*-models by implementing simplicity constraints, here in the form of edge projectors, into the partition function. This will in general change the topological character of the *BF*-models to nontopological ones.

We will then discuss in Section 6 a canonical description for lattice gauge theories and show that for the *BF*-theories the transfer operators are given by projectors onto the so-called physical Hilbert space. These are known as stabilizer conditions in quantum computing and in the description of string nets. Here, we will discuss the possibility to alter the projectors in order to break some subset of the translation symmetries of the *BF*-theories.

Next, in Section 7, we give an overview of group field theories for finite groups. Group field theories are quantum field theories on these finite groups and generate spin foams as the Feynman diagrams. Finally, in Section 8, we discuss the possibility of nonlocal spin foams. We will end with an outlook in Section 9.

#### 2. The Ising-Like Models

We will review the construction of the models [77, 78] that are dual to the Ising-like models. Henceforth, they will be referred to as *dual models*. As we will see, similar techniques lead to the spin foam representation for (Ising-like) gauge theories. The main ingredient comprises of a Fourier expansion of the couplings, which can be understood as functions of group variables. We will concentrate on the groups , with for the Ising models proper.

##### 2.1. Spin Foam Representation

Consider a generic lattice, of regular or random type,^{3} in which the edges are oriented. (In principle, this last condition is not necessary for the Ising models; that is, ). One defines an Ising model on such a lattice, with spins associated to the vertices and nearest-neighbour interactions, by utilizing the following partition function:
where , is a coupling constant, which is inversely proportional to temperature, and is the number of vertices in the lattice. Note that the group product is the standard product on , that is, addition modulo .^{4} The action describes the coupling between the spin at the source or starting vertex and the spin at the target or terminating vertex of every edge.

More generally, one can assume that (giving the various vector Potts models). The group product is the standard one on , that is, addition modulo . A further generalization of (1) is to allow the edge weights to be (even locally varying) functions of the two group elements ( and ) associated to the edge : Now, these group-valued functions can be expanded with respect to a basis given by the irreducible representations of the group [79]. For the groups , this is just the usual discrete Fourier transform: where are the group characters. Characters for the Abelian groups are multiplicative; that is, and .

Applying the character expansion to the weights in (2), one obtains (using the multiplicative property of the characters)
where if is the source of and if is the target of . Here, we simply resorted the product so that in the last line the factor involves all of the appearances of . We can now sum over , which according to the Fourier inversion theorem gives a delta function involving the of the edges adjacent to :
where is the -periodic delta function.^{5}

We have represented the partition function in a form where the group variables are replaced by variables taking values in the set of representation labels. This set inherits its own group product via this transform and is known as the Pontryagin dual of . For , the Pontryagin dual and the group itself are isomorphic, as both . We will call a representation of the form of (5), where the sum over group elements is replaced by a sum over representation labels, a *spin foam representation*.^{6} Generically speaking, representations meeting at vertices must satisfy a certain condition. For the Abelian models above, this condition stated that the oriented sum of the representation labels, meeting at a given vertex, had to vanish. In other words, the tensor product of the representations meeting at an edge has to be trivial. In other models, one finds other characteristics as the following. (i) For the non-Abelian groups, this triviality condition generalizes to the choice of a projector from the tensor product of all representations meeting at a vertex into the trivial representation. This choice is encoded as additional information attached to the vertices. (ii) For the gauge theory models we will examine later, the representations are attached to faces, while the triviality condition involves faces sharing a given edge.

Note that the number of representation labels differs from the number of group variables , since the number of edges is generally larger than the number of vertices. In addition, however, these representation labels are subject to constraints. These are referred to as the Gauss constraints, since they demand that the analogue of the electric field (the representation labels) be divergence free (the sum of the representation labels meeting at a vertex has to be zero modulo ).

##### 2.2. Dual Models

To finalize the construction of the dual models [78], one solves the Gauss constraints and replaces the representation labels associated to the edges of the lattice by variables defined on the dual lattice.(2*d*)A particularly clear example arises in two dimensions. On a differentiable manifold, a divergence-free vector field can be constructed from a scalar field via , . On the lattice, a similar construction leads to a dual model with variables associated to the vertices of the dual lattice.^{7} Indeed, for , one finds again the Ising model, just that the high- (low-) temperature regime of the original model is mapped to the low- (high-) temperature regime of the dual model.(3*d*)In three dimensions, a divergence-free vector field can be constructed from another vector field by taking its curl: . However, adding a gradient of a scalar field to does not change : leaves . This translates into a local gauge symmetry for the dual lattice theory. For the dual lattice model, variables reside on the edges of the dual lattice, and the weights define couplings among the edges bounding the faces (that is, the two-dimensional cells) of the dual lattice. Moreover, the aforementioned local gauge symmetry is realized at the vertices of the dual lattice (since the continuum gauge symmetry is parametrized by a scalar field). As a result, the dual model is a gauge theory, a class of theories that we will discuss in more detail in Section 3.(4*d*)In four dimensions, a similar argument to that just made in three dimensions leads to a class of dual models that are “higher-gauge theories” [80, 81]: the variables are associated to the faces of the dual lattice, and the weights define couplings among the faces bounding three-dimensional cells.(1*d*)Finally, let us examine one-dimensional models. Every vertex is adjacent to two edges. Hence, the Gauss constraints, expressed in (5), force the representation labels to be equal (we assume periodic boundary conditions). We can easily find that
For the second equality above, we are restricted to couplings that are homogenous on the lattice; that is, they do not depend on the position or orientation of the edge.

##### 2.3. High-Temperature Expansion and Spin Nets

This spin foam representation is well adapted to describe the high-temperature expansion [82]. Indeed, if one considers the Ising model with couplings as in (1), one obtains the following coefficients in the character expansion: One can now expand the partition function (5) in powers of the ratio , that is, into terms distinguished by the number of times the representation labelled appears.

The “ground” state is the configuration with for all edges. Edges with are said to be “excited”. The Gauss constraints in (5) enforce that the number of “excited” edges incident at each vertex must be even. On a cubical lattice, therefore, the lowest-order contribution arises at from those configurations with excited edges on the boundary of a single plaquette. The full amplitude at this order involves a combinatorial factor, namely, the number of ways one can embed the square (with four edges) into the lattice.

This reasoning extends to all Ising-like models and to higher-order terms in the perturbation series, which can be represented by closed graphs, known as polymers, embedded into the lattice. The expansion in terms of these graphs can then be organized in different ways [82]. Consider one particular contribution to the expansion, that is, a configuration where a some subset of edges carries nontrivial representation labels. In general, one may decompose such a subset into connected components, where one defines two sets of edges as disconnected if these sets do not share a vertex.

We will call such a connected component a *restricted spin net,* where restricted refers to the condition that all edges must carry a nontrivial representation. The Kronecker-s in (5) forbid the spin net to have open ends; that is, all vertices are bivalent or higher. (For the Ising model, with its lone nontrivial representation, only even valency is allowed.) Furthermore, edges meeting at a bivalent vertex, whether there is a “change of direction” or not, must carry the same representation label. (On this point, we have assumed that the couplings depend neither on the positions nor on the directions of the edges; otherwise, one must equip the spin nets with more embedding data.) The representation labels may only change at vertices that are trivalent or higher, so called branching points. This allows one to define *geometric (restricted) spin nets*,^{8} which are specified by two quantities: (i) the *connectivity* of the spin net vertices, that is, the valency of the branching points within the spin net and (ii) the *length* of the edges of the spin net, that is, the number of lattice edges that make up a spin net edge. The contribution of such a spin net to the free energy^{9} can be split into two parts, one is specific to the spin net itself and, from (4), one records that it is given by
The other part is a combinatorial factor determined by the number of embeddings of the spin net into the ambient lattice. This could be summarized as a measure (or entropic) factor, which carries information about the lattice in question, including its dimension.

##### 2.4. Some Considerations for Quantum Gravity

To reiterate, one notes from (8) that the spin net-specific contribution factorizes with respect to the edges of the spin net; for an edge with representation , it is given by where is the length of the spin net edge. Meanwhile, the combinatorial factors are determined by the ambient lattice and its structure. To describe quantum gravity, one is interested rather more in “background-independent” models, where such factors should not play a role. In this context, one can define an alternative partition function over geometric spin nets, whose edges carry both a representation label and an edge-length label. This length variable encodes the geometry of space time as a dynamical variable on the geometric spin net.

From that point of view, “background independence” motivates the search for models, within which the spin net weights are independent of these edge lengths. These can be termed *abstract spin nets* [83]. (Note that, with this definition, the abstract spin nets generally still remember some features of the ambient lattice; for instance, the valency of the spin net cannot be higher than the valency of the lattice.) For the Ising model, such abstract spin nets arise in the limit , that is, for zero temperature. In contrast to the high-temperature expansion, which is a sum ordered according to the length of the excited edges, abstract spin networks arise in a regime where this length does not play a role.

A similar structure arises in string-net models [55], which were designed to describe condensed phases of (scale-free) branching strings. Indeed, these string-net models are closely related to (the canonical formulation of) topological field theories, such as the *BF*-theories, which we will describe in Section 6. In this canonical formulation, spin net(work)s appear as one choice of basis for the Hilbert space of the *BF*-theories. String-net models assign amplitudes (corresponding to the so-called physical wave functions in the *BF*-theory) to spin networks that satisfy certain (gauge) symmetry properties. This implies that these spin networks can be freely deformed on the lattice without changing their amplitude.

In the following section, we will consider gauge systems, for which we will devise once again a spin foam representation. For these gauged models, rather than labelling the terms in the expansion using spin nets, we will utilize a different, yet similar structure known as a *spin foam*. These have a similar structure to spin nets, except that the various roles are played by simplices one dimension higher. To clarify, the roles of the vertices and edges are assumed by the edges and faces, respectively. With this proviso, other concepts introduced above translate nicely. Spin foam edges and faces are generally comprised of several edges and faces in the ambient lattice. Moreover, one can define both geometric and abstract spin foams [83]. Both are branched surfaces [1–6, 84] whose constituent faces carry a nontrivial representation label.^{10} The amplitudes for geometric spin foams will most often depend on the area of their surfaces, that is, the number of plaquettes in the ambient lattice constituting the spin foam surface in question, as well as the lengths of their branches, that is, the number of edges in the ambient lattice constituting the spin foam branch in question. For instance, this is the case in the Yang-Mills-like theories, where geometric spin foams arise in the strong coupling or high-temperature expansion.

The conditions for *abstract spin foams*, that is, amplitudes independent of the surface area, may be understood as requiring that the amplitudes be invariant under trivial (face and edge) subdivision. This can be used to specify measure factors for the spin foam models [1–5, 83, 85–87]. Abstract spin foams are also generated by the tensor models discussed in Section 7; however, the spin foams thus generated are not “restricted” in the manner we defined earlier; that is, the trivial representation label is allowed. Furthermore, the spin foams need not be embeddable into a given lattice. See also [88] for one possible connection between partition functions with restricted and unrestricted spin foams.

This discussion on geometric and abstract spin nets and spin foams assumes that the couplings do not depend on the position or direction of the edges or plaquettes in the lattice, otherwise, one would need to introduce more decorations for the spin nets and foams. However, one can show [89], for instance, that the Ising gauge model in 4D with varying couplings is universal. In other words, one can encode all other lattice models. Here, it would be interesting to develop the equivalent spin foam picture and to see how a spin foam with additional decorations could encode other spin foam models.

#### 3. Lattice Gauge Theories

The spin foam representation originated as a representation of partition functions for gauge theories [1–6, 11]. Such gauge theories can also be formulated with discrete groups [10, 52, 90], and for , they give the Ising gauge models [10]. A spin foam representation for the Yang-Mills theories with continuous Lie groups has been presented in [91], see also [92, 93], and the results can be easily adapted to discrete groups. For the Abelian (discrete and continuous) groups, these are again closely related to the construction of dual models [78] and the high-temperature (or strong coupling) expansion [84].

Given a lattice, one associates group variables to its (oriented) edges. To be more precise, a lattice in this context is an orientable 2-complex , within which one can uniquely specify the 2-cells, called faces or plaquettes. The aim is to construct models with gauge symmetries at the vertices, where the gauge action is given by and the are the gauge symmetry parameters. They are associated to the vertices of the lattice; remember that is the source vertex of and is the target vertex. Gauge-invariant quantities are given by the Wilson loops.^{11} Expressing the action as a function of these quantities ensures that it in turn is gauge invariant, as can be easily checked. A Wilson loop is the trace (in some matrix representation) of the oriented product of group variables associated to a loop of edges in the lattice.^{12} The “smallest” Wilson loops are those constructed from edges around a single lattice face. The associated face variables are then , where denotes the oriented product.

Given a unitary (finite-dimensional) matrix representation of a group , a Wilson-styled action is given by Here, is a coupling constant and is a matrix representation of . As we saw earlier, for the groups , the irreducible unitary representations are 1 dimensional and labelled by so that the representations matrices are realized by .

As before, let us generalize to a class of partition functions of the form where the weights are class functions; that is, just like the trace, they are invariant under conjugation . Moreover, we included a measure normalization factor for every summation over the group .

##### 3.1. Spin Foam Representation

The face weights , as class functions, can be expanded in characters. For the Abelian groups , this again just amounts to the discrete Fourier transform. Thus, the expansion takes the initial form while the derivation of the spin foam representation proceeds as in Section 2 as follows: where if the orientation of the edge coincides with the one induced by the “adjacent” face and otherwise. In the last step, we simply replace the Kronecker- weighting each edge by its square and associate one factor to each vertex. The reason is that in the spin foam representation one usually works with vertex amplitudes, which in this case is . Having said that, these amplitudes and the splitting are more complicated for both the non-Abelian groups (see Section 5) and continuous groups.

##### 3.2. Dual Models

The construction of the dual models [78] proceeds by solving for the Gauss constraints associated to the edges in (12), that is, by solving the Kronecker-s. These enforce that the oriented sum of representation labels, associated to the faces incident at a given edge, vanishes.(2*d*) For the lowest-dimensional case, namely, two dimensions, one obtains a “trivial” dual theory. Since every edge lies on just two faces, all of the representation labels coincide: , and one obtains (assuming that the do not depend on the position and orientation of the plaquettes, )
(3*d*) For the three-dimensional case, we have already seen that lattice models without any gauge symmetry are dual to lattice gauge models. For example, the Ising gauge model in is dual to the Ising model.(4*d*) In four dimensions, lattice gauge models are dual to lattice gauge models. In the case that the cohomology of the lattice is nontrivial, topological models might appear as the dual model [94].

##### 3.3. High-Temperature Expansion and Spin Foams

As for the Ising models, there is a high-temperature expansion, which for the Ising gauge model is also an expansion in powers of . The discussion is very similar to the one for the Ising models, with the difference being that the perturbative contributions are now represented by closed surfaces rather than closed polymers.

The “ground state” is the configuration with for all faces. The Gauss constraints demand that the number of excited plaquettes (those with ) incident at every edge must be even. In particular, this means that the excited plaquettes must form closed surfaces. In the case of the Ising gauge model on a hypercubical lattice, the lowest-order contribution arises at from those configurations with excited faces on the boundary of a single cube. As one might expect, there is a combinatorial factor arising from the number of embeddings of this cube into the ambient lattice.

More generally, the perturbative contributions are now described by closed, possibly branched surfaces. (For the free energy, one has only to consider connected components [95].) Here, a proper branching is a sequence of edges where at least three excited plaquettes meet. Excited plaquettes make up the elementary surfaces or spin foam faces. In other words, at the inner edges of these spin foam faces, there are always only two excited plaquettes meeting. Again the Gauss constraints demand that the representation labels of all of the plaquettes in a given spin foam face agree. For homogeneous and isotropic couplings, the spin foam amplitude (ignoring the combinatorial embedding factors) depends on the number of plaquettes inside each spin foam face, but not on how the spin foam is embedded into the lattice. The contribution of a spin foam face with label is . (In general, for models with the non-Abelian groups, the amplitudes might also depend on the length of the branching or spin foam edges.) This definition of geometric and abstract spin foams [83] parallels the one given for spin nets. Note that the condition for abstract spin foams invalidates the conditions of the high-temperature expansion, which is an expansion in the number of excited plaquettes of the underlying lattice.

In Section 4, we will discuss the *BF*-theory. This is a topological model that satisfies the conditions necessary to class it as an abstract spin foam. In other words, the amplitude does not depend on the areas of the spin foam faces. The reason is that for the *BF*-models. For the Ising gauge model, this is the case at zero temperature.^{13}

##### 3.4. Expansion around a Topological Sector

If one starts from the first line in (11) and keeps the summation over both group and representation labels, one obtains
This corresponds to a first-order representation of the Yang-Mills theory where the are the connection variables and the represent the dual (electric field or flux) variables. The model where all is a topological *BF*-model, whose partition function can be solved exactly. Hence, the Yang-Mills theory can be understood as a deformed *BF*-theory [92, 93, 96, 97]. The deformation appears here as the term; in the continuum, it is .^{14}

The expansion of models with propagating degrees of freedom, for instance, those modelling gravity and the Yang-Mills theories, around topological theories has been discussed, for instance, in [83, 98, 99].

In the case of the Ising model, one expands around the *BF*-theory partition function by introducing an expansion parameter − , which is small for low temperatures:
Note that to get to the second equality we have used the fact that . Thus, for a given configuration , the coefficient of records the number of excited faces. The coefficient of gives the number of ordered pairs of excited faces. The next coefficient records the number of ordered triples of faces and so on. The last coefficient of is one for the configuration and zero for all other configurations. The terms in this expansion can be seen as expectation values of observables in a *BF*-model, where the observables are the numbers of ordered -tuples of excited plaquettes, for each .

#### 4. Topological Models

In the previous section, we saw that theories of the Yang-Mills type can be understood as a deformation of the *BF*-theories, which we will discuss here in more detail. From the conventional standpoint, represents the curvature of a connection, usually taking values in the Lie group, while is a Lie algebra-valued -form. Hence, in constructing these models for discrete groups, the following question arises: with what should one replace the variables? We have seen that the -field corresponds to the representation labels. Examples for topological models with discrete groups are discussed in several texts [51, 53, 100, 101]. See also [102, 103] for the relation between the Chern-Simons-like theories and the *BF*-like theories in with discrete group . In this section, we will restrict to the Abelian groups, in particular . The models for the non-Abelian groups will be presented in Section 5.

The *BF*-theory is a topological field theory in which the equations of motion demand that the “local” curvature vanishes. One often starts with the partition function encoding this requirement. It is defined in a manner very similar to that of the previous sections (see, for instance, [6, 104]) as follows:
where again means the oriented product^{15} of edges around a face. Note that this is just a special case of the gauge theory partition functions (10), obtained by setting .

The partition function can be easily evaluated:
For the groups , an expansion of the type detailed in (14) yields
One may interpret the term in the exponential as a *BF*-theory action, where the representation labels represent the -field on the lattice, while the products represent curvature.

The spin foam representation can be derived in the same way as for the gauge theories in Section 3.1:
This recasts the partition function as
that is, as the number of assignments satisfying the Gauss constraint for every edge. In the following section, we will discuss the local symmetries of the *BF*-theory partition function. We will find that, on the set of all configurations , there acts a translation symmetry. This is realized at the 3-cells of the lattice. Hence, reflects the orbit volume of this translational gauge symmetry. It is this symmetry that leads to divergences in the *BF*-theory partition functions for the “continuous” Lie groups [105, 106].^{16}

##### 4.1. Symmetries of the Partition Function

We will now discuss the gauge symmetries of the partition function (16). As its form coincides with that of the gauge theories (10), the partition function is invariant under the usual gauge transformations , where are gauge parameters associated to the vertices. This is manifest in the representation given in (16).

There is a further symmetry, the so-called translation symmetry, which is easiest to see in the spin foam representation. This symmetry is based on the -cell of the lattice, that is, the cubes for a hypercubical lattice.^{17} Consider a field of gauge parameters associated to the cubes and define the gauge transformations
where if the orientation of agrees with that induced by and otherwise. If is a configuration that satisfies all of the Gauss constraints, then so does and vice versa. The reason stems from the following fact: if an edge is in the boundary of a -cell , then there are exactly two faces, say and , that are both in the boundary of and adjacent to . The orientation factors are such that the contributions of to the two faces and are of opposite signs in the Gauss constraint associated to . Hence, the stated result follows and the contributions of these two configurations, and , to the partition function are equal.

The above symmetry is a realization of the well-known cohomological principle that “the boundary of a boundary is zero” () and underlies the Bianchi identity for the curvature. With the Bianchi identity, one can also explain the translation symmetry; see [39, 49].

In , the *BF*-theory with as its gauge group has a gravitational interpretation. This translation symmetry corresponds to the diffeomorphism symmetry of the theory. The gauge field corresponds in the dual lattice to a field associated to the dual vertices, and one can interpret the gauge transformation as a translation^{18} of these dual vertices (which in the gravity models are the vertices in a triangulation of space time).

In , the -cells are dual to edges in the dual lattice, so this symmetry translates the dual edges rather than the dual vertices. For gravity-like models, on the other hand, one would like to break down this translation symmetry, where it is based on the dual edges, to symmetries based on the dual vertices. (Correspondingly, in , diffeomorphism symmetry of the action leads to the Noether charges encoding the contracted Bianchi identities and not the Bianchi identity itself.) In the case of gravity, one strategy is to impose the so-called simplicity constraints [56–61, 64] within the partition function. For models with finite groups, the question arises as to whether there is a class of gauge models with “translation-like” symmetries based on the dual vertices. Such models would be nontopological and feature propagating degrees of freedom, as a symmetry group based on dual vertices is much smaller than that for BF, where it is based on dual edges. See also the discussion in Sections 5 and 6.

#### 5. Gauge Theories for the Non-Abelian Groups

In the following, we will consider how to generalize the construction performed so far for the finite, but non-Abelian, groups . Details about representations can be found in [79].

Again, we denote by the irreducible, unitary representations of on , where . Note that for the non-Abelian groups can be larger than . Every function can be decomposed into matrix elements of representations; that is, where the sum ranges over all “equivalence classes” of irreducible representations of . The are given by where denotes the number of elements in the group . Introducing an inner product between functions by then the matrix elements of are orthogonal; that is, Furthermore, we denote by the character of ; then the -function on the group is given by which satisfies

##### 5.1. The -Gauge Theory on a Two-Complex

Consider an oriented two-complex^{19}. Denote the set of edges and the set of faces ; then by a *connection* we mean an assignment of group elements to edges . For every face , we denote the *curvature* of this connection by
which is the ordered product of group elements of edges in the boundary of , where we take depending on the relative orientation of and . In the non-Abelian case, the ordering of the group elements around a face/plaquette actually matters, and we choose here the convention as depicted in Figure 1.

As before, we define a gauge-invariant partition function by choosing a collection of weight-functions invariant under conjugation. These encode the action “and the path integral measure” of the system. The partition function is defined as
where is the number of elements in . Since can be expanded into characters via (22), the path integral can be written as
In sum, there is for each edge a representation appearing for every face adjacent to , . The corresponding sum results in (assuming that the orientations of all agree with those of for the moment, in order not to overburden the notation)
It is not hard to see that the operator , called the *Haar-intertwiner*, given by
which maps the edge-Hilbert space
to itself, is actually an orthogonal projector on the gauge-invariant subspace of .^{20} Note that, while in the Abelian case one has to sum over all representations over faces such that at all edges the “oriented” sum adds up to zero (as in (12)), in the non-Abelian generalization one has to sum over all representations such that at each edge the representations of the faces meeting at have to contain, in their tensor product, the trivial representation. Also, one obtains a nontrivial tensor for each edge, which in the case of the Abelian gauge groups is just the -number . Since the representations for the non-Abelian groups can be more than dimensional, in general the tensor has indices which are contracted at each vertex in the two-complex .

Choosing an orthonormal base , for the invariant subspace of , we get
In case the orientations of and do not agree for some , then is essentially replaced by in (32), which leads to the appearance of the dual^{21} representation in the tensor product (34) of the . In this case, labels a basis of intertwiner maps between the tensor product of all representations associated to faces with and the tensor product of all representations associated to the faces with .

In (35), one regards and as attached to, respectively, the endpoint and the beginning of . For each vertex in , the associated and can be contracted in a canonical way. For every face which touches , there are exactly two edges , in the boundary of that meet at . The definitions above are exactly such that, if one chooses a basis in each and the dual basis in each , in and the two indices associated to are in opposite position, so they can be contracted. See Figure 2 for an example.

Therefore, contracting all the appropriate at one vertex leaves one with the vertex-amplitude , which depends on the representations and intertwiners associated to the faces and edges that meet at (and the orientations of these). The vertex amplitude can be computed by evaluating the so-called *neighbouring spin network function*, living on graph which results from a dimensional reduction of the neighbourhood of . Construct a spin network, where there is a vertex for each edge touching , and a line between any two vertices for each face between two edges. Assigning the to the lines of the spin network and the intertwiners or to the vertices, depending on whether the corresponding edge is incoming or outgoing of , results in a spin network, whose evaluation gives .

With this, the state sum can, using (31) and (32), be written in terms of vertex amplitudes via

##### 5.2. Examples

The first example we consider is the *BF*-theory, which corresponds to an integral over only flat connections, that is, the choice . Since the -function can be decomposed into characters with , we get
If is the two-complex dual to a triangulation of a -dimensional manifold, then the vertex amplitude is essentially the analogue of the -symbol for ^{22}.

The second example that one usually considers is the Yang-Mills theory. The Wilson action can be specified, similar to the Abelian case, by choosing a unitary representation , so that where is the coupling constant. The weights are then given by .

##### 5.3. Constrained Models

There is a generalization of the state-sum models (36), coming from the desire to obtain nontopological generalizations of the *BF*-theory. It amounts to choosing the intertwiners not to span all of the invariant subspace, but only a proper subspace . This originates in the attempt to define a state-sum model for general relativity, which can be written as a constrained *BF*-theory. The subspace is viewed as the space of intertwiners satisfying the so-called simplicity constraints; see, for example, [56, 64, 107]. Examples for such models are the Barrett-Crane model [57] and the EPRL and FK models [58–61]. These models can also be written in the form of a path integral over connections [91, 108, 109], but we will not concern ourselves with this here.

In the following, we will introduce a class of models which can be seen as the generalization of the Barrett-Crane models to finite groups. Given a finite group , the model is a state-sum model for the group . Irreducible representations of are then pairs of irreducible representations of . In the edge-Hilbert space , there is a specific element of the space of gauge-invariant elements, called the *Barrett-Crane intertwiner*. It is only nonzero if and only if the representations and are dual to each other. In that case, for an edge with attached faces , consider the Hilbert space
Again, we have assumed that the orientations of all faces and agree, in order not to overburden the notation; otherwise, replace by its dual or, equivalently, exchange . If we consider the projector onto the subspace of elements invariant under the action of , then can be seen as an element of , which is naturally isomorphic to , because are dual to each other. The corresponding element in is invariant under , as one can readily see, and it is our distinguished element onto which is projecting.

Since the projection space for each edge is (at most) dimensional, the vertex amplitude for the BC-model does not depend on intertwiners, but only on the representations associated to the faces. In fact, since the representations are unitary, every representation is dual to itself, so we effectively have only one representation attached to each face. Using diagrammatical calculus [110], it is not hard to show that the vertex amplitude for the BC model is given by a sum over squares of the *BF*-theory model on the same vertex: that is,
where the sum ranges over an orthonormal basis of -intertwiners in for every edge at the vertex .

Let us shortly discuss the case of , the group of permutations in three elements. For , which is generated by , the permutation of the first two elements, and , which is the cyclic permutation, subject to the relation , the representations theory is well known [79]: There are three irreducible representations, called , , and . The first is the trivial representation and the second one maps a permutation to its sign . The third one is two dimensional, and The nontrivial tensor products of the representations decompose as Therefore, the intertwiner space in is, for example, three dimensional, when there are four faces attached to carrying the representation . Hence, for , the BC-model is an example for a constrained version of an -state-sum model, as defined in the last section, because projects onto a proper subspace of the invariant subspace of .

This class of models is an example for an abstract spin foam [83], since its amplitudes only depend on combinatorial data, and the topology of the two-complex . In particular, it leads to a model which is invariant under refinement of the two-complex by trivial subdivisions of edges or faces. Hence, these models provide potential examples for models which are background independent, but not topological (as is the case in the full theory).

Note that there is a wealth of different models, which come from different choices of nontrivial subspaces of , onto which is projecting. In particular, one could consider EPRL-like models where or and consider the subspace of intertwiners for , which are in the image of a boosting map , given by the fusion coefficients (see, e.g., [58–61] for ). Here, () is the group of “even” permutation of four elements. Another, possibly more geometric way would be to consider embeddings and thus construct proper subspaces of -intertwiners.^{23} Such models can be considered as truncations of the EPRL models as () is the “chiral” symmetry group of the tetrahedron and a subgroup of the rotation group. Here, it will be interesting to investigate the relation to the full models in more detail.

The models can serve to test many proposals for the full theory, for instance, how the different choices of edge projectors determine the physical degrees of freedom or particle content of the given theory. This is related to the problem of implementing the simplicity constraints into spin foam models. We are planning to investigate the features of these models, in particular their symmetries [65a, 65b] and behaviours under coarse graining, in future work.

#### 6. The Hilbert Space, the Transfer Operators, and Constraints

We intend to derive the *transfer operators* [111, 112] for the Yang-Mills-like gauge theories, having partition functions of the forms (10) and (30) and incorporating a generic finite group of cardinality . We will permit to be non-Abelian. The first part of the discussion is of a similar nature to that found in [112] for the Yang-Mills theory. However, we will also include the *BF*-theory as a special case. From the transfer operators, one can obtain, via a limiting procedure, the Hamiltonian operators. However, for the *BF*-theory, we will see that this limiting procedure is not necessary. Rather, the transfer operators can be understood as projectors onto the space of the so-called physical states. These can be characterized as lying in the kernel of the “Hamiltonian” constraints. This illustrates the general principle that a path integral with local symmetries acts a projector onto a constrained subspace [113, 114]. Conversely, one can construct a projector onto the constrained subspace as a path integral. See [115], in which this is performed for gravity in its formulation as an *BF*-theory.

The transfer operator is defined on a Hilbert space , so that the partition function can be written as
We assume for simplicity^{24} that the lattice is hypercubical. One lattice direction is designated as a time direction, in which there are periodic boundary conditions. The trace is the trace over states in the following Hilbert space. It is associated to the spatial lattice and is built as a direct product of the Hilbert spaces associated to the spatial edges . More succinctly written as , the “configuration” variable attached to any given edge is a group element, so the associated Hilbert space is the space of complex functions on the group, equipped with an inner product,^{25} that is, .

We will work in the representation in which a basis of eigenstates is given by . This is known as the holonomy or connection representation. For any unitary matrix representation, of the group , these are simultaneous eigenstates of the so-called edge holonomy operators : These basis states are orthonormal: where is the delta function on the group such that . We have also the completeness relation We use this resolution of identity times in (43) to rewrite the partition function as where we introduced to label the “constant time hypersurfaces” in the lattice. On the other hand, we will assume that our partition function is of the form Here, we introduced the weight functions of the holonomies around plaquettes , which in the form of the right-hand side can also be chosen to differ for spatial plaquettes and timelike plaquettes . (This is necessary if one wants to obtain the Hamiltonian in the limit of small lattice constant in timelike direction.) Since we are dealing with a gauge theory, the weights are class functions, that is, invariant under conjugation; furthermore, we will assume that ; that is, the weights are independent of the orientation of the face. A weight function can therefore be expanded into characters, which are linear combinations of matrix elements of representations (in fact, characters are just the trace). Hence, the weight functions will be quantized as edge holonomy operators.

On the right-hand side of (48), we have just split the product into factors labelled by the time parameter . Here, we assume that the plaquettes are the ones between time and . Comparing the two forms of the partition functions (47) and (48), we can conclude that The two outer factors can be easily quantized as multiplication operators, so that we can write with

To tackle the operator , note that the group elements associated to the timelike edges appearing in the plaquette weights can be understood as acting as gauge transformations associated to the vertices of the spatial lattice; see Figure 3.

That is, we define operators by where and are the source and target vertices of the edge , respectively. The operators generate gauge transformations as Now, we can see the plaquette weight in (52) as either a wave function at time or a wave function at time + on which the group elements associated to the timelike edges act as gauge transformations. The sum in (51) over the group elements then induces an averaging over all gauge transformations, that is, a projection onto the space of gauge-invariant states. Hence, we can write where The operator is a projector; that is, . One can easily show that for any ; hence, it projects onto the space of gauge-invariant functions. As the operator is made up from gauge-invariant plaquette couplings, it commutes with , so that we can write Here, we see an example for a general mechanism, namely, that the transfer operator for a partition function with local gauge symmetries acts as a projector onto the space of gauge-invariant states. We can characterize such gauge-invariant states as being annihilated by constraint operators. In the case of the usual lattice gauge symmetries, these constraints are the Gauss constraints, which we will discuss later on.

To discuss the action of , we introduce first the spin network basis, in which will be diagonal.

##### 6.1. Spin Network Basis

So far, we encountered the holonomy operators, that is, the multiplication operators in the configuration representation. The conjugated operators are the electric fields or fluxes, which act as “matrix” multiplication operators in the spin network basis. This is a convenient basis to discuss the remaining part of the transfer operator.

To obtain the spin network basis [11, 12], we just need to use that every function on the group can be expanded as a sum over the matrix elements of all irreducible unitary representations . For the Abelian groups, all irreducible representations are one dimensional; hence, and we obtain the discrete Fourier transform (3). In the general case of the non-Abelian groups, we define spin network states by which are orthonormal; that is, .

In the spin network basis, we can easily express the left and right translation operators . These are the conjugated operators to the holonomy operators (44). To simplify notation, we will just consider the Hilbert space and states associated to one edge and omit the edge subindex. We define In the spin network basis, where we sum over repeated indices. That is, The remaining factor of the transfer operator factorizes over the edges, and it is straightforward to check that it acts on one edge as (Here, one has to use that is a class function and invariant under inversion of the argument.) On a spin network state, we obtain where in the second line we used that a class function can be expanded into characters and in we used the third the orthogonality relation The expansion coefficients are given by That is, acts as a diagonal in the spin network basis, and as the eigenvalues only depend on the representation labels , it commutes with left and right translations. For the Yang-Mills theory (with Lie groups) in the limit of continuous time, one obtains for the Laplacian on the group [111, 112].

##### 6.2. Gauge-Invariant Spin Nets

Given a spin network function , which can be regarded as function of holonomies associated to edges, where is a connection, and a gauge transformation , an assignment of group elements to vertices of the network, then the gauge transformed spin network function is given by where and are the source and the target vertices of the edge . A gauge transformation can therefore be expressed as a linear operator in terms of and , as shown in the last section. Decomposing the spin network function into the basis , that is, one can readily see that the condition for a function to be invariant under all gauge transformations translates into a condition for the coefficients ; namely, where the product ranges over vertices , and each tensor , which has the indices for each edge outgoing of and for each edge incoming , is an invariant element of the tensor representation space that is, an intertwiner between the tensor product of incoming and outgoing representations for each vertex. An orthonormal basis on the space of gauge-invariant spin network functions therefore corresponds to a choice of orthonormal intertwiner in each tensor product representation space for each vertex, where the inner product is the tensor product of the inner products on each , .

Note that neither the holonomies nor left and right translations are gauge invariant; therefore, they map gauge-invariant functions to nongauge-invariant ones. However, they can be combined into gauge-invariant combinations, such as the holonomy of a Wilson loop within a net.

##### 6.3. Local Constraint Operators and Physical Inner Product

We have seen that for a general gauge partition function of the form of (48) the transfer operator (57) incorporates a projection onto the space of gauge-invariant states. This is a general feature of partition functions with gauge symmetries [113, 114]. Indeed in quantum gravity, one attempts to construct a projector onto gauge-invariant states by constructing an appropriate partition function.

As has been discussed in Section 4.1, the *BF*-models enjoy a further gauge symmetry, the translation symmetries (21). (A generalization of these symmetries holds also for the non-Abelian groups.) Hence, we can expect that another projector will appear in the transfer operator. Indeed, it will turn out that the transfer operator for the *BF*-models is just a combination of projectors.

This is straightforward to see as for the *BF*-models that the plaquette weights are given by , so that
Here, one might be worried by taking the square roots of the delta functions; however, we are on a finite group. The operator is diagonal in the basis and has only two eigenvalues: on states where all plaquette holonomies vanish and zero otherwise. The square root of is defined by taking the square root of these eigenvalues.

Hence, is proportional to another projector , which projects onto the states for which all the plaquette holonomies are trivial, that is, states with zero (local) curvature.

The final factor in the transfer operator is which, according to the matrix element of in (56) (putting ), is proportional to the identity
where the numerical factor arises due to our normalization convention for the group delta functions, which amounts to . Hence, the transfer operator for the *BF*-theory is given by
and since for the projectors we have , , the partition function simplifies to
The projection operators in the transfer operator ensure that only the so-called physical states are contributing to the partition function . These are states in the image of the projectors (we will call the corresponding subspace ) and can be equivalently described to be annihilated by constraint operators, which we will discuss below.

The objective in canonical approaches to quantum gravity [11, 12] is to characterize the space of physical states according to the constraints given by general relativity, which follow from the local symmetries of the theory. (The quantization of these constraints in a consistent way is highly complicated [116–118].) Indeed, also in general relativity, as well as in any theory which is invariant under reparametrizations of the time parameter, one would expect that the transfer operator (or the path integral) is given by a projector, so that the property would hold.^{26} This would also imply a certain notion of discretization independence, namely, the independence of transition amplitudes on the number of time steps used in the discretization; see also the discussion in [38] on the relation between reparametrization symmetry in the time parameter and discretization independence. For a lattice based on triangulations, one can introduce a local notion of evolution, the so-called tent moves [41, 119–121]. These moves evolve just one vertex and the adjacent cells. Local transfer operators can be associated to these tent moves. These would evolve the spatial hypersurface not globally, but locally, and would allow to choose some arbitrary order of the vertices to evolve. In this way, one can generate different slicings of a triangulation as well as different triangulations. The condition that the transfer operator associated to a tent move is given by a projector would then implement an even stronger notion of triangulation independence.

However, reparametrization invariance, which would lead to such transfer operators, is usually broken by discretizations already for one-dimensional toy models. For “coarse graining and renormalization” methods to regain this symmetry, see [36–38, 122]. In the case of gravity, or more generally nontopological field theories, these methods will however lead to nonlocal couplings for the partition functions [122], for which one would have to modify the transfer operator formalism.

Furthermore, one has to construct a physical inner product on the space of the physical sates. This process is quite trivial in the case considered here, as we are considering finite groups and all of the projectors are proper projectors. Hence, the physical states are proper states in the “so-called” kinematical Hilbert space , and the inner product of this space can be taken over for the physical states. In contrast, already for the *BF*-theory with “compact” Lie groups, the translation symmetries lead to noncompact orbits. This leads to physical states which are not normalizable with respect to the inner product of the kinematical Hilbert space. For the intricacies of this case (with ), see, for instance, [115]. A further complication for gravity is the very complicated non-commutative structure of the constraints [123–125].

Let us summarize the projectors onto . Also in the case of generalized projectors, a physical inner product can be defined [12, 126] between equivalence classes of kinematical states labelled by representatives through Kinematical states which project to the same (physical) state define the same equivalence class. This leads, for instance, to an identification of the two states (for the group ) in Figure 4. This is analogous to the action of the spatial diffeomorphism group in “ and ” loop quantum gravity, see the discussion in [115], and reflects the concept of abstract spin foams, whose amplitude does not depend on the embedding into the lattice, in Section 2. Indeed any two states which can be deformed into each other by applying the so-called stabilizer operators introduced below in (76) and (82) are equivalent to each other. The reason is that the physical Hilbert space is the common eigenspace to the eigenvalue with respect to all of the stabilizer operators. Hence, any part of a state that undergoes a nontrivial change if a stabilizer is applied will be projected out in the physical inner product.

Another problem in quantum gravity is then to find “quantum” observables [127–131] that are well defined on the physical Hilbert space. In the case of gravity, these Dirac observables are very hard to obtain explicitly; even classically there are almost none known. In particular, such Dirac observables have to be nonlocal [132]. If one interprets a quantum gravity model as a statistical system, a related task is to find a well-defined order parameter characterizing the phases. Expectation values of gauge-variant observables may be projected to zero under the physical inner product [133]: As in our case, we work with finite-dimensional Hilbert spaces, and the physical Hilbert space is just a subspace of the kinematical one; there is a construction principle for physical observables. For any operator on the kinematical Hilbert space, one can consider , which preserves the physical Hilbert space and corresponds to a “gauge averaging” of . However, many observables constructed in this way will turn out to be constants. For the *BF*-theory, observables can be constructed by considering a product of the stabilizer operators discussed below along non-contractable loops [54, 134–136]. The resulting observables are nonlocal, and the cardinality of a set of independent operators (equivalently the dimension of the physical Hilbert space) just depends on the topology of the lattice and not on its size.

We will now discuss the stabilizer operators which characterize the physical states. In topological quantum computing, these are known as stabilizer conditions [54].

We have considered the operators in (53) that implement a gauge transformation according to the assignments of group elements to the vertices of the “spatial sub-” lattice. These can be easily localized to the vertices of the spatial lattice by choosing the gauge parameters (an assignment of one group element to every vertex in the spatial lattice) such that all group elements are trivial except for one vertex . We will call the corresponding operators where now denotes an element in the gauge group (and not in the direct product ). These can be expressed by the left and right translation operators (59)
Physical states have to satisfy^{27}
for all vertices and all group elements . As the operators define a representation of the group, we can restrict to the elements of a generating set of the group.

This suggests considering the action of these “star” operators in the spin network basis: where on the right-hand side edges are the ones starting at , are those ending at and are all of the remaining edges in the spatial lattice. From this expression, one can conclude that the spin network states transform at in a tensor product of representations where is the dual representation to . Gauge-invariant states can be constructed by considering the projections of these representations to the trivial ones, or equivalently by contracting the matrix indices of the states with the intertwiners between the tensor product of “incoming” representations and the tensor product of “outgoing” representations at the vertices; see Section 6.2.

For the Abelian groups , the irreducible representations are one-dimensional; hence, we can omit the indices in the spin network basis. The tensor product of representations (79) is also one dimensional and equal to the trivial representation provided that for the representation labels of the outgoing and incoming edges at , respectively. Hence, we recover the Gauss constraints from the spin foam representation (12). Here, these Gauss constraints appear as based on vertices as opposed to edges as in (12). But the spatial vertices represent the timelike edges, and the Gauss constraints in the canonical formalism are just the Gauss constraints associated to the timelike edges in the spin foam representation (12).

Let us turn to the other projector . It maps to states for which all of the spatial plaquettes have trivial holonomies. Hence, the conditions on physical states can be written as where we have introduced the plaquette operators (corresponding to the flatness constraints) Here, should be a faithful representation; otherwise, one might find that also states with local non-trivial holonomies satisfy the conditions (81); see, for instance, the discussion in [137]. On the other hand, this can be taken as one possible generalization of the model. For the non-Abelian groups, the plaquette operators additionally depend on the choice of a vertex adjacent to the face, at which the holonomy around the face starts and ends. However, the conditions (81) do not depend on this choice of vertex: if a holonomy around a face is trivial for one choice of vertex, then it will be trivial for all other vertices in this face, as these holonomies just differ by a conjugation.

The *BF*-theory in any dimension is a topological field theory; that is, there are only finitely many physical degrees of freedom depending on the topology of space. Consequently the number of physical states, or equivalently of equivalence classes of kinematical states in the sense of (75), is finite and does not scale with the lattice volume.

There are a number of generalizations one could consider. First, one can couple matter in the form of defects to the models. In 3D, the *BF*-theory corresponds to a first-order formulation of 3D gravity. Point-like particles can be coupled to the model, see, for instance, [134–136], and lead to the violation of the flatness constraint (82) (coupling to the mass of the particles) and to the violation of the Gauss constraints (77) (coupling to the spin of the particles) at the position of the particle. The defects can be understood as changing the topology of the spatial lattice, that is, changing its first fundamental group. To have the same effect in, say, four dimensions, one needs to couple strings; see [138–140].

In the context of quantum computing [54], one does not necessarily impose the conditions (77) and (82) as constraints but as characterizations for the ground states of the system. Elementary excitations or quasi-particles are states in which either the flatness or the Gauss constraints are violated. These excitations appear in pairs (at least for the Abelian models [141]) and can be created by so-called ribbon operators [54, 141–143], which in the context of the proper particle models in 3D gravity are related to gauge-invariant “Dirac” observables [144]. For further generalizations based on symmetry breaking from to a subgroup of , see, for instance, [141].

In 4D gravity is not topological, rather there are propagating degrees of freedom. Similar to the discussion for the partition functions in Section 4.1, one can try to construct new models which are nearer to gravity by breaking down the symmetries of the *BF*-theory. In 3D, the flatness constraints are defined on the plaquettes of the lattice. Constraints act also as generators of gauge transformations (indeed the conditions (77) and (82) impose that physical states have to be gauge invariant), and the flatness constraints can be interpreted as translating the vertices of the dual lattice in space time [39, 40, 145, 146], which can be seen as an action of a diffeomorphism. In 4D, the same interpretation holds only for some exceptional cases corresponding to special lattices that do not lead to space time curvature; see the discussion in [107]. In general, the flatness constraints rather generate translations of the edges of the dual lattice [147].

A possible generalization leading to gravity-like models is therefore to replace the flatness conditions (81) based on plaquettes with some contractions of these conditions, such that these new conditions are based on the -cells, that is, cubes in a hyper-cubical lattice. This would correspond to the contractions of the form (the Hamiltonian constraints) and (diffeomorphism constraints) in the “complex” Ashtekar variables of the curvature with the electric flux variables [11, 12, 148–150]. The difficulty in constructing such models is consistency of the constraint algebra, that is, to find constraints that form a closed algebra. However to consider such models for finite groups should be much simpler than in the case of full gravity. Even if such consistent constraint algebras cannot be found, the physical Hilbert spaces could be constructed, for instance, with the techniques in [123–125, 151].

Furthermore, it will be illuminating to derive the transfer operators for the constrained models introduced in Section 5.3; as in the full theory, the relation between the covariant models and the Hamiltonians in the canonical formalism is still open [152, 153]. To this end, a connection representation [91, 109] can be employed.

#### 7. Tensor Models and Tensor Group Field Theories

Tensor models are theories of random topological spaces in the fashion of matrix models [154, 155], of which they may be seen as a superset.

Matrix models are a well-studied theory that can describe a multitude of different scenarios including gravity with cosmological constant [156–159] (optionally coupled to the Ising/Potts matter [160, 161] or the Yang-Mills matters theories [162]), certain string theories [163], the enumeration of virtual knots and links [164–167], and the list goes on. Moreover, they have inspired the development of a plethora of useful techniques to solve and analyse statistical ensembles of random matrices [168] such as the topological expansion [169–172], the eigenvalue method [173], the double-scaling limit, the method of orthogonal polynomials [174], and the character expansion [175, 176], to name just a few. Tensor models are a natural extension of this idea to higher dimensions. The ambition is to develop a similar technology, with similar applications, for tensor models in general.

Despite some initial pessimism [177], a recent spike of optimism occurred within the growing tensor model community, when it was discovered that a large class of such theories [18, 178–182] possessed a -expansion [183–189]; that is, their Feynman graphs could be partitioned into manageable subsets, using some parameter . In fact, it was shown that the leading order sector, in the large- limit, contained an infinite but manageable number of the Feynman graphs with the topology of the -sphere ( being the rank of the tensor). This leading order sector was analysed for a variety of models, which included the plain model [190, 191], placing the Ising/Potts [192] and dimer [193, 194] matter interactions on these -dimensional topological spaces, as well as dually weighing the spaces [195]. Critical exponents were extracted, indicating that this sector of the theory displayed identical physical properties to those of the branched polymer phase of the Euclidean dynamical triangulations [196, 197]. This likelihood was further confirmed when it was shown that their respective Hausdorff and spectral dimensions coincided [198]. Interestingly, aspects of universality have also been displayed by this leading order sector [199], while the quantum symmetries have been catalogued at all orders and take the suggestive form of a Virasoro-like algebra [200–202].

While the majority of the results stated above has been explicitly detailed for the simplest class of tensor models, the independent identically distributed IID tensor models, the initial result on the -expansion applies to many more classes of tensor models, including certain topological and quantum-gravity-inspired tensor models. Thus, a current aim is to develop the above formalism for these broader tensor model scenarios. As stated in the introduction, spin foams for quantum gravity generically involved the Lie groups, and so more structure, than that provided by tensor models, is needed to describe them.

Tensor group field theories (TGFTs) [19–21, 203] are to quantum field theory as tensor models are to quantum mechanics. Spin foam diagrams naturally arise as the Feynman graphs of TGFTs [1–5], and indeed various quantum-gravity-inspired TGFTs exist in the literature, such as the Boulatov-Ooguri theories [204, 205], the EPRL and FK theories [22–24, 206], and the Baratin-Oriti theories [25, 26]. As a result, they have garnered increasing attention from the quantum gravity community, since, inherently, all current spin foam models for 4d quantum gravity involve a truncation of the degrees of freedom (due to the use of a single lattice structure) and a manifest loss of diffeomorphism symmetry. With their explicit summation over -dimensional topological spaces, the hope is that TGFTs recapture these lost degrees of freedom. Indeed, some evidence has already been found at the level of the TGFT action for a seed of diffeomorphism symmetry [207]. Moreover, many techniques may be applied to these field theories that are idiosyncratic to quantum field theories (as opposed to quantum mechanics). Perhaps the most striking is the study of their renormalisation group properties [208–218], but also work has commenced on the study of mean field theory properties [219], matter coupling [220–222], various symmetry analyses [223, 224], and instantonic field theory solutions [225–228] (including cosmological applications [229]).

Having said all that, the aim of this paper is not to detail spin foam models with the Lie groups, but rather spin foam models with finite groups. This places us back under the purview of tensor models, and it is these that we will describe in the coming section.

To commence, let us detail some of the general setup.

##### 7.1. General Setup

Let us construct the class of *independent identically distributed* (IID) *model*s. We will attempt to be precise without being especially detailed. We refer the reader to [230] for a more thorough explanation. The fundamental variable is a complex rank- tensor (), which may be viewed as a map , where the are complex vector spaces of dimension . It is a tensor, so it transforms covariantly under a change of basis of each vector space independently. Its complex-conjugate is its contravariant counterpart.

One refers to their components in a given basis by and , where , , and the bar (−) distinguishes contravariant from covariant indices. We have purposefully denoted the indices by as they may be viewed as elements of respective .

As one might imagine, with these ingredients, one can build objects that are invariant under changes of bases. These so-called *trace invariants* are a subset of -dependent monomials that are built by pairwise contracting covariant and contravariant indices until all indices are saturated. It emerges readily that the pattern of contractions for a given trace invariant is associated to a unique *closed *-*colored graph*, in the sense that, given such a graph, one can reconstruct the corresponding trace invariant and vice versa.^{28}

In Figure 5, we illustrate the graph , the unique closed -colored graph with two vertices, which represents the unique quadratic trace invariant .

More generally, we denote the trace invariant corresponding to the graph by .

From now on, we will make two restrictions: (i) all of the vector spaces have the same dimension and (ii) we consider only *connected* trace invariants, that is, trace invariants corresponding to graphs with just a single connected component.

Given these provisos, the most general invariant action for such tensors is
where is the set of connected closed -colored graphs with vertices, is the set of coupling constants, and is the *degree* of (see [18] for its definition and properties). This defines the IID class of models.

##### 7.2. -Expansion

The central objects for further investigation are the *free energy* (per degree of freedom) associated to these models:
along with the various other -point Green functions. When facing such a quantity, the standard procedure is to expand it in a Taylor series with respect to the coupling constants and to evaluate the resulting Gaussian integrals in terms of the Wick contractions. It transpires that the Feynman graphs contributing to are none other than *connected closed **-colored graphs* with weight:
where
is a symmetry factor, and runs over the subgraphs with colors of , while and are the vertex and face sets of , respectively. For , the *degree * coincides with the genus of a surface specified by . A few more words of explanation are most definitely in order here. The graphs arise in the following manner. One knows that a given term in the Taylor expansion is a product of trace invariants upon which one performs Wick contractions. For such a term, one indexes these trace invariants by ; that is, we index their associated -colored graphs . A single Wick contraction pairs a tensor , lying somewhere in the product, with a tensor lying somewhere else. One represents such a contraction by joining the black vertex representing to the white vertex representing with a line of color . Thus, a complete set of the Wick contractions results in a connected (as one is dealing with the free energy) closed -colored graph. A particular Wick contraction is drawn in Figure 6.

It requires a bit more work to reconstruct the amplitude explicitly; see [230]. Importantly, is a nonnegative integer, and so one can order the terms in the Taylor expansion of (84) according to their power of . Quite evidently, therefore, one has a *-expansion*.

##### 7.3. Interpretation

Strikingly, -colored graphs represent -dimensional simplicial pseudomanifolds. We will give a rough presentation here. One distinguishes the -bubbles of species as those maximally connected subgraphs with the colors . These -bubbles are identified with the -simplices of the associated simplicial complexes. Moreover, one can see clearly that the -bubbles of species are nested within -bubbles of species , where . This set of nested relationships encodes the gluing of the simplices within the simplicial complex. This argument may be made rigorously. Thus, at the very least, rank- tensor models capture a sum over -dimensional topological spaces.

Moreover, a geometrical interpretation may be attached most readily to the amplitudes of the IID model by setting for all , where is some coupling constant. Then, one may rewrite the amplitudes as where denotes the number of simplices in the simplicial complex associated to and Interestingly, the discretization of the Einstein-Hilbert action on an equilateral -dimensional simplicial complex takes the form where and are the bare Newton and cosmological constants, respectively, and is the volume of the equilateral -simplex, while denotes the deficit angles associated to the -simplices. If one further associates to in the following fashion: then one finds that the Feynman amplitudes for the IID model are coincide with those in a Euclidean dynamical triangulations approach. We will return to this later.

##### 7.4. Large- Limit

In the large- limit, only one subclass of graphs survives, containing those graphs with . At this point, there is a marked difference between two and higher dimensions.(i)In the 2-dimensional model, is the genus of the graph in question. Thus, the graphs surviving in this limit are all graphs with the topology of the -sphere.(ii)In higher dimensions, that is , it was shown in [190, 191] that the only graphs surviving this limiting procedure are the melonic graphs (with this name stemming from their distinctive structure). While these have the topology of the -sphere, they do not constitute all possible -colored graphs with this topology.

The series constituting the leading order sector has a finite radius of convergence. This indicates that the theory displays critical behaviour, for some values of the coupling constants, characterised by some critical exponents. There are various scenarios that one might investigate. One simple, yet interesting case is to set for all , where is some coupling constant. In this scenario, the series display the following critical behaviour: where Multicritical behaviour can be extracted by tuning several coupling constants independently.

Given the description provided in the previous subsection, one notices that the large- limit corresponds to . Moreover, tuning , one enters the regime dominated by graphs with large numbers of vertices (or equivalently, simplicial complexes with large numbers of simplices). As it stands, this corresponds to a large-volume limit. However, with some more work, one can interpret it as a continuum limit. To begin, the average volume is To obtain a continuum limit, one should tune while keeping the average volume finite. This may be achieved in the following fashion: while keeping where is a renormalized cosmological constant.

##### 7.5. Coupling to the Ising and Potts Matter

The coupling to the Ising/Potts matter in tensor models is a direct generalisation of that scenario in matrix models [231, 232]. For the Ising matter, one considers a model with two complex tensors and action: where and is a coupling constant. This class of models has been examined in [192], where critical exponents have been extracted, which coincide with those of branched polymers.

This matter model may be extended to the -state Potts matter by increasing the number of tensors along with a suitable alteration of the propagator.

##### 7.6. Non-IID Models

The IID models are but the simplest of a much larger class of models possessing a -expansion, defined by actions of the form
The difference resides in the propagator and the scalings, , of the coupling constant, which can be chosen to reproduce any local (quantum-gravity-inspired) spin foam model (in this case, with the finite rather than the Lie group information). As an example, let us briefly detail the choice that yields the Boulatov-Ooguri class of tensor models. Firstly, the propagator is chosen to be
The can be specified so that the resulting amplitudes for the free energy take the form
where
Given that and are the edge and face sets of , respectively, while is the respective orientation of the edge lying in the boundary of the face , it is clear that is the holonomy around the face . As a result, the amplitude within square brackets is the *BF*-theory amplitude associated to the topological manifold represented by