We show that localization in quasiperiodically modulated, two-dimensional systems is stable to the presence of a finite density of ergodic grains. This contrasts with the case of randomly modulated systems, where such grains seed thermalizing avalanches. These results are obtained within a quantitatively accurate, self-consistent entanglement mean-field theory, which analytically describes two level systems connected to a central ergodic grain. The theory predicts the distribution of entanglement entropies of each two level system across eigenstates, and the late time values of dynamical observables. In addition to recovering the known phenomenology of avalanches, the theory reproduces exact diagonalization data, and predicts the spatial profile of the thermalized region when the avalanche fails.
We show that localization in quasiperiodically modulated, two-dimensional systems is stable to the presence of a finite density of ergodic grains. This contrasts with the case of randomly modulated systems, where such grains seed thermalizing avalanches. These results are obtained within a quantitatively accurate, self-consistent entanglement mean-field theory, which analytically describes two level systems connected to a central ergodic grain. The theory predicts the distribution of entanglement entropies of each two level system across eigenstates, and the late time values of dynamical observables. In addition to recovering the known phenomenology of avalanches, the theory reproduces exact diagonalization data, and predicts the spatial profile of the thermalized region when the avalanche fails.
We propose a multiscale diagonalization scheme to study disordered one-dimensional chains, in particular, the transition between many-body localization (MBL) and the ergodic phase, expected to be governed by resonant spots. Our scheme focuses on the dichotomy of MBL versus validity of the eigenstate thermalization hypothesis. We show that a few natural assumptions imply that the system is localized with probability one at criticality. On the ergodic side, delocalization is induced by a quantum avalanche seeded by large ergodic spots, whose size diverges at the transition. On the MBL side, the typical localization length tends to the inverse of the maximal entropy density at the transition, but there is a divergent length scale related to the response to an inclusion of large ergodic spots. A mean-field approximation analytically illustrates these results and predicts a power-law distribution for thermal inclusions at criticality.
We study bath-induced thermalization in the many-body localized XXZ spin chain subjected to a quasiperiodic magnetic field. We engineer a thermal inclusion by setting the field strength in one part of the chain below the localization threshold and analyze thermalization via the avalanche mechanism. To study the nature of such avalanches, we use two complementary observables, namely the two-point connected correlation function and the particle-number entropy. Surprisingly, the correlations alone show no signatures of avalanches. Instead, they display a logarithmic growth of the correlation length throughout the dynamics and predict localization lengths of the local integrals of motion that remain well below the avalanche threshold. In contrast, the particle number entropy shows clear signatures of the thermal avalanche, with the avalanche front progressively propagating deeper into the localized subsystem for sufficiently large baths. The discrepancy between the two observables shows that two-point correlations are unreliable in identifying thermal avalanches in quasiperiodic systems, as opposed to the random case. On one hand, qualitative results are consistent with the analytical predictions of the standard avalanche theory. On the other hand, significant quantitative deviations persist, which could be due to short-range resonances generated by the quasiperiodic potential. Our results suggest that the standard avalanche framework requires revision to account for the short-range correlations of quasiperiodic potentials.
Recent work by De Roeck et al. [Phys. Rev. B 95, 155129 (2017)] has argued that many-body localization is unstable in two and higher dimensions due to a thermalization avalanche triggered by rare regions of weak disorder. To examine these arguments, we construct several models of a finite ergodic bubble coupled to an Anderson insulator of noninteracting fermions. We first describe the ergodic region using a Gaussian orthogonal ensemble random matrix and perform an exact diagonalization study of small systems. The results are in excellent agreement with a refined theory of the thermalization avalanche that includes transient finite-size effects, lending strong support to the avalanche scenario. We then explore the limit of large system sizes by modeling the ergodic region via a Hubbard model with all-to-all random hopping: the combined system, consisting of the bubble and the insulator, can be reduced to an effective Anderson impurity problem. We find that the spectral function of a local operator in the ergodic region changes dramatically when coupling to a large number of localized fermionic states; this occurs even when the localized sites are weakly coupled to the bubble. In principle, for a given size of the ergodic region, this may arrest the avalanche. However, this back-action effect is suppressed and the avalanche can be recovered if the ergodic bubble is large enough. Thus, the main effect of the back-action is to renormalize the critical bubble size.
Many-body localized (MBL) systems lie outside the framework of statistical mechanics, as they fail to equilibrate under their own quantum dynamics. Even basic features of MBL systems, such as their stability to thermal inclusions and the nature of the dynamical transition to thermalizing behaviour, remain poorly understood. We study a simple central spin model to address these questions: a two-level system interacting with strength J with N ≫1 localized bits subject to random fields. On increasing J , the system transitions from an MBL to a delocalized phase on the vanishing scale J c ( N )∼1/ N , up to logarithmic corrections. 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 heterogeneous discontinuous transition out of the localized phase in single-site observables. We obtain these results by perturbative mapping to bond percolation on the hypercube at small J and by numerical exact diagonalization of the full many-body system. Our results support the arguments that the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions d that are high but finite. This article is part of the themed issue ‘Breakdown of ergodicity in quantum systems: from solids to synthetic matter’.
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.
Strong disorder in many-body quantum systems leads to many-body localization (MBL), manifested as long-time memory of the initial state. Spatial regions of anomalously weak disorder, arising due to fluctuations in the disordered potential, may act as ergodic inclusions that seed quantum avalanches and destabilize the MBL. To investigate this mechanism in a controlled setting, the authors plant ergodic inclusions by engineering regions of a weak disorder and compare the system's time evolution with predictions of the avalanche theory.