Abstract A novel method has been devised to compute the Local Integrals ofMotion~(LIOMs) for a one-dimensional many-body localized system. Inthis approach, a class of optimal unitary transformations is deducedin a tensor-network formalism to diagonalize the Hamiltonian of thespecified system. To construct the tensor network, we utilize theeigenstates of the subsystems' Hamiltonian to attain the desiredunitary transformations. Subsequently, we optimize the eigenstatesand acquire appropriate unitary localized operators that willrepresent the LIOMs tensor network. The efficiency of the method wasassessed and found to be both fast and almost accurate. In frameworkof the introduced tensor-network representation, we examine how theentanglement spreads along the considered many-body localized systemand evaluate the outcomes of the approximations employed in thisapproach. The important and interesting result is that in theproposed tensor network approximation, if the length of the blocksis greater than the length of localization, then the entropy growthwill be linear in terms of the logarithmic time. Also, it has beendemonstrated that, the entanglement can be calculated by onlyconsidering two blocks next to each other, if the Hamiltonian hasbeen diagonalized using the unitary transformation made by theprovided tensor-network representation.
Abstract A novel method has been devised to compute the Local Integrals ofMotion~(LIOMs) for a one-dimensional many-body localized system. Inthis approach, a class of optimal unitary transformations is deducedin a tensor-network formalism to diagonalize the Hamiltonian of thespecified system. To construct the tensor network, we utilize theeigenstates of the subsystems' Hamiltonian to attain the desiredunitary transformations. Subsequently, we optimize the eigenstatesand acquire appropriate unitary localized operators that willrepresent the LIOMs tensor network. The efficiency of the method wasassessed and found to be both fast and almost accurate. In frameworkof the introduced tensor-network representation, we examine how theentanglement spreads along the considered many-body localized systemand evaluate the outcomes of the approximations employed in thisapproach. The important and interesting result is that in theproposed tensor network approximation, if the length of the blocksis greater than the length of localization, then the entropy growthwill be linear in terms of the logarithmic time. Also, it has beendemonstrated that, the entanglement can be calculated by onlyconsidering two blocks next to each other, if the Hamiltonian hasbeen diagonalized using the unitary transformation made by theprovided tensor-network representation.
We propose a tensor network encoding the set of all eigenstates of a fully many-body localized system in one dimension. Our construction, conceptually based on the ansatz introduced in Phys. Rev. B 94, 041116(R) (2016), is built from two layers of unitary matrices which act on blocks of \nℓ \n contiguous sites. We argue that this yields an exponential reduction in computational time and memory requirement as compared to all previous approaches for finding a representation of the complete eigenspectrum of large many-body localized systems with a given accuracy. Concretely, we optimize the unitaries by minimizing the magnitude of the commutator of the approximate integrals of motion and the Hamiltonian, which can be done in a local fashion. This further reduces the computational complexity of the tensor networks arising in the minimization process compared to previous work. We test the accuracy of our method by comparing the approximate energy spectrum to exact diagonalization results for the random-field Heisenberg model on 16 sites. We find that the technique is highly accurate deep in the localized regime and maintains a surprising degree of accuracy in predicting certain local quantities even in the vicinity of the predicted dynamical phase transition. To demonstrate the power of our technique, we study a system of 72 sites, and we are able to see clear signatures of the phase transition. Our work opens a new avenue to study properties of the many-body localization transition in large systems.
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.
Subsystems of strongly disordered, interacting quantum systems can fail to thermalize because of the phenomenon of many-body localization (MBL). In this paper, we explore a tensor network description of the eigenspectra of such systems. Specifically, we will argue that the presence of a complete set of local integrals of motion in MBL implies an efficient representation of the entire spectrum of energy eigenstates with a single tensor network, a spectral tensor network. Our results are rigorous for a class of idealized systems related to MBL with integrals of motion of finite support. In one spatial dimension, the spectral tensor network allows for the efficient computation of expectation values of a large class of operators (including local operators and string operators) in individual energy eigenstates and in ensembles.
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 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
We consider the representation of operators in terms of tensor networks and their application to the ground-state approximation and time evolution of systems with long-range interactions. We provide an explicit construction to represent an arbitrary many-body Hamilton operator in terms of a one-dimensional tensor network (i.e., as a matrix product operator). For pairwise interactions, we show that such a representation is always efficient and requires a tensor dimension growing only linearly with the number of particles. For systems obeying certain symmetries or restrictions we find optimal representations with minimal tensor dimension. We discuss the analytic and numerical approximation of operators in terms of low-dimensional tensor operators. We demonstrate applications for time evolution and the ground-state approximation, in particular for long-range interaction with inhomogeneous couplings. The operator representations are also generalized to other geometries such as trees and two-dimensional lattices, where we show how to obtain and use efficient tensor network representations respecting a given geometry.