Rare regions with weak disorder (Griffiths regions) have the potential to spoil localization. We describe a non-perturbative construction of local integrals of motion (LIOMs) for a weakly interacting spin chain in one dimension, under a physically reasonable assumption on the statistics of eigenvalues. We discuss ideas about the situation in higher dimensions, where one can no longer ensure that interactions involving the Griffiths regions are much smaller than the typical energy-level spacing for such regions. We argue that ergodicity is restored in dimension d >1, although equilibration should be extremely slow, similar to the dynamics of glasses. This article is part of the themed issue ‘Breakdown of ergodicity in quantum systems: from solids to synthetic matter’.
Rare regions with weak disorder (Griffiths regions) have the potential to spoil localization. We describe a non-perturbative construction of local integrals of motion (LIOMs) for a weakly interacting spin chain in one dimension, under a physically reasonable assumption on the statistics of eigenvalues. We discuss ideas about the situation in higher dimensions, where one can no longer ensure that interactions involving the Griffiths regions are much smaller than the typical energy-level spacing for such regions. We argue that ergodicity is restored in dimension d >1, although equilibration should be extremely slow, similar to the dynamics of glasses. This article is part of the themed issue ‘Breakdown of ergodicity in quantum systems: from solids to synthetic matter’.
Many-body localization in disordered systems in one spatial dimension is typically understood in terms of the existence of an extensive number of (quasi) local integrals of motion (LIOMs) which are thought to decay exponentially with distance and interact only weakly with one another. By contrast, little is known about the form of the integrals of motion in disorder-free systems which exhibit localization. Here, we explicitly compute the LIOMs for disorder-free localized systems, focusing on the case of a linearly increasing potential. We show that while in the absence of interactions, the LIOMs decay faster than exponentially, the addition of interactions leads to the formation of a slow-decaying plateau at short distances. We study how the localization properties of the LIOMs depend on the linear slope, finding that there is a significant finite-size dependence, and present evidence that adding a weak harmonic potential does not result in typical many-body localization phenomenology. By contrast, the addition of disorder has a qualitatively different effect, dramatically modifying the properties of the LIOMS.
Many-body localization plays an increasing role in condensed matter theory, both because it challenges the fundaments of statistical physics, and because it allows us to engineer several new, exotic, stable phases of matter. In this paper, the authors address the issue of the stability of a many-body localized material in contact with an ergodic grain, i.e., an imperfect bath made of a few interacting degrees of freedom. Thanks to detailed microscopic analysis and numerics, they conclude that such an ergodic grain eventually destabilizes the localized phase in the following cases: if the spatial dimension is higher than one, or if the spatial dimension is one but the localization length of the localized material is larger than a fixed threshold value. In realistic materials, these ergodic grains are always present as Griffiths regions where the disorder is anomalously small, and hence, the authors conclude that the localized phase in such materials is unstable, strictly speaking. Transport and thermalization are however exponentially suppressed in the distance between ergodic grains.
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.
It is widely believed that many-body localization in one dimension is fragile and can be easily destroyed by thermal inclusions; however, there are still many open questions regarding the stability of the localized phase and under what conditions it breaks down. Here I construct models with dilute disorder, which interpolate between translationally invariant and fully random models, in order to study the breakdown of localization. This opens up the possibility to controllably increase the density of thermal regions and examine the breakdown of localization as this density is increased. At strong disorder, the numerical results are consistent with commonly used diagnostics for localization even when the concentration of thermalizing regions is high. At moderate disorder, I present evidence for slow dynamics and subdiffusive transport across a large region of the phase diagram, suggestive of a ``bad metal'' phase. This suggests that dilute disorder may be a useful effective model for studying Griffiths effects in many-body localization, and perhaps also in a wider class of disordered systems.
Two of the most active areas in quantum many-particle dynamics involve systems with an unusually large number of conservation laws. Many-body-localized systems generalize ideas of Anderson localization by disorder to interacting systems. While localization still exists with interactions and inhibits thermalization, the interactions between conserved quantities lead to some dramatic differences from the Anderson case. Quantum integrable models such as the XXZ spin chain or Bose gas with delta-function interactions also have infinite sets of conservation laws, again leading to modifications of conventional thermalization. A practical way to treat the hydrodynamic evolution from local equilibrium to global equilibrium in such models is discussed. This paper expands upon a presentation at a discussion meeting of the Royal Society on 7 February 2017. The work described was carried out with a number of collaborators, including Jens Bardarson, Vir Bulchandani, Roni Ilan, Christoph Karrasch, Siddharth Parameswaran, Frank Pollmann and Romain Vasseur.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.
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.