# Topological Order and Degenerate Singular Value Spectrum in Two-Dimensional Dimerized Quantum Heisenberg Model

## Abstract

We study the connection between topological order and degeneracy of the singular value spectrum by explicitly solving the two-dimensional dimerized quantum Heisenberg model in the form of tensor product state ansatz. Based on the ground state solution, we find non-zero topological entanglement entropy at the frustrated regime. It indicates a possible topological phase. Furthermore, we find that the singular value spectrum associated with each link in tensor product state is doubly degenerate only in this phase. Degeneracy of the singular value spectrum is robust against various types of perturbations, in accordance with our expectation for topological order. Our results support the connection among topological order, long range entanglement and the dominant degenerate singular values. In the context of tensor product state ansatz, the numerical evaluation of singular value spectrum costs far less computation power than the one for topological entanglement entropy. Our results provide a more viable way to numerically identify the topological order for the generic frustrated systems.

^{3}

## I Introduction

Investigation of topological orders for spin models is under intensive study in the recent years, partly due to its relation to the ground state of high Tc superconductor, the so called spin liquid phase (1), and partly because of the enlightenment of quantum information to characterize the ground states by quantum entanglement (2). These studies lead to new classification schemes (5); (6); (7); (3); (4); (8) of the quantum phases beyond the usual Landau-Ginzburg-Wilson paradigm. The new quantum phases such as the spin liquids are called the topological phases and cannot be characterized by a local order parameter. Instead, it could be characterized by the ground state degeneracy, quasiparticle statistics, existence of edge states, topological entanglement entropy (9); (10), and entanglement spectrum (11), i.e., eigen-spectrum of the reduced density matrix. The latter twos are motivated by the quantum information studies on the entangled nature of ground state, which could be efficiently approximated by the projected entangled pair state ansatz (12). One typical example for all the above characteristics of topological ordered state is the fractional quantum Hall states (13); (14).

For one-dimensional (1D) spin chain the quantum fluctuation could be more significant than in the higher dimensional ones so that it will melt the system into liquid, and thus becoming topological. The simplest example of the “topological phase” is the Haldane phase. The term with the quotation remark is to remind that the topological nature of the ground state is protected by some symmetry of the spin chain (3); (4). Once the symmetry is broken, the ground state wave function is equivalent to a trivial product state up to some local unitary transformation (15). The local unitary transformations can remove the short range entanglement, i.e., entanglement between the neighboring sites. In contrast, the long range entanglement is the entanglement among the distant sites, and cannot be removed by the local unitary transformations. Due to its robustness against the local operations, the long range entanglement causes the topological order of the nontrivial ground state.

On the other hand, for two-dimensional (2D) spin systems the topological phase is usually expected to exist in the frustrated systems since the large ground state degeneracy of the system is highly quantum and cannot have simple classical order. However, suffering from the lack of the reliable numerical tools in solving the 2D frustrated spin systems, only very few topological phases have been identified, except for some particularly exactly solvable models (16); (17); (18). The connection between the topological order and the quantum entanglement is less understood than in the 1D cases.

Despite of the recent progress in (5); (6); (7); (3); (4); (8) on the general scheme of classifying the topological orders caused by the long range entanglement, it is still not easy in practical to dynamically determine if exists a topological phase. It is then desirable to find a more numerical viable characterization for the topological order. Since the topological order is believed to be linked to the long range entanglement, we think the degeneracy pattern of the entanglement spectrum can be such an index. The intuitive picture is that the entanglement spectrum is the Schmidt values of the bi-partition of the whole system (4); (8), so the dominant degenerate singular values implies that the wave function is close to the maximally entangled state (Bell state). In this case, the two parts of the system is coherently long-range entangled.

The entanglement spectrum is easy to evaluate numerically for 1D spin systems in the ansatz of matrix product states (MPS) because it coincides with the Schmidt value spectrum of the bi-partition, and can be obtained by performing the singular value decomposition (SVD) for the matrix of the MPS. Indeed, the connection between entanglement spectrum and topological order is established for the MPS numerically (4). However, this is not the case for the 2D tensor product state (TPS) because the left and right states of the bi-partition are now the 1D spin chains (8), and the corresponding entanglement spectrum is not the same as the singular value spectrum. The latter is obtained by merging two neighboring 1-site tensors in TPS and then performing the SVD on the merging tensor. Despite that, these twos should be related, especially for the translationally invariant states. Therefore, if the entanglement spectrum can be used to characterize the 2D topological order, so does the singular value spectrum. Indeed, we will see this is the case in this work.

