We introduce a multi-scale diagonalization scheme to study the transition between the many-body localized and the ergodic phase in disordered quantum chains. The scheme assumes a sharp dichotomy between subsystems that behave as localized and resonant spots that obey the Eigenstate Thermalization Hypothesis (ETH). We establish a set of microscopic principles defining the diagonalization scheme, and use them to numerically study the transition in very large systems. To a large extent the results are in agreement with an analytically tractable mean-field analysis of the scheme: We find that at the critical point the system is almost surely localized in the thermodynamic limit, hosting a set of thermal inclusions whose sizes are power-law distributed. On the localized side the {\em typical} localization length is bounded from above. The bound saturates upon approach to criticality, entailing that a finite ergodic inclusion thermalizes a region of diverging diameter. The dominant thermal inclusions have a fractal structure, implying that averaged correlators decay as stretched exponentials throughout the localized phase. Slightly on the ergodic side thermalization occurs through an avalanche instability of the nearly localized bulk, whereby rare, supercritically large ergodic spots eventually thermalize the entire sample. Their size diverges at the transition, while their density vanishes. The non-local, avalanche-like nature of this instability entails a breakdown of single parameter scaling and puts the delocalization transition outside the realm of standard critical phenomena.
We introduce a multi-scale diagonalization scheme to study the transition between the many-body localized and the ergodic phase in disordered quantum chains. The scheme assumes a sharp dichotomy between subsystems that behave as localized and resonant spots that obey the Eigenstate Thermalization Hypothesis (ETH). We establish a set of microscopic principles defining the diagonalization scheme, and use them to numerically study the transition in very large systems. To a large extent the results are in agreement with an analytically tractable mean-field analysis of the scheme: We find that at the critical point the system is almost surely localized in the thermodynamic limit, hosting a set of thermal inclusions whose sizes are power-law distributed. On the localized side the {\em typical} localization length is bounded from above. The bound saturates upon approach to criticality, entailing that a finite ergodic inclusion thermalizes a region of diverging diameter. The dominant thermal inclusions have a fractal structure, implying that averaged correlators decay as stretched exponentials throughout the localized phase. Slightly on the ergodic side thermalization occurs through an avalanche instability of the nearly localized bulk, whereby rare, supercritically large ergodic spots eventually thermalize the entire sample. Their size diverges at the transition, while their density vanishes. The non-local, avalanche-like nature of this instability entails a breakdown of single parameter scaling and puts the delocalization transition outside the realm of standard critical phenomena.
We propose a multiscale diagonalization scheme to study disordered one-dimensional chains, in particular, the transition between many-body localization (MBL) and the ergodic phase, expected to be governed by resonant spots. Our scheme focuses on the dichotomy of MBL versus validity of the eigenstate thermalization hypothesis. We show that a few natural assumptions imply that the system is localized with probability one at criticality. On the ergodic side, delocalization is induced by a quantum avalanche seeded by large ergodic spots, whose size diverges at the transition. On the MBL side, the typical localization length tends to the inverse of the maximal entropy density at the transition, but there is a divergent length scale related to the response to an inclusion of large ergodic spots. A mean-field approximation analytically illustrates these results and predicts a power-law distribution for thermal inclusions at criticality.
We propose a scaling theory for the many-body localization (MBL) phase transition in one dimension, building on the idea that it proceeds via a ``quantum avalanche.'' We argue that the critical properties can be captured at a coarse-grained level by a Kosterlitz-Thouless (KT) renormalization group (RG) flow. On phenomenological grounds, we identify the scaling variables as the density of thermal regions and the length scale that controls the decay of typical matrix elements. Within this KT picture, the MBL phase is a line of fixed points that terminates at the delocalization transition. We discuss two possible scenarios distinguished by the distribution of rare, fractal thermal inclusions within the MBL phase. In the first scenario, these regions have a stretched exponential distribution in the MBL phase. In the second scenario, the near-critical MBL phase hosts rare thermal regions that are power-law-distributed in size. This points to the existence of a second transition within the MBL phase, at which these power laws change to the stretched exponential form expected at strong disorder. We numerically simulate two different phenomenological RGs previously proposed to describe the MBL transition. Both RGs display a universal power-law length distribution of thermal regions at the transition with a critical exponent ${\ensuremath{\alpha}}_{c}=2$, and continuously varying exponents in the MBL phase consistent with the KT picture.
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.
The random field S=21 Heisenberg chain exhibits a dynamical many body localization transition at a critical disorder strength, which depends on the energy density. At weak disorder, the eigenstate thermalization hypothesis (ETH) is fulfilled on average, making local observables smooth functions of energy, whose eigenstate-to-eigenstate fluctuations decrease exponentially with system size. We demonstrate the validity of ETH in the thermal phase as well as its breakdown in the localized phase and show that rare states exist which do not strictly follow ETH, becoming more frequent closer to the transition. Similarly, the probability distribution of the entanglement entropy at intermediate disorder develops long tails all the way down to zero entanglement. We propose that these low entanglement tails stem from localized regions at the subsystem boundaries which were recently discussed as a possible mechanism for subdiffusive transport in the ergodic phase.
While the prethermal regime of random many-body localized (MBL) systems is dominated by accidental many-body resonances, another class of resonances, originating from the underlying potential structure, is expected in large-size deterministic systems. It is known that this class of resonances can lead to single-particle critical phases that are neither localized nor extended, but the consequences in interacting systems remain unclear. In this work, we construct an asymptotically solvable model of a one-dimensional nearest-neighbor interacting spin chain, whose spatial structure induces a hierarchy of mirror-like many-body resonances. We derive two phases in the thermodynamic limit, characterized by the satisfaction and violation of a version of the weak eigenstate thermalization hypothesis (ETH). While these two phases are similar to the usual MBL and ETH phases, there exist rare eigenstates that behave like the opposite phase, interpreted as many-body scars and inverted scars. Surprisingly, the two phases can be separated by a finite-temperature phase transition, corresponding to a thermodynamic many-body mobility edge, which was often believed to be impossible. Our results also suggest the existence of delocalized rare regions in an otherwise-localized interacting Aubry-André model, even if there are no low-disorder regions like those in random systems. This challenges the common belief that there is no avalanche instability in quasiperiodic MBL.
Closed, interacting, quantum systems have the potential to transition to a many-body localized (MBL) phase under the presence of sufficiently strong disorder, hence breaking ergodicity and failing to thermalize. In this work we study the distribution of correlations throughout the ergodic-MBL phase diagram. We find the typical correlations in the MBL phase decay as a stretched exponential with range r eventually crossing over to an exponential decay deep in the MBL phase. At the transition, the stretched exponential goes as e−Ar, a decay that is reminiscent of the random singlet phase. While the standard deviation of the log(QMI) has a range dependence, the log(QMI) converges to a range-invariant distribution on all other moments (i.e., the skewness and higher) at the transition. The universal nature of these distributions provides distinct phenomenology of the transition different from both the ergodic and MBL phenomenologies. In addition to the typical correlations, we study the extreme correlations in the system, finding that the probability of strong long-range correlations is maximal at the transition, suggesting the proliferation of resonances there. Finally, we analyze the probability that a single bit of information is shared across two halves of a system, finding that this probability is non-zero deep in the MBL phase but vanishes at moderate disorder well above the transition.