We present a detailed analysis of the length-and timescales needed to approach the critical region of MBL from the delocalised phase, studying both eigenstates and the time evolution of an initial state.For the eigenstates we show that in the delocalised region there is a single length, which is a function of disorder strength, controlling the finite-size flow.Small systems look localised, and only for larger systems do resonances develop which restore ergodicity in the form of the eigenstate thermalisation hypothesis.For the transport properties, we study the time necessary to transport a single spin across a domain wall, showing how this grows quickly with increasing disorder, and compare it with the Heisenberg time.For a sufficiently large system the Heisenberg time is always larger than the transport time, but for a smaller system this is not necessarily the case.We conclude that the properties of the MBL transition cannot be explored using the system sizes or times available to current numerical and experimental studies.
We present a detailed analysis of the length-and timescales needed to approach the critical region of MBL from the delocalised phase, studying both eigenstates and the time evolution of an initial state.For the eigenstates we show that in the delocalised region there is a single length, which is a function of disorder strength, controlling the finite-size flow.Small systems look localised, and only for larger systems do resonances develop which restore ergodicity in the form of the eigenstate thermalisation hypothesis.For the transport properties, we study the time necessary to transport a single spin across a domain wall, showing how this grows quickly with increasing disorder, and compare it with the Heisenberg time.For a sufficiently large system the Heisenberg time is always larger than the transport time, but for a smaller system this is not necessarily the case.We conclude that the properties of the MBL transition cannot be explored using the system sizes or times available to current numerical and experimental studies.
Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.
"Under certain conditions, an interacting system can defy the concept of thermalization, a keystone in our understanding of physical processes. Many-body localization (MBL) is a phase of matter in which thermalization does not apply and ergodicity is broken. This striking behavior, which can appears on closed, interacting, quantum systems subject to strong disorder, has been the focus of a large body of theoretical, numerical, and experimental work in recent years. In this thesis we numerically study several aspects of MBL and the ergodic-MBL transition. Chapter 1 introduces the concept of thermalization in quantum systems, which relies on the idea of the eigenstate thermalization hypothesis (ETH). I then present localization as a phenomenon that breaks ETH, both in its single-particle as in its many-body versions. I discuss some of the main aspects that are known about MBL, as well as some of its open questions, some of which we tackle in later chapters. Chapter 2 discusses the main numerical methods used in this thesis for the study of MBL. I provide both theoretical background as well as discuss some more practical matters, such as their advantages and disadvantages, or details on the software developed/used in my studies. In Chapter 3 I present our work on MBL from the point of view of single-particle orbitals. In this work we access a complete set of approximate integrals of motion of a one-dimensional system by computing the one-particle orbitals (OPO) of highly-excited MBL energy eigenstates, which are obtained through the shift-and-invert matrix product state (SIMPS) algorithm. We then study the properties of the OPOs over large systems, up to L = 64. We find that the OPOs drawn from eigenstates at different energy densities have high overlap and their occupations are correlated with the energy of the eigenstates. Moreover, the standard deviation of the inverse participation ratio of these orbitals is maximal at the nose of the mobility edge. Also, the OPOs decay exponentially in real space, with a correlation length that increases at low disorder. In addition, we find that the probability distribution of the strength of the large-range coupling constants of the number operators generated by the OPOs approach a log-uniform distribution at strong disorder. In Chapter 4 I present our work on the hybridization of eigenstates in either the MBL or the ergodic phase, as well as at the transition. We do so by adiabatially evolving highly excited eigenstates of the Hamiltonian and measuring their their hybridization process with other eigenstates of the system. This hybridization, which dresses the eigenstate and has the potential of bringing it out of the MBL phase through the transition, is a consequence of the ""collisions"" of eigenstates in energy, which avoid level crossings every time their energy gap is small. The hybridization of eigenstates with each other involves only local regions of the system in the MBL phase, ignores locality in the ergodic phase, and is range-independent at the transition. This range independence suggests the proliferation of long-range resonances at the transition, as well as the divergence of a localization length. In Chapter 5 I present our studies on the typical and extreme (atypically strong) correlations across a one-dimensional system in the ergodic-MBL phase diagram. While typical correlations decay exponentially with range in the MBL phase, in the ergodic phase they are constant and independent of the range. Surprisingly, we identify a moderate region of the phase in which typical correlations decay as a stretched exponential with range r, and in particular as exp[-r^(1/2)] at the transition, a decay that is reminiscent of the random singlet phase. Moreover, at the transition the distribution of the logarithm of the correlations show vanishing even excess moments and non-zero range-invariant odd excess moments. This distinct behavior at the transition is in contrast with ergodic and MBL phenomenologies. In addition, we study the extreme correlations in the system. Our results suggest that strong long-range correlations proliferate at the transition, in contrast with a decay with range in the MBL phase and the lack of strong long-range correlations in the ergodic phase. Finally, we analyze the probability that a single bit of information is shared across two halves of a system, which has been proposed as a robust order parameter in the ergodic-MBL phase diagram. We find that this probability is non-zero deep in the MBL phase, but vanishes at moderate disorder, well above the transition, thus not providing a proper order parameter. Chapter 6 I summarize the work presented with an emphasis on context and perspective."
We study time dynamics of 1D disordered Heisenberg spin-1/2 chain focusing on a regime of large system sizes and a long time evolution. This regime is relevant for observation of many-body localization (MBL), a phenomenon that is expected to freeze the dynamics of the system and prevent it from reaching thermal equilibrium. Performing extensive numerical simulations of the imbalance, a quantity often employed in the experimental studies of MBL, we show that the regime of a slow power-law decay of imbalance persists to disorder strengths exceeding by at least a factor of 2 the current estimates of the critical disorder strength for MBL. Even though we investigate time evolution up to few thousands tunneling times, we observe no signs of the saturation of imbalance that would suggest freezing of system dynamics and provide a smoking gun evidence of MBL. We demonstrate that the situation is qualitatively different when the disorder is replaced by a quasiperiodic potential. In this case, we observe an emergence of a pattern of oscillations of the imbalance that is stable with respect to changes in the system size. This suggests that the dynamics of quasiperiodic systems remain fully local at the longest time scales we reach provided that the quasiperiodic potential is sufficiently strong. Our study identifies challenges in an unequivocal experimental observation of the phenomenon of MBL.
An interacting quantum system can transition from an ergodic to a many-body localized (MBL) phase under the presence of sufficiently large disorder. Both phases are radically different in their dynamical properties, which are characterized by highly excited eigenstates of the Hamiltonian. Each eigenstate can be characterized by the set of quantum numbers over the set of (local, in the MBL phase) integrals of motion of the system. In this work we study the evolution of the eigenstates of the disordered Heisenberg model as the disorder strength, W, is varied adiabatically. We focus on the probability that two `colliding' eigenstates hybridize as a function of both the range R at which they differ as well as the strength of their hybridization. We find, in the MBL phase, that the probability of a colliding eigenstate hybridizing strongly at range R decays as Pr(R)∝exp[−R/η], with a length scale η(W)=1/(Blog(W/Wc)) which diverges at the critical disorder strength Wc. This leads to range-invariance at the transition, suggesting the formation of resonating cat states at all ranges. This range invariance does not survive to the ergodic phase, where hybridization is exponentially more likely at large range, a fact that can be understood with simple combinatorial arguments. In fact, compensating for these combinatorial effects allows us to define an additional correlation length ξ in the MBL phase which is in excellent agreement with previous works and which takes the critical value 1/log(2) at the transition, found in previous works to destabilize the MBL phase. Finally, we show that deep in the MBL phase hybridization is dominated by two-level collisions of eigenstates close in energy.
The many-body localised (MBL) to thermal crossover observed in exact diagonalisation studies remains poorly understood as the accessible system sizes are too small to be in an asymptotic scaling regime. We develop a model of the crossover in short 1D chains in which the MBL phase is destabilised by the formation of many-body resonances. The model reproduces several properties of the numerically observed crossover, including an apparent correlation length exponent \nu=1 ν = 1 , exponential growth of the Thouless time with disorder strength, linear drift of the critical disorder strength with system size, scale-free resonances, apparent 1/\omega 1 / ω dependence of disorder-averaged spectral functions, and sub-thermal entanglement entropy of small subsystems. In the crossover, resonances induced by a local perturbation are rare at numerically accessible system sizes L L which are smaller than a \lambda λ . For L \gg \sqrt{\lambda} L ≫ λ (in lattice units), resonances typically overlap, and this model does not describe the asymptotic transition. The model further reproduces controversial numerical observations which Refs. claimed to be inconsistent with MBL. We thus argue that the numerics to date is consistent with a MBL phase in the thermodynamic limit.
The space of one-dimensional disordered interacting quantum models displaying a many-body localization (MBL) transition seems sufficiently rich to produce critical points with level statistics interpolating continuously between the Poisson statistics of the localized phase and the Wigner–Dyson statistics of the delocalized phase. In this paper, we consider the strong disorder limit of the MBL transition, where the level statistics at the MBL critical point is close to the Poisson statistics. We analyze a one-dimensional quantum spin model, in order to determine the statistical properties of the rare extensive resonances that are needed to destabilize the MBL phase. At criticality, we find that the entanglement entropy can grow with an exponent 0 < α < 1 anywhere between the area law α = 0 and the volume law α = 1 , as a function of the resonances properties, while the entanglement spectrum follows the strong multifractality statistics. In the MBL phase near criticality, we obtain the simple value ν = 1 for the correlation length exponent. Independently of the strong disorder limit, we explain why, for the many-body localization transition concerning individual eigenstates, the correlation length exponent ν is not constrained by the usual Harris inequality ν ≥ 2 / d , so that there is no theoretical inconsistency with the best numerical measure ν = 0 . 8 ( 3 ) obtained by Luitz et al. (2015).