# Renormalization of Group Field Theories for Quantum Gravity: New Computations and Some Suggestions

^{1}Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Potsdam-Golm, Germany^{2}Institute for Physics, Humboldt-Universität zu Berlin, Berlin, Germany^{3}Arnold-Sommerfeld-Center for Theoretical Physics, Ludwig-Maximilians-Universität, München, Germany

We discuss motivation and goals of renormalization analyses of group field theory models of simplicial 4d quantum gravity, and review briefly the status of this research area. We present some new computations of perturbative Group field theories amplitudes, concerning in particular their scaling behavior, and the numerical techniques employed to obtain them. Finally, we suggest a number of research directions for further progress.

## 1 Introduction

Group field theories (GFT) [1–3] are quantum field theories which aim at describing the fundamental quantum structures that constitute spacetime. They are quantum field theories *of* spacetime, rather than *on* spacetime. They are defined on group manifolds (hence the name), with an associated phase space given by the cotangent bundle of the same group.

Equivalently, in conjugate variables, the basic field maps *d* copies of the Lie algebra of the same group to the complex numbers, and should be understood in general as a non-commutative function (i.e., an element of a non-commutative algebra of functions), since the Lie algebra is in general a non-commutative manifold and this reflects on the algebra of functions defined on it. Depending on the specific model one is considering, various restriction can be imposed on the field, its domain, its target, and of course the choice of group manifold and ‘dimension’*d* are also model-dependent. What is general, in all current GFT models, is that the basic quanta of the theory, corresponding to the basic field excitations, can be represented as abstract cells or polyhedra with the *d* algebraic data forming the domain of the field associated to their (boundary) faces. When *d* is chosen as the dimension of the spacetime to be reproduced in some approximation, the corresponding GFT quanta can be understood as *d*

and then by its partition function (assumed here as being of statistical form):

Beside a quadratic local term (defined by an integral kernel convoluted with the two fields), the GFT action is determined by interactions (just one in the above example) that possess a characteristic ‘combinatorial non-locality’in that the interaction kernels pair non-locally the field arguments (*d* variables being contributed by each field entering the given interaction term). Interaction kernels of order *n* can be associated with possible ways of gluing together *n**d*-dimensional cell. The specific combinatorial patterns (i.e., the specific cells being associated to each fundamental interaction) and precise form of the interaction and kinetic kernels are part of the definition of each particular GFT model. However, from this generic aspect of the formalism follows one key fact: GFT Feynman diagrams

obtained gluing interaction vertices (*d*-cells) along their *d*-dimensional cellular complexes, of arbitrary topology (since a priori there is no restriction on the allowed gluing). The Feynman amplitudes assign a probability amplitude for each such cellular complex, seen as an elementary interaction process of the fundamental GFT quanta.

The above can be taken as a sketchy definition of the TGFT formalism, but it is of course the specification of particular models which gives tentative physical meaning to it and this meaning will therefore change in different contexts. In particular, when the physical interpretation of a given model is grounded in its perturbative expansion, it will affect what we expect about the properties of the Feynman amplitudes

Before doing so, let us spend some words to clarify the choice of a statistical form for the GFT partition function. The foundations of the GFT formalism have received some attention only recently, and much remains to be understood. Both statistical and complex weights for GFT fields summed over in the definition of the partition function can be found in the literature, but in fact most of the literature until recently has been focusing only on the perturbative expansion of the models, where the non-perturbative definition is less directly relevant, thus avoiding the need to investigate it. In particular, for the class of models to be analyzed in the next subsection, the interpretation in terms of quantum sums over discrete geometries (and topologies) is guaranteed by the form of the Feynman amplitudes (which are either complex or real with oscillating behavior, in line with their interpretation as defining a quantum discrete path integral for gravity on a lattice, as we are going to discuss). As for the non-perturbative definition of the GFT dynamics, instead, we have less guidelines, especially for the general definition of the formalism, before considering specific models. In the context of tensor models and tensorial group field theories treated as a generalization of usual local QFT framework, all non-perturbative analyses have simply assumed from the start a statistical definition, taken to be primary and in no special need of further justification; in particular, no attempt to understand it as the counterpart of an operatorial definition of the same model has been made. In non-quantum gravity applications, this is simply a choice of a ‘classical’statistical field theory context, usually adopted for purely mathematical reasons (better chance of making sense of the path integral).

For quantum gravity applications, one may be interested also in a better conceptual foundation of the formalism. In this context, a complex weight involving the GFT action as a pure phase is of course a possibility, and maybe more in line with the intuitive idea of a quantum gravity path integral (even though the quantum gravity interpretation of the GFT field itself is not straightforward, while as we remarked already, this interpretation be consistently associated to the discrete structures appearing in the GFT perturbative expansion). However, we should remind ourselves that these quantum gravity models should not be expected to encode any global unitary evolution, as it is in fact true for any fundamental quantum gravity dynamics, and this removes one strong motivation to insist on this formulation of the path integral. Nor we have a complete derivation of the GFT path integral from a canonical quantum gravity dynamics or, alternatively, from some formal field theoretic gravitational path integral, that could dictate one choice over another. To date, the only tentative derivation of a GFT partition function ‘from first principles’, was given in ref. 4 where it was seen to arise from a *quantum* statistical definition of equilibrium for a system of quantized simplices (indeed, the basic quanta of such GFT models), under a requirement of maximization of entropy and a choice of macroscopic conditions to be imposed on average (this choice concurs to the specification of the resulting GFT action). If among these constraints one includes some appropriate counterpart of the Hamiltonian constraint of canonical gravity (adapted to the discrete setting), then the GFT partition function can be seen as a sort of grandcanonical partition function relaxing the imposition of such constraint (in the sense that configurations satisfying the constraint are assigned greater weight, but fluctuations off the constraint surface are allowed) as compared with a ‘microcanonical’ensemble in which only solutions of the constraint are allowed. From a canonical quantum gravity perspective the latter, more restricted case would correspond to a definition of the physical inner product between quantum gravity states (with appropriate insertions of observables inside the GFT partition function), while in general the GFT formalism deals then with a broader class of quantum amplitudes. This scenario was also anticipated (more formally) in ref. 5, discussing the relation between GFT and canonical loop quantum gravity. This perspective also resonates (with many details still to be clarified, though) with the presentation of GFT from a quantum gravity perspective in ref. 6, where the ‘tree level’(thus dominant, from the perturbative point of view) GFT amplitudes were suggested to define the physical inner product of canonical quantum gravity, while the remaining GFT configurations were associated with topology changing processes (off-shell, from the canonical quantum gravity perspective).

### 1.1 Group Field Theories, Spin Foams and Other Quantum Gravity Formalisms

We recognize in this brief outline the straightforward generalization of how 2d surfaces are generated in the perturbative expansion of random matrix models. Indeed, GFTs can be seen as group-theoretic enrichment of random *tensor* models [7–9], to which they reduce if the Lie group domain is replaced by any finite set of *N* elements. The Feynman amplitudes become purely combinatorial, but the type of diagrams remains the same. Seen as tensors, GFT fields admit a natural action of unitary (and orthogonal) groups on their arguments. If one requires GFT interactions to be invariant under such unitary transformations, they can be fully classified, and we speak of *tensorial* GFT models. Most of the literature on GFT renormalization [10, 11] concerns these tensorial GFT models. More generally, focusing on tensorial aspects of GFTs allows to gain a greater control over the combinatorial structures of their states, diagrams and amplitudes and many of the results obtained in the simpler context of tensor models apply also to GFTs: the use of colors to encode the topology of Feynman diagrams, the large-N expansion, double scaling limits, universality results etc. The first two, in particular, are crucial for GFT renormalization.

In this contribution, we focus on GFT models which are ‘quantum geometric’: their fundamental quanta are quantized tetrahedra with a quantum geometry encoded in group-theoretic data. More precisely, the classical phase space of a single Lorentzian tetrahedron in 4d is chosen to be the cotangent bundle of 4 copies *G* is one of the chosen groups mentioned above. More precisely, the natural Hilbert space for a single tetrahedron in this class of models would be

For this class of geometric models, the GFT formalism benefits from direct links to other modern quantum gravity approaches, which can, viceversa, benefit from GFT tools and results.

First, of all, when

Next. for this class of models, the GFT Feynman amplitudes take the form of (non-commutative) simplicial gravity path integrals [16, 17], when written in (non-commutative) Lie algebra variables, which encode the discrete metric. The group variables, on the other hand, are understood as encoding the discrete gravity connection. They correspond indeed to discretizations of a classical formulation of gravity as a topological BF theory with added geometrician constraints, on the simplicial complex dual to the GFT Feynman diagrams. The specific way in which the BF action is discretized depend on the quantization map applied to Lie algebra variables, and different models correspond to different strategies for the imposition of the constraints and path integral measures.