The main motivation of this paper is to establish such an connection in 2D among topological order, long range entanglement and degenerate singular value spectrum by numerically studying a particular frustrated 2D spin model with its ground state in the form of the TPS ansatz. In the meantime we explore the properties of the singular value spectrum under various types of the local unitary transformations to test its robustness, which is the main characteristics of the topological order.

The model we will study is the staggered dimer spin 1/2 model on the square lattice, and will be referred as the J-J’ model form now on. It is an interesting 2D spin model with various issues about quantum phase transition. For the antiferromagnetic -bond regime, i.e., , this model has been studied by using the perturbation theory (19); (20), exact diagonalization (21), the coupled cluster method (21), iPEPS (22) and quantum Monte Carlo (23); (24). In this regime, there exists a quantum phase transition by tuning with a critical point at (19), separating the classical Néel ordered phase and a finite-gap disordered phase (23); (24).

In this paper we are more interested in the regime of the ferromagnetic -bond, i.e., because the plaquettes are now frustrated. The quantum Monte Carlo method cannot be used for the frustrated system because of the sign problem. Instead, some mean-field based approximation methods such as the renormalized spin wave theory (20), exact diagonalization (21) and coupled cluster method (21) have been used to study this regime. These studies found a phase transition at separating the Néel phase from a helical phase for classical spins, and the maximal frustration occurs around region. However, due to the frustrations we suspect that a quantum spin liquid phase with topological order may appear around the region of maximal frustration. We indeed find this is the case from our study.

The paper is organized as follows. In the next section we describe the Iterative Projection method of solving the 2D TPS ground state (25), in the literatures this method is also referred as the “simplified update” (26). This is the 2D generalization of the infinite time evolving block decimation (iTEBD) (27). We also review the method of the the tensor-network renormalization group (TRG) (28); (29), which are used to find the overlap of the wave functions in a polynomially efficient way. Based on these, we can obtain the singular value spectrum and topological entanglement entropy, and also perform the quantum state renormalization. We also briefly review the above concepts. In section 3 we present the phase diagram of the J-J’ model by evaluating the order parameter, entanglement entropy, dimer strength and the ground state energy. In section 4 we evaluate the topological entanglement entropy for the J-J’ model and a toric code like state to characterize the topological phase. In section 5 we focus on various properties of the singular value spectrum and its degeneracy to establish its connection with the topological order. Finally, we conclude our paper in section 6.

## Ii The Model and the Method

The model we study in this paper is the staggered dimer model on square lattice, also known as the J-J’ model. It is a quantum Heisenberg antiferromagnetic model with two kinds of nearest-neighbor exchange couplings, and its Hamiltonian is

(1) |

where the summations over and represent sums over the nearest-neighbor bonds as shown in Fig.1. Each square lattice consists of three bonds (the thin ones) and one bond (the thick one). We fix the bonds to be antiferromagnetic, i.e., , and consider as the parameter of this model.

### ii.1 Solving the model with Iterative Projection and TRG

We will use the Iterative Projection algorithm (25) (a.k.a. simplified update (26)) to find the translationally invariant ground state of the system in the form of the tensor product state ansatz

(2) |

where labels the physical spin of dimension at site for , ’s are rank-five tensors with the bond indices , and is to sum over all indices of tensors. We call the bond dimension and the physical dimension. This ansatz is able to handle an area-law by construction, but that doesn’t mean that the true ground state of the system necessarily obeys it. Depending on the situation, we should tune to fully capture the essential feature of the ground state.

The basic idea of Iterative Projection method is to evolve the system along the imaginary time in a very small time step so that the time evolution operator can be expanded as a sequence of two-site unitary operations through the Suzuki-Trotter decompositions. For long enough time evolution, the system will go to the ground state we are solving for. We further assume translational invariance of the ground state ansatz (2), then we only need to update 4 tensors as shown in Fig. 1 inside a unit cell for each time step.

Based on the solved ground state from Iterative Projection method, we can calculate some quantities to characterize the phase diagrams. For the symmetry breaking phase we can evaluate the order parameters which are the vacuum expectation value (vev) of some physical operators. However, when we insert a physical operator at some particular site such that the translational invariance is lost, one needs some efficient method to contract the exponentially large number of bounds. Here we will adopt the TRG method (28); (29) to do that. The basic idea of TRG is to renormalize the tensors of the TPS by keeping the relevant entanglement when coarse graining. In this way, we can reduce the size of the system while retaining the essential quantum correlations of the original ground state. Finally the whole system will be reduced to a unit cell with the renormalized tensor, we can then evaluate the vev’s faithfully enough within the renormalized unit cell. Besides, there are alternative ways to implement polynomially efficient evaluation of the vev. One interesting method is to contract the tensors in the 1D scheme vertically and horizontally (30), and the other one is to use the Monte Carlo sampling to enhance the rate of contracting the tensors (31).

