We construct an analytic theory of many-body localization (MBL) in random spin chains. The approach is based on a first quantized perspective in which MBL is understood as a localization phenomenon on the high dimensional lattice defined by the discrete Hilbert space of the clean system. We construct a field theory on that lattice and apply it to discuss the stability of a weak disorder (`Wigner-Dyson) and a strong disorder (`Poisson) phase.
Polynomially filtered exact diagonalization method (POLFED) for large sparse matrices is introduced. The algorithm finds an optimal basis of a subspace spanned by eigenvectors with eigenvalues close to a specified energy target by a spectral transformation using a high order polynomial of the matrix. The memory requirements scale better with system size than in the state-of-the-art shift-invert approach. The potential of POLFED is demonstrated examining many-body localization transition in 1D interacting quantum spin-1/2 chains. We investigate the disorder strength and system size scaling of Thouless time. System size dependence of bipartite entanglement entropy and of the gap ratio highlights the importance of finite-size effects in the system. We discuss possible scenarios regarding the many-body localization transition obtaining estimates for the critical disorder strength.
We investigate a many-body localization transition based on a Boltzmann transport theory. Introducing weak localization corrections into a Boltzmann equation, Hershfield and Ambegaokar re-derived the Wolfle-Vollhardt self-consistent equation for the diffusion coefficient [Phys. Rev. B {bf 34}, 2147 (1986)]. We generalize this Boltzmann equation framework, introducing electron-electron interactions into the Hershfield-Ambegaokar Boltzmann transport theory based on the study of Zala-Narozhny-Aleiner [Phys. Rev. B {bf 64}, 214204 (2001)]. Here, not only Altshuler-Aronov corrections but also dephasing effects are taken into account. As a result, we obtain a self-consistent equation for the diffusion coefficient in terms of the disorder strength and temperature, which extends the Wolfle-Vollhardt self-consistent equation in the presence of electron correlations. Solving our self-consistent equation numerically, we find a many-body localization insulator-metal transition, where a metallic phase appears from dephasing effects dominantly instead of renormalization effects at high temperatures. Although this mechanism is consistent with that of recent seminal papers [Ann. Phys. (N. Y). {bf 321}, 1126 (2006); Phys. Rev. Lett. {bf 95}, 206603 (2005)], we find that our three-dimensional metal-insulator transition belongs to the first order transition, which differs from the Anderson metal-insulator transition described by the Wolfle-Vollhardt self-consistent theory. We speculate that a bimodal distribution function for the diffusion coefficient is responsible for this first order phase transition.
Recent developments in matrix-product-state (MPS) investigations of many-body localization (MBL) are reviewed, with a discussion of benefits and limitations of the method. This approach allows one to explore the physics around the MBL transition in systems much larger than those accessible to exact diagonalization. System sizes and length scales that can be controllably accessed by the MPS approach are comparable to those studied in state-of-the-art experiments. Results for 1D, quasi-1D, and 2D random systems, as well as 1D quasi-periodic systems are presented. On time scales explored (up to $t approx 300$ in units set by the hopping amplitude), a slow, subdiffusive transport in a rather broad disorder range on the ergodic side of the MBL transition is found. For 1D random spin chains, which serve as a standard model of the MBL transition, the MPS study demonstrates a substantial drift of the critical point $W_c(L)$ with the system size $L$: while for $L approx 20$ we find $W_c approx 4$, as also given by exact diagonalization, the MPS results for $L = 50$--100 provide evidence that the critical disorder saturates, in the large-$L$ limit, at $W_c approx 5.5$. For quasi-periodic systems, these finite-size effects are much weaker, which suggests that they can be largely attributed to rare events. For quasi-1D ($dtimes L$, with $d ll L$) and 2D ($Ltimes L$) random systems, the MPS data demonstrate an unbounded growth of $W_c$ in the limit of large $d$ and $L$, in agreement with analytical predictions based on the rare-event avalanche theory.
We discuss the problem of localization in two dimensional electron systems in the quantum Hall (single Landau level) regime. After briefly summarizing the well-studied problem of Anderson localization in the non-interacting case, we concentrate on the problem of disorder induced many-body localization (MBL) in the presence of electron-electron interactions using numerical exact diagonalization and eigenvalue spacing statistics as a function of system size. We provide evidence showing that MBL is not attainable in a single Landau level with short range (white noise) disorder in the thermodynamic limit. We then study the interplay of topology and localization, by contrasting the behavior of topological and nontopological subbands arising from a single Landau level in two models - (i) a pair of extremely flat Hofstadter bands with an optimally chosen periodic potential, and (ii) a Landau level with a split-off nontopological impurity band. Both models provide convincing evidence for the strong effect of topology on the feasibility of many-body localization as well as slow dynamics starting from a nonequilibrium state with charge imbalance.
We propose a method for detecting many-body localization (MBL) in disordered spin systems. The method involves pulsed, coherent spin manipulations that probe the dephasing of a given spin due to its entanglement with a set of distant spins. It allows one to distinguish the MBL phase from a non-interacting localized phase and a delocalized phase. In particular, we show that for a properly chosen pulse sequence the MBL phase exhibits a characteristic power-law decay reflecting its slow growth of entanglement. We find that this power-law decay is robust with respect to thermal and disorder averaging, provide numerical simulations supporting our results, and discuss possible experimental realizations in solid-state and cold atom systems.