In the long-standing quest to reconcile gravity with quantum mechanics, profound connections have been unveiled between concepts traditionally pertaining to a quantum information theory, such as entanglement, and constitutive features of gravity, like holography. Developing and promoting these connections from the conceptual to the operational level unlock access to a powerful set of tools which can be pivotal toward the formulation of a consistent theory of quantum gravity. Here, we review recent progress on the role and applications of quantum informational methods, in particular tensor networks, for quantum gravity models. We focus on spin network states dual to finite regions of space, represented as entanglement graphs in the group field theory approach to quantum gravity, and illustrate how techniques from random tensor networks can be exploited to investigate their holographic properties. In particular, spin network states can be interpreted as maps from bulk to boundary, whose holographic behavior increases with the inhomogeneity of their geometric data (up to becoming proper quantum channels). The entanglement entropy of boundary states, which are obtained by feeding such maps with suitable bulk states, is then proved to follow a bulk area law with corrections due to the entanglement of the bulk state. We further review how exceeding a certain threshold of bulk entanglement leads to the emergence of a black hole-like region, revealing intriguing perspectives for quantum cosmology.

## I. INTRODUCTION

Holography has been a driving theme of research in quantum gravity since the discovery of the Bekenstein–Hawking area law for black hole entropy^{1,2} and the discussion on information loss and Hawking radiation.^{3,4} Aspects and realizations of the holographic principle, originally proposed by 't Hooft^{5} and later developed by Susskind^{6} and Bousso,^{7} have been extensively studied at both classical and quantum levels. Relevant instances include, out of a very wide range of contributions, early work on the microscopic interpretation of the black hole entropy;^{8–11} on the recovering of gravitational dynamics from the thermodynamics of boundaries;^{12,13} on the duality between the gravitational theory of asymptotically anti de Sitter (AdS) spacetime and a conformal field theory (CFT) leaving on its boundary, known as AdS/CFT correspondence^{14–17} (and within the latter, the Ryu–Takayanagi formula^{18,19} relates the boundary entanglement entropy to the area of a bulk surface); and on holography in loop quantum gravity.^{20–22}

In recent years, an intriguing connection among gravity, holography, and quantum entanglement has come to light. On the one hand, several results point to entanglement as the “glue” of spacetime;^{23–25} on the other hand, entanglement turns out to be intimately tied to holography in quantum many-body systems,^{26} and quantum spacetime can indeed be understood in several background-independent approaches to quantum gravity as a collection of (fundamental, “pre-geometric”) quantum entities,^{27} i.e., as a (background-independent) quantum many-body system.^{28} Understanding the origin of the gravity/holography/entanglement threefold connection would, therefore, be a major step toward the formulation of a theory of quantum gravity.^{29}

The main aim of this article is to review recent results^{30–32} that stand out for investigating holography directly at the level of *quantum gravity states* in a *quasi-local context* and via a *quantum information language*. The focus is on finite regions of 3D quantum space modeled by spin networks, i.e., graphs decorated by quantum geometric data (a formalism originally proposed by Penrose^{33}), which enter, as kinematical states, various background-independent approaches to quantum gravity.^{34–36} Crucially, such states are understood as arising from the entanglement of the quantum entities (“atoms of space”) composing the spacetime microstructure in the group field theory (GFT) framework,^{36,37} that is, as *graphs of entanglement*. This formalism has the remarkable property of realizing, directly at the level of the quantum microstructure of spacetime, the interrelation between entanglement and space connectivity supported by several results in quantum gravity contexts and beyond.^{18,19,23–25,38} Moreover, as entanglement graphs, the spin network states are put in correspondence with tensor networks,^{39} a quantum information language that efficiently encodes entanglement in quantum many-body systems. Such an information-theoretic perspective on spin network states is then exploited to investigate the role of entanglement (and quantum correlations more generally) in the holographic features of quantum spacetime via tensor network techniques.

As further reviewed in this article, the aforementioned approach is shared by a rich body of work at the interface of quantum gravity, quantum information, and condensed matter physics, which looked at entanglement on spin networks as a tool for probing and reconstructing geometry. It includes modeling of quantum black holes and the computation of the horizon entropy,^{40–43} the reconstruction of a notion of distance on spin networks from entanglement,^{44,45} the characterization of the entanglement entropy between an arbitrary region of a spin network and its complement,^{46–49} the use of entanglement to glue quantum polyhedra dual to spin network vertices,^{50,51} and the study of the holographic properties of spin network states.^{51–53}

The results of Refs. 30–32 reviewed here investigate holography in finite regions of quantum space from two different perspectives: (i) by studying the flow of information from the bulk to the boundary and (ii) by analyzing the information content of the boundary and its relationship with the bulk. The idea behind perspective (i) is the possibility (pointed out for the first time in Ref. 54) to interpret every spin network state as a bulk-to-boundary map, and the holographic character of the latter is traced back to how close it comes to being an isometry. The impact of combinatorial structure and geometric data of spin network states (matching random tensor networks) on the “isometry degree” of the corresponding bulk-to-boundary maps is then studied by relying on random tensor network methods. Perspective (ii) focuses on the entanglement entropy content of boundary states, obtained by feeding the aforementioned bulk-to-boundary maps with a bulk input state, upon varying the latter. The result is twofold: on the one hand, a bulk area law for the boundary entropy with corrections due to the bulk entanglement and on the other hand, the emergence of horizon-like surfaces when increasing the entanglement content of the bulk.

This focused review is structured in four sections. The first one is dedicated to the quantum gravity framework: Subsections II A and II B show the logical path from a quantized, elementary portion of space (a tetrahedron) to extended discrete quantum geometries and the dual spin network description; Subsection II C presents group field theories, quantum gravity models in which spin networks can be readily understood as graphs of entanglement, and as kinematic quantum gravity states; finally, Subsection II D illustrates the tensor network perspective on spin network states. In Sec. III, we give an overview of earlier results on the study of entanglement on spin network states and its role in reconstructing geometry. Section IV is dedicated to random tensor network techniques adapted to the considered quantum gravity framework; more specifically, it shows how to compute the Rényi-2 entropy of a certain class of spin network states via a statistical model. Section V contains the aforementioned results on the holographic features of spin network states matching random tensor networks from the perspective of bulk-to-boundary maps (Subsection V A) and of the entanglement entropy of boundary states (Subsection V B).

## II. QUANTUM GRAVITY STATES AS ENTANGLEMENT GRAPHS

Several approaches to quantum gravity, e.g., loop quantum gravity,^{34} spinfoam models,^{35} and group field theories^{36,37} (GFT), describe regions of 3D space via *spin networks*, graphs decorated by quantum geometric data. We review how spin networks can be constructed from elementary portions of space (e.g., small tetrahedra) quantized and glued together to form extended (discrete) spatial geometries; crucially, the gluing derives from entanglement, and spin networks can, thus, be regarded as the entanglement structure of many-body states for the set of elementary tetrahedra. We then introduce group field theories,^{36,37} quantum gravity models where the above picture is realized, and spin networks from many-body entanglement can be understood as kinematical quantum gravity states. We conclude by reviewing recent results^{30} on the formal correspondence between spin network states and tensor networks.

