Topology and many-body localization (MBL) have opened new avenues for preserving quantum information at finite energy density. Resonant delocalization plays a crucial role in destabilizing these phenomena. In this paper, we study the statistical properties of many-body resonances in a disordered interacting Ising model---which can host symmetry-protected topological order---using a Clifford circuit encoding of the real-space renormalization group which allows the resonant properties of the wave functions to be efficiently characterized. Our findings show that both the trivial and topologically ordered MBL phases remain stable in the presence of the resonances but, in the vicinity of the transition between them, localization is destabilized by resonance proliferation. Diverging susceptibility towards the development of an avalanche instability suggests an intervening ergodic phase. We are also able to access the local integrals of motion in the MBL phases and identify the topological edge-mode operators in the ordered phase. Our results have important implications for the stability of MBL and phase transitions between distinct MBL phases with and without symmetries.
Topology and many-body localization (MBL) have opened new avenues for preserving quantum information at finite energy density. Resonant delocalization plays a crucial role in destabilizing these phenomena. In this paper, we study the statistical properties of many-body resonances in a disordered interacting Ising model---which can host symmetry-protected topological order---using a Clifford circuit encoding of the real-space renormalization group which allows the resonant properties of the wave functions to be efficiently characterized. Our findings show that both the trivial and topologically ordered MBL phases remain stable in the presence of the resonances but, in the vicinity of the transition between them, localization is destabilized by resonance proliferation. Diverging susceptibility towards the development of an avalanche instability suggests an intervening ergodic phase. We are also able to access the local integrals of motion in the MBL phases and identify the topological edge-mode operators in the ordered phase. Our results have important implications for the stability of MBL and phase transitions between distinct MBL phases with and without symmetries.
Models of many-body localization (MBL) exhibit slow numerical drifts towards delocalization with increasing system size, for which no satisfactory theory exists. Numerical results indicate that these drifts are driven by the proliferation of many-body resonances at intermediate disorder strengths. We develop a statistical method to predict the distribution of resonance oscillation frequencies, accounting for how resonances formed at larger frequency scales affect the subsequent formation of resonances at lower frequencies. Working within the statistical Jacobi approximation (SJA), we derive a flow equation for a power-law exponent θ(w) characterizing the density of resonances at frequency scale w. A localized phase is described by a line of fixed points with θ(w)>0, while an instability of the flow signals resonance proliferation and the onset of thermalization. The predicted θ(w) agrees with numerical results for the Anderson model on random regular graphs and the Lévy--Rosenzweig--Porter random matrix ensemble, both of which host resonance-driven delocalization transitions. We further connect the flow to eigenstate properties such as the participation ratio and to dynamical observables such as the return probability. The predicted θ(w) also agrees with numerical measurements in real-space models of MBL at intermediate disorder strengths, providing a step towards explaining the finite-size drifts observed in MBL.
The effect of rare system-wide resonances in the many-body localization (MBL) transition has recently attracted significant attention. They are expected to play a prominent role in the stability of the MBL phase, prompting the development of new theoretical frameworks to properly account for their statistical weight. We employ a method based on an analogy with mean-field disordered glassy systems to characterize the statistics of transmission amplitudes between distant many-body configurations in Hilbert space, and apply it to the random-field XXZ spin chain. By introducing a Lagrange multiplier, which formally plays the role of an effective temperature controlling the influence of extreme outliers in the heavy-tailed distribution of propagators, we identify three distinct regimes: (i) an ergodic phase with uniform spreading in Hilbert space, (ii) an intermediate regime where delocalization is driven by rare, disorder-dependent long-range resonances, and (iii) a robust MBL phase where such resonances cannot destabilize localization. We derive a finite-size phase diagram in the disorder--interaction plane both in the spin and in the Anderson basis that quantitatively agrees with recent numerical results based on real-space spin-spin correlation functions. We further demonstrate that even infinitesimal interactions can destroy the Anderson insulator at finite disorder, with the critical disorder remaining finite down to small interaction strengths. By visualizing resonant transmission pathways on the Hilbert space graph, we provide a complementary perspective to real-space and spectral probes, revealing how the destabilization of the MBL phase at finite sizes stems from the emergence of resonant paths that become progressively rarer and shorter-ranged deep in the localized phase.
We examine the stability of marginally Anderson localized phase transitions between localized phases to the addition of many-body interactions, focusing in particular on the spin-glass to paramagnet transition in a disordered transverse field Ising model in one dimension. We find evidence for a perturbative instability of localization at finite energy densities once interactions are added, i.e. evidence for the relevance of interactions - in a renormalization group sense - to the non-interacting critical point governed by infinite randomness scaling. We introduce a novel diagnostic, the "susceptibility of entanglement", which allows us to perturbatively probe the effect of adding interactions on the entanglement properties of eigenstates, and helps us elucidate the resonant processes that can cause thermalization. The susceptibility serves as a much more sensitive probe, and its divergence can detect the perturbative beginnings of an incipient instability even in regimes and system sizes for which conventional diagnostics point towards localization. We expect this new measure to be of independent interest for analyzing the stability of localization in a variety of different settings.
We develop a real-space renormalization-group (RSRG) scheme by appropriately inserting the long-range hopping t∼r^{-α} with nearest-neighbor interaction to study the entanglement entropy and maximum block size for the many-body localization (MBL) transition. We show that for α<2 there exists a localization transition with renormalized disorder that depends logarithmically on the finite size of the system. The transition observed for α>2 does not need a rescaling in disorder strength. Most important, we find that even though the MBL transition for α>2 falls in the same universality class as that of the short-range models, the transition for α<2 belongs to a different universality class. Because of the intrinsic nature of the RSRG flow toward delocalization, MBL phase for α>2 might suffer an instability in the thermodynamic limit while the underlying systems support algebraic localization. Moreover, we verify these findings by inserting microscopic details to the RSRG scheme where we additionally find a more appropriate rescaling function for disorder strength; we indeed uncover a power-law scaling with a logarithmic correction and a distinctly different stretched exponential scaling for α<2 and α>2, respectively, by analyzing system with finite size. This finding further suggests that microscopic RSRG scheme is able to give a hint of instability of the MBL phase for α>2 even considering systems of finite size.
Strongly disordered systems in the many-body localized (MBL) phase can exhibit ground state order in highly excited eigenstates. The interplay between localization, symmetry, and topology has led to the characterization of a broad landscape of MBL phases ranging from spin glasses and time crystals to symmetry protected topological phases. Understanding the nature of phase transitions between these different forms of eigenstate order remains an essential open question. Here, we conjecture that no direct transition between distinct MBL orders can occur in one dimension; rather, an ergodic phase always intervenes. Motivated by recent advances in Rydberg-atom-based quantum simulation, we propose an experimental protocol where the intervening ergodic phase can be diagnosed via the dynamics of local observables.
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).