Understanding many-body quantum systems is one of the most challenging problems in contemporary condensed-matter physics. Tensor network representation of quantum states and operators are taking central stage in this pursuit and beyond. They prove to be a powerful numerical and conceptual tool, and indeed a new language altogether. This thesis investigates various aspects of these representations by focusing on two specific problems: the first half of the thesis is devoted to examining how 'stable' a tensor network representation is for two-dimensional quantum states with topological order, and the second half explores the representability of various unitary loop operators with tensor networks. In the numerical usage of the tensor networks, the tensor is varied as to find the representation of the ground states of the given Hamiltonian. In chapter two and three of this thesis we show that such a numerical program for topological phases can be 'ill-posed'. We show that tensor network can be an unstable representation for a topological phase: even an infinitesimal variation in the representation results in the loss of topological order, completely or partially. We diagnose this problem by identifying the exact causes of this instability, and find that it is only tensor variations in certain directions that result in instability, because they result in the condensation of bosonic quasi-particles of the phase. Such unstable variations are characterized by two properties: (1) they can replace a tensor in the tensor network without making the network collapse, and (2) their presence in the network represents the presence of a non-trivial topological charge. We prove that the general tensor representation of all string-net models suffer with such instabilities. We propose an exact mathematical operator to project out all such unstable variations and show its efficacy for a few models by direct calculations. Such an operator can be useful in numerical programs involving such tensor representations. We also point out that such variations play a crucial role in simulating topological phase transitions and their presence can be vital in an accurate simulation. In chapter four and five of this thesis we focus on the representability of unitary loop operators by tensor networks. Such operators not only provide an important tool in the study of dynamical process in one-dimensional systems, but also in understanding and classification of symmetry protected topological phases in two dimensions. To characterize all such operators, we find a necessary and sufficient condition for any loop tensor network operator of a given length to represent a unitary operator. In particular, it is shown that all unitary operators that map local operators to local operators (locality-preserving) can always be represented by a tensor network. Locality-preserving unitary loop operators are classified by a rational index called the GNVW index defined in Ref. [1] which measures how much information 'flows' along the loop. We define Rank-Ratio index for tensor network operators and show that it is completely equivalent to the GNVW index. Therefore, GNVW index of a unitary operator can be easily extracted from its tensor network representation. We find that, other than representing locality-preserving unitary maps, tensor networks can also represent unitary operators that map local operators to global (non-local) operators. These tensor network operators are found to have a long-ranged order similar to tensors that represent topological tensor network states in two dimensions.
Understanding many-body quantum systems is one of the most challenging problems in contemporary condensed-matter physics. Tensor network representation of quantum states and operators are taking central stage in this pursuit and beyond. They prove to be a powerful numerical and conceptual tool, and indeed a new language altogether. This thesis investigates various aspects of these representations by focusing on two specific problems: the first half of the thesis is devoted to examining how 'stable' a tensor network representation is for two-dimensional quantum states with topological order, and the second half explores the representability of various unitary loop operators with tensor networks. In the numerical usage of the tensor networks, the tensor is varied as to find the representation of the ground states of the given Hamiltonian. In chapter two and three of this thesis we show that such a numerical program for topological phases can be 'ill-posed'. We show that tensor network can be an unstable representation for a topological phase: even an infinitesimal variation in the representation results in the loss of topological order, completely or partially. We diagnose this problem by identifying the exact causes of this instability, and find that it is only tensor variations in certain directions that result in instability, because they result in the condensation of bosonic quasi-particles of the phase. Such unstable variations are characterized by two properties: (1) they can replace a tensor in the tensor network without making the network collapse, and (2) their presence in the network represents the presence of a non-trivial topological charge. We prove that the general tensor representation of all string-net models suffer with such instabilities. We propose an exact mathematical operator to project out all such unstable variations and show its efficacy for a few models by direct calculations. Such an operator can be useful in numerical programs involving such tensor representations. We also point out that such variations play a crucial role in simulating topological phase transitions and their presence can be vital in an accurate simulation. In chapter four and five of this thesis we focus on the representability of unitary loop operators by tensor networks. Such operators not only provide an important tool in the study of dynamical process in one-dimensional systems, but also in understanding and classification of symmetry protected topological phases in two dimensions. To characterize all such operators, we find a necessary and sufficient condition for any loop tensor network operator of a given length to represent a unitary operator. In particular, it is shown that all unitary operators that map local operators to local operators (locality-preserving) can always be represented by a tensor network. Locality-preserving unitary loop operators are classified by a rational index called the GNVW index defined in Ref. [1] which measures how much information 'flows' along the loop. We define Rank-Ratio index for tensor network operators and show that it is completely equivalent to the GNVW index. Therefore, GNVW index of a unitary operator can be easily extracted from its tensor network representation. We find that, other than representing locality-preserving unitary maps, tensor networks can also represent unitary operators that map local operators to global (non-local) operators. These tensor network operators are found to have a long-ranged order similar to tensors that represent topological tensor network states in two dimensions.
The tensor network representation of many-body quantum states, given by local tensors, provides a promising numerical tool for the study of strongly correlated topological phases in two dimensions. However, the representation may be vulnerable to instabilities caused by small variations in the local tensors. For example, the topological order in the tensor network representations of the toric code ground state has been shown in Chen, Zeng, Gu, Chuang, and Wen, Phys. Rev. B 82, 165119 (2010) to be unstable if the variations break certain Z2 symmetry of the tensor. In this work, we ask whether other types of topological orders suffer from similar kinds of instability and if so, what is the underlying physical mechanism and whether we can protect the order by enforcing certain symmetries on the tensor. We answer these questions by showing that the tensor network representations of all string-net models are indeed unstable, but the matrix product operator (MPO) symmetries of the tensors identified in \ifmmode \mbox{\c{S}}\else \c{S}\fi{}ahino\ifmmode \breve{g}\else \u{g}\fi{}lu, Williamson, Bultinck, Mari\"en, Haegeman, Schuch, and Verstraete, arXiv:1409.2150 can help to protect the order. In particular, we show that a subset of variations that break the MPO symmetries lead to instability by inducing the condensation of bosonic quasiparticles, which destroys the topological order in the wave function. Therefore such variations must be forbidden for the encoded topological order to be reliably extracted from the local tensors. On the other hand, if a tensor network based variational algorithm is used to simulate the phase transition due to boson condensation, such variation directions may prove important to access the continuous transition correctly.
Strongly correlated systems exhibit phenomena -- such as high-T_c superconductivity or the fractional quantum Hall effect -- that are not explicable by classical and semi-classical methods. Moreover, due to the exponential scaling of the associated Hilbert space, solving the proposed model Hamiltonians by brute-force numerical methods is bound to fail. Thus, it is important to develop novel numerical and analytical methods that can explain the physics in this regime. Tensor Network states are quantum many-body states that help to overcome some of these difficulties by defining a family of states that depend only on a small number of parameters. Their use is twofold: they are used as variational ansatzes in numerical algorithms as well as providing a framework to represent a large class of exactly solvable models that are believed to represent all possible phases of matter. The present thesis investigates mathematical properties of these states thus deepening the understanding of how and why Tensor Networks are suitable for the description of quantum many-body systems. It is believed that tensor networks can represent ground states of local Hamiltonians, but how good is this representation? This question is of fundamental importance as variational algorithms based on tensor networks can only perform well if any ground state can be approximated efficiently in such a way. While any state can be written as a tensor network state, the number of parameters needed for the description might be too large. This is not the case for one-dimensional systems: only a few parameters are required to have a good approximation of their ground states; that, in turn, allows for numerical algorithms based on tensor networks performing well. The situation in two dimensions is somewhat more complicated, but it is known that ground states of local Hamiltonians can be expressed as tensor networks with sub-exponentially many parameters. In the present thesis, we improve on these existing bounds strengthening the claim that the language of tensor networks is suitable to describe many-body systems. Another central question is how symmetries of the system such as translational invariance, time-reversal symmetry or local unitary symmetry can be reflected in tensor networks. This question is important as systems appearing in nature might intrinsically possess certain symmetries; on one hand, understanding these symmetries simplifies the description of these systems. On the other hand, the presence of symmetries leads to the appearance of novel phases -- symmetry-protected topological (SPT) order, -- and tensor networks provide the right language to classify these phases. In one dimension and for certain classes of two-dimensional tensor networks (states generated by so-called injective tensors) it is well understood how symmetries of the state can be described. A general framework, however, has yet to be developed. In the present thesis, we contribute to the development of the theory in two ways. We first investigate the question for injective tensors, and generalize the existing proof for any geometry including the hyperbolic geometry used in the AdS/CFT correspondence. Second, we introduce a class of tensor network states that include previously known examples of states exhibiting SPT order. We show how symmetries are reflected in these states thus deepening the understanding of SPT order in two dimensions.
The classical simulation of many-body quantum systems is an essential tool in understanding many fundamental aspects of condensed matter physics. But a major obstacle arises from the number of degrees of freedom involved in describing such systems, which is exponential in the system size. Recently, however, a class of numerical techniques based on structures called "tensor networks" has emerged, which allows many "typical" quantum states (such as the ground states of gapped, local Hamiltonians) to be represented much more efficiently. In this work we extend and apply these techniques to consider several central topics in quantum many-body physics. After reviewing the relevant background material from the field of tensor networks and tensor network states, we demonstrate a method for computing high order moments and cumulants of operators with respect to such states, including the so-called "Binder cumulant," a powerful tool for detecting phase transitions. Next, we employ tensor network algorithms to characterize the ground state phase diagram of a quantum spin model, including both symmetry-breaking phases and symmetry protected topological order, and find a signicant variety of phases and phase transitions. Finally, we consider the entanglement properties of quantum states exhibiting many-body localization, using a combination of exact diagonalization and tensor network techniques. | 231 pages
This thesis is devoted to the study of quantum many-body systems. This investigation is performed in the framework of Matrix Product States (MPS) and their generalization to higher dimensions, Projected Entangled Pair States (PEPS). In the first part of the work, we discuss the mathematical properties of such tensor network states in depth. In the first chapters we deal with one-dimensional systems, for which we use the connections between MPS and completely positive maps to prove several new results, such as the quantum version of the Wielandt's inequality or the construction of locally invariant states under a symmetry group. Chapters 5-7 are dedicated to higher-dimensional systems, for which we provide the conditions to construct invariant PEPS. The second part of the thesis is dedicated to applications in condensed matter. In chapter 8 we provide several methods of constructing quasi-solvable Hamiltonians with two-body interactions, while in chapter 9 we show that these tensor network states can be used as a laboratory for theoretical condensed matter in, for instance, the characterization of the string-order, the generalization of the Lieb-Schultz-Mattis theorem, the demonstration of new theorems relating entanglement to magnetisation or to long-range interactions, etc.
It appears that the tensor network ansatz, while not quite complete, is an efficient coordinate system for the tiny subset of a many-body Hilbert space which can be realized as a low energy state of a local Hamiltonian.However, we don't fully understand precisely which phases are captured by the tensor network ansatz, how to compute their physical observables (even numerically), or how to compute a tensor network representation for a ground state given a microscopic Hamiltonian. These questions are algorithmic in nature, but their resolution is intimately related to understanding the nature of quantum entanglement in many-body systems.For this reason it is useful to compute the tensor network representation of various `model' wavefunctions representative of different phases of matter; this allows us to understand how the entanglement properties of each phase are expressed in the tensor network ansatz, and can serve as test cases for algorithm development.Condensed matter physics has many illuminating model wavefunctions, such as Laughlin's celebrated wave function for the fractional quantum Hall effect, the Bardeen-Cooper-Schrieffer wave function for superconductivity, and Anderson's resonating valence bond ansatz for spin liquids.This thesis presents some results on exact tensor network representations of these model wavefunctions. In addition, a tensor network representation is given for the time evolution operator of a long-range one-dimensional Hamiltonian, which allows one to numerically simulate the time evolution of power-law interacting spin chains as well as two-dimensional strips and cylinders.
The focus of our investigations in this thesis is quantum spin systems in two dimensions. We examine phase transitions between topologically distinct phases of quantum matter using the framework of tensor network states. We study the phenomena of anyon condensation and confinement in the context of tensor network states where these two notions offer a robust probe to characterize the universal features of quantum phase transitions. We do a comprehensive study of the numerical methods enabled by the tensor network formalism for the study of quantum phase transitions. We map out the phase diagram of certain exotic phases of quantum matter (namely the D(mathbbZ4) quantum double, the toric code, and the double semion model) and identify the order and the universality classes of the phase transitions between these distinct topological phases. We analyze the phase boundaries of the toric code and the double semion model which exhibit first and second order phase transitions. We investigate the class of mathbbZ2-invariant tensor network states while taking in to account the spin rotation and lattice symmetries. The resulting tensor network states allow us to map the phase diagram of mathbbZ2 spin liquids. The behavior of the system at the phase boundaries of mathbbZ2 spin liquid is governed by the condensation of spinons and visions. We also present our preliminary findings regarding the mathbbZ4-invariant tensors with SU(2) symmetry. This approach enables us to study the spin liquid states of the toric code and the double semion model in a unified framework. We give a local tensor description for the approximate ground states of the Heisenberg antiferromagnet on the kagome lattice by using the variational manifolds of just three and five parameters in the gapped mathbbZ2 spin liquid phase. The approximation of the ground states we construct have an energy density that is remarkably close to the results from the state-of-the-art density matrix renormalization group and the exact diagonalization method. By analyzing the deconfinement of anyonic excitations, we also present our findings regarding the vicinity of the variational ground state in gapped mathbbZ2 spin liquid phase.