We use exact diagonalization to study the breakdown of many-body localization in a strongly disordered and interacting system coupled to a thermalizing environment. We show that the many-body level statistics cross over from Poisson to GOE, and the localized eigenstates thermalize, with the crossover coupling decreasing with the size of the bath in a manner consistent with the hypothesis that an infinitesimally small coupling to a thermodynamic bath should destroy localization of the eigenstates. However, signatures of incomplete localization survive in spectral functions of local operators even when the coupling to the environment is non-zero. These include a discrete spectrum and a gap at zero frequency. Both features are washed out by line broadening as one increases the coupling to the bath.
We use exact diagonalization to study the breakdown of many-body localization in a strongly disordered and interacting system coupled to a thermalizing environment. We show that the many-body level statistics cross over from Poisson to GOE, and the localized eigenstates thermalize, with the crossover coupling decreasing with the size of the bath in a manner consistent with the hypothesis that an infinitesimally small coupling to a thermodynamic bath should destroy localization of the eigenstates. However, signatures of incomplete localization survive in spectral functions of local operators even when the coupling to the environment is non-zero. These include a discrete spectrum and a gap at zero frequency. Both features are washed out by line broadening as one increases the coupling to the bath.
We use numerical exact diagonalization to analyze which aspects of the many-body localization phenomenon survive in an imperfectly isolated setting, when the system of interest is weakly coupled to a thermalizing environment. We show that widely used diagnostics (such as many-body level statistics and expectation values in exact eigenstates) cease to show signatures of many-body localization above a critical coupling that is exponentially small in the size of the environment. However, we also identify alternative diagnostics for many-body localization, in the spectral functions of local operators. Diagnostics include a discrete spectrum and a hierarchy of energy gaps, including a universal gap at zero frequency. These alternative diagnostics are shown to be robust, and continue to show signatures of many-body localization as long as the coupling to the bath is weaker than the characteristic energy scales in the system. We also examine how these signatures disappear when the coupling to the environment becomes larger than the characteristic energy scales of the system.
We study many-body-localized (MBL) systems that are weakly coupled to thermalizing environments, focusing on the spectral functions of local operators. These spectral functions carry signatures of localization even away from the limit of perfectly isolated systems. We find that, in the limit of vanishing coupling to a bath, MBL systems come in two varieties, with either discrete or continuous local spectra. Both varieties of MBL systems exhibit a ``soft gap'' at zero frequency in the spatially averaged spectral functions of local operators, which serves as a diagnostic for localization. We estimate the degree to which coupling to a bath broadens these spectral features, and we find that some characteristics of incipient localization survive as long as the system-bath coupling is much weaker than the characteristic energy scales of the system. We discuss the crossover to localization that occurs as the coupling to the external bath is tuned to zero. Since perfect isolation is impossible, we expect the ideas discussed in this paper to be relevant for experiments on many-body localization.
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, using the steady-state value of magnetization as a probe. Focusing on nonequilibrium dynamics and steady-state properties, our results indicate that within system sizes amenable to exact diagonalization, a strongly localized system interacting in a junction configuration retains remnant information on its initial state at long times despite coupling to a finite ergodic environment. In contrast, in a ladder configuration, strong dependencies on the initial state and coupling strength are observed, which can lead to either the loss or retention of information. Finally, we highlight and discuss discrepancies that can arise when similar methodologies are employed to infer localization or thermalization, revealing the need for careful interpretation.
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.
Coupling a many-body localized system to a thermal bath breaks local conservation laws and washes out signatures of localization. When the bath is nonthermal or when the system is also weakly driven, local conserved quantities acquire a highly nonthermal stationary value. We demonstrate how this property can be used to study the many-body localization phase transition in weakly open systems. Here, the strength of the coupling to the nonthermal baths plays a similar role as a finite temperature in a T=0 quantum phase transition. By tuning this parameter, we can detect key features of the many-body localization (MBL) transition: the divergence of the dynamical exponent due to Griffiths effects in one dimension and the critical disorder strength. We apply these ideas to study the MBL critical point numerically. The possibility to observe critical signatures of the MBL transition in an open system allows for new numerical approaches that overcome the limitations of exact diagonalization studies. Here, we propose a scalable numerical scheme to study the MBL critical point using matrix-product operator solution to the Lindblad equation.
Since the seminal work of Anderson, localization has been recognized as a standard mechanism allowing quantum many-body systems to escape ergodicity. This idea has acquired even more prominence in the last decade as it has been argued that localization---dubbed many-body localization (MBL) in this context---can sometimes survive local interactions in the presence of sufficiently strong disorder. A conventional signature of localization is in the statistical properties of the spectrum---spectral statistics---which differ qualitatively from those in the ergodic phase. Although features of the spectral statistics are routinely used as numerical diagnostics for localization, their derivation from first principles in the presence of nontrivial interactions has been lacking. Here we provide the example of a simple class of quantum many-body systems---which we dub strongly localized quantum circuits---that are interacting and localized and where the spectral statistics can be characterized exactly. Furthermore, we show that these systems exhibit a cascade of three different regimes for spectral correlations depending on the energy scale: at small, intermediate, and large scales they behave as disconnected patches of three decreasing sizes. We argue that these features appear in generic MBL systems, with the difference that only at the smallest scale do they become Poissonian.