We investigate the stability of the many-body localized phase for a system in contact with a single ergodic grain modeling a Griffiths region with low disorder. Our numerical analysis provides evidence that even a small ergodic grain consisting of only three qubits can delocalize a localized chain as soon as the localization length exceeds a critical value separating localized and extended regimes of the whole system. We present a simple theory, consistent with De Roeck and Huveneers's arguments in [Phys. Rev. B 95, 155129 (2017)PRBMDO2469-995010.1103/PhysRevB.95.155129] that assumes a system to be locally ergodic unless the local relaxation time determined by Fermi's golden rule is larger than the inverse level spacing. This theory predicts a critical value for the localization length that is perfectly consistent with our numerical calculations. We analyze in detail the behavior of local operators inside and outside the ergodic grain and find excellent agreement of numerics and theory.
We investigate the stability of the many-body localized phase for a system in contact with a single ergodic grain modeling a Griffiths region with low disorder. Our numerical analysis provides evidence that even a small ergodic grain consisting of only three qubits can delocalize a localized chain as soon as the localization length exceeds a critical value separating localized and extended regimes of the whole system. We present a simple theory, consistent with De Roeck and Huveneers's arguments in [Phys. Rev. B 95, 155129 (2017)PRBMDO2469-995010.1103/PhysRevB.95.155129] that assumes a system to be locally ergodic unless the local relaxation time determined by Fermi's golden rule is larger than the inverse level spacing. This theory predicts a critical value for the localization length that is perfectly consistent with our numerical calculations. We analyze in detail the behavior of local operators inside and outside the ergodic grain and find excellent agreement of numerics and theory.
Many-body localization plays an increasing role in condensed matter theory, both because it challenges the fundaments of statistical physics, and because it allows us to engineer several new, exotic, stable phases of matter. In this paper, the authors address the issue of the stability of a many-body localized material in contact with an ergodic grain, i.e., an imperfect bath made of a few interacting degrees of freedom. Thanks to detailed microscopic analysis and numerics, they conclude that such an ergodic grain eventually destabilizes the localized phase in the following cases: if the spatial dimension is higher than one, or if the spatial dimension is one but the localization length of the localized material is larger than a fixed threshold value. In realistic materials, these ergodic grains are always present as Griffiths regions where the disorder is anomalously small, and hence, the authors conclude that the localized phase in such materials is unstable, strictly speaking. Transport and thermalization are however exponentially suppressed in the distance between ergodic grains.
Recent work by De Roeck et al. [Phys. Rev. B 95, 155129 (2017)] has argued that many-body localization is unstable in two and higher dimensions due to a thermalization avalanche triggered by rare regions of weak disorder. To examine these arguments, we construct several models of a finite ergodic bubble coupled to an Anderson insulator of noninteracting fermions. We first describe the ergodic region using a Gaussian orthogonal ensemble random matrix and perform an exact diagonalization study of small systems. The results are in excellent agreement with a refined theory of the thermalization avalanche that includes transient finite-size effects, lending strong support to the avalanche scenario. We then explore the limit of large system sizes by modeling the ergodic region via a Hubbard model with all-to-all random hopping: the combined system, consisting of the bubble and the insulator, can be reduced to an effective Anderson impurity problem. We find that the spectral function of a local operator in the ergodic region changes dramatically when coupling to a large number of localized fermionic states; this occurs even when the localized sites are weakly coupled to the bubble. In principle, for a given size of the ergodic region, this may arrest the avalanche. However, this back-action effect is suppressed and the avalanche can be recovered if the ergodic bubble is large enough. Thus, the main effect of the back-action is to renormalize the critical bubble size.
When pushed out of equilibrium, generic interacting quantum systems equilibrate locally and are expected to evolve towards a locally thermal description despite their unitary time evolution. Systems in which disorder competes with interactions and transport can violate this expectation by exhibiting many-body localization. The strength of the disorder with respect to the other two parameters drives a transition from a thermalizing system towards a nonthermalizing one. The existence of this transition is well established both in experimental and numerical studies for finite systems. However, the stability of many-body localization in the thermodynamic limit is largely unclear. With increasing system size, a generic disordered system will contain with high probability areas of low disorder variation. If large and frequent enough, those areas constitute ergodic grains which can hybridize and thus compete with localization. While the details of this process are not yet settled, it is conceivable that if such regions appear sufficiently often, they might be powerful enough to restore thermalization. We set out to shed light on this problem by constructing potential landscapes with low disorder regions and numerically investigating their localization behavior in the Heisenberg model. Our findings suggest that many-body localization may be more stable than anticipated in other recent theoretical works.
Many-body localization provides a mechanism to avoid thermalization in isolated interacting quantum systems. The breakdown of thermalization may be complete, when all eigenstates in the many-body spectrum become localized, or partial, when the so-called many-body mobility edge separates localized and delocalized parts of the spectrum. Previously, De Roeck et al. [Phys. Rev. B 93, 014203 (2016)] suggested a possible instability of the many-body mobility edge in energy density. The local ergodic regions---so-called ``bubbles''---resonantly spread throughout the system, leading to delocalization. In order to study such instability mechanism, in this work we design a model featuring many-body mobility edge in particle density: the states at small particle density are localized, while increasing the density of particles leads to delocalization. Using numerical simulations with matrix product states, we demonstrate the stability of many-body localization with respect to small bubbles in large dilute systems for experimentally relevant timescales. In addition, we demonstrate that processes where the bubble spreads are favored over processes that lead to resonant tunneling, suggesting a possible mechanism behind the observed stability of many-body mobility edge. We conclude by proposing experiments to probe particle density mobility edge in the Bose-Hubbard model.
Closed generic quantum many-body systems may fail to thermalize under certain conditions even after long times, a phenomenon called many-body localization (MBL). Numerous studies support the stability of the MBL phase in strongly disordered one-dimensional systems. However, the situation is much less clear when a small part of the system is ergodic, a scenario which also has important implications for the existence of many-body localization in higher dimensions. Here we address this question experimentally using a large-scale quantum simulator of ultracold bosons in a two-dimensional optical lattice. We prepare two-component mixtures of varying relative population and implement a disorder potential which is only experienced by one of the components. The second non-disordered ''clean'' component plays the role of a bath of adjustable size that is collisionally coupled to the ''dirty'' component. Our experiments show how the dynamics of the dirty component, which, when on its own, show strong evidence of localization, become affected by the coupling to the clean component. For a high clean population, the clean component appears to behave as an effective bath for the system which leads to its delocalization, while for a smaller clean population, the ability of the bath to destabilize the system becomes strongly reduced. Our results reveal how a finite-sized quantum system can bring another one towards thermalization, in a regime of complex interplay between disorder, tunneling and intercomponent interactions. They provide a new benchmark for effective theories aiming to capture the complex physics of MBL in the weakly localized regime.
Closed generic quantum many-body systems may fail to thermalize under certain conditions even after long times, a phenomenon called many-body localization (MBL).Numerous studies support the stability of the MBL phase in strongly disordered one-dimensional systems.However, the situation is much less clear when a small part of the system is ergodic, a scenario which also has important implications for the existence of many-body localization in higher dimensions.Here we address this question experimentally using a largescale quantum simulator of ultracold bosons in a two-dimensional optical lattice.We prepare two-component mixtures of varying relative population and implement a disorder potential which is experienced only by one of the components.The second nondisordered "clean" component plays the role of a bath of adjustable size that is collisionally coupled to the "dirty" component.Our experiments show how the dynamics of the dirty component, which, when on its own, show strong evidence of localization, become affected by the coupling to the clean component.For a high clean population, the clean component appears to behave as an effective bath for the system which leads to its delocalization, while for a smaller clean population, the ability of the bath to destabilize the system becomes strongly reduced.Our results reveal how a finite-sized quantum system can bring another one towards thermalization, in a regime of complex interplay between disorder, tunneling, and intercomponent interactions.They provide a new benchmark for effective theories aiming to capture the complex physics of MBL in the weakly localized regime.