In fact, when the same Feynman amplitudes are recast as functions of group representations, using Peter-Weyl or Plancherel decomposition, they take the form of spin foam models [18]. Spin foam models have been introduced as a covariant language for computing spin network dynamics, so they can be understood as a covariant counterpart of canonical loop quantum gravity. A second perspective is to see spin foam amplitudes as a purely algebraic version of lattice gravity path integrals, or state sum models. In GFT, they arise as Feynman amplitudes. The correspondence is generic: for any given set of spin foam amplitudes associated to simplicial complexes (and admitting a local decomposition with respect to the complex), one can find a GFT action such that the perturbative expansion of the quantum partition function will produce the given amplitudes as Feynman amplitudes (and viceversa, any GFT action corresponds to a set of spin foam amplitudes). A complete definition of a spin foam model requires a prescription for the amplitudes to be associated to all possible cellular complexes (in some specified class) and an organization principle for them, i.e. one way of comparing, composing or selecting them, to obtain a single number for any observable one wants to compute. The GFT embedding provide one such clear organizing principle, by summing them in a QFT perturbative expansion. In addition, it provides a whole set of QFT tools that can be applied to study their mathematical foundations as well as for extracting physics. GFT renormalization can be seen, indeed, from this spin foam perspective.

### 1.2 Simplicial Group Field Theory Models for 4D Quantum Gravity

The starting point for the construction of simplicial GFT (and spin foam) models of 4d quantum gravity is the quantum geometry of a single tetrahedron in 4d [19].

