Recent numerical works have revealed the instability of many-body localized (MBL) phase in disordered quantum many-body systems with finite system sizes and over finite timescales. This instability originates from the rare Griffiths regions in the thermodynamic limit, which rapidly thermalize and induce an avalanche mechanism that destabilizes surrounding MBL regions. Here, we consider the [Formula: see text]-preserving interacting Ising Majorana chain, whose phase diagram is enriched by the presence of an intermediate ergodic phase separating two distinct MBL phases characterized by different long-range orders. We investigate the dynamic characteristics of the model when coupled to an infinite bath under perturbation, and through the scaling behavior of the slowest thermalization rate, we find how the critical disorder strength in finite-size systems is affected by the avalanche mechanism. We also implement the embedded inclusion model and use the time evolution of mutual information between each spin and the artificial Griffiths region to probe the diffusion of the thermal bubble. We observe that in finite-sized systems, the critical disorder strength gradually drifts away from the center. Our results show that both the MBL paramagnetic phase and MBL spin-glass phase are unstable at finite sizes.
Quantum avalanches in Z2-preserving interacting Ising Majorana chain
Recent numerical works have revealed the instability of many-body localized (MBL) phase in disordered quantum many-body systems with finite system sizes and over finite timescales. This instability originates from the rare Griffiths regions in the thermodynamic limit, which rapidly thermalize and induce an avalanche mechanism that destabilizes surrounding MBL regions. Here, we consider the [Formula: see text]-preserving interacting Ising Majorana chain, whose phase diagram is enriched by the presence of an intermediate ergodic phase separating two distinct MBL phases characterized by different long-range orders. We investigate the dynamic characteristics of the model when coupled to an infinite bath under perturbation, and through the scaling behavior of the slowest thermalization rate, we find how the critical disorder strength in finite-size systems is affected by the avalanche mechanism. We also implement the embedded inclusion model and use the time evolution of mutual information between each spin and the artificial Griffiths region to probe the diffusion of the thermal bubble. We observe that in finite-sized systems, the critical disorder strength gradually drifts away from the center. Our results show that both the MBL paramagnetic phase and MBL spin-glass phase are unstable at finite sizes.
Many-body localization is a profound phase of matter affecting the entire spectrum which emer ges in the presence of disorder in interacting many-body systems. Recently, the stability of many-body localization has been challenged by the avalanche mechanism, in which a small thermal region can spread, destabilizing localization and leading to global thermalization of the system. A key unresolved question is the critical competition between the thermal region's influence and the disorder strength required to trigger such an avalanche. Here, we numerically investigate many-body localization stability in an isolated Heisenberg spin chain of size L subjected to a disordered magnetic field. By embedding a tunable thermal region of size P, we analyze the system's behavior in both static and dynamical regimes using entanglement entropy and the gap ratio. Our study yields two main findings. Firstly, for strong disorder, the avalanche only occurs if the thermal region scales with system size, specifically when P/L exceeds a threshold value. Secondly, at strong disorder, we identify an intermediate phase between many-body localization and ergodic behavior as P increases. This intermediate phase leaves its fingerprint in both static and dynamic properties of the system and tends to vanish in the thermodynamic limit. Although our simulations are restricted to finite system sizes, the analysis suggests that these results hold in the thermodynamic limit for isolated many-body systems.
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.
We investigate the robustness of the many-body localized (MBL) phase to the quantum-avalanche instability by studying the dynamics of a localized spin chain coupled to a $T=\ensuremath{\infty}$ thermal bath through its leftmost site. By analyzing local magnetizations we estimate the size of the thermalized sector of the chain and find that it increases logarithmically slowly in time. This logarithmically slow propagation of the thermalization front allows us to lower-bound the slowest thermalization time, and find a broad parameter range where it scales fast enough with the system size that MBL is robust against thermalization induced by avalanches. The further finding that the imbalance---a global quantity measuring localization---thermalizes over a timescale that is exponential both in disorder strength and system size is in agreement with these results.
We numerically study quantum avalanches in one-dimensional disordered spin systems by attaching two XXZ spin chains. One chain has low disorder representing a rare Griffith's region, or thermal inclusion, and the second has larger disorder, i.e., disorder larger than the observed finite-size crossover. Comparing dynamics of this system to identical systems with uniformly large disorder, we find evidence for exponentially slow thermalization (in disorder) within the many-body localized regime when the rare region is present. We observe a decay of the spin imbalance in the bulk of the large disorder region that persists to long times $(\ensuremath{\sim}{10}^{4})$ and find a universal behavior of the spectral function.
Abstract The connection between entanglement dynamics and non-equilibrium statistics in isolated many-body quantum systems has been established both theoretically and experimentally. Many-body localization (MBL), a phenomenon where interacting particles in disordered (i.e. random) chains fail to thermalize, exemplifies this connection. However, the systematic proof of critical phenomena such as MBL remains challenging due to the lack of robust methods for analyzing many-body entanglement dynamics. In this paper, we identify MBL through quasi-periodic dynamics in the entanglement evolution of subsystems in a disordered Heisenberg chain. This new form of characterizing MBL, through stable quasi-periodic dynamics of entanglement-where ‘stable’ means they persist in the thermodynamic limit—concretely distinguishes between two competing scenarios: fully localized behavior of subsystems or slowly, exponentially slow in disorder, thermalizing subsystems—a heated controversy in the literature. Utilizing perturbation theory, we derive the entanglement dynamics of single spins through an infinite perturbative series, while also modeling rare Griffiths regions (locally thermal inclusions). Our results prove that in regimes of sufficiently strong disorder, the entanglement evolution of individual subsystems remains quasi-periodic in the thermodynamic limit, thereby providing concrete evidence for the stability of MBL dynamics in disordered Heisenberg chains. This behavior contrasts with the widely reported logarithmic growth of subsystem entanglement in the MBL phase. We show that the logarithmic growth observed in prior studies arises from statistical ensemble averaging, which is prohibited due to the intrinsic non-ergodic dynamics characteristic of MBL systems, rooted in their quasi-periodic features.
Strong disorder in many-body quantum systems leads to many-body localization (MBL), manifested as long-time memory of the initial state. Spatial regions of anomalously weak disorder, arising due to fluctuations in the disordered potential, may act as ergodic inclusions that seed quantum avalanches and destabilize the MBL. To investigate this mechanism in a controlled setting, the authors plant ergodic inclusions by engineering regions of a weak disorder and compare the system's time evolution with predictions of the avalanche theory.