### A. Quantum tetrahedron and the dual spin network vertex

Consider an elementary portion of 3D space, a tetrahedron, whose faces are labeled by an index $i=1,2,3,4$. The (classical) geometry of the tetrahedron can be described by four vectors ${L\u2192i}i=14$, with $L\u2192i$ normal to the *i*th face and having length equal to the face area, which satisfy the *closure constraint*^{55}

The equivalence class of the four vectors ${L\u2192i}i=14$ under global rotations encodes a geometrical configuration of the tetrahedron. Note that, as the vectors ${L\u2192i}i=14$ are elements of the *su*(2) Lie algebra, we can equivalently describe the geometry of the tetrahedron via the dual *SU*(2) group elements ${gi}i=14$ [more precisely, via the equivalence class of ${gi}i=14$ under global *SU*(2) action].

The quantization of the phase space of geometries of a tetrahedron^{56–58} leads to the Hilbert space $h=L2(G4/G)$, where $G=SU(2)$; i.e., the quantum state of geometry of a tetrahedron is described by a wave-function $f(g\u2192)$, where $g\u2192={g1,g2,g3,g4}$, that satisfies

with $hg\u2192:={hg1,hg2,hg3,hg4}$.

By the Peter–Weyl theorem, the wave-function $f(g\u2192)$ can be decomposed into irreducible representations $j\u2208\mathbb{N}/2$ of *SU*(2),^{59}

where we used a vector notation for a set of variables attached to the four faces of the tetrahedron, e.g., $j\u2192={j1,j2,j3,j4}$; the *magnetic index m ^{i}* (

*n*) labels a basis of the

^{i}*j*-representation space $Vji$ (its dual $Vji*$); and $Dminiji(gi)$ is the Wigner matrix representing the group element

^{i}*g*. When taking into account the gauge symmetry [see Eq. (2)], both the expansion coefficients and the Wigner matrices end up contracted with a

^{i}*SU*(2)-invariant tensor, i.e., an

*intertwiner ι*, pertaining to the Hilbert space

and ensuring the gauge invariant recoupling of the four spins ${ji}i=14$. Equation (3) then becomes^{59}

where

is the generic element of the *spin network basis* and $fn\u2192\iota j\u2192:=fm\u2192n\u2192j\u2192\iota m\u2192$. We denote by $dj:=2j+1$ the dimension of the representation space *V ^{j}*, and by $Dj\u2192$ the dimension of the intertwiner space $ij\u2192$.

The spin network basis ${|j\u2192n\u2192\iota \u27e9}$ diagonalizes the area and volume operators^{60–63} and, thus, possesses a clear geometrical interpretation; more specifically, the *SU*(2) spin *j ^{i}* determines the area of the

*i*-face of the tetrahedron, while the intertwiner

*ι*determines its volume.

The quantum tetrahedron can be graphically represented as a vertex with four *edges*, each one identified by a *color i*, where the *i*th edge (denoted by *e ^{i}*) is dual to the

*i*th face of the tetrahedron and carries the corresponding quantum data (see Fig. 1): in the group basis, the edge

*e*carries a group variable

^{i}*g*; in the spin network basis, the edge

^{i}*e*carries a spin

^{i}*j*and, at the free endpoint, the magnetic index

^{i}*n*, while the intertwiner quantum number

^{i}*ι*is attached to the vertex itself. This structure is called the

*spin network vertex*.

At the level of the Hilbert space of the quantum tetrahedron, the spin network decomposition performed via the Peter–Weyl theorem reads

where the intertwiner space $ij\u2192$ is defined in Eq. (4).

The above construction can be easily generalized to any elementary polyhedron. In particular, the quantum version of a $(d\u22121)$-simplex [which is the simplest possible $(d\u22121)$-polytope] is dual to a *d*-valent vertex and described by the Hilbert space

In the following, we take into account this generalization and adopt, for the $j\u2192$-spin sector, the notation

Also, to clarify the role of the different degrees of freedom of spin network vertices, for some equations, we write the basis element $|j\u2192n\u2192\iota \u27e9$ of $hj\u2192$ in the form

i.e., as an explicit tensor product of the basis states of the intertwiner and representation spaces $|j\u2192\iota \u27e9\u2208ij\u2192$ and $|jini\u27e9\u2208Vji$, respectively.

### B. Gluing tetrahedra: Spin networks for 3D quantum geometries

A region of 3D space can be arbitrary well approximated by a collection of (suitably small) polyhedra adjacent to each other. As we are going to show, the quantum geometry of such a discrete space can be described by a set of interconnected spin network vertices corresponding to the single polyhedra;^{64} the result is a *spin network graph*,^{35} i.e., a graph *γ* dual to the space partition and decorated by quantum geometric data, as shown in Fig. 2.

More precisely, a spin network graph represents the quantum version of a *twisted geometry*.^{65–67} The latter is a collection of polyhedra in which adjacent faces possess the same area but have, in general, different shapes and/or orientations. That is, only “neighboring relations” are present in a twisted geometry: the planes of adjacent faces are not necessarily parallel. Twisted geometries, thus, differ from standard Regge triangulations, in which faces of neighboring polyhedra, having the same area, shape, and orientation, perfectly adhere to each other. In the following, we explain how spin networks can arise from the gluing of vertices dual to polyhedra, where by “gluing” we mean establishing an adjacency relationship between them as defined in twisted geometries.

Consider a set $v=1,\u2026,N$ of open spin network vertices of valence *d*, which is described by the Hilbert space $hN=L2(Gd\xd7N/GN)$. We illustrate the gluing of vertices with an example. Given two vertices *v* and *w*, we want to glue the *i*th edge of *v* (denoted by $evi$), which carries the group variable $gvi$ with the *j*th edge of *w* (denoted by $ewj$), which carries $gwj$. As both edges are outgoing, the resulting link from *v* to *w* (denoted by $\u2113vwij$) carries the group element $gvi(gwj)\u22121$. Once connected, the two vertices are, thus, invariant under the simultaneous right action of the group on the edges $evi$ and $ewj$, as $gvih(gwjh)\u22121=gvi(gwj)\u22121$ $\u2200h\u2208SU(2)$. Starting from the set of open vertices in the state $\psi \u2208hN$, such a symmetry (that is, the gluing of edges) can be implemented via the following group averaging:^{64,68}

which, in fact, causes the resulting $\psi \gamma $ to depend on $gvi$ and $gwj$ only through the product $gvi(gwj)\u22121$. The wave-function $\psi \gamma $ is then associated with a graph *γ* involving the link $\u2113vwij$. In the group basis, the geometric data attached to a spin network graph, thus, consist in a group element on every edge of the graph with gauge invariance at each vertex. This structure is, therefore, described by the Hilbert space $h\gamma =L2(GL/GN)$, where *L* is the number of links of *γ*.^{35}

#### 1. Spin networks as entanglement graphs

