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.
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.
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.
"Under certain conditions, an interacting system can defy the concept of thermalization, a keystone in our understanding of physical processes. Many-body localization (MBL) is a phase of matter in which thermalization does not apply and ergodicity is broken. This striking behavior, which can appears on closed, interacting, quantum systems subject to strong disorder, has been the focus of a large body of theoretical, numerical, and experimental work in recent years. In this thesis we numerically study several aspects of MBL and the ergodic-MBL transition. Chapter 1 introduces the concept of thermalization in quantum systems, which relies on the idea of the eigenstate thermalization hypothesis (ETH). I then present localization as a phenomenon that breaks ETH, both in its single-particle as in its many-body versions. I discuss some of the main aspects that are known about MBL, as well as some of its open questions, some of which we tackle in later chapters. Chapter 2 discusses the main numerical methods used in this thesis for the study of MBL. I provide both theoretical background as well as discuss some more practical matters, such as their advantages and disadvantages, or details on the software developed/used in my studies. In Chapter 3 I present our work on MBL from the point of view of single-particle orbitals. In this work we access a complete set of approximate integrals of motion of a one-dimensional system by computing the one-particle orbitals (OPO) of highly-excited MBL energy eigenstates, which are obtained through the shift-and-invert matrix product state (SIMPS) algorithm. We then study the properties of the OPOs over large systems, up to L = 64. We find that the OPOs drawn from eigenstates at different energy densities have high overlap and their occupations are correlated with the energy of the eigenstates. Moreover, the standard deviation of the inverse participation ratio of these orbitals is maximal at the nose of the mobility edge. Also, the OPOs decay exponentially in real space, with a correlation length that increases at low disorder. In addition, we find that the probability distribution of the strength of the large-range coupling constants of the number operators generated by the OPOs approach a log-uniform distribution at strong disorder. In Chapter 4 I present our work on the hybridization of eigenstates in either the MBL or the ergodic phase, as well as at the transition. We do so by adiabatially evolving highly excited eigenstates of the Hamiltonian and measuring their their hybridization process with other eigenstates of the system. This hybridization, which dresses the eigenstate and has the potential of bringing it out of the MBL phase through the transition, is a consequence of the ""collisions"" of eigenstates in energy, which avoid level crossings every time their energy gap is small. The hybridization of eigenstates with each other involves only local regions of the system in the MBL phase, ignores locality in the ergodic phase, and is range-independent at the transition. This range independence suggests the proliferation of long-range resonances at the transition, as well as the divergence of a localization length. In Chapter 5 I present our studies on the typical and extreme (atypically strong) correlations across a one-dimensional system in the ergodic-MBL phase diagram. While typical correlations decay exponentially with range in the MBL phase, in the ergodic phase they are constant and independent of the range. Surprisingly, we identify a moderate region of the phase in which typical correlations decay as a stretched exponential with range r, and in particular as exp[-r^(1/2)] at the transition, a decay that is reminiscent of the random singlet phase. Moreover, at the transition the distribution of the logarithm of the correlations show vanishing even excess moments and non-zero range-invariant odd excess moments. This distinct behavior at the transition is in contrast with ergodic and MBL phenomenologies. In addition, we study the extreme correlations in the system. Our results suggest that strong long-range correlations proliferate at the transition, in contrast with a decay with range in the MBL phase and the lack of strong long-range correlations in the ergodic phase. Finally, we analyze the probability that a single bit of information is shared across two halves of a system, which has been proposed as a robust order parameter in the ergodic-MBL phase diagram. We find that this probability is non-zero deep in the MBL phase, but vanishes at moderate disorder, well above the transition, thus not providing a proper order parameter. Chapter 6 I summarize the work presented with an emphasis on context and perspective."
The phenomenon of many-body localization generalizes Anderson localization to interacting systems. Understanding this conceptually novel phenomenon requires a study of many-body eigenstates at finite energy densities. This is a very challenging task since the most efficient numerical methods such as, e.g., the density matrix renormalization group method, can only access the ground state and low lying excitations. In this work, the authors introduce a unitary tensor network based variational method that approximately finds all many-body eigenstates of fully localized Hamiltonians and scales polynomially with system size. The usefulness of their approach is demonstrated by considering the Heisenberg chain in a strongly disordered magnetic field.
Since the seminal work of Anderson, localization has been recognized as a standard mechanism allowing quantum many-body systems to escape ergodicity. This idea has acquired even more prominence in the last decade as it has been argued that localization---dubbed many-body localization (MBL) in this context---can sometimes survive local interactions in the presence of sufficiently strong disorder. A conventional signature of localization is in the statistical properties of the spectrum---spectral statistics---which differ qualitatively from those in the ergodic phase. Although features of the spectral statistics are routinely used as numerical diagnostics for localization, their derivation from first principles in the presence of nontrivial interactions has been lacking. Here we provide the example of a simple class of quantum many-body systems---which we dub strongly localized quantum circuits---that are interacting and localized and where the spectral statistics can be characterized exactly. Furthermore, we show that these systems exhibit a cascade of three different regimes for spectral correlations depending on the energy scale: at small, intermediate, and large scales they behave as disconnected patches of three decreasing sizes. We argue that these features appear in generic MBL systems, with the difference that only at the smallest scale do they become Poissonian.