No Arabic abstract
It is known that strong disorder in closed quantum systems leads to many-body localization (MBL), and that this quantum phase can be destroyed by coupling to an infinitely large Markovian environment. However, the stability of the MBL phase is less clear when the system and environment are of finite and comparable size. Here, we study the stability and eventual localization properties of a disordered Heisenberg spin chain coupled to a finite environment, and extensively explore the effects of environment disorder, geometry, initial state and system-bath coupling strength. Our numerical results indicate that in most cases, the system retains its localization properties despite the coupling to the finite environment, albeit to a reduced extent. However, in cases where the system and environment is strongly coupled in the ladder configuration, the eventual localization properties are highly dependent on the initial state, and could lead to either thermalization or localization.
We use low-depth quantum circuits, a specific type of tensor networks, to classify two-dimensional symmetry-protected topological many-body localized phases. For (anti-)unitary on-site symmetries we show that the (generalized) third cohomology class of the symmetry group is a topological invariant; however our approach leaves room for the existence of additional topological indices. We argue that our classification applies to quasi-periodic systems in two dimensions and systems with true random disorder within times which scale superexponentially with the inverse interaction strength. Our technique might be adapted to supply arguments suggesting the same classification for two-dimensional symmetry-protected topological ground states with a rigorous proof.
We numerically study both the avalanche instability and many-body resonances in strongly-disordered spin chains exhibiting many-body localization (MBL). We distinguish between a finite-size/time MBL regime, and the asymptotic MBL phase, and identify some landmarks within the MBL regime. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition, in infinite-length systems, occurs much deeper in the MBL regime than has been suggested by most previous studies. Our other landmarks involve system-wide resonances. We find that the effective matrix elements producing eigenstates with system-wide resonances are enormously broadly distributed. This means that the onset of such resonances in typical samples occurs quite deep in the MBL regime, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regime in to three subregimes: (i) at strongest disorder, typical samples do not have any eigenstates that are involved in system-wide many-body resonances; (ii) there is a substantial intermediate regime where typical samples do have such resonances, but the pair of eigenstates with the minimum spectral gap does not; and (iii) in the weaker randomness regime, the minimum gap is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain non-resonant and the system thus still appears MBL in many respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i).
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.
We theoretically study the response of a many-body localized system to a local quench from a quantum information perspective. We find that the local quench triggers entanglement growth throughout the whole system, giving rise to a logarithmic lightcone. This saturates the modified Lieb-Robinson bound for quantum information propagation in many-body localized systems previously conjectured based on the existence of local integrals of motion. In addition, near the localization-delocalization transition, we find that the final states after the local quench exhibit volume-law entanglement. We also show that the local quench induces a deterministic orthogonality catastrophe for highly excited eigenstates, where the typical wave-function overlap between the pre- and post-quench eigenstates decays {it exponentially} with the system size.
We introduce techniques for analysing 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 analyse 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 non-equilibrium dyanmics 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.