We construct and solve a classical percolation model with a phase transition that we argue acts as a proxy for the quantum many-body localization transition. The classical model is defined on a graph in the Fock space of a disordered, interacting quantum spin chain, using a convenient choice of basis. Edges of the graph represent matrix elements of the spin Hamiltonian between pairs of basis states that are expected to hybridize strongly. At weak disorder, all nodes are connected, forming a single cluster. Many separate clusters appear above a critical disorder strength, each typically having a size that is exponentially large in the number of spins but a vanishing fraction of the Fock-space dimension. We formulate a transfer matrix approach that yields an exact value $\ensuremath{\nu}=2$ for the localization length exponent, and also use complete enumeration of clusters to study the transition numerically in finite-sized systems.
We use exact diagonalization to explore the many-body localization transition in a random-field spin-1/2 chain. We examine the correlations within each many-body eigenstate, looking at all high-energy states and thus effectively working at infinite temperature. For weak random field the eigenstates are thermal, as expected in this nonlocalized, ``ergodic'' phase. For strong random field the eigenstates are localized with only short-range entanglement. We roughly locate the localization transition and examine some of its finite-size scaling, finding that this quantum phase transition at nonzero temperature might be showing infinite-randomness scaling with a dynamic critical exponent $z\ensuremath{\rightarrow}\ensuremath{\infty}$.
We study spectral and wave-function statistics for many-body localization transition in systems with long-range interactions decaying as $1/{r}^{\ensuremath{\alpha}}$ with an exponent $\ensuremath{\alpha}$ satisfying $d\ensuremath{\le}\ensuremath{\alpha}\ensuremath{\le}2d$, where
d is the spatial dimensionality. We refine earlier arguments and show that the system undergoes a localization transition as a function of the rescaled disorder ${W}^{*}=W/{L}^{2d\ensuremath{-}\ensuremath{\alpha}}lnL$, where
W is the disorder strength and
L the system size. This transition has much in common with that on random regular graphs. We further perform a detailed analysis of the inverse participation ratio (IPR) of many-body wave functions, exploring how ergodic behavior in the delocalized phase switches to fractal one at the critical point and on the localized side of the transition. Our analytical results for the scaling of the critical disorder
W with the system size
L and for the scaling of IPR in the delocalized and localized phases are supported and corroborated by exact diagonalization of spin chains.
The space of one-dimensional disordered interacting quantum models displaying a many-body localization (MBL) transition seems sufficiently rich to produce critical points with level statistics interpolating continuously between the Poisson statistics of the localized phase and the Wigner–Dyson statistics of the delocalized phase. In this paper, we consider the strong disorder limit of the MBL transition, where the level statistics at the MBL critical point is close to the Poisson statistics. We analyze a one-dimensional quantum spin model, in order to determine the statistical properties of the rare extensive resonances that are needed to destabilize the MBL phase. At criticality, we find that the entanglement entropy can grow with an exponent 0 < α < 1 anywhere between the area law α = 0 and the volume law α = 1 , as a function of the resonances properties, while the entanglement spectrum follows the strong multifractality statistics. In the MBL phase near criticality, we obtain the simple value ν = 1 for the correlation length exponent. Independently of the strong disorder limit, we explain why, for the many-body localization transition concerning individual eigenstates, the correlation length exponent ν is not constrained by the usual Harris inequality ν ≥ 2 / d , so that there is no theoretical inconsistency with the best numerical measure ν = 0 . 8 ( 3 ) obtained by Luitz et al. (2015).
We consider the one-dimensional lattice model of interacting fermions with disorder studied previously by Oganesyan and Huse [Phys. Rev. B 75, 155111 (2007)]. To characterize a possible many-body localization transition as a function of the disorder strength
W, we use an exact renormalization procedure in configuration space that generalizes the Aoki real-space renormalization procedure for Anderson localization one-particle models [H. Aoki, J. Phys. C 13, 3369 (1980)]. We focus on the statistical properties of the renormalized hopping
VL between two configurations separated by a distance
L in configuration space (distance being defined as the minimal number of elementary moves to go from one configuration to the other). Our numerical results point toward the existence of a many-body localization transition at a finite disorder strength
Wc. In the localized phase
W>Wc, the typical renormalized hopping ${V}_{L}^{typ}\ensuremath{\equiv}{e}^{\overline{\text{ln}\text{ }{V}_{L}}}$ decays exponentially in
L as $(\text{ln}\text{ }{V}_{L}^{typ})\ensuremath{\simeq}\ensuremath{-}\frac{L}{{\ensuremath{\xi}}_{loc}}$ and the localization length diverges as ${\ensuremath{\xi}}_{loc}(W)\ensuremath{\sim}{(W\ensuremath{-}{W}_{c})}^{\ensuremath{-}{\ensuremath{\nu}}_{loc}}$ with a critical exponent of order ${\ensuremath{\nu}}_{loc}\ensuremath{\simeq}0.45$. In the delocalized phase
W<Wc, the renormalized hopping remains a finite random variable as $L\ensuremath{\rightarrow}\ensuremath{\infty}$ and the typical asymptotic value ${V}_{\ensuremath{\infty}}^{typ}\ensuremath{\equiv}{e}^{\overline{\text{ln}\text{ }{V}_{\ensuremath{\infty}}}}$ presents an essential singularity $(\text{ln}\text{ }{V}_{\ensuremath{\infty}}^{typ})\ensuremath{\sim}\ensuremath{-}{({W}_{c}\ensuremath{-}W)}^{\ensuremath{-}\ensuremath{\kappa}}$ with an exponent of order $\ensuremath{\kappa}\ensuremath{\sim}1.4$. Finally, we show that this analysis in configuration space is compatible with the localization properties of the simplest two-point correlation function in real space.
We numerically study the possibility of many-body localization transition in a disordered quantum dimer model on the honeycomb lattice. By using the peculiar constraints of this model and state-of-the-art exact diagonalization and time evolution methods, we probe both eigenstates and dynamical properties and conclude on the existence of a localization transition, on the available time and length scales (system sizes of up to N=108 sites). We critically discuss these results and their implications.