The many-body localised (MBL) to thermal crossover observed in exact diagonalisation studies remains poorly understood as the accessible system sizes are too small to be in an asymptotic scaling regime. We develop a model of the crossover in short 1D chains in which the MBL phase is destabilised by the formation of many-body resonances. The model reproduces several properties of the numerically observed crossover, including an apparent correlation length exponent \nu=1 ν = 1 , exponential growth of the Thouless time with disorder strength, linear drift of the critical disorder strength with system size, scale-free resonances, apparent 1/\omega 1 / ω dependence of disorder-averaged spectral functions, and sub-thermal entanglement entropy of small subsystems. In the crossover, resonances induced by a local perturbation are rare at numerically accessible system sizes L L which are smaller than a \lambda λ . For L \gg \sqrt{\lambda} L ≫ λ (in lattice units), resonances typically overlap, and this model does not describe the asymptotic transition. The model further reproduces controversial numerical observations which Refs. claimed to be inconsistent with MBL. We thus argue that the numerics to date is consistent with a MBL phase in the thermodynamic limit.
The many-body localised (MBL) to thermal crossover observed in exact diagonalisation studies remains poorly understood as the accessible system sizes are too small to be in an asymptotic scaling regime. We develop a model of the crossover in short 1D chains in which the MBL phase is destabilised by the formation of many-body resonances. The model reproduces several properties of the numerically observed crossover, including an apparent correlation length exponent \nu=1 ν = 1 , exponential growth of the Thouless time with disorder strength, linear drift of the critical disorder strength with system size, scale-free resonances, apparent 1/\omega 1 / ω dependence of disorder-averaged spectral functions, and sub-thermal entanglement entropy of small subsystems. In the crossover, resonances induced by a local perturbation are rare at numerically accessible system sizes L L which are smaller than a \lambda λ . For L \gg \sqrt{\lambda} L ≫ λ (in lattice units), resonances typically overlap, and this model does not describe the asymptotic transition. The model further reproduces controversial numerical observations which Refs. claimed to be inconsistent with MBL. We thus argue that the numerics to date is consistent with a MBL phase in the thermodynamic limit.
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.
We numerically explore the many body localization (MBL) transition through the lens of the entanglement spectrum . While a direct transition from localization to thermalization is believed to be obtained in the thermodynamic limit (the exact details of which remain an open problem), in finite system sizes there exists an intermediate 'quantum critical' regime. Previous numerical investigations have explored the crossover from thermalization to criticality, and have used this to place a numerical lower bound on the critical disorder strength for MBL. A careful analysis of the high energy part of the entanglement spectrum (which contains universal information about the critical point) allows us to study the crossover from criticality to MBL, and we find evidence for such a crossover which could allow us to place a numerical upper bound on the critical disorder strength for MBL.
The space of one-dimensional disordered interacting quantum models displaying a many-body localization (MBL) transition seems sufficiently rich to produce critical points with level statistics interpolating continuously between the Poisson statistics of the localized phase and the Wigner–Dyson statistics of the delocalized phase. In this paper, we consider the strong disorder limit of the MBL transition, where the level statistics at the MBL critical point is close to the Poisson statistics. We analyze a one-dimensional quantum spin model, in order to determine the statistical properties of the rare extensive resonances that are needed to destabilize the MBL phase. At criticality, we find that the entanglement entropy can grow with an exponent 0 < α < 1 anywhere between the area law α = 0 and the volume law α = 1 , as a function of the resonances properties, while the entanglement spectrum follows the strong multifractality statistics. In the MBL phase near criticality, we obtain the simple value ν = 1 for the correlation length exponent. Independently of the strong disorder limit, we explain why, for the many-body localization transition concerning individual eigenstates, the correlation length exponent ν is not constrained by the usual Harris inequality ν ≥ 2 / d , so that there is no theoretical inconsistency with the best numerical measure ν = 0 . 8 ( 3 ) obtained by Luitz et al. (2015).
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.
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 present a detailed analysis of the length-and timescales needed to approach the critical region of MBL from the delocalised phase, studying both eigenstates and the time evolution of an initial state.For the eigenstates we show that in the delocalised region there is a single length, which is a function of disorder strength, controlling the finite-size flow.Small systems look localised, and only for larger systems do resonances develop which restore ergodicity in the form of the eigenstate thermalisation hypothesis.For the transport properties, we study the time necessary to transport a single spin across a domain wall, showing how this grows quickly with increasing disorder, and compare it with the Heisenberg time.For a sufficiently large system the Heisenberg time is always larger than the transport time, but for a smaller system this is not necessarily the case.We conclude that the properties of the MBL transition cannot be explored using the system sizes or times available to current numerical and experimental studies.