In the spin network basis, gluing edges correspond to entangling the degrees of freedom attached to their free ends.^{30} We clarify this point with the following example. Consider the gluing of two four-valent spin network vertices described by the wavefunction *ψ*,

The integral of the spin network basis elements, which is the factor implementing the gluing, is represented in Fig. 3 and performed in the following. To simplify the notation, label 4 is removed from all quantum numbers (e.g., *j*^{4} is denoted just as *j*); we also adopt the notation $n123={n1,n2,n3}$. By substituting Eq. (12) in the expression of Eq. (6), one obtains

The integral of the Wigner matrices, sketched in Fig. 3(a), is well-known in the representation theory and yields

where $Ikk\u2032$ is a bivalent intertwiner in the space $Vj\u2297Vj$ attached to the free ends of the to-be-glued edges

In the expression above, both the state coefficient and the spin network basis elements are contracted with a bivalent intertwiner. Crucially, this is equivalent to projecting $|\psi \u27e9$ on the following state of $Vj\u2297Vj$:

which is a singlet state. The entanglement between the two edges composing the link can be quantified via the von Neumann entropy. Denoting by *ρ _{s}* and

*ρ*the reduced density matrices of the edges attached to the source and target vertex of the link $\u2113$, respectively, one can easily check that

_{t}i.e., the entanglement entropy of the two subsystems reaches its maximum possible value: the state $|\u2113\u27e9$ is maximally entangled.

Therefore, starting from a set of open spin network vertices, the gluing of pairs of their edges is performed by entangling, in a singlet states, the spins on the corresponding free ends. The connectivity pattern of a set of vertices can, thus, be understood as an entanglement pattern among the degrees of freedom attached to the free ends of their open edges.

Finally note that, in the spin network basis, the graph resulting from the gluing of spin network vertices is decorated as follows: every link $\u2113$ carries a spin $j\u2113$, and every vertex *v* carries an intertwiner *ι _{v}*. Open edges, when present, carry an additional quantum number: the (non-contracted) magnetic index at their free end (which causes the spin network to transform non-trivially under gauge transformations acting on its boundary).

#### 2. Constructing entanglement graphs of arbitrary connectivity

The construction of spin network states of arbitrary connectivity *γ* from many-body states associated with *N* open vertices (where *N* is the number of vertices in *γ*) has been rigorously defined in Ref. 30. The first ingredient is a description of the combinatorial structure of graphs in terms of individual colored vertices. In the graph theory, the connectivity pattern of a set of *N* vertices (whose edges are not distinguished by a color) is encoded in the *adjacency matrix*, i.e., a *N *×* N* symmetric matrix *A* defined as follows: the generic element *A _{xy}* takes value 1 if vertices

*x*and

*y*are connected and 0 otherwise. This encoding can be easily generalized to the case in which edges departing from vertices are distinguished by a color

*i*, as it happens with spin network vertices. Assuming the absence of one-vertex loops, the generalized adjacency matrix takes the form

where *A _{vw}* is now a

*d*×

*d*matrix (and $0d\xd7d$ stands for the null

*d*×

*d*matrix) with element $(Avw)ij$ equal to 1 if vertices

*v*and

*w*are connected along edges of color

*i*and

*j*, respectively (i.e., $evi$ and $ewj$ are glued together) and 0 otherwise. To simplify the notation, and since the edge coloring does not play any particular role, one usually assumes that vertices can be connected only along edges of the same color. The matrix

*A*then takes a diagonal form

_{vw}with $avwi$ equal to 1 (0) if vertices *v* and *w* are connected (not connected) along their edges of color *i*; a link formed by $evi$ and $ewi$ is denoted as $\u2113vwi$.

The generalized adjacency matrix defined in Eqs. (19) and (20), thus, encodes the connectivity pattern *γ* of a set of *N* vertices; that is, “who is glued to whom.” The next ingredient for the implementation of *γ* on a set of open vertices is the operator performing the gluing of edges, defined as follows. The operator $P\u2113vwi$ creating the link $\u2113vwi$ acts on the edges $evi$ and $ewi$ by projecting their state onto the subspace characterized by the gluing symmetry (invariance under simultaneous right action of the group)

A spin network state associated with the generic graph *γ* can then be obtained from a set of open vertices in the state $\psi \u2208hN$ by applying to the latter a set of gluing operators according to the adjacency matrix *A* of *γ*,

As follows from Eq. (16), in the spin network basis, the gluing operator is a projection of edge spins onto maximally entangled states. The graph *γ* of the spin network state of Eq. (22) is, thus, realized as a *pattern of entanglement* of a set of vertices. Spin networks regarded as arising from the entanglement structure of states describing a collection of spin network vertices that are also referred to as *entanglement graphs*.

### C. Group field theories

A group field theory^{36,37} (GFT) is a theory of a quantum field $\varphi $ defined on *d* copies of a group manifold *G*. In the GFT model of simplicial quantum gravity, $\varphi $ is a bosonic field whose fundamental excitation is an elementary polyhedron, specifically the $(d\u22121)$-simplex dual to the *d*-valent spin network vertex introduced in Sec. II A. The action of the model takes the following form:

where $g\u2192={g1,\u2026,gd}$, $k(gi(qi)\u22121)$ is the kinetic kernel, responsible for the gluing of polyhedra (spin network vertices), which gives rise to extended spatial geometries (spin network graphs), *λ* is a coupling constant, and $v(gij(gji)\u22121)$ is the interaction kernel, which determines the interaction processes of polyhedra that generate *d*-dimensional spacetime manifolds of arbitrary topology. In particular, due to the simplicial interpretation of field quanta, the Feynman amplitudes of the theory are given by simplicial path integrals (a characteristic shared with simplicial approaches to quantum gravity^{69}) or, equivalently, spin foam models^{35} (representing “histories” of spin networks).

The GFT Fock space is constructed from the Hilbert space $h$ of the $(d\u22121)$-simplex (equivalently, the dual *d*-valent vertex) defined in Eq. (8)

It includes the spin network states in the form of Eq. (22), symmetrized over the vertex labels. Crucially, the symmetry under relabeling of vertices can be understood as a discrete version of diffeomorphism invariance^{30} (which is a necessary condition for background independence), as the vertex labels behave like “coordinates” over the spatial manifold described by the spin network.

Let us finally remark that spin networks arise, in this context, from the entanglement properties of many-body states describing a set of (indistinguishable) spin network vertices. More specifically, the entanglement structure of the many-body state can be identified with the graph formed by the vertices. In the following, we present the correspondence between spin network states and tensor networks, a quantum information language that realizes an analogous graphical encoding of many-body entanglement.

### D. The tensor network perspective

Consider a many-body system composed of *N d*-dimensional spins $s1,\u2026,sN$. A generic state for the system

is described by *d ^{N}* complex coefficients $Cs1,\u2026,sN$. The computational cost of this description can, however, be reduced by considering a

*tensor network decomposition*

^{39}of the state. It consists in replacing the tensor $Cs1,\u2026,sN$ with a collection of smaller tensors $Tisi$ interconnected via auxiliary indices $a\u2192i=ai1,\u2026,air$ (for simplicity, we assume each one having dimension

*D*),

