arxiv: v2 [cond-mat.str-el] 23 Sep 2017

Transcription

1 Holograpic encoding of universality in corner spectra Cing-Yu Huang, 1 Tu-Cie Wei, 1 and Román Orús 2 1 C. N. Yang Institute for Teoretical Pysics and Department of Pysics and Astronomy, tate University of New York at tony Brook, NY , UA 2 Institute of Pysics, Joannes Gutenberg University, Main, Germany arxiv: v2 [cond-mat.str-el] 23 ep 2017 In numerical simulations of classical and quantum lattice systems, 2d corner transfer matrices (CTMs and 3d corner tensors (CTs are a useful tool to compute approximate contractions of infinite-sie tensor networks. In tis paper we sow ow te numerical CTMs and CTs can be used, additionally, to extract universal information from teir spectra. We provide examples of tis for classical and quantum systems, in 1d, 2d and 3d. Our results provide, in particular, practical evidence for a wide variety of models of te correspondence between d-dimensional quantum and (d + 1-dimensional classical spin systems. We sow also ow corner properties can be used to pinpoint quantum pase transitions, topological or not, witout te need for observables. Moreover, for a ciral topological PEP we sow by examples tat corner tensors can be used to extract te entanglement spectrum of alf a system, wit te expected symmetries of te U(2 k Wess- Zumino-Witten model describing its gapless edge for k = 1, 2. We also review te teory beind te quantum-classical correspondence for spin systems, and provide a new numerical sceme for quantum state renormaliation in 2d using CTs. Our results sow tat bulk information of a lattice system is encoded olograpically in efficiently-computable properties of its corners. I. INTRODUCTION Corner Transfer Matrices (CTMs were introduced by Baxter in te context of exactly solvable models in 2d. In is 1968 paper [1] e laid, witout noticing it, some of te basics of CTMs, togeter wit tose of te density matrix renormaliation group (DMRG and matrix product states (MP, wen dealing wit te dominant eigenvector of a 1d transfer matrix. CTMs are a key ingredient in te exact solution of several statistical-mecanical models [2], and ave also inspired many advances in te study of quantum many-body entanglement [3 5]. CTMs ave also been important for te numerical simulation of lattice systems, bot classical and quantum. In retrospect, Baxter proposed in 1978 a variational metod over CTMs [6], inspired by an earlier numerical metod from 1941 by Kramers and Wannier [7]. Tis, in turn, was one of te inspirations of Nisino and Okunisi s CTM Renormaliation Group metod (CTMRG [8]. (For te avid reader, a good source of information about tis istory can be found in Ref.[51]. imilar numerical CTM tecniques are also currently used in te calculation of low-energy properties of infinite-sie quantum lattice systems in 2d [4, 9, 10], for wic tey ave become one of te standard tools in te approximate calculation of pysically-relevant quantities suc as expectation values of local observables and low energy excitations. CTMs and teir algoritms ave also been generalied to 3d by te so-called Corner Tensors (CTs [5, 11], in turn allowing to explore iger-dimensional systems wit Tensor Network (TN metods. till, CTMs and CTs contain a great amount of olograpic information about te bulk properties of te system wic, a bit surprisingly, as not yet been fully exploited in te context of numerical simulations. Apart from being a useful object in te calculation of observables, te corner objects also contain, by temselves, information about te universal properties of te simulated model, providing a nice instance of te bulk-boundary correspondence for Tensor Networks (TNs [12]. Bulk information is encoded olograpically at te boundary corners, in a way similar to te study of te so-called entanglement spectrum and entanglement Hamiltonians [13]. For instance, Pescel, Kaulke and Legea [3] sowed tat te entanglement spectrum of a quantum spin cain (w.r.t. a partition into two semi-infinite segments is identical, up to some normaliation constant, to te spectrum of some CTM in 2d, wic could be computed exactly in some cases. Tis was te case of te Ising and Heisenberg quantum spin cains in a transverse field, for wic tey were able to compute suc a entanglement spectrum exactly as eigenvalues of a corner Hamiltonian, wic ere we call corner energies. Neverteless, and in spite of tese results, te study of te pysical information encoded olograpically in CTMs and CTs as been traditionally overlooked in numerical simulations, especially in te case of 2d quantum lattice systems, in spite of te fact tat tis is, indeed, a quite natural ting to do. In tis paper we explore te fingerprints of universal pysics tat are encoded olograpically in numerical CTMs and CTs. We do tis by studying te eigenvalue spectra of tese objects or, more precisely, of contractions of tese objects, togeter wit its associated entropy, in a way to be explained later. We provide several examples of tis bot for classical and quantum systems, including classical and quantum Ising, XY, XXZ and N-state Potts models, as well as several instances of 2d Projected Entangled Pair tates (PEP [14] describing perturbed Z 2, Z 3, symmetry-protected, and ciral topological orders [15 22]. To acieve tis goal we use a variety of TN metods for CTMs and CTs. For te case of ground-state properties of a quantum Hamiltonian H q in d dimensions, we set up a corner metod for a d + 1 dimensional

2 2 TN as described in Ref. [5] via te imaginary-time evolution operator e τhq for large enoug τ. From a broad perspective, some of our results can be understood as a generaliation of te work by Pescel, Kaulke and Legea [3] to 2d quantum systems. Additionally, wenever we ave direct access to te ground-state wavefunction ψ G in te form of a TN (e.g., a PEP, we can also study te CTMs originating from te TN for te norm ψ G ψ G, wic can be regarded as te partition function of some fictitious 2d classical model wit complex weigts. Trougout tis paper we sall refer to tis setup as reduction CTM (rctm, since it is a sceme tat reduces te wavefunction to a partition function. uc CTMs are, in fact, readily available in several TN algoritms (suc as te full update and fast full update for infinite PEP [9, 23]. Along te way, we also compare different scemes for te classical-quantum correspondence, and provide some pedagogical derivations. Wen te quantum state ψ G is explicitly given by a TN, we can directly obtain its associated CTs. To do tis we propose a new sceme for quantum state renormaliation. In tis case, te entanglement spectrum of a partition (of infinite sie can be readily obtained by diagonaliing a contraction of CTs, as we sall explain. First we use te Ising PEP in te disorder pase as an example to demonstrate ow to obtain te CTs entanglement spectrum. Ten we also apply tis quantum state renormaliation to two cases of ciral topological ordered states, wit U(2 k edge modes (for k = 1, 2, and find te degeneracy pattern in te entanglement spectrum matces tat in te corresponding conformal tower for te vacuum of te U(2 k WZW model. Our work is organied as follows: in ec.(ii we provide a reminder on CTMs, CTs, some of teir properties, as well as a summary of previous relevant results. In ec.(iii we provide a summary of te TN numerical metods used to study te 1d, 2d and 3d classical and quantum lattice systems explored in tis paper. Moreover, we also provide a new numerical sceme for quantum state renormaliation in 2d using CTs. In ec.(iv we analye, as a first test, several models in te universality class of te quantum Ising spin cain in a transverse field. In ec.(v we sow ow te quantum-classical correspondence can be identified from corner properties, for 1d quantum vs 2d classical and 2d quantum vs 3d classical models. In tis section we also review te teory beind several approaces for te quantum-classical correspondence, namely, te partition function approac, Pescel s approac, and uuki s approac for te XY model [24]. Ten, in ec.(vi we provide furter examples were te calculation of corner properties is useful. In particular, we sow ow corner properties can be used to pinpoint pase transitions in quantum systems almost for free in common tensor network numerical algoritms, witout te need to compute observables explicitly. We sow tis for several PEP wit topological order, including symmetry-protected, as well as for te 2d XXZ model. In ec.(vii we sow ow CTs can be used to compute te entanglement spectrum of several bipartitions of an infinite 2d system. In particular, we apply te idea to ciral topological PEP [20], sowing tat te obtained spectra encode te expected symmetries of te ciral conformal field teory (CFT describing its gapless edge, specifically, U(2 k WZW models for k = 1, 2. Finally, in ec.(viii we wrap up wit a summary of te results, conclusions and perspectives. II. CORNER OBJECT A. Corner transfer matrices CTMs are objects tat can be defined for any 2d tensor network. Here, for simplicity, we assume te case of a 2d TN on a square lattice. uc a TN could be, e.g., te partition function of a classical lattice model, te time-evolution of a 1d quantum system, or te norm of a 2d PEP. To define wat a CTM is, we notice tat te contraction of te 2d TN can be obtained, at least teoretically, by multiplying four matrices C 1, C 2, C 3 and C 4, one for eac corner (see Fig. 1a. Terefore, one as tat Z = tr (C 1 C 2 C 3 C 4, (1 were Z is te scalar resulting from te contraction. Matrices C 1, C 2, C 3 and C 4 are te Corner Transfer Matrices of te system. Tey correspond to te (sometimes approximate contraction of all te tensors in eac one of te four corners of te 2d TN. In some cases, wen te appropriate lattice symmetries are present, te four CTMs are equal, i.e., C C 1 = C 2 = C 3 = C 4. For te sake of simplicity, in tis section we sall assume tat tis is te case, toug in te following sections te four CTMs are different wen computed numerically. It is also convenient to define diagonal CTMs C d = P CP 1. Depending on te symmetries of te system (and tus of C, matrix P may be arbitrary, unitary or ortogonal. Let us call te eigenvalues ν α, wit α = 1, 2,..., χ, and χ te bond dimension of te CTM. Ten, te contraction of te full TN reads Z = tr ( Cd 4 χ = να. 4 (2 α=1 In fact, one can understand tis as te trace of te exponential of a corner Hamiltonian H C, i.e., wit Z = tr ( e H C, (3 H C log ( C 4 d. (4 Notice tat a similar Hamiltonian can also be defined individually for eac one of te corners. Depending on te symmetries of te CTMs, H C may be a Hermitian operator or not. From te point of view of

3 3 Z C 1 C 2 C 4 C 3 C 1 C 2 C 4 C 3 FIG. 1: [Color online] Te contraction of a 2d square lattice of tensors results in a scalar Z, understood as te trace of te product of four CTMs, one for eac corner. A reduced density matrix ρ of a system wit a CTM at every corner. quantum states of 1d quantum lattice systems, it is well known [5] tat operator e H C is related to te reduced density matrix of alf an infinite cain (wit H C Hermitian in tis case, see Fig. 1b. In fact, te spectrum of cmidt coefficients λ α of alf an infinite quantum cain in its ground state is given by λ α = να. 2 Tese cmidt coefficients are related to te eigenvalues of te reduced density matrix of alf an infinite quantum system (te so-called entanglement spectrum [13] by = λ 2 α = να, 4 wic are known to codify universal information about te system wen close enoug to criticality [3]. In terms of, te contraction of te 2d TN reads Z = χ α=1. Aditionally, te eigenvalues ε α of te corner Hamiltonian H C read ε α log. (5 In tis paper we call tese eigenvalues ε α s corner energies. B. Corner Tensors imilarly to CTMs for 2d TNs, one can define corner objects for iger dimensions, wic we generically call Corner Tensors (CT. Formally speaking, a CT is te (sometimes approximate contraction of all te tensors at one of te corners of a TN. For instance, for a TN on a 3d cubic lattice, one would ave tat its contraction Z is equivalent to te contraction of eigt CTs, i.e., Z = f(c 1, C 2, C 3, C 4, C 5, C 6, C 7, C 8, (6 FIG. 2: [Color online] 3d corner tensors wic correspond to tracing over, respectively, tree quarters and alf of a given 2d quantum system. wit C i (i = 1,..., 8 eigt tree-index tensors (te CTs, and f( a function specifying te contraction pattern, see Fig. 2. For te case of systems wit CTs it is also possible to define corner Hamiltonians. For instance, contractions suc as te ones in Fig. 2 correspond, for te case of a 2d quantum lattice system, to tracing over tree quarters or alf of te infinite system. For quantum systems described by a 2d PEP, it is possible to obtain tese types of contractions by using te quantum state renormaliation sceme from ec.(iii. In suc cases, tese contractions correspond to te reduced density matrices ρ of eiter one quarter or alf an infinite 2d system, wit eigenvalues, α = 1,..., χ (entanglement spectrum. Te contraction of te full 3d TN tus amounts to Z = χ α=1, as in te lower-dimensional case of CTMs. Again, it is possible to define a corner Hamiltonian H C and corner energies ε α in an analogous way as for CTMs. C. Previous results CTMs and CTs ave proven to be important in a variety of contexts, bot for teory and numerics. In statistical mecanics tey were used to solve te ard exagon model and many oters [1, 2]. From te perspective of quantum information, it is well known tat te corner Hamiltonian H C is related to a quantum system wic, in some cases, can be diagonalied exactly [3]. Numerically, Baxter developed a variational metod to approximate te partition function per site of a 2d classical lattice model by truncating in te eigenvalue spectrum of te CTM [6]. Tis was later refined by Nisino and Okunisi, wo developed te Corner Transfer Matrix Renormaliation Group metod (CTMRG [8]. Alternative truncation scemes for CTMRG ave also been studied, based on a directional approac and wit a direct application in infinite-pep algoritms [4, 23]. In fact, CTMs ave been applied extensively in te calculation of effective environments in infinite-pep simulations [25]. Moreover, tey ave been used as well in te generaliation to 2d

4 4 C1 i= T1 T4 C 1 T2 T3 APPROACH AND METHOD A. (e C1 C1 U W0 T4 = T4 = C4 T4 U C 4 (d T3 C4 U W 00 C 4 C4 (f C 1 of te time-dependent variational principle [10], wic is also useful in te calculation of 2d excitations. As for generaliations, CTMs ave also been used in oter 2d geometries, including lattice discretiations of Ad manifolds [26]. Numerical metods wit CTMs were also implemented in systems wit periodic boundary conditions [27] as well as stocastic models [28]. Metods targeting directly te corner Hamiltonian ave also been considered [29, 30]. Finally, te iger-dimensional generaliation to corner tensors as also been used to develop new numerical simulation algoritms [5, 11]. = T 4 C3 FIG. 3: [Color online] 2d PEP on a square lattice and its renormalied version wit CTs III. T1 C1 C 1 U (c T 4 C4 W C2 C 1 + T = T = T C 4 (C 4 T C 1 (C 1 U DU U (g C 1 C 1 + W DW W FIG. 4: [Color online] 2d quantum state renormaliation wit corner tensors: a left move, were one column is absorbed to te left. Te procedure is te same as in te directional CT approac from Ref. [5], but on a single layer of PEP tensors instead of two layers. Consequently, at every step we need to renormalie wit isometries not just te bond indices, but also te pysical indices, wic proliferate at every iteration. everal prescriptions are possible for te calculation of te isometries, e.g, one could consider iger-order singular value decompositions of te resulting tensors [34], or compute te reduced density operators of te indices to be truncated [8]. Generalities In te following sections we sall sow ow te spectrum of eigenvalues ωα, or equivalently te spectrum of corner energies εα, encodes useful universal information wen computed numerically for a variety of classical and quantum lattice systems. Tis is also true for te corner entropy given by X ωα log ωα. (7 α In particular, we will sow explicitly ow te spectrum as well as te entropy exactly coincide if compared between some d-dimensional quantum and (d + 1-dimensional classical spin systems, as expected from te quantumclassical correspondence. Moreover, we will also study tem for a variety of oter models, including several instances of topologically-ordered states. We will see tat tis can be useful to pinpoint pase transitions as well as to study edge pysics of ciral topological states. Concerning numerical algoritms, in our simulations we ave used te following, depending on te nature of te system to be studied: 1. For 1d quantum: te infinite Time-Evolving Block Decimation (itebd [31] to approximate ground states. Te spectrum ωα obtained from CTMs is easily related [3] to te cmidt coefficients λα of a bipartition, readily available from itebd or idmrg [32], as ωα = λ2α. In some instances we also use te simplified one-directional 1d metod from Ref. [5]. 2. For 2d classical: 2d directional CTM approac [4]. 3. For 2d quantum: if a quantum Hamiltonian is given, ten we use te 3d directional CTM approac [4] to compute properties of CTs, as well as infinite-pep (ipep [23] to approximate ground states. If te ground state ψg i is given, ten we use te directional CTM approac for te doublelayer tensors of te norm [4] to compute te re(r duced spectrum ωα from rctm. Moreover, we also use te 2d quantum state renormaliation described in te next section to compute properties of CTs. As we sall see, tis metod is single-layer and targets directly te quantum state. 4. For 3d classical: metod [5]. B. simplified one-directional 2d 2d quantum state renormaliation wit CTs Te procedure of quantum state renormaliation is important in 2d to obtain te contractions from Fig. 2 in te quantum case, wic give te reduced density matrix by tracing spins in tree quadrants or alf-infinite plane. Te entanglement spectrum can ten be obtained

5 5 from te eigenvalues of suc reduced density matrix. We ave implemented our own approac for te case of a 2d PEP, using CTs and single-layer contractions. Tis procedure, wic is an independent algoritm by itself, is explained in detail in wat follows. Te quantum state renormaliation group (QRG transformation acts directly on a quantum state and aims to extract a fixed-point wave function encoding universal properties [33]. Te basic idea is to remove non-universal sort-range entanglement related to te microscopic details of te system. After many rounds of QRG, te original ground state flows to a simpler fixed-point state, from wic one can identify to wic pase te system belongs to. In order to determine te fixed-point wave function we make use of CTs, see Fig. 3. Te distinction from te usual QRG is tat ere te fixed-point wave function will be encoded in tese CTs. Te procedure is similar to te directional CTM approac from Ref. [4], but tis time acting directly on te PEP, wic is single-layer, and not on te TN for te norm, wic is double-layer. An example of a left-move is in Fig. 4, were we sow also a simple option to obtain te isometrics needed for te coarse-grainings. We follow tis procedure by absorbing rows and columns towards te left, up, rigt and down directions until convergence is reaced. In te end, te corner tensors C represent te renormaliation of one quadrant of te 2d PEP, and te alf-row/alf-column tensors T to te renormaliation of alf an infinite row or column of tensors in te PEP. One ten follows te contractions in Fig. 2 to obtain te corresponding reduced density matrix and ence te entanglement spectrum. one-directional 1d metod from Ref. [5] to obtain te corner spectrum. As argued in Ref. [3] we expect and verify tat {λ 2 α} agrees wit { }. (ii 2d classical Ising: te partition function is given by Z c = {s} wit classical Hamiltonian e βhc({s}, (9 H c {s} = i,j s s [j], (10 were β = 1/T is te inverse temperature, s = ±1 is a classical spin variable at site i, {s} is a spin configuration, and te sum in te Hamiltonian runs over nearest neigbours on te square lattice. Te model is exactly solvable, and te critical point satisfies β c = 1 2 log ( It is well known [31] tat te partition function Z c can be written as an exact 2d tensor network wit tensors on te sites of a square lattice. Te approximate contraction is terefore amenable to tensor network metods. We use te directional CTM approac to compute te corner spectra and corner entropy from te tensors defining te partition function of te model. (iii 2d Ising PEP: as explained in Ref. [22], it is actually possible to write an exact Projected Entangled Pair tate (PEP [14] wit bond dimension D = 2 wose expectation values are te ones of te 2d classical Ising model. Te way to construct tis PEP is simple: one starts by considering te quantum state IV. FIRT TET: THE 1D QUANTUM IING UNIVERALITY CLA ψ(β = 1 Z c e ( β 2 i,j σ σ[j] +, +,, +, (11 In order to build some intuition about te numerical information contained in te spectrum ε α of corner energies, we ave first performed a series of numerical tests in systems belonging to te universality class of te 1d quantum Ising model in a transverse field. Te analyed models undergo a 2nd order quantum or classical pase transition, wit te critical point being described by an effective (1 + 1-dimensional CFT of a free fermion [35]. Te models and metods considered are: (i 1d quantum Ising: te quantum Hamiltonian is given by H q = i σ x σ x [i+1] i σ, (8 wit σ α te corresponding α-pauli matrices at site i, and te transverse magnetic field, wit critical point at c = 1. We used itebd to approximate te ground state by a Matrix Product tate (MP [36] and ere te square of te cmidt coefficients λ 2 α (ence te entanglement spectrum is obtained. We also use te simplified wit β some inverse temperature and + te +1 eigenstate of σ x. It is easy to see tat te expectation values of tis quantum state matc te ones of te 2d classical Ising model, e.g., ψ(β σ σ [j] ψ(β = 1 s s [j] e βhc({s} = s s [j] β, Z c {s} (12 wit H c ({s} te classical Hamiltonian in Eq.(10, and β te expectation value in te canonical ensemble at inverse temperature β. For a square lattice, one can also see [22] tat te state ψ(β can be written exactly as a 2d PEP wit bond dimension D = 2. If A is te tensor defining te PEP, its non-ero coefficients are given by A A 0010 = (cos(β/24 = (cos(β/23 (sin(β/2 A = (cos(β/22 (sin(β/2 2 A 1110 = (cos(β/2 (sin(β/23 A = (sin(β/24 (13

6 6 and permutations tereof. In te above equations, te convention for te PEP indices is A i αβγδ, wit α, β, γ, δ te left, up, rigt and down indices, and i te pysical index (tis time in te +/ basis. By construction, tis PEP is critical at te same critical β c tan te classical Ising model, and belongs also to te same universality class. For te numerical simulations it is sometimes convenient to parametrie te PEP in terms of g = 1 2 arcsin(e β, and terefore g c For tis state, we computed te corner spectra and entropy from te double-layer TN defining its norm, using te directional CTM approac [4]. For tese tree models and te metods mentioned we ave computed te spectrum as a function of te relevant parameter (magnetic field, inverse temperature, perturbation..., as well as te corner entropy = α log. Te results are sown in Fig. 5. Te differences between models correspond to rescalings in te defining variables and parameters tat map te different models among tem. More specifically, we can rescale te parameters and g using te 2d classical Ising reduced temperature t = T/T c as te basic variable, wic is related to te magnetic field of te 1d quantum model by t = T c / arcsin 1/, and to te parameter g of te 2d Ising PEP by t = T c / log (sin(2g. As sown in te plots, in all cases one can see tat te entropy tends to ave te same type of divergence. Concerning te corner spectra, we see tat all te models reproduce te same type of brances on bot te symmetric and te symmetry-broken pases. As expected, all spectra matc perfectly between te different calculations, since te different models can be mapped into eac oter exactly. V. BENCHMARKING THE QUANTUM-CLAICAL CORREPONDENCE In tis section we consider te corner energies for a variety of quantum and classical systems, wic allows us to study in good detail te correspondence between quantum spin systems in d dimensions and classical systems in d + 1 dimensions. Tere are several approces and ere we focus mainly on tree of tem, wic we sall refer to as te partition-function metod [37], Pescel s metod [38], and uuki s metod [24], respectively. We will give pedagoical treatment, specialiing to a few models and sow numerical results for a variety of 1d and 2d quantum and 2d and 3d classical models. A. Partition-function approac We now review te standard procedure beind te partition-function approac for quantum-classical mapping and ten examine suc correspondence in terms of entanglement and corner spectra. Te main idea is tat, for a d-dimensional quantum Hamiltonian H q at inverse λ 2 α (r t( t( (c t (d t(g FIG. 5: [Color online] Entanglement spectra λ 2 α and te entanglement entropy obtained from itebd of 1d quantum Ising model wit parameter t( as te function of transverse field. Te corner spectra and te corner entropy of: also te 1d quantum Ising model wit parameter t( as te function of transverse field, but computed wit te simplified one-directional 1d metod [5]; (c 2d classical Ising model wit temperature t = T/T c; (d 2d quantum Ising PEP wit parameter t(g as te function of g. In (c,d te corner tensors are obtained from te rctm setting, see also examples in ec. VI. In all cases, te bond dimension of te CTMs - equivalent to te bond dimension of te MP in case - is χ = 40. temperature β, te canonical quantum partition function Z q = tr(e βhq can be evaluated by writing it as a pat integral in imaginary time, i.e., Z q = tr ( e βhq = m m e βhq m, (14 wit m a given basis of te Hilbert space. Introducing resolutions of te identity at intermediate steps in

7 7 imaginary time one as Z q = {m} m 0 U m L 1 m L 1 U m L 2 m 1 U m 0, (15 wit U e δτhq, δτ β/l 1 (smaller tan all time scales of H q, and were te sum is for all te configurations of m α, α = 0, 1,..., L 1 wit m L = m 0, i.e., periodic boundary condition in imaginary time. As suc, tis way of writing te partition function can be interpreted in some cases as te one of a classical model wit some variables m α along an extra dimension emerging from te imaginary-time evolution. In wat follows, we make tis specific for te quantum Ising and Potts models, and bencmark te teory wit numerical simulations using CTMs and CTs computing te corner spectra and corner entropy. 1. Transverse field quantum Ising model in d dimensions (i Mapping via te partition function: let us consider te quantum Ising model wit a transverse field in d dimensions for L spins. For convenience, we use now te following notation for its Hamiltonian: H q = J i,j σ σ [j] J x i σ x = H + H x, (16 were σ α is te αt Pauli matrix on site i, J is te interaction coupling, J x te field strengt, and te sum i, j runs over nearest-neigbors. Te canonical quantum partition function of tis model is given by Z q = tr ( e βhq = {η } e βhq {η }, (17 η wit {η } η [1], η [2],, η [L] te diagonal -basis of te N spins, so tat η = ±1, i = 1, 2,..., L. plitting te imaginary time β into infinitesimal time steps δτ we obtain {η (τ + δτ} e δτhq {η (τ} {η {η (τ + δτ} e δτhx e δτh (τ} = e δτh({η(τ} {η {η (τ + δτ} e δτhx (τ}, (18 were in te first line we performed a first-order Trotter approximation wit O(δτ 2 error. Next, we consider te term wit Hamiltonian H x. In te single-site -basis tis can be written as η (τ + δτ e δτjxσ x η (τ = η (τ + δτ e δτjxσ x η x η x η (τ η x =±1 = η x =±1 e δτjxη x η (τ + δτ η x η x η (τ. (19 We can now use te overlap relation ( η x η = 1 e iπ 2 and define η notation, we now ave η = e δτjxσ x η η x =±1 ( 1 η x 2 η (τ + δτ, η e δτjxη x 1 2 eiπ = 1 ( e δτjx + e δτjx η 2 = 1 ( 1 + e 2 eδτjx 2δτJx η 1 η 2, (20 η (τ. Using tis ( ( 1 η x 1 η η 2 η η. (21 Moreover, we ave te alternative representation η e δτjxσ x η Jτ η = Ce η ( = C cos(j τ + sin(j τ η ( = C cos(j τ 1 + tan(j τ η η η, (22 wit C a normaliation constant. Comparing Eqs. (21 and (22, we obtain te relation tan(j τ = e 2δτJx. Finally, te partition function Z q of te transverse-field quantum Ising model can be written as Z q {η} C e Js e Jτ α, i,j η (ταη[j] (τα α,i η (τα+1η (τα, (23 were te coupling constants along te imaginary-time (τ and space (s directions are given by J τ = tan 1 ( e 2δτJx J s = J δτ. (24 Terefore, te canonical quantum partition function of a d-dimensional quantum Ising model wit a transverse field at inverse temperature β can be approximately represented by te classical partition function of a (d + 1-dimensional classical Ising model of sie β in te imaginary-time direction. Te exact correspondence arrives if we take te number of sites L in te imaginary

8 Te case of a 3d classical Ising model on a cubic lattice, analogous to a 2d quantum Ising model in a transverse field on te square lattice, only introduces one more relation in additional to tose Eq. (27 for an extra coupling along a spatial direction. i.e., 8 (J K y K x (J s FIG. 6: [Color online] Coupling constants for a 2d classical Ising model. In connection wit te quantum-classical correspondence, te vertical direction corresponds to imaginarytime FIG. 7: Te diagonal transfer matrix of square lattice. time drection to be infinity, giving δ = β/l 0, and ten te corresponding classical model as te couplings J s 0 and J τ. In Monte Carlo simulations, tricks can be used to deal wit suc as a limit [39]. For our simulations using correspondence from suc a partitionfunction approac, we ave to take δ increasingly small to obtain te exact correspondence of te spectrum. Re-parametriing te derived classical 2d anisotropic Ising model (see Fig. 6 we ave βh c = ( K x s [i,j] s [i,j+1] + K y s [i,j] s [i+1,j], (25 i,j were K x, K y are respectively te oriontal and vertical couplings, s [i,j] = ±1 are classical spins at site [i, j], and te sum runs over nearest neigbors on a square lattice. Te classical canonical partition function of tis model is given by Z c = {s} D2 D1 e ( i,j Kxs[i,j] s [i,j+1] +K ys [i,j] s [i+1,j]. (26 Comparing Eq. (23 wit Eq. (26 we ten ave te relations K x = J s = J δτ, K y = J τ = tan 1 (e 2δτJx, (27 were we can set J = 1 and J x =. We tus obtain te relation between and K x, K y, tan K y = e 2Kx. (28 Te exact mapping is obtained in te limit K x 0 and K y 0. K x = J s = J δτ, K y = J s = J δτ, K = J τ = tan 1 (e 2δτJx. (29 uc a d-dimensional quantum Ising model is mapped to a corresponding (d+1-dimensional classical Ising model, wic as omogeneous couplings along d spatial dimensions, and is anisotropic in te extra (imaginary temporal dimension. (ii Pescel s mapping in 2d: In a work by Pescel [38], it was sown tat a 2d classical Ising model wit an isotropic coupling K is in exact correspondence to a 1d quantum spin cain wit Hamiltonian L 1 H q = i=1 σ x δσ [L] x L 1 λ σ σ [i+1], (30 i=1 were δ = cos 2K and λ = sin 2 K, by using a transfer matrix tecnique. Te transverse field labeled as δ at te rigt end can be neglected for large L. Ten one arrives at te usual omogeneous cain. Let us briefly review ow tis is derived. Consider te classical Hamiltonian of te 2d isotropic Ising model given by βh c = i,j K(s [i,j] s [i,j+1] + s [i,j] s [i+1,j], (31 were s [i,j] = ±1 is a classical spin at site [i, j] and β is te inverse temperature. Te partition function is given by Z c = {s} e (K i,j (s[i,j] s [i,j+1] +s [i,j] s [i+1,j]. (32 Firstly, by drawing te lattice diagonally (i.e., rotate te square lattice by 45 degrees, te sites can form a row as sown in Fig. 7, and tese rows can be classified into two types: open circles and solid circles. Tis means tat te number of rows must be even. Let now N be te number of rows and M is te number of sites in eac row. Moreover, let φ r denote all spins in row r wit 2 M possible values. In particular, te partition function can be represented by te diagonal-to-diagonal transfer matrix W and V as follows: Z c = (D 1 φ1,φ 2 (D 2 φ2,φ 3 (D 1 φ3,φ 4 φ 1 φ 2 φ N (D 1 φn 1,φ N (D 2 φn,φ 1. (33 Here, (D 1 φj,φ j+1 contains all Boltmann weigt factors of te spins (from open circles to solid circles in te

9 9 adjacent rows j and j +1. imilarly, (D 2 φj,φ j+1 contains te oter type of spins (from solid circles to open circles. We now consider tree rows labeled as φ, φ, φ, were φ = {s 1, s 2,..., s M } are te spins in te lower row and similarly for φ and φ. Ten te diagonal-to-diagonal transfer matrix is given by (D 1 φ,φ = e K( M j=1 (sj+1s j +sjs j, (D 2 φ,φ = ek( M j=1 (s j s j +s j s j+1. (34 Te partition function can tus be written as Z c = tr(d 1 D 2...D 1 D 2 = tr(d 1 D 2 N/2 = tr(v N/2. One can verify tat [H q, V ] = 0 if te couplings are cosen to satisfy δ = cos 2K, λ = sin 2 K. (35 If te lattice sie is large enoug, ten te single term σl x can be neglected. In tis case te Hamiltonian can be written as a 1d quantum Ising cain wit transverse field, H q /λ = L 1 i=1 σ x L 1 i=1 σ σ [i+1] wit = 1/λ = 1/ sin 2 K. It is wort mentioning tat te mapping is exact in te sense tat no limit in any parameter needs to be taken (in contrast to, e.g., te partition-function approac, were we ad δτ 0. (iii Numerical results: according to te mapping described above, we ave computed te corner spectra and te associated corner entropy for Ising models, first comparing te 1d quantum and 2d classical, and ten te 2d quantum and 3d classical, using te numerical tecniques mentioned earlier. On te one and, te comparison of 1d quantum vs 2d classical is sown in Fig. 8, were we also include in te second panel te mapping to te isotropic classical Ising model by Pescel [38]. Regarding te anisotropic classical model, te mapping becomes more and more precise as δτ 0, i.e. as K x becomes smaller. In our results, wen plotted wit respect to te same variables, we see a remarkably perfect agreement for all te numerical values of and among all te models. On te oter and, we sow in Fig. 9 our results comparing te 2d quantum vs 3d classical (anisotropic case. Te matc in tis case is not as perfect as in te 1d vs 2d case, but neverteless, it is still quite remarkable, especially considering te inner workings and associated errors of te iger-dimensional numerical algoritms tat we used. To understand furter te data obtained from te corners, we sow in Fig. 10 te corner entropies in more detail, as well as difference between te 2d classical corner entropy and te one for te 1d quantum case (wic equates te entanglement entropy. One can see in a more precise way tat te entropy in te classical anisotropic case tends to te quantum one as te coupling K x tends to ero, as expected from Eq. (27 for δτ 1. Te effect of a finite K x is better appreciated in Fig. 11, were one can see clearly ow te classical value tends to matc as a limiting case te quantum value as K x 0. Finally, in Fig. 12 we sow a comparison of te λ 2 α (K (c (K y (d (K y FIG. 8: [Color online] Entanglement spectra and entanglement entropy of te 1d quantum Ising model in a transverse field as obtained wit itebd. (b,c,d Corner spectra and corner entropy of: te 2d classical isotropic Ising model, as a function of, wit isotropic coupling K satisfying 1/ = sin 2 K; (c,d 2d anisotropic classical Ising model wit fixed K x = (c, K x = 1 (d, and K y as a function of satisfying tan K y = e 2Kx. Te corner bond dimension is χ = 20 in all cases. entropies for te 2d quantum vs 3d classical case. Again, as expected, te agreement between te quantum and te classical case improves as K x gets closer to ero. 2. Transverse field quantum N-Potts model in 1 dimension (i Mapping: we now consider te 1d quantum N-state Potts model in 1d for L sites. Te corresponding 1d

10 (K y (c (K y (d (K y FIG. 9: [Color online] Corner spectra and corner entropy of: te 2d quantum Ising model in a transverse field by using te simplified one-directional 2d metod [5]; (b,c,d 3d anisotropic classical Ising model (also wit te same metod wit fixed K x = K y =, K x = K y = 5 (c, and K x = K y = 1 (d, and K as a function of satisfying tan K = e 2Kx. Te corner bond dimension for te CTs is χ = 4 in all cases. quantum Potts Hamiltonian is given by L 1 H q = i=1 ( N 1 = H + H x, n=1 (Z Z [i+1] n L i=1 ( N 1 were operators Z and X at every site satisfy n=1 (X n (36 Z q = ω q q, X q = q 1, (37 wit ω = e i2π/n and q Z N. imilar to te case of te Ising model, te quantum 1d quantum 2d isotropic K x = K x =1 K x =01 FIG. 10: [Color online] Corner entropy of te 1d quantum Ising model, 2d classical isotropic Ising model (1/ = sin 2 K, and 2d classical anisotropic Ising model wit fixed K x =, K x = 1, and K x = 01 (tan K y = e 2Kx as a function of te transverse field wit bond dimension χ = 20. In te inset we sow te difference between te 2d corner entropies and te 1d entanglement entropy. canonical partition function is given again by Z q = tr ( e βhq = {η } e βhq {η }, (38 η but tis time {η } η [1], η [2],..., η [L] is te diagonal basis of Z for te L spins, so tat η = 0, 1, 2,..., N 1, i = 1, 2,..., L. Proceeding as for te Ising model in te previous section, now we ave a similar expression as in Eq. (18, but wit H and H x being te ones in Eq. (36. For te Hamiltonian term H we find η η [i+1] e δτ( N 1 = e δτϑ δ η η n=1 (Z Z [i+1] n η δ η [i+1] η [i+1] η [i+1], (39 were η η (τ + δτ and η η (τ. Te coefficient ϑ is ϑ = N 1 if η = η [i+1], and ϑ = 1 oterwise. Additionally, for te term H x one as η = e δτ( N 1 n=1 (X n η η N 1 cos(δτi + sin(δτ( (X n { cos(δτ = sin(δτ if η = η oterwise. n=1 η (40 For te classical case, te Hamiltonian of te 2d classical N-state Potts model on a square lattice is defined

11 K x d quantum K x = K x =5 K x = FIG. 12: [Color online] Corner entropy of te 2d quantum Ising model and 3d classical anisotropic Ising model wit fixed K x =, K x = 5, K x = 1 (tan K y = e 2Kx as a function of te transverse field wit corner dimension χ = 4. FIG. 11: [Color online] Corner entropy of te 2d classical anisotropic Ising model wit fixed (upper = 0.8 and (lower = 1.2 as a function of K x wit corner dimension χ = 20. Te blue dased lines sow te entanglement entropy of te ground state of te corresponding 1d quantum Ising model obtained by using te itebd metod. by βh c = i,j ( Kx δ s [i,j],s [i,j+1] + K yδ s [i,j],s [i+1,j], (41 wit Potts spin variables s [i,j] = 0, 1, 2,..., N 1 at eac site. Te classical partition function is ten Z c = {s} ( i,j e Kxδ s [i,j],s [i,j+1] +K yδ s [i,j],s [i+1,j]. (42 From Eqs. (39, (40, and (42 one finds te relations K x = Nδτ, tan(δτ = e Ky, (43 wic establis te quantum-classical mapping. (ii Numerical results: as we did for te case of te Ising model, now we ave bencmarked te quantumclassical correspondence by computing numerically te corner spectra and teir associated corner entropy for several quantum and classical Potts models. Our results are summaried in Fig. 13, were we sow te corner spectra and corner entropy for te 1d quantum and 2d classical N-state Potts models for N=2,3,4 and 5. Again, we find a remarkable almost-perfect matc for te corner properties as computed wit different metods for 1d quantum and 2d classical systems, once te parameters in te models are rescaled according to te relations found in te previous section. Te spectrums for te 2- state Potts model coincide wit tose of te Ising model, as expected. As N increases, we find small variations in te corner for different values of N, even toug te brances corresponding to te lowest corner spectra seem to be very similar for all te computed N. B. uuki s approac for te quantum XY model In a work by uuki [24] it was proven tat a 2d classical Ising model in te absence of a magnetic field and wit anisotropic couplings is equivalent, in te sense of aving te same expectation values and pysical properties, to te ground state of a XY quantum spin cain. Unlike te partition function approac, wic maps a quantum model to a classical model in one dimension iger, uuki s approac works from te oter direction: it starts from te (d + 1-dimensional classical partition function, and ten builds a d-dimensional quantum model wit te same pysics. We note tat te mapping is exact and does not involve te limit. However, if one uses te quantum XY model to study te transverse-field Ising model, ten a similar limit needs to be taken. Morever, it was known tat tere is a range of couplings in te quantum XY model tat tere is no valid classical correspondence (see te O region in Fig. 14. (i Te mapping: let us review te teory beind tis approac by considering first te classical Hamiltonian of te anisotropic 2d XY model, i.e., βh c = ( K x s [i,j] s [i,j+1] + K y s [i,j] s [i+1,j], (44 i,j were indices i, j denote respectively rows and columns,

12 λ 2 α (K y (c (K y (e (d (f ( (K y (g (K y 12 FIG. 13: [Color online] Entanglement spectra and entanglement entropy for te 1d quantum N-state Potts model in transverse field for N=2, (c N=3, (e N=4, and (g N=5 by using itebd metod. Corner spectra and corner entropy for te 2d classical N-state Potts model as a function of, were is a function of K y as in Eq. (43, wit K x = 1 for N=2, (d N=3, (f N=4, and ( N=5 computed wit te 2d directional CTM metod. Te corner bond dimension is χ = 20 in all cases. K x, K y are te oriontal and vertical couplings, s [i,j] = ±1 are classical spin variables at eac site, and β is te inverse temperature. For concreteness let us imagine tat we ave a finite periodic square lattice wit N M sites. Te canonical partition function is given by Z c = e (Kx i,j s[i,j] s [i,j+1] +K y i,j s[i,j] s [i+1,j]. (45 {s} Te first sum inside te brackets in te exponential is over oriontal edges, and te second over te vertical ones. Let N be te number of rows in te lattice and M te number of sites in eac row. Now let φ r denote all spins in row r, so tat φ r as 2 M possible values. Te partition function can tus be tougt of as a function of φ 1,...,φ N, and can be rewritten as Z c = φ 1... φ N T φ1,φ 2...T φn 1,φ N T φn.φ 1, (46 Here T φi,φ i+1 is te 1d transfer matrix of te system, wic contains all te Boltmann weigt factors of te spin in te adjacent rows. Let φ = {s 1 s 2,...s M } be te spins in a given row, and φ = {s 1, s 2,...s M } te ones in te following row. Ten te transfer matrix is given by T φ,φ = e (Kx i i sisi+1+ky sis i = e Kx i sisi+1 e i Ky sis i V 1 V 2. (47 Here V 1 can be decomposed as a product of 2 2 matrices, ( e K x e (V 1 si,s i+1 = Kx, (48 e Kx e Kx wic can also be written as (V 1 si,s i+1 = e Kx I + e Kx σ x = e Kx (I + e 2Kx σ x = (2 sin 2K x 1/2 e K x σx V 1 (i, (49 wit I te 2 2 identity matrix, σ x te x-pauli matrix, and were we define tan Kx e 2Kx (and tan K x e 2K x, as well as use te relation sin 2K x sin 2Kx = 1. Moreover, one as te 4 4 matrix ( (V 2 si,s j;s δsi,s i,s j δ i s j,s given by j (V 2 si,s j;s i,s j (+1, +1 (+1, 1 ( 1, +1 ( 1, 1 e Ky e = Ky e Ky e Ky = exp(k y σ i σ i+1 = cos K y I i I i+1 + sin K y σ i σ i+1 V 2 (i, i + 1. (50 It is clear tat te partition function is te trace of a matrix product, given by Z c = tr(v 1 V 2...V 1 V 2 = tr(v 1 V 2 N. (51 Tus, Z c can also be written as Z c = tr(v 1/2 2 V 1 V 1/2 2 N = tr(v N, (52

13 2 + 2 =1 1 F 13 Te 1d quantum Hamiltonian is a sum of two-body operators H q = i (i, i + 1. Tus, te commutator reads [V, H q ] = i [V, (i, i + 1] O P (/ =1 O F FIG. 14: [Color online] Pase diagram of 1d [40] and 2d [41] quantum XY model. Tere are tree pases: oscillatory (O, ferromagnetic (F, and paramagnetic (P. Te equation on top is te Barouc-McCoy circle [42] tat sets te boundary between te oscillatory and non-oscillatory ferromagnetic regions (wic is only a crossover. Te separation between F and P in is at = 1 and in te exact location is not known and only indicated scematically. or were and Z c = tr(v 1/2 1 V 2 V 1/2 1 = tr(v N, (53 V 1 = (2 sin 2K x M/2 e (K x m i=1 σi, (54 V 2 = e (Ky M i=1 σi σi+1. (55 Te next step is to sow tat V and te quantum Hamiltonian H q for te 1d quantum XY model can commute, and terefore ave common eigenvectors. Te usual XY quantum spin cain is defined by te Hamiltonian H q = ( J x σ x σ x [i+1] + J y σ y σ y [i+1] + σ, i i (56 were γ = (J x J y is te anisotropy, and te magnetic field. Te pase diagram of te model is well known [40] and is sketced in Fig. 14a. To prove tat te commutator of V and H q can sometimes be ero, we first define V 2 (i, i v 2 (i, i P = ( [v2 (i 1, iv 2 (i, i + 1v 2 (i + 1, i + 2 i V 1 (iv 1 (i + 1v 2 (i 1, iv 2 (i, i + 1v 2 (i + 1, i + 2, (i, i + 1] = 0. (57 Te last equality imposes a constraint on te couplings of te classical and quantum models in order for te commutator to vanis. One can see tat tis implies te relations between te couplings J y J x = e 4Kx, J x = 2e 2Kx cot(2k y, (58 wic make explicit te quantum-classical mapping. Importantly, for fixed, J x and J y, tese equations do not ave a real solution in te oscillatory pase of Fig. 14, so tat te mapping is only valid outside of tat pase. Finally, te mapping can also be extended easily to te 3d classical vs 2d quantum case, by considering a 2d omogeneous coupling K x = K y = K and adding an extra equation for K, i.e., J y J x = e 4K, J x = 2e 2K cot(2k. (59 (ii Numerical results: we ave explicitly cecked tis equivalence by computing numerically te corner spectra and te associated corner entropy for te quantum and classical XY models in 1d, 2d and 3d. For te 1d quantum vs 2d classical case, tis is sown in Fig. 15 for different values of te anisotropy in te quantum XY model. Te expressions in Eq. (58 ave only a real solution for K x, K y if te value of is outside of te oscillatory pase, as sown in te plots. We can see tat te agreement between te quantum and classical corner spectra and corner entropy is remarkably good, bot qualitatively and quantitatively, wit a sligtly larger error around te critical region = 1. Te comparison between 2d quantum vs 3d classical can be found in Fig. 16. Again in tis case te matc between te numerically-computed classical and quantum values is quite remarkable, considering te different numerical tecniques tat were used in tis case. VI. 2D CORNER PHAE TRANITION We now sow ow te study of corner properties can provide oter useful information wen studying a quantum or classical many-body system. In particular, we sow ow te corner spectra and corner entropy from 2d rctms (i.e., te CTMs obtained from te 2d TN for te

14 14 λ 2 α (K y (c (d (K y (e (f (K y FIG. 15: [Color online] Corner spectra and entropy for te 1d XY quantum spin cain, wit γ =, (c γ = 0.9, and (e γ = 0.99, togeter wit te corresponding anisotropic 2d classical Ising model, (d, (f, respectively, wit K x and K y being functions of as described in Eq.(58. Te corner bond dimension (equivalent to te MP bond dimension in te 1d quantum case is χ = 40 in bot cases. Te correspondence of parameters as a solution only for values of larger tan 0.85, (d, (f, and terefore te left and side of eac plot in te lower panel is empty. norm are useful in determining pase transitions witout te need to compute pysical observables. Te usual way to study quantum and classical pase transitions is troug te study of observables, wic ave specific properties at te transition point (e.g., te singular beavior of te observable. Te study of entanglement and correlations in many-body systems as sown us tat it is actually possible to study tese transitions from properties of te state only, suc as entanglement entropy, fidelities [43], entanglement spectra [13], and similar quantities. Following tis trend, in tis section we sow tat one can assess pase transitions from properties of te corners only, in particular te rctm tat we introduced in ec. I. Tis is very useful in te context of numerical simulations of, e.g., 2d quantum many-body systems, since suc corner objects are produced for free (e.g., in te infinite-pep metod wit a full or fast-full update [9, 23]. In wat follows we sow tree practical examples were pase transitions, bot topological and non-topological, can be clearly pinpointed by looking only at te corner objects. is given by H q = i,j ( σ x σ x [j] + σ y σ y [j] σ σ [j] i σ, (60 were as usual te sum i, j runs over nearest neigbors on te 2d square lattice, and is te anisotropy. In te large > 1 limit, it as been sown [44] tat a first-order transition takes place at some point 1 from a Néel pase to a spin-flipping pase. As te field increases furter, anoter pase transition at 2 = 2(1+ occurs towards te fully polaried pase. Here we consider te case wit = 1.5. We ave approximated te ground state of te model using te ipep algoritm wit simple update and bond dimension D = 2 [45], and ten computed te reduced corner spectra (r and entropy of te double-layer tensor defining te norm via te directional CTM approac, as a function of. Our results are sown in Fig. 17, were one can clearly see tat te two pase transitions are clearly pinpointed by te spectrum and te entropy. In particular, we observe te first transition appening at 1 1.8, and te second one at 2 = 5.0. A. 2d quantum XXZ model First we consider te 2d quantum XXZ model for spin- 1/2 on an infinite square lattice, under te effect of a uniform magnetic field along te -axis. Its Hamiltonian B. Perturbed Z N topological order Here we consider exact wavefunctions tat exibit topological pase transitions for Z 2 and Z 3 topological order.

15 FIG. 16: [Color online] Corner spectra and corner entropy of: 2d quantum XY model wit γ = 0.7 in a transverse field by using te simplified one-directional 2d metod [5]; te corresponding 3d anisotropic classical Ising model as a function of satisfying Eq. (58. Te corner bond dimension is χ = 4 in all cases. Te correspondence of parameters as a solution only for values of larger tan 1.45, and terefore te left and side of te plot in te lower panel is empty. (i 2d perturbed Z 2 Toric Code PEP: we consider te 2d PEP on a square lattice for te Toric Code ground state [22, 46], perturbed by a string tension g. Tis can be represented by a tensor A i,j,k,l αβγδ wit wit four pysical indices i, j, k, l = 0, 1 and four virtual indices α, β, γ, δ = 0, 1. Te coefficients of te tensor are given by A i,j,k,l i,j,k,l = { g i+j+k+l, if i + j + k + l = 0 mod 2, 0, oterwise. (61 Te norm of tis state can be described by a double-layer 2d TN on a square lattice, were at every site one as te tensor T ijkl ijkl T[ijkl], wit coefficients T[0000] = 1, T[1111] = g 8, T[0011] = T[0110] = T[1100] = T[1001] = g 4 T[0101] = T[1010] = g 4. (62 Parameter g is used to tune a crossover from a topological to a trivial pase. For g = 1 te state reduces to te ground state of te Toric Code model wit Z 2 topological order. For g = 0 it reduces to te polaried state (r FIG. 17: [Color online] Corner spectra (r for te norm of te numerical D = 2 PEP for te XXZ model in a field, at = 1.5, on te square lattice wit χ = 40, togeter wit te corner entropy computed from te corner spectra. 0, 0,, 0. Tere is a quantum pase transition between tese two pases wic, as sown in Ref. [47], occurs at g c One can see, moreover, tat te double tensor T consists of two copies of te partition function of te 2d classical Ising model in Eq. (10. In fact, one also finds te relation g = (sin(β 1/4, wit g te perturbation parameter of te Toric Code and β te inverse temperature of te Ising model. Bot models, terefore, belong to te same universality class. In tis case we ave implemented te directional CTM metod on te norm tensor T [4] to study te corner properties. Tis is sown in Fig. 18a, were one can see tat te corner spectrum and its associated entropy clearly pinpoint te quantum pase transition. (ii 2d perturbed Z 3 topological order: furtermore, we consider a 2d PEP wit Z 3 topological order under perturbations described by deformations {q 0, q 1, q 2 }. Te PEP is given by a tensor A i,j,k,l α,β,γ,δ wit wit four pysical indices i, j, k, l = 0, 1, 2 and four virtual indices α, β, γ, δ = 0, 1, 2, wit coefficients { q n i,j,k,l = 0 0 qn1 1 qn2 2, if i + j + k + l = 0 mod 3, 0, oterwise, (63 A i,j,k,l were n 0, n 1, n 2 means te number of te inner indices in 0, 1, and 2 respectively. We first study te case q 0 = 1, n 1 = 0, q 2 = g. In suc a case, te bond indices of te wavefunction live in an effective 2d Hilbert space spanned by 0 and 2. At g = 0, te remaining tensor represents a product state of all state 0. Terefore, te region near g = 0 is a trivial pase tat is adiabatically connected to a product state. At g > 0, te nonero components of