A Bulk topological order and isometric tensors

# Holographic spin networks from tensor network states

## Abstract

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.

## I Introduction

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 dimensional anti-deSitter (AdS) spacetime and is equivalent to a conformal field theory (CFT) that lives on the dimensional boundary of the spacetime (1); (2). The extra 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.

### i.1 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 many-body 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.

### i.2 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 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 (17); (16). 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 (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 , 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); (11); (12).

### i.3 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 , 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.

### i.4 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 many-body 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); (11); (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.

### i.5 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 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 and a compact, completely reducible symmetry group . Each site of is described by a Hilbert space on which the group acts by means of a unitary representation

 ^Vg:V→V,   ^Vg^V†g=^V†g^Vg=^I,

for all . Also consider a local Hamiltonian that acts on the lattice and has a global symmetry , namely,

 [^H, ⨂i^V(i)g]=0,   for all g∈G, (1)

where is a unitary representation of the symmetry group on site . We assume that the ground state of also has a global symmetry , namely,

 |Ψbound⟩=(⨂i^V(i)g)|Ψbound⟩. (2)

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 by means of an infinite symmetry-protected MERA tensor network. The tensor network is depicted in Fig. 1. An open index of the MERA labels an orthonormal basis on site of the lattice . State can be formally expanded as

 |Ψbound⟩=∑o1,o2,…^Ψo1,o2,…|o1⟩⊗|o2⟩⊗⋯ (3)

where the probability amplitudes 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 with sites to a coarse-grained lattice with 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: , where 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 (and the Hamiltonian ) is translation invariant and scale-invariant. Specifically, is a RG fixed point in a gapped or critical phase. Subsequently, a MERA representation of 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 and . The symmetry acts as on an incoming index (input space) and as on an outgoing index (output space). Tensors and remain invariant under the action of the symmetry, namely,

 ^u=(^Vg⊗^Vg) ^u (^V†g⊗^V†g),^w=(^Vg) ^w (^V†g⊗^V†g⊗^V†g), (4)

for all group elements , as depicted in Fig. 1(b). For brevity, we say that tensors and are -symmetric. The choice of -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

 ∑rs(^u)pqrs(^u†)rsp′q′=δpp′δqq′, ∑qrs(^w)pqrs(^w†)qrsp′=δpp′, (5)

depicted in Fig. 1(c).

## Iii Dual bulk state

In this section, we introduce a holographic description of the 1D state 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 on the manifold by locating two sites—each of which is described by the vector space —on every bond of the tensor network. Lattice is simply a collation of the degrees of freedom that appear in the RG flow of the ground state , and inherits the hyperbolic geometry of the tensor network.

Next, let us insert the 4-index copy tensor on each bond of the MERA, as depicted in Fig. 2, such that indices and 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 and . 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 respectively (as described in Sec. III.1). 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:

 ∑x(^Vg)lx (^c)xrl′r′=∑x(^Vg)xl′ (^c)lrxr′,∑x(^V†g)rx (^c)lxl′r′=∑x(^V†g)xr′ (^c)lrl′x, (6)

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 of the bulk lattice , where the probability amplitudes of are (formally) obtained by contracting all its bond indices, analogous to how the MERA encodes the state . Thus, we have ‘lifted’ the MERA representation of a quantum state of the 1D lattice to a quantum state of the 2D lattice . We refer to the bulk tensor network, comprised of copies of the ground state tensors and the copy tensor , as the lifted MERA.

### iii.1 Local gauge symmetry

An immediate consequence of the symmetry conditions Eq. (6) is that the bulk state has a local gauge symmetry . Let us introduce gauge transformations on the bulk lattice as follows. The symmetry acts on the two sites located on a bond differently, namely, as and respectively. (This choice corresponds to the action of ‘left’ and ‘right’ gauge transformations in lattice gauge theory.) Elementary gauge transformations act on the 4 bulk sites (counting clockwise starting at the top left in the graphical representation) immediately surrounding tensor as , and on the 4 bulk sites immediately surrounding tensor as . General gauge transformations act on larger regions of the lattice by composing these elementary gauge transformations.