where $Trn$ symbolizes the trace over the auxiliary indices performed according to a combinatorial pattern $n$ of the physical spins, *A* is the adjacency matrix describing the network $n$, and repeated indices are summed over. Note that the number of parameters needed to describe the tensor network has a polynomial scaling in the system size *N*, instead of an exponential one;^{39} in the case we considered, it is given by *NdD ^{r}*.

Spin networks regarded as entanglement graphs (according to the discussion of Sec. II B) formally correspond^{30} to a particular class of tensor networks, called *projected entangled pair states*^{70,71} (PEPSs). A PEPS is a collection of maximally entangled states $|\varphi \u27e9=\u2211a=1D|a\u27e9|a\u27e9$ of pairs of auxiliary systems projected locally onto physical systems $s1,\u2026,sN$ with the entangled pairs corresponding to the *links* of the resulting network $n$. Let $|\varphi \u2113\u27e9$ be the maximally entangled state corresponding to link $\u2113$ of $n$, and let *Q _{i}* be the operator at site

*i*projecting the auxiliary systems onto the physical one

*s*; then

_{i}where the tensor $Tisi$ has elements $(Tisi)ai1ai2\u2026=\u27e8si|Qi|ai1ai2,\u2026\u27e9$. The network $n$, thus, corresponds to the *pattern of entanglement* of the physical spins $s1,\u2026,sN$; in particular, the connectivity of $n$ is realized by pairs of auxiliary degrees of freedom in a maximally entangled state.

Similarly, spin networks can be understood as arising from the entanglement structure of a many-body system, as explained in Sec. II B. The degrees of freedom encoding the connectivity of the spin network are the ones living in the representation spaces attached to the edge free-ends. Specifically, a pair of edges $evi$ and $ewi$ is glued into a link $\u2113vwi$ when the spins living on $Vjvi$ and $Vjwi$ are in singlet states. The spin network counterpart to the link state $|\varphi \u2113\u27e9$ is, thus, the one of Eq. (17) that we rewrite here for clarity

where $jvi=jwi=j$.

Therefore, tensor networks and (completely generic) spin networks have in common the interpretation of links of the graph/network as maximally entangled pairs of systems (auxiliary degrees of freedom for the first, edge-spins for the latter). However, the spin network wave-function $\psi \gamma $ is not, in general, a tensor network; that is, it does not necessarily factorize over single-vertex tensors.

Nevertheless, spin network states obtained from the gluing of open vertices in the factorized state,

do formally correspond to tensor networks. In particular, they can be understood as PEPS, as the gluing procedure is effectively a projection of link states onto single-vertex states

where $\u2113$ is a short notation for the generic link $\u2113vwi,\u2009|\u2113\u27e9$ is the link state defined in Eq. (28), and $f\u2208hj\u2192$ is a fixed-spins vertex state.

Note that, when regarding Eq. (30) as a tensor network, the spins *j* on the graph *γ* correspond to “bond dimensions” of the tensor network indices. However, in the spin network formalism, the spins are not fixed parameters (as are the tensor network bond dimensions), but *dynamical variables*. Therefore, only the “fixed-spins case” given in Eq. (30) formally corresponds to an ordinary tensor network. The generalized case with link and vertex wave-functions spreading over all possible spins, thus, qualifies as a superposition of tensor networks. Furthermore, given the bosonic nature of the discrete entities the individual tensors are associated with, spin networks obtained from factorized many-body states correspond to (superpositions of) *symmetric* tensor networks (for more details, see Ref. 30).

## III. ENTANGLEMENT AND CORRELATIONS ON SPIN NETWORKS TO PROBE AND RECONSTRUCT GEOMETRY

In this section, we give an overview of a series of works on the study of correlations and entanglement entropy on spin networks which, in the spirit of the results presented in Sec. IV, are based on the interplay between quantum gravity and quantum information and/or condensed matter physics. The results are grouped by theme and presented mainly in chronological order.

### A. On the horizon surface: Correlations and bulk entropy

We start with early results on spin networks describing finite regions of 3D space bounded by a causal horizon. On the one hand, these results deal with the computation of the horizon entropy and the recognition of correlations between horizon subregions as responsible for corrections to the entropy area law;^{40} on the other hand, they concern the introduction of the concept of bulk entropy and its relationship with the boundary area.^{41}

#### 1. A black point model for the computation of the horizon entropy

In Ref. 40, Livine and Terno modeled the horizon of a static black hole (at the kinematic level) as a two-sphere made by 2*n* elementary patches, each one punctured by an edge carrying the spin $1/2$. (The argument is as follows: since any representation space *V ^{j}* can be decomposed into the symmetrized product of 2

*j*spin-$1/2$ representations, the spin-$1/2$ patch can be considered as the “elementary patch.”) We denote by

*R*the black hole region, so that its boundary $\u2202R$ corresponds to the horizon two-sphere. The Hilbert space $h\u2202R$ describing the set of boundary edges can be decomposed as

where $Dnj$ is the degeneracy space of states with spin *j*. The gauge-invariant subspace associated with the horizon is then given by the intertwiner space

where the superscript (0) is used to denote the presence of gauge-invariance. In this description, the bulk is, thus, coarse-grained to a single point (hence, the name “black point model”), as depicted in Fig. 4(a). The assumption that the surface is a causal horizon implies complete ignorance of the bulk geometry, and the boundary state is, therefore, given by

where ${|\iota r\u27e9}$ is a basis of the intertwiner space $h\u2202R0$ and *N* is the dimension of the latter. Note that the Boltzmann entropy of such a state coincides with its von Neumann entropy, both being equal to $log\u2009N$. The intertwiner-space dimension *N* is computed via random walk techniques, and the result for the entropy in the asymptotic limit $n\u2192\u221e$ is an area law with a logarithmic correction. The latter is shown to be given by the total amount of correlations between two halves of the horizon surface. Let us show the methodology, as this will be useful for subsequent discussion and for comparing the results presented in Sec. IV to the loop quantum gravity literature.

Consider the splitting of the boundary into a set $\u2202A$ of 2*k* qubits and a complementary set $\u2202B$ of $2(n\u2212k)$ qubits [see Fig. 4(a)]. Then

where $h\u2202A=(V1/2)\u22972k$ and $h\u2202B=(V1/2)\u22972(n\u2212k)$. (Note that such a factorization does not hold for the gauge-invariant subspace $h\u2202R(0)$, see the discussion in Ref. 44.) When decomposing each subspace into a direct sum over irreducible representations *j*, e.g., $h\u2202A=\u2295j=0V\u2202Aj\u2297D\u2202Aj$, the intertwiner states of $h\u2202R(0)\u2282h\u2202R$ turn out to be singlet states on $V\u2202Aj\u2297V\u2202Bj$ with extra indices *a _{j}* and

*b*labeling basis of the degeneracy spaces $D\u2202Aj$ and $D\u2202Bj$, respectively. This corresponds to unfolding the intertwiner as illustrated in Fig. 4(c). The horizon state then becomes as follows:

