Localization marks the breakdown of thermalization in subregions of quantum many-body systems in the presence of sufficiently large disorder. In this paper, we use numerical techniques to study thermalization and localization in a many-body system of coupled quantum harmonic oscillators obtained by discretizing a scalar quantum field theory in Minkowski spacetime. We consider a Gaussian initial state, constructed through a global mass quench, with a quadratic Hamiltonian, and solve for the system's exact dynamics without and with disorder in one and two spatial dimensions. We find that finite-size systems localize for sufficiently large disorder in both cases, such that the entanglement entropy of subregions retains its initial area-law behavior, and the system no longer develops long-range correlations. To probe the thermalization-to-localization transition further, we define a frequency gap ratio that measures adjacent gaps in the phase space eigenvalues of the Hamiltonian and study how it varies with disorder strength and system size. We find signatures of a chaotic regime at intermediate disorder in two spatial dimensions and argue that it is a finite-size effect, such that the system would localize for arbitrarily small disorder in the continuum in both one and two spatial dimensions, consistent with Anderson localization. Lastly, we use the frequency gap ratio to argue that in three spatial dimensions, on the other hand, the system would only localize for disorder strengths above a critical value in the continuum, again consistent with Anderson localization.
Localization marks the breakdown of thermalization in subregions of quantum many-body systems in the presence of sufficiently large disorder. In this paper, we use numerical techniques to study thermalization and localization in a many-body system of coupled quantum harmonic oscillators obtained by discretizing a scalar quantum field theory in Minkowski spacetime. We consider a Gaussian initial state, constructed through a global mass quench, with a quadratic Hamiltonian, and solve for the system's exact dynamics without and with disorder in one and two spatial dimensions. We find that finite-size systems localize for sufficiently large disorder in both cases, such that the entanglement entropy of subregions retains its initial area-law behavior, and the system no longer develops long-range correlations. To probe the thermalization-to-localization transition further, we define a frequency gap ratio that measures adjacent gaps in the phase space eigenvalues of the Hamiltonian and study how it varies with disorder strength and system size. We find signatures of a chaotic regime at intermediate disorder in two spatial dimensions and argue that it is a finite-size effect, such that the system would localize for arbitrarily small disorder in the continuum in both one and two spatial dimensions, consistent with Anderson localization. Lastly, we use the frequency gap ratio to argue that in three spatial dimensions, on the other hand, the system would only localize for disorder strengths above a critical value in the continuum, again consistent with Anderson 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.
Strong long-range interaction leads to localization in a closed quantum system without disorders. Employing the exact diagonalization method, I numerically investigate thermalization and many-body localization in translational invariant quantum chains with finite Coulomb interactions. In the computational basis, excluding all trivial degeneracies, the interaction-induced localization is well demonstrated in aspects of level statistics, eigenstate expectation values, and Anderson localization on graphs constructed of many-body bases. The nature of localization for generic eigenstates is attributed to the quasidisorder from the power-law interactions, and the full localization in the Hilbert space is similar to that in the disorder case. However, due to real-space symmetries, the long-time dynamics is dominated by the degenerated eigenstates and eventually reaches homogeneity in real space. On the other hand, the entanglement entropy exhibits size dependence beyond the area law for the same reason, even deep in the localized state, indicating an incomplete localization in real space.
It is believed that the two-dimensional (2D) Anderson model exhibits localization for any nonzero disorder in the thermodynamic limit and it is also well known that the finite-size effects are considerable in the weak disorder limit. Here we numerically study the quantum chaos to localization transition in the finite 2D Anderson model using standard indicators used in the modern literature such as the level spacing ratio, spectral form factor, variances of observable matrix elements, participation entropy, and the eigenstate entanglement entropy. We show that many features of these indicators may indicate emergence of robust single-particle quantum chaos at weak disorder. However, we argue that a careful numerical analysis is consistent with the single-parameter scaling theory and predicts the breakdown of quantum chaos at any nonzero disorder value in the thermodynamic limit. Among the hallmarks of this breakdown are the universal behavior of the spectral form factor at weak disorder, and the universal scaling of various indicators as a function of the parameter $u={(WlnV)}^{\ensuremath{-}1}$ where W is the disorder strength and V is the number of lattice sites.
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).
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.