### ii.2 Order parameters for Néel and dimer phases

For the J-J’ model considered here, we will evaluate the Néel order parameter i.e., given by to characterize the Néel ordered phase. We also evaluate the spin-spin correlations and the dimerization (32); (33); (34):

(3) |

and

(4) |

to characterize the dimmer strength of the disordered “dimer phase”. Note that the Hamiltonian of the model does not possess the dimerized order by construction. Besides, we also use TRG to evaluate the ground-state energy per site as .

### ii.3 Characterizing topological phase by topological entanglement entropy

On the other hand, it is more difficult to characterize the topological phase since there is no symmetry breaking order parameter for it. On typical quantity for the 2D systems is the topological entanglement entropy (9); (10). It is the sub-leading constant term in the entanglement entropy

(5) |

where is the boundary size of the block for which the entanglement entropy is evaluated by tracing out the degrees of freedom outside it.

Since the constant term is topological and universal, we can extract it by appropriately subtracting the entanglement entropies of different blocks. There are two subtraction schemes, one is proposed by Kitaev and Preskill (9), and the other is by Levin and Wen (10). For the square lattice it is easier to implement the numerical evaluation of the entanglement entropy for a square or rectangular block. In this case, we find that it is convenient to adopt Kitaev-Preskill’s scheme as shown in the left of Fig. 2 for the square block with even number of sites on each side of the square block. The extraction for the topological entanglement entropy then goes as follows

(6) |

where denotes the von Neumann entropy of the density matrix in the region .

### ii.4 Characterizing topological phase by degeneracy of singular value spectrum

Another quantity in characterizing the topological phase is the degeneracy of the entanglement spectrum. It has been used to characterize the topological order for some quantum Hall states (11), and recently to characterize some symmetry-protected topological ordered Haldane phase in spin 1 chain (4). The entanglement spectrum is refereed to the spectrum of the Schmidt values for the bi-partition of the whole system. By definition, these Schmidt values are the square roots of the eigenvalues of the reduced density matrix of either of the two partitions. Basically, the degeneracy of the entanglement spectrum implies the spins at different sites are almost in the maximally entangled state, hence it implies that the ground state has topological order. In 1D, the entanglement spectrum is the same as the singular value spectrum of the matrix product state, and its power in characterizing the topological order is fully demonstrated (5); (7); (3); (4). However, for general 2D spin system there is no consensus on the power of entanglement spectrum in characterizing the topological order, see (8) for recent study on this issue ^{4}

Unlike the 1D case, the 2D entanglement spectrum is not the same as the singular value spectrum of the 1-site tensor in TPS though these twos should be related, especially for translationally invariant states. Numerical task in evaluating the 2D entanglement spectrum is almost the same as for the topological entanglement entropy, and thus difficult. Instead, it is far more easier and straightforward to evaluate singular value spectrum for TPS by merging the tensors of neighboring sites along x- or y-direction and then doing SVD. In this way, we can obtain four kinds of the singular value spectrum, which are associated with the four different links out of a site. For example, in the J-J’ model there are three -bonds and one -bond, and they may have different singular spectrum. This can be done by just using Iterative Projection method without further invoking TRG.

Besides, from the singular value spectrum one can straightforwardly evaluate the bipartite entanglement measure per length as well as the single-site von Neumann entropy (1-tangle). Both can be used to characterize the quantum critical point (35).

### ii.5 Quantum state renormalization

Since the topological order is closely related to the long range entanglement, which however, is usually contaminated by the short range entanglement. In order to make the characteristic of the topological order more explicit, one can implement the quantum state renormalization method (15); (37); (38); (39), which can remove the short range entanglement by performing local unitary transformation. The successive quantum state renormalization will then flow the original state to a simpler fixed-point state, which, however, has the same topological order of the original state. We can then classify the topologically ordered phases by the TPS’s tensors of the fixed-point states.

We now briefly describe the method as follows, and for more details see (39). We first form a positive double tensor by merging two layers of tensors and with the physical indices contracted. Note that the resultant double tensor is constructed so that it is invariant under the local unitary transformation, and then we can spectrally decompose it again back into two tensors with physical indices. The spectral weights encode the relevance of the entangled components. In this step we have removed the short range entanglement by the local unitary transformation taking the tensor into according to the significance of the spectral weights. Next, we will coarse grain the lattice labelled by the tensor by implementing one step of block decimation in the TRG method.

