SU(N) symmetry is incompatible with the many-body localized (MBL) phase, even when strong disorder is present. However, recent studies have shown that finite-size
SU(2) systems exhibit nonergodic, subthermal behavior, characterized by the breakdown of the eigenstate thermalization hypothesis, and by the excited eigenstates entanglement entropy that is intermediate between area and volume law. In this paper, we extend previous studies of the
SU(2)-symmetric disordered Heisenberg model to larger systems, using the time-dependent density matrix renormalization group (tDMRG) method. We simulate quench dynamics from weakly entangled initial states up to long times, finding robust subthermal behavior at stronger disorder. Although we find an increased tendency towards thermalization at larger system sizes, the subthermal regime persists at intermediate time scales, nevertheless, and therefore should be accessible experimentally. At weaker disorder, we observe signatures of thermalization; however, entanglement entropy exhibits slow sublinear growth, in contrast to conventional thermalizing systems. Furthermore, we study dynamics of the
SU(3)-symmetric disordered Heisenberg model. Similarly, strong disorder drives the system into subthermal regime, albeit thermalizing phase is broader compared to the
SU(2) case. Our findings demonstrate the robustness of the subthermal regime in spin chains with non-Abelian continuous symmetry, and are consistent with eventual thermalization at large system sizes and long time scales, suggested by previous studies.
Previous studies reveal a crucial effect of symmetries on the properties of a single particle moving in a disorder potential.More recently, a phenomenon of many-body localization (MBL) has been attracting much theoretical and experimental interest.MBL systems are characterized by the emergence of quasilocal integrals of motion and by the area-law entanglement entropy scaling of its eigenstates.In this paper, we investigate the effect of a non-Abelian SUð2Þ symmetry on the dynamical properties of a disordered Heisenberg chain.While SUð2Þ symmetry is inconsistent with conventional MBL, a new nonergodic regime is possible.In this regime, the eigenstates exhibit faster than area-law, but still strongly subthermal, scaling of the entanglement entropy.Using extensive exact diagonalization simulations, we establish that this nonergodic regime is indeed realized in the strongly disordered Heisenberg chains.We use the realspace renormalization group (RSRG) to construct approximate excited eigenstates by tree tensor networks and demonstrate the accuracy of this procedure for systems of sizes up to L ¼ 26.As the effective disorder strength is decreased, a crossover to the thermalizing phase occurs.To establish the ultimate fate of the nonergodic regime in the thermodynamic limit, we develop a novel approach for describing many-body processes that are usually neglected by the RSRG.This approach is capable of describing systems of size L ≈ 2000.We characterize the resonances that arise due to such processes, finding that they involve an ever-growing number of spins as the system size is increased.Crucially, the probability of finding resonances grows with the system's size.Even at strong disorder, we can identify a large length scale beyond which resonances proliferate.Presumably, this proliferation would eventually drive the system to a thermalizing phase.However, the extremely long thermalization timescales indicate that a broad nonergodic regime will be observable experimentally.Our study demonstrates that, similar to the case of single-particle localization, symmetries control dynamical properties of disordered, many-body systems.The approach introduced here provides a versatile tool for describing a broad range of disordered manybody systems, well beyond sizes accessible in previous studies.
Many-body localization (MBL) enables new nonergodic dynamical phases of matter that avoid the fate of thermalization. Do other ergodicity-breaking phases beyond the MBL paradigm exist? Here, the authors introduce a method for constructing and probing the stability of candidate nonergodic phases. As an example, they consider disordered spin systems, in which a non-Abelian SU(2) symmetry precludes the conventional MBL phase that is characterized by quasilocal integrals of motion and area-law entanglement entropy of eigenstates. Instead, the authors construct a candidate nonergodic phase characterized by eigenstates having logarithmic scaling of entanglement entropy. They find that the phase is only marginally unstable and, therefore, nonthermal yet non-MBL dynamics can be observed in such disordered systems.
We study the infinite temperature dynamics of a prototypical one-dimensional system expected to exhibit many-body localization. Using numerically exact methods, we establish the dynamical phase diagram of this system based on the statistics of its eigenvalues and its dynamical behavior. We show that the nonergodic phase is reentrant as a function of the interaction strength, illustrating that localization can be reinforced by sufficiently strong interactions even at infinite temperature. Surprisingly, within the accessible time range, the ergodic phase shows subdiffusive behavior, suggesting that the diffusion coefficient vanishes throughout much of the phase diagram in the thermodynamic limit. Our findings strongly suggest that Wigner-Dyson statistics of eigenvalue spacings may appear in a class of ergodic but subdiffusive systems.
We study many-body localization (MBL) in a one-dimensional system of spinless fermions with a deterministic aperiodic potential in the presence of random interactions
Vij decaying as power-law ${V}_{ij}/{({r}_{ij})}^{\ensuremath{\alpha}}$ with distance
rij. We demonstrate that MBL survives even for $\ensuremath{\alpha}<1$ and is preceded by a broad nonergodic subdiffusive phase. Starting from parameters at which the short-range interacting system shows an infinite temperature MBL phase, turning on random power-law interactions results in many-body mobility edges in the spectrum with a larger fraction of ergodic delocalized states for smaller values of $\ensuremath{\alpha}$. Hence, the critical disorder
hcr, at which ergodic to nonergodic transition takes place, increases with the range of interactions. Time evolution of the density imbalance
I(t), which has power-law decay $I(t)\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\gamma}}$ in the intermediate to large time regime, shows that the critical disorder
hcI, above which the system becomes diffusionless (with $\ensuremath{\gamma}\ensuremath{\sim}0$) and transits into the MBL phase, is much larger than
hcr. In between
hcr and
hcI there is a broad nonergodic subdiffusive phase, which is characterized by the Poissonian statistics for the level spacing ratio, multifractal eigenfunctions, and a nonzero dynamical exponent $\ensuremath{\gamma}\ensuremath{\ll}1/2$. The system continues to be subdiffusive even on the ergodic side (
h<hcr) of the MBL transition, where the eigenstates near the mobility edges are multifractal. For
h<h0<hcr, the system is superdiffusive with $\ensuremath{\gamma}>1/2$. The rich phase diagram obtained here is unique to the random nature of long-range interactions. We explain this in terms of the enhanced correlations among local energies of the effective Anderson model induced by random power-law interactions.
An important challenge in the field of many-body quantum dynamics is to identify nonergodic states of matter beyond many-body localization (MBL). Strongly disordered spin chains with non-Abelian symmetry and chains of non-Abelian anyons are natural candidates, as they are incompatible with standard MBL. In such chains, real space renormalization group methods predict a partially localized, nonergodic regime known as a quantum critical glass (a critical variant of MBL). This regime features a treelike hierarchy of integrals of motion and symmetric eigenstates with entanglement entropy that scales as a logarithmically enhanced area law. We argue that such tentative nonergodic states are perturbatively unstable using an analytic computation of the scaling of off-diagonal matrix elements and accessible level spacing of local perturbations. Our results indicate that strongly disordered chains with non-Abelian symmetry display either spontaneous symmetry breaking or ergodic thermal behavior at long times. We identify the relevant length and timescales for thermalization: Even if such chains eventually thermalize, they can exhibit nonergodic dynamics up to parametrically long timescales with a nonanalytic dependence on disorder strength.