We construct a complete set of local integrals of motion that characterize the many-body localized (MBL) phase. Our approach relies on the assumption that local perturbations act locally on the eigenstates in the MBL phase, which is supported by numerical simulations of the random-field XXZ spin chain. We describe the structure of the eigenstates in the MBL phase and discuss the implications of local conservation laws for its nonequilibrium quantum dynamics. We argue that the many-body localization can be used to protect coherence in the system by suppressing relaxation between eigenstates with different local integrals of motion.

Certain wave functions of non-interacting quantum chaotic systems can exhibit "scars" in the fabric of their real-space density profile. Quantum scarred wave functions concentrate in the vicinity of unstable periodic classical trajectories. We introduce the notion of many-body quantum scars which reflect the existence of a subset of special many-body eigenstates concentrated in certain parts of the Hilbert space. We demonstrate the existence of scars in the Fibonacci chain-the onedimensional model with a constrained local Hilbert space realized in the 51 Rydberg atom quantum simulator [H. Bernien et al., arXiv:1707.04344]. The quantum scarred eigenstates are embedded throughout the thermalizing many-body spectrum, but surprisingly lead to direct experimental signatures such as robust oscillations following a quench from a charge-density wave state found in experiment. We develop a model based on a single particle hopping on the Hilbert space graph, which quantitatively captures the scarred wave functions up to large systems of L = 32 atoms. Our results suggest that scarred many-body bands give rise to a new universality class of quantum dynamics, which opens up opportunities for creating and manipulating novel states with long-lived coherence in systems that are now amenable to experimental study.

Recent numerical work by Bardarson et. al.[1] revealed a slow, logarithmic in time, growth of the entanglement entropy for initial product states in a putative many-body localized phase. We show that this surprising phenomenon results from the dephasing due to exponentially small interactioninduced corrections to the eigenenergies of different states. For weak interactions, we find that the entanglement entropy grows as ξ ln(V t/ ), where V is the interaction strength, and ξ is the single-particle localization length. The saturated value of the entanglement entropy at long times is determined by the participation ratios of the initial state over the eigenstates of the subsystem. Our work shows that the logarithmic entanglement growth is a universal phenomenon characteristic of the many-body localized phase in any number of spatial dimensions, and reveals a broad hierarchy of dephasing time scales present in such a phase. Introduction. While it is well-known that arbitrarily weak disorder localizes all single-particle quantummechanical states in 1 and 2 dimensions, the effect of disorder potential on the states of interacting systems largely remains an open problem. Works [2,3] conjectured that localization in a many-body system survives in the presence of weak interactions. When the strength of the interactions is increased, at some critical value a transition to the delocalized phase -a "many-body localization" transition -takes place, as observed in the numerical simulations [4][5][6][7][8][9][10][11][12][13][14].

We consider disordered many-body systems with periodic time-dependent Hamiltonians in one spatial dimension. By studying the properties of the Floquet eigenstates, we identify two distinct phases: (i) a many-body localized (MBL) phase, in which almost all eigenstates have area-law entanglement entropy, and the eigenstate thermalization hypothesis (ETH) is violated, and (ii) a delocalized phase, in which eigenstates have volume-law entanglement and obey the ETH. The MBL phase exhibits logarithmic in time growth of entanglement entropy when the system is initially prepared in a product state, which distinguishes it from the delocalized phase. We propose an effective model of the MBL phase in terms of an extensive number of emergent local integrals of motion, which naturally explains the spectral and dynamical properties of this phase. Numerical data, obtained by exact diagonalization and time-evolving block decimation methods, suggest a direct transition between the two phases.