Abstract The connection between entanglement dynamics and non-equilibrium statistics in isolated many-body quantum systems has been established both theoretically and experimentally. Many-body localization (MBL), a phenomenon where interacting particles in disordered (i.e. random) chains fail to thermalize, exemplifies this connection. However, the systematic proof of critical phenomena such as MBL remains challenging due to the lack of robust methods for analyzing many-body entanglement dynamics. In this paper, we identify MBL through quasi-periodic dynamics in the entanglement evolution of subsystems in a disordered Heisenberg chain. This new form of characterizing MBL, through stable quasi-periodic dynamics of entanglement-where ‘stable’ means they persist in the thermodynamic limit—concretely distinguishes between two competing scenarios: fully localized behavior of subsystems or slowly, exponentially slow in disorder, thermalizing subsystems—a heated controversy in the literature. Utilizing perturbation theory, we derive the entanglement dynamics of single spins through an infinite perturbative series, while also modeling rare Griffiths regions (locally thermal inclusions). Our results prove that in regimes of sufficiently strong disorder, the entanglement evolution of individual subsystems remains quasi-periodic in the thermodynamic limit, thereby providing concrete evidence for the stability of MBL dynamics in disordered Heisenberg chains. This behavior contrasts with the widely reported logarithmic growth of subsystem entanglement in the MBL phase. We show that the logarithmic growth observed in prior studies arises from statistical ensemble averaging, which is prohibited due to the intrinsic non-ergodic dynamics characteristic of MBL systems, rooted in their quasi-periodic features.
Abstract The connection between entanglement dynamics and non-equilibrium statistics in isolated many-body quantum systems has been established both theoretically and experimentally. Many-body localization (MBL), a phenomenon where interacting particles in disordered (i.e. random) chains fail to thermalize, exemplifies this connection. However, the systematic proof of critical phenomena such as MBL remains challenging due to the lack of robust methods for analyzing many-body entanglement dynamics. In this paper, we identify MBL through quasi-periodic dynamics in the entanglement evolution of subsystems in a disordered Heisenberg chain. This new form of characterizing MBL, through stable quasi-periodic dynamics of entanglement-where ‘stable’ means they persist in the thermodynamic limit—concretely distinguishes between two competing scenarios: fully localized behavior of subsystems or slowly, exponentially slow in disorder, thermalizing subsystems—a heated controversy in the literature. Utilizing perturbation theory, we derive the entanglement dynamics of single spins through an infinite perturbative series, while also modeling rare Griffiths regions (locally thermal inclusions). Our results prove that in regimes of sufficiently strong disorder, the entanglement evolution of individual subsystems remains quasi-periodic in the thermodynamic limit, thereby providing concrete evidence for the stability of MBL dynamics in disordered Heisenberg chains. This behavior contrasts with the widely reported logarithmic growth of subsystem entanglement in the MBL phase. We show that the logarithmic growth observed in prior studies arises from statistical ensemble averaging, which is prohibited due to the intrinsic non-ergodic dynamics characteristic of MBL systems, rooted in their quasi-periodic features.
We study time dynamics of 1D disordered Heisenberg spin-1/2 chain focusing on a regime of large system sizes and a long time evolution. This regime is relevant for observation of many-body localization (MBL), a phenomenon that is expected to freeze the dynamics of the system and prevent it from reaching thermal equilibrium. Performing extensive numerical simulations of the imbalance, a quantity often employed in the experimental studies of MBL, we show that the regime of a slow power-law decay of imbalance persists to disorder strengths exceeding by at least a factor of 2 the current estimates of the critical disorder strength for MBL. Even though we investigate time evolution up to few thousands tunneling times, we observe no signs of the saturation of imbalance that would suggest freezing of system dynamics and provide a smoking gun evidence of MBL. We demonstrate that the situation is qualitatively different when the disorder is replaced by a quasiperiodic potential. In this case, we observe an emergence of a pattern of oscillations of the imbalance that is stable with respect to changes in the system size. This suggests that the dynamics of quasiperiodic systems remain fully local at the longest time scales we reach provided that the quasiperiodic potential is sufficiently strong. Our study identifies challenges in an unequivocal experimental observation of the phenomenon of MBL.
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.
Isolated, many-body quantum systems, evolving under their intrinsic dynamics, exhibit a multitude of exotic phenomena and raise foundational questions about statistical mechanics. A flurry of theoretical work has been devoted to understanding how these systems reach thermal equilibrium in the absence of coupling to an external bath and, when thermalization does not occur, investigating the emergent non-equilibrium phases of matter. With the advent of synthetic quantum systems, such as ultra-cold atoms in optical lattices or trapped ions, these questions are no longer academic and can be directly studied in the laboratory. This dissertation explores the non-equilibrium phenomena that stem from the interplay between interactions, disorder, symmetry, topology, and external driving. First, we study how strong disorder, leading to many-body localization, can arrest the heating of a Floquet system and stabilize symmetry-protected topological order that does not have a static analogue. We analyze its dynamical and entanglement properties, highlight its duality to a discrete time crystal, and propose an experimental implementation in a cold-atom setting.Quenched disorder and the many-body localized state are crucial ingredients in protecting macroscopic quantum coherence. We explore the stability of many-body localization in two and higher dimensions and analyze its robustness to rare regions of weak disorder.We then study a second example of non-thermal behavior, namely integrability. We show that a class of random spin models, realizable in systems of atoms coupled to an optical cavity, gives rise to a rich dynamical phase diagram, which includes regions of integrability, classical chaos, and of a novel integrable structure whose conservation laws are reminiscent of the integrals of motion found in a many-body localized phase.The third group of disordered, non-ergodic systems we consider, spin glasses, have fascinating connections to complexity theory and the hardness of constraint satisfaction. We define a statistical ensemble that interpolates between the classical and quantum limits of such a problem and show that there exists a sharp boundary separating satisfiable and unsatisfiable phases.
Many-body localization is a profound phase of matter affecting the entire spectrum which emer ges in the presence of disorder in interacting many-body systems. Recently, the stability of many-body localization has been challenged by the avalanche mechanism, in which a small thermal region can spread, destabilizing localization and leading to global thermalization of the system. A key unresolved question is the critical competition between the thermal region's influence and the disorder strength required to trigger such an avalanche. Here, we numerically investigate many-body localization stability in an isolated Heisenberg spin chain of size L subjected to a disordered magnetic field. By embedding a tunable thermal region of size P, we analyze the system's behavior in both static and dynamical regimes using entanglement entropy and the gap ratio. Our study yields two main findings. Firstly, for strong disorder, the avalanche only occurs if the thermal region scales with system size, specifically when P/L exceeds a threshold value. Secondly, at strong disorder, we identify an intermediate phase between many-body localization and ergodic behavior as P increases. This intermediate phase leaves its fingerprint in both static and dynamic properties of the system and tends to vanish in the thermodynamic limit. Although our simulations are restricted to finite system sizes, the analysis suggests that these results hold in the thermodynamic limit for isolated many-body systems.
Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.
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.