Tensor network methods have progressed from variational techniques based on matrix-product states able to compute properties of one-dimensional condensed-matter lattice models into methods rooted in more elaborate states, such as projected entangled pair states aimed at simulating the physics of two-dimensional models. In this work, we advocate the paradigm that for two-dimensional fermionic models, matrix-product states are still applicable to significantly higher accuracy levels than direct embeddings into one-dimensional systems allow for. To do so, we exploit schemes of fermionic mode transformations and overcome the prejudice that one-dimensional embeddings need to be local. This approach takes the insight seriously that the suitable exploitation of both the manifold of matrix-product states and the unitary manifold of mode transformations can more accurately capture the natural correlation structure. By demonstrating the residual low levels of entanglement in emerging modes, we show that matrix-product states can describe ground states strikingly well. The power of the approach is exemplified by investigating a phase transition of spinless fermions for lattice sizes up to $10\ifmmode\times\else\texttimes\fi{}10$.
Tensor network methods have progressed from variational techniques based on matrix-product states able to compute properties of one-dimensional condensed-matter lattice models into methods rooted in more elaborate states, such as projected entangled pair states aimed at simulating the physics of two-dimensional models. In this work, we advocate the paradigm that for two-dimensional fermionic models, matrix-product states are still applicable to significantly higher accuracy levels than direct embeddings into one-dimensional systems allow for. To do so, we exploit schemes of fermionic mode transformations and overcome the prejudice that one-dimensional embeddings need to be local. This approach takes the insight seriously that the suitable exploitation of both the manifold of matrix-product states and the unitary manifold of mode transformations can more accurately capture the natural correlation structure. By demonstrating the residual low levels of entanglement in emerging modes, we show that matrix-product states can describe ground states strikingly well. The power of the approach is exemplified by investigating a phase transition of spinless fermions for lattice sizes up to 10\ifmmode×\else\texttimes\fi10.
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.
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.
Tensor network states, and in particular projected entangled pair states (PEPS) have been a strong ansatz for the variational study of complicated quantum many-body systems, thanks to their built-in entanglement entropy area law. In this work, we use a special kind of PEPS---gauged Gaussian fermionic PEPS (GGFPEPS)---to find the ground state of 2+1 dimensional pure Z2 lattice gauge theories for a wide range of coupling constants. We do so by combining PEPS methods with Monte-Carlo computations, allowing for efficient contraction of the PEPS and computation of correlation functions. Previously, such numerical computations involved the calculation of the Pfaffian of a matrix scaling with the system size, forming a severe bottleneck; in this work we show how to overcome this problem. This paves the way for applying the method we propose and benchmark here to other gauge groups, higher dimensions, and models with fermionic matter, in an efficient, sign-problem-free way.
This thesis treats the classical simulation of strongly-interacting many-body quantum-mechanical systems in more than one dimension using matrix product states and the more general tensor product states. Contrary to classical systems, quantum many-body systems possess an exponentially larger number of degrees of freedom, thereby significantly complicating their numerical treatment on a classical computer. For this thesis two different representations of quantum many-body states were employed. The first, the so-called matrix product states (MPS) form the basis for the extremely successful density matrix renormalization group (DMRG) algorithm. While originally conceived for one-dimensional systems, MPS are in principle capable of describing arbitrary quantum many-body states. Using concepts from quantum information theory it is possible to show that MPS provide a representation of one-dimensional quantum systems that scales polynomially in the number of particles, therefore allowing an efficient simulation of one-dimensional systems on a classical computer. One of the key results of this thesis is that MPS representations are indeed efficient enough to describe even large systems in two dimensions, thereby enabling the simulation of such systems using DMRG. As a demonstration of the power of the DMRG algorithm, it is applied to the Heisenberg antiferromagnet with spin S=1/2 on the kagome lattice. This model's ground state has long been under debate, with proposals ranging from static spin configurations to so-called quantum spin liquids, states where quantum fluctuations destroy conventional order and give rise to exotic quantum orders. Using a fully SU(2)-symmetric implementation allowed us to handle the exponential growth of entanglement and to perform a large-scale study of this system, finding the ground state for cylinders of up to 700 sites. Despite employing a one-dimensional algorithm for a two-dimensional system, we were able to compute the spin gap (i.e. the energy gap to the first spinful excitation) and study the ground state properties, such as the decay of correlation functions, the static spin structure factors, and the structure and distribution of the nearest-neighbor spin-spin correlations. Additionally, by applying a new tool from quantum information theory, the topological entanglement entropy, we could also with high confidence demonstrate the ground state of this model to be the elusive gapped Z2 quantum spin liquid with topological order. To complement this study, we also considered the extension of MPS to higher dimensions, known as tensor product states (TPS). We implemented an optimization algorithm exploiting symmetries for this class of states and applied it to the bilinear-biquadratic-bicubic Heisenberg model with spin S=3/2 on the z=3 Bethe lattice. By carefully analyzing the simulation data we were able to determine the presence of both conventional and symmetry-protected topological order in this model, thereby demonstrating the analytically predicted existence of the Haldane phase in higher dimensions within an extended region of the phase diagram. Key properties of this symmetry-protected topological order include a doubling of the levels in the entanglement spectrum and the presence of edge spins, both of which were confirmed in our simulations. This finding simultaneously validated the applicability of the novel TPS algorithms to the search for exotic order.
Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.