arXiv:1805.08258 · Physical review. B./Physical review. B
Exploration of the existence of a distinct quasi many-body localized phase: Numerical study of a translationally invariant system in the thermodynamic limit
We consider a quench in an infinite spin ladder describing a system with two species of bosons in the limit of strong interactions. If the heavy bosonic species has infinite mass, the model becomes a spin chain with quenched binary disorder which shows true Anderson localization (AL) or many-body localization (MBL). For finite hopping amplitude ${J}^{\ensuremath{'}}$ of the heavy particles, on the other hand, we find an exponential polarization decay with a relaxation rate which depends monotonically on ${J}^{\ensuremath{'}}$. Furthermore, the entanglement entropy changes from a constant (AL) or logarithmic (MBL) scaling in time $t$ for ${J}^{\ensuremath{'}}=0$ to a subballistic power law, ${S}_{\mathrm{ent}}\ensuremath{\sim}{t}^{\ensuremath{\alpha}}$ with $\ensuremath{\alpha}<1$, for finite ${J}^{\ensuremath{'}}$. We do not find a distinct regime in time where the dynamics for ${J}^{\ensuremath{'}}\ensuremath{\ne}0$ shows the characteristics of an MBL phase. Instead, we discover a time regime with distinct dephasing and entanglement times, leading to dynamics in this regime which is different both from a localized and from a fully ergodic phase.
We consider a quench in an infinite spin ladder describing a system with two species of bosons in the limit of strong interactions. If the heavy bosonic species has infinite mass, the model becomes a spin chain with quenched binary disorder which shows true Anderson localization (AL) or many-body localization (MBL). For finite hopping amplitude J\ensuremath′ of the heavy particles, on the other hand, we find an exponential polarization decay with a relaxation rate which depends monotonically on J\ensuremath′. Furthermore, the entanglement entropy changes from a constant (AL) or logarithmic (MBL) scaling in time t for J\ensuremath′=0 to a subballistic power law, Sent\ensuremath∼t\ensuremathα with \ensuremathα<1, for finite J\ensuremath′. We do not find a distinct regime in time where the dynamics for J\ensuremath′\ensuremath=0 shows the characteristics of an MBL phase. Instead, we discover a time regime with distinct dephasing and entanglement times, leading to dynamics in this regime which is different both from a localized and from a fully ergodic phase.
We investigate the phase transition between an ergodic and a many-body localized phase in infinite anisotropic spin-1/2 Heisenberg chains with binary disorder. Starting from the N\'eel state, we analyze the decay of antiferromagnetic order ms(t) and the growth of entanglement entropy Sent(t) during unitary time evolution. Near the phase transition we find that ms(t) decays exponentially to its asymptotic value ${m}_{s}(\ensuremath{\infty})\ensuremath{\ne}0$ in the localized phase while the data are consistent with a power-law decay at long times in the ergodic phase. In the localized phase, ${m}_{s}(\ensuremath{\infty})$ shows an exponential sensitivity on disorder with a critical exponent $\ensuremath{\nu}\ensuremath{\sim}0.9$. The entanglement entropy in the ergodic phase grows subballistically, ${S}_{\mathrm{ent}}(t)\ensuremath{\sim}{t}^{\ensuremath{\alpha}}$, $\ensuremath{\alpha}\ensuremath{\le}1$, with $\ensuremath{\alpha}$ varying continuously as a function of disorder. Exact diagonalizations for small systems, on the other hand, do not show a clear scaling with system size and attempts to determine the phase boundary from these data seem to overestimate the extent of the ergodic phase.
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 localized (MBL) systems are characterized by the absence of transport and thermalization and, therefore, cannot be described by conventional statistical mechanics. In this paper, using analytic arguments and numerical simulations, we study the behavior of local observables in an isolated MBL system following a quantum quench. For the case of a global quench, we find that the local observables reach stationary, highly nonthermal values at long times as a result of slow dephasing characteristic of the MBL phase. These stationary values retain the local memory of the initial state due to the existence of local integrals of motion in the MBL phase. The temporal fluctuations around stationary values exhibit universal power-law decay in time, with an exponent set by the localization length and the diagonal entropy of the initial state. Such a power-law decay holds for any local observable and is related to the logarithmic in time growth of entanglement in the MBL phase. This behavior distinguishes the MBL phase from both the Anderson insulator (where no stationary state is reached) and from the ergodic phase (where relaxation is expected to be exponential). For the case of a local quench, we also find a power-law approach of local observables to their stationary values when the system is prepared in a mixed state. Quench protocols considered in this paper can be naturally implemented in systems of ultracold atoms in disordered optical lattices, and the behavior of local observables provides a direct experimental signature of many-body localization.
Many-body localization of a disorder interacting boson system in one dimension is studied numerically at the filling factor being one-half, in terms of level statistics, local compressibility, correlation function and entanglement entropies. The von Neumann entanglement entropy is decoupled into a particle number entropy and a configuration entropy. An anomalous volume-law behavior is found for the configuration entanglement entropy to confirm a recent experimental observation [A. Lukin, M. Rispoli, R. Schittko, et al., Science 364, 256 (2019)] for sufficient strong disorder, while the particle number entropy fulfills an area-law corresponding to the total entropy for disordered spin chain. The localization length are extracted from a two-body correlation function for many-body localization states and corresponding time-evolutions states as well. A phase diagrams is established with consisting of an ergodic thermalized region and a many-body-localization region in a parameter space of the disorder strength and the energy density. Two regions are separated by a many-body mobility edge deducted from the standard deviation of the particle-number entanglement entropy, which appears consistent with that based on the localization length. Slow dynamics characterized by a logarithmic time-dependence is explicitly shown for both the particle number entropy and the configuration entropy in an intermediate regime of their time-evolutions, which does not show up in the Anderson localization case, i.e. non-interacting disorder systems.
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.
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)].