In this thesis, quantum many-body systems are studied by means of tensor network methods. After introducing the quantum many-body problem and discussing paradigmatic models as well as established methods for their solution, the focus is set on a practical approach to tensor network methods, always keeping in mind the realistic numerical applicability of the tools devised. In order to achieve this objective, we focus on loop-free network geometries, which allow to boost the computational performance to the greatest possible extent. The developed algorithms are described and benchmarked in detail, starting from their elementary constituents, namely the multidimensional arrays called tensors which are at the heart of any tensor network algorithm. We then apply the developed simulation techniques to a diverse collection of quantum many-body models of current research interest. We are interested, in particular, in a characterization of the ground-state properties of these models, but we also investigate problems beyond the zero-temperature equilibrium setting, such as time evolution or finite-temperature behavior. Specifically, novel findings for the phase diagrams of the disordered Bose-Hubbard chain and of the Heisenberg chain with additional biquadratic interactions are presented, as well as compelling new evidence for the existence of fractional quantum Hall states in the interacting Harper-Hofstadter model.
In this thesis, quantum many-body systems are studied by means of tensor network methods. After introducing the quantum many-body problem and discussing paradigmatic models as well as established methods for their solution, the focus is set on a practical approach to tensor network methods, always keeping in mind the realistic numerical applicability of the tools devised. In order to achieve this objective, we focus on loop-free network geometries, which allow to boost the computational performance to the greatest possible extent. The developed algorithms are described and benchmarked in detail, starting from their elementary constituents, namely the multidimensional arrays called tensors which are at the heart of any tensor network algorithm. We then apply the developed simulation techniques to a diverse collection of quantum many-body models of current research interest. We are interested, in particular, in a characterization of the ground-state properties of these models, but we also investigate problems beyond the zero-temperature equilibrium setting, such as time evolution or finite-temperature behavior. Specifically, novel findings for the phase diagrams of the disordered Bose-Hubbard chain and of the Heisenberg chain with additional biquadratic interactions are presented, as well as compelling new evidence for the existence of fractional quantum Hall states in the interacting Harper-Hofstadter model.
The study of interacting quantum many-body systems is one of the central challenges in modern physics. Interestingly, by understanding the physical properties of such systems from a microscopic description, one can explain many complex quantum phenomena such as the high-temperature superconductivity, fractional quantum Hall effect, etc. However, due to high complexity of these systems, theoretical investigations such as analytical and perturbation theories are limited to studying only a small subset of problems. In the past few decades, numerical simulation has emerged as an indispensable method for studying many phenomena, and has led to immense success and progress in recent years. The main obstacle in addressing quantum many-body systems is that the number of parameters required to represent a state of the system scales exponentially with the system size. However, one hardly needs the full Hilbert space but only a small subspace to describe physical properties of the system. From the point of view in quantum information theory, quantum entanglement is one of the main factors that is responsible for the complexity of studying quantum many-body systems. In particular, the amount of quantum entanglement in such system quantifies the correlation between constituents and therefore it determines how big the subspace is required to represent the system. In this thesis, we are interested in using tensor network states, a new emerging numerical technique to study dynamical and ground state properties of the system. The key ingredient is that a tensor network state which consists of a fairly small number of parameters can represent very well the complete wave function of a system. In one spatial dimension, Matrix Product State (MPS) has become a standard tool, combined with related algorithms such as Time Evolving Block Decimation algorithm (TEBD) to study dynamical properties of a locally-perturbed quantum lattice system. In this thesis, we introduce a technique called ``infinite boundary conditions'' that allows us to use a finite-size MPS to represent the state of an infinite-size system. As a result, one can simulate dynamics of a system for a longer time with a significant computational cost reduction compared to traditional methods. In two spatial dimensions, the tensor network known as Projected Entangled Pair States (PEPS) is emerging as a useful approach. However algorithms for optimizing the matrix elements of these tensors are complex and not very efficient. In particular, the efficiency of the optimization depends critically on the quality of the so-called ``environment'', which is expensive to compute. We have introduced a new method to accelerate the convergence of calculating the PEPS ground state by recycling the environment tensors for many optimization steps, vastly reducing the number of times that the full environment needs to be calculated.
Emergent phenomena of interacting quantum many-body systems arguably pose some of the greatest challenges to modern theoretical physics. Many of these systems are characterized by drastic changes in their properties as they are driven into a regime where quantum effects and interactions become relevant. Examples are the fractional quantum Hall effect and high-temperature superconductors, to name just two. New theoretical tools are necessary to study such systems, as even simple model Hamiltonians are not fully understood after decades of intense research. Numerical simulations are emerging as an indispensable tool alongside more traditional analytic approaches and have seen immense progress in recent years. In this thesis, we study tensor network states, a new class of numerical methods that have been developed in recent years drawing on insights from quantum information theory. Combining renormalization group methods with knowledge about the entanglement structure of ground states of quantum systems, tensor network states aim to efficiently describe ground states of strongly correlated quantum systems. Center stage is taken by the projected entangled-pair states. This ansatz for two-dimensional systems is discussed in detail and its capabilities are assessed by comparison to established methods and by exploring a frustrated system. We then develop a general formalism to exploit Abelian symmetries in a tensor network algorithm and demonstrate its application to projected entangled-pair states. A large part of this thesis is devoted to the discussion of a class of supersymmetric models for interacting lattice fermions. After introducing the concept of supersymmetry in quantum mechanics and an extensive overview of conformal field theory and its relation to numerical simulations, we study the critical theory of the supersymmetric model on the square ladder. We use two tensor network state algorithms, the density-matrix renormalization group and multi-scale entanglement renormalization along with exact diagonalization to explore this challenging system. In the final two chapters, we discuss two applications of tensor network states, namely simulations of the SU(3) Heisenberg model in two dimensions using two different tensor network state algorithms with the goal of illuminating the nature of its ground state, and the indistinguishability, a measure to detect phase transitions and classify wave functions which is particularly suitable for calculations based on tensor network states.
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
Strongly correlated quantum systems are one of the most interesting and challenging phenomena in modern condensed matter physics. Unfortunately, the theoretical access even to their low energy properties is limited by the exponential growth of the Hilbert space in a many-body setting. To mitigate this computational barrier, tensor network methods are one of the most competitive approaches to classically simulate complex quantum systems. In this thesis we will discuss recent developments in the field of tensor networks to efficiently and accurately simulate two-dimensional quantum many-body states in the thermodynamic limit. Hereby, we introduce the variational update globally optimizing the wave-function, encoded as tensor network in the form of an infinite projected entangled pair state (iPEPS), in order to obtain ground states of many-body Hamiltonian models. We focus in the discussion on the well-established method to approximately contract this two-dimensional structure, the corner transfer matrix renormalization group (CTMRG), and discuss its native application on square lattice and triangular lattice based models. Furthermore, we report the development of a refined numerical algorithm, the split-CTMRG, to significantly reduce its computational complexity, pushing the limits of feasible simulations within modern hardware capabilities. Finally, the method advancements are applied to a currently widely studied model of frustrated magnetism, the Heisenberg maple-leaf antiferromagnet. We are able to find a large phase of vanishing magnetization in contradiction to previous works and, to the best of our knowledge, the lowest reported variational ground state energy at the isotropic point of the model. As an outlook, we describe a potential adaption of the previously developed methods to variationally simulate two-dimensional many-body thermal states.
We present a compendium of numerical simulation techniques, based on tensor network methods, aiming to address problems of many-body quantum mechanics on a classical computer. The core setting of this anthology are lattice problems in low spatial dimension at finite size, a physical scenario where tensor network methods, both Density Matrix Renormalization Group and beyond, have long proven to be winning strategies. Here we explore in detail the numerical frameworks and methods employed to deal with low-dimensional physical setups, from a computational physics perspective. We focus on symmetries and closed-system simulations in arbitrary boundary conditions, while discussing the numerical data structures and linear algebra manipulation routines involved, which form the core libraries of any tensor network code. At a higher level, we put the spotlight on loop-free network geometries, discussing their advantages, and presenting in detail algorithms to simulate low-energy equilibrium states. Accompanied by discussions of data structures, numerical techniques and performance, this anthology serves as a programmer’s companion, as well as a self-contained introduction and review of the basic and selected advanced concepts in tensor networks, including examples of their applications.
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.