Abstract
The lowtemperature dynamics of quantum systems are dominated by the lowenergy eigenstates. For twodimensional systems in particular, exotic phenomena such as topological order and anyon excitations can emerge. While a complete lowenergy description of strongly correlated systems is hard to obtain, essential information about the elementary excitations is encoded in the eigenvalue structure of the quantum transfer matrix. Here we study the transfer matrix of topological quantum systems using the tensor network formalism and demonstrate that topological quantum order requires a particular type of ‘symmetry breaking’ for the fixed point subspace. We also relate physical anyon excitations to domainwall excitations at the level of the transfer matrix. This formalism enables us to determine the structure of the topological sectors in twodimensional gapped phases very efficiently, therefore opening novel avenues for studying fundamental questions related to anyon condensation and confinement.
Introduction
One of the major recent advances in the understanding of strongly correlated quantum many body systems has been the investigation of quantum entanglement in terms of area laws^{1,2}, the entanglement spectrum^{3} and the associated entanglement Hamiltonian^{4,5}. The structure of entanglement in gapped quantum systems has resulted in the powerful parameterization of quantum ground states in terms of socalled tensor network states, such as matrix product states (MPS)^{6,7} or their higher dimensional analogues, projected entangled pair states (PEPS)^{8,9}. In a translation invariant tensor network state, the entanglement features can be extracted from the leading eigenvector of the socalled (quantum) transfer matrix, which naturally appears whenever the quantum state is mapped to a classical partition function (using for example, a Trotter decomposition) and then sliced along the virtual direction, that is, the Trotter direction or direction of imaginary time^{10}. In a very recent publication^{11}, we have observed that the other eigenvalues of the transfer matrix also contain useful information, thst is, they contain crucial information about the elementary excitations and the corresponding dispersion relations of the system. This is, in many ways, surprising as this information is completely encoded in the groundstate description of the system, a priori without any information about the original Hamiltonian for which it was the exact or approximate ground state.
This relation is extremely useful for the case of twodimensional systems, where systematic methods to extract the dispersion relation of elementary excitations are virtually nonexistent. This was illustrated in ref. 11 by studying the Affleck–Kennedy–Lieb–Tasaki model^{12,13}. However, for the case of twodimensional quantum spin systems, exhibiting topological quantum order, the excitation spectrum can be much richer. In contrast to trivial phases, where the elementary excitations can provably be created as Bloch waves of localized perturbations^{14}, elementary excitations in topological phases are typically anyons that come with ‘strings attached’. Note that topologically nontrivial excitations also appear in onedimensional systems, such as in the Lieb–Liniger model^{15} or as a domainwall excitations in systems with discrete symmetry breaking. In terms of MPSs, those excitations can be represented with a Bloch like ansatz where an extra halfinfinite string operator is attached to the local perturbation^{16,17,18,19}.
The PEPS representation of topologically ordered ground states, either as variational approximation^{20,21} or as exact description of the Levin–Wen model wave functions^{22,23,24}, has been well established. Here we study the full spectrum of the transfer matrix of topologically ordered PEPS^{5} on the infinite plane or cylinder. Much like the prisoners in Plato’s cave, we observe onedimensional domain walls in the spectrum of the transfer matrix as shadows of the true anyons in the twodimensional world. By clarifying how the different anyon sectors are manifested at the virtual level, we can probe the dispersion relation of single anyon states. We discuss how the presence of anyons, and thus of topological order, requires a particular type of symmetry breaking in the fixed point subspace of the transfer matrix and how anyon condensation or confinement^{25,26} is reflected in these virtual description. We illustrate our results by studying the PEPS description of the toric code model with string tension and the resonating valence bond state. Our results also confirm that we can construct approximate eigenvectors of PEPS transfer matrices using the matrix product ansatz developed for onedimensional quantum Hamiltonians in refs 17, 27.
Results
Topological order in PEPS
The convience of the tensor network description of quantum states is that the global, topological properties of the state are reflected in the symmetries of the local tensors. Since topological phases are not characterized by local order parameters, these symmetries act purely on the virtual levels of the tensors. In particular, it was recently established that topological order in PEPS can be characterized by the existence of matrix product operators (MPO), which can be pulled through the lattice at the virtual level (see Fig. 1)^{28,29,30}. Closed MPO loops around a topologically trivial region define the invariant subspace on which the PEPS tensors are supported and in this way characterizes the topological properties of the state, such as the topological corrections to the entanglement entropy. They act as virtual operators O_{i} and satisfy a fusion algebra . Indeed, as shown in ref. 30, for the case of the Levin–Wen string net models, the different MPOs O_{i} can be associated to and labelled by the different string types i=1,…,N of the input category that defines the string net model. A PEPS in a trivial phase is characterized by a single MPO O_{1} that acts as the identity in the relevant subspace. Another interesting case is that of the quantum double models, which can be described using Ginjective PEPS^{28}. This is a special case of the formalism of ref. 30 where the MPOs are labelled by the group elements g∈G and correspond to representations O_{g}=U(g) of the group action at the virtual level. The pullingthrough condition is satisfied since the tensors are only supported on the invariant subspace defined by the projector P=∑_{g}O_{g}.
While the pullingthrough conditions ensure that the presence of an MPO string cannot be detected locally, noncontractible MPO loops can have global effects, such as adding a nontrivial flux in the system, and can, therefore, be used to map one ground state to another one. The relevance of these virtual MPOs is that also away from the the renormalization group fixed point—where the physical string operators are spread out and not exactly known^{31}—the MPOs at the virtual level remain strictly local and the ‘pullingthrough’ symmetry of the PEPS tensor is exactly preserved.
Anyon excitations in the PEPS picture
Having a translation invariant PEPS description of the ground state of a topological phase, one can easily argue that a suitable ansatz to model single anyon excitations is obtained by modifying the groundstate tensors in a local region (for example, a single site) and attaching a halfinfinite string to it, which is exactly given by this MPO at the virtual level. The MPO will give rise to the nontrivial braiding statistics of these excitations, while the ‘pullingthrough’ assures that the bulk of the string is locally unobservable, so that the energy density is left at its groundstate value sufficiently far away from the end point. Away from the renormalization group fixed point, these excitations will disperse and a proper eigenstate can be obtained by building a momentum superposition with the momentum k_{x} and k_{y} in the x and y direction.
Note that the topological quantum numbers of the anyon excitations are not completely specified by the string type i, but are determined by the structure of the excitation tensor in the ansatz of Fig. 2. For the case of the quantum doubles, it has been shown that the string type corresponds to the magnetic flux, whereas the charge quantum number is determined by the representation space on which the local tensor is supported^{28}. A complete characterization of the different anyon sectors in the PEPS formalism would take us to far and is presented elsewhere (N. Bultinck, M. Mariën, D. Williamson, J. Haegeman, F. Verstraete, manuscript in preparation).
Transfer matrix: symmetries and domain walls
As in ref. 11, one can now argue that the dominant contribution to the variational dispersion relation is coming from the normalization of these states, which is given by the sum of overlaps of ket and bra with string end points at different positions (x_{0}, y_{0}) and (x_{0}+Δx, y_{0}+Δy), as illustrated in Fig. 2. If we orient the strings along the x direction and first contract the tensor network along the y direction, the central object will be the transfer matrix in the x direction, as defined in Fig. 3. The pullingthrough condition of Fig. 1 ensures that commutes with infinite MPO strings along the x axis in the ket and bra level separately. We thus obtain , ∀i=1,…, N where O_{i} now denotes an infinite MPO of type i along the x axis. Normalization of the PEPS ground state requires that the largest eigenvalue of is 1, and the infinite power of the transfer matrix in the overlap of Fig. 2 in the regions y<y_{0} and y>y_{0}+Δy can be replaced by its left and right fixed point σ and ρ, which we represent as an infinite MPS with matrices A_{σ} and A_{ρ} (Fig. 3).
The transfer matrix can have a degenerate fixed point structure, since for a given right fixed point ρ, one can build other fixed points for all i, j=1,…, N. One could expect that this generally gives rise to an N^{2}dimensional fixed point subspace. However, at the renormalization group fixed point of the topological phase, we can easily check that the fixed point subspace of is exactly spanned by ρ_{k}=O_{k} for k=0,…, N, and is thus only Ndimensional. The degeneracy and labelling of the fixed point subspace remains intact throughout the topological phase, even though ρ_{k} will no longer exactly equal O_{k}. This implies, in particular, that can be expanded into a linear combination of ∑_{l}c_{kl}ρ_{l} with c_{k,k}=1, compatible with the fact that fusing will have a fusion channel k. We now argue why this property is required to support anyonic excitations with a halfinfinite virtual string of type i.
Contracting the tensor network in Fig. 2 from right to left up to position y_{0}+Δy gives rise to some boundary state ρ in the fixedpoint subspace, whose precise choice is set by the boundary conditions at y=+∞. The topological invariance ensures that this choice has no effect on local expectation values. As we now further contract from right to left, we pass the position y_{0}+Δy containing the excitation in the bra level. Here the boundary state is perturbed locally at x=x_{0}+Δx. In addition, it will be acted on by a halfinfinite MPO string of type i, which has the effect of changing the MPS tensors from A_{ρ} to A_{ρ′} with x<x_{0}+Δx. At the level of the transfer matrix, the boundary state now takes the form of a domain wall interpolating between the two different fixed points. Because of translation invariance, all overlaps corresponding to the momentum superposition in the x direction can be summed and the resulting state takes the form of a topologically nontrivial state with momentum k_{x}, similar to the domainwall excitation ansatz used for onedimensional Hamiltonians in ref. 17. We thus have to consider the spectrum of topologically nontrivial eigenstates ξ_{j} (Fig. 3) of with momentum k_{x}. Indeed, further contracting up to y_{0} yields the k_{x}dependent eigenvalues λ_{j}(k_{x}) of ξ_{j} to some power Δy, which dictates the k_{y} dependence of the dispersion relation of the corresponding physical excitation.
At point y_{0}, the boundary state is acted on with a second halfinfinite string (also with momentum k_{x} in the x direction), now in the ket level. After that, the state is collapsed onto the topologically trivial left fixed point σ. To have a nonvanishing overlap, the corresponding boundary state—with halfinfinite strings in both the ket and bra level—should have a contribution in the trivial sector. This will be true if the property stated above is satisfied, that is, for all ρ.
As we perturb the state out of the topological phase, there are two ways in which the existence of topologically nontrivial excitations can break down. If the transfer matrix has a unique (maximally ‘symmetric’) fixed point, which is invariant under the action of any string O_{i}, then the states ξ interpolate between the same fixed point left and right and are thus indistinguishable from topologically trivial local perturbations (without string). This scenario is realized when the corresponding anyon has condensed into the ground state^{25,26}. A second possibility is that the fixedpoint structure of has an even larger ‘symmetry breaking’, in such a way that the state with halfinfinite strings in ket and bra level is still completely topologically nontrivial and has zero overlap with the left fixed point σ. More specifically, the fidelity per site, defined as the largest eigenvalue λ of with , satisfies λ=e^{−t}<1 and the normalization of the state goes down as e^{−tL} with L the length of the strings. Hence, the nonzero value of t acts as a string tension and only bound states of two excitations connected by a finite string can exist, corresponding to the mechanism of confinement.
In the generic case, there can of course be several transitions corresponding to, for example, the condensation of only some anyonic sectors and the induced confinement of other anyon sectors. Note that, since topological phase transitions correspond to symmetry breaking phase transitions at the virtual level of the PEPS description, we can also find virtual order parameters, as illustrated in the toric code example. A more indepth study of these aspects of anyon condensation within the framework of PEPS will be provided elsewhere (J. Haegeman, N. Schuch, F. Verstraete, manuscript in preparation).
Mixed transfer matrix and momentum fractionalization
In the above discussion, we have explained how information about anyon excitations in topological phases can be obtained from the topologically nontrivial excitations of the translation invariant transfer matrix, which has a degenerate fixed point subspace in the case of topological order. Fig. 4 motivates an alternative approach. By using the pullingthrough property of the MPO, we can rewrite the eigenvalue equation for a topologically nontrivial excitation of as a normal (topologically trivial) eigenvalue problem for a socalled mixed transfer matrix. The latter is threaded by an MPO string and is thus defined on a larger vector space corresponding to the presence of additional MPO indices. Physically, we are effectively rotating the MPO strings attached to the anyon excitation to lie along the y direction.
By doing so, we can make the x direction finite and periodic, which allows to work on a cylinder with finite circumference. The fixed points of these mixed transfer matrices were first studied in ref. 5. In this case, physical translations in the x direction have a representation as modified translation operators at the virtual level of the transfer matrix with a nontrivial action on the extra MPO indices (see Fig. 5). This results in a momentum label that can have a fractional discretization in the circumference of the cylinder, similar to what is observed in the case of momentum polarization^{32}. We elaborate on this aspect in Supplementary Note 1 and Supplementary Figs 1 and 2.
Toric code model with string tension
Let us now illustrate this approach using the toric code ground state, to which we apply a local filtering
with Ω〉 the fully polarized spin state . The filtering induces dynamics to the elementary excitations of the toric code ground state and can drive the system into a trivial phase. Along the coordination axes β_{x}=0 or β_{z}=0, it can be interpreted as string tension in either the group or representation basis and can be solved exactly^{33,34}. The full twodimensional phase diagram as a function of β_{x} and β_{z} was studied in ref. 35 using a fidelity approach^{36,37,38}.
The PEPS representation of the toric code ground state was first constructed in ref. 39 and its properties were discussed at length in ref. 28. The PEPS tensors satisfy the property of Ginjectivity, where in this case . For this particular case, it means that the MPO projector P is given as O_{0}+O_{1}=^{⊗L}+Z^{⊗L}, with L the length of the MPO string and Z a representation of the nontrivial element of . Correspondingly, the transfer matrix along the x direction has the global symmetry with N_{x} the number of sites in the x direction. The filtering operation is applied at the physical level of the PEPS and has no effect on any of these properties. It does, however, influence the manifestation of the symmetry in the fixed point subspace of the transfer matrix.
At the virtual level of the PEPS, we can use local order parameters X⊗ and X⊗X to detect the symmetry breaking of Z⊗Z and Z⊗, respectively, in the fixedpoint subspace, where X is an operator such that XZ=−ZX. Fig. 6 shows the structure of the fixed point subspace. The topological ‘toric code’ phase is characterized by a doubly degenerate fixed point, which is invariant under the subgroup but breaks the symmetry under the action of . Creating a physical excitation with flux quantum number 0 (string of identities) or 1 (string of Z's) is manifested at the virtual level as a boundary state in the topologically trivial or, respectively, topologically nontrivial sector. That is, excitations with nontrivial flux correspond to domainwall excitations at the level of the transfer matrix. The charge of the physical excitation can be measured as a charge difference between the ket and bra level of the boundary state using the operator (which is a symmetry, that is, eigenvalue 1, for the fixed point subspace).
Following ref. 11, we can now probe the dispersion relation of the elementary excitations of the model by inspecting the spectrum of (minus logarithm of) the eigenvalues of the transfer matrix in the different sectors, for which we use the onedimensional excitation ansatz of ref. 17. The result is illustrated for various values of β_{z} and β_{x} inside and outside the topological phase in Fig. 6. In the topological phase (plots (a), (b), (e) and (f)), the eigenvectors can be labelled by the charge difference between ket and bra (corresponding to the physical charge) and the absence or presence of a halfinfinite string (corresponding to the physical flux). We can then relate the spectrum of the transfer matrix to the dispersion relation within the four physical topological sectors and recognize the charge and flux as elementary excitations (isolated branch). In the first column of Fig. 7, the gap in the charge sector closes, resulting in charge condensation and a phase transition to the trivial phase (plots (c) and (d)). Here the fixed point subspace breaks the full symmetry of the transfer matrix and the charge differences between ket and bra is no longer well defined. Correspondingly, we can create new topologically nontrivial excitations with a halfinfinite string of Z's in both ket and bra. That this sector has a gap Δ_{ZZ} indicates that physical flux excitations can no longer exist in isolation and must be confined to pairs, since their normalization goes down as exp(−Δ_{ZZ}L) with L the length of the string. In the right column of Fig. 7, the gap in the flux sector closes, corresponding to physical flux condensation. This triggers a phase transition to the trivial phase and results at the virtual level in a unique fixed point with the full symmetry of the transfer matrix. Correspondingly, there are no more topologically nontrivial excitations and we can now measure individual charge numbers of the ket and the bra level. That there is a gap in the (−,−) sector with negative charge in ket and bra indicates that physical charge excitations can no longer exist in isolation and must be confined.
Resonating valence bond state on the hexagonal lattice
Finally, in Fig. 8 we present the spectrum of the transfer matrix for the resonating valence bond state^{40} on the Kagome lattice, for which the PEPS is also injective^{41,42}. The Kagome lattice was blocked as illustrated in Fig. 9b, and the eigenvalues of the transfer matrix along the lattice vector were computed, giving full access to the momentum k_{1}. The phase of these eigenvalues can then be interpreted as momentum k_{2}, which allows to map them to the Brillouin zone according to Fig. 9c. This clearly allows to extract the physical elementary excitations^{43,44}, namely the spinons (S=1/2, no string) and visons (S=0, string), but also vison–spinon bound states (S=1/2, string), which seem to occur at energies lower than the vison energies. Finally, there are also an isolated branch of trivial spinless excitations (S=0, no string), which could be a bound state of either a vison pair or a spinon pair. For this model, the transfer matrix is a nonhermitian MPO and one can question the validity of an excitation ansatz based on a local perturbation on top of the fixed point (which is approximated as MPS), as this ansatz is only provably justified for the case of local Hamiltonians^{14}. Numerical evidence for the validity of our results is provided in Supplementary Note 2 and Supplementary Figs 2 and 3.
Discussion
We have illustrated how the eigenvalue spectrum of the onedimensional (quantum) transfer matrix provides a holographic description of the dispersion relations of elementary excitations in the full twodimensional quantum system. This holds true even in systems with topological order, where the elementary excitations are anyons. The presense of topological order gives rise to particular (virtual) symmetries of the transfer matrix. By carefully studying the manifestation of these symmetries in the fixedpoint subspace, we were able to relate the different topological sectors of the physical excitations to corresponding topologically nontrivial symmetry sectors (domain walls) at the virtual level. This shows, in particular, that the existence of anyon excitations requires a particular type of symmetry breaking of the doubled virtual symmetry in the fixedpoint subspace of the transfer matrix, whereas topological phase transitions give rise to a fixedpoint subspace with a larger or smaller degeneracy.
While these results might be reminiscent of the closely related bulkedge correspondence observed in chiral topological phases^{45,46}, we would like point out the subtle differences. In the PEPS formalism, the properties of the edge states are determined by the fixed points of the transfer matrix^{47}, whereas here we explicitly consider the complete (longdistance) spectrum of the transfer matrix. In addition, the framework for characterizing topological order in PEPS using MPOs, which is of central importance for our results, has so far only been made explicit for the nonchiral string net models, and it remains to be clarified how the recently discovered chiral PEPS^{48,49,50} fit within this framework.
This technique holds a powerful potential for studying fundamental questions of topological order and topological phase transitions. While we have studied transfer matrices originating from a tensor network representation of the ground state, the results presented in this paper should generalize to the fullquantum transfer matrix obtained from representing the ground state as an imaginary timepath integral. Whereas the exact path integral representation can have a gauge theory as virtual boundary, the PEPS truncation will eliminate the gauge degrees of freedom. Correspondingly, the local order parameter measuring the symmetry breaking transitions at the PEPS virtual level will map to a string operator in the temporal direction of the fullpath integral.
Additional information
How to cite this article: Haegeman, J. et al. Shadows of anyons and the entanglement structure of topological phases. Nat. Commun. 6:8284 doi: 10.1038/ncomms9284 (2015).
References
 1
Hastings, M. B. An area law for onedimensional quantum systems. J. Stat. Mech. 08, 24 (2007).
 2
Eisert, J., Cramer, M. & Plenio, M. B. Area laws for the entanglement entropy—a review. Rev. Mod. Phys. 82, 277–306 (2010).
 3
Li, H. & Haldane, F. D. M. Entanglement spectrum as a generalization of entanglement entropy: Identification of topological order in nonabelian fractional quantum hall effect states. Phys. Rev. Lett. 101, 010504 (2008).
 4
Cirac, J. I., Poilblanc, D., Schuch, N. & Verstraete, F. Entanglement spectrum and boundary theories with projected entangledpair states. Phys. Rev. B 83, 245134 (2011).
 5
Schuch, N., Poilblanc, D., Cirac, J. I. & PérezGarca, D. Topological order in the projected entangledpair states formalism: Transfer operator and boundary hamiltonians. Phys. Rev. Lett. 111, 090501 (2013).
 6
Fannes, M., Nachtergaele, B. & Werner, R. F. Finitely correlated states on quantum spin chains. Comm. Math. Phys. 144, 443–490 (1992).
 7
Östlund, S. & Rommer, S. Thermodynamic limit of density matrix renormalization. Phys. Rev. Lett. 75, 3537–3540 (1995).
 8
Verstraete, F. & Cirac, J. I. Renormalization algorithms for QuantumMany Body Systems in two and higher dimensions. Preprint at http://arxiv.org/abs/condmat/0407066 (2004).
 9
Verstraete, F., Murg, V. & Cirac, J. I. Matrix product states, projected entangled pair states, and variational renormalization group methods for quantum spin systems. Adv. Phys 57, 143–224 (2008).
 10
Suzuki, M. Quantum transfermatrix method and thermoquantum dynamics. Physica A 321, 334–339 (2003).
 11
Zauner, V. et al. Transfer matrices and excitations with matrix product states. New. J. Phys. 17, 053002 (2015).
 12
Affleck, I., Kennedy, T., Lieb, E. H. & Tasaki, H. Rigorous results on valencebond ground states in antiferromagnets. Phys. Rev. Lett. 59, 799–802 (1987).
 13
Affleck, I., Kennedy, T., Lieb, E. H. & Tasaki, H. Valence bond ground states in isotropic quantum antiferromagnets. Comm. Math. Phys. 115, 477–528 (1988).
 14
Haegeman, J. et al. Elementary excitations in gapped quantum spin systems. Phys. Rev. Lett. 111, 080401 (2013).
 15
Lieb, E. H. & Liniger, W. Exact analysis of an interacting bose gas. I. the general solution and the ground state. Phys. Rev. 130, 1605–1616 (1963).
 16
Mandelstam, S. Soliton operators for the quantized sineGordon equation. Phys. Rev. D 11, 3026 (1975).
 17
Haegeman, J. et al. Variational matrix product ansatz for dispersion relations. Phys. Rev. B 85, 100408 (2012).
 18
Draxler, D. et al. Particles, holes, and solitons: a matrix product state approach. Phys. Rev. Lett. 111, 020402 (2013).
 19
Milsted, A., Haegeman, J. & Osborne, T. J. Matrix product states and variational methods applied to critical quantum field theory. Phys. Rev. D 88, 085030 (2013).
 20
Dusuel, S., Kamfor, M., Orús, R., Schmidt, K. P. & Vidal, J. Robustness of a perturbed topological phase. Phys. Rev. Lett. 106, 107203 (2011).
 21
Schulz, M. D., Dusuel, S., Orús, R., Vidal, J. & Schmidt, K. P. Breakdown of a perturbed topological phase. New. J. Phys. 14, 025005 (2012).
 22
Levin, M. A. & Wen, X.G. Stringnet condensation: A physical mechanism for topological phases. Phys. Rev. B 71, 45110 (2005).
 23
Buerschaper, O., Aguado, M. & Vidal, G. Explicit tensor network representation for the ground states of stringnet models. Phys. Rev. B 79, 85119 (2009).
 24
Gu, Z.C., Levin, M., Swingle, B. & Wen, X.G. Tensorproduct representations for stringnet condensed states. Phys. Rev. B 79, 85118 (2009).
 25
Bais, F. A. & Slingerland, J. K. Condensateinduced transitions between topologically ordered phases. Phys. Rev. B 79, 045316 (2009).
 26
Burnell, F. J., Simon, S. H. & Slingerland, J. K. Phase transitions in topological lattice models via topological symmetry breaking. New. J. Phys. 14, 015004 (2012).
 27
Pirvu, B., Haegeman, J. & Verstraete, F. Matrix product state based algorithm for determining dispersion relations of quantum spin chains with periodic boundary conditions. Phys. Rev. B 85, 35130 (2012).
 28
Schuch, N., Cirac, I. & PerezGarcia, D. PEPS as ground states: Degeneracy and topology. Ann. Phys 325, 2153–2192 (2010).
 29
Buerschaper, O. Twisted injectivity in projected entangled pair states and the classification of quantum phases. Ann. Phys 351, 447–476 (2014).
 30
Burak şahinoğlu, M. et al. Characterizing Topological Order with Matrix Product Operators. Preprint at http://arxiv.org/abs/1409.2150 (2014).
 31
Hastings, M. B. & Wen, X.G. Quasiadiabatic continuation of quantum states: The stability of topological groundstate degeneracy and emergent gauge invariance. Phys. Rev. B 72, 45141 (2005).
 32
Tu, H.H., Zhang, Y. & Qi, X.L. Momentum polarization: An entanglement measure of topological spin and chiral central charge. Phys. Rev. B. 88, 195412 (2013).
 33
Papanikolaou, S., Raman, K. S. & Fradkin, E. Topological phases and topological entropy of twodimensional systems with finite correlation length. Phys. Rev. B 76, 224421 (2007).
 34
Castelnovo, C. & Chamon, C. Quantum topological phase transition at the microscopic level. Phys. Rev. B 77, 054433 (2008).
 35
Haegeman, J., Van Acoleyen, K., Schuch, N., Cirac, J. I. & Verstraete, F. Gauging quantum states: from global to local symmetries in manybody systems. Phys. Rev. X 5, 011024 (2015).
 36
Zhou, H.Q., Orus, R. & Vidal, G. Ground state fidelity from tensor network representations. Phys. Rev. Lett. 100, 80601 (2008).
 37
Zhou, H.Q. & Barjaktarevič, J. P. Fidelity and quantum phase transitions. J. Phys. A 41, 2001 (2008).
 38
Abasto, D. F., Hamma, A. & Zanardi, P. Fidelity analysis of topological quantum phase transitions. Phys. Rev. A 78, 010301 (2008).
 39
Verstraete, F., Wolf, M. M., PerezGarcia, D. & Cirac, J. I. Criticality, the area law, and the computational power of projected entangled pair states. Phys. Rev. Lett. 96, 220601 (2006).
 40
Anderson, P. The resonating valence bond state in La2CuO4 and superconductivity. Science 235, 1196 (1987).
 41
Schuch, N., Poilblanc, D., Cirac, J. I. & PérezGarca, D. Resonating valence bond states in the peps formalism. Phys. Rev. B 86, 115108 (2012).
 42
Poilblanc, D. & Schuch, N. Simplex spin liquids on the kagome lattice with projected entangled pair states: Spinon and vison coherence lengths, topological entropy, and gapless edge modes. Phys. Rev. B 87, 140407 (2013).
 43
Read, N. & Chakraborty, B. Statistics of the excitations of the resonatingvalencebond state. Phys. Rev. B 40, 7133–7140 (1989).
 44
Kivelson, S. Statistics of holons in the quantum hardcore dimer gas. Phys. Rev. B 39, 259–264 (1989).
 45
Wen, X. G. Chiral luttinger liquid and the edge excitations in the fractional quantum hall states. Phys. Rev. B 41, 12838–12844 (1990).
 46
Wen, X.G., Wu, Y.S. & Hatsugai, Y. Chiral operator product algebra and edge excitations of a fractional quantum hall droplet. Nucl. Phys. B 422, 476–494 (1994).
 47
Yang, S. et al. Edge theories in projected entangled pair state models. Phys. Rev. Lett. 112, 036402 (2014).
 48
Dubail, J. & Read, N. Tensor network trial states for chiral topological phases in two dimensions. ArXiv eprints (2013).
 49
Wahl, T. B., Tu, H.H., Schuch, N. & Cirac, J. I. Projected entangledpair states can describe chiral topological states. Phys. Rev. Lett. 111, 236805 (2013).
 50
Yang, S., Wahl, T. B., Tu, H.H., Schuch, N. & Cirac, J. I. Chiral projected entangledpair state with topological order. Phys. Rev. Lett. 114, 106803 (2015).
Acknowledgements
This work was initiated during the programme on ‘Quantum Hamiltonian Complexity’ held at the Simons Institute for the Theory of Computing. We acknowledge discussions with Nick Bultinck, Ignacio Cirac, Michaël Mariën, Burak Sahinoglu, Steve Simon, Joost Slingerland, Laurens Vanderstraeten, Karel Van Acoleyen and Dominic Williamson. We gratefully acknowledge support by EU grants SIQS, QUERG and WASCOSYS, the Austrian FWF SFB grants FoQuS and ViCoM. We further acknowledge the support from the Research Foundation Flanders (J.H., F.V.) and the Alexander von Humboldt foundation (N.S.), and computational resources provided by JARAHPC via grants jara0084 and jara0092.
Author information
Affiliations
Contributions
All authors contributed equally to all aspects of this work.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures 13, Supplementary Notes 12 and Supplementary References (PDF 462 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Haegeman, J., Zauner, V., Schuch, N. et al. Shadows of anyons and the entanglement structure of topological phases. Nat Commun 6, 8284 (2015). https://doi.org/10.1038/ncomms9284
Received:
Accepted:
Published:
Further reading

NonHermitian effects of the intrinsic signs in topologically ordered wavefunctions
Communications Physics (2020)

Donor polymer design enables efficient nonfullerene organic solar cells
Nature Communications (2016)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.