No Arabic abstract
We use the stochastic series expansion quantum Monte Carlo method, together with the eigenstate-to-Hamiltonian mapping approach, to map the localized ground states of the disordered two-dimensional Heisenberg model, to excited states of a target Hamiltonian. The localized nature of the ground state is established by studying the spin stiffness, local entanglement entropy, and local magnetization. This construction allows us to define many body localized states in an energy resolved phase diagram thereby providing concrete numerical evidence for the existence of a many-body localized phase in two dimensions.
Characterizing states of matter through the lens of their ergodic properties is a fascinating new direction of research. In the quantum realm, the many-body localization (MBL) was proposed to be the paradigmatic ergodicity breaking phenomenon, which extends the concept of Anderson localization to interacting systems. At the same time, random matrix theory has established a powerful framework for characterizing the onset of quantum chaos and ergodicity (or the absence thereof) in quantum many-body systems. Here we numerically study the spectral statistics of disordered interacting spin chains, which represent prototype models expected to exhibit MBL. We study the ergodicity indicator $g=log_{10}(t_{rm H}/t_{rm Th})$, which is defined through the ratio of two characteristic many-body time scales, the Thouless time $t_{rm Th}$ and the Heisenberg time $t_{rm H}$, and hence resembles the logarithm of the dimensionless conductance introduced in the context of Anderson localization. We argue that the ergodicity breaking transition in interacting spin chains occurs when both time scales are of the same order, $t_{rm Th} approx t_{rm H}$, and $g$ becomes a system-size independent constant. Hence, the ergodicity breaking transition in many-body systems carries certain analogies with the Anderson localization transition. Intriguingly, using a Berezinskii-Kosterlitz-Thouless correlation length we observe a scaling solution of $g$ across the transition, which allows for detection of the crossing point in finite systems. We discuss the observation that scaled results in finite systems by increasing the system size exhibit a flow towards the quantum chaotic regime.
We show that the magnetization of a single `qubit spin weakly coupled to an otherwise isolated disordered spin chain exhibits periodic revivals in the localized regime, and retains an imprint of its initial magnetization at infinite time. We demonstrate that the revival rate is strongly suppressed upon adding interactions after a time scale corresponding to the onset of the dephasing that distinguishes many-body localized phases from Anderson insulators. In contrast, the ergodic phase acts as a bath for the qubit, with no revivals visible on the time scales studied. The suppression of quantum revivals of local observables provides a quantitative, experimentally observable alternative to entanglement growth as a measure of the `non-ergodic but dephasing nature of many-body localized systems.
Entanglement is usually quantified by von Neumann entropy, but its properties are much more complex than what can be expressed with a single number. We show that the three distinct dynamical phases known as thermalization, Anderson localization, and many-body localization are marked by different patterns of the spectrum of the reduced density matrix for a state evolved after a quantum quench. While the entanglement spectrum displays Poisson statistics for the case of Anderson localization, it displays universal Wigner-Dyson statistics for both the cases of many-body localization and thermalization, albeit the universal distribution is asymptotically reached within very different time scales in these two cases. We further show that the complexity of entanglement, revealed by the possibility of disentangling the state through a Metropolis-like algorithm, is signaled by whether the entanglement spectrum level spacing is Poisson or Wigner-Dyson distributed.
We show that the one-particle density matrix $rho$ can be used to characterize the interaction-driven many-body localization transition in closed fermionic systems. The natural orbitals (the eigenstates of $rho$) are localized in the many-body localized phase and spread out when one enters the delocalized phase, while the occupation spectrum (the set of eigenvalues of $rho$) reveals the distinctive Fock-space structure of the many-body eigenstates, exhibiting a step-like discontinuity in the localized phase. The associated one-particle occupation entropy is small in the localized phase and large in the delocalized phase, with diverging fluctuations at the transition. We analyze the inverse participation ratio of the natural orbitals and find that it is independent of system size in the localized phase.
We present a fully analytical description of a many body localization (MBL) transition in a microscopically defined model. Its Hamiltonian is the sum of one- and two-body operators, where both contributions obey a maximum-entropy principle and have no symmetries except hermiticity (not even particle number conservation). These two criteria paraphrase that our system is a variant of the Sachdev-Ye-Kitaev (SYK) model. We will demonstrate how this simple `zero-dimensional system displays numerous features seen in more complex realizations of MBL. Specifically, it shows a transition between an ergodic and a localized phase, and non-trivial wave function statistics indicating the presence of `non-ergodic extended states. We check our analytical description of these phenomena by parameter free comparison to high performance numerics for systems of up to $N=15$ fermions. In this way, our study becomes a testbed for concepts of high-dimensional quantum localization, previously applied to synthetic systems such as Cayley trees or random regular graphs. We believe that this is the first many body system for which an effective theory is derived and solved from first principles. The hope is that the novel analytical concepts developed in this study may become a stepping stone for the description of MBL in more complex systems.