After repeating the above steps the original TPS will then flow into a fixed-point state with all the short range entanglement removed. If we perform quantum state renormalization while keeping the symmetry or gauge symmetry of the lattice, the final fixed-point tensors can then be used to classify the topological phases according to the symmetry structure.

## Iii Numerical results for the phase diagram of J-J’ model

In this section, we will present the numerical results for the phase diagrams of the J-J’ model of system size .

We first evaluate the vev of staggered magnetization in ground state as a function of by using the TRG method, and the result is shown in Fig. 2. In our calculation we consider the bond dimension up to and keep to ensure the accuracy of the TRG calculation. Here, is the cutoff on the bond dimension of the merging lattice during the coarse graining in the TRG method. We find that the numerical results converge quickly along range entanglementady for and the numerical program works efficiently. Therefore, we will take and keep in our numerical study in this paper.

From Fig. 3, it is observed that drops to zero both at and . This indicates two quantum critical points at the corresponding critical values of . It was known that a second-order transition from a Néel-ordered phase to a finite-gap disorder phase occurs at , which is obtained from the unbiased quantum Monte Carlo simulation (23); (24). Our result instead finds a critical point at . The discrepancy could be due to the not large enough used in our numerical calculation. Furthermore, we fit the critical exponent for the magnetization, i.e., . Our result is , which is close to the exponents of the 3D classical Heisenberg model, i.e., ^{5}

The other critical point at is frustrated since and , and there is possibility for the existence of topologically ordered phase due to the underlying degeneracy of the frustrated ground state. In Fig. 4 we show how the critical behavior in this regime changes when we change the bond dimension. Our result showing for is consistent with this expectation for a topologically ordered phase.

Since the system is frustrated in the regime of , our method based on TPS and TRG will be more reliable than the others. Despite that, we still compare with the results from the other methods. It was shown that classically there is a second order phase transition from the Néel phase to a helical phase at based on renormalized spin wave theory (RSWT) (20) or exact diagonalization (ED) (21). Furthermore, it was claimed in (21) that the critical point shifts to by using the coupled cluster method (CCM) to take the quantum fluctuations into account. However, we do not see such a shift in our results. Instead, we evaluate the ground state energy by TRG method and find the value is very close to the one obtained by the method of RSWT for a helical phase, especially in the Néel (classical) regime. The result is shown in Fig. 5. We see that for the agreement in energy between ours and RSTW’s starts to deviate, and our translationally invariant ground state is energetically more favored than the non-translationally invariant helical state. Similar deviation occurs for , this again reflects the relevance of quantum effect when the system is away from the classical Néel ordered regime.

For curiosity, we also calculate the dimerization and defined in (3) and (4) though they are not the order parameters of our dimerized model. The results is shown in Fig. 6. We see that at as expected, otherwise they are nonzero and reach some constant values at the disordered (dimer) regime. To visualize the distribution of the dimer strength, in Fig. 7 we explicitly show the values of the spin-spin correlation of the neighboring sites for three different values.

Beside using the order parameter such as in this case to characterize the quantum phase transitions, the entanglement measure can also do the job, see (35) for some examples. Especially, for the TPS ansatz it is quite straightforward to evaluate the bipartite entanglement per bond , where labels the two different types of the bonds, and the ’s are the singular values obtained from the SVD of the site tensors in performing the Iterative Projection algorithm. Similarly, one can also use TRG to obtain the one-site reduced density matrix, and then evaluate the von Neumann entropy, i.e., 1-tangle denoted as . The results are summarized in Fig. 8. Note that the derivatives of both and are discontinuous at two quantum critical points. Moreover, it shows that the classical Néel regime has lower entanglement compared to the other two phases dictated by the significant quantum effect. For the disordered dimer phase with quite amount of the quantum entanglement is due to the formation of the dimers, i.e., nearest-neighboring sites form the maximally entangled states (Bell states). Once the dimer is formed, the sites on the both ends of the dimer will not be correlated with the other sites, i.e., the monogamy of the entanglement.

The local nature of short range entanglement for the dimer phase is different from the long range entanglement for the topological phase in the regime . Later, we will see that the difference between these twos will be reflected in the degeneracy of the singular value spectrum.

## Iv Topological Entanglement Entropy

In this section, we will evaluate the topological entanglement entropy based on TPS ansatz. We will consider first the J-J’ model, and then a toric code like toy TPS state to demonstrate the power of TRG method in identifying the topological phase.

### iv.1 J-J’ model

In order to make sure that the phase in the regime is the topological phase, we evaluate the topological entanglement entropy (i.e., ) according to the prescriptions mentioned before. We also find that the entanglement entropy increases linearly with the block boundary size, this gives a consistent check for our numerical codes.