_{j}It is found that, for $2k=n$ (symmetric splitting of the horizon surface), the mutual information $I\rho (\u2202A:\u2202B)$, amounting to three times the entanglement between $\u2202A$ and $\u2202B$ (quantified, e.g., by the entanglement of formation), equals the logarithmic correction to the horizon entropy.

A possible relationship of the entanglement between $\u2202A$ and $\u2202B$ (for $\u2202A\u226a\u2202B$) with the evaporation process is also suggested, as the case *j *=* *0 corresponds to the detachment of the surface patch $\u2202A$ from the rest of the horizon.

#### 2. Bulk-topology contribution to the boundary entropy

In Ref. 41, Livine and Terno generalized the computation of the horizon entropy performed in Ref. 40 by taking into account the non-trivial structure of the bulk graph. In particular, they promoted the boundary state counting of Ref. 40 to a bulk state counting performed by gauge-fixing the holonomies on internal loops to avoid over-estimating the number of states seen by an external observer. (It is showed that, because of gauge invariance, the bulk degrees of freedom are truly carried by internal loops.) The horizon entropy (evaluated as the logarithm of the number of states supported by a bulk flower-graph with fixed boundary conditions) then turned out to depend on the topology of the graph through its number of loops.

### B. Distance from entanglement

#### 1. Correlations between disjoint regions of a spin network

In Ref. 44, Livine and Terno explored the correlations induced between two disjoint regions *A* and *B* of a spin network from the “outside geometry” *R* [i.e., the region of the spin network complementary to $A\u222aB$, see Fig. 4(b)]. Since $\u2202A\u222a\u2202B=\u2202R$, the gauge invariant state induced on the boundary of the two regions can be regarded as the result of coarse-graining *R* to a single intertwiner, as in the model of the previous paragraph. The presence of correlations between *A* and *B* can then be traced back to the fact that, because of the requirement of gauge invariance of $\u2202R$, the Hilbert space $h\u2202R(0)$ is not isomorphic to $h\u2202A(0)\u2297h\u2202B(0)$. The intertwiner on $\u2202R$ can, in fact, be unfolded into two vertices connected by a “fictitious link” as in Fig. 4(c), and $h\u2202A(0)\u2297h\u2202B(0)$ is recovered as the subspace with an internal link labeled by the trivial representation *j *=* *0 (which effectively corresponds to the absence of connection). The internal link, thus, encodes the entanglement between regions *A* and *B*, induced from the complementary region *R*. This entanglement is then related to a notion of distance between parts of the spin network, building on the idea that, in the absence of a background geometry, such a notion can only be defined in terms of correlations between the quantum degrees of freedom and is expected to be induced from the algebraic and combinatorial structure of the “outside geometry.”

In the same spirit, Ref. 45 by Feller and Livine shows how a notion of distance can be reconstructed from spin network states whose correlations map onto the standard Ising model.

### C. Entanglement entropy and holographic spin networks

Gauge-invariant degrees of freedom are non-local: the Hilbert space of a spin network graph $h\gamma $, indeed, does not factorize into the tensor product of Hilbert spaces describing the subgraphs into which *γ* can be split. We mentioned that in Refs. 40 and 44, this issue is overcome by embedding the intertwiner space into the tensor product of the Hilbert spaces, which are *not* gauge-invariant. (Each one being the tensor product of representations attached to a subset of boundary edges and to a “fictitious” internal-link.) Likewise, in Refs. 46 and 47, Donnelly showed how the entanglement entropy between an arbitrary region *R* of a spin network graph and its complement $R\xaf$ can be computed by embedding $h\gamma $ into an extended Hilbert space that factorizes over *R* and $R\xaf$ with the gauge symmetry broken at the interface of the two regions. More specifically, Ref. 46 takes the complete graph in a spin network basis state: the reduced density matrix *ρ _{R}* is, therefore, completely mixed, and the entropy is given by

An explanation of agreement with the result obtained from the isolated horizon framework in the limit of a large number of punctures is then provided: the spin network states representing the purification of *ρ _{R}* in the two frameworks have a Schmidt decomposition of the same rank. [Note that the result holds only asymptotically: the isolated-horizon entropy is less than Eq. (36), as it includes the gauge-invariance constraint on the boundary $\u2202R$.] Reference 47, instead, takes the whole graph in a completely generic state. The entropy of region

*R*then turns out to be given by the sum of three positive terms: the Shannon entropy of the distribution of boundary representations, the weighted average of $log\u2009(2j+1)$ overall boundary representations

*j*, and a term representing non-local correlations.

An alternative definition of entanglement entropy of regions of a spin network, similarly derived from the embedding of the Hilbert space of gauge-invariant states into an extended Hilbert space, is provided in Ref. 48 and relies on an extension procedure that is based on the excitation content of the theory instead of the underlying graph.

The computation of the entanglement entropy of spin network states and the study of a holographic regime via models and techniques from condensed matter physics is the methodology underlying the results on random spin networks to which this review is dedicated. It has been adopted in earlier work: in Ref. 52, Feller and Livine introduced a class of states inspired by Kitaev's toric code model, which satisfy an area law for entanglement entropy and whose correlation functions between distant spins are non-trivial.

We close this subsection with a general result on boundaries in quantum gravity: in Ref. 72, Bianchi *et al.* showed that boundary states associated with finite portions of spacetime, representing local gravitational processes with certain initial and final data, are mixed, pointing out that such a feature can be regarded as the consequence of tracing over the correlations between the region and its exterior.

### D. Gluing adjacent faces with entanglement

One of the main points of this review concerns the role of entanglement in the connectivity of space. Here, we recall recent results by Bianchi and collaborators on entanglement as a tool for gluing (in the sense specified below) elementary portions of space (spin network vertices). Crucially, this will allow us to differentiate between the various notions of gluing of spin network vertices and to clarify which degrees of freedom are involved in the corresponding entangling procedures. The variety we refer to stems from the distinction between vector geometries, which are defined below, and twisted geometries, of which tensor networks provide a quantum version. In fact, as explained in Sec. II B, a spin network describes a quantum twisted geometry in which “neighboring relations” of quantum polyhedra are codified by links: two intertwiners connected by a link $\u2113$ represent neighboring polyhedra, whose adjacent faces have equal area (determined by the spin $j\u2113$) but different shape and/or orientation, in general. Note that the absence of correlation between the polyhedra of a twisted geometry is translated, at the quantum level of the spin network, to the un-entangled nature of the intertwiner degrees of freedom (i.e., the quantum geometry of neighboring polyhedra has uncorrelated fluctuations). In a vector geometry, the normals to the adjacent faces of neighboring polyhedra are instead anti-parallel, i.e., the two faces adhere to each other, despite the possibly different shape.

