{'fraction_non_alphanumeric': 0.036162988115449914, 'fraction_numerical': 0.007541029994340691, 'mean_word_length': 4.349757795943082, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 20, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Holographic spin networks from tensor network states

In the holographic correspondence of quantum gravity, a global onsite symmetry at the boundary generally translates to a local gauge symmetry in the bulk. We describe one way how the global boundary onsite symmetries can be gauged within the formalism of the multi-scale renormalization ansatz (MERA), in light of the ongoing discussion between tensor networks and holography. We describe how to"lift"the MERA representation of the ground state of a generic one dimensional (1D) local Hamiltonian, which has a global onsite symmetry, to a dual quantum state of a 2D"bulk"lattice on which the symmetry appears gauged. The 2D bulk state decomposes in terms of spin network states, which label a basis in the gauge-invariant sector of the bulk lattice. This decomposition is instrumental to obtain expectation values of gauge-invariant observables in the bulk, and also reveals that the bulk state is generally entangled between the gauge and the remaining ("gravitational") bulk degrees of freedom that are not fixed by the symmetry. We present numerical results for ground states of several 1D critical spin chains to illustrate that the bulk entanglement potentially depends on the central charge of the underlying conformal field theory. We also discuss the possibility of emergent topological order in the bulk using a simple example, and also of emergent symmetries in the non-gauge ("gravitational") sector in the bulk. More broadly, our holographic model translates the MERA, a tensor network state, to a superposition of spin network states, as they appear in lattice gauge theories in one higher dimension.

The holographic principle, an anticipated feature of quantum gravity, asserts that at least certain theories of gravity can be described as quantum field theories that live in one less spacetime dimension. For example, in the AdS/CFT correspondence-a concrete realization of the holographic principle-the gravity system lives in a d + 1 dimensional anti-deSitter (AdS) spacetime and is equivalent to a conformal field theory (CFT) that lives on the d dimensional boundary of the spacetime [1,2]. The ex-tra dimension in the bulk spacetime is identified with the length scale of the boundary system, and the renormalization group equations essentially generalize the equations that describe gravity.
Recently, it has been proposed that the multi-scale entanglement renormalization ansatz (MERA) [3]-an efficient representation of ground states of local Hamiltonians on a lattice [4]-realizes at least some features of the AdS/CFT correspondence [5][6][7]. For example, the MERA representation of the ground state of a one dimensional (1D) quantum lattice system is a two dimensional (2D) hyperbolic tensor network, which also describes the RG flow of the ground state. Specifically, the MERA is based on a real space RG transformation, known as entanglement renormalization, that removes local entanglement before coarse-graining the state [8]. In particular, the extra dimension of the tensor network corresponds to length scale of the 1D system.
In Ref. 9 one of us introduced a toy model for holography based on the MERA representation of ground states of 1D local Hamiltonians. The model, dubbed tensor network state correspondence, illustrates a possible way in which the MERA could encode a dual 2D bulk description of a 1D ground state, however, without paying attention to the presence of onsite symmeties in the boundary theory. In this paper, we generalize the model in such a way that an onsite symmetry at the boundary is gauged in the bulk. In particular, this generalization allows us to establish a connection between the MERA and spin networks as they appear in lattice gauge theories in one higher dimension, while also realizing another important feature of the AdS/CFT correspondence using the MERA.

A. Tensor network state correspondence
The basic idea behind tensor network state correspondence is that a tensor network (with open indices) can be viewed as a representation of two different quantum many-body states (belonging to two different Hilbert spaces) depending on how a many-body Hilbert space is associated with the tensor network. We refer the reader for details to Ref. 9. Below, we only briefly summarize the main idea.
The open indices of the tensor network can be associated to sites of a quantum many-body system, specifically, we can use each open index to label an orthonormal basis on a different site of the many-body system. Subsequently, the tensor network defines a quantum manybody state of the system such that the probability amplitude of a given configuration of the sites is obtained by fixing the value of the open indices to correspond to that configuration, and contracting all the tensor networks together by summing over the bond indices. The quantum many-body state obtained following this prescription is referred to as a tensor network state. Examples of tensor network states include the MERA, matrix product states (MPS) [13], and projected entangled pair states [14] (PEPS).
Alternatively, one can associate both the open and bond indices of the tensor network with sites of a larger quantum many-body system, namely, by using each index (open or bond) in the tensor network to label an orthonormal basis on a different site of the system. Subsequently, the tensor network defines a different quantum many-body state whose amplitudes are obtained by fixing the value of all the indices of the tensor network and multiplying together the resulting tensor coefficients, one selected from each tensor. We refer to the many-body state obtained from the tensor network in this way as a tensor network bond state.
Thus, a generic tensor network (with open indices) can be viewed representing either as a tensor network state or as a tensor network bond state. In Ref. 9 we illustrated that the properties of these two states, which are obtained from the same tensor network, are related together in a systematic way. Thus, a tensor network may be viewed as a 'correspondence' between these two quantum many-body states.
Note that a tensor network bond state may be regarded as a regular tensor network state (where degrees of freedom are associated only with open indices and bond indices are summed over) by modifying the tensor network in a particular way. Namely, by inserting a three index copy tensor on each bond of the tensor network, see Ref. 9.
In Ref. 9, by applying this tensor network state correspondence to the MERA we obtained a toy model for holography (without considering symmetries). The tensor network state and the tensor network bond state obtained from a MERA correspond to the boundary and dual bulk state respectively. The bulk states obtained from the MERA in this way exhibit some interesting features. First, the bulk states satisfy an area law entanglement scaling [15]. Second, the bulk entanglement and correlations are organized according to holographic screens. And third, given the MERA representation of a critical boundary state, the boundary correlators of scaling operators (of the underlying CFT) can be obtained from the expectation value of extended bulk operators in certain dual bulk states. Some of these results caricature certain features of the AdS/CFT correspondence as described in Ref. 9.
In this paper, we further develop the toy model. We present a generalized construction of bulk states that retains the three features listed above, but also exhibits new features that result from the presence of a global onsite symmetry in the boundary description.
B. This paper: generalized holographic correspondence in the presence of onsite symmetries In the AdS/CFT correspondence, a global onsite symmetry of the boundary system generally translates to a local gauge symmetry in the dual bulk description [2]. Consequently, the bulk description generally consists, in addition to gravitational and matter degrees of freedom, gauge fields that are described by the boundary global symmetry group. In the quantum gravity regime, the bulk state is expected to be entangled between all these degrees of freedom. In this paper, we describe how these features of the AdS/CFT correspondence can be realized within the framework introduced in Ref. 9.
Symmetries must be properly accounted for in the RG description of a quantum many-body system, in order to reproduce the large-scale properties effectively. For example, consider 1D local, gapped Hamiltonians that have a global onsite Z 2 × Z 2 symmetry corresponding to π rotations about two orthogonal axes. These Hamiltonians can be partitioned into two different equivalence classes or quantum phases, each with distinct large length scale properties: the Haldane phase and the trivial phase [16,17]. More specifically, ground states belonging to the Haldane phase cannot be disentangled to a product state along the RG flow, as long as the RG (entanglement renormalization) transformations protect the symmetry. (Since product states are representative of the trivial phase.) One way to ensure this is to impose that the tensors that implement entanglement renormalization commute with the symmetry. The resulting symmetry-protected entanglement renormalization generates a MERA representation that captures both the expected RG flow of the ground state, and also its global symmetry exactly [18].
In this paper, we consider a local 1D Hamiltonian that has a global onsite symmetry G (which is not broken in the ground state). We represent its ground state by a symmetry-protected MERA and obtain a 'dual' 2D bulk state, by extending the construction of Ref. 9. The bulk state is decribed by a 2D tensor network that is obtained by inserting a 4-index, symmetric copy tensor on every bond of the MERA. Each copy tensor has two open indices, which correspond to bulk degrees of freedom that carry 'left' and 'right' gauge transformations. Thus, our construction, which takes into account the boundary symmetry G, leads to a bulk state in which the symmetry is gauged, thus realizing the holographic translation of a boundary global symmetry to a local gauge symmetry in the bulk. One may, in retrospect, view the (bulk) gauging of boundary symmetries as an underlying motivation for associating the dual bulk degrees of freedom with the bonds of the tensor network, as opposed to e.g. associating them with the tensors as in previous bulk descriptions of the MERA presented in Refs. 10-12. C. Connection to lattice gauge theories and spin networks More broadly, in this paper, we establish a connection between the MERA and lattice gauge theories (on a hyperbolic lattice) in one higher dimension. In a lattice gauge theory, the degrees of freedom are placed on the edges of the lattice, and elementary gauge transformations act on the sites located immediately around a vertex. In our bulk construction, the bulk lattice overlays the MERA tensor network, after it is embedded in a manifold, and the dual bulk degrees of freedom live on the edges of the bulk lattice (that is, the bonds of the tensor network). Elementary gauge transformations act on the bulk sites located immediately around a tensor. In particular, the two open indices of a copy tensor carry the 'left' and 'right' gauge transformations respectively. Mimicking this basic setup of a lattice gauge theory allows us to manifest a bulk gauge symmetry, which is seen to be dual to the global symmetry at the boundary.
We show how the bulk states decompose as a superposition of spin network states, as they appear in a 2D lattice gauge theory with gauge group G, where they span the gauge invariant subspace of the Hilbert space [19]. The spin network decomposition allows us to explore further parallels with holography. One, it reveals entanglement and correlations between the gauge degrees of freedom and the remaining bulk degrees of freedom that are not constrained by the symmetry. And second, by exposing the gauge degrees of freedom in the bulk, the spin network decomposition also allows one to calculate expectation values of gauge-invariant observables in the bulk. We also construct a local, gauge-invariant parent Hamiltonian for the MERA bulk states, see Appendix C.

D. Differences from previous work
A local symmetry also manifests simply in the bulk of a symmetry-protected MERA tensor network representation of a (1D) quantum many-body state with a global onsite symmetry, without reference to a bulk state [18,27,28]. However, we emphasize that in this paper we implement the holographic gauging of a global boundary symmetry more manifestly by means of boundary and bulk quantum states, while Ref. 18 describes the bulk gauging of the boundary symmetry only at the level of the tensor network. Having access to a bulk quantum state, in which the boundary symmetry appears gauged, allows us to probe interesting features in the bulk to explore further connections with holography. For example, we explore the entanglement and correlations between the gauge and the non-gauge degrees of freedom in the bulk, and emergent symmetries in the non-gauge sector of the bulk. On the other hand, no such notions can be defined when the symmetry is gauged only at the level of the tensor network.
In our construction, the local gauge symmetry is hardwired into the bulk tensor network ansatz. Explicit tensor network representations of quantum many-body states with a local gauge symmetry have been presented by other authors [20][21][22][23][24][25][26]. The 2D bulk states that we construct here indeed belong to the gauge-invariant tensor network ansatz e.g. presented in [21]. However, in this paper we focus on the construction of a gauge-invariant bulk state from the MERA representation of a 1D ground state, instead of, say, variationally minimizing the energy of a 2D gauge-invariant bulk Hamiltonian.
Ostensibly, our lifting procedure appears similar to the prescription to gauge quantum states presented in Ref [24]. There the authors describe how to gauge the global symmetry of a tensor network state. Specifically, they consider a 2D quantum many-body state represented by a PEPS tensor network and translate it to another 2D quantum many-body state with a gauged symmetry. In contrast, our construction produces a quantum manybody state in one higher dimension. Moreover, our approach is aimed at building a higher dimensional bulk description of symmetric ground states, whereas Ref [24] is not concerned with applications related to holography.
Some of the previous proposals for drawing a bulk description from the MERA, those presented in Refs. 10-12, associate the bulk degrees of freedom with the tensors of the MERA. In contrast, here we present a bulk description of the MERA by associating bulk degrees of freedom to the bonds of the tensor network, which is closer in spirit to the organization of the degrees of freedom in lattice gauge theories and allows for a more natural introduction of gauge transformations in the bulk and gauging of boundary onsite symmetries, which appears as a general rule of thumb in the AdS/CFT correspondence.

E. Organization of the paper
The paper is organized as follows. In Sec. II we briefly review the symmetry-protected MERA representation of a 1D ground state that has an onsite global symmetry. In Sec. III we describe how to lift the MERA representation to a 2D dual bulk state. In Sec. IV we describe how the bulk states decompose as a superposition of spin network states. In Sec. V, we present numerical results pertaining to the entanglement and correlations in bulk states dual to the ground states of several critical spin chains of interest. For example, we find evidence for a dependence of bulk entanglement on the central charge of the boundary critical system. We conclude with a brief summary and outlook in Sec. VI. The appendices contain some technical discussions and proofs. In Appendix A, we discuss the possibility of emergent topological order in the bulk using a simple example of a system with Z 2 symmetry. In Appendix B, we derive the Schmidt decomposition of the bulk state, which can be used to deduce area law entanglement in the bulk and also to construct a gauge-invariant parent Hamiltonian for the bulk state.
The latter is described in Appendix C.

II. BOUNDARY STATE WITH A GLOBAL SYMMETRY
Consider an infinite 1D lattice L and a compact, completely reducible symmetry group G. Each site of L is described by a Hilbert space V on which the group G acts by means of a unitary representation for all g ∈ G. Also consider a local HamiltonianĤ that acts on the lattice L and has a global symmetry G, namely, ∼ =Vg is a unitary representation of the symmetry group G on site i. We assume that the ground state |Ψ bound ofĤ also has a global symmetry G, namely, The superscript 'bound' appears in anticipation that the ground state plays the role of the boundary state in our holographic correspondence.
In this paper, we represent |Ψ bound by means of an infinite symmetry-protected MERA tensor network. The tensor network is depicted in Fig. 1. An open index o i of the MERA labels an orthonormal basis {|o i } on site i of the lattice L. State |Ψ bound can be formally expanded as where the probability amplitudesΨ o1,o2,... are obtained by contracting the tensor network, which involves summing over all the bond indices-indices that connect the tensors in the network. The MERA representation also describes the RG flow of the ground state. Each layer of tensors of the MERA, separated by dotted lines in Fig. 1, implements a real space RG transformation-known as entanglement renormalization-that maps a lattice L k with L (→ ∞) sites to a coarse-grained lattice L k+1 with L/3 sites. The MERA tensors are chosen so that the renormalization preserves the ground subspace at each step. Subsequent renormalization steps generate a sequence of increasingly coarse-grained lattices: L 0 → L 1 → L 2 · · · , where L 0 ∼ = L is the ultraviolet lattice. Thus, the extra dimension of the tensor network corresponds to length scale, in the sense that the residual tensor network obtained by discarding one or more bottom layers is a representation of the ground state on a coarse-grained lattice. For simplicity, and without loss of generality, in this paper we assume that the ground state |Ψ bound (and the HamiltonianĤ) is translation invariant and scaleinvariant. Specifically, |Ψ bound is a RG fixed point in a gapped or critical phase. Subsequently, a MERA representation of |Ψ bound can be composed from copies of the same two tensors,û andŵ, throughout the tensor network [29], see Fig. 1.
We decorate the indices of the MERA tensors with arrows, as depicted in Fig. 1(a), which indicate how the symmetry acts on the tensors. Tensorsû andŵ are linear transformations from input spaces (incoming indices) to output spaces (outgoing indices) asû : The symmetry acts asV g on an incoming index (input space) and asV † g on an outgoing index (output space). Tensorsû andŵ remain invariant under the action of the symmetry, namely, for all group elements g ∈ G, as depicted in Fig. 1(b). For brevity, we say that tensorsû andŵ are G-symmetric.
The choice of G-symmetric tensors captures the global symmetry, Eq. (2), exactly and also generates a symmetry protected RG flow [18,27]. The tensorsû andŵ are also isometric, namely, they satisfy depicted in Fig. 1(c).

III. DUAL BULK STATE
In this section, we introduce a holographic description of the 1D state |Ψ bound by extending the construction presented in Ref. 9 to the presence of symmetries. We refer the reader to Ref. 9 for a discussion about how the construction is inspired by and implements certain general features of the AdS/CFT correspondence.
Let us embed the MERA in a 2D manifold with a boundary, such that the open indices of the MERA are located at the boundary of the manifold and all the bond indices are located inside the bulk of the manifold. Construct a 2D lattice M on the manifold by locating two sites-each of which is described by the vector space Von every bond of the tensor network. Lattice M is simply a collation of the degrees of freedom that appear in the RG flow of the ground state |Ψ bound , and inherits the hyperbolic geometry of the tensor network.
Next, let us insert the 4-index copy tensor (ĉ) lr l r on each bond of the MERA, as depicted in Fig. 2, such that indices l and r are left open. We will define the components of the copy tensor in the next section (see Eqs. [13][14], but here it suffices to say that, colloquially, tensorĉ copies the basis states on a bond index of the MERA to each of the two open indices l and r. These indices label an orthonormal basis on the two sites located on that bond respectively, and in analogy to lattice gauge theory, we require that these indices carry the left and right gauge transformations on the bulk lattice M respectively (as described in Sec. III A). To this end, we demand that the copy tensorĉ fulfill the following equations that involve the action of symmetry on a single index of the tensor: see Fig. 3.
The new tensor network-the MERA with a copy tensor inserted on every bond-can be viewed as a representation of a quantum state |Ψ bulk of the bulk lattice M, where the probability amplitudes of |Ψ bulk are (formally) obtained by contracting all its bond indices, analogous to how the MERA encodes the state |Ψ bound . Thus, we have 'lifted' the MERA representation of a quantum state |Ψ bound of the 1D lattice L to a quantum state |Ψ bulk of the 2D lattice M. We refer to the bulk tensor network, comprised of copies of the ground state tensorsû,ŵ and the copy tensorĉ, as the lifted MERA.

A. Local gauge symmetry
An immediate consequence of the symmetry conditions Eq. (6) is that the bulk state |Ψ bulk has a local gauge symmetry G. Let us introduce gauge transformations on the bulk lattice M as follows. The symmetry G acts on the two sites located on a bond differently, namely, asV g andV † g respectively. (This choice corresponds to the action of 'left' and 'right' gauge transformations in lattice gauge theory.) Elementary gauge transformations (a) Graphical representation of the copy tensor (ĉ) lr l r . The symmetry acts asVg (blue solid circle) on an incoming index and asV † g (red solid circle) on an outgoing index for all g ∈ G. (b,c) The action of a symmetry operator on index l (index r) transfers to the index l (index r ). The left hand side of each equality depicts the action of the symmetry on the copy tensor according to index arrows, while the right hand side depicts an equivalent action of the symmetry on the tensor (not necessarily according to the arrows). (d) Tensorĉ remains invariant under the action of the symmetry according to the index arrows. Here we illustrate that the lifted MERA, and thus |Ψ bulk , remains invariant under the action of two elementary gauge transformations on the bulk lattice M, corresponding to two different group elements g, g ∈ G respectively. One gauge transformation acts on the 4 bulk sites (counting clockwise starting at the top left in the graphical representation) immediately surrounding tensorŵ (highlighted green) asVg ⊗V † g ⊗V † g ⊗V † g , and the other acts on the 4 bulk sites immediately surrounding tensorû (higlighted yellow) asV g ⊗V g ⊗V † g ⊗V † g . This is shown by means of two equalities: (a)=(b), which results from the symmetry properties of the copy tensor [ Fig. 3 , and (b)=(c), which results from the fact that the tensorŝ u andŵ are G-symmetric [ Fig. 1 act on the 4 bulk sites (counting clockwise starting at the top left in the graphical representation) immediately surrounding tensorû asV g ⊗V g ⊗V † g ⊗V † g , and on the 4 bulk sites immediately surrounding tensorŵ asV g ⊗V † g ⊗ V † g ⊗V † g . General gauge transformations act on larger regions of the lattice M by composing these elementary gauge transformations.
Let us consider the result of applying an elementary gauge transformation on a dual bulk state |Ψ bulk that is represented by a lifted MERA. In the lifted tensor network representation, the action of an elementary gauge transformation corresponds to contracting the symmetry operators on the open indices l, r of the bond tensors located immediately around anû orŵ tensor. The symmetry acts asV g andV † g on the open indices l ('left') and r ('right') respectively. Owing to Eq. (6) [ Fig. 3], oper-atorV g that is applied on an open index of a copy tensor 'slides' through to a bond index of the lifted MERA. Consequently, the action of a gauge transformation on the bulk state translates to contracting the symmetry operators with the tensors around which they are applied, see Fig. 4. However, the tensors are G-symmetric [Eq. (4) and Fig. 1(b)], which eliminates the symmetry operators. Thus, the lifted MERA, and therefore state |Ψ bulk , remains invariant under the action of local gauge transformations.
The gauge invariance of the bulk state is manifest in the same way as appears in lattice gauge theory as originally formulated by Kogut. There one introduces basis states on edges with the local degrees of freedom split into three subspaces as |j, n L , n R where j labels an irrep of the group, and the other two are labelled by matrix elements of these representations (one on the left and one on the right side of the edge). This is a Fourier basis conjugate to the group element labelled basis the two bases being related by the Peter-Weyl theorem. As in lattice gauge theory, the physical states in the bulk are those that are invariant under gauge transformations on a vertex that act with the same group element on all the neighbouring carrier spaces on the edges incident to that site. For Abelian models, all irreps are one dimensional so only one bulk degree of freedom would be needed per edge, but for non-Abelian gauge groups, two labels are needed since the irreps are matrices with components labelled by left and right pairs. We also remark that Elitzurs theorem applies to our bulk state in the sense that the expectation value of non-gauge invariant quantities are trivial by construction.

IV. SPIN NETWORK DECOMPOSITION OF BULK STATES
Let us now introduce a basis in the vector space V, which describes each site of the boundary lattice L and also each site of the bulk lattice M. Under the action of the symmetry, V generally decomposes as where the symmetry acts on space S j by means of the irreducible representation (irrep) of G labelled by quantum number (or charge) j, and D j is the degeneracy space of irrep j. Accordingly, the symmetry operatorsV g decompose asV In particular, note that the symmetry operators act trivially on the degeneracy spaces, namely, as the d j × d j We denote by {|j, m j } an orthonormal basis in the irrep space S j , by {|j, t j } an orthonormal basis in the degeneracy space D j , and by {|j, t j , m j ≡ |j, t j ⊗ |j, m j } the basis on the total space V. For example, if G = SU (2), then the symmetry charge j ∈ {0, 1 2 , 1, 3 2 , · · · } is the total spin, m ∈ {−j, −j +1, · · · , j} is the spin projection along the z-axis. (For simplicity, we assume that G is multiplicity-free.) The description simplifies considerably for an Abelian symmetry, for example G = Z n , U (1), since all the irreps of an Abelian group have dimension 1, that is, dim(S j )=1 for all j.
According to the Wigner-Eckart theorem, the Gsymmetric tensorsû andŵ [Eq. (4)] decompose in terms of the intertwiners of G. If the components of tensorsû andŵ are denoted as (û) pq rs and (ŵ) p qrs respectively, then in the irrep basis where a, b, c, d denote symmetry charges, the tensors decompose asû Here where, for example, e, m e |a, m a ; b, m b ≡ e, m e | · (|a, m a ⊗ |b, m b ) are the Clebsch-Gordan coefficients that describe the change of basis from the tensor product basis |a, m a ⊗ |b, m b to the total charge basis |e, m e . Analogously, we have [see Fig. 5 Finally, (û deg e ) ab cd and (ŵ deg f ) a bcd in Eq. (10) are degeneracy tensors, namely, multi-linear maps between the degeneracy spaces, , and represent the part of the tensor that is not fixed by the symmetry. We refer the reader to Refs [27] for a more detailed exposition on such decompositions of Gsymmetric tensors.
The copy tensorĉ is G-symmetric [ Fig. 3(d)] and therefore also decomposes according to the Wigner-Eckart theorem. The equalities Eq. (6) imply that only the trivial interwiner appears in the decomposition, namely, an intertwiner with trivial intermediate charge. E.g., for G = SU (2) the trivial charge corresponds to the spin 0 irrep. Specifically, tensorĉ decomposes as (see Fig. 5 whereĉ deg a,a is a 4-index degeneracy tensor, a * denotes the conjugate charge of a (namely, charges a and a * fuse to the trivial charge), andÎ 0 is the intertwiner defined according to Eq. (11) for the trivial intermediate charge e = 0. The intertwinerÎ 0 is, in fact, equal to the tensor product of the identityÎ a and the identityÎ a on the irrep spaces S a and S a respectively [Eq. (7)].
Denote the components of a degeneracy tensorĉ deg a,a by (ĉ deg a,a ) The only non-zero components are given by (ĉ deg a,a ) tata tata = 1.
Equations (13)- (14) define the copy tensor. By decomposing tensors according to Eq. (10) and Eq. (13), the entire lifted MERA tensor network decomposes as shown in Fig. 6. Here the sum is over the symmetry charges carried by all the indices of the tensor network, and also the internal intertwining charges that appear in the decomposition of each tensor. The tensor networks appearing on the left in the figure are composed only of the degeneracy tensors, and represent the support of the bulk state |Ψ bulk on the sector of the Hilbert space that is not constrained by the symmetry. On the other hand, the tensor networks appearing on the right in Fig. 6 are composed only of intertwiners of G, and are thus completely fixed by the symmetry. These tensor networks are nothing but spin network states, which here label an orthonormal basis in the support of the bulk state within the gauge-invariant subspace of the bulk lattice M, analogous to their role in lattice gauge theories [19].
In order to make an analogy with the AdS/CFT correspondence, we interpret the degeneracy degrees of freedom as possibly including 'gravitational' degrees of freedom. (Or more generally, 'emergent' gauge degrees of freedom, see Sec. V D.) Thus, in the context of holography, the bulk state |Ψ bulk may be interpreted as an entangled state of gauge fields (described by the spin net-works) living on a 2D quantum geometry (described by the degeneracy tensors). We remark that the bulk construction described here may be readily generalized by also exposing and lifting the internal intertwining charges (that is, the e i 's that appear in Fig. 6 also appear as open indices in the lifted MERA), which allows to incorporate 'gauge matter' in the model. However, we do not pursue this here.
We remark that in a lattice gauge theory, based on a continous gauge group G, one often has to truncate the irreps that appear on the bonds of the spin networks, in order to make calculations tractable. In our bulk construction, the irreps that appear on the bonds of the holographic spin networks are also truncated, since they are carried over from the MERA representation of the ground state. However, the trucation here results from practical considerations in MERA simulations. One systematically assigns only a finite number of irreps on the bonds of the MERA in the variational energy minimization for a given G-symmetric Hamiltonian. Bond irreps are selected with the aim of obtaining the smallest energy possible, within the constraints imposed by the available computational resources.
The spin network decomposition separates the gauge degrees of freedom from the remaining (degeneracy) degrees of freedom in the bulk. This leads to three interesting applications. First, the decomposition allows one to introduce meaningful gauge-invariant observables in the bulk, since it exposes quantum numbers in the bulk (within the gauge-invariant sector). Second, it reveals correlations between the gauge and the remaining degrees of freedom. And third, it allows one to trace out the gauge degrees of freedom and thus probe the nature of the degeneracy degrees of freedom. We make some remarks pertaining to the first application in Sec. IV A below. The second and third applications are explored in Sec. V.

A. Gauge-invariant bulk operators
As mentioned above, the spin network decomposition of the lifted MERA allows one to introduce gaugeinvariant operators in the bulk. Simple examples are operators that act non-trivially on the spin network states and as the identity on the degeneracy degrees of freedom. Figure 7 illustrates a tensor network contraction that equates to the expectation value of such a gaugeinvariant (wilson) loop operator in the bulk. See Appendix D for examples of interesting gauge-invariant loop operators. In the context of holography, it may also be possible to infer some information about the curvature of the ambient space in which the gauge field lives [40]. For a pure gauge theory on a flat space the vacuum state is described as having a flat connection everywhere. However, for curved space, the expectation value will differ in general. Thus, we could expect to infer metric curvature by measuring local holonomies.
Gauge-invariant loop operators may also be used to detect topological order in the bulk. In Appendix A, we explore the topological order of the bulk state for the simple case of Z 2 symmetry. (The discussion readily generalizes to Z n symmetry.)

B. Some properties of the bulk state
The boundary state |Ψ bound is recovered from a bulk state by projecting every pair of sites located on a bond to the state |+ ∈ (V ⊗ V) defined as |+ ≡ a,ta,ma |a, t a , m a ⊗ |a * , t a * , m a * .
State |+ is isomorphic to the identity matrix after the identification |a * , t a * , m a * ↔ a, t a , m a |. Thus, applying the projectorP k ≡ |+ +| on the two bulk sites located on bond k is equivalent to contracting the identity |+ with the copy tensor located on the bond. This contraction results in the identity, as depicted in Fig. 8(a). Thus, the action of the projectorP k eliminates the copy tensor located on bond k of the lifted MERA. By applying the projector on all the bonds of the lifted MERA, all the copy tensors are eliminated and we recover the MERA, and thus the boundary state |Ψ bound .
It is readily checked that the G-symmetric copy tensor is an isometry satisfying the equality depicted in Fig. 8(b). This, along with the fact that the MERA tensorsû andŵ are isometries, ensures that a bulk state is normalized (see Ref. 9, Appendix A), and also exhibits the bulk features of the simpler lifted MERA described in Ref. 9, namely: (i) the presence of holographic screens, (ii) a simple dictionary that translates boundary correlators to expectation values of extended bulk operators, and (iii) a causal cone structure that can be exploited to compute bulk expectation values efficiently. These properties essentially rely on the fact that tensorsû andŵ are isometries.

V. BULK ENTANGLEMENT
Given a subsystem of the bulk lattice M, we define its perimeter and area as the number of sites that are located at the boundary and inside the subsystem respectively. For a generic state belonging to the lattice M, subsystem entanglement entropy is expected to scale as the subsystem's area. In contrast, the subsystem entanglement entropy in a bulk state scales at most as the perimeter of the subsystem, see Appendix B. Such an entanglement scaling is commonly exhibited by ground states of local Hamiltonians in condensed matter physics, where it is often called 'area law entanglement' [15]. In fact, given a lifted MERA, which represents a bulk state |Ψ bulk , one can construct a local, gauge-invariant bulk Hamiltonian whose ground state is |Ψ bulk , as described in Appendix C.
In the remainder of this section we consider bulk states dual to 1D critical ground states, and explore any potential dependence of the bulk entanglement on the central charge of the CFT that describes the critical system in the continuum. We are motivated by the fact that in the AdS/CFT correspondence, the leading order of quantum fluctuations in the bulk is O(1/c) where c 1 is the central charge of the CFT [34].
However, in order to compare bulk properties corresponding to different critical boundary states one has to address an ambiguity that arises from the fact that the MERA representation of a 1D ground state is not unique as discuss in the next section. Eq. (14), selects out a particular MERA representation of the ground state-the one expressed in a bond basis in which the degeneracy tensors {ĉ deg a } have these components. On the other hand, the degeneracy copy tensors {ĉ deg a } 'commute' only with diagonal matrices. Namely, a contraction ofĉ deg a with a diagonal matrix on any index is equal to a contraction of the tensor with the same diagonal matrix on any other index. This implies that the bulk states obtained by lifting different MERA representations of the same ground state are not generally related to each other by one-site unitary transformations on the bulk lattice, and therefore they have different entanglement. Thus, our bulk construction generally relates a given ground state to a set of bulk states with different entanglement.
However, in this paper, we restrict attention to MERA representations that are made of G-symmetric and isometric tensors. While G-symmetric tensors ensure that the bulk state-obtained by lifting the MERA by inserting copies of the G-symmetric copy tensor-has a local gauge symmetry G (as described in Sec. III A), the choice of isometric tensors leads to the desirable bulk features listed (i)-(iii) in Sec. IV B.
To this end, we restrict {M k : V → V} k to unitary matrices that commute with the symmetry, namely, [M k ,V g ] = 0 for all g ∈ G. Since V decomposes as Eq. (8), Schur's lemma (a special case of the Wigner-Eckart decomposition) implies that matrixM k decomposes asM k = a (M k,a ⊗Î ηa ). Thus, the bond transformations are restricted to act as the identityÎ ηa on the bonds of the spin networks, which also restricts the set of the dual bulk states. In particular, one can exploit this restriction on the bond transformations to partially fix a basis on the total bond space V, in the different MERA representations of |Ψ bound . Specifically, we fix the irrep basis {|a, m a } on the bonds of the spin networks, while a basis on the bonds of the degeneracy tensor networks corresponds to a choice of the bond transformationsM k,a (with respect to a given MERA representation).
Therefore, here we probe for any statistical dependence of the bulk entanglement on the boundary central charge, by randomly sampling from the set of all allowed dual bulk states. Recall that we only consider bulk states that are obtained by lifting MERA tensor networks composed of G-symmetric and isometric tensors. (This corresponds to restricting the intrinsic bond transformations {M k } to unitary matrices that commute with the symmetry G.)

B. Critical spin chains
To this end, we considered the ground states of the following 1D critical spin models: where i labels sites of a 1D infinite lattice on which the Hamiltonian acts,σ x ,σ y ,σ z are Pauli matrices, the oper-atorŜ α is the α component of the spin−1 representation of su (2) Table I. We determined a symmetry-protected MERA representation of the ground state of each of these models using the variational energy minimization algorithm for the scale-invariant MERA de-scribed in Ref. 29, adapted to the presence of symmetries [27,35]. We considered only Abelian symmetries here, which appear either as the total symmetry or as subgroup symmetry. Specifically, we obtained a Z 2 -symmetric MERA representation for the ground state of the Ising model and the Blume-Capel model, and a Z 3 -symmetric MERA representation of the ground state of the Potts model.   (15). We use a compact notation a(da) to denote an irrep a and its degeneracy da that appears in the irrep decomposition, Eq. (7), of the Hilbert space of one site of the lattice. The two irreps of Z2 are labelled by 0 and 1 respectively. The three irreps of Z3 are labelled by 0, 1 and 2 respectively. We label the two irreps of U (1) that appear on each site of the XXZ model by −1 and 1. For example, for the Blume-Capel model 0(2) ⊕ 1(1) denotes that each site of the lattice decomposes as the direct sum of two copies of Z2 irrep 0 and one copy of Z2 irrep 1.  Table II.) The total bond dimension (namely, the dimension of the total bond representation) is equal to 12 for all the simulations. For the XXZ model, we obtained both a Z 2 -symmetric MERA representation of the ground states for ∆ ∈ {0, 0.71, 0.81, 0.87, 1} (corresponding to a Z 2 subgroup symmetry), and also an U (1)-symmetric MERA representation for ∆ ∈ {0, 1}. The U (1) symmetry of the XXZ models with ∆ ∈ {0, 1} corresponds to the total symmetry for ∆ = 0 and a subgroup symmetry for ∆ = 1.

model symmetry site representation
The representation of the symmetry on each site of the lattice for the various models is listed in Table II. Table  III lists the symmetry representation that we assigned to the MERA bonds for the ground state simulations. We tried a few different charge and degeneracy combinations and the choices listed in the table III resulted in the smallest error in the ground state energy density as determined from the resulting MERA.
The error in the estimated ground state energy density for the Ising model was O(10 −8 ) and the relative error in the estimated central charge was 0.6%. For the remaining models, the error in the estimated ground state energy density was at most O(10 −4 ), and the relative error in the estimated central charge was at most 3%. For all the models, the relative error in the estimated smallest six D 'tower' of bulk sites FIG. 9. The bulk sites (red squares) located in the region highlighted yellow were considered to obtain the plots shown in Fig. 10 and Fig. 11. These sites are located along an infinitely long tower of theŵ tensors. The bulk sites located in the region highlighted blue were considered to obtain the entanglement negativities listed in Table IV. These consist of sites located around a loop of tensors and two sites located at the bottom of the loop.
The models listed in Eq. (15) are not scale-invariant, but flow to a scale-invariant fixed point after possibly several RG (entanglement renormalization) steps. We discarded the non-scale invariant part of the MERA before lifting it to a bulk state. (That is, we considered the renormalized scale-invariant ground state of each model.)

C. Bulk entanglement vs boundary central charge
Before proceeding to our results, we remark that defining entanglement entropy in gauge-invariant states is subtle since a gauge-invariant Hilbert space does not usually have a tensor product structure. A possible approach, one that we have followed here, is to embed the Hilbert space into a larger tensor product space-the tensor product of the Hilbert spaces on each of the links of a lattice gauge theory. See, for example, a recent work  For the ground state of each of the models listed in Eq. (15), we randomly selected 10 5 dual bulk states (restricting the corresponding bond transformations {M k } to unitary matrices that commute with the respective symmetry), and computed the second Renyi entanglement entropy R tower per site, R tower ≡ −log 2 (Tr(ρ tower ) 2 ). (16) Hereρ tower is the reduced density matrix of all the bulk sites located along the infinitely long tower ofŵ tensors (highlighted yellow in Fig. 9). We partitioned the Renyi entanglement entropy density values in to 100 equally spaced bins. Figure 10 shows the probability distribution of the Renyi entanglement entropy R tower per site for the critical Ising model, which illustrates that the different bulk states indeed have different entanglement.
In Fig. 11 we plot the maximum, minimum and the mode of the probability distributions of R tower per site for all the critical models. The plot indicates a statistical trend that the Renyi entropy density generally increases with increase in the boundary central charge. Note also the clustering of data for different XXZ models, which have the same central charge. From these results it appears that the bulk entanglement entropy R tower depends predominantly on the central charge, as compared to other microscopic details of these models.
The plot in Fig. 11 also suggests that a MERA representation based on a larger subgroup symmetry may correlate with a decrease in the Renyi entropy R tower . Specifically, for ground states of the two XXZ models with ∆ = 0, 1 the maximum, mode and minimum Renyi entropies obtained from the U (1)-symmetric MERA representation were found to be smaller than those obtained from the Z 2 -symmetric MERA representation. On the other hand, the two representations gave approximately equal estimates for the ground state energies, central charges, and few lowest scaling dimensions.
Note that a Z 2 -symmetric and a U (1)-symmetric MERA representation of a given ground state are expected to correspond to two bulk states with different entanglement respectively. This is because a U (1)symmetric MERA representation can be converted to a Z 2 -symmetric MERA representation by applying bond transformations to change the bond basis from a U (1) irrep basis to the subgroup Z 2 basis listed in Table III, which likely alter the bulk entanglement. However, we do not know how to account for the decrease in entropy when the larger symmetry was considered here, and whether this behaviour is more general than illustrated by these results.

D. Entanglement between gauge and degeneracy degrees of freedom
Finally, we probed the bulk entanglement between the gauge and degeneracy degrees of freedom for the case of Z 2 and Z 3 symmetry. We considered a small region of the bulk lattice and obtained a reduced density matrix by tracing out all degrees of freedom outside the region, and also the gauge degrees of freedom inside the region. This was achieved by using the spin network decomposition of the bulk state, which exposes separate open indices in the lifted MERA corresponding to the gauge and non-gauge degrees of freedom respectively. In order to trace out the gauge degrees of freedom in a region, one also contracts the open indices of the spin networks that are located with the region but not the corresponding degeneracy indices.
The smallest region for which we found non-zero entanglement negativity, a measure of quantum entanglement, is depicted as region D in Fig. 9. Letρ [D] denote the reduced density matrix of region D by tracing out all bulk sites outside D, and also the gauge degrees of freedom inside D. We computed the entanglement negativity where λ i are the eigenvalues of the matrix obtained by taking the partial transpose ofρ [D] with respect to some of the sites in D. (A non-zero value of the negativity indicates that the state has quantum entanglement.) We selected two different bulk states dual to the ground state of each critical model listed in Eq. (15), and computed the value of n(ρ [D] ) for both these bulk states. These values are listed in Table IV. The fact that this entanglement negativity is positive for these models indicate that (at least) these dual bulk states have quantum entanglement between the gauge and degeneracy degrees of freedom.
The plot in Fig. 12 shows the spectrum of the reduced density matrixρ deg k , obtained by tracing out all bulk sites except the bulk site k and also tracing out the gauge degrees of freedom on the site k. (That is,ρ deg k has support only on the non-gauge sector of site k.) Notice the appearance of approximate degeneracies in the spectrum. (We sampled a few different bulk states, this plot is illustrative of the typical degeneracies that we observed.) One possible way to account for these degeneracies is the emergence of a non-Abelian symmetry in the nongauge ('gravitational') sector of the bulk. For instance, if the bulk state has an (emergent) non-Abelian on-site symmetry, say which acts only on the non-gauge sector of the Hilbert space, then the reduced density matrixρ deg k must commute with this symmetry. Consequently, by applying Schur's lemma,ρ deg k must decompose aŝ Here j is an irrep of the emergent symmetry,ρ deg k,j is a density matrix that acts on the degeneracy space of charge j, andÎ ηj is the η j × η j identity matrix that acts on the irrep j. The spectrum ofρ deg k is clearly degenerate, in accordance with this decomposition, specifically, the degeneracy of an eigenvalue ofρ deg k,j is at least η j . In the scenario of an emergent symmetry, the degeneracies in the spectrum ofρ deg k can be used to infer a possible set of emergent symmetry charges, which can be used to decorate the bonds of the degeneracy tensor networks that appear in Fig. 6 (analogous to how the spin networks are decorated with the symmetry charges). Broadly speaking, in this case, it may be possible to further decompose the degeneracy tensor networks in terms of spin networks composed of intertwiners of the emergent symmetry, thus refining the bulk construction presented in this paper. We leave further exploration of any emergent bulk symmetries for future work.

VI. SUMMARY AND OUTLOOK
In this paper, we described a toy model for constructing a holographic description of a 1D quantum lattice system, equipped with the action of a local Hamiltonian that has a global onsite symmetry G. Specifically, we lifted a MERA representation of the ground state, which also has the global symmetry, to a tensor network representation of a quantum state of a 2D lattice on which the symmetry appears gauged. This was achieved by embedding the MERA in a 2D manifold, and inserting 4-index tensors on the bonds of the tensor network. The 1D ground state and the dual 2D quantum state are seen to live on the boundary and in the bulk of the manifold respectively. In order to manifest a gauge symmetry in the bulk, it was essential to use G-symmetric tensors, which compose the MERA representation and generate a symmetry protected RG flow, and require that the copy tensors, which were used to lift the MERA, fulfill particular symmetry properties, those depicted in Fig. 3.
In this way, our toy model translates a 1D boundary state with a global onsite symmetry to a 2D bulk state in which the symmetry appears gauged. In the AdS/CFT correspondence, a global onsite symmetry at the boundary is also gauged in the bulk as a general rule of thumb. In light of the ongoing discussion between the MERA and holography, we take the view that any legitimate bulk de-scription of the MERA must implement the holographic gauging of global boundary symmetries. In particular, making this demand may narrow the choices for the possible bulk degrees of freedom. As we have shown in this paper, the holographic gauging of boundary symmetries is very conveniently realized by associating bulk degrees of freedom with the bonds of the MERA, as opposed to its tensors as has been considered in some of the previous works [10][11][12], since it allows us to introduce gauge transformations as in lattice gauge theory.
We further showed how the bulk states decompose as a superposition of spin network states, which label a basis in the gauge-invariant sector of the bulk Hilbert space. Thus, our bulk construction brings together tensor network states and spin network states, as they appear in lattice gauge theories with gauge group G. The spin network decomposition of the bulk state allows one to introduce meaningful gauge-invariant observables in the bulk, since it exposes quantum numbers in the bulk (within the gauge-invariant sector). It also allows us to explore further parallels with holography. For example, the decomposition reveals correlations between the gauge and the remaining ('gravitational') degrees of freedom, and is also instrumental to probe any emergent symmetries in the non-gauge ('gravitational') degrees of freedom (since the spin network decomposition allows one to trace out the gauge degrees of freedom in the bulk).
Spin networks also appear in various quantum gravity models where they label a gauge-invariant basis in the kinematic Hilbert space of the theory, for example, in loop quantum gravity [37]. Towards the completion of this work, we found recent papers which also explore connections between tensor network states and spin network states, specifically as they appear in loop quantum gravity [38], and in the context of group field theory [39].
This work demonstrates a useful toy model for exploring basic features of holography using tensor networks. Beyond holography, our formalism may be viewed as a general correspondence between a 1D ground state with a global symmetry G and a 2D many-body state with a local symmetry G, which may also be useful in characterizing and relating together different types of quantum phases of matter as illustrated in Appendix A.

Appendix A: Bulk topological order and isometric tensors
In this appendix, we discuss the topological order of the bulk states obtained from the MERA as described in this paper, applied to the case of Z 2 symmetry. (The discussion can be readily generalized to Z n symmetry.) In particular, we illustrate an interesting interplay between bulk (Abelian) topological order and the choice of isometric tensors. More specifically, the bulk state obtained by lifting a MERA made of Z n -symmetric and isometric tensors does not have a non-trivial Z n topological order.

Example using Z2 symmetry
For the purpose of this section, we specialize the notation introduced in Sec. II to the case of a Z 2 symmetry. Let L here denote an infinite 1D lattice, each site of which is described by vector space V ∼ = C 2 and is equipped with the action of the group Z 2 = {Î,Ẑ}. The group acts on the space V by means of the unitary rep-resentationÎ = 1 0 0 1 ,Ẑ = 1 0 0 −1 . Under the action of the symmetry, the space V decomposes as where V e and V o are the two irreps of Z 2 . We denote by |e ≡ 0 and |o ≡ 1 a basis in the one dimensional vector spaces V e and V o respectively.
Let |Φ bound denote the (unnormalized) GHZ state belonging to the lattice L, where |± = (|e ± |o ) and e.g. | + + · · · ≡ (|+ ⊗ |+ ⊗ · · · ). State |Φ bound has a global Z 2 symmetry since |Φ bound = (Ẑ ⊗Ẑ ⊗ . . .)|Φ bound . Consider a Z 2 -symmetric MERA representation T of |Φ bound comprised of copies of two simple tensorŝ which replace copies of the tensorsû andŵ in Fig. 1 respectively. Tensorû GHZ is simply the identity, and the components ofŵ GHZ are: Note that tensorŵ GHZ is an isometry satisfying It is readily checked that tensorŵ GHZ is also In order to verify that the MERA tensor network T indeed represents the state |Φ bound , Eq. (A1), we need to to contract all the tensors of T to obtain the probability amplitudesΦ i1i2··· , where i 1 i 2 · · · ∈ {e, o} denote the open indices of T . The simple tensor network T can be contracted algebraically to obtain That is, |Φ bound is an (unnormalized) equal superposition of kets |o 1 o 2 · · · ≡ |o 1 ⊗ |o 2 ⊗ · · · labelled by bit strings o 1 o 2 · · · with even number of 1's. In the basis |± on each site, this state is simply the GHZ state |Φ bound . Let us lift the tensor network T by inserting the Z 2symmetric copy tensorĉ [Z 2 ] on each bond of the tensor network. The copy tensorĉ [Z 2 ] is defined by specializing Eq. (13) and Eq. (14) to Z 2 , namely, where i, j, k, l = {e, o}. By construction, the bulk state |Φ bulk represented by this lifted tensor network has a local Z 2 gauge symmetry, as described in Sec. III. We now ask whether the bulk state |Φ bulk has Z 2 topological order? In order for the state |Φ bulk to have a Z 2 topological order it must be invariant under the action of two non-commuting, deformable loop operators [30,31]. One considers two types of loops, namely, (A) paths comprised of a closed sequence of the tensor network bonds, and (B) closed paths in the ambient manifold (outside the tensor network) that intersects only the bonds of the tensor network, such that the two bulk sites associated with the intersected bonds are located inside and outside of the loop respectively.
State |Φ bulk is invariant under any type B loop ofẐ's, which follows simply from the fact that |Φ bulk has a local Z 2 gauge symmetry. See Lemma 1 in Appendix A 3. In addition, |Φ bulk must be invariant under the action of type A loops ofX operatorŝ X ≡ V → V,X ≡ |e o| + |o e|.
(A8) However, it can be shown that the bulk expectation value of any type A loop ofX's is identically zero, see Lemma 2 in Appendix A 3. This implies that the state |Φ bulk does not have Z 2 topological order. However, since the MERA representation of a quantum many-body state is generally not unique (see, for example, the discussion in Sec. V A). Is there then a different MERA representation, composed of isometric and Z 2 -symmetric tensors, of the GHZ state that perhaps has Z 2 topological order? The answer is still no, since the proof presented in Appendix A 3 applies to any lifted MERA that is obtained by replacingû GHZ andŵ GHZ with arbitrary isometric and Z 2 -symmetric tensors.
Next, consider a different MERA tensor network T topo comprised of copies of tensorsŵ topo ≡ŵ GHZ , andû topo defined as: Analogous to T , the tensor network T topo also represents the GHZ state |Φ bound . Namely, by contracting all the tensors of T topo and changing the site basis to |± one obtains the probability amplitudes in Eq. (A1). On the other hand, tensor network T topo cannot be obtained from T by applying bond transformations. The bulk state |Φ bulk topo remains invariant under the action ofX operators applied on bulk sites located around a plaquette on the bulk lattice, here illustrated for the action ofX's on the highlighted plaquette (red). This is shown by means of two equalities. The first equality results from applying the equality depicted in Fig. 14(a) to all the copy tensors located around the plaquette. The second equality results from using the equalities depicted in Fig. 14(c)-(d), which eliminates theX operators.
Tensorû topo is Z 2 -symmetric sinceû topo = (Ẑ ⊗ Z)û topo (Ẑ † ⊗Ẑ † ). However,û topo is neither an isometry nor a unitary tensor. Instead,û topo is a projector, [That is, (û topo ) 2 = (û topo ).] Let |Φ bulk topo denote the bulk state obtained by lifting T topo , by inserting copies of the Z 2 -symmetric copy tensor c [Z 2 ] on the bonds of T topo . We show that |Φ bulk topo is the ground state of the Z 2 surface code Hamiltonian, here defined on a hyperbolic lattice [30]. The ground state of the Z 2 surface code is known to have Z 2 topological order. Let v and p denote the vertices and plaquettes of the bulk hyperbolic lattice, and s(v) and s(p) denote the set of bulk sites immediately surrounding vertex v and those located around plaquette p respectively. The Z 2 surface code Hamiltonian is defined here aŝ The bulk state |Φ bulk topo is the ground state ofĤ topo because it remains invariant under the action of the vertex terms in Eq. (A11) (which are simply Z 2 gauge transformations)-by virtue of our bulk constructionand also under the action of the plaquette terms in Eq. (A11), as illustrated in Fig. 15.
Since each pair of bulk sites associated with a bond are effectively supported only on a qubit subspace, by virtue of the gauge symmetry, our bulk state is indeed stabilised by a set of vertex and plaquette operators equivalent to the surface code state. As the state |Φ bulk has no further symmetries or additional degrees of freedom then this is the same phase as the surface code: a phase with Z 2 topological order. We also refer the reader to Ref. 32 where the topological properties of the lifted tensor network, which represents state |Φ bulk topo , are explicitly demonstrated. (In Ref. 32 the authors consider a 3-index copy tensor, which is isomorphic to the Z 2 -symmetric copy tensorĉ [Z 2 ] after each pair of bond sites is projected to an effective qubit space.) The bulk state represented by the lifted T topo tensor network turned out to have topological order because the tensor network contains non-isometric tensors, and thus avoids the argument for the absence of bulk topological order presented in Appendix A 3. Thus, this example reveals an interesting interplay between the presence of an Abelian topological order in the bulk and the isometric constraints that are usually imposed on the MERA tensors. On the other hand, for non-Abelian symmetries, isometric tensors may be compatible with the presence of a bulk topological order. Viewed as ground states of local Hamiltonians, the quantum many-body states |Φ bound and |Φ bulk topo , described above, belong to two different quantum phases of matter: the boundary state |Φ bound (a GHZ state) belongs to a quantum phase with spontaneously broken Z 2 symmetry, while the bulk state |Φ bulk topo belongs to a quantum phase with Z 2 topological order. (In a Z 2 symmetry broken phase, the ground state is 2-fold degenerate. The only Z 2 -symmetric ground states are GHZ type states dressed with local entanglement, see e.g. Appendix C in Ref. 33.) The Z 2 -symmetric ground states at the RG fixed point in the phase are exactly GHZ states. Recall that the MERA representation of a ground state describes the RG flow of the ground state to a fixed point wavefunction, which is characteristic of the quantum phase to which the ground state belongs [18]. If the ground state belongs to a Z 2 symmetry broken phase, and if we target the Z 2 -symmetric ground subspace by protecting the symmetry along the RG flow (employing Z 2 -symmetric tensors) then the ground state flows to a GHZ state.
Subsequently, if the fixed point tensors in the MERA representation of a Z 2 -symmetry broken ground state are, in fact,û topo andŵ topo (which represent a GHZ state) it is tempting to conclude that our bulk construction leads to an emergent Z 2 topological order at the RG fixed point in a 1D Z 2 symmetry broken phase. Here by emergent we mean that we can systematically obtain a quantum state with Z 2 topological order from the MERA representationof a Z 2 symmetry broken ground state (by lifting the MERA). However, this argument does not quite work since the MERA representation of the GHZ state is not unique. For example, as mentioned previously, the GHZ state may also be represented by the MERA T , which does not lift to a topologically ordered bulk state.
On the other hand, we have isolated a condition under which our holographic correspondence could lead to an emergent Abelian topological order in the bulk. Namely, if non-isometric tensors are permitted, and in fact preferred, in the MERA representation of a ground state, while protecting the symmetry along the RG flow. In this case, the fixed point tensors in a Z 2 -broken phase are indeed given by Eq. (A4), and therefore one can argue for an emergent topological order as described above.

Proofs
In this section, we prove two lemmas that were used in the discussion presented in this Appendix.
Consider a vector space V ∼ = C 2 that is equipped with the action of the group Z 2 = {Î,Ẑ}. The group acts on the space V by means of the unitary representation The tensor network contraction resulting from simplifying the tensors located below the loop. First, the tensors located at the bottom are contracted with their adjoints and thus cancel out. Consequently, pairs of copy tensors are seen to be contracted together and also cancel out, thanks to the equality depicted Fig. 8(b). (c) Zoom in to the contraction around tensor v4. Tensorû cancels with its adjoint, resulting in two separate contractions. Each of these is identically zero since the trace ofX is 0.
Under the action of the symmetry, the space V decomposes as V ∼ = V e ⊕V o where V e and V o are the two irreps of Z 2 . Denote by |e and |o a basis in the one dimensional vector spaces V e and V o . Also, define the Z 2 irrep flip operatorX ≡ |e o| + |o e|. Consider a MERA tensor network T composed from arbitrary isometric and Z 2 -symmetric tensorsû : Let T denote the lifted MERA obtained by inserting the Z 2 -symmetric copy tensor, Eq. (14), on the bonds of T .
Lemma 1. Consider a loop C in the ambient manifold, in which the (lifted) MERA is embedded, (i) that intersects only copy tensors, and (ii) the two bulk sites associated with the open indices of an intersected copy tensor are located inside and outside of the loop respectively. Also, consider the loop operatorẐ C ≡ iẐ i that acts on all bulk sites i located immediately inside loop C. The lifted MERA T , and thus the bulk state it represents, is invariant under the action ofẐ C .
Proof. As an illustration of the general proof, consider the specific loop operator applied along the loop enclos- where µ α > 0 are the Schmidt coefficients, and {|Ω [A] α } and {|Ω [B] α } is an orthonormal basis in spaces V (A) and V (B) respectively. The decomposition B1 is useful since the reduced density matrix of the parts A and B is diagonal in the Schmidt basis, namely, In particular, the rank of the reduced density matriceŝ ρ [A] andρ [B] is equal to n.
In this appendix, we derive a Schmidt decomposition for a bulk state |Ψ bulk . We will use this Schmidt decomposition to derive two results: (i) the bulk state exhibits an area law scaling of entanglement, and (ii) |Ψ bulk can be viewed as the ground state of a local, gauge-invariant Hamiltonian.
For our purposes, we decompose the copy tensorĉ, Eq. (14), in terms of a 3-index, G-symmetric tensorx depicted in Fig. 17. In the irrep basis, |e ≡ |a, t a , m a , |f ≡ |b, t b , m b , |g ≡ |c, t c , m c , |h ≡ |d, t d , m d , the index values that correpond to non-zero components ofĉ satisfy The intermediate index γ ≡ (0, t γ ) carries only the trivial charge and takes a d a number of values. We establish a one-to-one correspondence between the index γ and the degeneracy index (a, t a ) and denote it as γ ↔ γ(a, t a ). In the irrep basis, tensorx decomposes as (Wigner-Eckart theorem)x where the only non-zero components of the degeneracy tensorsx deg a are (x deg a ) ta γ(a,ta),ta = 1 4 √ η a , for all t a ∈ {1, 2, . . . , d a }. (B7) Next, consider a closed path on the ambient manifold-into which the (lifted) MERA is embeddedthat intersects only the new bonds resulting from decomposing the copy tensors, as illustrated by the dashed contour depicted in Fig. 17. Such a path bipartitions the lifted MERA, and the bulk lattice M, into parts P and R. The two bulk sites that are associated with each intersected copy tensor are split amongst parts P and R respectively.
LetP andR denote the tensors obtained by contracting all the tensors located in parts P and R respectively. Also, let p ≡ (p 1 , p 2 , . . . , p N , p N +1 , . . .) and r ≡ (r 1 , r 2 , . . . , r N , r N +1 , . . .) denote the tuple of all open indices located in P and R respectively.
The bulk state |Ψ bulk can be expressed as This can be understood as follows. The one-to-one map γ(a, t a ) also corresponds to a one-to-one map between the kets |γ ↔ |a, t a . Two different tuples α and α differ in some entry γ i , and therefore also correspond to different elements of the tensor product basis |p and |p (and also |r and |r ). However, p|p = δ p,p , thus leading to Eq. (B10). labelled by the conjugate irrep j * . In the spin network basis the holonomy operator is given bŷ where D (j) (h) is the irrep j of the group element h ∈ G, and the sum is taken over all irrep labels. The operator B j p inserts a flux of type j into the plaquette and its matrix elements in the spin network basis can be determined using recoupling formulae, see e.g. Ref. 41.
Generically, consider a plaquette which is an n-gon with boundary edges carrying irreps j 1 , j 2 , . . . j n , and incident edges to the vertices labelled by k 1 , k 2 , . . . k n . The spin network for this state can be labelled |Φ; k 1 , . . . k n ; j 1 , . . . j n where Φ indicates the configurations of all the other edges not touching the plaquette. Assuming for simplicity that all boundary edges have the same orientation as p and all incident edges are directed toward the vertices of p, the matrix elements for the plaquette flux operator are: Φ; k 1 , . . . k n ; j 1 , . . . j n |B s p |Φ; k 1 , . . . k n ; j 1 , . . . j n =   n l=1,r=1 where the F 's are recoupling coefficients (6-j symbols) that describe coupling the three irreps (labelled by superscripts) to a fourth total irrep (labelled by subscript). For other orientations of the edges simply replace the labels of the reverse oriented edges in the expressions above by their conjugates. Let us consider specific examples. where o(e, p) = 1(−1) if the edge e is oriented the same direction (opposite direction) to the orientation of the plaquette boundary ∂p. Plaquettes can be give a uniform orientation which we choose to be counterclockwise. Herê X j = d−1 k=0 |k ⊕ d j k|, where ⊕ d is addition modulo d. The group Z 2 is particularly simple having self conjugate irreps meaning the only allowed gauge field configurations are non intersecting loops yieldingB 0 p = 1 and B 0 p = e∈∂pX . The group U (1).-For the gauge group G = U (1), the irreps are labelled by integers, j ∈ Z, and the conjugate irrep satisfies j * = −j. Consider a vertex of the spin network with two incoming edges labeled by irreps j 1 and j 2 and an outgoing edge labelled j 3 . These must satisfy the branching rule j 1 + j 2 − j 3 = 0 implying only closed strings of irreps or sums thereof appear on the network. The plaquette operator iŝ whereL j = k∈Z |k + j k|. These operators are infinite dimensional but when the irreps appearing in the spin network are truncated so that j min ≤ j ≤ j max , then one can use finite dimensional truncations ofL j .

Non-Abelian gauge groups
The group SU (2).-The irreps j are labelled in the set j ∈ {0, 1 2 , 1, 3 2 , 2, . . . .} and the irreps are self dual, j * = j. The branching rule at vertex with two incoming edges of irreps j 1 and j 2 and one outgoing edge of irrep j 3 satisfies |j 1 − j 2 | ≤ j 3 ≤ j 1 + j 2 , implying branching strings nets can occur on the network. One must use the full expression for the matrix elements ofB j p where the F matrices are proportional to the Wigner 6 − j symbols: