No Arabic abstract
We show how the thermodynamic properties of large many-body localized systems can be studied using quantum Monte Carlo simulations. To this end we devise a heuristic way of constructing local integrals of motion of very high quality, which are added to the Hamiltonian in conjunction with Lagrange multipliers. The ground state simulation of the shifted Hamiltonian corresponds to a high-energy state of the original Hamiltonian in case of exactly known local integrals of motion. We can show that the inevitable mixing between eigenstates as a consequence of non-perfect integrals of motion is weak enough such that the characteristics of many-body localized systems are not averaged out in our approach, unlike the standard ensembles of statistical mechanics. Our method paves the way to study higher dimensions and indicates that a full many-body localized phase in 2d, where (nearly) all eigenstates are localized, is likely to exist.
The Loschmidt echo, defined as the overlap between quantum wave function evolved with different Hamiltonians, quantifies the sensitivity of quantum dynamics to perturbations and is often used as a probe of quantum chaos. In this work we consider the behavior of the Loschmidt echo in the many body localized phase, which is characterized by emergent local integrals of motion, and provides a generic example of non-ergodic dynamics. We demonstrate that the fluctuations of the Loschmidt echo decay as a power law in time in the many-body localized phase, in contrast to the exponential decay in few-body ergodic systems. We consider the spin-echo generalization of the Loschmidt echo, and argue that the corresponding correlation function saturates to a finite value in localized systems. Slow, power-law decay of fluctuations of such spin-echo-type overlap is related to the operator spreading and is present only in the many-body localized phase, but not in a non-interacting Anderson insulator. While most of the previously considered probes of dephasing dynamics could be understood by approximating physical spin operators with local integrals of motion, the Loschmidt echo and its generalizations crucially depend on the full expansion of the physical operators via local integrals of motion operators, as well as operators which flip local integrals of motion. Hence, these probes allow to get insights into the relation between physical operators and local integrals of motion, and access the operator spreading in the many-body localized phase.
Strongly disordered systems in the many-body localized (MBL) phase can exhibit ground state order in highly excited eigenstates. The interplay between localization, symmetry, and topology has led to the characterization of a broad landscape of MBL phases ranging from spin glasses and time crystals to symmetry protected topological phases. Understanding the nature of phase transitions between these different forms of eigenstate order remains an essential open question. Here, we conjecture that no direct transition between distinct MBL orders can occur; rather, a thermal phase always intervenes. Motivated by recent advances in Rydberg-atom-based quantum simulation, we propose an experimental protocol where the intervening thermal phase can be diagnosed via the dynamics of local observables.
Recent experiments in quantum simulators have provided evidence for the Many-Body Localized (MBL) phase in 1D and 2D bosonic quantum matter. The theoretical study of such bosonic MBL, however, is a daunting task due to the unbounded nature of its Hilbert space. In this work, we introduce a method to compute the long-time real-time evolution of 1D and 2D bosonic systems in an MBL phase at strong disorder and weak interactions. We focus on local dynamical indicators that are able to distinguish an MBL phase from an Anderson localized one. In particular, we consider the temporal fluctuations of local observables, the spatiotemporal behavior of two-time correlators and Out-Of-Time-Correlators (OTOCs). We show that these few-body observables can be computed with a computational effort that depends only polynomially on system size but is independent of the target time, by extending a recently proposed numerical method [Phys. Rev. B 99, 241114 (2019)] to mixed states and bosons. Our method also allows us to surrogate our numerical study with analytical considerations of the time-dependent behavior of the studied quantities.
Using numerically exact methods we study transport in an interacting spin chain which for sufficiently strong spatially constant electric field is expected to experience Stark many-body localization. We show that starting from a generic initial state, a spin-excitation remains localized only up to a finite delocalization time, which depends exponentially on the size of the system and the strength of the electric field. This suggests that bona fide Stark many-body localization occurs only in the thermodynamic limit. We also demonstrate that the transient localization in a finite system and for electric fields stronger than the interaction strength can be well approximated by a Magnus expansion up-to times which grow with the electric field strength.
Many aspects of many-body localization (MBL), including dynamic classification of MBL phases, remain elusive. Here, by performing real-space renormalization group (RSRG) analysis we propose that there are two distinct types of MBL phases: strong MBL induced by quasiperiodic (QP) potential and weak MBL induced by random potential. Strong and weak MBL phases can be distinguished by their different probability distributions of thermal inclusion and entanglement entropy: exponential decay in strong MBL phases but power-law decay in weak MBL. We further discuss underlying mechanisms as well as experimental implications of having two distinct types of MBL phases. Strong MBL induced by QP potential may provide a more robust and promising platform for quantum information storage and processing.