Despite enormous efforts devoted to the study of the many-body localization (MBL) phenomenon, the nature of the high-energy behavior of the Heisenberg spin chain in a strong random magnetic field is lacking consensus. Here, we take a step back by exploring the weak interaction limit starting from the Anderson localized (AL) insulator. Through shift-invert diagonalization, we find that, below a certain disorder threshold h^{*}, weak interactions necessarily lead to an ergodic instability, whereas at strong disorder the AL insulator directly turns into MBL, in agreement with a simple interpretation of the avalanche theory for restoration of ergodicity. We further map the phase diagram for the generic XXZ model in the disorder h-interaction Δ plane. Taking advantage of the total magnetization conservation, our results unveil the remarkable behavior of the spin-spin correlation functions: in the regime indicated as MBL by standard observables, their exponential decay undergoes an inversion of orientation ξ_{z}>ξ_{x}. We find that the longitudinal length ξ_{z} is a key quantity for capturing ergodic instabilities, as it increases with system size near the thermal phase, in sharp contrast to its transverse counterpart.
Despite enormous efforts devoted to the study of the many-body localization (MBL) phenomenon, the nature of the high-energy behavior of the Heisenberg spin chain in a strong random magnetic field is lacking consensus. Here, we take a step back by exploring the weak interaction limit starting from the Anderson localized (AL) insulator. Through shift-invert diagonalization, we find that, below a certain disorder threshold h^{*}, weak interactions necessarily lead to an ergodic instability, whereas at strong disorder the AL insulator directly turns into MBL, in agreement with a simple interpretation of the avalanche theory for restoration of ergodicity. We further map the phase diagram for the generic XXZ model in the disorder h-interaction Δ plane. Taking advantage of the total magnetization conservation, our results unveil the remarkable behavior of the spin-spin correlation functions: in the regime indicated as MBL by standard observables, their exponential decay undergoes an inversion of orientation ξ_{z}>ξ_{x}. We find that the longitudinal length ξ_{z} is a key quantity for capturing ergodic instabilities, as it increases with system size near the thermal phase, in sharp contrast to its transverse counterpart.
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.
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.
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 (SiSi+rz)-S-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 (SiSi+rz)-S-z correlator exhibits striking differences between the MBL and Anderson insulator cases.
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.
Despite a very good understanding of single-particle Anderson localization in one-dimensional (1D) disordered systems, many-body effects are still full of surprises; a famous example is the interaction-driven many-body localization (MBL) problem, about which much has been written, and perhaps the best is yet to come. Interestingly enough the noninteracting limit provides a natural playground to study nontrivial multiparticle physics, offering the possibility to test some general mechanisms with very large-scale exact diagonalization simulations. In this paper, we first revisit the 1D many-body Anderson insulator through the lens of extreme value theory, focusing on the extreme polarizations of the equivalent spin chain model in a random magnetic field. A many-body-induced chain breaking mechanism is explored numerically, and compared to an analytically solvable toy model. A unified description, from weak to large disorder strengths W emerges, where the disorder-dependent average localization length $\ensuremath{\xi}(W)$ governs the extreme events leading to chain breaks. In particular, tails of the local magnetization distributions are controlled by $\ensuremath{\xi}(W)$. Remarkably, we also obtain a quantitative understanding of the full distribution of the extreme polarizations, which is given by a Fr\'echet-type law. In a second part, we explore finite interaction physics and the MBL question. For the available system sizes, we numerically quantify the difference in the extreme value distributions between the interacting problem and the noninteracting Anderson case. Strikingly, we observe a sharp ``extreme-statistics transition'' as W changes, which may coincide with the MBL transition.
In this paper, we report our numerical analysis of energy level spacing statistics for one-dimensional spin-1/2 XXZ model with random on-site longitudinal magnetic fields Bi. We concentrate on the strong disorder limit where the system is expected to be in either a paramagnetic many-body localized (MBL) state or a spin-glass MBL state. By analyzing the energy-level spacing statistics as a function of strength of random magnetic field h/J⊥ (−h≤Bi≤h), interaction strength Jz/J⊥, energy of the many-body state E and the number of spin-↑ particles in the system M=∑i(siz+21), we show that there exists a small region between the paramagnetic and spin-glass MBL phases where an ergodic phase emerges. The ergodic phase emerges at the middle of the many-body energy spectrum when M∼2N and is not adiabatically connected to the ergodic phase that occurs in the weak-disorder, weak-interaction limit.