Our result is shown in Fig. 9. The one for 2 by 2 block is evaluated by Kitaev-Preskill’s scheme, and the one for 3 by 3 block by Levin-Wen’s scheme. For larger blocks, it requires far more power of computation in TRG for evaluating the block entanglement entropy. For example, if we would like to evaluate the entanglement entropy for the 4 by 4 block, then the number of entries for the reduced density matrix will be about times of the one for 3 by 3 block. That means it will take about 100 times of the computing time than for the 3 by 3 block, and it is beyond our computational capacity. This is the disadvantage when trying to identify the topological phase with the topological entanglement entropy numerically. Despite that, we still spent months to evaluate the topological entanglement entropy for the 4 by 4 block at as shown in Fig. 9. The values of these two data points are quite consistent with the ones for the 3 by 3 block. More closely, for the 3 by 3 block at and at , and for the 4 by 4 block at and at , respectively.

Our result shows that is consistent with zero for the non-topological phases (Néel and dimer phases). However, we find that is negative and unphysical in the topological phase for 2 by 2 block. This could be due to the small boundary size effect from the next sub-leading term of with . On the other hand, for the 3 by 3 block we expect the small boundary size effect could be suppressed so that in the topological phase should be positive. We find that this is indeed the case, and this is also confirmed by the two data points for the 4 by 4 block. In Fig. 9, there is a long crossover for as we vary between and . Then, approaches around and then remains at that value for more negative as expected from its topological nature. Naively, we expect a quantum critical point separating topological and non-topological phase, and the crossover could be due to the short range entanglement in our TPS solution. We expect that it should be removed by suitable quantum state renormalization.

Though our result is subjected to the errors from the small size effect and also from the approximation in our numerical method of Iterative Projection and TRG, it does indicate a nonzero and thus a topological phase in the regime with . We may expect the value of will be improved if the above errors could be reduced.

### iv.2 Toric code like state

Due to the limitation of our computation power, it is hard to evaluate the topological entanglement entropy to high precision by increasing the bond dimension in the TPS ansatz for the J-J’ model. Instead, we now consider a one-parameter toric code like TPS state proposed in (39) with . In this case we can evaluate the topological entanglement entropy more precisely to characterize the topological order, and demonstrate the power of TRG method on this issue. Later on, we will compare the topological entanglement entropy of this state with the pattern of degenerate singular value spectrum of the same state, and demonstrate that the latter can be also used to characterize the topological order without invoking heavy computation capacity. Though this one-parameter state has no relation to the J-J’ model, however, it yields the supporting evidence for the connection between topological order and degenerate singular values in the topological phase of the J-J’ model, and possibly for more other frustrated spin systems.

This toy TPS state on the square lattice motivated by the toric code model is characterized by the TPS tensor (39) on each vertex in Fig. 10 in the form of with four physical indices and four bond indices (i.e., bond dimension ). Its entries are given by

(8) |

The parameter is used to tune the property of this state from the topological phase to the trivial phase. Note that this state is not solved from any model Hamiltonian but is devised on purpose to demonstrate the phase transition between topological and trivial phases. For , it reduces to the ground state of the toric code model with topological order. For , it reduces to a product state of all zeros, i.e., a zero state. Therefore, as we vary , the state will go through a phase transition, and in (39) it showed that the quantum critical point occurs at which separates the topological phase from the phase of product state.

With the above TPS state, we then use the Kitaev-Preskill’s scheme to numerically evaluate the topological entanglement entropy for 2 by 2 block indicated by the blue dash line in Fig. 10. Since this state has only and thus requires not much computation power to carry out the numerical calculation up to high precision. Again, we verify that the entanglement entropy obeys the area law. Our result is shown in Fig. 11, and it indicates a crossover around separating the trivial phase with and the topological phase with . We may observe a sharp jump at if we increase the block size, however, it is then beyond our computation capacity.

## V Degeneracy of Singular Value Spectrum and Topological Order

In general, there are some ways of classifying the possible topological phases, however, there is no specific method with consensus in identifying the topological phase of a particular dynamical model. Topological entanglement entropy (9); (10) is usually thought to be good index in characterizing the topological order, but it is difficult to evaluate numerically to high precision due to the limitation of the computation capacity as discussed in the previous section.

