We study the properties of excited states in one-dimensional many-body localized (MBL) systems using a matrix product state algorithm. First, the method is tested for a large disordered noninteracting system, where for comparison we compute a quasiexact reference solution via a Monte Carlo sampling of the single-particle levels. Thereafter, we present extensive data obtained for large interacting systems of $L\ensuremath{\sim}100$ sites and large bond dimensions $\ensuremath{\chi}\ensuremath{\sim}1700$, which allows us to quantitatively analyze the scaling behavior of the entanglement $S$ in the system. The MBL phase is characterized by a logarithmic growth $S(L)\ensuremath{\sim}log(L)$ over a large scale separating the regimes where volume and area laws hold. We check the validity of the eigenstate thermalization hypothesis. Our results are consistent with the existence of a mobility edge.
We study the properties of excited states in one-dimensional many-body localized (MBL) systems using a matrix product state algorithm. First, the method is tested for a large disordered noninteracting system, where for comparison we compute a quasiexact reference solution via a Monte Carlo sampling of the single-particle levels. Thereafter, we present extensive data obtained for large interacting systems of L\ensuremath∼100 sites and large bond dimensions \ensuremathχ\ensuremath∼1700, which allows us to quantitatively analyze the scaling behavior of the entanglement S in the system. The MBL phase is characterized by a logarithmic growth S(L)\ensuremath∼log(L) over a large scale separating the regimes where volume and area laws hold. We check the validity of the eigenstate thermalization hypothesis. Our results are consistent with the existence of a mobility edge.
A prime characterization of many-body localized (MBL) systems is the entanglement of their eigenstates; in contrast to the typical ergodic phase whose eigenstates are volume law, MBL eigenstates obey an area law. In this paper, we show that a spin-disordered Hubbard model has both a large number of area-law eigenstates as well as a large number of eigenstates whose entanglement scales logarithmically with system size (log-law). This model provides a microscopic Hamiltonian which is neither ergodic nor MBL. We establish these results through a combination of analytic arguments based on the eta-pairing operators combined with a numerical analysis of eigenstates. In addition, we describe and simulate a dynamic time evolution approach starting from product states through which one can separately probe the area-law and log-law eigenstates in this system.
A key property of many-body localized Hamiltonians is the area law entanglement of even highly excited eigenstates. Matrix product states (MPS) can be used to efficiently represent low entanglement (area law) wave functions in one dimension. An important application of MPS is the widely used density matrix renormalization group (DMRG) algorithm for finding ground states of one-dimensional Hamiltonians. Here, we develop two algorithms, the shift-and-invert MPS (SIMPS) and excited state DMRG which find highly excited eigenstates of many-body localized Hamiltonians. Excited state DMRG uses a modified sweeping procedure to identify eigenstates, whereas SIMPS applies the inverse of the shifted Hamiltonian to a MPS multiple times to project out the targeted eigenstate. To demonstrate the power of these methods, we verify the breakdown of the eigenstate thermalization hypothesis in the many-body localized phase of the random field Heisenberg model, show the saturation of entanglement in the many-body localized phase, and generate local excitations.
We consider disordered many-body systems with periodic time-dependent Hamiltonians in one spatial dimension. By studying the properties of the Floquet eigenstates, we identify two distinct phases: (i) a many-body localized (MBL) phase, in which almost all eigenstates have area-law entanglement entropy, and the eigenstate thermalization hypothesis (ETH) is violated, and (ii) a delocalized phase, in which eigenstates have volume-law entanglement and obey the ETH. The MBL phase exhibits logarithmic in time growth of entanglement entropy when the system is initially prepared in a product state, which distinguishes it from the delocalized phase. We propose an effective model of the MBL phase in terms of an extensive number of emergent local integrals of motion, which naturally explains the spectral and dynamical properties of this phase. Numerical data, obtained by exact diagonalization and time-evolving block decimation methods, suggest a direct transition between the two phases.
Many-body localization (MBL) addresses the absence of thermalization in interacting quantum systems, with nonergodic high-energy eigenstates behaving as ground states, only area-law entangled. However, computing highly excited many-body eigenstates using exact methods is very challenging. Instead, we show that one can address high-energy MBL physics using ground-state methods, which are much more amenable to many efficient algorithms. We find that a localized many-body ground state of a given interacting disordered Hamiltonian H0 is a very good approximation for a high-energy eigenstate of a parent Hamiltonian, close to H0 but more disordered. This construction relies on computing the covariance matrix, easily achieved using density-matrix renormalization group for disordered Heisenberg chains up to L=256 sites.
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.
Many-body localization of a disorder interacting boson system in one dimension is studied numerically at the filling factor being one-half, in terms of level statistics, local compressibility, correlation function and entanglement entropies. The von Neumann entanglement entropy is decoupled into a particle number entropy and a configuration entropy. An anomalous volume-law behavior is found for the configuration entanglement entropy to confirm a recent experimental observation [A. Lukin, M. Rispoli, R. Schittko, et al., Science 364, 256 (2019)] for sufficient strong disorder, while the particle number entropy fulfills an area-law corresponding to the total entropy for disordered spin chain. The localization length are extracted from a two-body correlation function for many-body localization states and corresponding time-evolutions states as well. A phase diagrams is established with consisting of an ergodic thermalized region and a many-body-localization region in a parameter space of the disorder strength and the energy density. Two regions are separated by a many-body mobility edge deducted from the standard deviation of the particle-number entanglement entropy, which appears consistent with that based on the localization length. Slow dynamics characterized by a logarithmic time-dependence is explicitly shown for both the particle number entropy and the configuration entropy in an intermediate regime of their time-evolutions, which does not show up in the Anderson localization case, i.e. non-interacting disorder systems.