We propose a new approach to probing ergodicity and its breakdown in quantum many-body systems based on their response to a local perturbation. We study the distribution of matrix elements of a local operator between the systems eigenstates, finding
a qualitatively different behaviour in the many-body localized (MBL) and ergodic phases. To characterize how strongly a local perturbation modifies the eigenstates, we introduce the parameter ${cal G}(L)=langle ln (V_{nm}/delta) rangle$, which represents a disorder-averaged ratio of a typical matrix element of a local operator $V$ to the energy level spacing, $delta$; this parameter is reminiscent of the Thouless conductance in the single-particle localization. We show that the parameter ${cal G}(L)$ decreases with system size $L$ in the MBL phase, and grows in the ergodic phase. We surmise that the delocalization transition occurs when ${cal G}(L)$ is independent of system size, ${cal G}(L)={cal G}_csim 1$. We illustrate our approach by studying the many-body localization transition and resolving the many-body mobility edge in a disordered 1D XXZ spin-1/2 chain using exact diagonalization and time-evolving block decimation methods. Our criterion for the MBL transition gives insights into microscopic details of transition. Its direct physical consequences, in particular logarithmically slow transport at the transition, and extensive entanglement entropy of the eigenstates, are consistent with recent renormalization group predictions.
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.
