No Arabic abstract
We enquire into the quasi-many-body localization in topologically ordered states of matter, revolving around the case of Kitaev toric code on ladder geometry, where different types of anyonic defects carry different masses induced by environmental errors. Our study verifies that random arrangement of anyons generates a complex energy landscape solely through braiding statistics, which suffices to suppress the diffusion of defects in such multi-component anyonic liquid. This non-ergodic dynamic suggests a promising scenario for investigation of quasi-many-body localization. Computing standard diagnostics evidences that, in such disorder-free many-body system, a typical initial inhomogeneity of anyons gives birth to a glassy dynamics with an exponentially diverging time scale of the full relaxation. A by-product of this dynamical effect is manifested by the slow growth of entanglement entropy, with characteristic time scales bearing resemblance to those of inhomogeneity relaxation. This setting provides a new platform which paves the way toward impeding logical errors by self-localization of anyons in a generic, high energy state, originated in their exotic statistics.
It is typically assumed that disorder is essential to realize Anderson localization. Recently, a number of proposals have suggested that an interacting, translation invariant system can also exhibit localization. We examine these claims in the context of a one-dimensional spin ladder. At intermediate time scales, we find slow growth of entanglement entropy consistent with the phenomenology of many-body localization. However, at longer times, all finite wavelength spin polarizations decay in a finite time, independent of system size. We identify a single length scale which parametrically controls both the eventual spin transport times and the divergence of the susceptibility to spin glass ordering. We dub this long pre-thermal dynamical behavior, intermediate between full localization and diffusion, quasi-many body localization.
We construct a family of many-body wave functions to study the many-body localization phase transition. The wave functions have a Rokhsar-Kivelson form, in which the weight for the configurations are chosen from the Gibbs weights of a classical spin glass model, known as the Random Energy Model, multiplied by a random sign structure to represent a highly excited state. These wave functions show a phase transition into an MBL phase. In addition, we see three regimes of entanglement scaling with subsystem size: scaling with entanglement corresponding to an infinite temperature thermal phase, constant scaling, and a sub-extensive scaling between these limits. Near the phase transition point, the fluctuations of the Renyi entropies are non-Gaussian. We find that Renyi entropies with different Renyi index transition into the MBL phase at different points and have different scaling behavior, suggesting a multifractal behavior.
The presence of flat bands is a source of localization in lattice systems. While flat bands are often unstable with respect to interactions between the particles, they can persist in certain cases. We consider a diamond ladder with transverse hopping that possesses such stable flat bands and show that many-body localization appears in the presence of interactions. We demonstrate that the eigenstate thermalization hypothesis is violated and verify localization by time evolution of local observables, revival probabilities, and participation ratios. Thus, this system appears to be an example for many-body localization without disorder.
We theoretically study correlations present deep in the spectrum of many-body-localized systems. An exact analytical expression for the spectral form factor of Poisson spectra can be obtained and is shown to agree well with numerical results on two models exhibiting many-body-localization: a disordered quantum spin chain and a phenomenological $l$-bit model based on the existence of local integrals of motion. We also identify a universal regime that is insensitive to the global density of states as well as spectral edge effects.
Polynomially filtered exact diagonalization method (POLFED) for large sparse matrices is introduced. The algorithm finds an optimal basis of a subspace spanned by eigenvectors with eigenvalues close to a specified energy target by a spectral transformation using a high order polynomial of the matrix. The memory requirements scale better with system size than in the state-of-the-art shift-invert approach. The potential of POLFED is demonstrated examining many-body localization transition in 1D interacting quantum spin-1/2 chains. We investigate the disorder strength and system size scaling of Thouless time. System size dependence of bipartite entanglement entropy and of the gap ratio highlights the importance of finite-size effects in the system. We discuss possible scenarios regarding the many-body localization transition obtaining estimates for the critical disorder strength.