Let us consider the result of applying an elementary gauge transformation on a dual bulk state 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 of the bond tensors located immediately around an or tensor. The symmetry acts as and on the open indices (‘left’) and (‘right’) respectively. Owing to Eq. (6) [Fig. 3], operator 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 -symmetric [Eq. (4) and Fig. 1(b)], which eliminates the symmetry operators. Thus, the lifted MERA, and therefore state , 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 where 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 Elitzurâs 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 , which describes each site of the boundary lattice and also each site of the bulk lattice . Under the action of the symmetry, generally decomposes as

 V≅⨁jDj⊗Sj, (7)

where the symmetry acts on space by means of the irreducible representation (irrep) of labelled by quantum number (or charge) , and is the degeneracy space of irrep . Accordingly, the symmetry operators decompose as

 ^Vg=⨁j(^Idj⊗^Vg,j),   ∀g∈G. (8)

In particular, note that the symmetry operators act trivially on the degeneracy spaces, namely, as the identity on the degeneracy space , where is the dimension of the space .

We denote by an orthonormal basis in the irrep space , by an orthonormal basis in the degeneracy space , and by the basis on the total space . For example, if , then the symmetry charge is the total spin, is the spin projection along the -axis. (For simplicity, we assume that is multiplicity-free.) The description simplifies considerably for an Abelian symmetry, for example , since all the irreps of an Abelian group have dimension 1, that is, dim()=1 for all .

According to the Wigner-Eckart theorem, the -symmetric tensors and [Eq. (4)] decompose in terms of the intertwiners of . If the components of tensors and are denoted as and respectively, then in the irrep basis

 |p⟩≡|a,ta,ma⟩,   |q⟩≡|b,tb,mb⟩,|r⟩≡|c,tc,mc⟩,    |s⟩≡|d,td,md⟩, (9)

where denote symmetry charges, the tensors decompose as

 ^u≡⨁abcd;e (^u\tiny dege)abcd⊗(^Ie)abcd,^w≡⨁abcd;f (^w% \tiny degf)abcd⊗(^Jf)abcd, (10)

depicted in Fig. 5(a)-(b). Here

 (^Ie)abcd:(Sa⊗Sb)→(Sc⊗Sd)(^Jf)abcd:Sa→(Sb⊗Sc⊗Sd)

are 4-index intertwiners of the symmetry group , whose components are completely fixed by the properties of the group representations. (The intermediate charges and label a basis in a vector space of intertwiners.) The components of are given by [see Fig. 5(c)]

 [(^Ie)abcd]mambmcmd≡∑me⟨e,me|a,ma;b,mb⟩⟨e,me|c,mc;d,md⟩, (11)

where, for example, are the Clebsch-Gordan coefficients that describe the change of basis from the tensor product basis to the total charge basis . Analogously, we have [see Fig. 5(d)]

 [(^Jf)abcd]mambmcmd≡∑mf⟨f,mf|b,mb;cmc⟩⟨a,ma|f,mf;dmd⟩. (12)

Finally, and in Eq. (10) are degeneracy tensors, namely, multi-linear maps between the degeneracy spaces,

 (^u\tiny dege)abcd:(Da⊗Db)→(Dc⊗Dd)(^w\tiny degf)abcd:Da→(Db⊗Dc⊗Dd),

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 -symmetric tensors.

The copy tensor is -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 the trivial charge corresponds to the spin 0 irrep. Specifically, tensor decomposes as (see Fig. 5(e))

 ^c≡⨁a,a′ ^c\tiny dega,a′⊗(^I0)aa∗a′a′∗, (13)

where is a 4-index degeneracy tensor, denotes the conjugate charge of (namely, charges and fuse to the trivial charge), and is the intertwiner defined according to Eq. (11) for the trivial intermediate charge . The intertwiner is, in fact, equal to the tensor product of the identity and the identity on the irrep spaces and respectively [Eq. (7)].

Denote the components of a degeneracy tensor by where . The only non-zero components are given by

 (^c\tiny dega,a)tatatata=1. (14)

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 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 , 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 , 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.4.) Thus, in the context of holography, the bulk state may be interpreted as an entangled state of gauge fields (described by the spin networks) 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 ’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 , 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 -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.1 below. The second and third applications are explored in Sec. V.

### iv.1 Gauge-invariant bulk operators

