The simulation and numerical study of large, strongly correlated quantum systems containing Fermions or using real-time evolution in finite dimensions is still an essentially unsolved problem, primarily due to the exponential growth of the Hilbert state space with system size and the occurrence of the so-called sign problem in Monte Carlo studies. In this area, the use of tensor-network methods, for one-dimensional systems chief among them the density matrix renormalisation group (DMRG) and matrix-product states (MPS), has grown in importance in recent years. This thesis first recapitulates the use of non-abelian symmetries such as SU(2)-Spin in arbitrary tensor networks with an extensive review of the published literature including detailed algorithms and implementation hints. Implementing such symmetries can lead to a considerably more efficient representation of states in the tensor network. This part is intended to be suitable as an implementation-oriented introduction to tensor networks in general and the implementation of non-abelian symmetries in particular. Second, it introduces a series of technical improvements for the MPS methods. These improvements include a faster convergence scheme for MPS-DMRG, a systematic approach to the construction of matrix-product operators and an improved Krylov time evolution method as well as the combination of several well-known techniques into a single tensor network toolkit, SYTEN. The effectiveness of these improvements is demonstrated in numerical examples. Third, the toolkit is applied to the study of two models of current research interest: A one-dimensional spin chain in a staggered external magnetic field is studied and confinement of the elementary spinon excitations, as predicted by analytical arguments, found numerically using real-time evolution and evaluation of the dynamical structure factor. Additionally, the Hubbard model in two dimensions is studied extensively at various system sizes, geometries, interaction strengths U and filling factors n using up to 30'000 SU(2)-Spin-symmetric states equivalent to approx. 100'000 states in other MPS-DMRG implementations. Hints of a possible phase coexistence in the region 0.85 < n < 0.95 are found at intermediate interaction strengths U = 4 and U = 6 as well as a consistently striped ground state in the region n ≈ 0.875.
The simulation and numerical study of large, strongly correlated quantum systems containing Fermions or using real-time evolution in finite dimensions is still an essentially unsolved problem, primarily due to the exponential growth of the Hilbert state space with system size and the occurrence of the so-called sign problem in Monte Carlo studies. In this area, the use of tensor-network methods, for one-dimensional systems chief among them the density matrix renormalisation group (DMRG) and matrix-product states (MPS), has grown in importance in recent years. This thesis first recapitulates the use of non-abelian symmetries such as SU(2)-Spin in arbitrary tensor networks with an extensive review of the published literature including detailed algorithms and implementation hints. Implementing such symmetries can lead to a considerably more efficient representation of states in the tensor network. This part is intended to be suitable as an implementation-oriented introduction to tensor networks in general and the implementation of non-abelian symmetries in particular. Second, it introduces a series of technical improvements for the MPS methods. These improvements include a faster convergence scheme for MPS-DMRG, a systematic approach to the construction of matrix-product operators and an improved Krylov time evolution method as well as the combination of several well-known techniques into a single tensor network toolkit, SYTEN. The effectiveness of these improvements is demonstrated in numerical examples. Third, the toolkit is applied to the study of two models of current research interest: A one-dimensional spin chain in a staggered external magnetic field is studied and confinement of the elementary spinon excitations, as predicted by analytical arguments, found numerically using real-time evolution and evaluation of the dynamical structure factor. Additionally, the Hubbard model in two dimensions is studied extensively at various system sizes, geometries, interaction strengths U and filling factors n using up to 30'000 SU(2)-Spin-symmetric states equivalent to approx. 100'000 states in other MPS-DMRG implementations. Hints of a possible phase coexistence in the region 0.85 < n < 0.95 are found at intermediate interaction strengths U = 4 and U = 6 as well as a consistently striped ground state in the region n ≈ 0.875.
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.
The research work in this thesis is based on strongly interacting quantum lattice systems. The biggest part of research was conducted using state-of-the-art tensor network simulations. Tensor networks provide efficient and highly accurate representations of quantum states when any simply connected patch of the quantum state is slightly entangled. Matrix Product States (MPS) is a tensor network representation for quantum states which is quasi-exact for one-dimensional systems when entanglement entropy of any bipartition of the state follows ``area-law". Projected Entangled Pair States (PEPS) is the extension of MPS for higher dimensional systems, where entanglement entropy area-law is non-trivial. These formalisms and their relevant ground state optimization techniques, Density Matrix Renormalization Group (DMRG) for MPS and Simple Update (SU) with Tensor Renormalization Group (TRG) for PEPS are thoroughly analysed. Entanglement entropy area-law is usually obeyed by the ground state of local Hamiltonians, while generic highly excited states follow a volume-law (they span a finite part of the Hilbert space). Recently, a class of interacting system was shown to undergo a dynamical phase transition (Many Body Localization) where the entropy of every eigenstate follows an area-law. This transition is achieved when the system is highly disordered and the quantum many body state becomes localized similarly to the free particle Anderson localization. The area-law property makes highly excited eigenstates of Many Body Localized (MBL) systems efficiently represented by the MPS ansatz. We develop a highly optimized algorithm (eDMRG) which goes beyond ground state optimization and successfully target eigenstates in any part of the spectrum. This algorithm is used together with analytical calculations, based on local integrability of MBL systems, to identify the universal behaviour of the entanglement spectrum of highly excited eigenstates in MBL systems. In the second part we study interacting bosons in 2D optical lattices. Strongly interacting bosons are simulated using the Bose-Hubbard (BH) model when the interactions are strictly local and the Extended Bose-Hubbard (EBH) model, when additional dipolar interactions are present. The ground states of BH/EBH Hamiltonians in a hexagonal lattice are studied. The phase spectrum includes various insulating and critical phases which are studied in detail using infinite-PEPS ansatz. Additional results on entanglement entropy scaling with respect to the filling of the lattice are presented. Finally, a strongly interacting Harper-Hofstadter Hamiltonian is realized by combining synthetic fields with the BH model. The homogeneity of the system is then broken by a parabolic trapping potential, similar to the ones used in cold atom experiments. Using time-dependent Gutzwiller ansatz (GA) the expansion dynamics of the cloud in large square lattices are studied. In contrast to the expansion dynamics of the BH model, it is found that the synthetic fields generate a self-trapping effect. Using phenomenology and simulations the dynamics are studied in the hard-core and soft-core limit.
The field of tensor networks, kicked off in 1992 by Steve White's invention of the spectacularly successful density matrix renormalization group (DMRG) algorithm, has exploded in popularity in recent years. Tensor networks are poised to play a role in helping us solve some of the greatest open physics problems of our time, such as understanding the nature of high-temperature superconductivity and illuminating a theory of quantum gravity. DMRG and extensions based on a class of variational states known as tensor network states have been indispensable tools in helping us understand both numerically and theoretically the properties of complicated classical and quantum many-body systems. However, practical challenges to these techniques still remain, and algorithmic developments are needed before tensor network algorithms can be applied to more physics problems. In this thesis we present a variety of recent advancements to tensor network algorithms. First we describe a DMRG-like algorithm for noninteracting fermions. Noninteracting fermions, naturally being gapless and therefore having high levels of entanglement, are actually a challenging setting for standard DMRG algorithms, and we believe this new algorithm can help with tensor network calculations in that setting. Next we explain a new algorithm called the variational uniform matrix product state (VUMPS) algorithm that is a DMRG-like algorithm that works directly in the thermodynamic limit, improving upon currently available MPS-based methods for studying infinite 1D and quasi-1D quantum many-body systems. Finally, we describe a variety of improvements to algorithms for contracting 2D tensor networks, a common problem in tensor network algorithms, for example for studying 2D classical statistical mechanics problems and 2D quantum many-body problems with projected entangled pair states (PEPS). One is a new variant of the corner transfer matrix renormalization group (CTMRG) algorithm of Nishino and Okunishi that improves the numerical stability for contracting asymmetric two-dimensional tensor networks compared to the most commonly used method. Another is the application of the VUMPS algorithm to contracting 2D tensor networks. The last is a new alternative to CTMRG, where the tensors are solved for with eigenvalue equations instead of a power method, which we call the fixed point corner method (FPCM). We present results showing the transfer matrix VUMPS algorithm and FPCM significantly improve upon the convergence time of CTMRG. We expect these algorithms will play an important role in expanding the set of 2D classical and 2D quantum many-body problems that can be addressed with tensor networks.
The simulation of strongly correlated quantum systems on classical computers is an essential tool for investigating and understanding quantum materials.Due to the inherent "curse of dimensionality" of quantum systems, such simulations pose a challenging task, and have spurred the development of several computational approaches and approximations.For quasi one-dimensional lattice systems, the density matrix renormalization group (DMRG) framework (Schollwöck, 2005;White, 1992White, , 1993) has emerged as one of the most successful methods.Modern formulations are based on matrix product states (MPS) and operators (MPO) (Schollwöck, 2011): briefly speaking, a matrix product state is an Ansatz for the quantum many-body wavefunction amendable to efficient numerical manipulation.The favorable approximation properties of matrix product states have been theoretically proven for ground states of gapped Hamiltonians with local interactions (Hastings, 2007;Verstraete & Cirac, 2006), and recent work has extended this to quantum systems at finite temperature (Barthel, 2017).Over the years many variations and generalizations of the DMRG framework have been devised (Barthel, 2013;Stoudenmire & White, 2010;Verstraete, Murg, & Cirac, 2008;Vidal, 2003Vidal, , 2008)), and further developments and improvements of associated algorithms remain an active field of research.
This thesis is focused on many-body localization (MBL) and the development of algorithms using the tensor networks representation for many-body localized systems. Chapter 1 is a detailed discussion on the MBL phenomenon. Chapter 2 gives an introduction to the numerical techniques that are frequently used in studying many-body quantum systems. Chapter 3 studies the shift-and-invert matrix product state method for obtaining highly excited states of MBL systems. Chapter 4 studies an efficient full diagonalization method using Wegner unitary gates. Chapter 5 studies the distribution of the entanglement entropy in the MBL-ergodic transition region. Chapter 6 studies the the influence of pseudo-spin SU (2) symmetry on the MBL phase. Chapter 7 is a side project that studies the influence of the duality-twisted boundary conditions on the critical twist defect chains. The MBL-ergodic transition is a dynamical phase transition which happens at finite energy density for a disordered and isolated many-body interacting system. It is a natural extension of Anderson localization where the localized phase survives the inter-particle interaction. In Chapter 1, we describe its relation to quantum thermalization, typical lattice models that harbor a MBL phase, MBL system’s phenomenological description, the characteristics of the MBL phase, experimental studies of MBL phenomenon in 1D and higher dimensions, the latest developments and the open questions on the nature of the MBL-ergodic transition. In Chapter 2, we provide a short but self-contained discussion of basic numerical methods that are frequently used in studying MBL physics. We mainly focus on the exact diagonalization methods and the matrix product state (MPS) representation, and explain their advantages and disadvantages. In Chapter 3, because obtaining the interior eigenstates of MBL Hamiltonians is a key step to studying the MBL phases and the transitions, we develop algorithms that can capture individual excited states to high fidelity. Using the fact that eigenstates of many-body localized systems have area-law entanglement and can therefore be efficiently represented as a matrix product state (MPS), we designed two algorithms that can generate excited states in the MPS representation, and use them to test the basic properties of MBL in the regime of large one-dimensional random field Heisenberg chains that were previously inaccessible due to the limitations of exact diagonalization (ED). In Chapter 4, we investigate an efficient way to diagonlize a fully many-body localized (FMBL) Hamiltonian. Using the Wegner flow technique for subsystems, we constructively build unitary gates that can greedily transform a local part of the FMBL Hamiltonian with the aim of lowering its average energy variance. Applying the Wegner unitary successively can make the system “flow” towards the diagonalized form. We find that the performance of this method is mainly controlled by the length of each individual Wegner unitary gate, and gets better with increasing length of the gate. We compare the performance of this constructive method with previous methods based on quasi-Newtonian optimization methods, and find that they have similar performance, while our constructive method is far more efficient. We also comment the possibility of constructing a tensor network state that resembles a multi-scale entanglement renormalization ansatz (MERA) using this constructive method. In Chapter 5, using strong subadditivity (SSA) theorem, we develop two order parameters – cut-averaged entanglement entropy (CAEE) and its slope (SCAEE), and show that they can be used to directly identify the volume or area law scaling in single eigenstates. We study the distribution of the SCAEE over disorder realizations. This distribution appears Gaussian at weak disorder, while on the MBL side, the distribution of the SCAEE is peaked at zero slope and has an exponential tail. In the critical regime we find that the distribution of the SCAEE is bimodal both over multiple disorder realizations as well as for single disorder realizations. The variance of the SCAEE distribution in the transition region seems to grow with system size when considered over disorder realizations. Our system sizes are too small to pin down its maximal value, but they are consistent with (among other possibilities) the maximal variance possible which would lead to half the states having zero SCAEE and half having maximal SCAEE. This scenario would lead to an entanglement entropy at the transition which scales as a volume law with half its thermal value. In Chapter 6, we consider the strongly disordered one-dimensional Hubbard model with spin disorder, which, as we argue, cannot be categorized as a conventional MBL phase. Though the spin rotation SU (2) symmetry is broken by the spin disorder, this model still preserves the pseudo-spin SU (2) symmetry under periodic boundary conditions. On the theoretical side, using the pseudo-spin algebra we show that a significant number of the excited states at any disorder strength have logarithmic correction to their von Neumann entanglement entropy, which violates the area-law entanglement for a typical MBL phase. A common feature of this group of excited states is that they usually have far more double occupancies than single occupancies. On the numerical side, we studied the time evolution after quantum quench from two extreme cases of product states – all single occupancies at quarter filling, and all double occupancies at half filling. We find in the former case entanglement entropy behaves as in a typical MBL system, while the latter is clearly delocalized, which suggests the existence of incomplete set of local integrals of motion. With the above evidence, it is convincing that this model provides a playground for studying a non-ergodic, non-MBL phase, and its relation to continuous SU (2) symmetry. Chapter 7 deviates from the common scene of MBL physics. In this chapter, we consider a series of one-dimensional critical twist defect models with various boundary conditions. These twist defects are embedded in the background of a generalization of Kitaev Z 2 toric code model, and can also be understood as a discrete Z k gauge theory in its deconfined phase. We analytically construct the duality transformation of these critical models, and use numerical methods to extract information on the central charges, sound velocities, and conformal dimensions. With the numerical data, we show the relation between the CFT contents of the critical chains with periodic boundary conditions and the critical chains with duality-twisted boundary conditions.
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.