No Arabic abstract
The many-body localised (MBL) to thermal crossover observed in exact diagonalisation studies remains poorly understood as the accessible system sizes are too small to be in an asymptotic scaling regime. We develop a model of the crossover in short 1D chains in which the MBL phase is destabilised by the formation of many-body resonances. The model reproduces several properties of the numerically observed crossover, including an apparent correlation length exponent $ u=1$, exponential growth of the Thouless time with disorder strength, linear drift of the critical disorder strength with system size, scale-free resonances, apparent $1/omega$ dependence of disorder-averaged spectral functions, and sub-thermal entanglement entropy of small subsystems. In the crossover, resonances induced by a local perturbation are rare at numerically accessible system sizes $L$ which are smaller than a emph{resonance length} $lambda$. For $L gg sqrt{lambda}$, resonances typically overlap, and this model does not describe the asymptotic transition. The model further reproduces controversial numerical observations which Refs. [v{S}untajs et al, 2019] and [Sels & Polkovnikov, 2020] claimed to be inconsistent with MBL. We thus argue that the numerics to date is consistent with a MBL phase in the thermodynamic limit.
Many-body localized (MBL) systems lie outside the framework of statistical mechanics, as they fail to equilibrate under their own quantum dynamics. Even basic features of MBL systems such as their stability to thermal inclusions and the nature of the dynamical transition to thermalizing behavior remain poorly understood. We study a simple model to address these questions: a two level system interacting with strength $J$ with $Ngg 1$ localized bits subject to random fields. On increasing $J$, the system transitions from a MBL to a delocalized phase on the emph{vanishing} scale $J_c(N) sim 1/N$, up to logarithmic corrections. In the transition region, the single-site eigenstate entanglement entropies exhibit bi-modal distributions, so that localized bits are either on (strongly entangled) or off (weakly entangled) in eigenstates. The clusters of on bits vary significantly between eigenstates of the emph{same} sample, which provides evidence for a heterogenous discontinuous transition out of the localized phase in single-site observables. We obtain these results by perturbative mapping to bond percolation on the hypercube at small $J$ and by numerical exact diagonalization of the full many-body system. Our results imply the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions $d$ that are high but finite.
We numerically study both the avalanche instability and many-body resonances in strongly-disordered spin chains exhibiting many-body localization (MBL). We distinguish between a finite-size/time MBL regime, and the asymptotic MBL phase, and identify some landmarks within the MBL regime. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition, in infinite-length systems, occurs much deeper in the MBL regime than has been suggested by most previous studies. Our other landmarks involve system-wide resonances. We find that the effective matrix elements producing eigenstates with system-wide resonances are enormously broadly distributed. This means that the onset of such resonances in typical samples occurs quite deep in the MBL regime, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regime in to three subregimes: (i) at strongest disorder, typical samples do not have any eigenstates that are involved in system-wide many-body resonances; (ii) there is a substantial intermediate regime where typical samples do have such resonances, but the pair of eigenstates with the minimum spectral gap does not; and (iii) in the weaker randomness regime, the minimum gap is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain non-resonant and the system thus still appears MBL in many respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i).
The transition between ergodic and many-body localized phases is expected to occur via an avalanche mechanism, in which emph{ergodic bubbles} that arise due to local fluctuations in system properties thermalize their surroundings leading to delocalization of the system, unless the disorder is sufficiently strong to stop this process. We propose an algorithm based on neural networks that allows to detect the ergodic bubbles using experimentally measurable two-site correlation functions. Investigating time evolution of the system, we observe a logarithmic in time growth of the ergodic bubbles in the MBL regime. The distribution of the size of ergodic bubbles converges during time evolution to an exponentially decaying distribution in the MBL regime, and a power-law distribution with a thermal peak in the critical regime, supporting thus the scenario of delocalization through the avalanche mechanism. Our algorithm permits to pin-point quantitative differences in time evolution of systems with random and quasiperiodic potentials, as well as to identify rare (Griffiths) events. Our results open new pathways in studies of the mechanisms of thermalization of disordered many-body systems and beyond.
Some interacting disordered many-body systems are unable to thermalize when the quenched disorder becomes larger than a threshold value. Although several properties of nonzero energy density eigenstates (in the middle of the many-body spectrum) exhibit a qualitative change across this many-body localization (MBL) transition, many of the commonly-used diagnostics only do so over a broad transition regime. Here, we provide evidence that the transition can be located precisely even at modest system sizes by sharply-defined changes in the distribution of extremal eigenvalues of the reduced density matrix of subsystems. In particular, our results suggest that $p* = lim_{lambda_2 rightarrow ln(2)^{+}}P_2(lambda_2)$, where $P_2(lambda_2)$ is the probability distribution of the second lowest entanglement eigenvalue $lambda_2$, behaves as an order-parameter for the MBL phase: $p*> 0$ in the MBL phase, while $p* = 0$ in the ergodic phase with thermalization. Thus, in the MBL phase, there is a nonzero probability that a subsystem is entangled with the rest of the system only via the entanglement of one subsystem qubit with degrees of freedom outside the region. In contrast, this probability vanishes in the thermal phase.
We theoretically study the response of a many-body localized system to a local quench from a quantum information perspective. We find that the local quench triggers entanglement growth throughout the whole system, giving rise to a logarithmic lightcone. This saturates the modified Lieb-Robinson bound for quantum information propagation in many-body localized systems previously conjectured based on the existence of local integrals of motion. In addition, near the localization-delocalization transition, we find that the final states after the local quench exhibit volume-law entanglement. We also show that the local quench induces a deterministic orthogonality catastrophe for highly excited eigenstates, where the typical wave-function overlap between the pre- and post-quench eigenstates decays {it exponentially} with the system size.