On the other hand, the degeneracy of the entanglement spectrum has recently implemented to characterize the topological orders for some 2D quantum Hall states (11) and some 1D symmetry protected topological phase (4). For 2D TPS, it is more easier to consider the singular value spectrum, which is related to the entanglement spectrum in some way. Despite there is no rigorous proof on its power of characterizing the 2D topological orders, for the translationally invariant states the degeneracy of the singular value spectrum by itself indicates the existence of the long range entanglement. Recently it is argued that the long range entanglement is the characteristics of the topological order from the point of view of the quantum information (5); (6); (7); (4). The tensor of the TPS encodes the information about the entanglement for the bi-partitions of the system. The singular value spectrum obtained from it weighs different components of quantum correlation between bi-partition regions in the ground state wave function. Therefore, if there is a double degeneracy of the singular value spectrum, with the help of the translational invariance it may imply that the bi-partitions are almost in the maximally entangled states. It is the long range entanglement since it involves the the collective coherences of the spins on half of the space to be strongly correlated with all the spins in the other half.

By considering the frustrated regime of the J-J’ model, we demonstrate the above picture is on the right track to uncover the topological nature of the long range entanglement. It needs more study on both numerical and analytical sides to solidify the above connection. From our on-going work (36), we find more positive evidences.

### v.1 Some examples

Instead of providing the proof of the above picture, we give examples. A well-studied one is the symmetry-protected Haldane phase, which is described by the Afflect-Kennedy-Lieb-Tasaki (AKLT) state. The AKLT state can be put into the form of the MPS represented by the matrices . By merging two matrices and then performing SVD, we find that the singular values are doubly degenerate. This is also proved in more rigorous way in (4). Since the topological phase is robust again small perturbations, we consider the spin-1 model Hamiltonian

(9) |

At the large , there is a trivial phase where all the spins are in the eigenstates of . On the other hand, for small the system is in the gapped Haldane phase, which is related to the AKLT state by local unitary transformation. We use the iTEBD method to solve the ground states in the MPS ansatz and then find the singular value spectrum. We find that the phase transition between Haldane phase and the trivial phase is characterized by the degeneracy of the singular value (entanglement) spectrum.

The relation between the degeneracy of the singular value spectrum and the topological order for the 2D spin system is less explored. Here we consider the simplest toric code model (16), and the explicit TPS form of its ground state is represented by the tensor as discussed in the previous section. We can find the singular value spectrum by merging the tensors and then performing SVD. Again, we find that four kinds of the singular value spectrum are all two-fold degenerate.

### v.2 Topological order of J-J’ model by singular value spectrum and its robustness

The singular value spectrum for each bond in the J-J’ model from the TPS of our numerical ground state solution is given in Fig. 12(a). Each small circle with a middle bar represents a specific singular value of the singular value spectrum. Note that the singular value spectrum for the - and -bond may show different degeneracy patterns, however, if the patterns are similar we will only show one of them. We see that the singular values are doubly degenerate for the regime in which the Néel order also vanishes. Therefore, should be the critical point separating the topological and non-topological phases, along with the evidence of nonzero topological entanglement entropy discussed previously. On the other hand, for the dimmer phase in the regime with , we find the singular values for the -bond are doubly degenerate, but others aren’t. This reflects the peculiar feature of the strong dimmer bond but not the topological order because the topologically entanglement entropy is zero in this phase. Moreover, our results show that not just the dominant singular value but almost all the singular values for that bond are also doubly degenerate.

Recall that the precision for the topological entanglement entropy and the sharp quantum transition point is hard to achieve due to the limitation on the computation capacity. However, with the same capacity we can obtain the degenerate singular value spectrum quite easily with a sharp quantum transition point. In contrast with the topological entanglement entropy, this is the advantage to use the degenerate singular value spectrum to characterize the topological order numerically with TPS ansatz.

Since the topological order should be robust against the local perturbations, the degeneracy of the singular value spectrum should also does so if it can be used to characterize the topological orders. We first add a random perturbation at each step of Iterative Projection method, where the random 2-body Hamiltonian is a real symmetric matrix with 10 independent matrix elements given by the random numbers in the interval . The result is shown in Fig. 12(b). We see that the singular value spectrum is changed but the pattern of its degeneracy remains intact. This demonstrates the robustness of the degeneracy of the singular value spectrum, and thus the robustness of the topological order.

To further check the robustness, we add a perturbation of the form which breaks the spin SO(3) symmetry and the time-reversal symmetry. For small , the above perturbation does not destroy the double degeneracy of the singular value spectrum. Once is larger than the critical value , the double degeneracy will be lost. The results for and are shown in Fig. 13(a) and (b), respectively. On the other hand, in Fig. 13(c) we show the result by adding a perturbation of the form which breaks the translational symmetry in the direction perpendicular to the dimers. Again, we see that the perturbation with cannot destroy the double degeneracy.

With above generic types of perturbations, we see that the degeneracy pattern of the singular value spectrum is quite robust, this reflects the topological nature of the phase.

### v.3 Remove short range entanglement