In Ref. 50, it was shown that a quantum version of vector geometries can be obtained from a spin network graph by entangling the intertwiner degrees of freedom. They introduced a class of states, called *Bell-network states*, constructed by creating between intertwiners at nearest-neighbour nodes the analogous of the spin-spin correlations of Bell singlet states. These correlations ensure that the normals to the adjacent faces of the corresponding quantum polyhedra are always back-to-back, i.e., that the face planes are parallel. Then, exactly as a Bell singlet state can be understood as a uniform superposition of back-to-back spins over all space directions, a Bell-network state at fixed spins represents a uniform superposition over all vector geometries. In Ref. 51, it was further shown that the entanglement entropy of Bell-network states obeys an area law.

## IV. RANDOM SPIN NETWORKS AND DUAL STATISTICAL MODELS

So far, we introduced the spin network formalism (shared by several approaches to quantum gravity) to describe regions of quantum space (time) and pointed out that entanglement plays a crucial role in this description: it is at the origin of space connectivity. When facing the problem of extracting continuum gravitational physics from such a fundamental description, a crucial issue to be dealt with is the interplay between quantum correlations among the geometric data and global kinematic (and possibly dynamic) geometric features of the spacetime regions considered. Entanglement entropy turned out to be a key tool in this regard.^{18,19,24,25,73}

The computation of the entanglement entropy of spin network states can be highly simplified by the use of random tensor network techniques. This clearly requires to restrict the attention to spin network states given by (superpositions of) random tensor networks. We introduce such a class of states in Sec. IV A and dedicate Sec. IV B to illustrate how random tensor network techniques can be used to translate their Rényi entropies into partition functions of a classical Ising model.

### A. Random spin networks

Consider the class of states (introduced in Sec. II D) that are obtainable from the gluing of a set of vertices each one described by a state *f _{v}*,

where $j\gamma \u2192,\u2009n\gamma \u2192$, and $\iota \gamma \u2192$ are, respectively, spins, magnetic indices, and intertwiners attached to the graph *γ*, and *I* is the bivalent intertwiner introduced in Eq. (15). A coarse graining of these states is then implemented via uniform randomization over the geometric data. The randomization is performed on each vertex separately. This is a necessary requirement for the entropy calculation to be mapped into the evaluation of the free energy of a statistical model (see below). It is also assumed that the spin network states are peaked on specific values $j\gamma \u2192$ of the edge spins. This assumption allows us to work in a fixed spin-sector and, thus, largely simplifies the calculation. The attention is, therefore, restricted to states of the form

where each tensor $(f)n\u2192\iota j\u2192$ is picked randomly from its Hilbert space $hj\u2192$ [defined in Eq. (9)] according to the uniform probability distribution. In the following, we omit the rhs of Eq. (38) that is the explicit reference to the edge spins $j\u2192\gamma $; therefore, unless otherwise stated, $|\psi \gamma \u27e9$ refers to the fixed-spin state of Eq. (38).

### B. Rényi entropy from the Ising partition function

As pointed out in Sec. IV A, our main focus is on finite regions of quantum space described by spin networks corresponding to random tensor networks. In the following, we show that the entanglement content of these states can be conveniently computed via the Rényi entropies.

Given the spin network state $|\psi \gamma \u27e9$ of Eq. (38), consider the reduced state associated with a region *R* of the graph *γ*: $\rho R=TrR\xaf[\rho ]$, where $\rho =|\psi \gamma \u27e9\u27e8\psi \gamma |$ and $TrR\xaf$ is the trace over all degrees of freedom (magnetic indices and/or intertwiners) of the region $R\xaf$ complementary to *R*. The Rényi-2 entropy of *ρ _{R}* is a measure of entanglement given by

The computation of this is quantity that is performed via the *replica trick*, which is based on the possibility to express the trace of a reduced density matrix *ρ _{R}* as a trace over two copies of the density matrix

*ρ*associated with the entire system [here, we assumed

*ρ*to be normalized, i.e., $Tr(\rho )=1$]

where the operator *S _{R}*, called

*swap operator*, acts on the two copies of the Hilbert space $hR$ associated with

*R*as follows:

with $|r\u27e9$ and $|r\u2032\u27e9$ being elements of an orthonormal basis of $hR$. An illustration of the replica trick of Eq. (40) is given in Fig. 5.

By applying the replica trick, the Rényi-2 entropy of region *R* of the spin network described by the state $\rho =|\psi \gamma \u27e9\u27e8\psi \gamma |$ can be written as

where the presence of the denominator takes into account the possible non-normalization of *ρ*, and where the swap operator *S _{R}* acts on two copies of the Hilbert space

associated with the spin network region *R*.

Note that *Z*_{1} and *Z*_{0} are quadratic functions of the random vertex states $\rho v:=|fv\u27e9\u27e8fv|$, and their average is, therefore, easier to compute than the average of the entropy. This leads to the proposal^{74} of expanding the latter in powers of the fluctuations $\delta Z1=Z1\u2212Z1\xaf$ and $\delta Z0=Z0\u2212Z0\xaf$ (the overline is used to denote the average value under randomization of the vertex states),

In Ref. 74, Hayden *et al.* showed that for large enough bond dimensions, which in the present framework correspond to the edge spins, the fluctuations are suppressed, i.e.,

where $\u2243$ refers to asymptotic equality as the edge spins go to infinity. In particular, they proved that, for a tensor network with homogeneous bond dimensions equal to *D*, given an arbitrary small parameter $\delta >0$, it holds

with probability $P(\delta )=1\u2212Dc/D$, where *D _{c}* is a critical bond dimension depending on

*δ*and on the number

*N*of vertices as $Dc\u221d\delta \u22122ecN$, with

*c*being a constant factor.

Thanks to Eq. (45) the computation of the average entropy can, thus, be traced back to the computation of the average quantities $Z1\xaf$ and $Z0\xaf$. Let us focus on $Z1\xaf$, as $Z0\xaf$ is simply given by the latter upon reducing the swap operator to the identity operator. The quantity $Z1\xaf$ can be written as^{31}

where $\rho \u2113:=|\u2113\u27e9\u27e8\u2113|$ and $\rho v:=|fv\u27e9\u27e8fv|$. For each vertex *v*, the average over the two copies of the state *ρ _{v}* can be computed via the Schur's lemma,

^{75}which yields

where $dv$ is the dimension of the vertex Hilbert space $hjv\u2192$ [see Eq. (9)] and *S _{v}* is the swap operator on $hjv\u2192\u2297hjv\u2192$. This is a crucial step since it brings out, as we are going to explain, two-level variables (one for each vertex) corresponding to the spins of a Ising model living on the graph

*γ*. Note in fact that, once Eq. (48) is inserted into Eq. (47), the latter can be written as

where $\sigma v=\xb11$ is a two-level variable associated with vertex *v*, $\sigma \u2192={\sigma 1,\u2026,\sigma N}$, and

is a constant factor. That is, $Z1\xaf$ has been written as a sum of $2N$ terms involving the identity ($I$) or the swap operator (*S _{v}*) for each of the

*N*vertices, and the variable

*σ*encodes the presence of one or the other ($I$ for $\sigma v=+1$ and

_{v}*S*for $\sigma v=\u22121$) in every term of the sum.

