The dynamical phase diagram of interacting disordered systems has seen substantial revision over the past few years. Theory must now account for a large prethermal many-body localized regime in which thermalization is extremely slow, but not completely arrested. We derive a quantitative description of these dynamics in short-ranged one-dimensional systems using a model of successive many-body resonances. The model explains the decay timescale of mean autocorrelators, the functional form of the decay-a stretched exponential-and relates the value of the stretch exponent to the broad distribution of resonance timescales. The Jacobi method of matrix diagonalization provides numerical access to this distribution, as well as a conceptual framework for our analysis. The resonance model correctly predicts the stretch exponents for several models in the literature. Successive resonances may also underlie slow thermalization in strongly disordered systems in higher dimensions, or with long-range interactions.
The dynamical phase diagram of interacting disordered systems has seen substantial revision over the past few years. Theory must now account for a large prethermal many-body localized regime in which thermalization is extremely slow, but not completely arrested. We derive a quantitative description of these dynamics in short-ranged one-dimensional systems using a model of successive many-body resonances. The model explains the decay timescale of mean autocorrelators, the functional form of the decay-a stretched exponential-and relates the value of the stretch exponent to the broad distribution of resonance timescales. The Jacobi method of matrix diagonalization provides numerical access to this distribution, as well as a conceptual framework for our analysis. The resonance model correctly predicts the stretch exponents for several models in the literature. Successive resonances may also underlie slow thermalization in strongly disordered systems in higher dimensions, or with long-range interactions.
The Kohlrausch(-Williams-Watts) law of stretched exponential relaxation has been observed for over a century and a half in diverse complex classical systems. Here we show that this law describes relaxation quite generically in closed (executing Schrödinger dynamics), interacting disordered many-body systems, using interaction range and disorder strength as primary tuning parameters. This we observe for both time-independent and periodically driven (Floquet) systems. Finite-size analysis indicates the persistence of this nonthermal relaxation regime in the thermodynamic limit thus defining a distinct dynamical regime. This regime exhibits a peak in the time scale of the long-time relaxation (following a transient), upon crossing over from weak to strong disorder. We provide a simple picture of this behavior, which naturally accounts for its generic occurrence. Formation of spin-glass — one of the possible mechanisms for stretched relaxation appears incidental to the occurrence of Kohlrausch law in our context. Finally, we provide a simple non-Hermitian Hamiltonian formulation for the dynamics of a single spin embedded in the disordered chain. This provides an analytical formula that captures not only the Kohlrausch relaxation of the disorder averaged autocorrelation but also captures the largely diverse dynamics of an arbitrary target spin in the system. Our work hence also provides a concrete quantification of the “prethermal slowness” in many-body disordered system.
Models of many-body localization (MBL) exhibit slow numerical drifts towards delocalization with increasing system size, for which no satisfactory theory exists. Numerical results indicate that these drifts are driven by the proliferation of many-body resonances at intermediate disorder strengths. We develop a statistical method to predict the distribution of resonance oscillation frequencies, accounting for how resonances formed at larger frequency scales affect the subsequent formation of resonances at lower frequencies. Working within the statistical Jacobi approximation (SJA), we derive a flow equation for a power-law exponent θ(w) characterizing the density of resonances at frequency scale w. A localized phase is described by a line of fixed points with θ(w)>0, while an instability of the flow signals resonance proliferation and the onset of thermalization. The predicted θ(w) agrees with numerical results for the Anderson model on random regular graphs and the Lévy--Rosenzweig--Porter random matrix ensemble, both of which host resonance-driven delocalization transitions. We further connect the flow to eigenstate properties such as the participation ratio and to dynamical observables such as the return probability. The predicted θ(w) also agrees with numerical measurements in real-space models of MBL at intermediate disorder strengths, providing a step towards explaining the finite-size drifts observed in MBL.
Understanding how closed quantum systems dynamically approach thermal equilibrium presents a major unresolved problem in statistical physics. Generically, non-integrable quantum systems are expected to thermalize as they comply with the Eigenstate Thermalization Hypothesis. However, in the presence of strong disorder, the dynamics can possibly slow down to a degree that systems fail to thermalize on experimentally accessible timescales, as in spin glasses or many-body localized systems. In general, particularly in long-range interacting quantum systems, the specific nature of the disorder necessary for the emergence of a prethermal, metastable state--distinctly separating the timescales of initial relaxation and subsequent slow thermalization--remains an open question. We study an ensemble of Heisenberg spins with a tunable distribution of random coupling strengths realized by a Rydberg quantum simulator. We observe a drastic change in the late-time magnetization when increasing disorder strength. The data is well described by models based on pairs of strongly interacting spins, which are treated as thermal for weak disorder and isolated for strong disorder. Our results indicate a crossover into a pair-localized prethermal regime in a closed quantum system of thousands of spins in the critical case where the exponent of the power law interaction matches the spatial dimension.
Closed, interacting, quantum systems have the potential to transition to a many-body localized (MBL) phase under the presence of sufficiently strong disorder, hence breaking ergodicity and failing to thermalize. In this work we study the distribution of correlations throughout the ergodic-MBL phase diagram. We find the typical correlations in the MBL phase decay as a stretched exponential with range r eventually crossing over to an exponential decay deep in the MBL phase. At the transition, the stretched exponential goes as e−Ar, a decay that is reminiscent of the random singlet phase. While the standard deviation of the log(QMI) has a range dependence, the log(QMI) converges to a range-invariant distribution on all other moments (i.e., the skewness and higher) at the transition. The universal nature of these distributions provides distinct phenomenology of the transition different from both the ergodic and MBL phenomenologies. In addition to the typical correlations, we study the extreme correlations in the system, finding that the probability of strong long-range correlations is maximal at the transition, suggesting the proliferation of resonances there. Finally, we analyze the probability that a single bit of information is shared across two halves of a system, finding that this probability is non-zero deep in the MBL phase but vanishes at moderate disorder well above the transition.
We prove the existence of non-equilibrium phases of matter in the prethermal regime of periodically-driven, long-range interacting systems, with power-law exponent α>d, where d is the dimensionality of the system. In this context, we predict the existence of a disorder-free, prethermal discrete time crystal in one dimension -- a phase strictly forbidden in the absence of long-range interactions. Finally, using a combination of analytic and numerical methods, we highlight key experimentally observable differences between such a prethermal time crystal and its many-body localized counterpart.
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.