The quantum geometry of this basic building block, and the extended structured built from it, can be described in various parametrizations [20, 21], and a number of generalizations can also be defined [22, 23] and imported in the GFT framework. Classically, one can use two equivalent characterization of a tetrahedral geometry, leading immediately to an algebraic translation. First, one can start with assigning four vectors *V* (thus satisfying

An equivalent encoding of the classical geometry of a single tetrahedron uses directly the variables of discretized topological BF theory. All geometric quantities of a single tetrahedron can be computed starting from four bivectors *V*. The phase space of a single tetrahedron can be taken to be the cotangent bundle

A spin foam amplitude, that is a GFT Feynman amplitude written in representation variables, will be assigned to any given simplicial complex, dual to a GFT Feynman diagram. The basic building block is an assignment of a quantum amplitude to each 4-simplex, i.e., a ‘vertex’of the spin foam complex given by the GFT Feynman diagram, with this amplitude function of the algebraic data associated to the five tetrahedra on its boundary. These boundary data can be written as

One last comment about the discrete geometry of these models. The construction sketched above, at the classical level, leads to a full characterization of the discrete geometry of the 4d simplicial complex (to which the spin foam amplitudes are associated), equivalent to the more standard characterization in terms of edge lengths, as used in Regge calculus, even though it uses a different set of classical variables (it corresponds, indeed, to a formulation of classical simplicial geometry in terms of the discrete counterpart of the variables of BF theory, suitably constrained, or to so-called ‘area-angle’Regge calculus [25]). The translation of the same characterization at the quantum level, and in particular the correct imposition of the geometrician constraints on quantum states and amplitudes, is the crucial point for ensuring the correctness of the model from a discrete geometric point of view, and it is still subject to debate in the literature. In particular, one would expect to find back the Regge action for metric (edge length) variables, or an equivalent classical reformulation, in the semi-classical expansion of the spin foam amplitudes for a generic simplicial complex (or of the corresponding simplicial path integral). Many results are available (see the cited references) on this issue for a single 4-simplex and for extremely simple complexes, at least for the EPRL model, but the results are mixed, and the situation is especially unclear for larger complexes.

The general formula for the spin foam amplitudes, for all the models in this class, in the Riemannian setting, for given cellular complex

The coefficients *f* are matrix elements of the map between Spin

where we have a Spin*C* is the 3j-symbol, and the 15j-symbol is the one of the first type. We have indicated with *w* the function of group representations that characterizes the implementation of the simplicity constraints defining each model, depending on the representations of

The amplitude’s formula can be rewritten in terms of propagators as follows:

where the propagator

In order to derive the master Integral expression for a GFT Feynman

The models we will deal with in the following are the EPRL model [18] and the Duflo BO model [17], whose defining maps are:

where the *T* function is given by:

but it can also be given an expression purely in terms of representation labels. See ref. 17 for more details.

We note that the relative simplicity/complexity of these two models is highly dependent on the basis in which they are expressed, with the flux representation switching such relative complexity with respect to the spin representation given above.

We also point out that other spin foam models, obtained from alternative strategies of imposition of the same geometrician constraints and thus also belonging to the same general class we are considering, can be cast in principle in the same general form, and studied by the same method we will illustrate in the following. Beside difficulties, for some of them, in achieving an explicit and manageable expression for their corresponding *w* coefficients in representation variables, that makes the analysis more cumbersome, it would indeed be very interesting to perform the same scaling analysis of amplitudes and compare with our results.

More details about the construction of spin foam amplitudes, as well as all the ingredients we mentioned as entering in such construction, in a language well adapted to their GFT embedding, can be found in ref. 17.

## 2 Renormalization of Group Field Theories for 4D Quantum Gravity

Let us now discuss motivation and current status of renormalization of simplicial GFT (and spin foam) models for quantum gravity.

Beyond the connection to spin foam models and simplicial gravity path integrals, the general strategy for renormalization of GFT models [10, 11, 26] is to treat them as ordinary QFTs defined on a Lie group manifold, thus using the group structures (topology, Killing forms, etc) to define ‘scales’and mode integration. A natural notion of scale, to be used to label the RG flow, is provided by group representations, which index the spectrum of differential operators on the group, e.g., the Laplace-Beltrami operator, in turn often used to define the propagator of GFT models. Cut-offs imposed as part of a renormalization group scheme are then imposed on representation labels; for example, in the case of

For a proper renormalization group scheme, however, two more ingredients are needed: control over the theory space corresponding to a given GFT model, i.e. the space of allowed interactions; a detailed characterization of the combinatorics of (the cellular complexes dual to) GFT Feynman diagrams. On neither of these two points much is known for simplicial 4d gravity models. As a result, most work in the context of GFT renormalization has been done focusing on tensorial GFT models, where the above limitations are not present.

Before discussing the goals of GFT renormalization, we spend a few words of caution concerning the physical interpretation of the renormalization group scheme and derived flows. With scales associated to group representation labels, the natural cutoffs entering as UV cutoffs are for large representations. The associated RG then flows from large to small representations (from UV to IR). In LQG and simplicial quantum geometry representation labels identify eigenvalues of geometric operators (e.g., triangle areas or tetrahedral volumes). Large representation labels correspond to large values of such geometric quantities. Thus we have an apparent inversion of roles here, with large distances/volumes playing the role of UV scales in GFT. Caution however should be exercised. In both LQG and simplicial geometry, we know an area of a surface, say, to result roughly speaking from the sum of the individual areas of all elementary surfaces forming the one under consideration, so that one has

Also, let us comment on the importance of a better understanding of the symmetry properties of these 4d gravity models, both for the characterization of the corresponding theory space and for their relation to continuum gravity, which is of course crucially characterized by diffeomoprhism symmetry. The issue of symmetries in GFT models is very important but also very much open. At the general level, we do not know much about symmetries of 4d GFT models, beyond the Lorentz invariance of the kinetic and interaction kernels and of the Feynman amplitudes (implemented as in usual lattice gauge theories, since the amplitudes are in fact lattice gauge theories for a (constrained) Lorentz connection). Because of the simplicity constraints and also of the simplicial combinatorics that characterize them, moreover, even the tensorial symmetry typical of tensor models is not present (or at least not manifest). Moreover, even for the few symmetries we know of, in other models, the analysis of their consequences, for example in terms of conservation laws, is complicated by the non-local nature of the GFT interaction (see the analysis [27, 28]). Concerning diffeomorphisms, strictly speaking (being smooth transformations) they are not defined in a discrete context like that of GFT Feynman amplitudes, i.e. spin foam models and lattice gravity path integrals, and thus the question becomes whether we can identify some analogue of diffeomorphic symmetry that, in a continuum limit, could be then identified with the one characterizing GR. There are several analyses of such question for 3d (topological) models at the level of spin foam amplitudes [29], lattice gravity (see for example [30–32]) and corresponding GFT formulation [33], but nothing similar in the 4d gravity case (where the 4d counterpart of the symmetry identified in the 3d case is actually broken, at the discrete level [34]). When attempting a reconstruction of an effective dynamics of geometry in a continuum approximation, as done in the context of GFT cosmology, one has to proceed in terms of observables of the fundamental theory that have a chance to correspond to diffeomorphic invariant observables in GR, since all the structures of continuum GR on which diffeomorphisms act, e.g., manifold points, directions an coordinate functions, but also fields defined on the same manifold, are simply not present in the theory.

### 2.1 Quantum Consistency and Perturbative Renormalization

GFT models are first defined in perturbative expansion and it is in this perturbative formulation that spin foam amplitudes, and simplicial gravity path integrals, appear. The perturbative GFT amplitudes generically diverge and regularizations have to be imposed. It is this truncation that corresponds to working at a given ‘scale’. Is this definition of the quantum dynamics of GFT models consistent? is the spin foam description consistent? Here, consistency means first of all valid for all ranges of dynamical variables, under (controlled) removal of regulators. If not, the GFT model as defined in perturbative expansion, and thus the corresponding spin foam model (and simplicial path integral) cannot be trusted. In the GFT language, this is recognized immediately to be the issue of perturbative renormalizability of a given model. We should only trust, from the spin foam or lattice gravity point of view, only GFT models that turn out to be (perturbatively) renormalizable.

We note in passing that there should be no requirement that the model is finite (in the sense of presenting no divergence even before any renormalization); first, we have no obvious reason to expect it, if the model contains an infinite number of degrees of freedom; second, renormalizable models are usually more interesting, as QFTs, than finite ones, since they have a non-trivial RG flow and new effective physics at each scale.

Let us clarify further what we mean, here, to avoid possible misunderstandings. As a general point about field theories, we are saying that finiteness of the Feynman amplitudes associated to a given subset of diagrams, or even to all diagrams involved in a given “scattering process” is not so important, per se, and in fact not necessarily desirable. What is important is that the scattering amplitudes can be -made finite- by suitable renormalization procedure (at any order in perturbation theory), if originally divergent in terms of bare couplings, and after resumming all the diagrams involved in their computation. The final renormalized scattering (or transition) amplitudes are what is physically relevant. A theory that is instead simply finite in the sense of not requiring any renormalization, even if clearly easier to deal with, would be less interesting from a physical point of view because this finiteness would probably indicate that the quantum dynamics is not very rich and it does not change much across scales (i.e. when more of its quantum degrees of freedom are accounted for). The consequence would also be a less interesting phase diagram. This should explain our comments about finiteness of GFT Feynman amplitudes. As for the finiteness of “quantum gravity scattering amplitudes”, what we should expect or desire depends on how we interpret the terms. If we take a given GFT model to be a tentative definition of full quantum gravity, then for sure we should hope that its “transition or scattering amplitudes” be finite, in the end, i.e., after renormalization. If the relation between GFT and quantum GR is as in the first case discussed above, i.e., the two are “equivalent”, then this also implies that the transition or scattering amplitudes of quantized GR will be finite, when properly defined and after renormalization (even if the renormalization procedure as well as the observables expressing the amplitudes may look very different in the two formulations). If the relation is as in the second case, and thus “quantum GR” is just an effective theory, then we do not have to expect that its transition or scattering amplitudes are finite tout court, but only within the domain of validity of the approximations or truncations leading to it within the fundamental theory.

Thus, the requirement of perturbative renormalizability is an important constraint, which helps removing from consideration inconsistent constructions. Here, the GFT embedding proves potentially very important also for spin foam models (and loop quantum gravity). All known GFT and spin foam models present several ambiguities, some intrinsic to any quantization procedure, others specific to simplicial GFT (and spin foam) models of quantum gravity. Requiring perturbative renormalizability means constraining such ambiguities. To name one, we have little constraints of the face amplitudes of spin foam models, even though they can drastically affect the scaling behavior of the GFT and spin foam amplitudes, to the point of allowing to achieve perturbative finiteness easily by simply fixing them to this end (which also shows why finiteness *per se* cannot be a goal, without a proper physical understanding), while perturbative renormalizability is a much trickier requirement. For example, see how simple modifications of the Barrett-Crane model (which is also the limit of both EPRL model and Duflo model for infinite Immirzi parameter) affect the resulting amplitudes [35, 36].

Let us list some of them. A first one is combinatorial: why restricting to simplicial complexes? These are the ones for which we have a better understanding of the discrete geometry underlying our models, and in particular of the simplicity constraints that characterize them. But what other cellular complexes should be included in the theory for consistency, e.g., because corresponding to the counterterms required for taming the perturbative divergences? Others concern the underlying quantization and imposition of simplicity constraints. Being functions of the flux variables (which are non-commutative), they depend on which quantization map is chosen to quantize such variables. Different choices result in different discrete gravity actions and different simplicial path integral measures, thus different spin foam amplitudes. Also the very definition of the simplicity conditions as operator equations acting on quantum states depends on the chosen quantization map, from which follow thus different constraints on representation variables in the spin foam amplitudes. Further, the strategy by which simplicity constraints are imposed produces in general different models or versions of the same type of models (this is apparent in the Riemannian case, while in the Lorentzian one we only have experience of different versions of the EPRL model). These and other ambiguities are discussed, e.g., in ref. 17. Using

So, where do we stand, on this important issue? For simplicial GFT models of 4d quantum gravity the answer is, unfortunately, that we are only at the very beginning. The main reasons have been already mentioned. First, we do not know enough of their symmetries to characterize the relevant theory space. Second, the amplitudes for these models are very involved and technically challenging to compute, mostly due to the fact that the imposition of simplicity constraints makes them defined not simply on Lie group manifolds but on particular sub-manifolds of these (usually not even corresponding to homogeneous spaces). Third, dominant configurations (i.e., those giving the most divergent contribution to the amplitudes) are not just flat connections or similarly simple, but correspond to richer configurations from the point of view of simplicial geometry; possibly, they correspond to (or possibly include) the whole set of Regge geometries found as saddle-point configurations in the asymptotic analysis of spin foam amplitudes and corresponding simplicial path integrals. Therefore even power counting results are hard to obtain, and the brute force analysis of divergences is not advanced enough to indicate the needed counterterms, forming the theory space. If the theory space is hard to characterize also in the simpler 3d simplicial case (corresponding to topological BF theory), at least the amplitudes are manageable enough to obtain complete power counting theorems [38], identify some counterterms [39] and nice finiteness results [40].

So computational challenges are one big obstacle. It is on this aspect that we focus in the next section, presenting some new results in the Riemannian context. These new results should be added to other ones we have on the calculation of radiative corrections and basic divergences of both Riemannian and Lorentzian simplicial spin foam models, and on explicit evaluations of their building blocks (mainly the vertex amplitudes). For a partial list, see [41–49] and references therein. Future progress will build on these hard-won calculations. In turn these results build on the hard-won understanding, based on both analytical and numerical studies, regarding the asymptotic properties of SU

For a comparison, one has to look at the amount of knowledge we have accumulated on tensorial (thus colored) GFT models [10, 11, 26]. Here we know several (classes of) models which are rigorously proven to be perturbative renormalizable, comprising both abelian and non-abelian models, on homogeneous spaces, with or without gauge invariance (closure condition), in different dimensions. Divergences are associated to bubbles, i.e. cells of the complex dual to the cellular complex associated to a GFT Feynman diagram, and typically the most divergent diagrams that form the relevant theory space of renormalizable theories are melonic ones, also singled out in tensor models. However, we also know example of TGFT models which are renormalizable outside the melonic truncation [62], and these examples may be relevant also for the case of simplicial GFT models, since the structure of their divergences presents some aspects of the simplicial case.

### 2.2 Continuum Limit and Non-Perturbative Renormalization

GFT models of quantum gravity are bona fide QFTs, thus they possess infinite degrees of freedom, as we expect quantum gravity to do (at least thinking of it naively as a quantum theory of the gravitational field). Control over a very large number of degrees of freedom can only be achieved step by step, within some truncation scheme. With the inclusion of more and more degrees of freedom, we can expect a richer and richer set of new phenomena to be unraveled, simply because the physics of many (quantum, interacting) degrees of freedom is very different from that of few of them. In particular, we expect new phases to be revealed. Controlling the full quantum dynamics is controlling the continuum limit of GFT models, and this implies mapping out as best as we can the phase diagram of the same models. In practical terms, it means being able to evaluate the full GFT partition function, for given values of coupling constants. This is the problem of computing the full non-perturbative renormalization group flow of any given GFT model.

Given the mentioned structural connections, understanding the non-perturbative renormalization of a quantum gravity GFT model implies controlling the continuum limit of the corresponding lattice gravity path integral and spin foam model, and the full quantum dynamics of the corresponding canonical loop quantum gravity formulation. The characterization of the continuum quantum gravity phase diagram and the identification of one phase where an effective general relativistic dynamics of spacetime can be extracted is in fact the key outstanding open issue in the field [63–66].

This should already make clear why the precise relation between the (non-perturbative) renormalization of GFT models for 4d quantum gravity and the (non-perturbative) renormalization of continuum quantum GR treated as an ordinary field theory (as in the asymptotic safety approach) can only be envisaged in a very tentative manner. Let us give only some comment on our own tentative perspective on this. take a given GFT model that can be fully defined at the non-perturbative level, thus associated with a continuum phase diagram where RG flow trajectories are well-defined from the deep UV (in the GFT sense) to the full IR (still in the GFT sense), and thus accounting for all the (infinite) degrees of freedom of the model; to achieve this situation is the goal of non-perturbative GFT renormalization, as explained. A matching with GR requires that one can also compute, in the same model, observables which characterize fully a 4d geometry and that can be shown to satisfy the GR equations in a classical approximation. Now, we can envisage two possibilities. If the rewriting is, in the appropriate sense, exact, i.e., if one can in principle go from the GFT formulation of the theory to the geometric “quantum GR” one, in the same continuum limit, then the GFT model could be seen, in fact, as a definition of “quantum GR”, without any change in dynamical degrees of freedom. In this case, one could expect that there exist a translation of the RG picture of the given GFT model into the one obtained by a non-perturbative RG treatment of GR, for example as provided (ideally) by the asymptotic safety scenario, and an isomorphism between their corresponding phase diagrams and RG trajectories. If the rewriting requires, instead, some truncation of the dynamical degrees of freedom of the GFT model, is valid only for a subset of the GFT observables, or some other drastic approximation to be valid, i.e., if “quantum GR” turns out to be only an effective, emergent description of some sector of the full quantum GFT, then the situation is different. In this case, we should not expect that the GFT phase diagram matches the GR one, and we can only expect that it will reproduce a portion of it, for scales and regime of couplings where the needed approximations and truncations hold. This regime will probably be the one corresponding to “low energies” from the standard GR and effective QFT perspective. Of course, all the above is very much tentative and it is hard to envisage the precise relation at the current stage of development of GFT as well as of “quantum GR”, even though a number of features of GFT models (e.g., the fact that they include a sum over topologies and not just geometries, at least at the discrete level) would suggest that the second scenario is more likely.Where do we stand, at the non-perturbative renormalization level? Beside work on the non-perturbative RG flow of tensor models [67, 68], a lot of activity has focused on the analysis of GFT models proper [11, 26]. Two main strategies have been followed. One is based on constructive methods, mostly focusing on the resummation of the perturbative series, e.g., showing Borel summability. The other is based on functional renormalization group analysis, either (mostly) based on the Wetterich-Morris equation for the effective action, or the Polchinski equation for n-point functions. For the same reasons that limited work on perturbative GFT renormalization, little is known about the general RG flow of simplicial GFT models of 4d quantum gravity. Simplicial GFT models in 3d have been shown to be Borel summable [69, 70] and phase transitions for the GFT formulation of simplicial BF theory in any dimension has been shown to exist [71]. But no similar analysis has been carried over to the 4d gravity case, where, as mentioned, we even lack perturbative indications.

The tensorial GFT case, on the other hand, has been widely explored, mostly via functional renormalization techniques, with many results on a variety of models, again both abelian and non-abelian, with and without gauge symmetries, based on compact as well as non-compact Lie groups, in different dimensions. Concerning UV behavior, asymptotic freedom is found in many examples and asymptotic safety is found in others [72], in various truncations, and the perturbative results have been reproduced from a non-perturbative standpoint. More results on the relevance of truncation beyond the melonic sector have been found [73], and the use of Ward identities for studying the RG flow have been explored [26]. Concerning IR behavior (i.e., the actual continuum limit), work is more limited (and more difficult) at the analytic level, but hints have been found, in various truncations and for various models, of a non-trivial phase diagram. In particular, hints of the existence of Wilson-Fisher fixed points (often found alongside asymptotic freedom in the UV) and of broken (or condensate) phases have been obtained [62, 74], indirectly supporting parallel work on the extraction of continuum gravitational physics from such condensate phases [75–77].

Even if it is unclear, at this stage, which of these results holds also in the simplicial 4d models, with their additional quantum geometric intricacies, all this work on tensorial GFT models has certainly led to a better understanding of GFT renormalization group schemes (and flows). This will certainly turn out to be useful also for the analysis of full-blown quantum gravity models.

## 3 Bubble Divergences and Radiative Corrections in Group Field Theories: Some New Results

In this section we report on some recent results concerning the leading order radiative corrections to *N*-point functions (*Duflo model* and the EPRL model, whose amplitudes we have recalled above. We refer to the cited literature for more details on motivations, construction and features of these GFT (and spin foam) models. Also, we limit our presentation to a summary of results and procedures; a more detailed presentation with be left for a forthcoming publication.

### 3.1 A Warm-Up Example: The 2-Point Function of the Ooguri Group Field Theory Model

As a warm-up, we recall the general procedure to compute the degree of divergences of GFT amplitudes both in the holonomic and spin formulation of GFT models, using the simpler case of the Ooguri GFT model for 4d topological BF theory with local

*Analytical evaluation in group variables* - For models defined on the full group manifold or a corresponding homogeneous space (as it is the case in most tensorial GFT renormalization analyses), the evaluation is often conveniently done in the group representation. It can be carried out analytically, and it proceeds as follows. We compute first the bulk or amputated amplitude ^{1} (if any). The expression so obtained is the irreducible *Master Integral*

*Numerical evaluation in the spin basis* - When the analytic evaluation in group variables is not possible, it is often more convenient to pass to the equivalent expression in terms of group representations (like the ones given above for 4d gravity models), and then proceed numerically, along similar steps as in group variables. First we compute the bulk or amputated amplitude by setting to zero all the spins labeling the external faces and using the appropriate identities for degenerate recoupling coefficients. Next we compute the Master Integral

Let us illustrate the general procedure with an example. We consider the leading order radiative (melonic) correction to the two-point function of this GFT model. The associated Feynman diagram is shown in Figure 1. The diagram has four external and six internal faces, none of which is contractible.

**FIGURE 1**. The picture shows the LO (melonic) radiative correction to the simplicial Ooguri model two-point function

In the holonomic formulation, the amplitude can be written as follows:

In this case, with or without gauge fixing, we can perform all the integrations exactly. Neglecting the contributions from the external faces, without a regularization we would find the divergent result:

The master integral can be regularized either via a sharp cut-off or by heat kernels:

In both cases the amplitude’s degree of divergence reads:

The same result can be recovered by evaluating the amplitude in the spin basis. We have:

Since there are no contractible faces, by setting to zero all the spins labeling the external faces we immediately obtain the expression of the regularized master integral

Equation 3.6 is the result of applying the identity

Thus the degree of divergence can be obtained by evaluating the master integral’s exact formula, or approximately from the above asymptotic formula, by combining the volume factor (replacing the redundant summations) and the face weights with the large-*j* behavior of the Wigner

in agreement with the analytical result obtained in the group formulation.

### 3.2 Radiative Corrections in Simplicial Group Field Theories Models for Quantum Gravity

We now report on some recent results concerning the leading order radiative corrections to *N*-point functions (

We identify the relevant 1PI Feynman diagrams, compute the corresponding master integral formulae and use them to evaluate the master integrals’ scaling (i.e., the diagrams’ superficial degree of divergence) as a function of the cutoff. We also comment on the diagrams’ combinatorial properties and on the structure of the corresponding counterterms. Finally we show how these results can be applied ‘beyond perturbation theory’to characterize to all orders the scaling of the necklace graphs (an important subclass of diagrams appearing in the radiative correction, also identified as the relevant graphs for renormalizability in the tensorial model of ref. 62). We will derive first general formulae that applied to all models in the chosen class (again, simplicial models constructed from constraining those for topological BF theory) and then specialize to the models of interest by choosing the relevant form for the coefficients *w*, encoding the geometrician conditions characterizing them.

##### 3.2.1 Leading Order Corrections to the *N*-Point Functions

The relevant 1PI GFT diagrams, appearing in the perturbative expansion of the 2-point and 4-point functions at the leading order in the GFT coupling constant λ, are shown in Table 1 with the notation explained in the caption.

**TABLE 1**. The 1PI diagrams contributing to the LO expansion of the 2-point and 4-point functions

Let us make a few remarks before moving to the analysis of specific diagrams.

*Selected diagrams* - The diagrams showed in Table 1 are the only potentially divergent 1PI GFT diagrams at the leading order. All the other leading order corrections to the

*point diagrams* - The diagram ^{2}, if any, can be subtracted by mass renormalization (as done in ordinary QFT). The associated irreducible master integral will be denoted as

*The**point diagram* - The diagram *N*-point functions with

To summarize: in order to determine the scaling behavior and divergent structure of the leading corrections to the

##### 3.2.2 The 2-Point Function

To derive the expression of the master integral

where the integer *l* denotes the number of simplicity constraint insertions and the propagator^{3}*w* appearing in the expression *any* simplicial GFT model for constrained BF theory.

Before applying this formula to specific models, let us give a few more details on how it derived, the procedure being in fact the same for the other diagrams.

The melonic 2-point diagram

Hence after setting to zero the spins *type-A* NineJ symbols. Upon using the identity

For completeness, the full formula for the propagator is given by:

The formula for the degenerate propagator (with one vanishing triad) can be written as follows:

Upon setting

The above formulas are completely general. For the euclidean EPRL model they can be further simplified yielding the equations

Let us now focus on the *Duflo model*. In this case the expression (3.9) is still too complicated to be evaluated exactly, even numerically, as it stands. In order to simplify it, we use the asymptotic formula for the

we can rewrite the master integral as follows:

The coefficient ^{4} as a function of the cutoff for different values of the parameters β and *l*.

In the case of the EPRL model, instead, the general formula (3.9) simplifies rather drastically^{5}, yielding the following result

where the (degenerate) propagator

A sample of the results of the numerical evaluation of the relevant master integrals^{6} is in Table 2.

**TABLE 2**. **Top:** numerical evaluation of the expression **Bottom:** numerical evaluation of the EPRL model’s master integral (3.16) as a function of the cutoff

The degree of divergence ^{7} of the scaling exponent

**TABLE 3**. Summary of the estimated scaling exponent

In the case of the Duflo model, some further subtleties arise in the evaluation, due to the more involved nature of the simplicity or geometrician coefficients. These subtleties require additional care in the numerical evaluation of scaling exponents, which are worth emphasizing here, since they are of more general validity in this class of spin foam amplitudes. The master integral formula (3.15) relies on the use of the asymptotic formula for the *j*s are identified, which we cannot do in the Duflo model), the full expression (3.15) is expected to provide *only* a lower bound for

In order to derive an asymptotic formula for the melonic master integral

After appropriate simplifications, the general formula Equation 3.9 becomes.

where *j* regime the propagator (A.14) can be very well approximated by the following expressions:

The EPRL formula can be analytically derived, as shown in Supplementary Appendix A (see A.16, A.17). The corresponding formula for the Duflo model follows from a direct numerical evaluation of the propagator. It is worth noticing that also the Duflo propagator, containing the much more involved Duflo geometrician coefficients Equation 1.6, peak on the same configurations of

The computed values of the scaling exponent α can be found in Table 4.

**TABLE 4**. **Top panel:** numerical evaluation of the propagator **Bottom panel:** computed estimates of the propagator’s scaling exponent α for various values of β and *l*.

The values of ω, obtained by substituting the identities Equation 3.18 into the master integral formula Equation 3.1, are listed in Table 3.

To summarize: the leading order (melonic) correction to the self-energy Equation 3.9 appears to be convergent for the Duflo model and divergent for the EPRL one. The degree of divergence we computed for the EPRL model is in excellent agreement with known analytical results in the literature [41, 42]. Concerning the *Duflo model* for the case

##### 3.2.3 The 4-Point Function

We now focus on the leading (melonic) correction to the 4-point function. The corresponding GFT Feynman diagram

Once more, the above formula is completely general and valid for any simplicial GFT (spin foam) model for constrained BF theory.For the EPRL and Duflo models it specializes to:

where in the first expressions we used the same notations of Equation 3.15.

The degree of divergence of the Master Integral in the Equation 3.20 can be computed by fitting the data shown in Table 5. The resulting values of ω are reported in Table 6.

**TABLE 5**. Numerical evaluation of the expression *l* and β. Here

**TABLE 6**. Numerical values of the divergence’s degree for the 4-point amplitude

The scaling of the EPRL 4-point amplitude can be directly read off from the corresponding formula Equation 3.21.

To summarize: the leading order radiative correction to the 4-point function converges for the Duflo model while it diverges quadratically in the EPRL model. Neglecting possible ambiguities in the definition of both models (which, as we emphasized earlier, could affect the face amplitudes and thus the precise scaling behavior) it would then seem that the Duflo model does not require renormalization, at least at this order, while the EPRL model does. But of course higher orders are needed to establish the renormalizability of both models, thus it is hard to draw too many conclusions from this result. More than the divergence degree in itself, it is important to notice that, since the

Such counterterm is incompatible with a pure simplicial theory space (e.g., a strictly simplicial EPRL model would be non-renormalizable), and this signals the need to extend the theory space of geometric GFT models beyond the simplicial ansatz to include tensorial bubble interactions.

##### 3.2.4. Next-To-Leading Order Corrections to the *N*-point Functions

We now show how to generalize the master integral formulas for the melonic 2-point function

The master integral expressions Equations 3.9, 3.17 rely on the property that each internal link of the diagram is shared exactly by one external face and three internal ones. Hence, after setting to zero the spins labeling the external faces we are left with: i) a pair of

The above combinatorial property is true for any tadpole-free GFT diagram with one and only one external link for each simplicial vertex (here denoted as

with

**FIGURE 3**. Example of GFT diagrams of the type

Upon using the identities Equation 3.18, the asymptotic master integral Equation 3.24 takes the following form:

For the Duflo model the value of α must be computed on case by case basis (see Table 4), while for the EPRL model we have

The degree of divergence ω can also be written in terms of the number of vertices

For the diagrams in Figure 3 the above formula give us the following results, reported in Table 7: the Duflo model amplitudes for all four diagrams are finite; for the EPRL model^{8} the first and last diagrams might be logarithmically divergent and therefore require a deeper analysis [78]. The diagram

**TABLE 7**. Degree of divergence for the diagrams in Figure 3 computed from the power counting formula Equation 3.27.

Once more, this is in fact the crucial lesson we draw from this analysis of divergences, more important, we think, than the precise scaling of the amplitudes, for the reasons already explained.

##### 3.2.5. Beyond Perturbation Theory: The Necklace Diagrams

In the previous two subsections we discussed the all the leading order and a subclass of next-to-leading order radiative corrections to the *N*-point functions with

A (connected) GFT Feynman diagram belongs to the necklace class if (and only if) it consists of an open chain of vertices where each vertex (except the first and last ones) is connected only to its two closest neighbors. Here we restrict to the set of *k*-necklace diagrams with *k* denotes the number of external links, since all other higher order diagrams are either 1-particle reducible or manifestly convergent. Three examples of necklace diagrams are shown in Table 8.

The necklace diagrams share the following remarkable property: the set of internal faces of a *k*-necklace diagram can always be decomposed into the direct sum of two subsets

where the (three) faces in each set *k*-necklace diagram factorizes into *N* identical copies of the master integral

To summarize: the master integral *k*-necklace graph requires an extension of the theory space to include the appropriate tensor invariant interactions of order four and six. A more detailed analysis of the combinatorial structure of the require counterterms is left for future work [78].

## 4 The Road Ahead: Some Suggestions

We close with a set of suggestions for research directions to be pursued, toward a complete understanding the renormalization group flow of simplicial GFT models for 4d quantum gravity.

### 4.1 Scaling and More Scaling, and the Discrete Geometry of Divergent Configurations

The first suggestion is stating the obvious: compute, compute, compute^{9}. We need to (pardon us the word pun) scale up the effort in investigating the scaling of simplicial GFT amplitudes, at the same time trying to import insights from tensorial GFT amplitudes. We need to know much more about the divergent configurations and their dependence on the combinatorics of the underlying cellular complex. Lacking better tools, hard brute-force computations of spin foam amplitudes are the inevitable duty, and these in turn can only build on a better control of the relevant building block the vertex amplitude or, in GFT language, the vertex kernel (which plays an important role also in non-perturbative calculations). Brute force alone will not lead us far, however. On the one hand, we need to develop a more refined analytic understanding of these kernels and resulting amplitudes, and to identify simplified expressions that capture the relevant scaling properties, and their behavior under coarse-graining. On the other hand, where analytical methods do not reach, we need numerical ones to take over; numerical tools for the evaluation of GFT amplitudes are thus badly needed. Packages like have been developed for local QFT Feynmanology would be of course most welcome. On the analytic side, another important objective should be to characterize in detail the (simplicial) geometric meaning of the dominant, most divergent configurations. This is needed to understand the nature of the needed counterterms, but it may also provide insights on the physical features of such GFT models, even beyond their discrete formulation.

The goal here is not to so much to be able to compute GFT Feynman diagrams to arbitrary order (e.g., in vertex or loop expansions). Even in standard QFT, for the physical questions for which the perturbative expansion is the correct approximation scheme, we need to compute (very) many diagrams, but there is often no need nor possibility to go beyond some (usually low) order of approximation and beyond a certain (usually small) number of physical degrees of freedom in the chosen boundary states (of course, the two restrictions go together, since for highly populated boundary states, even the simplest diagrams are of high order). A clear physical picture behind this approximation scheme is as important as computational power. Moreover, coming to the specific quantum gravity case, we would argue that the perturbative GFT expansion, and the description of the dynamics in terms of elementary processes involving few of the fundamental quanta, ie. the usual spin foam language, is not the most convenient approximation to capture the effective continuum physics of quantum gravity, e.g., concerning early cosmology or quantum black holes.

Nor the goal of such analysis of perturbative GFT divergences is establishing that one specific GFT model is finite. Not only renormalizability is a more subtle and possibly interesting feature that finiteness, but one can imagine playing with the ambiguities entering the construction of any given model to modify its scaling behavior and turning it into a finite one. This could be a way to fix or constrain such ambiguities, of course, but it also shows that finiteness per se probably should not be a goal, and that physical conditions fixing the same ambiguities are needed. The main goal of actual computations of GFT amplitudes should then be to provide solid indications on the general power counting of divergences, on the way to a renormalizability proof, and, even more, to indicate the relevant counterterms to be added to the model, and thus the relevant theory space of the starting GFT model. More generally, the goal of such perturbative calculations should be to provide information and tools to be employed to go beyond the perturbative setting and dwell into non-perturbative GFT renormalization. It is only the latter that can provide us with the insights and the results we need to truly explore the formal solidity and effective continuum physics of GFT models of quantum gravity.

### 4.2 Group Field Theories Theory Space, Colors and Relation Between Simplicial and Tensorial Models

We have emphasized several times already the importance of defining the relevant theory space of simplicial GFT models, in order to set up a proper renormalization scheme (perturbative and non-perturbative). Much more work should be devoted to this issue, in particular understanding more about the symmetries of such quantum gravity models. One question is whether the yet to be identified theory space of simplicial GFT models relates to the one of tensorial GFTs. We speculate that they do and, in fact, some hints that the two may largely coincide are known. First, taking seriously the tensorial nature of GFT fields implies coloring (thus distinguishing and ordering) their arguments and, as a consequence, their Feynman diagrams. As we noted, this coloring allows a precise control over the topology of the cellular complexes dual to these Feynman diagrams [7] and, in turn, this greater control allowed for many results that are central in renormalization analyses of tensorial GFTs (e.g., large-N expansions) [10, 11]. A precise control over the topology of the Feynman diagrams, i.e., the cellular complexes on which spin foam amplitudes are based, is needed also in simplicial GFTs, if one aims at identifying the nature of divergences, leading to precise power counting results. It is also needed for identifying key symmetries, as we know already in the case of topological BF models [80], where the complete power counting also relied on the full topological information on the underlying cellular complex [38]. Thus we have a strong argument for relying on colors also in simplicial GFT models of 4d quantum gravity; the form of the corresponding spin foam amplitudes would remain unchanged, but they would now be defined on full 4d cellular complexes, rather than just their 2-skeleton. Assuming we work on colored simplicial GFT models, we then have two preliminary results that suggest a close relation with tensorial models. One is that using colors one can identify similar symmetries in the simplicial case than one finds in the tensorial one [28]. The other is that integrating out all fields except one in a colored simplicial GFT model (in any dimension, with trivial kinetic term) produces an equivalent tensorial GFT model for the remaining field (with the same coupling constant for all interactions) [81]. A third general fact pointing in the same direction is that divergences in simplicial GFT models for topological BF theory, which is the starting point of the construction of simplicial 4d gravity models, are associated to bubbles in the cellular complex, which are in fact the cells associated to allowed interactions in tensorial GFTs with the same base group manifold. These results, in our opinion, suggest that there could be a single theory space containing both (colored) simplicial models and tensorial ones, with interaction kernels in the tensorial directions yet to be identified.

### 4.3 Group Field Theories Models With Local Directions

The third suggestion for further research is to devote attention to the renormalization of GFT models which combine the combinatorially non-local pairing structure on geometric variables, in GFT interactions, with the presence of local directions. This includes both simplicial GFT models and tensorial ones, with the distinction referring to the pairing of geometric variables.

There are two main examples of such ‘mixed’models. One is the tensorial models used to describe SYK-like many-body systems [82], whose renormalization has been in fact studied in several cases. Here the non-local, tensorial indices are usually reduced to finite sets (we have thus simple tensor models, rather than full GFTs) and the single local direction is a time variable. The standard SYK models are indeed quantum mechanical models in 0 + 1 dimensions, with generalizations to higher dimensions (thus, with more local directions) having been proposed. The other class of mixed models is the extension of (simplicial) GFT quantum gravity models to include scalar fields coupled to gravitational degrees of freedom [83]. These extended models have been studied in particular in the context of GFT condensate cosmology [75–77, 84, 85], with the additional scalar fields playing (also) the role of clock and rods that allow to define relational, diffeo-invariant observables in terms of which an effective cosmological dynamics can be extracted from the GFT hydrodynamics.

The potential physical interest of these models, and of their renormalization analysis, is thus obvious. They present several interesting issues. The presence of both local and non-local directions may modify sensibly the renormalization flow and the structure of divergences, thus leading to different dominant diagrams and effective dynamics in both UV and IR sectors. One can also envisage setting up an altogether different renormalization group scheme, adopting a notion of scale tied to the scalar (local) directions, rather than the group manifold (or involving both), potentially producing very different results. Such focus on the flow parametrized by variables with a (tentative) physical interpretation as relational time/space variables may also allow a more direct physical interpretation of the renormalization flow itself, e.g., in a cosmological context (even though similar cautionary remarks as for the usual renormalization scheme would apply here).

### 4.4 Relation With Lattice Spin Foam Renormalization

We have emphasized how renormalizing a GFT model is tantamount to renormalizing (and studying the continuum limit of) the corresponding discrete gravity path integral and spin foam amplitudes, from a different standpoint. But the GFT formalism is only one way to provide a complete definition of spin foam models, the other being to view them as a peculiar (because background independent) lattice theory and setting up some appropriate refinement procedure. Therefore, it would be very important to compare results obtained in the context of GFT renormalization, especially for simplicial quantum gravity models, with the results and techniques developed for renormalizing spin foam amplitudes from a lattice gauge theory perspective [65, 66, 86–88].

In this lattice-focused approach to spin foam renormalization, a cut-off is also imposed on representation variables, but the notion of ‘scale’is rather given by the combinatorial complexity of the underlying lattice, and the renormalization group flow is driven by refinement/coarse-graining steps ordered by such complexity. Refinement/coarse-graining steps affect both bulk lattices and boundary graphs, and the flow of quantum amplitudes is constrained by the requirement of their consistency under restriction to coarser boundary states.

Despite their differences, the two renormalization schemes share several, since also GFT subtraction moves amount to lattice coarse-graining steps, and corresponding maps between associated amplitudes are also built-in in the (perturbative) QFT renormalization steps used in the GFT context. Still, a detailed work of translation between the two frameworks would be very useful. This work may require, on the GFT side, a combination of functional renormalization group techniques, since we are interested in the continuum limit of spin foam models, and perturbative expansions, given that spin foam models arise in such expansion. This comparison would be beneficial for both approaches; in particular, it would emphasize the role of combinatorial complexity of boundary states in the GFT renormalization flow. This work should be carried out for all models that have been studied in both settings (also in the lattice renormalization approach work has been confined mostly to highly simplified models), aiming of course at unraveling the continuum phase diagram of 4d quantum gravity from two perspectives at once.

### 4.5 Group Field Theories Renormalization Via Tensor Networks

One powerful set of techniques coming from the theory of quantum many-body systems, that have been already applied in the context of lattice-based renormalization of spin foam models, uses the language of tensor networks [89, 90]. This language is useful both for numerical studies and for emphasizing the role of entanglement in the renormalization group flow [91, 92]; in particular, it allows to unravel topological quantum phases of many-body systems.

In the case of GFT models, the interest in importing techniques from tensor networks goes beyond these general facts, and stems also from the fact that GFT states themselves can be seen as generalized tensor networks [93], and by the related fact that entanglement is responsible for the basic connectivity between GFT quanta that gives rise to extended discrete structures labeled by quantum geometric data. The many facets of the GFT formalism, moreover, would allow for a manifold application of tensor network techniques. On the one hand the basic GFT field is a tensor and its quantum states are tensor networks, as mentioned; on the other hand, it remains a QFT, calling for continuum tensor network techniques as employed, say, in standard scalar quantum field theory [94]. At the same time, its Feynman amplitudes are lattice gauge theories, to which a different set of tensor network techniques can be applied [95] (as developed in the context of spin foam lattice renormalization). And they remain quantum many-body systems, peculiar for their background independent nature, but still conventional enough to allow the deployment of tensor network methods taken from their natural context.

## Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

## Author Contributions

All authors have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

MF is supported by the International Max Planck Research School for Mathematical and Physical Aspects of Gravitation, Cosmology and Quantum Field Theory. DO gratefully acknowledges funding from Deutsche Forschungs Gemeinschaft.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2020.552354/full#supplementary-material.

^{1}An internal face of

^{2}The first melonic correction to the self energy might also have a subleading divergence as in the case of the Ooguri model and (Lorentzian) EPRL model. Such divergence is responsible for the wave function renormalization. This in turns requires a modification of the covariance with a second-order derivative term in order to account for the new counterterm.

^{3}According to the diagram’s connectivity, each internal edge of the graph

^{4}The main limitation on range of cutoff values we can test depends on the computational resources available. Here we choose

^{5}The *l*, since the simplicity constraints act as a projector.

^{6}They will be discussed in detail, alongside the results for other values of the various parameters, in a follow-up publication [78].

^{7}To derive an estimate for ω we took the statistical average of the values obtained by fitting the data points in the cutoff ranges

^{8}We point out again that we are studying both models under specific choices fixing the various ambiguities that enter the construction of the spin foam amplitudes. These ambiguities affect, in general, the scaling results.

^{9}To be clear: 1-loop and 2-loop calculations would be probably enough to extract a lot of interesting properties from the perturbative expansion of GFT models for 4d quantum gravity. They may even be sufficient, since we actually expect that most interesting physics should be looked for elsewhere, i.e., within different regimes and approximations of the fundamental quantum dynamics, and not in the perturbative (spin foam) expansion. Moreover, as we stressed, such perturbative calculations should have as main goal to identify the theory space within which the relevant GFT models should be placed, more than simply pushing the perturbative analysis for its own sake. We emphasize, however, that, in fact, a lot of computational effort is required to perform such 1-loop and 2-loop calculations because one needs to consider many and complicated GFT diagrams, and involved spin foam amplitudes, already at this order. These are the calculations we intend to encourage.

## References

1. Oriti D. *The microscopic dynamics of quantum space as a group field theory*. In: Proceedings, foundations of space and time: reflections on quantum gravity. 2009 Aug 10–14; Cape Town, South Africa. Cambridge, United Kingdom: Cambridge University Press (2011). Vol. 257.

4. Chirco G, Kotecha I, Oriti D. Statistical equilibrium of tetrahedra from maximum entropy principle. *Phys Rev* (2019) D99:086011. doi:10.1103/physrevd.99.086011

5. Oriti D. Group field theory as the 2nd quantization of loop quantum gravity, class. *Quant. Grav* (2016) 33:085005. doi:10.1088/0264-9381/33/8/085005

6. Freidel L. Group field theory: an Overview. *Int J Theor Phys* (2005) 44:1769–1783. doi:10.1007/s10773-005-8894-1

8. Rivasseau V. The tensor track, III. *Fortschritte der Physik* (2014) 62:81–107. doi:10.1002/prop.201300032

10. Carrozza S. *Tensorial methods and renormalization in group field theories*. [PhD thesis]. Paris: Orsay, LPT (2013).

11. Carrozza S. Flowing in group field theory space: a review. *Sigma* 12 (2016). 070. doi:10.3842/sigma.2016.070

12. Pereira R. *Spinfoams from simplicial geometry*. [PhD thesis]. Aix-Marseille II: U. Mediterranee (2010).

13. Dupuis M, Livine ER. Lifting SU(2) spin networks to projected spin networks. *Phys Rev D* (2010) 82:064044. doi:10.1103/physrevd.82.064044

16. Baratin A, Oriti D. Group field theory and simplicial gravity path integrals: a model for Holst-Plebanski gravity. *Phys. Rev* (2012) D85:044003. doi:10.1103/physrevd.85.044003

17. Finocchiaro M, Oriti D. Spin foam models and the Duflo map. *Classical Quant. Grav* (2020) 37:015010. doi:10.1088/1361-6382/ab58da

18. Perez A. The spin-foam approach to quantum gravity. *Living Rev. Relat* (2013) 16:3. doi:10.12942/lrr-2013-3

19. Dupuis M, Ryan JP, Speziale S. Discrete gravity models and Loop Quantum Gravity: a short review. *Sigma* (2012) 8:052. doi:10.3842/sigma.2012.052

20. Baez JC, Barrett JW. The Quantum tetrahedron in three-dimensions and four-dimensions. *Adv. Theor. Math. Phys*. (1999) 3:815

21. Freidel L, Speziale S. Twisted geometries: a geometric parametrisation of SU(2) phase space. *Phys. Rev* (2010) D82:084040. doi:10.1103/physrevd.82.084040

22. Freidel L, Livine ER. Bubble networks: framed discrete geometry for quantum gravity. *Gen. Relat. Gravit*. (2019) 51:9. doi:10.1007/s10714-018-2493-y

23. Freidel L, Livine ER, Pranzetti D. Gravitational edge modes: from Kac–Moody charges to Poincaré networks. *Classical Quant. Grav*. (2019) 36:195014. doi:10.1088/1361-6382/ab40fe

24. Reisenberger MP, Rovelli C. Space-time as a feynman diagram: the connection formulation, class. *Quant. Grav*. (2001) 18:121–140. doi:10.1088/0264-9381/18/1/308

25. Dittrich B, Speziale S. Area-angle variables for general relativity. *New J. Phys*. (2008) 10:083006. doi:10.1088/1367-2630/10/8/083006

26. Baloitcha E, Lahoche V, Ousmane Samary D. Flowing in discrete gravity models and Ward identities: a review. arXiv:2001.02631v2

27. Kegeles A, Oriti D. Generalized conservation laws in non-local field theories. *J. Phys*. (2016) A49:135401. doi:10.1088/1751-8113/49/13/135401

28. Kegeles A, Oriti D. Continuous point symmetries in group field theories. *J. Phys*. (2017) A50:125402. doi:10.1088/1751-8121/aa5c14

29. Freidel L, Louapre D. Diffeomorphisms and spin foam models. *Nucl. Phys*. (2003) B662:279–298. doi:10.1016/s0550-3213(03)00306-7

30. Dittrich B. Diffeomorphic symmetry in quantum gravity models. *Adv. Sci. Lett*. (2008) 2(10):151. doi:10.1166/asl.2009.1022

31. Bahr B, Dittrich B. (Broken) gauge symmetries and constraints in Regge calculus. *Classical Quant. Grav*. (2009) 26:225011. doi:10.1088/0264-9381/26/22/225011

32. Bahr B, Dittrich B. Improved and perfect actions in discrete gravity. *Phys. Rev. D* (2009) 80:124030. doi:10.1103/physrevd.80.124030

33. Baratin A, Girelli F, Oriti D. Diffeomorphisms in group field theories. *Phys Rev. D* (2011) 83:104051. doi:10.1103/physrevd.83.104051

34. Dittrich B, Kamiński W, Steinhaus S. Discretization independence implies non-locality in 4D discrete quantum gravity. *Classical Quant. Grav*. (2014) 31:245009. doi:10.1088/0264-9381/31/24/245009

35. De Pietri R, Freidel L, Krasnov K, Rovelli C. Barrett-Crane model from a Boulatov-Ooguri field theory over a homogeneous space. *Nucl. Phys*. (2000) B574:785–806. doi:10.1016/s0550-3213(00)00005-5

36. Perez A, Rovelli C. A Spin foam model without bubble divergences. *Nucl. Phys* (2001) B599:255–282. doi:10.1016/s0550-3213(01)00030-x

37. Perez A. On the regularization ambiguities in loop quantum gravity. *Phys. Rev*. (2006) D73:044007. doi:10.1103/PhysRevD.73.044007

38. Bonzom V, Smerlak M. Bubble divergences: sorting out topology from cell structure. *Ann Henri Poincaré* (2012) 13:185. doi:10.1007/s00023-011-0127-y

39. Ben Geloun J, Bonzom V. Radiative corrections in the Boulatov-Ooguri tensor model: the 2-point function. *Int. J. Theor. Phys*. (2011) 50:2819–2841. doi:10.1007/s10773-011-0782-2

40. Ben Geloun J. On the finite amplitudes for open graphs in Abelian dynamical colored Boulatov–Ooguri models. *J. Phys*. (2013) A46:402002. doi:10.1088/1751-8113/46/40/402002

41. Perini C, Rovelli C, Speziale S. Self-energy and vertex radiative corrections in LQG. *Phys. Lett*. (2009) B682:78. doi:10.1016/j.physletb.2009.10.076

42. Ben Geloun J, Gurau R, Rivasseau V. EPRL/FK group field theory. *Europhys. Lett*. (2010) 92:60008. doi:10.1209/0295-5075/92/60008

43. Bonzom V, Dittrich B. Bubble divergences and gauge symmetries in spin foams. *Phys. Rev* (2013) D88:124021. doi:10.1103/physrevd.88.124021

44. Riello A. Self-energy of the lorentzian engle-pereira-rovelli-livine and freidel-krasnov model of quantum gravity. *Phys. Rev*. (2013) D88:024011. doi:10.1103/physrevd.88.024011

45. Chen LQ. Bulk amplitude and degree of divergence in 4d spin foams. *Phys. Rev*. (2016) D94:104025. doi:10.1103/physrevd.94.104025

46. Donà P. Infrared divergences in the EPRL-FK spin foam model. *Classical Quant. Grav*. (2018) 35:175019. doi:10.1088/1361-6382/aad38f

47. Dona P, Sarno G. Numerical methods for EPRL spin foam transition amplitudes and Lorentzian recoupling theory. *Gen. Relat. Gravit*. (2018) 50:127. doi:10.1007/s10714-018-2452-7

48. Donà P, Gozzini F, Sarno G. Searching for classical geometries in spin foam amplitudes: a numerical method. arXiv:1909.07832

49. Donà P, Fanizza M, Sarno G, Speziale S. Numerical study of the Lorentzian Engle-Pereira-Rovelli-Livine spin foam amplitude. *Phys. Rev*. (2019) D100:106003. doi:10.1103/physrevd.100.106003

50. Haggard HM, Littlejohn RG. Asymptotic of the wigner 9j symbol. *Classical Quant. Grav*. (2010) 27:135010. doi:10.1088/0264-9381/27/13/135010

51. Yu L, Littlejohn RG. Semiclassical analysis of the Wigner 9j symbol with small and large angular momenta. arXiv:1104.1499

53. Bonzom V, Fleury P. Asymptotic of wigner 3nj-symbols with small and large angular momenta: an elementary method. *J. Phys* (2012) A45:075202. doi:10.1088/1751-8113/45/7/075202

54. Barrett JW, Fairbairn WJ, Hellmann F. Quantum gravity asymptotic from the SU(2) 15j symbol. *Int. J. Mod. Phys*. (2010) A25:2897–2916. doi:10.1142/s0217751x10049281

55. Barrett JW, Dowdall R, Fairbairn WJ, Gomes H, Hellmann F, Pereira R. Asymptotic of 4d spin foam models. *Gen. Relat. Gravit*. (2011) 43:2421. doi:10.1007/s10714-010-0983-7

56. Barrett JW, Dowdall R, Fairbairn WJ, Gomes H, Hellmann F. Asymptotic analysis of the EPRL four-simplex amplitude. *J. Math. Phys*. (2009) 50:112504. doi:10.1063/1.3284398

57. Johansson H, Forssén C. Fast and accurate evaluation of Wigner 3j, 6j, and 9j symbols using prime factorisation and multi-word integer arithmetic. *SIAM J. Sci. Stat. Comput*. (2016) 38:A376. doi:10.1137/15M1021908

58. Speziale S. Boosting Wigner’s nj-symbols. *J. Math. Phys*. (2017) 58:032501. doi:10.1063/1.4977752

59. Donà P, Fanizza M, Sarno G, Speziale S. SU(2) graph invariants, Regge actions and polytopes. *Classical Quant. Grav*. (2018) 35:045011. doi:10.1088/1361-6382/aaa53a

60. Dona P, Speziale S. Asymptotic of lowest unitary SL(2,C) invariants on graphs. *Phys. Rev*. (2020) D102:086016. doi:10.1103/physrevd.102.086016

61. Dona P, Fanizza M, Martin-Dussaud P, Speziale S. Asymptotic of SL(2,C) coherent invariant tensors. arXiv:2011.13909

62. Carrozza S, Lahoche V, Oriti D. Renormalizable Group Field Theory beyond melonic diagrams: an example in rank four. *Phys. Rev*. (2017) D96:066007. doi:10.1103/physrevd.96.066007

63. Oriti D. Group field theory as the microscopic description of the quantum spacetime fluid: a New perspective on the continuum in quantum gravity. PoS QG-PH. [arXiv:0710.3276]

65. Dittrich B, The continuum limit of loop quantum gravity–a framework for solving the theory In: A Ashtekar, and J Pullin, editors *Loop quantum gravity: the first 30 years* (2017). p 153.

66. Delcamp C, Dittrich B. Towards a phase diagram for spin foams. *Classical Quant. Grav*. (2017) 34:225006. doi:10.1088/1361-6382/aa8f24

67. Eichhorn A, Koslowski T, Pereira AD. Status of background-independent coarse-graining in tensor models for quantum gravity. *Universe* (2019) 5:53. doi:10.3390/universe5020053

68. Eichhorn A, Lumma J, Pereira AD, Sikandar A. Universal critical behavior in tensor models for four-dimensional quantum gravity. *J. High Energy Phys*. (2020) 02:110. doi:10.1007/jhep02(2020)110

69. Freidel L, Louapre D. Nonperturbative summation over 3-D discrete topologies. *Phys. Rev*. (2003) D68:104004. doi:10.1103/physrevd.68.104004

70. Magnen J, Noui K, Rivasseau V, Smerlak M. Scaling behaviour of three-dimensional group field theory, Class. *Quant. Grav*. (2009) 26:185012. doi:10.1088/0264-9381/26/18/185012

71. Baratin A, Carrozza S, Oriti D, Ryan J, Smerlak M. Melonic phase transition in group field theory. *Lett. Math. Phys*. (2014) 104:1003–1017. doi:10.1007/s11005-014-0699-9

72. Carrozza S, Lahoche V. Asymptotic safety in three-dimensional SU(2) Group Field Theory: evidence in the local potential approximation. *Classical Quant. Grav*. (2017) 34:115004. doi:10.1088/1361-6382/aa6d90

73. Ben Geloun J, Koslowski TA, Oriti D, Pereira AD. Functional Renormalization Group analysis of rank 3 tensorial group field theory: the full quartic invariant truncation. *Phys. Rev*. (2018) D97:126018. doi:10.1103/physrevd.97.126018

74. Ben Geloun J, Martini R, Oriti D. Functional renormalisation group analysis of tensorial group field theories onRd. *Phys. Rev*. (2016) D94:024017. doi:10.1103/physrevd.94.024017

75. Oriti D, Sindoni L, Wilson-Ewing E. Emergent Friedmann dynamics with a quantum bounce from quantum gravity condensates. *Classical Quant. Grav*. (2016) 33:224001. doi:10.1088/0264-9381/33/22/224001

76. Gielen S, Sindoni L. Quantum cosmology from group field theory condensates: a review. *Sigma* (2016) 12:082. doi:10.3842/sigma.2016.082

77. Oriti D. The universe as a quantum gravity condensate. *Compt. Rendus Phys* (2017) 18:235–245. doi:10.1016/j.crhy.2017.02.003

79. Celoria M, Finocchiaro M, Oriti D. *Numerical analysis of the fusion coefficients of the Duflo spin foam model for 4d Riemannian quantum gravity. To appear*.

80. Bonzom V, Smerlak M. Gauge symmetries in spinfoam gravity: the case for ‘cellular quantization’. *Phys. Rev. Lett*. (2012) 108:241303. doi:10.1103/physrevlett.108.241303

81. Gurau R. A generalization of the Virasoro algebra to arbitrary dimensions. *Nucl. Phys*. (2011) B852:592–614. doi:10.1016/j.nuclphysb.2011.07.009

82. Delporte N, Rivasseau V. The tensor track V: holographic tensors. In: Proceedings, 17th hellenic school and workshops on elementary particle physics and gravity (CORFU2017); 2017 Sep 2–28; Corfu, Greece. Trieste, Italy: SISSA (2017).

83. Li Y, Oriti D, Zhang M. Group field theory for quantum gravity minimally coupled to a scalar field. *Classical Quant. Grav*. (2017) 34:195001. doi:10.1088/1361-6382/aa85d2

84. Wilson-Ewing E. A relational Hamiltonian for group field theory. *Phys. Rev*. (2019) D99:086017. doi:10.1103/physrevd.99.086017

85. Gielen S, Oriti D. Cosmological perturbations from full quantum gravity. *Phys. Rev. D*. (2018) 98:106019. doi:10.1103/PhysRevD.98.106019

86. Dittrich B, Schnetter E, Seth CJ, Steinhaus S. Coarse graining flow of spin foam intertwiners. *Phys. Rev* (2016) D94:124050. doi:10.1103/physrevd.94.124050

87. Bahr B, Steinhaus S. Hypercuboidal renormalization in spin foam quantum gravity. *Phys. Rev*. (2017) D95:126006. doi:10.1103/physrevd.95.126006

88. Bahr B. On background-independent renormalization of spin foam models. *Classical Quant. Grav*. (2017) 34:075001. doi:10.1088/1361-6382/aa5e13

89. Dittrich B, Mizera S, Steinhaus S. Decorated tensor network renormalization for lattice gauge theories and spin foam models. *New J. Phys*. (2016) 18:053009. doi:10.1088/1367-2630/18/5/053009

90. Cunningham WJ, Dittrich B, Steinhaus S. Tensor network renormalization with fusion charges: applications to 3d lattice gauge theory. arXiv:2002.10472

92. Orus R. A practical introduction to tensor networks: matrix product states and projected entangled pair states. *Ann. Phys*. (2014) 349:117–158. doi:10.1016/j.aop.2014.06.013

93. Chico G., Oriti D., Zhang M., Group field theory and tensor networks: toward a Ryu–Takayanagi formula in full quantum gravity, *Class. Quant. Grav*.35 (2018) 115011

94. Hu Q, Franco-Rubio A, Vidal G. Continuous tensor network renormalization for quantum fields. arXiv:1809 05176

Keywords: quantum gravity, group field theory, spin foam models, renormalization, lattice gravity

Citation: Finocchiaro M and Oriti D (2021) Renormalization of Group Field Theories for Quantum Gravity: New Computations and Some Suggestions. *Front. Phys.* **8**:552354. doi: 10.3389/fphy.2020.552354

Received: 15 April 2020; Accepted: 10 December 2020;

Published: 25 February 2021.

Edited by:

Benjamin Bahr, University of Erlangen Nuremberg, GermanyCopyright © 2021 Finocchiaro and Oriti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marco Finocchiaro, marcfinocchiaro@gmail.com; Daniele Oriti, daniele.oriti@physik.lmu.de