Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2303.02748 · Physical review. B./Physical review. B

Many-body delocalization from embedded thermal inclusion

J. Clayton Peacock, Dries Sels

Published Jul 18, 2023• 20 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.108.l020201

Related papers

Quantum avalanches in $$\mathbb {Z}_2$$-preserving interacting Ising Majorana chain

L. M. Zhang, Kai Xu, Heng Fan

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.

arXiv:2501.06064•Scientific Reports•Dec 19, 2025•1 citation

Toy model for anomalous transport and Griffiths effects near the many-body localization transition

Marco Schirò, Marco Tarzia

We introduce and study a toy model for anomalous transport and Griffiths effects in one-dimensional (1D) quantum disordered isolated systems near the many-body localization (MBL) transition. The model is constituted by a collection of 1D tight-binding chains with on-site random energies, locally coupled to a weak perturbation, corresponding to a random matrix which belongs to the Gaussian orthogonal ensemble (GOE) and mimics the effect of thermal inclusions due to delocalizing interactions by providing a local broadening of the Poisson spectrum. While in the absence of such a coupling the model is localized as expected for the one-dimensional Anderson model, increasing the coupling with the GOE perturbation we find a transition to a conducting phase driven by the proliferation of quantum avalanches, which does not fit the standard paradigm of the Anderson transition. In particular an intermediate Griffiths region emerges, where exponentially distributed insulating segments coexist with a few rare resonances. Typical correlations decay exponentially fast, while average correlations decay as stretched exponentials and diverge with the length of the chain, indicating that the conducting inclusions have a fractal structure and that the localization length is broadly distributed at the critical point. This behavior is consistent with a Kosterlitz-Thouless-like criticality of the MBL transition. Transport and relaxation are dominated by rare resonances and rare strong insulating regions, and show anomalous behaviors strikingly similar to those observed in recent simulations and experiments in the bad metal delocalized phase preceding MBL. In particular, we find subdiffusive transport and power-law decay of the return probability at large times, with exponents that gradually change as one moves across the intermediate region. Concomitantly, the ac conductivity vanishes near zero frequency with an anomalous power law.

arXiv:1909.07160•Physical review. B./Physical review. B•Jan 17, 2020•15 citations

Thermal avalanches in isolated many-body localized systems

Muhammad Sajid, Rozhin Yousefjani, Abolfazl Bayat

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.

arXiv:2506.04834•Physical review. B./Physical review. B•Oct 1, 2025•2 citations

Griffiths effects and slow dynamics in nearly many-body localized systems

Sarang Gopalakrishnan, Kartiek Agarwal, Eugene Demler +2 more

The low-frequency response of systems near a many-body localization transition can be dominated by rare regions that are locally critical or ``in the other phase.'' It is known that in one dimension, these rare regions can cause the dc conductivity and diffusion constant to vanish even inside the delocalized thermal phase. Here, we present a general analysis of such Griffiths effects in the thermal phase near the many-body localization transition: we consider both one-dimensional and higher-dimensional systems, subject to quenched randomness, and discuss both linear response (including the frequency- and wave-vector-dependent conductivity) and more general dynamics. In all the regimes we consider, we identify observables that are dominated by rare-region effects. In some cases (one-dimensional systems and Floquet systems with no extensive conserved quantities), essentially all long-time local observables are dominated by rare-region effects; in others, generic observables are instead dominated by hydrodynamic long-time tails throughout the thermal phase, and one must look at specific probes, such as spin echo, to see Griffiths behavior.

arXiv:1511.06389•Physical review. B./Physical review. B•Apr 11, 2016•150 citations

Catching thermal avalanches in the disordered XXZ model

Tomasz Szołdra, Piotr Sierant, Maciej Lewenstein +1 more

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.

arXiv:2402.01362•Physical review. B./Physical review. B•Apr 8, 2024•21 citations

Bath-induced delocalization in interacting disordered spin chains

Dries Sels

