Discrete Gravity on Random Tensor Network and Holographic R\'enyi Entropy

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 $|\Psi\rangle$ using random tensor networks as the holographic mapping, applied to the Wheeler-deWitt wave function of bulk Euclidean discrete gravity in 3 dimensions. The entanglement R\'enyi entropy of $|\Psi\rangle$ is shown to holographically relate to the on-shell action of Einstein gravity on a branch cover bulk manifold. The resulting R\'enyi entropy $S_n$ of $|\Psi\rangle$ approximates with high precision the R\'enyi 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 AdS$_3$/CFT$_2$ correspondence on random tensor networks, and provide a new proposal to approximate CFT ground state.


I. INTRODUCTION
The tensor network is a quantum state of many-body system constructed by contracting tensors according to a network graph with nodes and links (FIG.1). It is originated in condense matter physics because tensor network states efficiently compute ground states of many-body quantum systems [1,2]. In addition, the tensor network has wide applications to quantum information theory by its relation to error correcting codes and quantum entanglement [3]. It recently relates to quantum machine learning [4], as well as neuroscience [5].
One of the fascinating developments of the tensor network is the recent relation to the AdS/CFT correspondence and emergent gravity program [6,7]. The AdS/CFT correspondence proposes that the quantum gravity theory on ddimensional Anti-de Sitter (AdS) spacetime is equivalent to a conformal field theory (CFT) living at the (d − 1)-dimensional boundary of AdS. It offers a dictionary between the observables of the d-dimensional bulk gravity theory and those of the (d − 1)-dimensional boundary CFT. Properties of the bulk gravity and geometry may be reconstructed or emergent from the boundary CFT, known as the emergent gravity program [8,9]. As an important ingredient of AdS/CFT, the bulk geometry relates holographically to the entanglement in the boundary CFT, via the Ryu-Takayanagi (RT) formula which identifies the entanglement entropy S EE (A) of a (d −1)dimensional boundary region A with the area Ar min of the bulk (d − 2)-dimensional minimal surface anchored to A [10][11][12][13][14][15]. G N is the Newton constant in d dimensions. S EE (A) satisfying RT formula is referred to as the holographic entanglement entropy (HEE) [16][17][18][19][20][21][22] The tensor network can be understood as a discrete version of the AdS/CFT correspondence [23][24][25][26][27][28][29][30]. The tensor network states approximates the CFT state at the boundary, while the structure of tensor network emerges an bulk dimension built by layers of tensors. Tensors in the tensor network correspond to local degrees of freedom in the bulk [31][32][33]. The architecture of the TN may be viewed as a process of real-space renormalization, such as multiscale entanglement renormalization ansatz (MERA), where the renormalization scale relates to the coordinate of the emergent dimension [2,7,34,35]. The feature of tensor network makes it an interesting tool for realizing the AdS/CFT correspondence constructively from many-body quantum states. Among many recent progress, an important result is to reproduce the HEE on the tensor network.
The optimized scale-invariant MERA tensor networks are able to approximate ground states of (discreted) CFTs e.g. the critical Ising and Potts models. It is observed that the entanglement entropy of optimized MERA is proportional to the number of links on a minimal cut through the tensor network, as a nice realization of the RT formula [6,7,36]. Moreover there has been two recent interesting exactly-solvable models that realize the RT formula on tensor networks, including [24] with tensor networks with perfect tensors, and [25] with random tensor networks. Given a boundary region A containing a number of open links of the tensor network, the entanglement entropy of the tensor network in each of the above approaches reproduces an analog of the RT formula (See e.g. [37][38][39][40][41][42][43][44][45] for some more recent developments) where Min(# cut ) is the minimal number of tensor network links cut by a surface anchored to A. D is the bond dimension (range of tensor index). Eq.(2) is in general an upper bound in the tensor network, which is saturated in the case of perfect or random tensor network. The random tensor network approach has been related to the quantum geometry in loop quantum gravity (LQG) ( [46][47][48][49] for reviews), which relates Eq.
(2) to the geometrical RT formula Eq.(1) [50]. Although the perfect and random tensor network approaches give exactly solvable models and have nice features, it is known that both approaches suffer the issue of flat entanglement spectrum. Although the entanglement (Von Neumann) entropy Eq.(2) is consistent with the RT formula, Rényi entropies S n (A) from both approaches are all identical to Eq. The tensor network state |Ψ is given by an expansion with certain basis in the Hilbert space of many-body system |Ψ = .. a i ... |a 1 , a 2 , · · · a N . The coefficients is constructed by distributing a rank-6 tensor T p at each 6-valent node p, such that each tensor index associates to a link adjacent to p. Connecting 2 nodes by a link means contracting the corresponding indices γ l . formula of Rényi entropy has a nontrivial n dependence since the CFT Rényi entropy does [15]. For instance, in any 2d CFT (CFT 2 ), the ground state has the universal Rényi entropy [51] S n (A) = 1 + 1 n where c is the central charge, l A is the length of the region A, and δ is a UV cut-off. The CFT Rényi entropy manifestly has a nontrivial n dependence. The mismatch of Rényi entropy implies that both tensor network states in [24,25] fail to approximate the CFT ground state. An idea of the reason behind it may be extracted from the AdS/CFT: In the AdS/CFT, the CFT ground state (at strong-coupling) is dual to the bulk semiclassical AdS spacetime geometry. However the tensor networks designed in [24,25] only consider the geometry of a spatial slice in AdS, without any input about time evolution. In the continuum AdS/CFT context, the correct entanglement spectrum are obtained by considering the spacetime geometry, and taking into account the bulk dynamics given by the Einstein equation [11,12,14,15], while [24,25] don't have any dynamical input. In this work, we take the above idea as a hint to improve the tensor network approach, in which we resolve the issue.
Note that the optimized MERA for the CFT can derive the correct Rényi entropies as well as the non-flat entanglement spectrum [36]. However the goal of this work is to finding analytically a generic family of tensor networks that provide the correct CFT Rényi entropies without going through any numerical optimization at all. Our results develop important analytic methods for the tensor network AdS/CFT.
In this work, we resolve the above Rényi entropy issue by having the bulk dynamics on random tensor network states. The dynamics relates tensor networks to the discrete gravity model known as Regge calculus. We construct the state |Ψ which is proposed as an approximation to the CFT ground state. As is anticipated by our proposal, |Ψ indeed reproduces correctly the RT formula and CFT ground state Rényi entropy S n with correct n dependence. Our results develop the framework of the AdS/CFT correspondence on random tensor networks, and provide a new proposal to approximate the large-N CFT ground state. As another interesting aspect of our work, we show that the random tensor network effectively relates the boundary Rényi entropy to the partition functions of bulk gravity on the branched cover spacetime. This relation has been a key assumption in the existing AdS/CFT Rényi entropy computations [11,12,15], while our work derives this relation in the context of AdS 3 /CFT 2 . We also derive the duality between the boundary Rényi entropy and bulk cosmic brane (in AdS 3 /CFT 2 ), which has been proposed as a conjecture in [15].
In this paper, we focus on 2d CFT and 3d bulk spacetime (AdS 3 /CFT 2 ) in Euclidean signature. The CFT state |Ψ is constructed by implementing bulk gravity dynamics to random tensor network states studied in [32]. The random tensor network constructed in [32] has random tensors at all nodes p, and has labels a p,p on links (p, p ). Each a p,p labels the non-maximal entangled state |a p,p on each link. The tensor network is dual to a tiling of 2d spatial slice Σ (FIG.2). The entanglement entropy of |a p,p relates to the length L of the edge intersecting (p, p ) in the tiling. Thus each random tensor network as boundary CFT state, denoted by | a , determines a set of edge lengths L in the bulk. When the tiling is a triangulation, edge lengths uniquely determines a discrete geometry on Σ which approximates the smooth geometry as the triangulation is refined [52]. We make linear combinations of | a to write CFT states as where the coefficients Φ( a) can be understood as a wave function of bulk geometry, by the relation between a and bulk edge lengths. Eq.(4) is a holographic mapping from the bulk state Φ to the boundary state Ψ. As mentioned above, the CFT ground state is expected dual to a physical state in the bulk corresponding to the semiclassical AdS spacetime geometry. But the link labels a in Φ( a) only relate to the geometry of 2d spatial slice Σ. To let Φ( a) encode spacetime geometry, we propose Φ( a) to be the Wheeler-deWitt wave function. Namely, Φ( a) is a path integral of (Euclidean) Einstein gravity on spacetime M whose boundary contains Σ (FIG.3(a)). The geometry on Σ determined by a is the boundary condition of the path integral. We construct Φ( a) to sum all possible bulk spacetime geometries satisfying the boundary condition, while in the semiclassical limit, it localizes at the classical AdS spacetime in 3d. The semiclassical limit relates to the large bond dimension of tensor network. The tensor network link variables a are the boundary data in the spatial slice Σ, while the boundary data at other boundaries (including the infinitely past/future) doesn't affect our discussion.
Since a are data of discrete geometry, Φ( a) is the discrete version of Wheeler-deWitt wave function: It is a path integral of Regge calculus. Regge calculus is a discretization of Einstein gravity by triangulating spacetime geometries [53]. The discrete spacetime geometry is given by the edge lengths in the triangulation, known as the Regge geometry. Regge geometries converge to smooth geometries as triangulations are refined [52]. Φ( a) is a sum over all (Euclidean) Regge geometries on the spacetime M, weighted by the exponentiated Einstein-Regge action [54]. The detailed explanations of Φ( a) and random tensor networks are presented in Section II.
The Rényi entropy S n (A) of |Ψ at arbitrary n ≥ 2 is computed in Section III. The computation involves averages of the random tensors at nodes p in the tensor networks [25]. It turns out that the random average effectively glues 2n copies 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 equivalent to large bond dimensions in the tensor networks, where I Bulk (M n ) is the on-shell gravity action evaluated at the bulk solution on the branch cover manifold M n (FIG.4).
The bulk solution has the Z n replica symmetry. Eq.(5) has been an assumption in AdS/CFT derivations of HEE in e.g. [11,12,14,15]. But it is now derived from |Ψ and random tensor networks. The derivation also gives the duality between the boundary Rényi entropy and bulk cosmic brane (in AdS 3 /CFT 2 ), which has been proposed as a conjecture in [15]. As a result, we show that S n (A) reproduces the RT formula for holographic Rényi entropy for 2d CFT (Hung-Myers-Smolkin-Yale formula in [14]) Here Ar min is the geodesic length in AdS 3 . The above result of S n (A) gives the Rényi entropy Eq.(3) of CFT 2 ground state with correct n dependence. Our result is valid in the regime of large CFT central charge. This work applies discrete geometry method such as Regge calculus to study tensor networks (see [55] for other application of discrete gravity in AdS/CFT). |Ψ encodes the dynamics of bulk geometries, which is given by the discrete Einstein equation. It is interesting to further understand how the bulk dynamics might relate to the dynamics of boundary CFT, and whether a boundary CFT Hamiltonian might be induced from the bulk dynamics. It is interesting to compare our proposal of CFT ground state |Ψ to the MERA approach, to extend our result to finite central charges and understand the relation between |Ψ and concrete many-body systems such as critical spin-chains. The understanding of these aspects should develop tensor network models to realize the AdS/CFT correspondence at the dynamical level. The research on these aspects is currently undergoing.

II. RANDOM TENSOR NETWORK AND WHEELER-DEWITT WAVE FUNCTION
In this work we consider trivalent random tensor network states. A tensor network is viewed as a discrete 2d spatial slice Σ of 3d bulk spacetime. It is made by a large number of trivalent random tensors |V p ∈ H ⊗3 ≡ H p at each tensor network node p. The Hilbert space H is of dimension D. We decompose H into a number of subspaces H ⊕ a V a and denote dim(V a ) ≡ d[a]. Each internal link (p, p ) of the tensor network associates with a maximal entangled state in V a ⊗ V a of certain a, where |µ p is a basis in H. It satisfies a p,p |b p,p = δ ab . A class of random tensor networks | a can be defined by the (partial) inner product between |V p at all p and |a p,p on all internal links The inner product takes place in H at each end point p or p of each link. | a is a state in the boundary Hilbert space The label a relates to the amount of entanglement on each internal link (p, p ). The entanglement entropy S (|a p,p ) of The above class of tensor network states is proposed in [32]. We consider the boundary state |Ψ as a linear combination Here we understand the trivalent tensor network to be dual to a triangulation of the spatial slice Σ (FIG.2). Namely, each node p located at the center of a triangle ∆ p in the triangulation. Each link (p, p ) intersects transversely an internal edge shared by 2 triangles ∆ p , ∆ p . Open links in tensor network intersect transversely the edges at the boundary of triangulation.
The label a is understood as the discrete geometry in the bulk of Σ [32], in the sense that the edge length L of intersecting (p, p ) is proportional to the entanglement entropy S (|a p,p ) on each link: Here P = G N is the Planck length in 3d. In this proposal, the bulk geometry is understood as emergent from the entanglement in tensor network state. The relation can be obtained from the recent proposal of understanding tensor networks as the effective theory from coarse graining quantum gravity at Planck scale [50], in which one derives that the bond dimension d[a p,p ] of tensor network | a satisfies d[a p,p ] e L /4 P . By the relation between a and bulk geometry, Eq.(9) is a boundary state by summing over all bulk spatial geometry on Σ, while Φ( a) is a wave function of bulk geometry. Eq.(9) defines a holographic mapping from the bulk states of geometry to the boundary states of CFT. We propose the following boundary state |Ψ whose bulk wave function Φ( a) (pre-image of the holographic mapping) is an Wheeler-deWitt wave function in 3d Euclidean gravity. Namely Φ( a) is a path integral of gravity on a 3d solid cylinder M, whose boundary includes Σ in addition to the boundary where CFT lives (FIG.3(a)). The geometry a is the boundary condition on Σ in the path integral. The path integral may also depend on the boundary conditions at other boundaries of M. But we make those boundary conditions implicit since they play no role in the following analysis.
Since a gives a discrete geometry with a set of edge lengths L , more precisely, Φ( a) is a discrete version of the Wheeler-deWitt wave function. Indeed, we consider a sufficiently refined triangulation of M, and impose discrete metrics on the triangulation. Namely each tetrahedron in the triangulation carries a 3d hyperbolic geometry with constant curvature −L −2 AdS . Tetrahedron edges are geodesics in the hyperbolic space, and have edge lengths L . The set of edge lengths {L } on the triangulation defines a discrete metric of Regge geometry [53,54,56]. We define Φ( a) to be a path integral of discrete gravity on the triangulation by summing over all L in the bulk of M Φ( a) := L e −S Regge (M) (11) The boundary condition at Σ is L ⊂Σ = 4 P ln d[a p,p ]. We use L instead of integration because L are assumed as discrete data, to be consistent with L ⊂Σ . S Regge (M) is the Regge action of Euclidean gravity on the triangulated 3-manifold M evaluated at the discrete metric {L }: L AdS relates to P and the central charge of CFT by L AdS = 2 3 c P . In the following, the large bond dimension limit d[a p,p ] 1 will be taken in the computation of Renyi entropies. The limit corresponds to the bulk semiclassical limit P → 0 and large central charge c 1 as both L and L AdS are kept fixed. ε is the bulk deficit angle hinged at the bulk edge . ε is a discretization of the bulk curvature. Each bulk edge is shared by a number of tetrahedra t. In each t, the dihedral angle between 2 faces joint at is denoted by θ(t, ). The deficit angle is defined by Each boundary edge is shared by 2 boundary triangles. Θ is the angle between their outward pointing normals, equivalently Θ relates to the boundary extrinsic curvature. The 1st term in S Regge is the discretization of Ricci scalar term of Einstein-Hilbert action, while the 2nd term is the discretization of Gibbons-Hawking boundary term [57]. The last term is the cosmological constant term where V(M) is the total volume of M. All quantities ε , Θ , and V(M) are determined by edge lengths L . In order to be the boundary condition of Regge geometry, L ⊂Σ = 4 P ln d[a p,p ] have to be the edge lengths of hyperbolic triangles, which triangulate Σ. L ⊂Σ have to be a discrete metric of Σ, which constrains the allowed data a entering the sum in Eq. (9).
Applying Φ( a) in Eq. (11) to the holographic mapping Eq.(9), we obtain a boundary CFT state |Ψ , and we propose the resulting |Ψ to be the ground state of the boundary large-N CFT, in the bulk semiclassical regime P L equivalent to the large bond dimension limit by Eq.(10). As P L , the path integral Φ( a) localizes at the solution of equation of motion (deriving equation of motion uses the Schläfli identity of hyperbolic tetrahedra −δV(t)/L 2 AdS = ⊂t L δθ(t, ), see e.g. [56]) which is the discretized Einstein equation in 3d. Vanishing ε everywhere means that the 3d Regge geometry is a smooth Euclidean AdS 3 . So Φ( a) is a semiclassical wave function of bulk AdS 3 geometry. The holographic mapping is expected to map the bulk semiclassical state of AdS 3 to the ground state of boundary CFT 2 .
In Section III, we check our proposal by computing the Rényi entropies S n of the state |Ψ . We show that |Ψ indeed reproduces correctly the Rényi entropies of CFT ground state with the correct n dependence, in the regime P L . Before we come to the discussion of the Rényi entropy.
There is a remark about the Regge calculus that we employ above: The Regge calculus with cosmological constant Λ has 2 formulations, depending on the curvature of tetrahedra in Regge geometries [56,58,59]. One may use flat tetrahedra in the same way as the Regge calculus without Λ, in order to have simpler Regge geometries. But the disadvantage is that the constant curvature space is an approximation rather than exact solution. From this perspective, it is more convenient for us to employ the second formulation with constant curvature tetrahedra, where the constant curvature is consistent with the cosmological constant Λ. Consequently, the Regge equation of motion Eq. (14) gives the smooth AdS 3 geometry as the exact solution, although we start with a discrete gravity theory.

III. RÉNYI ENTROPIES
We compute Rényi entropies S n of the state |Ψ by specifying a boundary region A ⊂ ∂Σ which contains a subset of open links. Recall that |Ψ is made by random tensors at nodes, the n-th Rényi entropy is given by an average over random tensors [25] The fluctuation away from the average is discussed in Appendix A. ρ A is the reduced density matrix by tracing out the degrees of freedom located in the complementĀ = ∂Σ \ A. tr(ρ n A ) can be conveniently written in terms of the pure density matrix ρ = |Ψ Ψ| where the trace is taken in n copies of boundary Hilbert space. C (n) A cyclicly permutes the states of region A, leaving the states ofĀ invariant, e.g. for n = 2, C (2) where |µ forms a basis in the boundary Hilbert space.
Define the pure state density matrix ρ P = |E a,Φ E a,Φ | where |E a,Φ = a Φ( a) ⊗ p,p |a pp . tr(ρ n A ) can be written as where the trace is taken in all H p ≡ H ⊗3 at all nodes. tr(ρ n A ) contains 2n copies of wave functions Φ( a) or Φ( a) * .
The random average trρ n A relates to average n copies of random tensors |V p V p | at each node p. Taking an arbitrary reference state |0 p ∈ H p , the random tensor |V p = U p |0 p with U p unitary transformation on H ⊗3 ≡ H p . The Haar random average is given by [25,60]: where dU p is the Haar measure on the group of all unitary transformations. g p ∈Sym n sums over all permutations g p acting on H ⊗n p . The overall constant C n,p = g p ∈Sym n trg p = Inserting this result in trρ n A , the average trρ n A becomes a sum over all permutations {g p } at all nodes p, where each term associates to a choice of g p at each p. It is straight-forward to compute that for large bond dimension D 1, the sum over {g p } is dominated by the contribution from {g p } satisfying the following boundary condition: g p = I if p is adjacent toĀ, while g p = (C (n) ) −1 if p is adjacent to A [25].
Given {g p } satisfying the boundary condition, {g p } contains different domains on Σ with different permutations. We denote by R g the closed region in which p ∈ R g are of constant g p = g. R g ∩R g ≡ S g,g denotes the domain wall shared by R g , R g with two different permutations g g .
Recall that the Wheeler-deWitt wave functions Φ( a) is a path integral on a triangulated manifold M, where a is the boundary condition. The product of 2n copies of Φ( a) or Φ( a) * in tr(ρ n A ) is the path integral on the product of 2n copies of manifolds M orM (FIG.3)  Here is a description of the manifold determined by each {g p }: Inside a domain R g , the i-th copy ofM is glued to the g(i)-th copy of M. In a neighboring domain R g , the i-th copy ofM is glued to the g (i)-th copy of M. Therefore the gluing in R g and R g results in a branch cut, where the domain wall S g,g is a 1d branch curve containing all branch points. Taking all domains with different permutations into account, each {g p } determines a manifold M {g p } made by gluing n copies of M and n copies ofM, with a number of branch cuts. FIG.4 illustrates a simple situation with n = 2, where there are 2 domains of the identity I and cyclic C (2) . The domain wall is a branch curve with Z 2 symmetry. This situation generalizes to a general domain wall S g,g and a cycle c in g −1 g . Indeed, each domain wall S g,g becomes χ(g −1 g ) branching curves S g,g (c), where χ(g −1 g ) is the number of cycles in g −1 g . Each branch curve S g,g (c) associates to a cycle c ∈ C(g −1 g ), and has a local Z n c symmetry, where n c is the length of the cycle c (the number of involved elements i ∈ c), satisfying c n c = n.
At each branch curve S g,g (c) where 2n c copies of M and M meet, Regge actions in the 2n c copies of Φ( a), Φ( a) * contribute boundary terms L Θ to each ⊂ S g,g (c). In addition, thanks to a normalization factor which contributes (1 − n c ) ln d[a p,p (c)] = (1 − n c ) L 4 P (see Appendix B), the total contributions at each S g,g (c) precisely make a new bulk term 1 8π P ⊂S L ε of Regge action. On the other hand, it is easy to see that when M glues toM inside a domain R g , a pair of boundary terms from S Regge (M) and S Regge (M) again makes a bulk term of Regge action on the glued manifold [57].
As a result, tr(ρ n A ) is written as a sum of discrete path integrals of Regge actions on different manifolds M {g p } : It allows us to translate the Rényi entropy of boundary CFT to the bulk geometry. As is mentioned above, we consider the regime P L . On each M {g p } , the dominant contribution comes from the solution of equation of motion It implies that as the leading contribution, the geometry on M {g p } is smooth AdS 3 everywhere. The on-shell action gives Focus on a given M {g p } , locally at each S g,g (c) of a given cycle c, the geometry has a local Z n c symmetry at both continuum level and discrete level, because we use the same triangulation on all copies of M andM. We cut a local neighborhood N n c at S g,g (c) from M {g p } , and make a Z n c quotient. The orbifold is denoted byN n c = N n c /Z n c . The geometry onN n c has a conical singularity at S g,g (c) with deficit angle 2π 1 − 1 n c . In the language of [15], the geometry we derive is back-reacted by a cosmic brane with tension T n c = n c −1 4n c P , located at S g,g (c). We may analytic continue n c by considering arbitrary conical singularity or brane tension.
The geometry of branch curve S g,g (c) is determined by the equation of motion as in [12]. Both on-shell geometries on N n c andN n c are AdS 3 , except the conical singularity ofN n c .N n c is a fundamental domain in N n c of Z n .N n c may be obtained by cutting N n c into n c identical pieces, pick up one piece, followed by identifying its 2 cut boundaries.N n c is AdS 3 away from the singularity. So the glued boundaries can be chosen to be identical hyperbolic surfaces intersecting at the singularity. The singularity S g,g (c) as the intersection has to be a geodesic (hyperbola) in the hyperbolic plane. The length L S g,g (c) (n c ) of S g,g (c) explicitly depends on n c . Since the triangulation of M has been fixed, we only consider S g,g (c) made by the edges in the triangulation. Otherwise the equation of motion cannot be satisfied and the domain wall S g,g (c) doesn't give leading order contribution.
Consider the volume ofN n c . We analytic continue n c and compute the derivative. By Schläfli identity of hyperbolic tetrahedra and keeping ε = 0 fixed in the bulk ∂ n c V N n c = AdS q 2 L S g,g (c) (q) dq, (22) where N 1 has no singularity at S g,g (c). Because the geometry is smooth AdS 3 onN n c away from S g,g (c), to compute L S g,g (c) (q), we use the metric onN q in the hyperbolic foliation [13,14]: The periodicity of τ is τ ∼ τ + 2π. r satisfies r ≥ L AdS /q. du is the geodesics length in the hyperbolic plane with unit curvature. S g,g (c) is located at the origin r = L AdS /q. Rewriting the metric by (2qL AdS ) −1 ξ 2 = r − L AdS /q and considering the limit ξ → 0 shows the length of S g,g (c) is given by L S g,g (c) (q) = q −1 L AdS S g,g (c) du ≡ q −1 L AdS l S g,g , where l S g,g is the geodesic length of S g,g evaluated in the hyperbolic plane with unit curvature. l S g,g is independent of c. As a result, The volume of N c : V(N n c ) = n c V N n c , i.e. n c times Eq. (22). When we glue back N n c in M {g p } . The first term in Eq. (22) gives n c V (N 1 ), and effectively replaces N n c by n c copies of N 1 , which resolves the branch curve S g,g (c) in M {g p } . When all branch curves are resolved, M {g p } reduces to n copies M 1 . Therefore when we sum all domain walls and all cycles, It is shown in Appendix C that the sum is dominant by {g p } with only a single domain wall S separating I in RĀ and (C (n) ) −1 in R A , where R A (or RĀ) is the region bounded by the boundary region A (orĀ) and the domain wall (FIG.3). We denote the corresponding The contribution of any other M {g p } is much less than Eq. (25), with the gap of order L / P = 4 ln d[a] 1. As a result, the dominant contribution of tr(ρ n A ) in Eq.(19) is given by The denominator tr(ρ A ) n in Eq. (15) can be computed in a very similar manner, since which different from Eq.(17) by removing C (n) A in the trace. We still use the Haar random average Eq. (18), and write tr(ρ A ) n as a sum over all permutations {g p } at all nodes. However because C (n) A is absent, as D 1 the dominant configurations of {g p } satisfy the boundary condition that g p = I at the entire boundary. Thus suppose {g p } has different domains with different g p , domain walls are detached from the boundary, and contain closed curves.
Using the same argument as the above, we can write tr(ρ A ) n as a sum of path integral of Regge action on different M {g p } , similar to Eq. (19). domain walls become the branch curves in M {g p } , which contain closed curves. However, since the intersection of two hyperbolic surfaces cannot give closed branch curves, M {g p } with closed branch curves doesn't admit AdS 3 geometry. Thus the equation of motion doesn't have any solution, except g p = I identically without any domain wall. As a result tr(ρ A ) n is dominant at the configuration that g p = I everywhere where the boundary terms are identical to the ones appearing in Eq. (26). We find that the average Rényi entropy is given by where is the on-shell action of Einstein gravity on 3-manifold M n . The relation Eq. (29) has been an assumption in the existing derivation of RT formula from AdS/CFT [11,12,15]. But it is now derived from the state Eq.(9) using random tensor networks.
Inserting Eq. (25) in Eqs (29) and (30), we obtain the RT formula of Rényi entropy for CFT 2 , which has the desired n dependence.
where L AdS l S corresponds to Ar min the geodesic length in AdS 3 in Eq.(1). The usual RT formula is recovered as n → 1. The above result reproduces the Renyi entropy computed by Hung-Myers-Smolkin-Yale in [14] using the AdS/CFT assumptions. To see it is indeed the right Rényi entropy of the boundary CFT 2 , recall the central charge of CFT relates to L AdS and P by c = 3L AdS 2 P , and l S relates to length l A of boundary interval A by l S 2 ln(l A /δ) in Poincaré patch, where δ is a UV cut-off. It gives which matches precisely the Rényi entropy Eq.(3) of CFT 2 with correct n dependence [51].

IV. CONCLUSION AND OUTLOOK
In this paper we explicitly construct a boundary state |Ψ as a linear combination of random tensor networks to understand the holographic duality at the level of many-body quantum states. |Ψ is shown to reproduce holographically the (1+1)d CFT Rényi entropies S n (A) with the correct n-dependence, which agrees with the known behavior of CFT ground states.
The boundary state |Ψ is given by the holographic mapping applied to the bulk discretized Wheeler-deWitt wave function Φ( a) as a path integral of (2+1)d gravity. The holographic CFT Rényi entropies are derived by the average of random tensors in the tensor networks. The random tensors effectively glue the path integrals and relate boundary Rényi entropies to partition functions of bulk gravity on branched cover spacetimes. The relation between boundary Rényi entropies and bulk gravity partition functions on branched cover spacetimes has been a key assumption in the existing holographic Rényi entropy computations [11,12,15], while our work derives this relation with random tensor networks in the context of AdS 3 /CFT 2 . The derivation also gives the duality between the boundary Rényi entropy and bulk cosmic brane (in AdS 3 /CFT 2 ), which has been proposed as a conjecture in [15].
As another key feature, our derivation doesn't assume the bulk geometry corresponding to the entropy calculation. The geometries on branched cover manifolds are emergent from the boundary CFT state |Ψ via the large bond dimension limit of tensor networks. |Ψ is proposed as a new way to approximate the ground state of large-N CFT.
The results in this paper suggest that |Ψ is a promising realization of the holographic duality in the tensor network, which deserves to be further developed. It is interesting to generalize the derivation to higher dimensional CFTs, which is expected to derive the conjectured duality between the boundary Rényi entropy and bulk cosmic brane in general context. A current research undergoing is to understand the conformal symmetry of |Ψ . |Ψ proposed as the CFT ground state is expected to be (approximately) invariant under boundary conformal transformations. The bulk Wheeler-deWitt wave function and its semiclassical behavior turns out to be a nice tool to understand the symmetry perspectives. It is also important to compare |Ψ with the optimized MERA, given that both approaches captures the behavior of CFT ground state. The comparison may require to extend our analysis to the regime of finite bond dimensions (the finite-N CFT).
The holographic mapping with random tensor networks should be applied to study the excited states in the boundary CFT. The states excited by certain boundary operators acting on |Ψ should be the image of the holographic mapping from a bulk wave function, created by a bulk operator acting on the present Wheeler-deWitt wave function Φ( a). It suggests a holographic relation between boundary and bulk operators and is a promising starting point to understand the bulk reconstruction in the AdS/CFT [61,62]. Since Φ( a) are the coefficients in the expansion of |Ψ with random tensor networks, the relation between boundary and bulk operators should be an analog of the standard duality in quantum mechanics between operators acting on Dirac kets and operators acting on wave functions. This paper for the first time relates Regge calculus, as a popular tool in general relativity, to the tensor network, which is extensively studied in quantum information, condensed matter physics, and the AdS/CFT. The framework developed in this paper is a promising platform for communications between different fields. Regge calculus closely relates to loop quantum gravity, especially its covariant formulation [59,[63][64][65][66][67]. This work is an interesting starting point to explore the relation between loop quantum gravity, AdS/CFT correspondence, and string theory.
Acknowledgments MH acknowledges Ling-Yan Hung and Yidun Wan at Fudan University in Shanghai, Yong-Shi Wu at University of Utah, Bei Zeng at University of Waterloo, and Jie Zhou at Perimeter Institute for Theoretical Physics for the hospitality during his visits, stimulating discussions, and collaborations. MH acknowledges the useful discussions with Xiao-liang Qi and Zhao Yang, and especially thanks Ling-Yan Hung for explaining the result in [14]. MH also acknowledges support from the US National Science Foundation through grant PHY-1602867, and the Start-up Grant at Florida Atlantic University, USA.

Appendix A: Bound on Fluctuation
In this appendix we examine the fluctuation of the Rényi entropy S n (A) from the average value S n (A) in Eq. (31), to qualify how well is the approximation. We show that in the regime P L the fluctuation is generically small. The method used in the following is similar to [25].
We denotes by Z(n) = tr(ρ n A ) and Z(n) ∞ the average value of Z(n) as P L (same as in Eq. (30)). We consider the following fluctuation of Z(n): Z(n) ≥ Z(n) ∞ because in the approximation we made as P L , the neglected terms in the sums are all non-negative. Z(n) 2 is computed in a similar way as the above, using the random average formula Eq.(18), changing n by 2n. It leads to that the dominant contribution of Z(n) 2 is again given by a sum over permutations {g p } at all p, whose boundary condition is g p = I inĀ and g p = (1 · · · n)(n + 1 · · · 2n) in A. Hence Z(n) 2 is also written as a sum of path integrals on different M {g p } . The situation of a single domain wall separating I in RĀ and (1 · · · n)(n + 1 · · · 2n) in R A gives again the dominant contribution. The 3-manifold M {g p } in this case is simply 2 copies of M n . As a result, which implies the following bound Bounding the fluctuation of Z(n) by ε/4 has the following probability by Markov inequality,: Prob Similar conclusion can be drawn for Z(1) n . Bounds on the fluctuations of Z(n), Z(1) n implies the bound on the fluctuation of S n (A). The probability of violating the following bound is where we have used that | ln(1±ε/4)| ≤ ε/2 for small ε. When P /L ε 2 , the above bound of fluctuation is satisfied with a high probability 1 − O P /ε 2 L . This appendix demonstrates a graphic computation of the random average tr(ρ n A ). It is sufficient to demonstrate the case of n = 2, while the generalization to n > 2 is obvious.
FIG .5 lists the basic ingredients of presenting graphically the computation of tr(ρ 2 A ). We have suppress 1 spatial dimension in the presentation where the tensor network | a is illustrated by a chain, and two ends of the chain illustrate regions A andĀ.
FIG.6 is a graphic illustration of tr(ρ n A ) in Eqs. (16) and (17) at n = 2. tr(ρ n A ) contains n copies of Φ * ( a) and n copies of Φ( a).
The random average Eq. (18) is performed for 2n copies of random tensors locally at each node p. Applying Eq.(18) to tr(ρ n A ) gives a sum over choices of permutations {g p } at all p. The computation of tr(ρ 2 A ) is illustrated in FIG.7, which shows a typical term corresponding to two domains with g p = I and g p = C (2) . The second step in FIG.7 uses the relations  in FIG.5 and FIG.8. The red lines illustrate the identification and summation of the boundary data a in path integrals Φ( a), Φ( a) * , which glues the path integrals. The resulting path integral is on the branch cover manifold in FIG.4, where the branch points are in the domain wall S. Each link (p, p ) crossing S associates a normalization factor d[a] −1 , generalized to d[a] 1−n c for an arbitrary pair g p , g p , where n c is the length of the cycle c in g −1 p g p .

Appendix C: Domain Wall Contributions
In this section, we prove that the configuration {g p } with a single domain wall indeed gives the leading contribution to {g p } . Let's consider a more generic case shown in FIG.9(a), where more than one domain-walls are created in the bulk of Σ. We are going to show that this configuration always contribute less than a single domain-wall.
Recall the computation of V(M {g p } ) and the discussion below Eq. (24). When we sum all domain walls and all cycles, in Eq.(C1) where l S is the geodesic length on the hyperbolic plane with unit curvature. Given trivalent intersection of domain-walls, which separate three domains with permutations g 1 , g 2 , g 3 (FIG.10), we have the following "triangle in-  1 − n 2 c n c l S 13 + min l S 12 , l S 23 which is less than a single domain wall contribution. For any intersection with e.g. 4 domain walls, one can always shift the end point of one domain wall away from the intersection and obtain a smaller l S (greater contribution to Eq.(C1)). It reduces the 4-valent intersection back to trivalent situation, which implies the contribution after the above shift is still smaller than the single domain wall configuration. The same argument applies to the intersection with larger number of domain walls.
Therefore we find that the contribution of the multidomain-wall configuration is less or equal to the single domain-wall configuration l S g,g c∈C(g −1 g ) 1 − n 2 c n c ≤ 1 − n 2 n l S .
FIG. 10: A trivalent intersection of domain walls has less contribution than a single domain wall.