Arxiv:1705.01964V2 [Hep-Th] 27 Jul 2017 Where Min(#Cut) Is the Minimal Number of Tensor Network Armin Links Cut by a Surface Anchored to A

Discrete Gravity on Random Tensor Network and Holographic Renyi´ Entropy

Muxin Han1, 2 and Shilin Huang3 1Department of Physics, Florida Atlantic University, 777 Glades Road, Boca Raton, FL 33431, USA 2Institut fürQuantengravitation, UniversitätErlangen-Nürnberg, Staudtstr. 7/B2, 91058 Erlangen, Germany 3Institute for Interdisciplinary Information Sciences, Tsinghua University, Beijing 100084, China In this paper we apply the discrete gravity and Regge calculus to tensor networks and Anti-de Sitter/conformal field theory (AdS/CFT) correspondence. We construct the boundary many-body quantum state Ψ using random | i tensor networks as the holographic mapping, applied to the Wheeler-deWitt wave function of bulk Euclidean discrete gravity in 3 dimensions. The entanglement Renyi´ entropy of Ψ is shown to holographically relate to | i the on-shell action of Einstein gravity on a branch cover bulk manifold. The resulting Renyi´ entropy S of Ψ n | i approximates with high precision the Renyi´ entropy of ground state in 2-dimensional conformal field theory (CFT). In particular it reproduces the correct n dependence. Our results develop the framework of realizing the AdS3/CFT2 correspondence on random tensor networks, and provide a new proposal to approximate the CFT ground state.

PACS numbers:

I. INTRODUCTION to local degrees of freedom in the bulk [31–33]. The architec- ture of the TN may be viewed as a process of real-space renor- The tensor network is a quantum state of many-body sys- malization, such as multiscale entanglement renormalization tem constructed by contracting tensors according to a network ansatz (MERA), where the renormalization scale relates to the graph with nodes and links (FIG.1). It is originated in con- coordinate of the emergent dimension [2, 7, 34, 35]. The fea- dense matter physics because tensor network states efficiently ture of tensor network makes it an interesting tool for realizing compute ground states of many-body quantum systems [1, 2]. the AdS/CFT correspondence constructively from many-body In addition, the tensor network has wide applications to quan- quantum states. Among many recent progress, an important tum information theory by its relation to error correcting codes result is to reproduce the HEE on the tensor network. and quantum entanglement [3]. It recently relates to quantum The optimized scale-invariant MERA tensor networks are machine learning [4], as well as neuroscience [5]. able to approximate ground states of (discreted) CFTs e.g. the One of the fascinating developments of the tensor net- critical Ising and Potts models. It is observed that the entan- work is the recent relation to the AdS/CFT correspondence glement entropy of optimized MERA is proportional to the and emergent gravity program [6, 7]. The AdS/CFT corre- number of links on a minimal cut through the tensor network, spondence proposes that the quantum gravity theory on d- as a nice realization of the RT formula [6, 7, 36]. Moreover dimensional Anti-de Sitter (AdS) spacetime is equivalent to a there has been two recent interesting exactly-solvable models conformal field theory (CFT) living at the (d 1)-dimensional that realize the RT formula on tensor networks, including [24] boundary of AdS. It offers a dictionary between− the observ- with tensor networks with perfect tensors, and [25] with ran- ables of the d-dimensional bulk gravity theory and those of dom tensor networks. Given a boundary region A containing a the (d 1)-dimensional boundary CFT. Properties of the bulk number of open links of the tensor network, the entanglement gravity− and geometry may be reconstructed or emergent from entropy of the tensor network in each of the above approaches the boundary CFT, known as the emergent gravity program reproduces an analog of the RT formula (See e.g. [37–45] for [8, 9]. As an important ingredient of AdS/CFT, the bulk geom- some more recent developments) etry relates holographically to the entanglement in the bound- S EE(A) = Min(#cut) ln D, (2) ary CFT, via the Ryu-Takayanagi (RT) formula ·

arXiv:1705.01964v2 [hep-th] 27 Jul 2017 where Min(#cut) is the minimal number of tensor network Armin links cut by a surface anchored to A. D is the bond dimension S EE(A) = , (1) 4GN (range of tensor index). Eq.(2) is in general an upper bound in the tensor network, which is saturated in the case of per- which identifies the entanglement entropy S EE(A) of a (d 1)- fect or random tensor network. The random tensor network − dimensional boundary region A with the area Armin of the bulk approach has been related to the quantum geometry in loop (d 2)-dimensional minimal surface anchored to A [10–15]. quantum gravity (LQG) ([46–49] for reviews), which relates − GN is the Newton constant in d dimensions. S EE(A) satisfy- Eq.(2) to the geometrical RT formula Eq.(1) [50]. ing RT formula is referred to as the holographic entanglement Although the perfect and random tensor network ap- entropy (HEE) [16–22] proaches give exactly solvable models and have nice fea- The tensor network can be understood as a discrete version tures, it is known that both approaches suffer the issue of of the AdS/CFT correspondence [23–30]. The tensor network flat entanglement spectrum. Although the entanglement (Von states approximates the CFT state at the boundary, while the Neumann) entropy Eq.(2) is consistent with the RT formula, structure of tensor network emerges an bulk dimension built Renyi´ entropies S n(A) from both approaches are all identical by layers of tensors. Tensors in the tensor network correspond to Eq.(2) with trivial n dependence. But the (generalized) RT 2

FIG. 2: A simple trivalent tensor network with 2 nodes and 5 links (4 open links). The tensor network is dual to a triangulation with 2 FIG. 1: An example of tensor network with rank-6 tensors. triangles. The tensor network state Ψ is given by an expansion with cer- | i tain basis in the Hilbert space of many-body system Ψ = P Q | i (T ) a , a , a . The coefficients is constructed ai,γl p p γ1 γ2 ... ai... 1 2 N As another interesting aspect of our work, we show that { } | ···T i p by distributing a rank-6 tensor p at each 6-valent node , such that the random tensor network effectively relates the boundary each tensor index associates to a link adjacent to p. Connecting 2 nodes by a link means contracting the corresponding indices γ . Renyi´ entropy to the partition functions of bulk gravity on l the branched cover spacetime. This relation has been a key assumption in the existing AdS/CFT Renyi´ entropy compu- formula of Renyi´ entropy has a nontrivial n dependence since tations [11, 12, 15], while our work derives this relation in the CFT Renyi´ entropy does [15]. For instance, in any 2d CFT the context of AdS3/CFT2. We also derive the duality be- ´ (CFT2), the ground state has the universal Renyi´ entropy [51] tween the boundary Renyi entropy and bulk cosmic brane (in AdS /CFT ), which has been proposed as a conjecture in [15]. ! ! 3 2 1 c lA In this paper, we focus on 2d CFT and 3d bulk spacetime S n(A) = 1 + ln , (3) n 6 δ (AdS3/CFT2) in Euclidean signature. The CFT state Ψ is constructed by implementing bulk gravity dynamics to| ran-i where c is the central charge, lA is the length of the region A, dom tensor network states studied in [32]. The random tensor and δ is a UV cut-off. The CFT Renyi´ entropy manifestly has network constructed in [32] has random tensors at all nodes a nontrivial n dependence. p, and has labels ap,p0 on links (p, p0). Each ap,p0 labels the The mismatch of Renyi´ entropy implies that both tensor non-maximal entangled state ap,p0 on each link. The ten- network states in [24, 25] fail to approximate the CFT ground sor network is dual to a tiling| of 2di spatial slice Σ (FIG.2). state. An idea of the reason behind it may be extracted The entanglement entropy of a relates to the length L of | p,p0 i ` from the AdS/CFT: In the AdS/CFT, the CFT ground state the edge ` intersecting (p, p0) in the tiling. Thus each random (at strong-coupling) is dual to the bulk semiclassical AdS tensor network as boundary CFT state, denoted by ~a , deter- | i spacetime geometry. However the tensor networks designed mines a set of edge lengths L` in the bulk. When the tiling is a in [24, 25] only consider the geometry of a spatial slice in triangulation, edge lengths uniquely determines a discrete ge- AdS, without any input about time evolution. In the contin- ometry on Σ which approximates the smooth geometry as the uum AdS/CFT context, the correct entanglement spectrum are triangulation is refined [52]. We make linear combinations of obtained by considering the spacetime geometry, and taking ~a to write CFT states as into account the bulk dynamics given by the Einstein equation | i X [11, 12, 14, 15], while [24, 25] don’t have any dynamical in- Ψ = Φ(~a) ~a . (4) | i | i put. In this work, we take the above idea as a hint to improve ~a the tensor network approach, in which we resolve the issue. Note that the optimized MERA for the CFT can derive the where the coefficients Φ(~a) can be understood as a wave func- correct Renyi´ entropies as well as the non-flat entanglement tion of bulk geometry, by the relation between ~a and bulk edge spectrum [36]. However the goal of this work is to finding lengths. Eq.(4) is a holographic mapping from the bulk state analytically a generic family of tensor networks that provide Φ to the boundary state Ψ. the correct CFT Renyi´ entropies without going through any As mentioned above, the CFT ground state is expected dual numerical optimization at all. Our results develop important to a physical state in the bulk corresponding to the semi- analytic methods for the tensor network AdS/CFT. classical AdS spacetime geometry. But the link labels ~a in In this work, we resolve the above Renyi´ entropy issue by Φ(~a) only relate to the geometry of 2d spatial slice Σ. To let having the bulk dynamics on random tensor network states. Φ(~a) encode spacetime geometry, we propose Φ(~a) to be the The dynamics relates tensor networks to the discrete gravity Wheeler-deWitt wave function. Namely, Φ(~a) is a path inte- model known as Regge calculus. We construct the state Ψ gral of (Euclidean) Einstein gravity on spacetime M whose which is proposed as an approximation to the CFT ground| i boundary contains Σ (FIG.3(a)). The geometry on Σ deter- state. As is anticipated by our proposal, Ψ indeed repro- mined by ~a is the boundary condition of the path integral. We duces correctly the RT formula and CFT ground| i state Renyi´ construct Φ(~a) to sum all possible bulk spacetime geometries entropy S n with correct n dependence. Our results develop satisfying the boundary condition, while in the semiclassical the framework of the AdS/CFT correspondence on random limit, it localizes at the classical AdS spacetime in 3d. The tensor networks, and provide a new proposal to approximate semiclassical limit relates to the large bond dimension of ten- the large-N CFT ground state. sor network. 3

⌃ ⌃

⇤

⌃ ⇤

FIG. 4: The (triangulated) manifold Mn (n = 2) made by gluing 2n copies of M.

FIG. 3: (Triangulated) 3-manifolds M in (a), M¯ in (b), and M1 in (c). The tensor network link variables ~a are the boundary data in the spatial slice Σ, while the boundary data at other boundaries (including Here Armin is the geodesic length in AdS3. The above result the infinitely past/future) doesn’t affect our discussion. of S n(A) gives the Renyi´ entropy Eq.(3) of CFT2 ground state with correct n dependence. Our result is valid in the regime of large CFT central charge. Since ~a are data of discrete geometry, Φ(~a) is the discrete This work applies discrete geometry method such as Regge version of Wheeler-deWitt wave function: It is a path integral calculus to study tensor networks (see [55] for other applica- of Regge calculus. Regge calculus is a discretization of Ein- tion of discrete gravity in AdS/CFT). Ψ encodes the dynam- | i stein gravity by triangulating spacetime geometries [53]. The ics of bulk geometries, which is given by the discrete Einstein discrete spacetime geometry is given by the edge lengths in equation. It is interesting to further understand how the bulk the triangulation, known as the Regge geometry. Regge ge- dynamics might relate to the dynamics of boundary CFT, and ometries converge to smooth geometries as triangulations are whether a boundary CFT Hamiltonian might be induced from refined [52]. Φ(~a) is a sum over all (Euclidean) Regge ge- the bulk dynamics. It is interesting to compare our proposal ometries on the spacetime M, weighted by the exponentiated of CFT ground state Ψ to the MERA approach, to extend | i Einstein-Regge action [54]. The detailed explanations of Φ(~a) our result to finite central charges and understand the relation and random tensor networks are presented in Section II. between Ψ and concrete many-body systems such as criti- | i The Renyi´ entropy S n(A) of Ψ at arbitrary n 2 is com- cal spin-chains. The understanding of these aspects should puted in Section III. The computation| i involves≥ averages of develop tensor network models to realize the AdS/CFT cor- the random tensors at nodes p in the tensor networks [25]. It respondence at the dynamical level. The research on these turns out that the random average effectively glues 2n copies aspects is currently undergoing. of path integrals Φ(~a)(Φ(~a)∗) on M, and relate S n(A) to the path integrals of gravity on branch cover spacetimes made by 2n copies of M. We derives that in the bulk semiclassical limit II. RANDOM TENSOR NETWORK AND equivalent to large bond dimensions in the tensor networks, WHEELER-DEWITT WAVE FUNCTION

1 In this work we consider trivalent random tensor network S n(A) [IBulk(Mn) nIBulk(M1)] , (5) ' 1 n − states. A tensor network is viewed as a discrete 2d spatial − slice Σ of 3d bulk spacetime. It is made by a large number where I (M ) is the on-shell gravity action evaluated at Bulk n of trivalent random tensors 3 at each tensor the bulk solution on the branch cover manifold M (FIG.4). p ⊗ p n network node p. The Hilbert|V spacei ∈ H is of≡H dimension D. We The bulk solution has the Z replica symmetry. Eq.(5) has n decompose into a number of subspacesH V and been an assumption in AdS/CFT derivations of HEE in e.g. a a denote dim(VH) d[a]. Each internal link (p,Hp ) 'of the⊕ tensor [11, 12, 14, 15]. But it is now derived from Ψ and ran- a 0 network associates≡ with a maximal entangled state in V V dom tensor networks. The derivation also gives| thei duality a a of certain a, ⊗ between the boundary Renyi´ entropy and bulk cosmic brane (in AdS /CFT ), which has been proposed as a conjecture d[a] 3 2 1 X in [15]. As a result, we show that S n(A) reproduces the RT a = µ µ (7) | p,p0 i | ip ⊗ | ip0 formula for holographic Renyi´ entropy for 2d CFT (Hung- √d[a] µ=1 Myers-Smolkin-Yale formula in [14]) where µ is a basis in . It satisfies a b = δ .A ! | ip H h p,p0 | p,p0 i ab 1 Armin class of random tensor networks ~a can be defined by the S n(A) 1 + . (6) | i ' n 8G (partial) inner product between p at all p and ap,p on all N |V i | 0 i 4 internal links L` on the triangulation defines a discrete metric of Regge geometry{ } [53, 54, 56]. We define Φ(~a) to be a path integral of ~a = p,p app p p . (8) | i ⊗ 0 h 0 | ⊗ |V i discrete gravity on the triangulation by summing over all L` The inner product takes place in at each end point p or p0 in the bulk of M H N∂ X of each link. ~a is a state in the boundary Hilbert space , S (M) | i H Φ(~a):= e− Regge (11) where N∂ is the number of open links. L The label ~a relates to the amount of entanglement on each ` internal link (p, p0). The entanglement entropy S ( ap,p ) of The boundary condition at Σ is L` Σ = 4`P ln d[ap,p ]. We use | 0 i P ⊂ 0 a is ln d[a ], where d[a ] is effectively the bond di- instead of integration because L` are assumed as discrete | p,p0 i p,p0 p,p0 L` mension on (p, p0) in ~a . data, to be consistent with L` Σ. S Regge(M) is the Regge ac- The above class of tensor| i network states is proposed in [32]. tion of Euclidean gravity on⊂ the triangulated 3-manifold M We consider the boundary state Ψ as a linear combination evaluated at the discrete metric L` : | i { } [68]   X 1  X X V(M) S Regge(M) =  L` ε` + L` Θ`  . Ψ = Φ(~a) ~a . (9) −8π`  − 2  | i | i P ` bulk(M) ` ∂M LAdS ~a ⊂ ⊂ Here we understand the trivalent tensor network to be dual LAdS relates to `P and the central charge of CFT by LAdS = 2 to a triangulation of the spatial slice Σ (FIG.2). Namely, each 3 c`P. In the following, the large bond dimension limit

p ∆ d[ap,p0 ] 1 will be taken in the computation of Renyi en- node located at the center of a triangle p in the triangulation. Each link (p, p0) intersects transversely an internal edge tropies. The limit corresponds to the bulk semiclassical limit ` shared by 2 triangles ∆ , ∆ . Open links in tensor network `P 0 and large central charge c 1 as both L` and LAdS p p0 → intersect transversely the edges at the boundary of triangula- are kept fixed. ε` is the bulk deficit angle hinged at the bulk tion. edge `. ε` is a discretization of the bulk curvature. Each bulk The label ~a is understood as the discrete geometry in the edge ` is shared by a number of tetrahedra t. In each t, the bulk of Σ [32], in the sense that the edge length L` of ` in- dihedral angle between 2 faces joint at ` is denoted by θ(t, `). tersecting (p, p0) is proportional to the entanglement entropy The deficit angle is defined by S ( ap,p ) on each link: X | 0 i ε = 2π θ(t, `), ` bulk. (12) ` − ⊂ L` 4`P ln d[ap,p ] (10) t, ` t ≡ 0 ⊂ Here `P = GN ~ is the Planck length in 3d. In this proposal, Each boundary edge ` is shared by 2 boundary triangles. Θ` the bulk geometry is understood as emergent from the entan- is the angle between their outward pointing normals, equiva- glement in tensor network state. The relation can be obtained lently from the recent proposal of understanding tensor networks as X ff Θ` = π θ(t, `), ` boundary. (13) the e ective theory from coarse graining quantum gravity at − ⊂ t, ` t Planck scale [50], in which one derives that the bond dimen- ⊂ L`/4` sion d[ap,p ] of tensor network ~a satisfies d[ap,p ] e P . 0 0 Θ` relates to the boundary extrinsic curvature. The 1st term in By the relation between ~a and| i bulk geometry,' Eq.(9) is a S Regge is the discretization of Ricci scalar term of Einstein- boundary state by summing over all bulk spatial geometry on Hilbert action, while the 2nd term is the discretization of Σ, while Φ(~a) is a wave function of bulk geometry. Eq.(9) Gibbons-Hawking boundary term [57]. The last term is the defines a holographic mapping from the bulk states of geom- cosmological constant term where V(M) is the total volume etry to the boundary states of CFT. We propose the following of M. All quantities ε`, Θ`, and V(M) are determined by edge boundary state Ψ whose bulk wave function Φ(~a) (pre-image | i lengths L`. of the holographic mapping) is an Wheeler-deWitt wave func- In order to be the boundary condition of Regge geometry, tion in 3d Euclidean gravity. Namely Φ(~a) is a path integral of L` Σ = 4`P ln d[ap,p ] have to be the edge lengths of hyper- gravity on a 3d solid cylinder M, whose boundary includes Σ ⊂ 0 bolic triangles, which triangulate Σ. L` Σ have to be a discrete in addition to the boundary where CFT lives (FIG.3(a)). The metric of Σ, which constrains the allowed⊂ data ~a entering the geometry ~a is the boundary condition on Σ in the path integral. sum in Eq.(9). The path integral may also depend on the boundary conditions Applying Φ(~a) in Eq.(11) to the holographic mapping at other boundaries of M. But we make those boundary condi- Eq.(9), we obtain a boundary CFT state Ψ , and we propose tions implicit since they play no role in the following analysis. the resulting Ψ to be the ground state of| thei boundary large- Since ~a gives a discrete geometry with a set of edge lengths | i N CFT, in the bulk semiclassical regime `P L` equivalent L`, more precisely, Φ(~a) is a discrete version of the Wheeler- to the large bond dimension limit by Eq.(10). As `P L`, the deWitt wave function. Indeed, we consider a sufficiently path integral Φ(~a) localizes at the solution of equation of mo- refined triangulation of M, and impose discrete metrics on tion (deriving equation of motion uses the Schlafli¨ identity of the triangulation. Namely each tetrahedron in the triangula- 2 P hyperbolic tetrahedra δV(t)/L = ` t L`δθ(t, `), see e.g. tion carries a 3d hyperbolic geometry with constant curva- − AdS ⊂ 2 [56]) ture L−AdS . Tetrahedron edges are geodesics in the hyperbolic− space, and have edge lengths L . The set of edge lengths ε = 0, ` bulk(M). (14) ` ` ∀ ⊂ 5

n which is the discretized Einstein equation in 3d. Vanishing The random average trρA relates to average n copies of ran- ε everywhere means that the 3d Regge geometry is a smooth dom tensors at each node p. Taking an arbitrary ref- ` |VpihVp| Euclidean AdS3. So Φ(~a) is a semiclassical wave function of erence state 0p p, the random tensor p = Up 0p with bulk AdS geometry. The holographic mapping is expected to U unitary transformation| i ∈ H on 3 .|V Thei Haar| randomi 3 p H ⊗ ≡ Hp map the bulk semiclassical state of AdS3 to the ground state average is given by [25, 60]: of boundary CFT2. Z n In Section III, we check our proposal by computing the n ⊗ p p ⊗ = dUp Up 0p 0p Up† Renyi´ entropies S n of the state Ψ . We show that Ψ indeed |V ihV | | ih | | i | i X reproduces correctly the Renyi´ entropies of CFT ground state 1 n n = gp p⊗ p∗⊗ (18) with the correct n dependence, in the regime `P L`. Cn,p ∈ H ⊗ H gp Sym Before we come to the discussion of the Renyi´ entropy. ∈ n There is a remark about the Regge calculus that we employ where dU is the Haar measure on the group of all unitary Λ p P above: The Regge calculus with cosmological constant has g transformations. gp Symn sums over all permutations p act- 2 formulations, depending on the curvature of tetrahedra in n ∈ P ing on p⊗ . The overall constant Cn,p = g Sym trgp = Regge geometries [56, 58, 59]. One may use flat tetrahedra H p∈ n (dim p + n 1)!/(dim p 1)!. in the same way as the Regge calculus without Λ, in order to H − H − Inserting this result in trρn , the average trρn becomes a sum have simpler Regge geometries. But the disadvantage is that A A over all permutations g at all nodes p, where each term as- the constant curvature space is an approximation rather than p sociates to a choice of{ g} at each p. It is straight-forward to exact solution. From this perspective, it is more convenient p compute that for large bond dimension D 1, the sum over for us to employ the second formulation with constant cur- g is dominated by the contribution fromg satisfying the vature tetrahedra, where the constant curvature is consistent p p following{ } boundary condition: g = I if p{ is} adjacent to A¯, with the cosmological constant Λ. Consequently, the Regge p while g = ( (n)) 1 if p is adjacent to A [25]. equation of motion Eq.(14) gives the smooth AdS geometry p − 3 Given g Csatisfying the boundary condition, g contains as the exact solution, although we start with a discrete gravity p p different domains{ } on Σ with different permutations.{ We} denote theory. by Rg the closed region in which p Rg are of constant gp = g. R R denotes the domain∈ wall shared by R , R with g∩ g0 ≡ Sg,g0 g g0 III. RENYI´ ENTROPIES two different permutations g , g0. Recall that the Wheeler-deWitt wave functions Φ(~a) is a path integral on a triangulated manifold M, where ~a is the We compute Renyi´ entropies S of the state Ψ by specify- n boundary condition. The product of 2n copies of Φ(~a) or ing a boundary region A ∂Σ which contains a| subseti of open Φ(~a) in tr(ρn ) is the path integral on the product of 2n copies links. Recall that Ψ is made⊂ by random tensors at nodes, the ∗ A of manifolds M or M¯ (FIG.3) with identical triangulations. n-th Renyi´ entropy| isi given by an average over random tensors The result of random average Eq.(18) permutes the boundary [25] conditions ~a and identifies boundary conditions ~a. The trace n in Eq.(17) effectively glue the 2n path integrals. As a result, 1 tr(ρA) S n(A) = ln . (15) the sum over gp becomes a sum of path integrals on different 1 n (trρ )n { } − A manifolds. Each gp determines a gluing of 2n copies of M ¯ { } The fluctuation away from the average is discussed in Ap- and M in certain manner. The path integral is defined on the = pendix A. ρ is the reduced density matrix by tracing out the resulting manifold. See FIG.4 for an example of n 2. A A graphic computation of tr(ρn ) is presented in Appendix B. degrees of freedom located in the complement A¯ = ∂Σ A. A tr(ρn ) can be conveniently written in terms of the pure density\ Here is a description of the manifold determined by each A g : Inside a domain R , the i-th copy of M¯ is glued to the matrix ρ = Ψ Ψ { p} g | ih | g(i)-th copy of M. In a neighboring domain Rg0 , the i-th copy h (n)i ¯ tr(ρn ) = tr (ρ ρ) , (16) of M is glued to the g0(i)-th copy of M. Therefore the gluing A ⊗ · · · ⊗ CA in Rg and Rg0 results in a branch cut, where the domain wall where the trace is taken in n copies of boundary Hilbert g,g0 is a 1d branch curve containing all branch points. Tak- (n) ingS all domains with different permutations into account, each space. A cyclicly permutes the states of region A, leaving C (2) (1) (1) gp determines a manifold M g made by gluing n copies of the states of A¯ invariant, e.g. for n = 2, ( µ µ ¯ p A ` A A { } ¯ { } (2) (2) (2) (1) (1) (2) C | i | i ⊗ M and n copies of M, with a number of branch cuts. FIG.4 µ µ ¯) = µ µ ¯ µ µ ¯ where µ forms a | iA| iA | ` iA| iA ⊗ | iA| iA | i illustrates a simple situation with n = 2, where there are 2 basis in the boundary Hilbert space. domains of the identity I and cyclic (2). The domain wall is Define the pure state density matrix ρP = E~a,Φ E~a,Φ where a branch curve with Z symmetry. ThisC situation generalizes P ~ ρn | ih | 2 E~a,Φ = ~a Φ(a) p,p0 app0 . tr( A) can be written as 1 | i ⊗ | i to a general domain wall g,g0 and a cycle c in g− g0. Indeed, S 1 n each domain wall g,g0 becomes χ(g− g0) branching curves n n ⊗ (n) tr(ρ ) = tr ρ⊗ p p p . (17) S1 1 A P A g,g0 (c), where χ(g− g0) is the number of cycles in g− g0. Each ⊗ |V ihV | C S 1 branch curve g,g0 (c) associates to a cycle c C(g− g0), and where the trace is taken in all 3 at all nodes. tr(ρn ) has a local Z Ssymmetry, where n is the length∈ of the cycle c p ⊗ A nc c P contains 2n copies of wave functionsH ≡ HΦ(~a) or Φ(~a) . (the number of involved elements i c), satisfying n = n. ∗ ∈ c c 6 ˆ At each branch curve g,g0 (c) where 2nc copies of M and tetrahedra and keeping ε` = 0 fixed in the bulk ∂nc V Nnc = ¯ S 2 M meet, Regge actions in the 2nc copies of Φ(~a), Φ(~a)∗ con- 2πLAdS 2 L g,g (c)(nc). Integrating this relation gives tribute boundary terms L Θ to each ` (c). In ad- nc S 0 ` ` ⊂ Sg,g0 dition, thanks to a normalization factor which contributes Z n 2 L` c 2πL (1 nc) ln d[ap,p (c)] = (1 nc) (see Appendix B), the total AdS 0 4`P V Nˆ n = V (N1) + L (c)(q) dq, (22) − − c q2 Sg,g0 contributions at each g,g0 (c) precisely make a new bulk term 1 1 P S ` L`ε` of Regge action. On the other hand, it is easy 8π`P where N has no singularity at (c). ⊂S ¯ 1 g,g0 to see that when M glues to M inside a domain Rg, a pair of S ˆ Because the geometry is smooth AdS3 on Nnc away from boundary terms from S Regge(M) and S Regge(M¯ ) again makes a g,g (c), to compute L (c)(q), we use the metric on Nˆ q in the 0 g,g0 bulk term of Regge action on the glued manifold [57]. hyperbolicS foliation [13,S 14]: n As a result, tr(ρA) is written as a sum of discrete path inte-  2  2 grals of Regge actions on different manifolds M gp :  r 1  dr { } ds2 =   L2 dτ2 + + r2du2. (23) X  L2 − q2  AdS r2 1 n S Regge(M gp ) AdS L2 q2 tr(ρ ) e− { } . (19) AdS − A ∼ L { `} The periodicity of τ is τ τ + 2π. r satisfies r L /q. R ∼ ≥ AdS It allows us to translate the Renyi´ entropy of boundary CFT to du is the geodesics length in the hyperbolic plane with unit the bulk geometry. curvature. g,g0 (c) is located at the origin r = LAdS /q. Rewrit- S 1 2 As is mentioned above, we consider the regime `P L`. ing the metric by (2qLAdS )− ξ = r LAdS /q and consider- ξ − On each M gp , the dominant contribution comes from the so- ing the limit 0 shows the length of g,g0 (c) is given by { } 1→ R 1 S lution of equation of motion L (c)(q) = q− LAdS du q− LAdS l , where l is g,g0 (c) g,g0 g,g0 S Sg,g0 ≡ S S the geodesic length of g,g0 evaluated in the hyperbolic plane ε` = 0, ` bulk(M gp ) (20) S ∀ ⊂ { } with unit curvature. l , is independent of c. As a result, Sg g0 It implies that as the leading contribution, the geometry on Z nc 2π 1 n2 M g is smooth AdS everywhere. The on-shell action gives c p 3 L (c)(q) dq = − πLAdS l (24) { } 2 g,g0 2 g,g0 − 1 q S nc S V M gp 1 ( { }) X 8π` 2 +boundary terms S Regge(M gp ) − P L e− { } e AdS (21) The volume of N : V(N ) = n V Nˆ , i.e. n times ∼ c nc c nc c L` { } Eq.(22). When we glue back Nnc in M gp . The first term { } in Eq.(22) gives ncV (N1), and effectively replaces Nnc by nc Focus on a given M gp , locally at each g,g0 (c) of a given { } S copies of N1, which resolves the branch curve g,g0 (c) in M gp . cycle c, the geometry has a local Znc symmetry at both contin- S { } When all branch curves are resolved, M g reduces to n copies uum level and discrete level, because we use the same triangu- { p} M . Therefore when we sum all domain walls and all cycles, lation on all copies of M and M¯ . We cut a local neighborhood 1 Z Nnc at g,g0 (c) from M gp , and make a nc quotient. The orb- X X 2 S { } 3 1 nc Nˆ = N /Z Nˆ V M gp = nV (M1) πL − l g,g ifold is denoted by nc nc nc . The geometry on nc has a { } AdS S 0 − 1 nc 1 , c C(g g ) conical singularity at g,g (c) with deficit angle 2π 1 . In Sg g0 ∈ − 0 S 0 − nc the language of [15], the geometry we derive is back-reacted nc 1 It is shown in Appendix C that the sum is dominant by gp by a cosmic brane with tension Tnc = −` , located at g,g0 (c). { } 4nc P S with only a single domain wall separating I in RA¯ and We may analytic continue nc by considering arbitrary conical (n) 1 S ( )− in RA, where RA (or RA¯) is the region bounded by the singularity or brane tension. boundaryC region A (or A¯) and the domain wall (FIG.3). We The geometry of branch curve g,g (c) is determined by the 0 denote the corresponding M g by Mn S { p} equation of motion as in [12]. Both on-shell geometries on Nnc ˆ ˆ ˆ and Nnc are AdS3, except the conical singularity of Nnc . Nnc 2 3 1 n is a fundamental domain in Nn of Zn. Nˆ n may be obtained V (Mn) = nV (M1) πL − l (25) c c − AdS n S by cutting Nnc into nc identical pieces, pick up one piece, fol- ˆ lowed by identifying its 2 cut boundaries. Nnc is AdS3 away The contribution of any other M g is much less than Eq.(25), { p} from the singularity. So the glued boundaries can be chosen to with the gap of order L`/`P = 4 ln d[a] 1. be identical hyperbolic surfaces intersecting at the singularity. n As a result, the dominant contribution of tr(ρA) in Eq.(19) The singularity (c) as the intersection has to be a geodesic Sg,g0 is given by (hyperbola) in the hyperbolic plane. The length L (c)(nc) of g,g0 S 1 V(Mn) +boundary terms g,g (c) explicitly depends on nc. Since the triangulation of M 8π`P L2 0 n − AdS S tr(ρA) e (26) has been fixed, we only consider g,g0 (c) made by the edges ∼ in the triangulation. Otherwise theS equation of motion cannot The denominator tr(ρ )n in Eq.(15) can be computed in a be satisfied and the domain wall (c) doesn’t give leading A g,g0 very similar manner, since order contribution. S ˆ Consider the volume of Nnc . We analytic continue nc and n n n tr(ρ ) = tr ρ⊗ ⊗ . (27) compute the derivative. By Schlafli¨ identity of hyperbolic A P ⊗p |VpihVp| 7

(n) which different from Eq.(17) by removing A in the trace. IV. CONCLUSION AND OUTLOOK We still use the Haar random average Eq.(18),C and write tr(ρ )n as a sum over all permutations g at all nodes. How- A { p} In this paper we explicitly construct a boundary state Ψ as ever because (n) is absent, as D 1 the dominant config- a linear combination of random tensor networks to understand| i CA urations of gp satisfy the boundary condition that gp = I at the holographic duality at the level of many-body quantum { } the entire boundary. Thus suppose gp has different domains states. Ψ is shown to reproduce holographically the (1+1)d { } | i with different gp, domain walls are detached from the bound- CFT Renyi´ entropies S n(A) with the correct n-dependence, ary, and contain closed curves. which agrees with the known behavior of CFT ground states. Using the same argument as the above, we can write tr(ρ )n The boundary state Ψ is given by the holographic mapping A | i as a sum of path integral of Regge action on different M g , applied to the bulk discretized Wheeler-deWitt wave function { p} similar to Eq.(19). domain walls become the branch curves in Φ(~a) as a path integral of (2+1)d gravity. The holographic M g , which contain closed curves. However, since the inter- CFT Renyi´ entropies are derived by the average of random { p} section of two hyperbolic surfaces cannot give closed branch tensors in the tensor networks. The random tensors effec- curves, M g with closed branch curves doesn’t admit AdS3 tively glue the path integrals and relate boundary Renyi´ en- { p} geometry. Thus the equation of motion doesn’t have any so- tropies to partition functions of bulk gravity on branched cover lution, except gp = I identically without any domain wall. As spacetimes. The relation between boundary Renyi´ entropies n a result tr(ρA) is dominant at the configuration that gp = I and bulk gravity partition functions on branched cover space- everywhere times has been a key assumption in the existing holographic Renyi´ entropy computations [11, 12, 15], while our work de- V(M ) n 1 +boundary terms 8π`P L2 rives this relation with random tensor networks in the con- n − AdS tr(ρA) e , (28) ∼ text of AdS3/CFT2. The derivation also gives the duality be- where the boundary terms are identical to the ones appearing tween the boundary Renyi´ entropy and bulk cosmic brane (in in Eq.(26). AdS3/CFT2), which has been proposed as a conjecture in [15]. We find that the average Renyi´ entropy is given by As another key feature, our derivation doesn’t assume the bulk geometry corresponding to the entropy calculation. The 1 geometries on branched cover manifolds are emergent from S n(A) [ln Z(n) n ln Z(1) ] (29) ' 1 n ∞ − ∞ the boundary CFT state Ψ via the large bond dimension limit − of tensor networks. Ψ |isi proposed as a new way to approxi- where mate the ground state| ofi large-N CFT. 1 V (Mn) The results in this paper suggest that Ψ is a promising real- ln Z(n) (30) | i ∞ π` 2 ization of the holographic duality in the tensor network, which ≡ −8 P LAdS deserves to be further developed. It is interesting to generalize is the on-shell action of Einstein gravity on 3-manifold Mn. the derivation to higher dimensional CFTs, which is expected The relation Eq.(29) has been an assumption in the existing to derive the conjectured duality between the boundary Renyi´ derivation of RT formula from AdS/CFT [11, 12, 15]. But entropy and bulk cosmic brane in general context. A current it is now derived from the state Eq.(9) using random tensor research undergoing is to understand the conformal symmetry networks. of Ψ . Ψ proposed as the CFT ground state is expected to Inserting Eq.(25) in Eqs (29) and (30), we obtain the RT be| (approximately)i | i invariant under boundary conformal trans- formula of Renyi´ entropy for CFT2, which has the desired n formations. The bulk Wheeler-deWitt wave function and its dependence. semiclassical behavior turns out to be a nice tool to understand ! the symmetry perspectives. It is also important to compare Ψ 1 LAdS | i S n(A) 1 + l (31) with the optimized MERA, given that both approaches cap- ' n 8`P S tures the behavior of CFT ground state. The comparison may require to extend our analysis to the regime of finite bond di- where LAdS l corresponds to Armin the geodesic length in S mensions (the finite-N CFT). AdS3 in Eq.(1). The usual RT formula is recovered as n 1. The above result reproduces the Renyi entropy computed→ by The holographic mapping with random tensor networks Hung-Myers-Smolkin-Yale in [14] using the AdS/CFT as- should be applied to study the excited states in the boundary sumptions. To see it is indeed the right Renyi´ entropy of the CFT. The states excited by certain boundary operators acting on Ψ should be the image of the holographic mapping from boundary CFT2, recall the central charge of CFT relates to | i 3LAdS a bulk wave function, created by a bulk operator acting on LAdS and `P by c = , and l relates to length lA of bound- 2`P S the present Wheeler-deWitt wave function Φ(~a). It suggests ary interval A by l 2 ln(lA/δ) in Poincare´ patch, where δ is a UV cut-off. It givesS ' a holographic relation between boundary and bulk operators and is a promising starting point to understand the bulk recon- ! ! 1 c l struction in the AdS/CFT [61, 62]. Since Φ(~a) are the coef- S (A) 1 + ln A (32) n ' n 6 δ ficients in the expansion of Ψ with random tensor networks, the relation between boundary| i and bulk operators should be which matches precisely the Renyi´ entropy Eq.(3) of CFT2 an analog of the standard duality in quantum mechanics be- with correct n dependence [51]. tween operators acting on Dirac kets and operators acting on 8 wave functions. copies of Mn. As a result, This paper for the first time relates Regge calculus, as a 2 " !# Z(n) Y C ,p ` C ,p popular tool in general relativity, to the tensor network, which = 2n + P , 2n 2 2 1 O 2 1 (A2) is extensively studied in quantum information, condensed Z(n) p Cn,p L` Cn,p ≤ matter physics, and the AdS/CFT. The framework developed ∞ in this paper is a promising platform for communications be- which implies the following bound tween different fields. Regge calculus closely relates to loop !2 ! Z(n) ` quantum gravity, especially its covariant formulation [59, 63– 1 O P (A3) 67]. This work is an interesting starting point to explore the Z(n) − ≤ L` relation between loop quantum gravity, AdS/CFT correspon- ∞ Bounding the fluctuation of Z(n) by ε/4 has the following dence, and string theory. probability by Markov inequality,:

2 ! Z(n) 1 ! Acknowledgments Z(n) ε Z(n) `P ∞ − . Prob 1 2 O 2 (A4) Z(n) − ≥ 4 ≤ ε ≤ ε L` ∞ 4 MH acknowledges Ling-Yan Hung and Yidun Wan at Fu- n dan University in Shanghai, Yong-Shi Wu at University of Similar conclusion can be drawn for Z(1) . Bounds on the n Utah, Bei Zeng at University of Waterloo, and Jie Zhou at fluctuations of Z(n), Z(1) implies the bound on the fluctuation Perimeter Institute for Theoretical Physics for the hospitality of S n(A). The probability of violating the following bound is 2 during his visits, stimulating discussions, and collaborations. of O `P/ε L` MH acknowledges the useful discussions with Xiao-liang Qi n ! and Zhao Yang, and especially thanks Ling-Yan Hung for ex- 1 Z(n) Z(1) S n(A) S n(A) ln + ln plaining the result in [14]. MH also acknowledges support − ≤ n 1 Z(n) Z(1)n − ∞ ∞ from the US National Science Foundation through grant PHY- ε (A5) 1602867, and the Start-up Grant at Florida Atlantic Univer- ≤ sity, USA. where we have used that ln(1 ε/4) ε/2 for small ε. When 2 | ± | ≤ `P/L` ε , the above bound of fluctuation is satisfied with a high probability 1 O ` /ε2L . − P ` Appendix A: Bound on Fluctuation

Appendix B: Compute tr(ρn ) In this appendix we examine the fluctuation of the Renyi´ en- A tropy S n(A) from the average value S n(A) in Eq.(31), to qual- ify how well is the approximation. We show that in the regime This appendix demonstrates a graphic computation of the random average tr(ρn ). It is sufficient to demonstrate the case `P L` the fluctuation is generically small. The method used A in the following is similar to [25]. of n = 2, while the generalization to n > 2 is obvious. n FIG.5 lists the basic ingredients of presenting graphically We denotes by Z(n) = tr(ρA) and Z(n) the average value ∞ 2 of Z(n) as `P L` (same as in Eq.(30)). We consider the the computation of tr(ρA). We have suppress 1 spatial dimen- following fluctuation of Z(n): sion in the presentation where the tensor network ~a is illustrated by a chain, and two ends of the chain illustrate| i regions !2     A and A¯. Z(n)  Z(n)2   Z(n)  1 =  1 2  1 FIG.6 is a graphic illustration of tr(ρn ) in Eqs.(16) and (17)  2    A Z(n) − Z(n) − − Z(n) − = ρn Φ ~ ∞ ∞ ∞ at n 2. tr( A) contains n copies of ∗(a) and n copies of   Φ(~a).  Z(n)2   1 (A1) The random average Eq.(18) is performed for 2n copies of ≤ Z(n)2 −  ∞ random tensors locally at each node p. Applying Eq.(18) to n tr(ρA) gives a sum over choices of permutations gp at all p. Z(n) Z(n) because in the approximation we made as `P { } ≥ ∞ ρ2 L , the neglected terms in the sums are all non-negative. The computation of tr( A) is illustrated in FIG.7, which shows ` = 2 a typical term corresponding to two domains with gp I Z(n) is computed in a similar way as the above, using the (2) and gp0 = . The second step in FIG.7 uses the relations random average formula Eq.(18), changing n by 2n. It leads in FIG.5 andC FIG.8. The red lines illustrate the identifica- to that the dominant contribution of Z(n)2 is again given by a tion and summation of the boundary data ~a in path integrals sum over permutations gp at all p, whose boundary condition Φ(~a), Φ(~a)∗, which glues the path integrals. The resulting path is g = I in A¯ and g {= }(1 n)(n + 1 2n) in A. Hence p p ··· ··· integral is on the branch cover manifold in FIG.4, where the 2 Z(n) is also written as a sum of path integrals on different branch points are in the domain wall . Each link (p, p0) cross- S 1 M gp . The situation of a single domain wall separating I in ing associates a normalization factor d[a]− , generalized to { } S1 nc RA¯ and (1 n)(n + 1 2n) in RA gives again the dominant d[a] − for an arbitrary pair gp, gp0 , where nc is the length of ··· ··· 1 contribution. The 3-manifold M gp in this case is simply 2 the cycle c in gp− gp . { } 0 9 a a

⇤⇤ = =

ap,p = ap,p0 bp,p0 = | 0 i p p0 h | i ab a tr(⇢2 ) = A¯ I I (2) A

b b A p C p0 +

··· p A¯ |Vpi |V 0 i A ~a = | i p p0 a1 a2 a3 a4 ⇤⇤ A¯ A 1

= = A¯ A + | i d[a] S

···

FIG. 5: Graphic presentation of the link state a , the inner prod- | p,p0 i uct ap,p bp,p , the tensor network ~a , and Ψ . Blue triangles in ~a ρ2 h 0 | 0 i | i | i | i FIG. 7: The graphical computation of the average tr( A). are random tensors located at nodes p. Connecting the tensor |Vpi Φ ~ ~ a a network to (a) (gray box) with red lines illustrate the sum over a in b b deﬁning Ψ . The gray box illustrates Φ(~a) as the path integral over | i 1 1 M. = = d[a] S d[a] ab bc cd c c d d ⇤⇤

2 tr(⇢A) = A¯ p p0 A FIG. 8: The graphical computation involving inner products among

4 link states ap,p0 , bp,p0 , cp,p0 , dp,p0 at a domain wall . The in-

| i | i | i | i S volved link data a, b, c, d are all identiﬁed at , which makes a S S branch point (curve) of the branch cover manifold. The inner prod- 1 ucts give a prefactor d[a]− .

FIG. 6: Graphic illustration of tr(ρ2 ). A equality” (see Appendix D for a proof) X 1 n2 X 1 n2 X 1 n2 − c − c + − c . (C3) Appendix C: Domain Wall Contributions nc ≥ nc nc c C(g 1g ) c C(g 1g ) c C(g 1g ) ∈ 1− 3 ∈ 1− 2 ∈ 2− 3 It implies each trivalent intersection gives the following con- In this section, we prove that the configuration g with a p tribution single domain wall indeed gives the leading contribution{ } to P 2 2 2 X 1 nc X 1 nc X 1 nc gp . Let’s consider a more generic case shown in FIG.9(a), { } l 13 − + l 12 − + l 23 − where more than one domain-walls are created in the bulk of S nc S nc S nc c C(g 1g ) c C(g 1g ) c C(g 1g ) Σ. We are going to show that this configuration always con- ∈ 1− 3 ∈ 1− 2 ∈ 2− 3 2 tribute less than a single domain-wall. X 1 nc − l 13 + min l 12 , l 23 Recall the computation of V(M ) and the discussion be- ≤ nc S S S gp c C(g 1g ) low Eq.(24). When we sum all domain{ } walls and all cycles, ∈ 1− 3 which is less than a single domain wall contribution. 2 For any intersection with e.g. 4 domain walls, one can al- V M gp nV (M ) X X 1 n L − { } = − 1 + − c AdS l (C1). ways shift the end point of one domain wall away from the 2 2 n 8` Sg,g0 8π`PLAdS 8π`PLAdS 1 c P intersection and obtain a smaller l (greater contribution to g,g0 c C(g− g0) S ∈ Eq.(C1)). It reduces the 4-valent intersectionS back to trivalent Thus each domain-wall carries the contribution proportional situation, which implies the contribution after the above shift to is still smaller than the single domain wall configuration. The same argument applies to the intersection with larger number 2 of domain walls. X 1 nc l g,g − (C2) Therefore we find that the contribution of the multi- S 0 nc c C(g 1g ) ∈ − 0 domain-wall configuration is less or equal to the single domain-wall configuration in Eq.(C1) where l is the geodesic length on the hyper- S X 2 2 bolic plane with unit curvature. Given trivalent intersection 1 nc 1 n l g,g − − l . (C4) S 0 nc ≤ n S of domain-walls, which separate three domains with permu- c C(g 1g ) ∈ − 0 tations g1, g2, g3 (FIG.10), we have the following “triangle in- 10

Let f be a generator of a norm on Symn. For every x ,... x 1, we have 1 k ≥ k k k Pk k X X f (x ) X f ( xi) X f (x ) = x i x i=1 = f ( x ). i i x ≤ i Pk i i=1 i=1 i i=1 i=1 xi i=1

Let f be a generator of a norm on Symn. For every Pk Pk x1,... xk 1, we have i=1 f (xi) f ( i=1 xi k + 1). In- deed, when≥ k = 2, we have ≥ −

R x1 R x2 f (x1) + f (x2) = f 0(x)dx + f 0(x)dx FIG. 9: shows the space Σ with boundary ∂Σ divided into regions A 1 1 R x1 R x1+x2 1 and A¯. Σ contains the domain-walls (1), (2), , (8), which divide f 0(x)dx + − f 0(x)dx = f (x1 + x2 1). ··· 1 x1 the bulk of Σ into regions I, II, , VI. Each bulk region associates a ≥ − ··· (n) 1 permutation gI,II, ,VI , with gI = I and gII = ( )− . ··· C Suppose the argument holds when k = k0. When k = k0 + 1, we have

kX0+1 Xk0 kX0+1 f (x ) f ( x k0 + 1) + f (x ) f ( x k0). i ≥ i − k0+1 ≥ i − i=1 i=1 i=1

Let d : Sym R be a norm on Sym whose generator is n → n f , g Symn be a permutation, c Symn be a cycle. Then we have∈d(cg) d(c) + d(g). This statement∈ can be proven by the following: ≤ FIG. 10: A trivalent intersection of domain walls has less contribu- Let A C(g) be the set of all cycles in C(g) that are disjoint tion than a single domain wall. with c, Then⊆ A C(cg). Let B = C(g) A, B = C(cg) A. ⊆ 1 \P 2 P \ We obtain d(c) + d(g) d(cg) = f (nc) + f (nr) f (nr). − r B − r B P P P∈ 1 ∈ 2 Then N = r B1 nc = r B2 nc leads to f (nr) f (N). By ∈ ∈ r B ≤ Appendix D: Proof of triangle inequality (C3) P ∈ 2 the fact that B1 nc, we have f (nr) f (N B1 + 1) | | ≤ r B ≥ − | | ≥ ∈ 1 For every permutation g Symn (in the following discus- f (N nc + 1). Therefore sion, we assume that n 2),∈ we can decompose g as − ≥ d(c) + d(g) d(cg) f (nc) + f (N nc + 1) f (N) 0. Yk − ≥ − − ≥ g = ci, Let d : Symn R be a norm on Symn whose generator is i=1 f . For g , g S→, we have d(g g ) d(g ) + d(g ). Indeed, 1 2 ∈ n 1 2 ≤ 1 2 Pk let C(g1) = c1,..., ck . We have where ci Symn are disjoint cycles such that i=1 nci = n. { } ∈ = ,..., Denote C(g) c1 ck as the set of disjoint cycles whose k k g { } Y Y product is . d(g g ) = d(( c )g ) d(c ) + d(( c )g ) Let d : Sym R be a function which satisﬁes that there 1 2 i 2 ≤ 1 i 2 n i=1 i=2 exists a function→f C2[1, + ) such that (i) f (x) 0 for ∈ ∞ 00 ≤ k f (x) 0 X x 1; (ii) x 0 for x 1; (iii) For each permutation ... d(c ) + d(g ) = d(g ) + d(g ). ≥ ≥ P≥ ≤ ≤ i 2 1 2 g Symn, we have d(g) = c C(g) f (nc). We say d is a norm i=1 ∈ ∈ on Symn, and f is the generator of d. 2 Let f be a generator of a norm d on Sym , then f (1) = 0. x 1 f (x) 0 2 n f (x) = − is a generator, since = 3 0 when Indeed let g = (1)(2 ... n), then d(g) = f (1) + f (n 1) = x x x ≥ x 1, and f (x) = 2 0. + = − 00 x3 f (1) d(g), hence f (1) 0. ≥ − ≤ 2 P 1 nc We know f 0(x) 0, because leting g(x) = f (x)/x, we have As a result the function d : Sym R, d(g) = c C g − R x n ( ) nc f (x) = [x g(x)] =≥ x g (x) + g(x) g(x) = g (x) 0. is a norm on Sym . → ∈ 0 · 0 · 0 ≥ 1 0 ≥ n

[1] J. C. Bridgeman and C. T. Chubb (2016), 1603.03039. [3] A. J. Ferris and D. Poulin, Phys. Rev. Lett. 113, 030501 (2014). [2] G. Evenbly and G. Vidal, Physical Review Letters 115, 180405 [4] E. Miles Stoudenmire and D. J. Schwab, ArXiv e-prints (2016), (2015), 1412.0732. 1605.05775. 11

[5] A. Novikov, D. Podoprikhin, A. Osokin, and D. P. Vetrov, [39] J. Cotler, P. Hayden, G. Salton, B. Swingle, and M. Walter CoRR abs/1509.06569 (2015), URL http://arxiv.org/ (2017), 1704.05839. abs/1509.06569. [40] G. Chirco, D. Oriti, and M. Zhang (2017), 1701.01383. [6] R. Orus, Eur. Phys. J. B87, 280 (2014), 1407.6552. [41] F. Pastawski and J. Preskill (2016), 1612.00017. [7] B. Swingle, Phys. Rev. D86, 065007 (2012), 0905.1317. [42] A. Peach and S. F. Ross, Class. Quant. Grav. 34, 105011 (2017), [8] M. Van Raamsdonk, in Proceedings, Theoretical Advanced 1702.05984. Study Institute in Elementary Particle Physics: New Fron- [43] D. A. Roberts and B. Yoshida, JHEP 04, 121 (2017), tiers in Fields and Strings (TASI 2015): Boulder, CO, USA, 1610.04903. June 1-26, 2015 (2017), pp. 297–351, 1609.00026, URL [44] P. Hosur, X.-L. Qi, D. A. Roberts, and B. Yoshida, JHEP 02, https://inspirehep.net/record/1484863/files/ 004 (2016), 1511.04021. arXiv:1609.00026.pdf. [45] L. Smolin (2016), 1608.02932. [9] E. P. Verlinde (2016), 1611.02269. [46] M. Han, W. Huang, and Y. Ma, Int.J.Mod.Phys. D16, 1397 [10] S. Ryu and T. Takayanagi, Phys. Rev. Lett. 96, 181602 (2006), (2007), gr-qc/0509064. hep-th/0603001. [47] A. Ashtekar and J. Lewandowski, Class.Quant.Grav. 21, R53 [11] A. Lewkowycz and J. Maldacena, JHEP 08, 090 (2013), (2004), gr-qc/0404018. 1304.4926. [48] T. Thiemann, Modern Canonical Quantum General Relativity [12] X. Dong, A. Lewkowycz, and M. Rangamani, JHEP 11, 028 (Cambridge University Press, 2007). (2016), 1607.07506. [49] C. Rovelli and F. Vidotto, Covariant Loop Quantum [13] H. Casini, M. Huerta, and R. C. Myers, JHEP 05, 036 (2011), Gravity: An Elementary Introduction to Quantum Grav- 1102.0440. ity and Spinfoam Theory, Cambridge Monographs on [14] L.-Y. Hung, R. C. Myers, M. Smolkin, and A. Yale, JHEP 12, Mathematical Physics (Cambridge University Press, 2014), 047 (2011), 1110.1084. ISBN 9781107069626, URL https://books.google.com/ [15] X. Dong, Nature Commun. 7, 12472 (2016), 1601.06788. books?id=4VjeBAAAQBAJ. [16] M. Rangamani and T. Takayanagi (2016), 1609.01287. [50] M. Han and L.-Y. Hung, Phys. Rev. D95, 024011 (2017), [17] J. de Boer and J. I. Jottar, JHEP 04, 089 (2014), 1306.4347. 1610.02134. [18] B. Chen and J.-q. Wu, JHEP 09, 015 (2016), 1605.06753. [51] P. Calabrese and J. L. Cardy, J. Stat. Mech. 0406, P06002 [19] X. Dong, Phys. Rev. Lett. 116, 251602 (2016), 1602.08493. (2004), hep-th/0405152. [20] W. Song, Q. Wen, and J. Xu, JHEP 02, 067 (2017), 1610.00727. [52] J. Barrett and P. Parker, Journal of Approximation Theory 76, [21] M. Ammon, A. Castro, and N. Iqbal, JHEP 10, 110 (2013), 107 (1994), ISSN 0021-9045. 1306.4338. [53] T. Regge, Nuovo Cim. 19, 558 (1961). [22] Y. Ling, P. Liu, C. Niu, J.-P. Wu, and Z.-Y. Xian, JHEP 04, 114 [54] R. Friedberg and T. Lee, Nucl.Phys. B242, 145 (1984). (2016), 1502.03661. [55] S. S. Gubser, M. Heydeman, C. Jepsen, M. Marcolli, S. Parikh, [23] Q. Hu and G. Vidal (2017), 1703.04798. I. Saberi, B. Stoica, and B. Trundy (2016), 1612.09580. [24] F. Pastawski, B. Yoshida, D. Harlow, and J. Preskill, JHEP 06, [56] B. Bahr and B. Dittrich, Phys.Rev. D80, 124030 (2009), 149 (2015), 1503.06237. 0907.4323. [25] P. Hayden, S. Nezami, X.-L. Qi, N. Thomas, M. Walter, and [57] J. Hartle and R. Sorkin, Gen.Rel.Grav. 13, 541 (1981). Z. Yang, JHEP 11, 009 (2016), 1601.01694. [58] B. Bahr and B. Dittrich, New J.Phys. 12, 033010 (2010), [26] M. Miyaji, T. Takayanagi, and K. Watanabe, Phys. Rev. D95, 0907.4325. 066004 (2017), 1609.04645. [59] H. M. Haggard, M. Han, W. Kaminski, and A. Riello, Nucl. [27] A. Bhattacharyya, Z.-S. Gao, L.-Y. Hung, and S.-N. Liu, JHEP Phys. B900, 1 (2015), 1412.7546. 08, 086 (2016), 1606.00621. [60] A. W. Harrow (2013), 1308.6595. [28] B. Czech, P. H. Nguyen, and S. Swaminathan, JHEP 03, 090 [61] A. Hamilton, D. N. Kabat, G. Lifschytz, and D. A. Lowe, Phys. (2017), 1612.05698. Rev. D74, 066009 (2006), hep-th/0606141. [29] Z.-X. Luo, E. Lake, and Y.-S. Wu, ArXiv e-prints (2016), [62] D. L. Jafferis (2017), 1703.01519. 1611.01140. [63] J. W. Barrett, R. Dowdall, W. J. Fairbairn, F. Hellmann, and [30] K. Li, M. Han, G. Long, Y. Wan, D. Lu, B. Zeng, and R. Pereira, Class.Quant.Grav. 27, 165009 (2010), 0907.2440. R. Laflamme (2017), 1705.00365. [64] F. Conrady and L. Freidel, Phys.Rev. D78, 104023 (2008), [31] X.-L. Qi (2013), 1309.6282. 0809.2280. [32] X.-L. Qi, Z. Yang, and Y.-Z. You (2017), 1703.06533. [65] M. Han and M. Zhang, Class.Quant.Grav. 30, 165012 (2013), [33] A. Bhattacharyya, L.-Y. Hung, Y. Lei, and W. Li (2017), 1109.0499. 1703.05445. [66] M. Han (2017), 1705.09030. [34] G. Vidal, Physical Review Letters 101, 110501 (2008), quant- [67] W. Kaminski, M. Kisielowski, and H. Sahlmann (2017), ph/0610099. 1705.02862. [35] G. Vidal, Phys. Rev. Lett. 99, 220405 (2007), cond- [68] The set of ~a are shown to form an overcomplete basis in the | i mat/0512165. boundary Hilbert space if there is no restriction on the link vari- [36] G. Evenbly and G. Vidal, Journal of Statistical Physics 145, 891 ables ~a [32]. Whether ~a with ~a as tiling of the spatial slice | i (2011), 1106.1082. still form an over-complete basis is an interesting question but [37] A. May (2016), 1611.06220. doesn’t affect the present analysis. [38] Y. Li, M. Han, M. Grassl, and B. Zeng (2016), 1612.04504.