As mentioned above, the spin network decomposition of the lifted MERA allows one to introduce gauge-invariant 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 gauge-invariant (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 symmetry. (The discussion readily generalizes to symmetry.)

### iv.2 Some properties of the bulk state

The boundary state is recovered from a bulk state by projecting every pair of sites located on a bond to the state defined as

 |+⟩≡∑a,ta,ma|a,ta,ma⟩⊗|a∗,ta∗,ma∗⟩.

State is isomorphic to the identity matrix after the identification . Thus, applying the projector on the two bulk sites located on bond 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 projector eliminates the copy tensor located on bond 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 .

It is readily checked that the -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 , 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 , 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 , one can construct a local, gauge-invariant bulk Hamiltonian whose ground state is , 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 where 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.

### v.1 Many bulk states dual to a ground state

Given a MERA representation of a ground state, one can obtain another equivalent MERA representation of the state by inserting a resolution of identity on bond , and multiplying the matrices and with the two tensors that are connected by the bond respectively. The two MERAs are an equivalent representation of the ground state, since the expectation value of any observable is the same in both the representations. (Obtaining an expectation value from the MERA involves contracting all the bond indices, and is multiplied with in the process.)

Clearly, inserting the copy tensor, defined according to Eq. (14), selects out a particular MERA representation of the ground state—the one expressed in a bond basis in which the degeneracy tensors have these components. On the other hand, the degeneracy copy tensors ‘commute’ only with diagonal matrices. Namely, a contraction of 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 -symmetric and isometric tensors. While -symmetric tensors ensure that the bulk state—obtained by lifting the MERA by inserting copies of the -symmetric copy tensor—has a local gauge symmetry (as described in Sec. III.1), the choice of isometric tensors leads to the desirable bulk features listed - in Sec. IV.2.

To this end, we restrict to unitary matrices that commute with the symmetry, namely, for all . Since decomposes as Eq. (8), Schur’s lemma (a special case of the Wigner-Eckart decomposition) implies that matrix decomposes as . Thus, the bond transformations are restricted to act as the identity 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 , in the different MERA representations of . Specifically, we fix the irrep basis 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 transformations (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 -symmetric and isometric tensors. (This corresponds to restricting the intrinsic bond transformations to unitary matrices that commute with the symmetry .)

### v.2 Critical spin chains

To this end, we considered the ground states of the following 1D critical spin models:

 ^H\small{ising}=∑i^σix^σi+1x+^σiz,^H\small{bc}=∑i−^Six^Si+1x+α(^Six)2+β^Siz,^H\small{potts}=−∑i^Pi(^PT)i+1+(^PT)i^Pi+1+^Mi,^H\small{xxz}=∑i^σix^σi+1x+^σiy^σi+1y+Δ^σiz^σi+1z, (15)

where labels sites of a 1D infinite lattice on which the Hamiltonian acts, are Pauli matrices, the operator is the component of the spin representation of , and and are Potts matrices:

 ^P=⎛⎜⎝010001100⎞⎟⎠;^M=⎛⎜⎝2000−1000−1⎞⎟⎠

The Blume-Capel model is critical for , and the XXZ model is critical for . The central charges and total symmetry groups of these models are listed in Table 1.

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 -symmetric MERA representation for the ground state of the Ising model and the Blume-Capel model, and a -symmetric MERA representation of the ground state of the Potts model.

For the XXZ model, we obtained both a -symmetric MERA representation of the ground states for (corresponding to a subgroup symmetry), and also an -symmetric MERA representation for . The symmetry of the XXZ models with corresponds to the total symmetry for and a subgroup symmetry for .

The representation of the symmetry on each site of the lattice for the various models is listed in Table 2. Table 3 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 3 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 and the relative error in the estimated central charge was . For the remaining models, the error in the estimated ground state energy density was at most , and the relative error in the estimated central charge was at most . For all the models, the relative error in the estimated smallest six scaling dimensions was between to .

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.)

### v.3 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 presented in Ref. (36) and references contained therein.

For the ground state of each of the models listed in Eq. (15), we randomly selected dual bulk states (restricting the corresponding bond transformations to unitary matrices that commute with the respective symmetry), and computed the second Renyi entanglement entropy per site,

 R\tiny tower≡−log2(Tr(^ρ\tiny tower% )2). (16)

Here 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 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 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 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 . Specifically, for ground states of the two XXZ models with the maximum, mode and minimum Renyi entropies obtained from the -symmetric MERA representation were found to be smaller than those obtained from the -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 -symmetric and a -symmetric MERA representation of a given ground state are expected to correspond to two bulk states with different entanglement respectively. This is because a -symmetric MERA representation can be converted to a -symmetric MERA representation by applying bond transformations to change the bond basis from a irrep basis to the subgroup basis listed in Table 3, 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.

### v.4 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 and 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 in Fig. 9. Let denote the reduced density matrix of region by tracing out all bulk sites outside , and also the gauge degrees of freedom inside . We computed the entanglement negativity given by

 n(^ρ[D])=∑i(|λi|−λi)/2, (17)

where are the eigenvalues of the matrix obtained by taking the partial transpose of with respect to some of the sites in . (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 for both these bulk states. These values are listed in Table 4. 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 , obtained by tracing out all bulk sites except the bulk site and also tracing out the gauge degrees of freedom on the site . (That is, has support only on the non-gauge sector of site .) 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 non-gauge (‘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 must commute with this symmetry. Consequently, by applying Schur’s lemma, must decompose as

 ^ρ\tiny degk=⨁j(^ρ\tiny degk,j⊗^Iηj). (18)

Here is an irrep of the emergent symmetry, is a density matrix that acts on the degeneracy space of charge , and is the identity matrix that acts on the irrep . The spectrum of is clearly degenerate, in accordance with this decomposition, specifically, the degeneracy of an eigenvalue of is at least .

In the scenario of an emergent symmetry, the degeneracies in the spectrum of 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 . 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 -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 description 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 . 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 and a 2D many-body state with a local symmetry , which may also be useful in characterizing and relating together different types of quantum phases of matter as illustrated in Appendix A.

Acknowledgements.—This research was largely completed while SS was employed at the Center for Engineered Quantum Systems in Macquarie University. We thank Guifre Vidal, John Baez, Sundance Bilson-Thompson, Yichen Shi, Giandemenico Palumbo and Rob Pfeifer for stimulating discussions. SS acknowledges the hospitality of the Perimeter Institute for Theoretical Physics where a part of this work was presented. We acknowledge support from the ARC via the Centre of Excellence in Engineered Quantum Systems (EQuS), project number CE110001013 and from DP160102426.

## 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 symmetry. (The discussion can be readily generalized to 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 -symmetric and isometric tensors does not have a non-trivial topological order.

### a.1 Example using Z2 symmetry

For the purpose of this section, we specialize the notation introduced in Sec. II to the case of a symmetry. Let here denote an infinite 1D lattice, each site of which is described by vector space and is equipped with the action of the group . The group acts on the space by means of the unitary representation . Under the action of the symmetry, the space decomposes as

 V≅Ve⊕Vo,

where and are the two irreps of . We denote by and a basis in the one dimensional vector spaces and respectively.

Let denote the (unnormalized) GHZ state belonging to the lattice ,

 |Φbound⟩≡|++⋯⟩+|−−⋯⟩, (19)

where and e.g. . State has a global symmetry since .

Consider a -symmetric MERA representation of comprised of copies of two simple tensors

 ^uGHZ:V⊗V→V⊗V,   ^wGHZ:V→V⊗V⊗V, (20)

which replace copies of the tensors and in Fig. 1 respectively. Tensor is simply the identity,

 (^uGHZ)ijkl=δikδjl,   i,j,k,l∈{e,f}, (21)

and the components of are:

 (^wGHZ)ijkl={1,if (i+j+k+l)%mod2=00,otherwise. (22)

Note that tensor is an isometry satisfying

 ∑ijk(^wGHZ)ijkl(^wGHZ∗)jkli′=δii′. (23)

It is readily checked that tensor is also -symmetric, since .

In order to verify that the MERA tensor network indeed represents the state , Eq. (19), we need to to contract all the tensors of to obtain the probability amplitudes , where denote the open indices of . The simple tensor network can be contracted algebraically to obtain

 ^Φi1i2⋯={1,if ∑kik is% even,0,otherwise. (24)

That is, is an (unnormalized) equal superposition of kets