Over a decade of work has culminated in the consensus that one-dimensional systems, subject to sufficiently large disorder, fail to thermalize and possess an extensive set of local integrals of motion. In this Letter, I will provide numerical evidence for the contrary. In particular, this work studies the dynamics of disordered spin chains which are weakly coupled to a Markovian bath. Within this approach, the critical disorder for stability to quantum avalanches exceeds ${W}^{*}\ensuremath{\gtrsim}20$ in the random field Heisenberg chain. In stark contrast to the Anderson insulator, the avalanche threshold drifts considerably with system size, with no evidence of saturation in the studied regime.

arXiv:2108.10796•Physical review. B./Physical review. B•Jul 27, 2022•140 citations
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∼104)(\ensuremath{\sim}{10}^{4})(\ensuremath∼104) and find a universal behavior of the spectral function.
Quantum avalanches in Z2\mathbb {Z}_2Z2​-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.
We introduce and study a toy model for anomalous transport and Griffiths effects in one-dimensional (1D) quantum disordered isolated systems near the many-body localization (MBL) transition. The model is constituted by a collection of 1D tight-binding chains with on-site random energies, locally coupled to a weak perturbation, corresponding to a random matrix which belongs to the Gaussian orthogonal ensemble (GOE) and mimics the effect of thermal inclusions due to delocalizing interactions by providing a local broadening of the Poisson spectrum. While in the absence of such a coupling the model is localized as expected for the one-dimensional Anderson model, increasing the coupling with the GOE perturbation we find a transition to a conducting phase driven by the proliferation of quantum avalanches, which does not fit the standard paradigm of the Anderson transition. In particular an intermediate Griffiths region emerges, where exponentially distributed insulating segments coexist with a few rare resonances. Typical correlations decay exponentially fast, while average correlations decay as stretched exponentials and diverge with the length of the chain, indicating that the conducting inclusions have a fractal structure and that the localization length is broadly distributed at the critical point. This behavior is consistent with a Kosterlitz-Thouless-like criticality of the MBL transition. Transport and relaxation are dominated by rare resonances and rare strong insulating regions, and show anomalous behaviors strikingly similar to those observed in recent simulations and experiments in the bad metal delocalized phase preceding MBL. In particular, we find subdiffusive transport and power-law decay of the return probability at large times, with exponents that gradually change as one moves across the intermediate region. Concomitantly, the ac conductivity vanishes near zero frequency with an anomalous power law.
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 LLL subjected to a disordered magnetic field. By embedding a tunable thermal region of size PPP, 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/LP/LP/L exceeds a threshold value. Secondly, at strong disorder, we identify an intermediate phase between many-body localization and ergodic behavior as PPP 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.
The low-frequency response of systems near a many-body localization transition can be dominated by rare regions that are locally critical or ``in the other phase.'' It is known that in one dimension, these rare regions can cause the dc conductivity and diffusion constant to vanish even inside the delocalized thermal phase. Here, we present a general analysis of such Griffiths effects in the thermal phase near the many-body localization transition: we consider both one-dimensional and higher-dimensional systems, subject to quenched randomness, and discuss both linear response (including the frequency- and wave-vector-dependent conductivity) and more general dynamics. In all the regimes we consider, we identify observables that are dominated by rare-region effects. In some cases (one-dimensional systems and Floquet systems with no extensive conserved quantities), essentially all long-time local observables are dominated by rare-region effects; in others, generic observables are instead dominated by hydrodynamic long-time tails throughout the thermal phase, and one must look at specific probes, such as spin echo, to see Griffiths behavior.
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.
Over a decade of work has culminated in the consensus that one-dimensional systems, subject to sufficiently large disorder, fail to thermalize and possess an extensive set of local integrals of motion. In this Letter, I will provide numerical evidence for the contrary. In particular, this work studies the dynamics of disordered spin chains which are weakly coupled to a Markovian bath. Within this approach, the critical disorder for stability to quantum avalanches exceeds ${W}^{*}\ensuremath{\gtrsim}20$ in the random field Heisenberg chain. In stark contrast to the Anderson insulator, the avalanche threshold drifts considerably with system size, with no evidence of saturation in the studied regime.