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, 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.
We propose a Gaussian disorder to study many-body localization (MBL) in one-dimensional quantum system and coupled system-environment models. Unlike uniform disorder, Gaussian disorder is controlled by the standard deviation σ (driving MBL transitions) and the mean μ (causing global energy shifts). By calculating the von Neumann entropy, mean gap ratio, and dynamical indicators, we elucidate the ergodic-to-MBL transition as σ , with Gaussian disorder exhibiting more pronounced finite-size critical point drift effects compared to conventional random disorder. Through finite-size scaling analysis, it shows that there are no significant differences between the two disorder types at our assumed system sizes. We construct two system-environment coupling configurations (ladder and staggered), showing that under weak coupling, the system and environment evolve independently, whereas strong coupling induces cooperative localization. The ladder configuration shows a strong dependence on the initial state, which may lead to either the loss or retention of information, while no such dependence is observed in the staggered configuration. This can be explained by the changes in energy density of the initial state caused by variations in different model structures. When treating the system-environment as a unified two-dimensional ladder, asymmetric disorder on the two chains can drive a global MBL phase transition, with slightly different behaviors compared to symmetric disorder. Additionally, we emphasize the necessity of correctly selecting the initial state when characterizing the phase transition through dynamic indicators.
We investigate the stability of the many-body localized phase against quantum avalanche instabilities in a one-dimensional Heisenberg spin chain with long-range power-law interactions (V∝r−α). By combining exact diagonalization of static properties with Lindblad master equation simulations of open-system dynamics, we systematically map the interplay between interaction range and disorder strength. Our finite-size scaling analysis of entanglement entropy identifies a critical interaction exponent αc≈2, which separates a fragile regime, characterized by an exponentially diverging critical disorder, from a robust short-range regime. To rigorously test the system's resistance to avalanches, we couple the boundary to an infinite-temperature bath and track the propagation of the thermalization front into the localized bulk. We find that the characteristic thermalization time follows a unified scaling law, Trth∼exp[κ(α)LW] (herein, L is the system size, and W is the disorder intensity), which diverges exponentially with the product of system size and disorder strength. This suppression enables the derivation of a quantitative stability criterion, Wstab(α), representing the minimum critical disorder strength required to maintain avalanche stability. Our results confirm that the MBL phase remains asymptotically stable in the thermodynamic limit when disorder exceeds an interaction-dependent threshold, bridging theoretical debates on long-range MBL and providing a roadmap for observing these dynamics in experimental platforms such as Rydberg atom arrays.
Closed generic quantum many-body systems may fail to thermalize under certain conditions even after long times, a phenomenon called many-body localization (MBL). Numerous studies have left no doubt about the stability of the MBL phase in strongly disordered one-dimensional systems. However, the situation is much less clear when a small part of the system is ergodic, a scenario which also has important implications for the existence of many-body localization in higher dimensions. Here we address this question experimentally using a large-scale quantum simulator of ultracold bosons in a two-dimensional optical lattice. We prepare two-component mixtures of varying relative population and implement a disorder potential which is only experienced by one of the components. The second non-disordered clean component plays the role of a small quantum bath of adjustable size that is collisionally coupled to the dirty component. Our experiments provide evidence that MBL is remarkably stable in 2D and can even survive in the presence of a small bath for our observed long timescales. A too large bath, however, ultimately leads to the destruction of MBL.
This paper introduces Gaussian disorder, characterized by two parameters:the expected value and the standard deviation.Studying this type of disorder enhances our understanding of how many-body localization (MBL) transition is influenced by the central value and breadth of disorder.We apply Gaussian disorder to a one-dimensional spin-1/2 chain to reveal its features.Our findings indicate that the crucial parameter driving localization transitions is the standard deviation; as it increases, the system transitions from the ergodic phase to the localized phase, while the expected value has no effect,so the critical factor influencing MBL transition is the breadth of disorder rather than its central value.We then explore Gaussian disorder in two system-environment models (hereafter referred to as the sye-env models, consisting of ladder configuration and staggered configuration),analyzing the mutual influence when they are coupled.We extensively explore the impacts brought about by system disorder,environment disorder,initial state,and the strength of the sys-env interaction.Finally,for theses two different models,we consider the system and environment as a whole,studying how the localization properties of the entire system are influenced by the disorder strength of each chain and the interaction strength between the two chains.In the ladder configuration, we observe a strong dependence on the initial state,which can lead to either the loss or retention of information,whereas no such dependence on the initial state is found in the staggered configuration.We emphasize that when using dynamic indicators to analyze system characteristics, it is essential to consider the influence of different initial states,as this may lead to significant differences under certain structural conditions,and not all initial states can accurately reflect the localization properties of the many-body system.
Despite the exponentially large amount of information required in the quantum description of many-body systems, finite size numerical simulations have paved the way for recent progress in the understanding of quantum phases of matter both in and out-of-equilibrium. \n \nIn this spirit, using exact diagonalization for a small number of spins ~ 20, first, we investigate the fate of both ergodic quantum many-body systems that thermalize at long times and systems with quenched disorder that exhibit many-body localization (MBL) and fail to thermalize (due to an extensive set of local conserved degrees of freedom) under periodic driving. Ergodic systems always delocalize in energy space and heat up to infinite temperature, for both local driving, in which the time dependence of the Hamiltonian is restricted to a few contiguous degrees of freedom, and global driving. On the other hand, MBL systems remain localized at finite energy density for local driving, while in the case of global driving these systems can either heat up to infinite temperature or remain localized depending on the parameters of the drive. Underlying the latter is the emergence of an effective time independent MBL Hamiltonian describing the long time dynamics of the system. These numerical predictions have recently been verified experimentally in a one dimensional quasi-disordered optical lattice of interacting fermions wherein unitary time evolution can be experimentally probed. \n \nSubsequently, we study the stability of MBL systems with respect to another perturbation, in this case by thermal inclusions, which remains poorly understood. We consider a simple model to address this question: a two level system interacting with 𝑁 localized bits subject to random fields. In the large 𝑁 limit, the localized bits are thermalized for any nonzero ratio of interactions to disorder strength. Moreover, at finite 𝑁, in the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either "on" (strongly entangled) or "off" (weakly entangled) in eigenstates. The clusters of "on" bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogenous discontinuous transition out of the localized phase in single site observables. These results also imply the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions 𝑑 that are high but finite. \n \nThe second part of this thesis is motivated by the recent breakthroughs of artificial intelligence in domains such as computer vision involving highly dimensional data of complexity not unlike the one found in the description of many-body systems. Using Monte Carlo simulations to generate synthetic data, we find that kernel methods, a class of machine learning algorithms, can learn the mathematical \nform of physical discriminators associated with phase transitions, such as order parameters and Hamiltonian constraints, for a set of two dimensional spin models: the ferromagnetic Ising model, a conserved-order-parameter Ising model, and the Ising gauge theory. \n \nOur work shows that the combination of numerical simulations with the ability of machine learning algorithms for automatic feature detection hold promise to become standard tools in the exploration and characterization of different phases of matter in condensed matter and other many-body systems.