As described previously, quantum state renormalization consists of two steps: the first step is to optimize the local unitary transformation to remove the short range entanglement, and the second step is to coarse grain the site tensors according to the TRG like method while keeping the physical bonds. After successive steps of quantum state renormalization, the TPS will flow to a fixed point with most of short range entanglement being removed. Therefore, one would expect the fixed-point TPS will be the representative of the universal class of topologically ordered states. This quantum informational way of classifying the topological orders are recently proposed in (5); (6); (7); (4); (39).

We then would like to study how the singular value spectrum evolves under the quantum state renormalization group (RG) flow procedure. Since we conjecture that the degeneracy pattern of the singular value spectrum encodes the long range entanglement, we would expect it is robust under the quantum state renormalization, which only removes the short range entanglement. This is indeed the case for the J-J’ model and the result is shown in Fig. 14. Moreover, we see that under the quantum state RG flow, the number of the dominant singular values decreases in the topological phase. We can think of the extreme case that there is only one dominant doubly degenerate singular value with all others being suppressed, then the whole lattice is in the GHZ state, i.e., maximal long range entanglement. So, the suppression of the singular values except the dominant ones corresponds to the removal of the short range entanglement under the quantum state RG flow, and the long range entanglement will be encoded in the degeneracy of the dominant singular values of the fixed-point TPS.

Finally, to get more supporting evidences about the topological nature of the singular value spectrum under quantum state RG flow, we consider the case for the toric code like TPS state (8) discussed before. To draw a parallel comparison with the topological entanglement entropy presented in Fig. 11, we evaluate the difference between the two largest singular values, denoted as . The result is shown in Fig. 15.

From Fig. 15, we see that before quantum state renormalization, it seems there is no degenerate singular value spectrum except at . Moreover, the degeneracy at is 2-fold, which may reflect the underlying symmetry apart from the long range entanglement. This is in sharp contrast to the crossover behavior of topological entanglement entropy around as shown in Fig. 11. However, after one step of quantum state renormalization, we obtain the crossover similar to the one for topological entanglement entropy. As we perform more steps of quantum state renormalization, the crossover become sharper and sharper around . It finally approaches to a sharp quantum critical point separating the topological phase () from the non-topological one (). Similar behavior for some scale invariant quantity of the same state (8) under quantum state RG flow is also found in (39).

Moreover, we can get some lesson as follows. Though the degeneracy of the singular value spectrum may indicate the existence of the long range entanglement, the reverse statement may not be true due to the contamination of the short range entanglement. However, the hidden long range entanglement and the underlying symmetry such as for the toric code ground state will emerge under the quantum state RG flow as shown in Fig. 15. This is consistent with the robustness check for the topological order of toric code done in (40).

## Vi Conclusion

In this paper, we identify a new topological phase for the J-J’ model by explicitly calculating the phase diagram of order parameters, topological entanglement entropy and the degeneracy of the singular value spectrum. This result demonstrates that the recent quantum information inspired methods such as TPS, Iterative Projection method and TRG are powerful enough in tackling the frustrated spin systems and identifying the topological orders. Especially, we demonstrate that the concepts of short range entanglement and long range entanglement introduced in (5); (6) are very useful in providing intuitive picture for the entangled nature of topological order. It is also helpful in interpreting the numerical results based on TPS ansatz.

We have also found an intriguing connection between the degeneracy of the singular value spectrum and topological order by explicitly solving the J-J’s model as a nontrivial example. The topological order may not always reflect in the degeneracy of the singular value spectrum due to the disturbing short range entanglement. However, as we show explicitly, we can remove short range entanglement by quantum state renormalization so that the degeneracy of the singular value spectrum after quantum state renormalization indicates the topological order. The implication of such a connection is far-reaching in numerical identification of the topological phase in 2D systems since the computation power of obtaining the singular value spectrum is far less than the one of evaluating the topological entanglement entropy.

Based on our observation, we believe that we can use the degeneracy of the singular value spectrum based on TPS ansatz to identify the topological phases of other interesting frustrated 2D systems. We will report more works along this direction in the near future.

## Acknowledgements

We thank H.Y. Chen, Z. C. Gu and A. Sandvik for helpful discussions. We also thank Yu-Cheng Lin for the generous support on computer facility. This work is also supported by Taiwan’s NSC grant 97-2112-M-003-003-MY3. We also thank the support of NCTS.

### Footnotes

- 896410093@ntnu.edu.tw
- linfengli@phy.ntnu.edu.tw
- preprint:
- In (8), the entanglement spectrum of a bulk region is associated with a the boundary excitation spectrum. Their numerical study revealed that the boundary Hamiltonians become non-local for the topologically ordered states.
- By increasing one may improve the above the results to be close to the results by the quantum Monte Carlo method, though it may require far large computational power. In this work, we are more interested in the identification of the topological phase, and will be satisfied by the above accuracy for the phase transition.

