**1. Introduction**

Two potentially revolutionary ideas have inspired much work in contemporary theoretical physics. Both ideas herald from the use of the general information theoretic approach to the problem of quantum gravity. The first idea is that—*the world is holographic*—with the physics of several semi-classical systems (by which we mean systems in which matter is treated quantum-mechanically while spacetime and geometry are treated classically) entirely captured on spacetime regions of one dimension lower.

The evidence for holography is already suggestive when considering classical gravitational systems like black holes, or more general causal horizons, and the semi-classical physics of quantum fields in their vicinity. This evidence is strengthened by numerous results in the context of the AdS/CFT correspondence [1,2], where, in particular, holography as found in semi-classical gravitational systems is put in correspondence with general properties of non-gravitational, purely quantum mechanic *dual* many-body systems. Indeed, holographic features deeply characterise condensed matter physics—hence the suggestion that there may be a *purely quantum mechanic origin of holography* that may in fact underlie classical, gravitational, geometric physics as studied in the general relativistic context.

The second idea is that *geometry itself originates from entanglement*. The recent quantum information-theoretic paradigm for gravity has provided a new vision of the cosmos wherein the universe, together with its topology, its geometry and its macroscopic dynamics, arise from the entanglement between the fundamental constituents of some exotic underlying quantum system.

In this direction, in particular, along with the increasing impact of condensed matter physics in string theory and gauge/gravity duality, people have started exploring the use of tensor network (TN) algorithms [3–12] from condensed matter theory in quantum gravity [13–20], providing interesting insights on holographic duality and its generalisation in terms of geometry/entanglement correspondence [21,22]. Approaches like the AdS/MERA [10,23], where the geometry of the auxiliary tensor network decomposition of the quantum many-body vacuum state is interpreted as a representation of the dual spatial geometry, are providing an intriguing constructive framework for investigating holography beyond AdS/CFT [13–16,24].

However, if everything is quantum at its root, and this is true not only for ordinary systems living in spacetime but for *space, time and geometry themselves*, then this implies that the very holographic behaviour of the universe is the result of purely quantum properties of the microscopic constituents of spacetime. In addition, this line of research, therefore, indicates a strong need for a quantum foundation of holography.

In fact, background independence naturally leads to a quantum description of the universe in terms of fundamental quantum many-body physics of discrete and purely algebraic microscopic constituents [25–34], from which spacetime emerges only at an effective, approximate level, out of a texture of quantum correlations [35–43]. This means that the two suggestions that holography has a purely quantum origin and that geometry itself comes from entanglement are extremely natural when seen from the perspective of quantum gravity formalisms, in which spacetime and geometry are ultimately emergen<sup>t</sup> notions. However, more than that, differently from the semiclassical framework, in the non-perturbative scheme, holography as detected in gravitational systems, as well as any macroscopic feature of our geometric universe, not only would result from purely quantum properties of the microscopic constituents of spacetime, but they can— *only*—be understood in this light.

This perspective is manifest in the *Group Field Theory* (GFT) formalism, a promising convergence of the insights and results from matrix models [44,45], loop quantum gravity and simplicial approaches into a background independent quantum field theory setup. The GFT approach to quantum gravity [30, 46–49] provides a very general quantum many-body formulation of the spacetime micro-structure, for instance of the spin networks and discrete quantum geometry states of Loop Quantum Gravity [25–27, 29], with a Fock space description where the quantum GFT fields create and annihilate elementary building blocks of space, interpreted as (*d* − 1)-simplices in *d* spacetime dimensions, organised in nontrivial combinatorial tensor network structures.

As a higher order generalisation of matrix models, the GFT formalism at the same time provides a field-theoretic and inherently covariant framework for generalising the *tensor networks* approach to the holographic aspects of quantum many-body systems in condensed matter and in the AdS/CFT context. This makes GFT a very effective framework to investigate how space-time geometry, together with its holographic behaviour and macroscopic dynamics, arise from entanglement between the fundamental constituents.

In this paper, we review a series of recent results [18–20] concerning the definition of entanglement entropy in the GFT framework and the characterisation of its holographic behaviour. In particular, we focus on the definition of the notion of entanglement entropy in the second quantised formalism of GFT setting and on the set of choices which eventually lead to a holographic behaviour for the entanglement entropy.

The manuscript is organised as follows: Section 2 shortly reviews the framework of group field theory while focussing the attention of the reader on those aspects of the GFT fields that play a major role in the forthcoming derivation. Section 3 introduces the second quantisation formalism for group field theory, defines the notion of multi-particle state observable for quantum geometry in the GFT Fock space and specifies a class of GFT coherent state basis, necessary for the definition of the entanglement entropy expectation value. In Section 4, the statistical derivation of the Rényi entropy for a bipartite GFT open spin network state is reviewed, with an emphasis on the Fock space setting. The formal mapping of the expectation values of the Rényi entropy to BF theory amplitudes is described and the resulting divergence degree and scaling of the entropy analysed. Section 5 briefly comments on the results on the holographic scaling for the interacting GFT case. A brief discussion and an appendix on the notion of coherent states over-completeness close the manuscript.

### **2. Group Field Theory**

Group field theories (GFT) are quantum field theories defined on *d* copies of a compact Lie group *G* with combinatorially non-local interactions. The dynamics of the GFT field

$$
\phi: G^{\times d} \to \mathbb{C}
$$

are specified by a probability measure

$$\mathrm{d}\mu\_{\mathbb{C}}(\phi,\overline{\phi})\exp\left(-S\_{\mathrm{int}}[\phi,\overline{\phi}]\right),\tag{1}$$

comprised of a *Gaussian measure* d*μC*, associated with a positive covariance kernel operator *C* defining the propagator of the theory [45],<sup>1</sup> and a *perturbation* around it, given by an interaction term *<sup>S</sup>*int[*φ*, *φ*], generically parametrized as

$$S\_{\mathrm{int}}[\boldsymbol{\phi}, \overline{\boldsymbol{\phi}}] = \sum\_{\substack{\mathcal{T} \\ p, q \ge 0}} \sum\_{\substack{\mathfrak{p} + q = \mathcal{T} \\ p, q \ge 0}} \int \prod\_{\mathfrak{p} = 1}^{p} \mathrm{d} \mathbf{g}\_{\mathfrak{p}}^{\prime} \, \overline{\boldsymbol{\phi}}(\mathbf{g}\_{\mathfrak{p}}^{\prime}) \prod\_{q = 1}^{q} \mathrm{d} \mathbf{g}\_{\mathfrak{q}} \, \boldsymbol{\phi}(\mathbf{g}\_{\mathfrak{q}}) \, \lambda\_{\mathcal{T}} V\_{\mathcal{T}}(\mathbf{g}\_{1 \prime}^{\prime} \cdot \cdot \cdot \cdot , \mathbf{g}\_{\mathfrak{p}}^{\prime}; \mathbf{g}\_{1 \prime} \cdot \cdot \cdot \cdot , \mathbf{g}\_{\mathfrak{q}}) \, \lambda\_{\mathcal{T}} V\_{\mathcal{T}}(\mathbf{g}\_{\mathfrak{q}}^{\prime}) \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \, \lambda\_{\mathcal{T}} \,$$

where I denotes a term in the set of elementary interactions, *V*I is the specific monomial in the fields associated with the interaction, and *λ*I is the respective coupling constant. Together with the field valence *d* and symmetry, the specific choice of the covariance and interaction kernels identifies the GFT model.

We are particularly concerned with three peculiar aspects of the GFT formalism, which will combine at the hearth of the following derivation. The first is that the dynamical *d*-valent GFT field *φ* combinatorially behaves as an infinite dimensional rank-*d* tensor, with indices labelled by elements of the compact Lie group *G* [18]. This is apparent in the combinatorially non-local structure of the interaction kernels *V*I , as functions on *Gd*×|I| for finite sets of interactions. The kernels *V*I do not impose coincidence of points in the group space *<sup>G</sup>*<sup>×</sup>*d*, but the whole set of the *d* × |I| field arguments is partitioned into pairs, convoluted "strandwise" by the kernel,

$$V\_{\mathcal{T}}(\{\mathbf{g}\}\_{\mathcal{T}}) = V(\{\mathbf{g}\_{\mathfrak{p}}\}\_{\mathfrak{l}}\{\mathbf{g}\_{\mathfrak{q}}^{-1}\}\_{\mathfrak{l}}).\tag{2}$$

Indeed, one can see GFTs as higher-rank, infinite dimensional generalizations of random matrix models [45]. For instance, if we take *G* as Z*N*, then group fields identically reduce to rank-*d* tensors,<sup>2</sup> where integrability with respect to the discrete Dirac measure *μ* is satisfied for all fields considered.

The second aspect is that specific GFT partition functions, where the group *G* is identified with the local gauge subgroup of gravity and the kernels properly chosen, define generating functions

$$
\overline{\varphi}\_1 \cdot \varphi\_2 = \int \mathbf{dg} \overline{\varphi}\_1(\mathbf{g}) \,\varphi\_2(\mathbf{g}) \,\varphi\_1
$$

for any two square-integrable functions *ϕ*1 and *ϕ*2 on *G<sup>d</sup>*.)

$$\mathbf{g} = (\mathfrak{g}\_1, \dots, \mathfrak{g}\_d) \in G^d, \qquad \mathbf{d}\mathbf{g} = \mathbf{d}\mathfrak{g}\_1 \cdots \mathbf{d}\mathfrak{g}\_d.$$

2 A rank-*d* tensor *T* with index cardinality *N* is a complex field on *d* copies of the cyclic group Z*N*:

$$\mathcal{T}: \mathcal{Z}\_N^{\times d} \to \mathbb{C}\_r$$

which defines a state in H*d*,*<sup>N</sup>* the space of

tensors with fixed rank *d* and index cardinality *N*. Neglecting the structure of the cyclic group, H*d*,*<sup>N</sup>* is reduced to C*<sup>N</sup><sup>d</sup>* .

