The many-body localization (MBL) problem continues to reveal its rich physical content. Herey, the authors investigate long-range correlations in the random-field Heisenberg chain and find rare, system-spanning ``cat states'' hidden in the high-energy many-body spectrum. These states are identified as nearly degenerate resonant pairs, well captured by an analytical ansatz. Remarkably, such atypical eigenstates persist even at very strong disorder, a regime in which conventional diagnostics suggest robust localization. These results offer fresh insight into the MBL problem and reopen fundamental questions about the mechanisms that might destabilize the MBL phase.
We investigate long-range resonances in quasiperiodic many-body localized (MBL) systems. Focusing on the Heisenberg chain in a deterministic Aubry-André potential, we complement standard diagnostics by analyzing the structure of long-distance pairwise correlations at high energy. Contrary to the expectation that the ergodic-MBL transition in quasiperiodic systems should be sharper due to the absence of Griffiths regions, we uncover a broad unconventional regime at strong quasiperiodic potential, characterized by fat-tailed distributions of longitudinal correlations at long distance. This reveals the presence of atypical eigenstates with strong long-range correlations in a regime where standard diagnostics indicate stable MBL. We further identify these anomalous eigenstates as quasidegenerate pairs of resonant cat states, which exhibit entanglement at long distance. These findings advance the understanding of quasiperiodic MBL and identify density-correlation measurements in ultracold atomic systems as a probe of long-range resonances.
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.
Motivated by recent debates around the many-body localization (MBL) problem, and in particular its stability against systemwide resonances, we investigate long-distance spin-spin correlations across the phase diagram of the random-field XXZ model, with a particular focus on the strong disorder regime. Building on state-of-the-art shift-invert diagonalization techniques, we study the high-energy behavior of transverse and longitudinal correlation functions, computed at the largest possible distance, for a broad range of disorder and interaction strengths. Our results show that while transverse correlations display a fairly stable exponential decay over the entire XXZ phase diagram, longitudinal correlations exhibit markedly different behavior, revealing distinct physical regimes. More precisely, we identify an intermediate disorder region where standard observables show well-converged MBL behavior [Colbois et al., Phys. Rev. Lett. 133, 116502 (2024)] while the distributions of longitudinal correlations reveal unexpected fat tails toward large values. These rare events strongly influence the average decay of longitudinal correlations, which we find to be algebraic in a broad region inside the supposed MBL phase, whereas the typical decay remains mostly exponential. At stronger disorder and weaker interactions, this intermediate regime is replaced by a more conventional exponential decay with short correlation lengths for both typical and average correlators, as expected for standard localization. Our findings shed light on the systemwide instabilities and raise important questions about the impact of such rare but large long-range correlations on the stability of the MBL phase. Finally, we discuss the possible fate of the intermediate region in the context of recent perspectives in the field.
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.
Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) – at least in one dimensional systems – leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder – in the thermal phase – we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in S_i^z S_{i+r}^z S i z S i + r z correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the S_i^zS_{i+r}^z S i z S i + r z correlator exhibits striking differences between the MBL and Anderson insulator cases.