Slow particle dynamics observed numerically even deep inside the many-body localized (MBL) regime has called the stability of MBL in the thermodynamic limit into question, and its microscopic origin remains unknown. Here, we show that this slow dynamics originates from many-body quantum beats that arise from the interaction-induced modulation of oscillations associated with single-particle hopping processes. We relate this mechanism to local integrals of motion (LIOMs) and show that it survives at arbitrarily large distances between LIOMs, is consistent with the stability of MBL in the thermodynamic limit, and is generic to many-body systems. We present new numerical results for quasiperiodic potentials and, based on the quantum beats mechanism, develop an analytical model that explains number entropy growth and quasiparticle spreading, as well as clarifies the distinct MBL phenomenology observed in quasiperiodic and random potentials. Lastly, we propose concrete signatures for observing many-body quantum beats in existing experimental platforms.
Slow particle dynamics observed numerically even deep inside the many-body localized (MBL) regime has called the stability of MBL in the thermodynamic limit into question, and its microscopic origin remains unknown. Here, we show that this slow dynamics originates from many-body quantum beats that arise from the interaction-induced modulation of oscillations associated with single-particle hopping processes. We relate this mechanism to local integrals of motion (LIOMs) and show that it survives at arbitrarily large distances between LIOMs, is consistent with the stability of MBL in the thermodynamic limit, and is generic to many-body systems. We present new numerical results for quasiperiodic potentials and, based on the quantum beats mechanism, develop an analytical model that explains number entropy growth and quasiparticle spreading, as well as clarifies the distinct MBL phenomenology observed in quasiperiodic and random potentials. Lastly, we propose concrete signatures for observing many-body quantum beats in existing experimental platforms.
We study the many-body localization (MBL) transition for interacting fermions subject to quasiperiodic potentials by constructing the local integrals of motion (LIOMs) in the MBL phase as time-averaged local operators. We study numerically how these time-averaged operators evolve across the MBL transition. We find that the norm of such time-averaged operators drops discontinuously to zero across the transition; as we discuss, this implies that LIOMs abruptly become unstable at some critical localization length of order unity. We analyze the LIOMs using hydrodynamic projections and isolating the part of the operator that is associated with interactions. Equipped with these data we perform a finite-size scaling analysis of the quasiperiodic MBL transition. Our results suggest that the quasiperiodic MBL transition occurs at considerably stronger quasiperiodic modulations and has a larger correlation-length critical exponent than previous studies had found.
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.
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.
Recent theoretical and numerical evidence suggests that localization can survive in disordered many-body systems with very high energy density, provided that interactions are sufficiently weak. Stronger interactions can destroy localization, leading to a so-called many-body localization transition. This dynamical phase transition is relevant to questions of thermalization in extended quantum systems far from the zero-temperature limit. It separates a many-body localized phase, in which localization prevents transport and thermalization, from a conducting (``ergodic'') phase in which the usual assumptions of quantum statistical mechanics hold. Here, we present numerical evidence that many-body localization also occurs in models without disorder but rather a quasiperiodic potential. In one dimension, these systems already have a single-particle localization transition, and we show that this transition becomes a many-body localization transition upon the introduction of interactions. We also comment on possible relevance of our results to experimental studies of many-body dynamics of cold atoms and nonlinear light in quasiperiodic potentials.
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.
Many-body localization (MBL) provides a mechanism by which interacting quantum systems evade thermalization, leading to persistent memory of initial conditions and slow entanglement growth. Probing these dynamical signatures in large systems and at long evolution times remains challenging for both classical simulations and current quantum devices. Here we experimentally investigate the ergodic-MBL crossover in quasiperiodic Floquet Ising systems using up to 144 qubits on an IBM Quantum processor. By implementing deep Floquet circuits reaching up to 5000 cycles, we access long-time many-body dynamics beyond the regime explored in previous quantum computing experiments. Measurements of autocorrelation functions reveal a smooth crossover from rapid thermalization at weak quasiperiodic potential strength to persistent correlations in the strong-disorder regime. Notably, in addition to the one-dimensional system, we also observe clear signatures consistent with localization behavior in the two-dimensional system. Furthermore, the quantum Fisher information exhibits logarithmic growth over thousands of Floquet cycles, providing evidence for slow entanglement spreading characteristic of the MBL regime. These results demonstrate that programmable quantum processors can serve as experimental platforms for probing nonergodic quantum many-body dynamics and exploring localization phenomena in regimes beyond the reach of classical simulations.