<sup>1</sup> We use a vector notation for the configuration space variables and its Haar measure (We will also use the short-hand notation: 

for the covariant quantization of Loop Quantum Gravity (LQG) in terms of spin foam models (for instance, [28]). In particular, LQG spin network states, describing three-dimensional discrete quantum geometries at the boundary of the spin foam transition amplitudes, can be expressed as expectation values of specific GFT operators. In a formalism of second quantisation for the GFT, spin network boundary states then become elements of the GFT Fock space and spin network vertices, intended as atoms of space that can be put in direct correspondence with fundamental GFT quanta that are created or annihilated by the field operators of GFT.

When *G* is set to correspond to the local gauge group of GR, e.g., the Lorentz group SO(1, *d* − 1) or its universal covering, or SU(2) in dimension *d* = 3, the gauge symmetry leads to an invariance of the GFT action under the (right) diagonal action of *G*. GFT fields are constrained to satisfy a *gauge invariance condition*, defined as a global symmetry of the GFT field under simultaneous translation of its group variables:

$$\forall h \in G, \qquad \phi(g\_1h, \dots, g\_dh) = \phi(g\_1, \dots, g\_d). \tag{3}$$

The gauge invariance condition, or *closure constraint*, is the main dynamical ingredient of GFT models for quantum BF theory in arbitrary dimension. In *d* = 3, SU(2) BF theory can be interpreted as a theory of Euclidean gravity, and therefore SU(2) GFT with closure constraint provides a natural arena in which to formulate 3D Euclidean quantum gravity models. A typical example is the Boulatov model [50] (which generates Ponzano–Regge spin foam amplitudes [51]), which will constitute an important ingredient of the following derivation.

The third aspect is that, for such (simplicial) GFT models, endowed with a geometric interpretation of the dynamical fields, the very field-theoretic nature of the GFT formalism provides a powerful tool to describe quantum geometry states as a peculiar *quantum many-body* systems in the formalism of *second quantisation*. For large systems in quantum mechanics, we know that the concept of a particle fades away and is replaced by the notion of an excitation of a given mode of the field representing the particle. Similarly, in the GFT description, we expect the solid graph description of spin networks quantum geometry to fade into a dynamical net of excitations of the GFT field over a vacuum. Given the tensorial behaviour of the GFT field, such a quantum many-body description turns quantum geometry states into collective, purely combinatorial and algebraic analogues of quantum *tensor networks states*.

These three ingredients together make the GFT second quantised formalism an specially convenient setting to *quantitatively* investigate the relation between geometry and entanglement in quantum gravity, taking advantage of the most recent techniques and tools of quantum statistical mechanics, information theory and condensed matter theory. In particular, quantum tensor network algorithms provide a constructive tool to investigate the roots of the holographic behaviour of gravity at the quantum level.

### **3. The GFT Fock Space**

In a second quantisation scheme [30], multi-particle states of the quantum GFT field *φ*(**g**) can be organised in a Fock space F generated by a Fock vacuum |0 and field operators

$$
\hat{\boldsymbol{\varphi}}(\mathbf{g}) \equiv \hat{\boldsymbol{\varphi}}(\mathbf{g}\_1, \dots, \mathbf{g}\_d), \qquad \hat{\boldsymbol{\varphi}}^\dagger(\mathbf{g}) \equiv \hat{\boldsymbol{\varphi}}^\dagger(\mathbf{g}\_1, \dots, \mathbf{g}\_d), \tag{4}
$$

assumed to be invariant under the diagonal action of the group *ϕ*<sup>+</sup>(**g***h*) ≡ *ϕ*<sup>+</sup>(*g*1*h*, ··· , *gVh*) = *ϕ*<sup>+</sup>(**g**), for *h* ∈ *G*, consistently with (3), and to obey canonical commutation relations (bosonic statistics)

$$\begin{aligned} \left[\hat{\boldsymbol{\varphi}}(\mathbf{g}), \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}')\right] &=& \int \mathrm{d}h \prod\_{1}^{d} \delta\left(\mathcal{g}\_{i} \mathrm{h}(\mathcal{g}\_{i}'^{-1})\right) = \mathbbm{1}\_{\mathbb{G}}(\mathcal{g}\_{i}, \mathbb{g}\_{i}'),\\ \left[\hat{\boldsymbol{\varphi}}(\mathbf{g}), \hat{\boldsymbol{\varphi}}(\mathbf{g}')\right] &=& \left[\hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}), \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}')\right] = 0. \end{aligned} \tag{5}$$

In these terms, the Fock vacuum is the state with no quantum geometrical or matter degrees of freedom, satisfying *ϕ*<sup>+</sup>(**g**)|<sup>0</sup> = 0 for all arguments. A generic single-particle state |*φ* with wavefunction *φ*, consisting of a *d*-valent node with links labelled by group elements **g** ≡ (*g*1, ..., *gd*), is written as

$$
\langle \Phi \rangle = \int\_{G^{\times d}} \mathbf{dg} \,\phi(\mathbf{g}) \, | \mathbf{g} \rangle \,\tag{6}
$$

where d**g** = ∏*<sup>d</sup> i*=1 d*gi* is the Haar measure, *φ* is an element of the single-particle Hilbert space H = *<sup>L</sup>*<sup>2</sup>(*G*<sup>×</sup>*<sup>d</sup>*), and |**g** = |*g*1 × ... × |*gd* a basis (of Dirac distributions) in H. For *d* = 4, this is the space of states of a quantum tetrahedron [52]. Notice then that one can think of (6) as the analogue of a quantum tensor state where each group element *g* corresponds to an index *i* variable in a continuous (∞-dim) index space H*i* = *<sup>L</sup>*<sup>2</sup>(*G*). 3

The complete Fock space is given by the direct sum of *n*-particle sectors H⊗*n*, restricted to states that are invariant under graph automorphisms of vertex relabelling in the spin network picture, in order to consider multi-particle states that only depend on the intrinsic combinatorial structure of their interaction pattern (a discrete counterpart of continuum diffeomorphisms). Therefore, one has

$$\mathbb{F} \equiv \bigoplus\_{n \ge 0} \text{sym}[\mathbb{H}^{\otimes n}].\tag{9}$$

The symmetry condition is consistent with the assumed *bosonic statistics* and it implies *indistinguishability* for the quanta of the quantum many-body system.

Generic GFT observables <sup>O</sup><sup>+</sup>[*ϕ*<sup>+</sup>, *ϕ*<sup>+</sup>†] in the Fock space are defined in terms of a series of many-body operators expressed as a function of the field operators (or of the basic creation/annihilation operators). For instance, a (*n*)-body operator O*n* acting on p vertices and resulting in q particles is written as [30]

$$\hat{\mathcal{O}}\_{\scriptscriptstyle\rm n}[\hat{\boldsymbol{\varrho}}, \hat{\boldsymbol{\varrho}}^{\dagger}] = \sum\_{n=2}^{\infty} \sum\_{\substack{p+q=n\\p,q\ge 0}} \int \prod\_{p=1}^{p} \mathbf{dg}\_{\textit{p}}^{\prime} \, \hat{\boldsymbol{\varrho}}^{\dagger}(\mathbf{g}\_{\textit{p}}^{\prime}) \prod\_{q=1}^{q} \mathbf{dg}\_{\textit{q}} \, \hat{\boldsymbol{\varrho}}(\mathbf{g}\_{\textit{q}}) \, \mathcal{O}\_{\textit{n}}(\mathbf{g}\_{1\prime}^{\prime}, \cdots, \mathbf{g}\_{p\prime}^{\prime}; \mathbf{g}\_{1\prime}, \cdots, \mathbf{g}\_{q\prime}), \tag{10}$$

where O*n*(**g** 1, ··· , **g** *p*; **g**1, ··· , **<sup>g</sup>***q*) denote the matrix elements of a corresponding first-quantized operator.

### *3.1. Multi-Particle State Observables*

Analogous with the case of the single particle state in (6), we can think of a quantum many-body system as a collective state generated by the action of a *multi-particle group-field operator* on the Fock vacuum. We can define a **product** *n*-particle state, comprising *n* disconnected nodes, by the multiple action of the creation field operator in the group representation of the Fock space, e.g.,

$$\left< \mathbf{g}\_1, \mathbf{g}\_2, \dots, \mathbf{g}\_n \right> = \frac{1}{\sqrt{n!}} \prod\_{a=1}^n \hat{\boldsymbol{\rho}}^\dagger(\mathbf{g}\_a) \left| 0 \right> . \tag{11}$$

$$
\langle j\_1, \ldots j\_d | i\_1, \ldots, i\_d \rangle = \delta\_{i\_1 i\_1} \cdot \ldots \cdot \delta\_{i\_d i\_d} \,. \tag{7}
$$

With respect to this basis, we decompose a tensor *T* into its components *Ti*1...*id*, which introduces an isomorphism to C*<sup>N</sup><sup>d</sup>* :

$$\left| \left| T \right> \right> = \mathop{\rm s.t.}\_{i\_1, \dots, i\_d \in \mathbb{Z}\_{\forall d}} T\_{i\_1 \dots i\_d} \left| i\_1, \dots, i\_d \right> \,. \tag{8}$$

<sup>3</sup> The analogy with a tensor state is apparent again for *G* = *ZN*. Let H*d*,*<sup>N</sup>* be the space of tensors with fixed rank *d* and index cardinality *N*. Neglecting the structure of the cyclic group, H*d*,*<sup>N</sup>* is reduced to C*<sup>N</sup><sup>d</sup>* . The linear structure, the scalar product and the completeness of C*<sup>N</sup><sup>d</sup>* establish H*d*,*<sup>N</sup>* to be a Hilbert space. A basis of H*d*,*<sup>N</sup>* is chosen by |*<sup>i</sup>*1, ..., *id* , defined as:

*Universe* **2019**, *5*, 211

Because the *<sup>ϕ</sup>*<sup>+</sup>†(**g***a*) commute with each other, the order of the particles does not affect the state and we have |..., **g***<sup>a</sup>*, ..., **g***b*, ... = |..., **g***b*, ..., **g***<sup>a</sup>*, ... . The *n*-particle state |**g**1, ..., **g***n* defines a multi-particle basis of the Fock space F, with orthogonality relation given by

$$\langle \mathbf{g}'\_{1'}, \dots, \mathbf{g}'\_{n} | \mathbf{g}\_{1'}, \dots, \mathbf{g}\_{n'} \rangle = \frac{1}{n!} \sum\_{\pi \in S\_n} \prod\_{a=1}^{n} \int \, \mathrm{d}h\_a \, \delta^4(\mathbf{g}'\_a h\_a \mathbf{g}\_{\pi(a)}^{-1}) \,. \tag{12}$$

In addition, because of the required symmetry of the fields *ϕ*<sup>+</sup>(**g***h*) = *ϕ*<sup>+</sup>(**g**), *h* ∈ *G*, the *n*-particle state is right invariant

$$\left| \cdot \cdot \right\rangle \left. \mathbf{g}\_a h\_{a\prime} \cdot \cdot \cdot \right\rangle = \left| \cdot \cdot \right\rangle \left. \mathbf{g}\_{a\prime} \cdot \cdot \right\rangle \cdot \tag{13}$$

Finally, the resolution of identity in the Fock space can be written in terms of the *n*-particle state as

$$\mathbb{1}\_{\mathcal{F}} = |0\rangle\langle 0| + \sum\_{n=1}^{\infty} \int \prod\_{a=1}^{n} \mathbf{dg}\_{a} \left| \mathbf{g}\_{1\prime} \cdots \mathbf{\prime}, \mathbf{g}\_{n} \right\rangle \left\langle \mathbf{g}\_{1\prime} \cdots \mathbf{\prime}, \mathbf{g}\_{n} \right| \,. \tag{14}$$

One can check immediately that 1F |**g**1, ··· , **g***n* = |**g**1, ··· , **g***n* .

A generic GFT **multi-particle state**, analogous to (6), will then be based on a specific configuration of *n* fields, characterized by a multi-particle wavefunction Ψ*<sup>n</sup>*, generated by a multi-particle operator in the GFT Fock space

$$\hat{\Psi}\_{\boldsymbol{n}}[\hat{\boldsymbol{\varphi}}^{\dagger}] \quad = \sum\_{\substack{\boldsymbol{p}+\boldsymbol{q}=\boldsymbol{n} \\ \boldsymbol{p}, \boldsymbol{q}\geq 0}} \int \prod\_{\mathbf{p}=1}^{\boldsymbol{p}} \mathrm{d} \mathbf{g}\_{\boldsymbol{p}}^{\prime} \, \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}\_{\boldsymbol{p}}^{\prime}) \prod\_{\mathbf{q}=1}^{\boldsymbol{q}} \mathrm{d} \mathbf{g}\_{\boldsymbol{q}} \, \hat{\boldsymbol{\varphi}}(\mathbf{g}\_{\boldsymbol{q}}) \, \Psi\_{\boldsymbol{n}}(\mathbf{g}\_{1}^{\prime}, \cdots, \mathbf{g}\_{\boldsymbol{p}}^{\prime}; \mathbf{g}\_{1}, \cdots, \mathbf{g}\_{\boldsymbol{q}}) \tag{15}$$

via the repeated action of the GFT field operators on the Fock space.

