No Arabic abstract
We compute and compare the decay lengths of several correlation functions and effective coupling constants in the many-body localized (MBL) phase. To this end, we consider the distribution of the logarithms of these couplings and correlators: in each case the log-coupling follows a normal distribution with mean and variance that grow linearly with separation. Thus, a localization length is asymptotically sharply defined for each of these quantities. These localization lengths differ numerically from one another, but all of them remain short up to the numerically observed MBL transition, indicating stability of the MBL phase against isolated ergodic inclusions. We also show how these broad distributions may be extracted using interferometric probes such as double electron-electron resonance (DEER) and the statistics of local spin precession frequencies.
Many-body localization is a striking mechanism that prevents interacting quantum systems from thermalizing. The absence of thermalization behaviour manifests itself, for example, in a remanence of local particle number configurations, a quantity that is robust over a parameter range. Local particle numbers are directly accessible in programmable quantum simulators, in systems of cold atoms even in two spatial dimensions. Yet, the classical simulation aimed at building trust in quantum simulations is highly challenging. In this work, we present a comprehensive tensor network simulation of a many-body localized systems in two spatial dimensions using a variant of an iPEPS algorithm. The required translational invariance can be restored by implementing the disorder into an auxiliary spin system, providing an exact disorder average under dynamics. We can quantitatively assess signatures of many-body localization for the infinite system: Our methods are powerful enough to provide crude dynamical estimates for the transition between localized and ergodic phases. Interestingly, in this setting of finitely many disorder values, which we also compare with simulations involving non-interacting fermions and for which we discuss the emergent physics, localization emerges in the interacting regime, for which we provide an intuitive argument, while Anderson localization is absent.
We investigate dynamical quantum phase transitions in disordered quantum many-body models that can support many-body localized phases. Employing $l$-bits formalism, we lay out the conditions for which singularities indicative of the transitions appear in the context of many-body localization. Using the combination of the mapping onto $l$-bits and exact diagonalization results, we explicitly demonstrate the presence of these singularities for a candidate model that features many-body localization. Our work paves the way for understanding dynamical quantum phase transitions in the context of many-body localization, and elucidating whether different phases of the latter can be detected from analyzing the former. The results presented are experimentally accessible with state-of-the-art ultracold-atom and ion-trap setups.
Strongly correlated systems can exhibit surprising phenomena when brought in a state far from equilibrium. A spectacular example are quantum avalanches, that have been predicted to run through a many-body--localized system and delocalize it. Quantum avalanches occur when the system is locally coupled to a small thermal inclusion that acts as a bath. Here we realize an interface between a many-body--localized system and a thermal inclusion of variable size, and study its dynamics. We find evidence for accelerated transport into the localized region, signature of a quantum avalanche. By measuring the site-resolved entropy we monitor how the avalanche travels through the localized system and thermalizes it site by site. Furthermore, we isolate the bath-induced dynamics by evaluating multipoint correlations between the bath and the system. Our results have fundamental implications on the robustness of many-body--localized systems and their critical behavior.
Quantum coherence quantifies the amount of superposition a quantum state can have in a given basis. Since there is a difference in the structure of eigenstates of the ergodic and many-body localized systems, we expect them also to differ in terms of their coherences in a given basis. Here, we numerically calculate different measures of quantum coherence in the excited eigenstates of an interacting disordered Hamiltonian as a function of the disorder. We show that quantum coherence can be used as an order parameter to detect the well-studied ergodic to many-body-localized phase transition. We also perform quantum quench studies to distinguish the behavior of coherence in thermalized and localized phases. We then present a protocol to calculate measurement-based localizable coherence to investigate the thermal and many-body localized phases. The protocol allows one to investigate quantum correlations experimentally in a non-destructive way, in contrast to measures that require tracing out a subsystem, which always destroys coherence and correlation.
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).