_{v}Given the form of the vertex Hilbert space $hjv\u2192$, the swap operator *S _{v}* factorizes as follows:

i.e., into a swap operator $Sv0$ for (the double copy of) the intertwiner Hilbert space $ij\u2192v$ and a swap operator $Svi$ for (the double copy of) the representation space $Vjvi$ on each edge $evi$, as shown in Fig. 5(d). Crucially, the same applies to the swap operator *S _{R}*,

Consequently, to every open edge $evi$ of the graph *γ*, one can attach a two-level variable $\mu vi=\xb11$ (also called *pinning spin*^{74}), encoding whether ($\mu vi=\u22121$) or not ($\mu vi=+1$) an additional swap operator acts on (the double copy of) its Hibert space; that is, whether or not it belongs to region *R*. The same holds true for the intertwiner on each vertex *v* of the graph for which the two-level variable $\nu v=\xb11$ is introduced.

By performing the trace in Eq. (49), one finally obtains the quantity $Z1\xaf$ that corresponds to the partition function of a classical Ising model

with $A1(\sigma \u2192)$ being the Ising action

where *d _{j}* is the dimension of the representation space

*V*and $Dj\u2192$ the dimension of the intertwiner space $ij\u2192$ (see Sec. II A). Note that the Ising model is defined on the graph

^{j}*γ*: Eq. (54) involves interactions between nearest neighbors Ising spins, where the adjacency relationship is determined by

*γ*(two Ising spins interact only if the corresponding vertices are connected by a link); every Ising spin also interacts with the pinning spins located at its vertex (e.g., the Ising spin

*σ*of a vertex

_{v}*v*on the boundary interacts with the pinning field

*ν*on the intertwiner of

_{v}*v*and with the pinning field $\mu vi$ on the open edge $evi$ of

*v*).

As far as $Z0\xaf$ is concerned, we pointed out that it corresponds to $Z1\xaf$ with $R=\u2205$ (in fact $S\u2205=I$). Therefore, it holds that $Z0\xaf=\u2211\sigma \u2192e\u2212A0(\sigma \u2192)$, where *A*_{0} is given in Eq. (54) with all pinning spins equal to +1,

Note also that, since $Z0\xaf$ and $Z1\xaf$ enter $S2(\rho R)\xaf$ only via their ratio, the computation of the entropy the constant factor in Eqs. (54) and (55) is irrelevant; we, therefore, omit it in the following.

To study the properties of the partition function $Z1\xaf$, it is useful to rewrite the Ising action $A1(\sigma \u2192)$ in the form $A1(\sigma \u2192)=\beta H1(\sigma \u2192)$, where $\beta :=dj$ with *j* being the average spin on *γ*, and

The parameter *β* then plays the role of inverse temperature of the Ising model. As we are working in the high spins regime, the partition function $Z1\xaf$ is dominated by the lowest energy configuration

The same applies to $Z0\xaf$ and, since $min\sigma \u2192H0=0$ [where *H*_{0} is given in Eq. (56) with $\mu vi=\nu v=+1$$\u2200v,evi\u2208\gamma $], it holds that

Therefore, the average entropy can be finally computed via the following formula:

with *β* being the average dimension of the edge spins and $H1(\sigma \u2192)$ being the Ising-like Hamiltonian defined in Eq. (56).

## V. HOLOGRAPHIC ENTANGLEMENT IN SPIN NETWORK STATES

We present recent works that explored the connection between holographic features of regions of quantum space and entanglement of their quantum geometric data for spin network states obtainable from the gluing of random vertex states.

### A. Bulk-to-boundary quantum channels: Isometric mapping of quantum-geometric data

Colafranceschi *et al.*^{31} analyzed the flow of information from the bulk to the boundary of regions of quantum space described by the class of spin network states defined in Eq. (38) to determine under which conditions such a flow can be holographic.

Let us start by providing the definitions of bulk and boundary of a spin network, as given in Ref. 31. Consider a spin network with combinatorial pattern *γ* and edge spins $j\gamma \u2192$. The *boundary* consists in the set of open edges of *γ* (denoted by $\u2202\gamma $) decorated by the respective spins and is described by the Hilbert space

let $|n\xaf\u27e9:=\u2297e\u2208\u2202\gamma |jene\u27e9$ be the basis element of the boundary space $h\u2202\gamma $. The *bulk* is the set of vertices of *γ* (denoted by $\gamma \u0307$) together with the intertwiners attached to them and is described by the Hilbert space

let $|\iota \xaf\u27e9:=\u2297v|j\u2192v\iota v\u27e9$ be the basis element of the bulk space $h\gamma \u0307$.

The flow of information from the bulk to the boundary is identified with the bulk-to-boundary map that every spin network state implicitly defines once regarding the bulk space as *input* and the boundary space as *output*. More specifically, every spin network state of the form

(to simplify the notation, we omitted the edge spins, as they are fixed) can be regarded as a *map*$m$ from the bulk to the boundary Hilbert space, having components

The map $m$ associated with $|\varphi \gamma \u27e9$, therefore, acts on a generic bulk state $|\zeta \u27e9\u2208h\gamma \u0307$ as follows:

i.e., by evaluating the spin network state on $|\zeta \u27e9$ or, in tensor network language, by feeding the bulk input with $|\zeta \u27e9$ (see Fig. 7).

The reduced (and normalized) bulk state takes the form

where $\rho \gamma =|\varphi \gamma \u27e9\u27e8\varphi \gamma |$ and $D\gamma \u0307$ is the dimension of the bulk Hilbert space $h\gamma \u0307$. It follows from Eq. (65) that if the reduced bulk state is maximally mixed, namely, $\rho \gamma \u0307=I/D\gamma \u0307$, the map $m$ is an isometry, i.e., $m\u2020m=I$. Moreover, the corresponding superoperator on the space of bulk operators, $\Lambda (\xb7):=m\xb7m\u2020$, is a completely positive trace preserving (CPTP) map with a Choi–Jamiołkowski state

where

is a maximally entangled state of two copies of the bulk (see Fig. 6).

Colafranceschi *et al.*^{31} studied the bulk-to-boundary map $m$ of a spin network state of the form of Eq. (38) to analyze the relationship between the combinatorial structure and geometric data of a spin network on the one hand, and the isometric character of the corresponding map on the other hand. The latter is quantified via the Rényi-2 entropy of the reduced bulk state [see Eq. (65)]. Thanks to the random nature of the vertex tensors, the entropy is computed via an Ising partition function according to the technique illustrated in Sec. IV B. In particular,

with $H1(\sigma \u2192)$ being the Ising-like Hamiltonian

It is found that spin network graphs made of four-valent vertices (dual to 3D spatial geometries) with an homogeneous assignment of edge spins do not realize an isometric mapping of data from the bulk to boundary. Coherently, increasing the inhomogeneity of the spins assigned to a spin network with four-valent vertices increases the “isometry degree” of the corresponding bulk-to-boundary map.