We are interested in the structure of quantum correlations of the multi-particle state <sup>Ψ</sup>*n*(**g**1, ··· , **g***p*; **g**1, ··· , **<sup>g</sup>***q*). As it is the case for any highly entangled quantum many-body system, disentangling the information on the quantum correlations of Ψ*n* is highly nontrivial. Therefore, we focus on a special class of multi-particle operators Ψ + Γ, where the wave-function is explicitly constructed via a pairwise contractions scheme of single node states, in correspondence with a given **network architecture** Γ. We write

$$\Psi\_{\Gamma}[\hat{\boldsymbol{\varphi}}, \hat{\boldsymbol{\varphi}}^{\dagger}](\mathbf{g}\_{\boldsymbol{\theta}}) = \sum\_{\substack{p+q=n\\p,q\ge 0}} \int \prod\_{\mathbf{p}=1}^{p} \mathbf{d} \mathbf{g}\_{\mathbf{p}}^{\prime} \, \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}\_{\mathbf{p}}^{\prime}) \prod\_{q=1}^{q} \mathbf{d} \mathbf{g}\_{\mathbf{q}} \, \hat{\boldsymbol{\varphi}}(\mathbf{g}\_{\mathbf{q}}) \prod\_{\ell\in\Gamma}^{L} \mathbf{d}h\_{\ell} \prod\_{\ell\in\Gamma} L\_{\ell} \left(\mathbf{g}\_{\mathbf{s}(\ell)}^{\prime} h\_{\ell} \mathbf{g}\_{\mathbf{t}(\ell)}^{-1}\right) \tag{16}$$

with *ϕ*+† operators generating the nodes connected by link kernels *L* to form an open network with **g***∂* dangling indices, via an overall integration over **g**∈Γ. The expression of the link convolution kernel connecting the nodes pairwise is left generic at this stage, with the only requirement to preserve the overall gauge invariance of the network state.

We shall see such a class of multi-particle operators as *tensor networks* operators, where we reduce the entanglement structure of the multi-particle state to local correlations induced by the generic link kernels *L*, propagated non-locally via nodes.

Notice that the most generic tensor network state of the theory would involve superpositions of both network architectures (combinatorial structures) and number of particles corresponding to the same number of boundary degrees of freedom,

$$\hat{\Psi}[\hat{\boldsymbol{\varphi}},\hat{\boldsymbol{\varphi}}^{\dagger}](\mathbf{g}\_{\mathbf{b}}) = \sum\_{\{\Gamma\}} \sum\_{\begin{subarray}{c} \mathbf{p}=\mathbf{p}\end{subarray}}^{\infty} \sum\_{\begin{subarray}{c} p+q=n\\p,q\geq 0\end{subarray}} \int \prod\_{\mathbf{p}=1}^{p} \mathbf{dg}\_{\mathbf{p}}^{\prime} \,\hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}\_{\mathbf{p}}^{\prime}) \prod\_{\mathbf{q}=1}^{q} \mathbf{dg}\_{\mathbf{q}} \,\hat{\boldsymbol{\varphi}}(\mathbf{g}\_{\mathbf{q}}) \,\mathcal{L}\_{p,q}^{\Gamma}(\mathbf{g}\_{1}^{\prime},\cdots,\mathbf{g}\_{p}^{\prime};\mathbf{g}\_{1},\cdots,\mathbf{g}\_{q}),\tag{17}$$

where, for simplicity, we indicate by L<sup>Γ</sup> *p*,*q* the set of pairwise link convolutions (gluing) functions associated with a given graph Γ, and some suitable symmetry quotient factor removing equivalent graph configurations is assumed.

### *3.2. GFT Coherent State Basis*

In the context of quantum gravity, we specify the GFT formalism to the case *G* = *SU*(2), the relevant local gauge subgroup of gravity, and we understand the group elements *g* as a generalisation of the embedded parallel transports (holonomies) of the gravitational *G*-connection of loop quantum gravity. The symmetric GFT *d*-valent fields as *d*-simplices (convex polyhedra, e.g., for *d* = 4, these are tetrahedra), with *d* number of faces labelled by dual Lie algebra-valued flux variables, become single "quanta" of twisted geometry states expressed in terms of quantum spin network basis [29].

In this setting, on the one hand, we are interested in working with states in the Fock space that can be eventually put in relation with extended macroscopic 3D geometries. To this aim, the natural choice consists of looking for a coherent state basis in F, defined by exponential operators providing desirable coherence properties, having *macroscopic* occupation numbers for given modes controlled by the wave-function [53–56]. More concretely, this choice will allow us to compute quantum averages of many-body systems in thermal equilibrium using functional integrals over group field configurations [18–20].

The simplest class of such states is given by the single-particle (condensate) coherent states

$$\begin{split} |\langle \boldsymbol{\varrho} \rangle| &\equiv \frac{1}{\mathcal{N}\_{\boldsymbol{\Phi}}} \exp \left[ \int \mathbf{dg} \,\boldsymbol{\Phi}(\mathbf{g}) \hat{\boldsymbol{\phi}}^{\dagger}(\mathbf{g}) \right] |0\rangle \equiv \frac{1}{\mathcal{N}\_{\boldsymbol{\Phi}}} \sum\_{n=0}^{\infty} \frac{1}{n!} \prod\_{a}^{n} \left[ \int \mathbf{dg}\_{\boldsymbol{\Phi}} \,\boldsymbol{\Phi}(\mathbf{g}\_{\boldsymbol{\Phi}}) \hat{\boldsymbol{\phi}}^{\dagger}(\mathbf{g}\_{\boldsymbol{\Phi}}) \right] |0\rangle \\ &\equiv \quad \frac{1}{\mathcal{N}\_{\boldsymbol{\Phi}}} \sum\_{n=0}^{\infty} \frac{1}{\sqrt{n!}} \int [\mathbf{dg}]^{n} \,\boldsymbol{\Phi}(\mathbf{g}\_{1}) \times \dots \times \boldsymbol{\Phi}(\mathbf{g}\_{n}) | \mathbf{g}\_{1}, \dots, \mathbf{g}\_{n} \rangle. \end{split} \tag{18}$$

For the last equality, we use the definition of the *n*-particle state (11). *φ*(**g**) is the field on H that has the same gauge symmetry as *<sup>ϕ</sup>*<sup>+</sup>†(**g**), namely *φ*(**g***h*) = *φ*(**g**), and N*φ* is the normalization 4

$$\mathcal{N}\_{\Phi}^{2} = \exp\left[\int \mathrm{dg}\,\overline{\phi(\mathbf{g})}\phi(\mathbf{g})\right]. \tag{22}$$

One can show that |*ϕ* is the eigenstate of the field operator *ϕ*<sup>+</sup>(**g**) such that

$$
\langle \hat{\boldsymbol{\varphi}}(\mathbf{g}) \mid \boldsymbol{\varphi} \rangle = \boldsymbol{\phi}(\mathbf{g}) \mid \boldsymbol{\varphi} \rangle \,. \tag{23}
$$

4 Define an operator *A*+ as

> The commutator between *A*+ and *A*<sup>+</sup>† is

$$
\widehat{A} = \int \mathrm{d}\mathbf{g} \,\,\phi(\mathbf{g}) \widehat{\phi}^{\dagger}(\mathbf{g}).\tag{19}
$$

$$[\hat{A}^\dagger, \hat{A}] = \int \text{d}\mathbf{g} \,\overline{\phi(\mathbf{g})} \phi(\mathbf{g}).\tag{20}$$

Then, *ϕ*|*ϕ* can be given as

$$\begin{split} 1 &= \quad \langle \boldsymbol{\varrho} | \boldsymbol{\varrho} \rangle = \mathcal{N}\_{\boldsymbol{\varrho}}^{-2} \langle 0 | \, \mathbf{e}^{\hat{A}^{\dagger}} \mathbf{e}^{\hat{A}} | 0 \rangle \\ &= \quad \mathcal{N}\_{\boldsymbol{\varrho}}^{-2} \langle 0 | \, \mathbf{e}^{\hat{A}} \mathbf{e}^{\hat{A}^{\dagger}} \mathbf{e}^{[\hat{A}^{\dagger}, \hat{A}]} | 0 \rangle \\ &= \quad \mathcal{N}\_{\boldsymbol{\varrho}}^{-2} \exp \left[ \int \, \mathbf{d} \mathbf{g} \, \overline{\boldsymbol{\phi}(\mathbf{g})} \boldsymbol{\phi}(\mathbf{g}) \right]. \end{split} \tag{21}$$

In the third equality, we use the Baker–Campbell–Hausdorff formula. Indeed, we have

$$\hat{\boldsymbol{\varphi}}(\mathbf{g}) \frac{1}{n!} \prod\_{a}^{n} \left[ \int \mathrm{d} \mathbf{g}\_{a} \, \boldsymbol{\phi}(\mathbf{g}\_{a}) \hat{\boldsymbol{\phi}}^{\dagger}(\mathbf{g}\_{a}) \right] \tag{24}$$

$$=\frac{1}{n!} \sum\_{k=1}^{n} \prod\_{a \neq k}^{n} \left[ \int \text{dg}\_{a} \, \phi(\mathbf{g}\_{a}) \hat{\phi}^{\dagger}(\mathbf{g}\_{a}) \right] \int \text{d}h \, \text{dg}\_{k} \, \phi(\mathbf{g}\_{k}) \delta(\mathbf{g} \text{lg}\_{k}^{-1}) \tag{25}$$

$$=\frac{1}{(n-1)!} \prod\_{a}^{n-1} \left[ \int \mathrm{d} \mathbf{g}\_{a} \, \phi(\mathbf{g}\_{a}) \hat{\phi}^{\dagger}(\mathbf{g}\_{a}) \right] \int \mathrm{d} \mathbf{h} \, \phi(\mathbf{g} \, \mathbf{h}) \tag{26}$$

$$=\left.\phi(\mathbf{g})\frac{1}{(n-1)!}\prod\_{a}^{n-1}\left[\int \mathrm{d}\mathbf{g}\_{a}\,\phi(\mathbf{g}\_{a})\hat{\phi}^{\dagger}(\mathbf{g}\_{a})\right].\tag{27}$$

In the first equality, we use the commutator (5) between *ϕ*+ and *ϕ*<sup>+</sup>†. In the last equality, we use the fact that *φ*(**g**) is right invariant. Thus, when *ϕ*<sup>+</sup>(**g**) acts on |*ϕ* , it gives (23). In particular, coherent states |*ϕ* provide an over-complete basis of the Fock space F (see Appendix A for details).

Via Equation (23), one immediately obtains the tensor fields *φ*(**g**) and *φ*(**g**) in terms of expectation values of *ϕ*<sup>+</sup>(**g**) and *ϕ*<sup>+</sup>†(**g**) with respect to |*ϕ* 

$$
\langle \boldsymbol{\varrho} | \widehat{\boldsymbol{\varrho}}(\mathbf{g}) | \boldsymbol{\varrho} \rangle = \boldsymbol{\Phi}(\mathbf{g}), \quad \langle \boldsymbol{\varrho} | \widehat{\boldsymbol{\varrho}}^{\dagger}(\mathbf{g}) | \boldsymbol{\varrho} \rangle = \overline{\boldsymbol{\Phi}(\mathbf{g})}.\tag{28}
$$

Accordingly, we can express the multi-particle state as the expectation value of a group-field network operator (16) in the 2nd quantised basis of eigenstates of the GFT quantum field operator. +

For the network operator Ψ <sup>Γ</sup>[*ϕ*<sup>+</sup>, *<sup>ϕ</sup>*<sup>+</sup>†](**g***∂*) defined in Equation (16), we ge<sup>t</sup>

