Within the Thouless energy scale, many-body Hamiltonians behave like random matrices. We study the Thouless energy of the disordered Heisenberg quantum spin chain. This model has become the standard model for the study of many-body localization and the transition to an ergodic phase for weak disorder. We focus on the ergodic side of the Many-Body Localization transition. The Thouless energy extracted from long range spectral statistics and from the power spectrum of the fluctuations of the full momentum distribution coincide. This scale behaves non-monotonically and is peaked at an intermediate value of the disorder. An analysis of the scaling behavior of the Thouless energy and its comparison with the number of populated eigenstates after a quench suggests that this scale is not large enough to guarantee thermalization in the thermodynamic limit.
Within the Thouless energy scale, many-body Hamiltonians behave like random matrices. We study the Thouless energy of the disordered Heisenberg quantum spin chain. This model has become the standard model for the study of many-body localization and the transition to an ergodic phase for weak disorder. We focus on the ergodic side of the Many-Body Localization transition. The Thouless energy extracted from long range spectral statistics and from the power spectrum of the fluctuations of the full momentum distribution coincide. This scale behaves non-monotonically and is peaked at an intermediate value of the disorder. An analysis of the scaling behavior of the Thouless energy and its comparison with the number of populated eigenstates after a quench suggests that this scale is not large enough to guarantee thermalization in the thermodynamic limit.
We study the ergodic side of the many-body localization transition in its standard model, the disordered Heisenberg quantum spin chain. We show that the Thouless energy, extracted from long-range spectral statistics and the power-spectrum of the full momentum distribution fluctuations, is not large enough to guarantee thermalization. We find that both estimates coincide and behave nonmonotonically, exhibiting a strong peak at an intermediate value of the disorder. Furthermore, we show that nonthermalizing initial conditions occur well within the ergodic phase with larger probability than expected. Finally, we propose a mechanism, driven by the Thouless energy and the presence of anomalous events, for the transition to the localized phase.
We study a one-dimensional XXZ spin chain in a random field on the metallic side of the many-body localization transition by level statistics. For a fixed interaction, and intermediate disorder below the many-body localization transition, we find that, asymptotically, the number variance grows faster than linear with a disorder-dependent exponent. This is consistent with the existence of an anomalous Thouless energy in the spectrum. In noninteracting disordered metals, this is an energy scale related to the typical time for a particle to diffuse across the sample. In the interacting case, it seems related to a more intricate anomalous diffusion process. This interpretation is not fully consistent with recent claims that for intermediate disorder, level statistics are described by a plasma model with power-law decaying interactions whose number variance grows slower than linear. As disorder is further increased, still on the metallic side, the Thouless energy is gradually washed out. In the range of sizes we can explore, level statistics are scale invariant and approach Poisson statistics at the many-body localization transition. Slightly below the many-body localization transition, spectral correlations, well described by critical statistics, are quantitatively similar to those of a high-dimensional, noninteracting, disordered conductor at the Anderson transition.
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.
Disordered interacting spin chains that undergo a many-body localization transition are characterized by two limiting behaviors where the dynamics are chaotic and integrable. However, the transition region between them is not fully understood yet. We propose here a possible finite-size precursor of a critical point that shows a typical finite-size scaling and distinguishes between two different dynamical phases. The kurtosis excess of the diagonal fluctuations of the full one-dimensional momentum distribution from its microcanonical average is maximum at this singular point in the paradigmatic disordered J_1 J 1 - J_2 J 2 model. For system sizes accessible to exact diagonalization, both the position and the size of this maximum scale linearly with the system size. Furthermore, we show that this singular point is found at the same disorder strength at which the Thouless and the Heisenberg energies coincide. Below this point, the spectral statistics follow the universal random matrix behavior up to the Thouless energy. Above it, no traces of chaotic behavior remain, and the spectral statistics are well described by a generalized semi-Poissonian model, eventually leading to the integrable Poissonian behavior. We provide, thus, an integrated scenario for the many-body localization transition, conjecturing that the critical point in the thermodynamic limit, if it exists, should be given by this value of disorder strength.
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 random field S=21 Heisenberg chain exhibits a dynamical many body localization transition at a critical disorder strength, which depends on the energy density. At weak disorder, the eigenstate thermalization hypothesis (ETH) is fulfilled on average, making local observables smooth functions of energy, whose eigenstate-to-eigenstate fluctuations decrease exponentially with system size. We demonstrate the validity of ETH in the thermal phase as well as its breakdown in the localized phase and show that rare states exist which do not strictly follow ETH, becoming more frequent closer to the transition. Similarly, the probability distribution of the entanglement entropy at intermediate disorder develops long tails all the way down to zero entanglement. We propose that these low entanglement tails stem from localized regions at the subsystem boundaries which were recently discussed as a possible mechanism for subdiffusive transport in the ergodic phase.