### References

- T. Senthil, A. Vishwanath, L. Balents, S. Sachdev and M. P. A. Fisher, Science 303, 1490 (2004).
- T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 032110 (2002).
- F. Pollmann, E. Berg, A. M. Turner, and M. Oshikawa, arXiv:0909.4059.
- F. Pollmann, E. Berg, A. M. Turner, and M. Oshikawa, Phys. Rev. B 81, 064439 (2010).
- Z.-C. Gu and X.-G. Wen, Phys. Rev. B 80, 155131 (2009); X. Chen, Z.-C. Gu, and X.-G. Wen, arXiv:1008.3745 and arXiv:1103.3323.
- X. Chen, B. Zeng, Z.-C. Gu, I. L. Chuang, and X.-G. Wen, Phys. Rev. B 82, 165119 (2010); Z.-C. Gu, Z. Wang, and X.-G. Wen, arXiv:1010.1517.
- N. Schuch, D. Pérez-García, and I. Cirac, arXiv:1010.3732.
- J. I. Cirac, D. Poilblanc, N. Schuch, F. Verstraete, arXiv:1103.3427.
- A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006).
- M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 (2006).
- H. Li and F. D. M. Haldane, Phys. Rev. Lett. 101, 010504 (2008).
- F. Verstraete, J. I. Cirac, arXiv:cond-mat/0407066; V. Murg, F. Verstraete, and J. I. Cirac, Phys. Rev. A. 75, 033605(2007).
- X. G. Wen and Q. Niu. Phys. Rev. B. 41, 9377 (1990).
- B. A. Friedman and G. C. Levine. Phys. Rev. B 78,035320 (2008).
- F. Verstraete, J. I. Cirac, J. I. Latorre, E. Rico, and M. M. Wolf, Phys. Rev. Lett. 94, 140601 (2005).
- A. Yu. Kitaev, Ann. Phys. 303, 2 (2003).
- M. Levin, X.-G. Wen, Phys. Rev. B 71, 045110 (2005).
- A. Kitaev, Ann. Phys. 321 , 2 (2006).
- R. R. P. Singh, M. P. Gelfand, and D. A. Huse, Phys. Rev. Lett. 61, 2484 (1988).
- N.B. Ivanov, S.E. Krüger, and J. Richter, Phys. Rev. B 53, 2633 (1996).
- S.E. Krüger, J. Richter, D. J. J. Farnell and R. F. Bishop, Phys. Rev. B 61, 14607 (2000).
- B. Bauer, G. Vidal, M. Troyer, J. Stat. Mech. P09006 (2009).
- S. Wenzel, L. Bogacz, and W. Janke1, Phys. Rev. Lett. 101, 127202 (2008).
- F.-J. Jiang, arXiv:0911.0653.
- H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008).
- P. Corboz, R. Orus, B. Bauer, G. Vidal, Phys. Rev. B 81, 165104 (2010).
- G. Vidal, Phys. Rev. Lett. 98, 070201(2007); R. Orús and G. Vidal, Phys. Rev. B. 78, 155117 (2008).
- M. Levin and C. P. Nave, Phys. Rev. Lett. 99, 120601 (2007).
- Z. C. Gu, M. Levin and X. G. Wen, Phys. Rev. B. 78, 205116 (2008).
- J. Jordan, R. Orus, G. Vidal, F. Verstraete, J. I. Cirac, Phys. Rev. Lett. 101, 250602 (2008)
- L. Wang, I. Pizorn, F. Verstraete, Phys. Rev. B 83, 134421 (2011).
- P. W. Leung and N.-W. Lam, Phys. Rev. B 53, 2213 (1996).
- R. R. P. Singh, Z. Weihong, C. J. Hamer,and J. Oitmaa, Phys. Rev. B 60, 7278 (1999).
- Y. Yu and G. Müller, V. S. Viswanath, Phys. Rev. B 54, 9242 (1996).
- C.-Y. Huang, F.-L. Lin, Phys. Rev. A 81, 032304 (2010).
- C.-Y. Huang and F.-L. Lin, in preparation.
- G. Vidal, Phys. Rev. Lett. 99, 220405 (2007).
- R. Orús, Phys. Rev. Lett. 100, 130502 (2008).
- X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B. 82, 155138 (2010).
- S. Dusuel, M. Kamfor, R. Orus, K. P. Schmidt, J. Vidal, Phys. Rev. Lett. 106, 107203 (2011).