We study many‐body localization (MBL) from the perspective of integrals of motion (IOMs). MBL can be understood phenomenologically through the existence of macroscopically many localized IOMs. We develop a systematic procedure based on IOM to calculate many‐body quantities. Displacement transformations made clear that any operator can be expanded in 1‐,2‐ ... n ‐particles terms. We use this property to develop a systematic procedure to approximately calculate IOMs and many‐body quantities. We characterize the decay with distance of the IOM's and their interactions through effective localization lengths. For all values of disorder the typical IOMs are localized, suggesting the importance of rare fluctuations in understanding the MBL‐to‐ergodic transition. image
We study many‐body localization (MBL) from the perspective of integrals of motion (IOMs). MBL can be understood phenomenologically through the existence of macroscopically many localized IOMs. We develop a systematic procedure based on IOM to calculate many‐body quantities. Displacement transformations made clear that any operator can be expanded in 1‐,2‐ ... n ‐particles terms. We use this property to develop a systematic procedure to approximately calculate IOMs and many‐body quantities. We characterize the decay with distance of the IOM's and their interactions through effective localization lengths. For all values of disorder the typical IOMs are localized, suggesting the importance of rare fluctuations in understanding the MBL‐to‐ergodic transition. image
Many-body localization provides a generic mechanism of ergodicity breaking in quantum systems. In contrast to conventional ergodic systems, many-body-localized (MBL) systems are characterized by extensively many local integrals of motion (LIOM), which underlie the absence of transport and thermalization in these systems. Here we report a physically motivated construction of local integrals of motion in the MBL phase. We show that any local operator (e.g., a local particle number or a spin-flip operator), evolved with the system's Hamiltonian and averaged over time, becomes a LIOM in the MBL phase. Such operators have a clear physical meaning, describing the response of the MBL system to a local perturbation. In particular, when a local operator represents a density of some globally conserved quantity, the corresponding LIOM describes how this conserved quantity propagates through the MBL phase. Being uniquely defined and experimentally measurable, these LIOMs provide a natural tool for characterizing the properties of the MBL phase, in both experiments and numerical simulations. We demonstrate the latter by numerically constructing an extensive set of LIOMs in the MBL phase of a disordered spin-chain model. We show that the resulting LIOMs are quasilocal and use their decay to extract the localization length and establish the location of the transition between the MBL and ergodic phases.
Many interesting properties of many-body localization (MBL) can be understood using a phenomenological model, involving local integrals of motion (LIOMs). Usually, LIOMs can only be computed approximately. Here, the authors explicitly compute a complete set of LIOMs nonperturbatively by maximizing their overlap with physical spin operators. The scheme is immune to accidental resonances and can be applied even at the phase transition point, providing a powerful tool to study microscopic features of MBL as characterized by various localization lengths.
Many-body localization (MBL) is a novel prototype of ergodicity breaking due to the emergence of local integrals of motion (LIOMs) in a disordered interacting quantum system. To better understand the role played by the existence of such LIOMs, we explore and study some of their structural properties across the MBL transition. We first consider a one-dimensional XXZ spin chain in a disordered magnetic field and introduce and implement a nonperturbative, fast, and accurate method of constructing LIOMs. In contrast to already existing methods, our scheme allows obtaining LIOMs not only in the deep MBL phase but, rather, near the transition point too. Then, we take the matrix representation of LIOM operators as an adjacency matrix of a directed graph whose elements describe the connectivity of ordered eigenbasis in the Hilbert space. Our cluster-size analysis for this graph shows that the MBL transition coincides with a percolation transition in the Hilbert space. By performing finite-size scaling, we compare the critical disorder and correlation exponent $\ensuremath{\nu}$ both in the presence and absence of interactions. Finally, we also discuss how the distribution of diagonal elements of LIOM operators in a typical cluster signals the transition.
"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."
Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.
Conservation laws are of paramount importance in our understanding of classical and quantum dynamics. Here, we present a general framework for constructing exact quantum integrals of motion with the desired locality and quantum numbers, which will be illustrated for the case of many-body-localization (MBL). The latter has been understood theoretically in terms of the existence of an extensive number of local integrals of motion (LIOMs). Using our approach, we show that one can express the task of finding LIOMs as an optimization problem. For some specifications, this problem surprisingly connects to the question of finding classical ground states of spin-glass models. Our work unifies previous results obtained in the MBL context and reveals intriguing connections between many-body localization, spin-glass physics, and constrained optimization problems.