We study the transitions between ergodic and many-body localized phases in spin systems, subject to quenched disorder, including the Heisenberg chain and the central spin model. In both cases systems with common spin lengths
1/2 and 1 are investigated via exact numerical diagonalization and random matrix techniques. Particular attention is paid to the sample-to-sample variance ${({\mathrm{\ensuremath{\Delta}}}_{s}r)}^{2}$ of the averaged consecutive-gap ratio $\ensuremath{\langle}r\ensuremath{\rangle}$ for different disorder realizations. For both types of systems and spin lengths we find a maximum in ${\mathrm{\ensuremath{\Delta}}}_{s}r$ as a function of disorder strength, accompanied by an inflection point of $\ensuremath{\langle}r\ensuremath{\rangle}$, signaling the transition from ergodicity to many-body localization. The critical disorder strength is found to be somewhat smaller than the values reported in the recent literature. Further information about the transitions can be gained from the probability distribution of expectation values within a given disorder realization.
We analyze the finite-size scaling of the average gap ratio and the entanglement entropy across the many-body localization (MBL) transition in one dimensional Heisenberg spin chain with quasiperiodic (QP) potential. By using the recently introduced cost-function approach, we compare different scenarios for the transition using exact diagonalization of systems up to 22 lattice sites. Our findings suggest that the MBL transition in the QP Heisenberg chain belongs to the class of the Berezinskii-Kosterlitz-Thouless transition, the same as in the case of uniformly disordered systems as advocated in recent studies. Moreover, we observe that the critical disorder strength shows a clear sublinear drift with the system size as compared to the linear drift seen in random disordered models, suggesting that the finite-size effects in the MBL transition for the QP systems are less severe than that in the random disordered scenario. Moreover, deep in the ergodic regime, we find an unexpected double-peak structure of distribution of on-site magnetizations that can be traced back to the strong correlations present in the QP potential.
We study the many-body localization problem in the non-Abelian SU(2)-invariant random antiferromagnetic exchange model in 1D. Exact and sparse matrix diagonalization methods are used to calculate eigenvalues and eigenvectors of the Hamiltonian matrix. We investigate the behavior of the energy level gap-ratio statistic, participation ratio, entanglement entropy and the entanglement spectral parameter as a function of disorder strength. Different distributions of random couplings are considered. We find, up to L = 24 , a clear distinction between our non-Abelian model and the more often studied random field Heisenberg model: the regime of seemingly localized behavior is much less pronounced in the random exchange model than in the field model case.
Many-body localization occurs in isolated quantum systems when Anderson localization persists in the presence of finite interactions. Despite strong evidence for the existence of a many-body localization transition, a reliable extraction of the critical disorder strength is difficult due to a large drift with system size in the studied quantities. In this Letter, we explore two entanglement properties that are promising for the study of the many-body localization transition: the variance of the half-chain entanglement entropy of exact eigenstates and the long time change in entanglement after a local quench from an exact eigenstate. We investigate these quantities in a disordered quantum Ising chain and use them to estimate the critical disorder strength and its energy dependence. In addition, we analyze a spin-glass transition at large disorder strength and provide evidence for it being a separate transition. We, thereby, give numerical support for a recently proposed phase diagram of many-body localization with localization protected quantum order [Huse et al., Phys. Rev. B 88, 014206 (2013).
We investigate the transition from the many-body localized phase to the ergodic thermalized phase at an infinite temperature in an
XY spin chain with
L spins, which experiences power-law decaying interactions in the form of ${V}_{ij}\ensuremath{\propto}1/{\left|i\ensuremath{-}j\right|}^{\ensuremath{\alpha}}$ (
i,j=1,...,L) and a random transverse field. By performing large-scale exact diagonalization for the chain size up to
L=18, we systematically analyze the energy gap statistics, the half-chain entanglement entropy, and the uncertainty of the entanglement entropy of the system at different interaction exponents $\ensuremath{\alpha}$. The finite-size critical scaling allows us to determine the critical disorder strength
Wc and the critical exponent $\ensuremath{\nu}$ at the many-body-localization phase transition, as a function of the interaction exponent $\ensuremath{\alpha}$ in the limit $L\ensuremath{\rightarrow}\ensuremath{\infty}$. We find that both
Wc and $\ensuremath{\nu}$ diverge when $\ensuremath{\alpha}$ decreases to a critical power of ${\ensuremath{\alpha}}_{c}\ensuremath{\simeq}1.16\ifmmode\pm\else\textpm\fi{}0.17$, indicating the absence of many-body localization for $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$. Our result differs from the analytical prediction of ${\ensuremath{\alpha}}_{c}=3/2$ found in a previous work based on a scaling argument [A. L. Burin, Phys. Rev. B 92, 104428 (2015)], but is close to the recent numerical result ${\ensuremath{\alpha}}_{c}\ensuremath{\approx}1$ obtained from quantum dynamics simulation [A. Safavi-Naini, M. L. Wall, O. L. Acevedo, A. M. Rey, and R. M. Nandkishore, Phys. Rev. A 99, 033610 (2019)].
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.