$$
\begin{split}
\langle\boldsymbol{\varrho}|\,^{\Psi}\boldsymbol{\Gamma}\_{\Gamma}[\widehat{\boldsymbol{\varrho}},\widehat{\boldsymbol{\varrho}}^{\dagger}]\,^{\mathop{\!}}\boldsymbol{\varrho}\rangle &=&\sum\_{\begin{subarray}{c}p+q=\boldsymbol{\imath}\\p,q\geq 0\end{subarray}}\int\prod\_{p=1}^{p}\mathrm{d}\mathbf{g}\_{p}^{\prime}\,\overline{\boldsymbol{\phi}(\mathbf{g}\_{p}^{\prime})}\prod\_{q=1}^{q}\mathrm{d}\mathbf{g}\_{q}\,\boldsymbol{\phi}(\mathbf{g}\_{q})\prod\_{\ell\in\Gamma}^{L}\mathrm{d}\boldsymbol{h}\_{\ell}\,\prod\_{\ell\in\Gamma}L\_{\ell}\left(\mathrm{g}\_{\mathbf{s}(\ell)}h\_{\ell}\mathbf{g}\_{\mathbf{t}(\ell)}^{-1}\right) \\
&=&\Psi\_{\Gamma}[\overline{\boldsymbol{\Phi}},\boldsymbol{\Phi}](\mathbf{g}\_{\boldsymbol{\Theta}}).
\end{split}
\tag{29}
$$

This is a group field tensor network state based on graph Γ with *n* nodes and *L* links. In particular, for *L hs*()*gh*−<sup>1</sup> *t*() = *δ hs*()*gh*−<sup>1</sup> *t*(), we can think of the expectation values of the associated operators as peculiar projected entangled-pairs tensor network states (PEPS)

$$\left|\Psi\_{\Gamma}\right\rangle \equiv \int \mathbf{dg}\_{\partial} \,\Psi\_{\Gamma}(\mathbf{g}\_{\partial}) \left|\mathbf{g}\_{\partial}\right\rangle \equiv \bigotimes\_{\ell \in \Gamma} \langle L\_{\ell} \vert \bigotimes\_{v \in \Gamma} \left|\phi\_{v}\right\rangle \tag{30}$$

obtained by the contraction of maximally entangled link states

$$|L\_{\ell}\rangle = \int \mathbf{dg}\_{s(\ell)} \mathbf{dg}\_{t(\ell)} \, \delta\left(\mathcal{g}\_{s(\ell)} h\_{\ell} \mathbf{g}\_{t(\ell)}^{-1}\right) |\mathcal{g}\_{s(\ell)}\rangle \otimes \left|\mathcal{g}\_{t(\ell)}\right\rangle \tag{31}$$

via tensor states |*φ* on some generically open graph architecture Γ.

The basis |**g***∂* labels the uncontracted dangling indices comprising the boundary of the auxiliary tensor network representation. Differently from standard PEPS, the GFT networks are further characterised by the inherently *random* character of the tensors |*φ* , induced by their field-theoretic statistical descritpion. In this light, in particular, we can see states like (29) as a generalisation of the random tensor network states (RTNs) recently introduced in [15], where the statistical characterisation of the standard RTN gets mapped in the momenta of the GFT partition function.<sup>5</sup>

<sup>5</sup> When Γ is closed, one further recognises ΨΓ[*φ*, *φ*](**g**∈<sup>Γ</sup>) to be equivalent to the cylindrical functions describing the quantum geometry of a closed spacial hypersurface in the kinematic Hilbert space of LQG [57].

### **4. Bipartite Entanglement of a GFT Network State**

Given a group field tensor network state |ΨΓ ∈ H*∂* = /∈*<sup>∂</sup> L*2 [*G*], a bipartition of the boundary degrees of freedom corresponds to a factorisation of the boundary Hilbert space H*∂* into two subspaces H*A* and H*B*, such that

$$\mathbb{H}^{\mathfrak{g}} = \mathbb{H}\_A \otimes \mathbb{H}\_B. \tag{32}$$

The entanglement of the boundary state |ΨΓ across the bipartition in H*A* and H*B* is measured by the von Neumann entropy

$$S(A) = -\text{Tr}\phi\_A \ln \phi\_{A\prime} \tag{33}$$

where

