We analyze the unusual slow dynamics that emerges in the bad metal delocalized phase preceding the many-body localization transition by using single-particle Anderson localization on the Bethe lattice as a toy model of many-body dynamics in Fock space. We probe the dynamical evolution by measuring observables such as the imbalance and equilibrium correlation functions, which display slow dynamics and power laws strikingly similar to those observed in recent simulations and experiments. We relate this unusual behavior to the nonergodic spectral statistics found on Bethe lattices. We discuss different scenarios, such as a true intermediate phase which persists in the thermodynamic limit versus a glassy regime established on finite but very large time and length scales only, and their implications for real-space dynamical properties. In the latter, slow dynamics and power laws extend on a very large time window but are eventually cut off on a time scale that diverges at the many-body localization transition.
We analyze the unusual slow dynamics that emerges in the bad metal delocalized phase preceding the many-body localization transition by using single-particle Anderson localization on the Bethe lattice as a toy model of many-body dynamics in Fock space. We probe the dynamical evolution by measuring observables such as the imbalance and equilibrium correlation functions, which display slow dynamics and power laws strikingly similar to those observed in recent simulations and experiments. We relate this unusual behavior to the nonergodic spectral statistics found on Bethe lattices. We discuss different scenarios, such as a true intermediate phase which persists in the thermodynamic limit versus a glassy regime established on finite but very large time and length scales only, and their implications for real-space dynamical properties. In the latter, slow dynamics and power laws extend on a very large time window but are eventually cut off on a time scale that diverges at the many-body localization transition.
Using the noninteracting Anderson tight-binding model on the Bethe lattice as a toy model for the many-body quantum dynamics, we propose a novel and transparent theoretical explanation of the anomalously slow dynamics that emerges in the bad metal phase preceding the many-body localization (MBL) transition. By mapping the time decorrelation of many-body wave functions onto directed polymers in random media, we show the existence of a glass transition within the extended regime of loopless Cayley trees separating a metalliclike phase at small disorder, where delocalization occurs on an exponential number of paths, from a bad metal-like phase at intermediate disorder, where resonances are formed on rare, specific, disorder dependent site orbitals on very distant generations. The physical interpretation of subdiffusion and nonexponential relaxation emerging from this picture is complementary to the Griffiths one, although both scenarios rely on the presence of heavy-tailed distribution of the escape times. We relate the dynamical evolution in the glassy phase to the depinning transition of directed polymers, which results in macroscopic and abrupt jumps of the preferred delocalizing paths when a parameter like the energy is varied, and produce a singular behavior of the overlap correlation function between eigenstates at different energies. By comparing the quantum dynamics on loopless Cayley trees and random regular graphs, we discuss the effect of loops, showing that in the latter slow dynamics and apparent power laws extend on a very large time window but are eventually cutoff on a timescale that diverges at the MBL transition.
We investigate dynamics of a single mobile impurity immersed in a bath of Anderson localized particles and focus on the regime of relatively strong disorder and interactions. In that regime, the dynamics of the system is particularly slow, suggesting, at short times, an occurrence of many-body localization. Considering longer timescales, we show that the latter is a transient effect and that, eventually, the impurity spreads subdiffusively and induces a gradual delocalization of the Anderson insulator. The phenomenology of the system in the considered regime of slow dynamics includes a subdiffusive growth of mean square displacement of the impurity, power-law decay of density correlation functions of the Anderson insulator, and a power-law growth of entanglement entropy in the system. We observe a similar regime of slow dynamics also when the disorder in the system is replaced by a sufficiently strong quasiperiodic potential.
We study the delocalization dynamics of interacting disordered hard-core bosons for quasi-1D and 2D geometries, with system sizes and timescales comparable to state-of-the-art experiments. The results are strikingly similar to the 1D case, with slow, subdiffusive dynamics featuring power-law decay. From the freezing of this decay we infer the critical disorder W_{c}(L,d) as a function of length L and width d. In the quasi-1D case W_{c} has a finite large-L limit at fixed d, which increases strongly with d. In the 2D case W_{c}(L,L) grows with L. The results are consistent with the avalanche picture of the many-body localization transition.
Characterizing the delocalization transition in closed quantum systems with a many-body localized phase is a key open question in the field of nonequilibrium physics. We exploit the fact that localization of particles as realized in Anderson localization and standard many-body localization (MBL) implies Fock-space localization in single-particle basis sets characterized by a real-space index. Using a recently introduced quantitative measure for Fock-space localization computed from the density distributions, the occupation distance, we systematically study its scaling behavior across delocalization transitions and identify critical points from scaling collapses of numerical data. Excellent agreement with literature results is found for the critical disorder strengths of noninteracting fermions, such as the one-dimensional Aubry-Andr\'e model and the three-dimensional Anderson model. We observe a distinctively different scaling behavior in the case of interacting fermions with random disorder consistent with a Kosterlitz-Thouless transition. Finally, we use our measure to extract the transition point as a function of filling for interacting fermions.
Disordered quantum systems undergoing a many-body localization (MBL) transition fail to reach thermal equilibrium under their own dynamics. Distinguishing between asymptotically localized or delocalized dynamics based on numerical results is, however, nontrivial due to finite-size effects. Numerical linked cluster expansions (NLCE) provide a means to tackle quantum systems directly in the thermodynamic limit but are challenging for models without translational invariance. Here we demonstrate that NLCE provides a powerful tool to explore MBL by simulating quench dynamics in disordered spin-1/2 two-leg ladders and Fermi-Hubbard chains. Combining NLCE with an efficient real-time evolution of pure states, we obtain converged results for the decay of the imbalance on long time scales and show that, especially for intermediate disorder below the putative MBL transition, NLCE outperforms direct simulations of finite systems with open or periodic boundaries. Furthermore, while spin is delocalized even in strongly disordered Hubbard chains with frozen charge, we unveil that an additional tilted potential leads to a drastic slowdown of the spin imbalance and nonergodic behavior on accessible times. Our work sheds light on MBL in systems beyond the well-studied disordered Heisenberg chain and emphasizes the usefulness of NLCE for this purpose.
We introduce and study a toy model for anomalous transport and Griffiths effects in one-dimensional (1D) quantum disordered isolated systems near the many-body localization (MBL) transition. The model is constituted by a collection of 1D tight-binding chains with on-site random energies, locally coupled to a weak perturbation, corresponding to a random matrix which belongs to the Gaussian orthogonal ensemble (GOE) and mimics the effect of thermal inclusions due to delocalizing interactions by providing a local broadening of the Poisson spectrum. While in the absence of such a coupling the model is localized as expected for the one-dimensional Anderson model, increasing the coupling with the GOE perturbation we find a transition to a conducting phase driven by the proliferation of quantum avalanches, which does not fit the standard paradigm of the Anderson transition. In particular an intermediate Griffiths region emerges, where exponentially distributed insulating segments coexist with a few rare resonances. Typical correlations decay exponentially fast, while average correlations decay as stretched exponentials and diverge with the length of the chain, indicating that the conducting inclusions have a fractal structure and that the localization length is broadly distributed at the critical point. This behavior is consistent with a Kosterlitz-Thouless-like criticality of the MBL transition. Transport and relaxation are dominated by rare resonances and rare strong insulating regions, and show anomalous behaviors strikingly similar to those observed in recent simulations and experiments in the bad metal delocalized phase preceding MBL. In particular, we find subdiffusive transport and power-law decay of the return probability at large times, with exponents that gradually change as one moves across the intermediate region. Concomitantly, the ac conductivity vanishes near zero frequency with an anomalous power law.