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.
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.
Many-body localization provides a generic mechanism of ergodicity breaking in quantum systems. In contrast to conventional ergodic systems, many-body-localized (MBL) systems are characterized by extensively many local integrals of motion (LIOM), which underlie the absence of transport and thermalization in these systems. Here we report a physically motivated construction of local integrals of motion in the MBL phase. We show that any local operator (e.g., a local particle number or a spin-flip operator), evolved with the system's Hamiltonian and averaged over time, becomes a LIOM in the MBL phase. Such operators have a clear physical meaning, describing the response of the MBL system to a local perturbation. In particular, when a local operator represents a density of some globally conserved quantity, the corresponding LIOM describes how this conserved quantity propagates through the MBL phase. Being uniquely defined and experimentally measurable, these LIOMs provide a natural tool for characterizing the properties of the MBL phase, in both experiments and numerical simulations. We demonstrate the latter by numerically constructing an extensive set of LIOMs in the MBL phase of a disordered spin-chain model. We show that the resulting LIOMs are quasilocal and use their decay to extract the localization length and establish the location of the transition between the MBL and ergodic phases.
The presence and character of local integrals of motion---quasilocal operators that commute with the Hamiltonian---encode valuable information about the dynamics of a quantum system. In particular, strongly disordered many-body systems can generically avoid thermalization when there are extensively many such operators. In this work, we explicitly construct local conserved operators in one-dimensional spin chains by directly minimizing their commutator with the Hamiltonian. We demonstrate the existence of an extensively large set of local integrals of motion in the many-body localized phase of the disordered XXZ spin chain. These operators are shown to have exponentially decaying tails, in contrast to the ergodic phase where the decay is (at best) polynomial in the size of the subsystem. We study the algebraic properties of localized operators and confirm that in the many-body localized phase, they are well described by ``dressed'' spin operators.
Many-body localized (MBL) systems are often described using their local integrals of motion, which, for spin systems, are commonly assumed to be a local unitary transform of the set of on-site spin-z operators.We show that this assumption cannot hold for topologically ordered MBL systems.Using a suitable definition to capture such systems in any spatial dimension, we demonstrate a number of features, including that MBL topological order, if present, (i) is the same for all eigenstates, (ii) is robust in character against any perturbation preserving MBL, and (iii) implies that on topologically nontrivial manifolds a complete set of integrals of motion must include nonlocal ones in the form of local-unitary-dressed noncontractible Wilson loops.Our approach is well suited for tensor-network methods and is expected to allow these to resolve highly excited finite-size-split topological eigenspaces despite their overlap in energy.We illustrate our approach on the disordered Kitaev chain, toric code, and X-cube model.
Conservation laws are of paramount importance in our understanding of classical and quantum dynamics. Here, we present a general framework for constructing exact quantum integrals of motion with the desired locality and quantum numbers, which will be illustrated for the case of many-body-localization (MBL). The latter has been understood theoretically in terms of the existence of an extensive number of local integrals of motion (LIOMs). Using our approach, we show that one can express the task of finding LIOMs as an optimization problem. For some specifications, this problem surprisingly connects to the question of finding classical ground states of spin-glass models. Our work unifies previous results obtained in the MBL context and reveals intriguing connections between many-body localization, spin-glass physics, and constrained optimization problems.
We introduce techniques for analyzing the structure of quantum states of many-body localized (MBL) spin chains by identifying correlation clusters from pairwise correlations. These techniques proceed by interpreting pairwise correlations in the state as a weighted graph, which we analyze using an established graph theoretic clustering algorithm. We validate our approach by studying the eigenstates of a disordered XXZ spin chain across the MBL to ergodic transition, as well as the nonequilibrium dynamics in the MBL phase following a global quantum quench. We successfully reproduce theoretical predictions about the MBL transition obtained from renormalization group schemes. Furthermore, we identify a clear signature of many-body dynamics analogous to the logarithmic growth of entanglement. The techniques that we introduce are computationally inexpensive and, in combination with matrix product state methods, allow for the study of large-scale localized systems. Moreover, the correlation functions we use are directly accessible in a range of experimental settings, including cold atoms.
Isolated quantum systems with quenched randomness exhibit many-body localization (MBL), wherein they do not reach local thermal equilibrium even when highly excited above their ground states. It is widely believed that individual eigenstates capture this breakdown of thermalization at finite size. We show that this belief is false in general and that a MBL system can exhibit the eigenstate properties of a thermalizing system. We propose that localized approximately conserved operators (l∗-bits) underlie localization in such systems. In dimensions d>1, we further argue that the existing MBL phenomenology is unstable to boundary effects and gives way to l∗-bits. Physical consequences of l∗-bits include the possibility of an eigenstate phase transition within the MBL phase unrelated to the dynamical transition in d=1 and thermal eigenstates at all parameters in d>1. Near-term experiments in ultracold atomic systems and numerics can probe the dynamics generated by boundary layers and emergence of l∗-bits.