$$\rho\_A \equiv \frac{\rho\_A}{\text{Tr}\rho'} \quad \rho\_A \equiv \text{Tr}\_B \rho = \text{Tr}\_B \left| \Psi\_\Gamma \right\rangle \left\langle \Psi\_\Gamma \right| \tag{34}$$

defines the (normalised) marginal on H*A* of the GFT multi-particle density matrix,

$$\rho = |\Psi\_{\Gamma}\rangle\langle\Psi\_{\Gamma}| = \text{Tr}\left[\bigotimes\_{\ell}|L\_{\ell}\rangle\langle L\_{\ell}|\bigotimes\_{v}|\phi\_{v}\rangle\langle\phi\_{v}|\right] \equiv \text{Tr}\bigotimes\_{\ell}\rho\_{\ell}\bigotimes\_{v}\rho\_{v}.\tag{35}$$

A representation of *ρ* for the case of a simple 2-vertices graph is given in Figure 1.

It is computationally convenient to derive the von Neumann entropy as the limit of the Rényi entropy *SN*(*A*), via standard replica trick. The Rényi entropy is defined as

$$S\_N(A) = \frac{1}{1 - N} \ln \text{Tr} \bar{\rho}\_A^N = \frac{1}{1 - N} \ln \frac{\text{Tr} \rho\_A^N}{(\text{Tr} \rho)^N} \equiv \frac{1}{1 - N} \ln \frac{\text{Tr} (\rho^{\otimes N} \mathbb{P}\_A)}{\text{Tr} (\rho^{\otimes N})},\tag{36}$$

where P*A* = <sup>P</sup>(*π*0*A*; *n*, *d*) is the 1-cycle permutation operator in *SN* acting on the reduced Hilbert space H*A*,

$$\mathbb{P}(\pi\_A^0; N, d) = \prod\_{s=1}^N \delta\_{\mu\_A^{([s+1]\_D)} \mu\_A^{(s)} \prime} \tag{37}$$

and *d* is the dimension of the Hilbert space in the same region *A* (see Figure 1). Explicitly, one has

$$\left| \mathbb{P}\_A \left| a\_1, b\_1 \right> \left| a\_2, b\_2 \right> \cdots \left| a\_N, b\_N \right> = \left| a\_2, b\_1 \right> \left| a\_3, b\_2 \right> \cdots \left| a\_1, b\_N \right> \tag{38}$$

with /*i* |*ai* ∈ H*A* and /*i* |*bi* ∈ H*B*.

> The Rényi entropy *SN*(*A*) coincides with the von Neumann entropy *S*(*A*) as *N* goes to 1, which is

$$S(A) = \lim\_{N \to 1} S\_N(A). \tag{39}$$

We are interested in carrying on this measurement directly in the Fock space of GFT.

**Figure 1.** Graphical representation of the density matrix state for a simple bipartite GFT network state |G , comprising two bivalent internal nodes, and trace of the *N* = 2 replica with the action of the cyclic permutation (swap) operator. Notice that the label *A* ¯ , generally indicating the complementary marginal, corresponds to the *B*-labelling in the main text.

### *4.1. Expected Rènyi Entropy in the Fock Space*

Due to the random character of the nodes {*φv*}, induced by their dynamical GFT description, the measure of the entanglement will be necessarily given in expectation value.

Now, by construction, the Rényi entropy *SN*(*A*) is a functional of the GFT fields *φ* and *φ*, hence it can be promoted to an observable *SN*(*A*)[*ϕ*<sup>+</sup>, *ϕ*<sup>+</sup>†] in the GFT Fock space. We shall then derive the expectation value of the Rényi entropy in the Fock space using single-particle coherent state basis, by inserting the resolution of identity 1F ( and set *K* = 0, see Appendix A), with a normal ordering **:** ··· **:** such that all *ϕ*+† is to the left of *ϕ*<sup>+</sup>:<sup>6</sup>

$$\mathbb{E}\left[\mathcal{S}\_N(A)[\hat{\boldsymbol{\varphi}}, \hat{\boldsymbol{\varphi}}^\dagger]\right] \equiv \frac{\mathbb{C}}{Z\_0} \text{Tr}\left(: \mathcal{S}\_N(A)[\hat{\boldsymbol{\varphi}}, \hat{\boldsymbol{\varphi}}^\dagger] \, \mathrm{e}^{-\bar{\mathcal{S}}[\hat{\boldsymbol{\varphi}}, \hat{\boldsymbol{\varphi}}^\dagger]} : \mathbb{1}\_\mathcal{F}\right),\tag{41}$$

where the very GFT action *S* + [*ϕ*+, *ϕ*<sup>+</sup>†] is constructed as a quantum many-body operator on the Fock space. By (36), the explicit form of the expectation value reads

$$\mathbb{E}[S\_N(A)[\hat{\boldsymbol{\wp}}, \hat{\boldsymbol{\wp}}^{\dagger}]] = \frac{1}{1 - N} \mathbb{E}\left[\ln \frac{\text{Tr}(\boldsymbol{\rho} \, [\hat{\boldsymbol{\wp}}, \hat{\boldsymbol{\wp}}^{\dagger}]^{\otimes N} \mathbb{P}\_A)}{\text{Tr}(\boldsymbol{\rho} \, [\hat{\boldsymbol{\wp}}, \hat{\boldsymbol{\wp}}^{\dagger}])^N}\right] \quad . \tag{42}$$

If we choose to work in perturbative GFT regime, such that *<sup>S</sup>*[*φ*, *φ*] = *S*0 + *<sup>λ</sup>S*int[*φ*, *φ*], with *λ* 1, then we deal with a polynomially perturbed generalised [45] Gaussian distribution for the random field, hence we can take advantage of the central limit theorem to ge<sup>t</sup> a good approximation of (42) in terms of a Taylor expansion around the mean values

$$\mathbb{E}[S\_N(A)[\hat{\boldsymbol{\wp}}, \hat{\boldsymbol{\wp}}^{\dagger}]] \approx \frac{1}{1-N} \ln \frac{\mathbb{E}[\text{Tr}(\rho \, [\hat{\boldsymbol{\wp}} \, \hat{\boldsymbol{\wp}}^{\dagger} \, ]^{\otimes N} \mathbb{P}\_A)]}{\mathbb{E}[\text{Tr}(\rho \, [\hat{\boldsymbol{\wp}} \, \hat{\boldsymbol{\wp}}^{\dagger} \, ]^{\otimes N})]} \quad , \tag{43}$$

$$\begin{split} \mathbf{Z}\_{0} &= \quad \mathsf{CTr} \left( \mathbf{:} \mathbf{e}^{-S[\hat{\boldsymbol{\phi}}, \hat{\boldsymbol{\phi}}^{\dagger}]} : \right) = \mathsf{CTr} \left( \mathbf{:} \mathbf{e}^{-S[\hat{\boldsymbol{\phi}}, \hat{\boldsymbol{\phi}}^{\dagger}]} : \mathbb{1}\_{\mathrm{F}} \right) \\ &= \quad \int \mathcal{D}\boldsymbol{\phi} \mathcal{D}\overline{\boldsymbol{\phi}} \, \delta \mathbb{C} \left[ \boldsymbol{\phi} - \boldsymbol{\phi}\_{\mathrm{CF}} \right] \, \mathbf{e}^{-S[\hat{\boldsymbol{\phi}}, \overline{\boldsymbol{\phi}}]} . \end{split} \tag{40}$$

<sup>6</sup> The GFT vacuum amplitude in the same basis reads

plus corrections due to fluctuations around the mean values that are suppressed whenever the dimensionality of our quantum system gets extremely large. For standard random tensor network states [15], such a *typical* regime is realised in the *large bond* limit, when the dimension of the random tensor index (bond) space H is extremely large, dim(H) 1.

In the GFT setting, one has to deal with infinite dimensional bond spaces, H = *<sup>L</sup>*<sup>2</sup>(*G*), which automatically set the derivation in the typicality regime of [15]. Nevertheless, link spaces are regularised via a cut-off on the group representation space, such that

$$<\langle \mathbf{g} | \mathbf{g}' \rangle = D(\Lambda) \delta\_{\mathbf{g}, \mathbf{g}'} \tag{44}$$

with *<sup>δ</sup>g*,*g* equal to 1 if *g* = *g* and 0, otherwise. Therefore, with *<sup>D</sup>*(Λ), the dimension of the bonds (links) of the network, one can consistently assume typicality to hold in the *D*(Λ) 1 regime.

Given (43), the expression of the *N*-th Rényi entropy is mapped into a ratio of averaged partition functions,

$$\frac{Z\_N}{Z\_0^N} \equiv \frac{\mathbb{E}\left[\text{Tr}(\rho\left[\hat{\boldsymbol{\varphi}}, \hat{\boldsymbol{\varphi}}^\dagger\right]^{\otimes N} \mathbb{P}\_A)\right]}{\mathbb{E}\left[\text{Tr}(\rho\left[\hat{\boldsymbol{\varphi}}, \hat{\boldsymbol{\varphi}}^\dagger\right]^{\otimes N})\right]},\tag{45}$$

where we removed the bars over *ZN*, *<sup>Z</sup>N*0 to simplify the notation.

Let us rewrite *ρ* [*ϕ*<sup>+</sup>, *ϕ*<sup>+</sup>†] ⊗*N* as a trace contraction of individual link and nodes density matrices, similarly to (35). We restrict for simplicity to tensor network observables (16) with *n* vertices and *L* links, and with *p* = *n*(*q* = 0)

$$\hat{\Psi}\_{\Gamma}[\hat{\boldsymbol{\varphi}}^{\dagger}](\mathbf{g}\_{\mathfrak{d}}) = \int \prod\_{\mathfrak{p}=1}^{n} \mathbf{d} \mathbf{g}\_{\mathfrak{p}} \prod\_{\mathfrak{p}=1}^{n} \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}\_{\mathfrak{p}}) \prod\_{\ell \in \Gamma}^{L} \mathbf{d}h\_{\ell} \prod\_{\ell \in \Gamma}^{L} L\_{\ell} \left(\mathbf{g}\_{\mathfrak{e}(\ell)} h\_{\ell} \mathbf{g}\_{\mathfrak{t}(\ell)}^{-1}\right) \tag{46}$$

where now the p label coincides with the *v*-label of the nodes of the graph Γ. The *Nth* replica of the density matrix operator describing a graph observable reads

$$\boldsymbol{\rho}^{\otimes N} \quad = \left( \int \prod\_{v=1}^{n} \mathbf{d} \mathbf{g}\_{v} \mathbf{d} \mathbf{g}\_{v}^{\prime} \prod\_{v=1}^{n} \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}\_{v}) \, \hat{\boldsymbol{\varphi}}(\mathbf{g}\_{v}^{\prime}) \prod\_{\ell}^{L} \mathbf{d} h\_{\ell} \, \mathbf{d} h\_{\ell}^{\prime} \prod\_{\ell}^{L} \mathbf{L}\_{\ell} L\_{\ell}^{\prime} \right)^{\otimes N} \tag{47}$$

$$=\int \left(\prod\_{\upsilon=1}^{n} \mathbf{d} \mathbf{g}\_{\upsilon} \mathbf{d} \mathbf{g}\_{\upsilon}^{\prime} \prod\_{\upsilon=1}^{n} \hat{\boldsymbol{\varphi}}^{\dagger}(\mathbf{g}\_{\upsilon}) \,\hat{\boldsymbol{\varphi}}(\mathbf{g}\_{\upsilon}^{\prime})\right)^{\otimes N} \left(\prod\_{\ell}^{L} \mathbf{d} \boldsymbol{h}\_{\ell} \,\mathbf{d} \mathbf{h}\_{\ell}^{\prime} \prod\_{\ell}^{L} \mathbf{L}\_{\ell} \boldsymbol{L}\_{\ell}^{\prime}\right)^{\otimes N} \tag{48}$$

$$= \left[ \text{Tr} \left[ (\bigotimes\_v \rho\_v)^{\otimes N} (\bigotimes\_\ell \rho\_\ell)^{\otimes N} \right] . \right] . \tag{49}$$

The linearity of the trace allows for moving the expectation operator inside the integral and letting it act on the N replicas of the products of fields. We then write (45) as

$$\frac{\text{Tr}\left[\bigotimes\_{\ell}\rho\_{\ell}^{N}\mathbb{E}[\left(\bigotimes\_{v}\rho\big[\hat{\boldsymbol{\rho}}^{\dagger}(\mathbf{g}\_{\nu})\big,\hat{\boldsymbol{\rho}}(\mathbf{g}\_{v}^{\prime})\big]\big|\big,\_{v}^{\otimes N}\right]\mathbb{P}\_{A}\right]}{\text{Tr}\left[\bigotimes\_{\ell}\rho\_{\ell}^{N}\mathbb{E}[\left(\bigotimes\_{v}\rho\big[\hat{\boldsymbol{\rho}}^{\dagger}(\mathbf{g}\_{\nu})\big,\hat{\boldsymbol{\rho}}(\mathbf{g}\_{v}^{\prime})\big]\big|\big}\_{v}\right]^{\otimes N}\right]}\tag{50}$$

and we focus on the calculation of <sup>E</sup>[/*v ρ* ⊗*Nv* ].

The derivation in [15,18] thereby proceeds by making two strongly simplifying assumptions. The first consists of the restriction to the case of a Gaussian (free) group field theory for describing the single node field statistics. The second assumes the expected value *φ* of the field operator at each node to be individually independently distributed (i.i.d.). The latter assumption corresponds, from a

physical viewpoint, to considering a non-interacting quantum many body system. The latter condition translates in particular into a *local averaging* condition, namely

$$\mathbb{E}\left[\left(\bigotimes\_{v} \rho\_{v}\right)^{\otimes N}\right] \to \bigotimes\_{v} \mathbb{E}\left[\rho\_{v}^{\otimes N}\right] \quad , \tag{51}$$

which allows for an explicit calculation of the expectation value in terms of a product of *n* 2*N*-point function of the free group field theory

$$\mathbb{E}\left[\rho\_v^{\otimes N}\right] \equiv \frac{\mathbb{C}}{Z\_0} \text{Tr}\left( : \rho[\hat{\boldsymbol{\phi}}^{\dagger}(\mathbf{g}\_v), \hat{\boldsymbol{\phi}}(\mathbf{g}\_v')]^{\otimes N} \, \mathbf{e}^{-\hat{S}\_0[\hat{\boldsymbol{\phi}}, \hat{\boldsymbol{\phi}}^{\dagger}]} : 1\_{\mathbb{F}}\right) \tag{52}$$

$$= \frac{1}{Z\_0} \int \mathcal{D}\boldsymbol{\phi} \mathcal{D}\overline{\boldsymbol{\phi}} \, \delta\_{\mathbb{C}}[\boldsymbol{\phi} - \boldsymbol{\phi}\_{\text{GF}}] \, \langle \boldsymbol{\uprho} | : [\hat{\boldsymbol{\phi}}^{\dagger}(\mathbf{g}\_v), \hat{\boldsymbol{\uprho}}(\mathbf{g}\_v')]^{\otimes N} \, \mathbf{e}^{-\hat{S}\_0[\hat{\boldsymbol{\phi}}, \hat{\boldsymbol{\uprho}}^{\dagger}]} : |\boldsymbol{\uprho}\rangle$$

$$= \frac{1}{Z\_0} \int \mathcal{D}\boldsymbol{\phi} \mathcal{D}\overline{\boldsymbol{\uprho}} \, \delta\_{\mathbb{C}}[\boldsymbol{\phi} - \boldsymbol{\uprho}\_{\text{GF}}] \, (\boldsymbol{\uprho}, \overline{\boldsymbol{\uprho}})^N \, \mathbf{e}^{-S\_0[\hat{\boldsymbol{\uprho}}, \hat{\boldsymbol{\uprho}}]} \tag{53}$$

$$\hat{\mathbf{E}} = \begin{bmatrix} \mathbf{E}\_0 \left[ (\phi\_v \vec{\phi}\_v)^N \right] & \mathbf{\bar{I}} \end{bmatrix} \tag{54}$$

In the free case, in particular, one can evaluate the 2*N*-point functions at each node directly via Wick's theorem

$$\begin{split} \mathbb{E}\_{0} \left[ \prod\_{a}^{N} \phi\_{\boldsymbol{v}}(\mathbf{g}\_{\boldsymbol{\pi}}) \overline{\phi\_{\boldsymbol{v}}(\mathbf{g}\_{\boldsymbol{\pi}}')} \right] &= \quad \mathscr{C} \sum\_{\pi\_{\boldsymbol{v}} \in \mathcal{S}\_{N}} \int \prod\_{a}^{N} \mathrm{d}h\_{a} \prod\_{a}^{N} \delta\_{\boldsymbol{v}} \left( h\_{d} \mathbf{g}\_{d} \mathbf{g}\_{\pi(a)}^{\dagger} \right) \\ &= \quad \mathscr{C} \sum\_{\pi\_{\boldsymbol{v}} \in \mathcal{S}\_{N}} \int \prod\_{a}^{N} \mathrm{d}h\_{a} \, \mathbb{P}\_{\mathbf{h}\_{\boldsymbol{v}}}(\boldsymbol{\pi}), \end{split}$$

where the permutation operator acts strandwise (locally on the link spaces)

$$\mathbb{P}\_{\mathbf{g}}(\pi) \equiv \prod\_{a}^{N} \delta \left( h\_{a} \mathbf{g}\_{a} \mathbf{g}\_{\pi(a)}^{\prime} \right) = \prod\_{s=1}^{4} \prod\_{a}^{N} \delta \left( h\_{a} \mathbf{g}\_{s a} \mathbf{g}\_{s \pi(a)}^{\prime} \right) \equiv \prod\_{s}^{4} \mathbb{P}\_{\mathbf{g}}^{s}(\pi),$$

with **g** independent from **g**, *a* labelling the replica order at each node, **h** denoting the set of *ha*, *a* = 1, ··· , *N*.

When *ha* = 1 for all *a* from 1 to *N*,

$$\mathbb{P}\_{\mathbb{1}}(\pi) \quad = \prod\_{a}^{N} \delta \left( \mathbf{g}\_{a} \mathbf{g}\_{\pi(a)}^{\prime} \right) \tag{55}$$
 
$$= \quad \mathbb{P}(\pi; N, D^{4}) = \prod\_{s}^{4} \mathbb{P}^{s}(\pi; N, D^{4}),$$

where <sup>P</sup>(*π*; *N*, *D*<sup>4</sup>) and <sup>P</sup>*<sup>s</sup>*(*π*; *N*, *D*<sup>4</sup>) are the representations of *π* ∈ S*N* on H⊗<sup>4</sup> and H, respectively. The averaged partition functions, *Z N* and *Z N* 0become

$$\begin{split} \mathcal{Z}\_{N} &\approx \mathcal{C}^{V\_{\Gamma}} \sum\_{\pi\_{\nu} \in \mathcal{S}\_{N}} \int \prod\_{v} \mathbf{dh}\_{v} \operatorname{Tr} \left[ \bigotimes\_{\ell} \rho\_{\ell}^{N} \bigotimes\_{v} \mathbb{P}\_{\mathbf{h}\_{v}} (\pi\_{v}) \mathbb{P} (\pi\_{A}^{0}; N, d) \right] \\ &\equiv \quad \mathcal{C}^{V\_{\Gamma}} \sum\_{\pi\_{\nu} \in \mathcal{S}\_{N}} \int \prod\_{v} \mathbf{dh}\_{v} \bigvee\_{\ell} \mathcal{N}\_{\ell} (\mathbf{h}\_{v}, \pi\_{v}), \\ \mathcal{Z}\_{0}^{N} &= \mathcal{C}^{V\_{\Gamma}} \sum\_{\pi\_{\nu} \in \mathcal{S}\_{N}} \int \prod\_{v} \mathbf{dh}\_{v} \operatorname{Tr} \left[ \bigotimes\_{\ell} \rho\_{\ell}^{N} \bigotimes\_{v} \mathbb{P}\_{\mathbf{h}\_{v}} (\pi\_{v}) \right] \\ &\equiv \quad \mathcal{C}^{V\_{\Gamma}} \sum\_{\pi\_{\nu} \in \mathcal{S}\_{N}} \int \prod\_{v} \mathbf{dh}\_{v} \bigvee\_{\ell} (\mathbf{h}\_{v}, \pi\_{v}), \end{split} \tag{57}$$

respectively corresponding to summations of Feynman graphs N*A*(**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*) and N0(**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*) labelled by permutation operators <sup>P</sup>**h***v* (*<sup>π</sup>v*), at each node *v*, contracted with the *ρ<sup>N</sup>* densities at each link (see Figure 2). The difference between the reduced and full (density matrix) networks is encoded in the boundary condition, as *Z N* is defined with <sup>P</sup>(*π*<sup>0</sup> *A*; *N*, *d*) on *A* of *∂*Γ and <sup>P</sup>(<sup>1</sup>; *N*, *d*) on *B* of *∂*Γ, while *Z N* 0 is defined with <sup>P</sup>(<sup>1</sup>; *N*, *d*) for all boundary region *∂*Γ.

**Figure 2.** Free Group Field Theory propagators acting among two sets of N replicas for each node of the given graph.

It is easy to see at this stage what the impact of the i.i.d. assumption is in reducing the complexity of the propagator.<sup>7</sup> In this case, the permutation operator P acts only among the N replicas of the same node, independently taken across the graph. Clearly, / *v* <sup>E</sup>(*ρ<sup>N</sup> v* ) = E( /*<sup>V</sup> v ρ<sup>N</sup> v* ) as the permutation group S2*N* is much smaller that S2*NV*, reducing to (2*N*!)*<sup>V</sup>* the number of permutation patterns with respect to the (2*NV*)! allowed patterns of the indistinguishable case. Such a strong truncation is nevertheless consistent with a restriction to the tree-like, or "melonic", sector of the Feynmann diagrams of the 2*NV*-point function, which we expect to provide the leading order contribution to the divergence degree.

From a qualitative point of view, a truly general result in GFT would require to consider a fully interacting GFT and to drop the i.i.d. assumption. Moreover, it would require to leave the *tensor*

$$\mathbb{E}\_{0}\left[\bigotimes\_{v} (\left|\phi\_{\mathbb{P}}\right>\left<\phi\_{\mathbb{P}}\right|)^{\otimes N}\right] = \mathcal{C} \sum\_{\pi \in S\_{N \mid \mathbb{P}}} \mathbb{P}(\pi)$$

$$= \sum\_{\pi \in S\_{N \mid \mathbb{P}}} \prod\_{v=1}^{N} \prod\_{n=1}^{V\_{\mathbb{P}}} \int \mathrm{d}h\_{\pi \mathsf{d}} \mathbf{dg}\_{\mathsf{p} \mathsf{d}} \mathrm{d}\mathsf{g}'\_{\pi(n\mathsf{x})} \,\delta\left(h\_{\pi \mathsf{d}} \mathbf{g}\_{\pi \mathsf{d}} \mathbf{g}'\_{\pi(n\mathsf{x})}\right)^{-1} |\left|\mathbf{g}\_{\mathsf{p} \mathsf{a}}\right\rangle\left<\mathbf{g}'\_{\pi(n\mathsf{x})}\right|,\tag{58}$$

<sup>7</sup> Note that, if we keep the indistinguishability condition, the calculation changes

where <sup>S</sup>*NV*Γ is the permutation group of *NV*Γ objects, which corresponds to the permutations of *NV*Γ nodes; *hna* comes from the required gauge symmetry of the propagator; C is a constant factor.

*network* representation of the GFT multi-particle operator, and derive the entropy for the generic operator given in (15).

### *4.2. Mapping to BF Theory Partition Function*

The partition functions *ZN* and *<sup>Z</sup>N*0 correspond to summations of two *auxiliary* networks N*A*(**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*) and N0(**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*). We shall proceed by giving the main ingredients of the derivation, while referring the reader to the original literature (see e.g., [18]) for a fully detailed description of the combinatorics of the calculation.

First of all, in opening the expressions for *ZN* and *<sup>Z</sup>N*0 , we shall notice that the action of <sup>P</sup>**h***v* (*<sup>π</sup>v*) at each node is decoupled among the incident legs. Due to the strandwise action of the propagator, the value of the networks N*A*(**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*) and N0(**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*) can be written as factorised products over internal (e.g., trivial propagators acting in the picture below) and boundary links

$$\mathcal{N}\_{A}(\mathbf{h}\_{\upsilon}, \boldsymbol{\pi}\_{\upsilon}) = \prod\_{\boldsymbol{\varepsilon} \in \Gamma} \mathcal{L}\_{\boldsymbol{\varepsilon}}(\boldsymbol{\pi}\_{\upsilon}, \boldsymbol{\pi}\_{\upsilon'}; \mathbf{h}\_{\upsilon}, \mathbf{h}\_{\upsilon'}) \prod\_{\boldsymbol{\varepsilon} \in A} \mathcal{L}\_{\boldsymbol{\varepsilon}}(\boldsymbol{\pi}\_{\upsilon}, \boldsymbol{\pi}\_{A}^{0}; \mathbf{h}\_{\upsilon}) \prod\_{\boldsymbol{\varepsilon} \in B} \mathcal{L}\_{\boldsymbol{\varepsilon}}(\boldsymbol{\pi}\_{\upsilon}, \mathbb{1}; \mathbf{h}\_{\upsilon}), \tag{59}$$

$$\mathcal{N}\_0(\mathbf{h}\_{\upsilon\prime}, \boldsymbol{\pi}\_{\upsilon}) = \prod\_{c \in \Gamma} \mathcal{L}\_{\varepsilon}(\boldsymbol{\pi}\_{\upsilon\prime}, \boldsymbol{\pi}\_{\upsilon\prime}; \mathbf{h}\_{\upsilon\prime}, \mathbf{h}\_{\upsilon\prime}) \prod\_{c \in \partial\Gamma} \mathcal{L}\_{\varepsilon}(\boldsymbol{\pi}\_{\upsilon\prime}, \mathbb{1}; \mathbf{h}\_{\upsilon}). \tag{60}$$

On the internal (dotted) links, L(*<sup>π</sup>*, *<sup>π</sup>*, **h**, **h**) can be written as a trace of a modified representation of a permutation group element ≡ (*π*)−1*π*, such that

$$\mathcal{L}(\boldsymbol{\pi}, \boldsymbol{\pi}'; \mathbf{h}, \mathbf{h}') = \text{Tr}\left[\mathbf{P\_h}(\boldsymbol{\pi}) \rho\_\varepsilon^N \mathbf{P\_{h'}}(\boldsymbol{\pi})\right] = \text{Tr}\left[\mathbf{P\_h}\left((\boldsymbol{\pi}')^{-1}\boldsymbol{\pi}\right)\right] \equiv \text{Tr}\left[\mathbf{P\_h}\left(\boldsymbol{\mathcal{O}}\right)\right],\tag{61}$$

where

$$\mathbf{h} = \left\{ H\_{\mathfrak{a}} \mid H\_{\mathfrak{a}} \equiv \left( h'\_{\mathfrak{a}\circ(a)} \right)^{\dagger} h\_{\mathfrak{a}\prime} \,\,\forall a = 1, \dots, N \right\}.$$

In particular, any element ∈ S*N* can be expressed as the product of disjoint cycles C*i*

$$
\varpi \equiv \prod\_{i}^{\chi(\mathcal{O})} \mathcal{C}\_{i\nu} \tag{62}
$$

leading to a specific simple form for the individual link contributions

$$\mathcal{L}(\pi, \pi'; \mathbf{g}, \mathbf{g}') = \quad \text{Tr}\left[\mathbf{P}\_{\mathbf{g}}\left(\mathcal{O}\right)\right] = \prod\_{i} \text{Tr}\left[\mathbf{P}\_{\mathbf{g}}\left(\mathcal{G}\_{i}\right)\right] \tag{63}$$

$$= \quad \prod\_{i} \int \prod\_{k=1}^{r\_{i}} \text{dg}\_{a\_{k}^{i}} \delta\left(H\_{a\_{k}^{i}} \mathcal{g}\_{a\_{[k]\_{k}}^{i}}^{\dagger}\right) = \prod\_{i}^{\chi(\mathcal{O})} \delta\left(\prod\_{k=1}^{\ell\_{r\_{i}}} H\_{a\_{k}^{i}}\right) \tag{64}$$

which is expressed as the product of the traces of the individual cycles C*i*. Indeed, one can eventually realise [18,20] that the integral of the pattern networks N (**<sup>h</sup>***<sup>v</sup>*, *<sup>π</sup>v*) on the gauge holonomies are equivalent to the amplitudes of a three-dimensional topological BF field theory [50], with given boundary condition. Such amplitudes are discretized on a specific 2-complex comprised by the *N* replicas of the networks, with each different pattern P corresponding to a different 2-complex. The simple form of the various functions entering the calculation of the entropy directly follows from the specific approximations used in the calculation of expectation values, namely the choice to neglect the GFT interactions and to consider tensor-network like observables with simple delta functions associated with the links' kernels.

In the local averaging setting, the expectation values of the partition functions therefore reduce to contractions of single node 2 *N*-point functions of the GFT model, expanded in series of Feynman amplitudes and corresponding to a local BF gravity spin-foams (see Figure 3). In particular, if one specifies to a 3-valent GFT field theory, with simplicial four-field interaction kernels, each domain amplitude will correspond to a Boulatov (or 3-*d* BF theory) spinfoam amplitude, whose semi-classical limit coincides with Ponzano–Regge gravity [50].

**Figure 3.** (Above) Simple open graph associated with a boundary tensor state (and its conjugate). (Down) Cartoon picture of the expectation value of the partition function defined by the trace of the boundary tensor state density matrix. The *N* replicas of the nodes of the graph define domain regions D, connected among them via (replicas of) the graph links. The boundary links are contracted among the conjugate graphs via a global trace. The expectation value generates a set of independent spin-foam channels N among conjugate domains (GFT 2 *N*-point functions), due to the i.i.d. assumption considered. Indistinguishability would instead merge the different spin-foams into a single one, raising the degree of correlations of the boundary degrees of freedom. Nevertheless, the local approximation seems to convey the leading order contribution of the quantum entanglement of the bipartite system.

=

### *4.3. Entanglement Scaling and Divergence Degree*

Given the expression for the Rényi entropy, *e*(<sup>1</sup>−*<sup>N</sup>*)*SN* = *ZN*/*ZN*0 , finding the scaling of the entanglement entropy amounts to identifying the most divergent terms of the partition functions *ZN* and *<sup>Z</sup>N*0 , which corresponds to the divergence degree of the BF theory amplitudes discretized on a lattice [58–62]. The recipe given in [18] consists of a combination of coarse graining and combinatorics.

Once given global boundary conditions for N0, *π* = 1 and *h* = 1, one can coarse-grain the boundary of N0 into a single node with *π* = 1 and *h* = 1. Accordingly, the boundary of N*A* gets coarse-grained into two nodes, one of which corresponds to *A* with *π* = C0, *h* = 1 and the other to *B* with *π* = 1 and *h* = 1. The corresponding closed graphs are denoted as Γ0 and Γ*AB*. Once we take the expectation value on the *N* replicas at each node, then a given pattern <sup>P</sup>(*<sup>π</sup>v*) divides Γ0 and Γ*AB* into regions (set of nodes), each one coloured with permutation group *πm* and *N* integrals over *hm* (see Figure 4 below).

**Figure 4.** Example of a permutation pattern <sup>P</sup>(*<sup>π</sup>v*) dividing Γ*AB* into regions (set of nodes), each one colored with permutation group element *π* and *N* integrals over holonomies **h**.

In full analogy with the *SN* Ising model [15], links connecting different regions identify boundaries which can be interpreted as domain walls. We are interested in finding the scaling behaviour of the Rényi entropy in the large bond regime *D*(Λ) 1, classically corresponding to the long-range ordered phase for such an Ising-like model, where the entropy of a boundary region is known to be directly related to the energy of the domain wall between domains of the order parameter [15].

Different local regions *R*, with uniform boundary conditions, can be further coarse-grained via gauge invariance to single block nodes. For a given region, the degree of divergence is counted by the number of links in the region, minus the number of trivialised links on a maximal spanning tree *TR*

$$\#\_{R} = \left(\#\_{c \in R} - \#\_{T\_{R}}\right) N. \tag{64}$$

One can then generally show that the pattern with no domain walls, Γ0, corresponding to the ordered phase in which all nodes are assigned with the same permutation group, has the highest degree of divergence

$$\#\_0 = (\#\_{\mathfrak{e}\in\Gamma\_0} - \#\_{T\_{\Gamma\_0}})N,\tag{65}$$

where #*e*∈Γ0is the number of links in the Γ0 graph.

For the bipartite graph Γ*AB*, defined by the assignment of different boundary conditions for *A* and *B*, patterns with a single domain wall have higher divergence degree than the multi-domain walls configurations. The coarse-grained graph contains only two block nodes and one finds [18]

$$\#\_{\Gamma\_{AB}(\pi\_m)} \le \quad \#\_{AB} = \#\_0 + (1 - N) \min(\#\_{\ell \in \tilde{\mathcal{O}}\_{AB}}).\tag{66}$$

*Universe* **2019**, *5*, 211

> The leading order divergence terms of the amplitudes are given by

$$\begin{array}{rcl} Z\_0^N &=& \mathbb{C}^{V\_\Gamma}[D(\Lambda)]^{\mathfrak{H}\_0}, \\ Z\_N &=& \mathbb{C}^{V\_\Gamma}[D(\Lambda)]^{\mathfrak{H}\_0 + (1 - N)\min(\mathfrak{H}\_{\mathfrak{e} \in \mathfrak{d}\_{AB}})} \left[1 + \mathcal{O}(D(\Lambda)^{-1})\right] \end{array}$$

and the *N*th order Rényi entropy *SN* eventually reads

$$\mathrm{e}^{(1-N)\mathbb{S}\_N} = \frac{Z\_N}{Z\_0^N} = \left[ D(\Lambda) \right]^{(1-N)\min(\mathfrak{h}\_{\varepsilon \ominus \partial\_{AB}})} \left[ 1 + \mathcal{O}(D(\Lambda)^{-1}) \right]. \tag{67}$$

When *N* goes to 1, the leading term of the entanglement entropy *S*EE is given by

$$S\_{\rm EE} = \min(\#\_{c \in \hat{\theta}\_{AB}}) \ln D(\Lambda),\tag{68}$$

,

where *D*(Λ) with Λ 1 reads as a regularisation of each BF bubble divergence *<sup>δ</sup>*(1).

The result, within a different formal setting, reproduces the universal behaviour typical of random tensor network states [15]. However, differently from the case of a standard random tensor network, the considered GFT states carry an inherent geometric characterisation. The graph Γ is dual to a 2D simplicial complex. Each node is dual to a triangle and each link is dual to an edge of this complex, and the specific GFT model endows the simplicial complex with dynamical geometric data. In particular, the proportionality of the entropy to the cardinality of the minimal domain wall *σ*min ≡ min(#*e*∈*∂AB* ) has a clear geometric interpretation, which becomes apparent in passing from a group element to a spin representation description of the dynamical fields, in the sense of discrete geometry. In this case, indeed, we have

$$\text{Area}(\sigma\_{\text{min}}) = \sum\_{\mathbf{c} \in \sigma\_{\text{min}}} \ell\_{\mathbf{c}}(j\_{\mathbf{c}}) = \langle \ell\_{j\_{\mathbf{c}}} \rangle |\sigma\_{\text{min}}|\_{\text{\textquotedblleft}} \tag{69}$$

where the length of each edge *j*, in any given eigenstate of the length operator [26], is a function of the irreducible representation *je* associated with it, and to the dual link. Therefore, the cardinality of the minimal global domain wall can be interpreted as the 1*D*-area (length) of a dual discrete minimal one-dimensional path and we can write |*<sup>σ</sup>*min<sup>|</sup> = Area(*<sup>σ</sup>*min)/ *je* . [20].

Equation (68) can then be understood as the discrete tensor network analogue of the Ryu–Takayanagi formula in the context of group field theory [21], if we consider the path integral averaging over the open network Γ as a simplified model of a bulk/boundary (spinfoam/network state) duality [18].

### **5. Holographic Scaling for Interacting with GFT**

A complete geometric characterisation of the result requires working with a fully interacting GFT model. Indeed, it is the simplicial character of the interaction kernels that actually provides a connection between GFT and quantum gravity spin-foams. In this sense, both consistency and robustness of the result in (68) require exploring the possible modifications to the RT formula induced by group field interactions.

At the perturbative level, this amounts to calculating the first order correction of E[*SN*(*A*)] in the GFT coupling constant *λ*, namely <sup>E</sup>[Tr(*ρ*<sup>⊗</sup>*N*P*A*)] and <sup>E</sup>[Tr(*ρ*<sup>⊗</sup>*<sup>N</sup>*)], with

$$\mathbb{E}[\text{Tr}(\rho^{\ominus N}\mathbb{P}\_A)] = \quad \mathbb{E}\_0[\text{Tr}(\rho^{\ominus N}\mathbb{P}\_A)] + \lambda\mathbb{E}\_0\left[\mathbb{S}\_{\text{int}}[\phi, \overline{\phi}]\text{Tr}(\rho^{\ominus N}\mathbb{P}\_A)\right],\tag{70}$$

$$\mathbb{E}[\text{Tr}(\rho^{\ominus N})] \quad = \quad \mathbb{E}\_0[\text{Tr}(\rho^{\ominus N})] + \lambda\mathbb{E}\_0\left[\mathbb{S}\_{\text{int}}[\phi, \overline{\phi}]\text{Tr}(\rho^{\ominus N})\right].$$

This is a combinatorially highly non-trivial problem, as the interaction processes correspond to further stranded diagrams that contribute to the expectation value of *Z*(*N*) *A*/0(see Figure 5). The behaviour of the amplitudes for rank-3 fields and Boulatov (simplicial) four fields interaction kernel <sup>V</sup>sym(**g**(1)**g**(2)**g**(1)**g**(2))

$$\int \prod\_{\ell=1}^{4} d\mathfrak{h}\_{\ell} \, \delta(\mathfrak{h}\_1 \mathfrak{g}\_1^{(1)}, \mathfrak{h}\_3 \underline{\mathfrak{g}}\_1^{(1)}) \delta(\mathfrak{h}\_1 \mathfrak{g}\_2^{(1)}, \mathfrak{h}\_4 \underline{\mathfrak{g}}\_2^{(2)}) \delta(\mathfrak{h}\_1 \mathfrak{g}\_3^{(1)}, \mathfrak{h}\_2 \mathfrak{g}\_3^{(2)})$$
 
$$\delta(\mathfrak{h}\_2 \mathfrak{g}\_1^{(2)}, \mathfrak{h}\_4 \underline{\mathfrak{g}}\_1^{(2)}) \delta(\mathfrak{h}\_2 \mathfrak{g}\_2^{(2)}, \mathfrak{h}\_3 \underline{\mathfrak{g}}\_2^{(1)}) \delta(\mathfrak{h}\_3 \underline{\mathfrak{g}}\_3^{(1)}, \mathfrak{h}\_4 \underline{\mathfrak{g}}\_3^{(2)})$$

is studied in [20]. Therein, the authors provide a series of theorems aimed at constraining the complexity of the combinatorial pattern. In particular, it was shown that patterns with a single interaction happening between two incoming and two outgoing fields of the same network node *v* leads at most to the same number of divergences as in a maximal case of the free theory, where the interaction at *v* is replaced by a free propagation [20].

**Figure 5.** Examples of insertions of Boulatov 4-fields interaction vertices V in the amplitude among replicas of the individual node of the graph. Multiple of interaction vertices generate new faces in the amplitudes.

Interestingly, the situation changes if one drops the assumption of gauge invariance (closure) for the GFT field. In this case, the expected *N*th Rényi entanglement entropy is estimated to be

$$\mathbb{E}\_{\text{n-sym}}[S\_N(\rho\_A)] \approx [\ln D(\Lambda) - \lambda N] \min(\#\_{e \in \hat{\mathcal{O}}\_{AB}}).\tag{71}$$

The linear order corrections do modify the asymptotical scaling of the Rènyi entanglement entropy with the area of a minimal surface. In particular, the results in [20] show that the proportionality factor is corrected by an additive term linear in the perturbed group field theory coupling constant *λ*. However, no additive leading order correction to the area scaling entropy formula emerges from the analysis. A more systematic analysis of such dynamical regularisation and its relation with the emergence of an effective gravitational coupling is an open interesting issue for future work.

### **6. Discussion and Conclusions**

Group field theories realise a kinematical description of quantum space geometry in terms of discrete, pre-geometric degrees of freedom of combinatorial and algebraic nature, described in terms of spin-network states. In loop quantum gravity (LQG), such theories play the role of *auxiliary* field theories, whose partition functions, for appropriate choices of the kernels, provide

generating functionals for the LQG spinfoams: a covariant path integral realisation of spacetime as a transition amplitude between boundary spin network states. More generally, GFTs provide a versatile field-theoretic tool to study the very emergence of space-time quantum geometry via path integral techniques and a quantum many-body approach associated with their second quantised formalism [57,63]. Moreover, as higher order generalisations of matrix models, group field theories can be put in direct formal relation with tensor network algorithms [10–12,17,18] in condensed matter theory.

In this review, we considered a setting *similar* to that of a gauge field theory on a lattice, which in the background independent quantum gravity context consists of nets of fundamental quanta of space that admit an interpretation as quantized fuzzy geometries.

The entanglement structure of the wave-function for such collections of quanta is encoded into a ∞-rank tensor of coefficients, which we interpret here as an *open* GFT network state. In [18–20], such an open GFT spin network is interpreted as a *tensor network representation* of some quantum geometry wave-function, with physical indices corresponding to boundary degrees of freedom of the "auxiliary" GFT network state, realised by a spin network.

In this setting, the degree of entanglement of a generic quantum region of space can be measured *holographically*, in terms of the entanglement entropy of the bipartite auxiliary spin network. Such entanglement is directly related to the topology of the internal network, which, differently from many similar derivations in quantum gravity, is only partially fixed by the choice of the auxiliary network architecture. Indeed, due to the random character of the nodes, expectation values imply a dynamical characterisation of the internal graph topology, which end up being directly specified by the choice of the GFT model. Therefore, different GFT models a priori induce different architecture for the internal network, with higher order interactions corresponding to a higher degree of connectivity of the graph.

The GFT formalism, along the lines first proposed by [15], allows for an explicit computation of the entanglement entropy, which in [18–20] is realised within a set of simplifying assumptions:

1. indistinguishability and *i.i.d.*: the symmetry of the group field naturally selects a bosonic quantum many-body characterisation for the quanta comprising our network states. This translates into indistinguishability for the nodes comprising the network. Furthermore, along with the results in [15], in the statistical evaluation of the averaged Rényi entropy, 'tensor' fields at the nodes are individually independently distributed. This leads to a sensible simplification of the auxiliary bulk averages, reduced to a set of *local* amplitudes defined among conjugate domains of replicas of the fields at each node. Such i.i.d. assumption is a natural choice in [15], where random tensor networks are associated with maximally mixed states. In our setting, accordingly, the *i.i.d.* assumption presupposes actual *absence of dynamical interaction* among the quanta sitting at the nodes of the graph. Effectively, quanta couples through the network only via local adjacent link entanglement. Such a non-interacting quantum many-body assumption is very strong. Nevertheless, it is consistent with the simplifying choice of dealing with *product* coherent states in the adopted 2nd quantised formalism. On the other hand, the i.i.d. assumption embodies some a priori knowledge of the position of the individual nodes (their replicas and conjugates) on the graph, hence somehow violating indistinguishability.


As a result, the set of assumptions above induces a structure of correlations that is essentially local in its *leading* contribution to the entropy, even when the interaction among fields at different nodes is considered [18,20]. This suggests that the emergence of the dominant holographic behaviour for the entropy is intertwined with the universality features of the large bond regime. In the presence of (weak) interactions, such area scaling behaviour is shown to remain a solid feature of the *quantum typical regime* [20]. 8

Many aspects can be tuned to further test such a holographic behaviour in the specific aproach. At the level of the free theory, it would be interesting to look at changes in the area scaling due to minimal modification of the very free propagator, via more general heat-kernel techniques, possibly expressed in terms of a different choice of mutual information among the half link states [14]. For higher order perturbations, where the diagrams are dominated by the bulk structure, induced by the creation of new interaction vertices [20], we may expect significant changes in the area scaling behaviour, even if the local averaging scheme induced by the i.i.d. assumption is preserved. Outside the typicality regime, it is natural to expect that fluctuations in the annealed average (43) become more and more relevant and the area scaling behaviour is eventually lost. A detailed study in this sense requires strong computational efforts, necessary to explore the properties and combinatorics of the dynamically induced bulk auxiliary networks.

More generally, concerning the potential of the GFT formalism in the study of holographic entanglement in quantum geometry, much still remain to be understood. We showed how the formalism of *second quantization* for GFT, and the use of a coherent state basis, allow for a new approach to the study of entanglement in quantum geometry. Standard quantum mechanical correlations among links and nodes of the quantum many-body network state are replaced by correlations among *excitations* of modes of the tensor field representing the particle. This is likely the right formalism for understanding the diverse roles played by entanglement in quantum gravity, from the local entanglement among single quanta responsible for the connectivity of a quantum geometry configuration, to the local and non-local entanglement among modes of the GFT field, intended as elementary excitations induced by the *collective* behaviour of the underlying many-body quantum system. This would effectively allow for maximally disregarding the auxiliary spin network structures and work in full indistinguishability with bosonic GFT quanta of geometry in the Fock space.

<sup>8</sup> A possibly deep relation between such universal regime and the universality proper of tensor models in the large *D* regime (see e.g., [67]) has been pointed out and discussed in [18].

However, harnessing entanglement in the 2nd quantisation formalism is still a fundamental open problem already in condensed matter and quantum information theory [68–70]. In our derivation, the second quantisation scheme was in fact only partially realised. Indistinguishability would have actually prevented us from defining an AB-factorisation of the Hilbert space. In this regard, we proceeded by selecting a single graph Hilbert space within the full Fock space, hence coming back to a first quantisation formalism. Indistinguishability was further violated by the local averaging approach, in Section 4.1, induced by the i.i.d. assumption. A GFT derivation that makes full use of the mode aspect of second quantization is therefore substantially left for future work. Success in this sense would be evocative of (and consistent with) the general perspective that sees continuum spacetime and geometry as emergen<sup>t</sup> from the collective, quantum many-body description of the fundamental GFT degrees of freedom, the same perspective motivating a large part of the literature and in particular the one concerned with GFT renormalization (both perturbative and non-perturbative) [71–82].

Most recent work in this direction focussed on entanglement among GFT modes induced by interaction, starting from a generalization of the Bogoliubov description of a weakly interacting Bose gas to the GFT framework [83]. Analogous insights on the inherently dynamic character of entanglement among GFT modes in second quantization appear in the series of recent results proposing a general procedure for constructing states that describe macroscopic, spatially homogeneous universes as Bose–Einstein condensates (see e.g., [53]). Strong results in this sense would provide a new quantitative tool to unravel the relation between entanglement, holography and emergence of spacetime geometry in quantum gravity.

For instance, an interesting direct correspondence between black holes and Bose–Einstein condensates (BECs) of gravitons at the point of maximal packing was recently proposed in [84], though within a standard *effective field theory framework*. In this case, the physics of maximally-packed gravitational systems is identified with the general behaviour of BECs at the critical point of the quantum phase transition, while collective nearly gapless excitations of the quantum condensate are shown to define the holographic degrees of freedom responsible for the known semiclassical holographic properties of black holes (BH) [84]. In this exemplary case, GFTs would provide a unique *non-perturbative* formalism to investigate the foundations of black hole entropy [85] and holography as a general quantum phenomenon of nature.

### **Funding:** This research received no external funding.

**Acknowledgments:** It is a pleasure to acknowledge Daniele Oriti, Mingyi Zhang and Alex Goeßmann for their precious contribution in the exploration of such a fascinating topic. I am particularly grateful to Hermann Nicolai and to the Albert Einstein Institute in Potsdam, for the support received during the preparation of the manuscript.

**Conflicts of Interest:** The author declare no conflict of interest.

### **Appendix A. Coherent States Over-Completeness**

Consider two coherent state |*ϕ* and |*ϕ* . One can prove that they are not orthogonal

$$
\langle \boldsymbol{\varrho} | \boldsymbol{\varrho}' \rangle = \frac{1}{\mathcal{N}\_{\boldsymbol{\varrho}'} \mathcal{N}\_{\boldsymbol{\varrho}'}} \exp \left[ \int \mathbf{dg} \,\overline{\boldsymbol{\varrho}(\mathbf{g})} \,\boldsymbol{\varrho}'(\mathbf{g}) \right]. \tag{A1}
$$

The state |*ϕ* can be decomposed by the *n*-particle state basis |··· , **g***<sup>a</sup>*, ··· as

$$<\langle \boldsymbol{\varrho} \rangle = \left| 0 \right\rangle \left\langle 0 \middle| \boldsymbol{\varrho} \right\rangle + \sum\_{n=1}^{\infty} \int \prod\_{a=1}^{n} \mathbf{dg}\_{a} \left| \mathbf{g}\_{1}, \dots, \mathbf{g}\_{n} \right\rangle \left\langle \mathbf{g}\_{1}, \dots, \mathbf{g}\_{n} \middle| \boldsymbol{\varrho} \right\rangle,\tag{A2}$$

where **g**1, ··· , **<sup>g</sup>***n*|*ϕ* is given as

$$
\begin{split}
\langle\langle\mathbf{g}\_{1},\cdots,\mathbf{g}\_{n}|\boldsymbol{\varrho}\rangle\ &=&\frac{1}{\mathcal{N}\_{\boldsymbol{\varrho}}}\frac{1}{\sqrt{n!}}\prod\_{a}^{n}\left[\int\!\!\!\mathbf{d}\mathbf{g}\_{a}^{\prime}\,\boldsymbol{\varrho}(\mathbf{g}\_{a}^{\prime})\right]\langle\mathbf{g}\_{1},\cdots,\mathbf{g}\_{n}|\mathbf{g}\_{1}^{\prime},\cdots,\mathbf{g}\_{n}^{\prime}\rangle\\&=&\frac{1}{\mathcal{N}\_{\boldsymbol{\varrho}}}\frac{1}{(n!)^{3/2}}\sum\_{\pi\in\mathcal{S}\_{\boldsymbol{a}}}\prod\_{a}^{n}\left[\int\!\!\!\mathbf{d}\mathbf{g}\_{a}^{\prime}\,\mathrm{d}\boldsymbol{h}\_{a}\,\boldsymbol{\varrho}(\mathbf{g}\_{a}^{\prime})\delta\left(\mathbf{g}\_{\pi(a)}\boldsymbol{h}\_{\boldsymbol{a}}(\mathbf{g}\_{a}^{\prime})^{\dagger}\right)\right]\\&=&\frac{1}{\mathcal{N}\_{\boldsymbol{\varrho}}}\frac{1}{\sqrt{n!}}\prod\_{a}^{n}\,\boldsymbol{\varrho}(\mathbf{g}\_{a}).
\end{split}
\tag{A3}
$$

In order to obtain the resolution of identity in terms of |*ϕ* , let us first introduce the gauge fixed field *ϕ*GF

$$\varphi\_{\rm GF}(\mathbf{g}) \equiv \varphi(\mathbf{g}h\_1^{-1}) = \varphi(1, \, h\_2h\_1^{-1}, \, h\_3h\_1^{-1}, \, h\_4h\_1^{-1}) \equiv \varphi([\mathbf{g}]).\tag{A4}$$

Then, we have the identities


where *K* is a parameter which we assume to be real, and *C* is a constant number

$$\mathcal{C} \equiv \int \mathcal{D}q \mathcal{D}\overline{\mathbf{p}} \,\delta\_{\mathcal{C}}[\varphi - \varphi\_{\rm GF}] \, \exp\left[ (-K - 1) \int \mathbf{dg} \,\overline{\varphi(\mathbf{g})} \,\varphi(\mathbf{g}) \right]. \tag{A6}$$

For the third equality, we use the Wick's theorem, and, for the fourth equality, we reintroduce the gauge symmetry such that the arguments of *ϕ* are equal in weight.

The resolution of identity in the Fock space F is in terms of |*ϕ* 

$$\mathbb{1}\_{\mathcal{F}} = \mathbb{C}^{-1} \int \mathcal{D}\boldsymbol{\varrho} \mathcal{D}\overline{\boldsymbol{\varrho}} \, \delta\_{\mathbb{C}} [\boldsymbol{\varrho} - \boldsymbol{\varrho}\_{\operatorname{GF}}] \, \mathrm{e}^{-K \int \mathrm{d} \mathbf{g} \, \overline{\boldsymbol{\varrho}(\mathbf{g})} \, \boldsymbol{\varrho}(\mathbf{g})} \, |\boldsymbol{\varrho}\rangle \, \langle \boldsymbol{\varrho}| \,. \tag{A7}$$
