The nature of the dynamical quantum phase transition between the many-body localized (MBL) phase and the thermal phase remains an open question, and one line of attack on this problem is to explore this transition numerically in finite-size systems. To maximize the contrast between the MBL phase and the thermal phase in such finite-size systems, we argue one should choose a Floquet model with no local conservation laws and rapid thermalization to ``infinite temperature'' in the thermal phase. Here we introduce and explore such a Floquet spin chain model and show that standard diagnostics of the MBL-to-thermal transition behave well in this model even at modest sizes. We also introduce a physically motivated space-time correlation function, which peaks at the transition in the Floquet model, but is strongly affected by conservation laws in Hamiltonian models.
The nature of the dynamical quantum phase transition between the many-body localized (MBL) phase and the thermal phase remains an open question, and one line of attack on this problem is to explore this transition numerically in finite-size systems. To maximize the contrast between the MBL phase and the thermal phase in such finite-size systems, we argue one should choose a Floquet model with no local conservation laws and rapid thermalization to ``infinite temperature'' in the thermal phase. Here we introduce and explore such a Floquet spin chain model and show that standard diagnostics of the MBL-to-thermal transition behave well in this model even at modest sizes. We also introduce a physically motivated space-time correlation function, which peaks at the transition in the Floquet model, but is strongly affected by conservation laws in Hamiltonian models.
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.
Coupling a many-body localized system to a thermal bath breaks local conservation laws and washes out signatures of localization. When the bath is nonthermal or when the system is also weakly driven, local conserved quantities acquire a highly nonthermal stationary value. We demonstrate how this property can be used to study the many-body localization phase transition in weakly open systems. Here, the strength of the coupling to the nonthermal baths plays a similar role as a finite temperature in a T=0 quantum phase transition. By tuning this parameter, we can detect key features of the many-body localization (MBL) transition: the divergence of the dynamical exponent due to Griffiths effects in one dimension and the critical disorder strength. We apply these ideas to study the MBL critical point numerically. The possibility to observe critical signatures of the MBL transition in an open system allows for new numerical approaches that overcome the limitations of exact diagonalization studies. Here, we propose a scalable numerical scheme to study the MBL critical point using matrix-product operator solution to the Lindblad equation.
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.
Isolated quantum systems with quenched randomness exhibit many-body localization (MBL), wherein they do not reach local thermal equilibrium even when highly excited above their ground states. It is widely believed that individual eigenstates capture this breakdown of thermalization at finite size. We show that this belief is false in general and that a MBL system can exhibit the eigenstate properties of a thermalizing system. We propose that localized approximately conserved operators (l∗-bits) underlie localization in such systems. In dimensions d>1, we further argue that the existing MBL phenomenology is unstable to boundary effects and gives way to l∗-bits. Physical consequences of l∗-bits include the possibility of an eigenstate phase transition within the MBL phase unrelated to the dynamical transition in d=1 and thermal eigenstates at all parameters in d>1. Near-term experiments in ultracold atomic systems and numerics can probe the dynamics generated by boundary layers and emergence of l∗-bits.
In isolated quantum many-body systems periodically driven in time, the asymptotic dynamics at late times can exhibit distinct behavior such as thermalization or dynamical freezing. Understanding the properties of and the convergence towards infinite-time (nonequilibrium) steady states, however, remains a challenging endeavor. We propose a physically motivated Krylov space perspective on Floquet thermalization which offers a natural framework to study rates of convergence towards steady states and, simultaneously, an efficient numerical algorithm to evaluate infinite-time averages of observables within the diagonal ensemble. The effectiveness of our algorithm is demonstrated by applying it to the periodically driven mixed-field Ising model, reaching system sizes of up to 30 spins. Our method successfully resolves the transition between the ergodic and dynamically frozen phases and provides insight into the nature of the Floquet eigenstates across the phase diagram. Furthermore, we show that the long-time behavior is encoded within the localization properties of the Ritz vectors under the Floquet evolution, providing an accurate diagnostic of ergodicity.
We use exact diagonalization to explore the many-body localization transition in a random-field spin-1/2 chain. We examine the correlations within each many-body eigenstate, looking at all high-energy states and thus effectively working at infinite temperature. For weak random field the eigenstates are thermal, as expected in this nonlocalized, ``ergodic'' phase. For strong random field the eigenstates are localized with only short-range entanglement. We roughly locate the localization transition and examine some of its finite-size scaling, finding that this quantum phase transition at nonzero temperature might be showing infinite-randomness scaling with a dynamic critical exponent $z\ensuremath{\rightarrow}\ensuremath{\infty}$.