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 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 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.
Isolated quantum systems at strong disorder can display many-body localization (MBL), a remarkable phenomena characterized by an absence of conduction even at finite temperatures. As the ratio of interactions to disorder is increased, one expects that an MBL phase will eventually undergo a dynamical phase transition to a delocalized phase. Here we constrain the nature of such a transition by exploiting the strong subadditivity of entanglement entropy, as applied to the many-body eigenstates close to the transition in general dimensions. In particular, we show that at a putative continuous transition between an MBL and an ergodic delocalized phase, the critical eigenstates are necessarily thermal, and therefore, the critical entanglement entropy equals the thermal entropy. We also explore a qualitatively different continuous localization-delocalization transition, where the delocalized phase is non-ergodic whose volume law entanglement entropy tends to zero as the transition is approached.
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.
The transition between many-body localized states and the delocalized thermal states is an eigen-state phase transition at finite energy density outside the scope of conventional quantum statistical mechanics. In this work we investigate the properties of the transition by studying the behavior of the entanglement entropy of a subsystem of size LA in a system of size L>LA near the critical regime of the many-body localization transition. The many-body eigenstates are obtained by exact diagonalization of a disordered quantum spin chain under twisted boundary conditions to reduce the finite-size effect. We present a scaling theory based on the assumption that the transition is continuous and use the subsystem size LA/ξ as the scaling variable, where ξ is the correlation length. We show that this scaling theory provides an effective description of the critical behavior and that the entanglement entropy follows the thermal volume law at the transition point. We extract the critical exponent governing the divergence of ξ upon approaching the transition point. We also study the participation entropy in the spin-basis of the domain wall excitations and show that the transition point and the critical exponent agree with those obtained from finite size scaling of the entanglement entropy. Our findings suggest that the many-body localization transition in this model is continuous and describable as a localization transition in the many-body configuration space.
We investigate the stability of the many-body localized phase against quantum avalanche instabilities in a one-dimensional Heisenberg spin chain with long-range power-law interactions (V∝r−α). By combining exact diagonalization of static properties with Lindblad master equation simulations of open-system dynamics, we systematically map the interplay between interaction range and disorder strength. Our finite-size scaling analysis of entanglement entropy identifies a critical interaction exponent αc≈2, which separates a fragile regime, characterized by an exponentially diverging critical disorder, from a robust short-range regime. To rigorously test the system's resistance to avalanches, we couple the boundary to an infinite-temperature bath and track the propagation of the thermalization front into the localized bulk. We find that the characteristic thermalization time follows a unified scaling law, Trth∼exp[κ(α)LW] (herein, L is the system size, and W is the disorder intensity), which diverges exponentially with the product of system size and disorder strength. This suppression enables the derivation of a quantitative stability criterion, Wstab(α), representing the minimum critical disorder strength required to maintain avalanche stability. Our results confirm that the MBL phase remains asymptotically stable in the thermodynamic limit when disorder exceeds an interaction-dependent threshold, bridging theoretical debates on long-range MBL and providing a roadmap for observing these dynamics in experimental platforms such as Rydberg atom arrays.
An interacting quantum system can transition from an ergodic to a many-body localized (MBL) phase under the presence of sufficiently large disorder. Both phases are radically different in their dynamical properties, which are characterized by highly excited eigenstates of the Hamiltonian. Each eigenstate can be characterized by the set of quantum numbers over the set of (local, in the MBL phase) integrals of motion of the system. In this work we study the evolution of the eigenstates of the disordered Heisenberg model as the disorder strength, W, is varied adiabatically. We focus on the probability that two `colliding' eigenstates hybridize as a function of both the range R at which they differ as well as the strength of their hybridization. We find, in the MBL phase, that the probability of a colliding eigenstate hybridizing strongly at range R decays as Pr(R)∝exp[−R/η], with a length scale η(W)=1/(Blog(W/Wc)) which diverges at the critical disorder strength Wc. This leads to range-invariance at the transition, suggesting the formation of resonating cat states at all ranges. This range invariance does not survive to the ergodic phase, where hybridization is exponentially more likely at large range, a fact that can be understood with simple combinatorial arguments. In fact, compensating for these combinatorial effects allows us to define an additional correlation length ξ in the MBL phase which is in excellent agreement with previous works and which takes the critical value 1/log(2) at the transition, found in previous works to destabilize the MBL phase. Finally, we show that deep in the MBL phase hybridization is dominated by two-level collisions of eigenstates close in energy.