Let us close this section by commenting on the comparison of this work with Ref. 54, where the idea of interpreting the spin network states as maps from the bulk to the boundary first appeared. Chen and Livine^{54} pointed out that spin network wavefunctions with support on an open graph can be regarded as linear forms on the boundary Hilbert space (the space of spin states living on the open edges of the spin network), and that coarse-graining the bulk, i.e., integrating over the bulk holonomies, then induces a probability distribution for the boundary degrees of freedom. Based on that, they proved the following: any boundary density matrix can be obtained via the bulk-to-boundary coarse-graining procedure from a pure bulk state with support on a graph composed of a single vertex connecting all boundary edges to a single bulk loop. A crucial difference between the map of Ref. 54 and $m$ is that the latter does not perform a coarse graining of the bulk (intended as tracing out the bulk holonomies); instead, it *evaluates* the (pure) spin network state on a given bulk configuration (specifically, a given state for the intertwiner degrees of freedom), thereby yielding a boundary state. Consequently, the latter is a pure state if the bulk input state is pure. In contrast, the boundary density matrix resulting from the bulk-to-boundary coarse-graining of Ref. 54 applied to a pure spin network state is typically mixed.

### B. Holographic states and black hole modeling

As illustrated in Sec. V, Ref. 31 investigated holography on spin network states having the form of Eq. (38), regarding them as maps from the bulk to the boundary. Inspired by similar questions, Ref. 32 studied the same class of states from a different perspective: it analyzed the boundary states returned by the bulk-to-boundary map on varying the bulk input state. That is,

where $m$ is the bulk-to-boundary map corresponding to the spin network state (and a random tensor network) $|\psi \gamma \u27e9$, defined in Eq. (38); $|\zeta \u27e9\u2208h\gamma \u0307$ is the input bulk state, and $|\eta \u27e9\u2208h\u2202\gamma $ is the output boundary state. In particular, it focused on the entanglement content of portion *A* of the output boundary state (see Fig. 7). Again, given the random character of the state, the entanglement measure considered is the Rényi-2 entropy, computed via the Ising model. The result is as follows:

where

with $\zeta \u2193$ being the bulk state reduced to the region with Ising spins pointed down. From Eq. (72), one can note that every misalignment between the Ising spins *σ _{v}* and

*σ*on a link $\u2113vwi$ carries a contribution to the entropy equal to $(log\u2009djvwi)/\beta $, i.e., to (the logarithm of) the dimension of that link, normalized by

_{w}*β*(the average value that quantity can take). The same holds for the pinning spin $\mu vi$ and the Ising spin

*σ*on a boundary edge $evi$. As a result, the first two terms of the rhs of Eq. (72) provide the “area” of the Ising domain wall, i.e., the surface separating the spin-down region (externally bounded by

_{v}*A*) from the spin-up region, where the area is given by a weighted sum of the links crossing it (with weights proportional to the logarithm of the link dimensions). Let $\Sigma (\sigma \u2192)$ be the aforementioned surface for the Ising configuration $\sigma \u2192$, and

its area, as defined above. The Ising Hamiltonian of Eq. (72) can then be written as follows:

Combining Eq. (71) with Eq. (74), one then finds that, for $S2(\zeta \u2193)\u226a\beta \Sigma (\sigma \u2192)$, the Rényi-2 entropy follows an area law with a small correction deriving from the bulk entanglement (see Fig. 8):

For $S2(\zeta \u2193)=O(\beta \Sigma (\sigma \u2192))$, instead, the Rényi-2 entropy follows an “area + volume law,”

In fact, $S2(\eta A)\xaf$ depends on a comparable extent on the entanglement content of the surface $\Sigma (\sigma )\u2192$ (link entanglement) and of the spin-down region bounded by it (intertwiner entanglement in *ζ*).

In Ref. 32, it was also showed that increasing the entanglement content of a region of the bulk can turn the boundary of that region into a horizon-like surface (see Fig. 9), as the Ising domain wall that determines the entropy cannot access it. Notably, this result can be regarded as a realization of the proposal made by Krasnov and Rovelli^{21} of defining a quantum black hole as the part of a spin network that does not influence observables at infinity.

## VI. DISCUSSION

We reviewed very recent work that contributes to the research effort, which looks at holography not as an asymptotic global property (as it was originally conceived), but as a quasi-local property entering the description of finite spatial boundaries (spacetime corners). The main novelty of the illustrated approach is the use of spin network states formally corresponding to (generalized) random tensor networks. The defining feature of this class of states is the randomness of the wavefunctions associated with the individual spin network vertices, which has the remarkable property of mapping the correlations of the spin network states to that of a classical Ising model living on the same graph. This enables to investigate the entanglement content of the spin network by relying on standard condensed matter and quantum information techniques. Moreover, the randomization over vertex wavefunctions can be understood as a local coarse graining on the vertex data and, thus, makes this type of states of immediate interest for GFT cosmology.^{76–79}

Reference 31 specifically studied the flow of information from the bulk to the boundary through the Choi–Jamiołkowski duality, computing the Rényi entropy of the Choi–Jamiołkowski state through a random tensor technique that traces it back to the evaluation of Ising partition functions. The result is a positive correlation between the inhomogeneity of the edge spins and the “isometry degree” of the bulk-to-boundary map. The same technique is applied in Ref. 32 to the computation of the Rényi entropy of boundary states and leads to the derivation of (an analogue of) the Ryu–Takayanagi formula.^{18,19} Interestingly, Ref. 32 also showed that the presence of a bulk region with high entanglement entropy can turn the boundary of that region into a horizon-like surface, hereby offering a concrete example of the definition of quantum black holes given in Ref. 21 with a picture that recalls the “quantum graphity” of Ref. 80.

The illustrated work paves the way to an extensive application of quantum information tools to the study of the spacetime microstructure and the modeling of quantum black holes. In particular, the superposition of graphs (which is necessary to bring the analysis at the dynamical level) may be implemented by enriching the spin network structure with data encoding the amount of link-entanglement between vertices and using such data to manipulate the combinatorial structure of the graph, analogously to what has been done for random tensor networks.^{81} As far as an information-theoretic characterization of black hole horizons is concerned, the illustrated techniques are, for example, expected to enable the derivation of a “threshold condition” for the emergence of horizon-like surfaces in finite regions of quantum space, analogously to the one obtained from the typicality approach to the study of the local behavior of spin networks.^{82}

While the present article covered only a particular corner of the burgeoning field at the crossroads between quantum information and gravity, it is hoped our focused review might inspire further research and continue to motivate fruitful cross fertilization of methods and concepts between these two cutting-edge areas of theoretical physics, ultimately leading to their unification or confluence within a more fundamental theory yet to be discovered.

## ACKNOWLEDGMENTS

The authors would like to thank Goffredo Chirco, Daniele Oriti, and Aron Wall for useful discussions and comments. E.C. acknowledges funding from the DAAD via the scholarship programme “Research Grants-Short-Term Grants, 2021” and thanks the Ludwig Maximilian University of Munich for the hospitality.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

## DATA AVAILABILITY

Data sharing is not applicable to this article as no new data were created or analyzed in this study.