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.
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.
Since the seminal work of Anderson, localisation has been recognised as a standard mechanism allowing quantum many-body systems to escape ergodicity. This idea acquired even more prominence in the last decade as it has been argued that localisation -- dubbed many-body localisation (MBL) in this context -- can sometimes survive local interactions in the presence of sufficiently strong disorder. A conventional signature of localisation 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 localisation, they have never been derived from first principles in the presence of non-trivial interactions. Here we fill this gap and provide the example of a simple class of quantum many-body systems -- which we dub strongly localised quantum circuits -- that are interacting, localised, and where the spectral statistics can be characterised 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 they become Poissonian.
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.
Recently it has been suggested that many-body localization (MBL) can occur in translationally-invariant systems, and candidate 1D models have been proposed. We find that such models typically exhibit strong finite-size effects due to the presence of two or more vastly different energy scales, including an artificial separation of density of states into bands separated by large energy gaps. By carefully choosing the parameters of the models such that finite-size effects are minimized, we perform extensive numerical studies of small systems using exact diagonalization. Based on the variety of diagnostics, including entanglement properties and the behavior of local observables, we find the systems exhibit thermal (ergodic), rather than MBL-like behavior. Our results suggest that MBL in translationally-invariant systems with two very different energy scales is less robust than perturbative arguments suggest, possibly pointing to the importance of non-perturbative effects which induce delocalization in the thermodynamic limit.
The study of many-body localized (MBL) phases intrinsically links spectral properties with eigenstate characteristics: localized systems exhibit Poisson level statistics and area-law entanglement entropy, while ergodic systems display volume-law entanglement and follow random matrix theory predictions, including level repulsion. Here, we introduce the interacting Anderson Quantum Sun model, which significantly deviates from these conventional expectations. In addition to standard localized and ergodic phases, we identify a regime that exhibits volume-law entanglement coexisting with intermediate spectral statistics. We also identify another nonstandard regime marked by Poisson level statistics, sub-volume entanglement growth, and rare-event-dominated correlations, indicative of emerging ergodic instabilities. These results highlight unconventional routes of ergodicity breaking and offer fresh perspectives on how Anderson localization may be destabilized.
Some interacting disordered many-body systems are unable to thermalize when the quenched disorder becomes larger than a threshold value. Although several properties of nonzero energy density eigenstates (in the middle of the many-body spectrum) exhibit a qualitative change across this many-body localization (MBL) transition, many of the commonly used diagnostics only do so over a broad transition regime. Here we provide evidence that the transition can be located precisely even at modest system sizes by sharply defined changes in the distribution of extremal eigenvalues of the reduced density matrix of subsystems. In particular, our results suggest that $p*={lim}_{{\ensuremath{\lambda}}_{2}\ensuremath{\rightarrow}ln{(2)}^{+}}{P}_{2}({\ensuremath{\lambda}}_{2})$, where ${P}_{2}({\ensuremath{\lambda}}_{2})$ is the probability distribution of the second lowest entanglement eigenvalue ${\ensuremath{\lambda}}_{2}$, behaves as an ``order parameter'' for the MBL phase: p∗>0 in the MBL phase, while p∗=0 in the ergodic phase with thermalization. Thus, in the MBL phase, there is a nonzero probability that a subsystem is entangled with the rest of the system only via the entanglement of one subsystem qubit with degrees of freedom outside the region. In contrast, this probability vanishes in the thermal phase.
We suggest that if a localized phase at nonzero temperature T>0 exists for strongly disordered and weakly interacting electrons, as recently argued, it will also occur when both disorder and interactions are strong and T is very high. We show that in this high-T regime, the localization transition may be studied numerically through exact diagonalization of small systems. We obtain spectra for one-dimensional lattice models of interacting spinless fermions in a random potential. As expected, the spectral statistics of finite-size samples cross over from those of orthogonal random matrices in the diffusive regime at weak random potential to Poisson statistics in the localized regime at strong randomness. However, these data show deviations from simple one-parameter finite-size scaling: the apparent mobility edge ``drifts'' as the system's size is increased. Based on spectral statistics alone, we have thus been unable to make a strong numerical case for the presence of a many-body localized phase at nonzero T.