Many-body Localized Phase Research Articles

Local resonances and parametric level dynamics in the many-body localized phase

By varying the disorder realisation in the many-body localised (MBL) phase, we investigate the influence of resonances on spectral properties. The standard theory of the MBL phase is based on the existence of local integrals of motion (LIOM), and eigenstates of the time evolution operator can be described as LIOM configurations. We show that smooth variations of the disorder give rise to avoided level crossings, and we identify these with resonances between LIOM configurations. Through this parametric approach, we develop a theory for resonances in terms of standard properties of non-resonant LIOM. This framework describes resonances that are locally pairwise, and is appropriate in arbitrarily large systems deep within the MBL phase. We show that resonances are associated with large level curvatures on paths through the ensemble of disorder realisations, and we determine the curvature distribution. By considering the level repulsion associated with resonances we calculate the two-point correlator of the level density. We also find the distributions of matrix elements of local observables and discuss implications for low-frequency dynamics.

Dynamical obstruction to localization in a disordered spin chain.

We analyze a one-dimensional XXZ spin chain in a disordered magnetic field. As the main probes of the system's behavior, we use the sensitivity of eigenstates to adiabatic transformations, as expressed through the fidelity susceptibility, in conjunction with the low-frequency asymptotes of the spectral function. We identify a region of maximal chaos-with exponentially enhanced susceptibility-which separates the many-body localized phase from the diffusive ergodic phase. This regime is characterized by slow transport, and we argue that the presence of such slow dynamics highly constrains any possible localization transition in the thermodynamic limit. Rather, the results are more consistent with absence of the localized phase.

Dynamics of Negativity of a Wannier–Stark Many‐Body Localized System Coupled to a Bath

An interacting system subjected to a strong linear potential can host a many‐body localized (MBL) phase when being slightly perturbed. This so‐called Wannier–Stark or “tilted‐field” MBL phase inherits many properties from the well‐investigated disordered MBL phase, and provides an alternative route to experimentally engineer interacting localized systems without quenched disorder. Herein, the dynamics of entanglement in a Wannier–Stark MBL system coupled to a dephasing environment is investigated. As an accessible entanglement proxy, the third Rényi negativity is used, which reduces to the third Rényi entropy in case the system is isolated from the environment. This measure captures the characteristic logarithmic growth of interacting localized phases in the intermediate‐time regime, where the effects of the coupling to the environment are not yet dominating the dynamics. Thus, it forms a tool to distinguish Wannier–Stark MBL from noninteracting Wannier–Stark localization up to intermediate time‐scales, and to quantify quantum correlations in mixed‐state dynamics.

Detecting ergodic bubbles at the crossover to many-body localization using neural networks

The transition between ergodic and many-body localized phases is expected to occur via an avalanche mechanism, in which \emph{ergodic bubbles} that arise due to local fluctuations in system properties thermalize their surroundings leading to delocalization of the system, unless the disorder is sufficiently strong to stop this process. We propose an algorithm based on neural networks that allows to detect the ergodic bubbles using experimentally measurable two-site correlation functions. Investigating time evolution of the system, we observe a logarithmic in time growth of the ergodic bubbles in the MBL regime. The distribution of the size of ergodic bubbles converges during time evolution to an exponentially decaying distribution in the MBL regime, and a power-law distribution with a thermal peak in the critical regime, supporting thus the scenario of delocalization through the avalanche mechanism. Our algorithm permits to pin-point quantitative differences in time evolution of systems with random and quasiperiodic potentials, as well as to identify rare (Griffiths) events. Our results open new pathways in studies of the mechanisms of thermalization of disordered many-body systems and beyond.

Real-time dynamics of one-dimensional and two-dimensional bosonic quantum matter deep in the many-body localized phase

Recent experiments in quantum simulators have provided evidence for the Many-Body Localized (MBL) phase in 1D and 2D bosonic quantum matter. The theoretical study of such bosonic MBL, however, is a daunting task due to the unbounded nature of its Hilbert space. In this work, we introduce a method to compute the long-time real-time evolution of 1D and 2D bosonic systems in an MBL phase at strong disorder and weak interactions. We focus on local dynamical indicators that are able to distinguish an MBL phase from an Anderson localized one. In particular, we consider the temporal fluctuations of local observables, the spatiotemporal behavior of two-time correlators and Out-Of-Time-Correlators (OTOCs). We show that these few-body observables can be computed with a computational effort that depends only polynomially on system size but is independent of the target time, by extending a recently proposed numerical method [Phys. Rev. B 99, 241114 (2019)] to mixed states and bosons. Our method also allows us to surrogate our numerical study with analytical considerations of the time-dependent behavior of the studied quantities.

Local density of states and scattering rates across the many-body localization transition

Characterizing the many-body localization (MBL) transition in strongly disordered and interacting quantum systems is an important issue in the field of condensed matter physics. We study the single particle Green's functions for a disordered interacting system in one dimension using exact diagnonalization in the infinite temperature limit and provide strong evidence that single particle excitations carry signatures of delocalization to MBL transition. In the delocalized phase, the typical values of the local density of states and the scattering rate are finite while in the MBL phase, the typical values for both the quantities become vanishingly small. The probability distribution functions of the local density of states and the scattering rate are broad log-normal distributions in the delocalized phase while the distributions become very narrow and sharply peaked close to zero in the MBL phase. We also study the eigenstate Green's function for all the many-body eigenstates and demonstrate that both, the energy resolved typical scattering rate and the typical local density of states, can track the many-body mobility edges.

Analyzing Nonequilibrium Quantum States through Snapshots with Artificial Neural Networks.

Current quantum simulation experiments are starting to explore nonequilibrium many-body dynamics in previously inaccessible regimes in terms of system sizes and timescales. Therefore, the question emerges as to which observables are best suited to study the dynamics in such quantum many-body systems. Using machine learning techniques, we investigate the dynamics and, in particular, the thermalization behavior of an interacting quantum system that undergoes a nonequilibrium phase transition from an ergodic to a many-body localized phase. We employ supervised and unsupervised training methods to distinguish nonequilibrium from equilibrium data, using the network performance as a probe for the thermalization behavior of the system. We test our methods with experimental snapshots of ultracold atoms taken with a quantum gas microscope. Our results provide a path to analyze highly entangled large-scale quantum states for system sizes where numerical calculations of conventional observables become challenging.

Interplay and competition between disorder and flat band in an interacting Creutz ladder

We clarify the interplay and competition between disorder and flat band in the Creutz ladder with interparticle interactions focusing on the system's dynamics. Without disorder, the Creutz ladder exhibits flat-band many-body localization (FMBL). In this work, we find that disorder generates drastic effects on the system, i.e., addition of it induces a thermal phase first and further increase of it leads the system to the conventional many-body-localized (MBL) phase. The competition gives novel localization properties and unconventional quench dynamics to the system. We first draw the global sketch of the localization phase diagram by focusing on the two-particle system. The thermal phase intrudes between the FMBL and MBL phases, the regime of which depends on the strength of disorder and interactions. Based on the two-particle phase diagram and the properties of the quench dynamics, we further investigate finite-filling cases in detail. At finite-filling fractions, we again find that the interplay and competition between the interactions and disorder in the original flat-band Creutz ladder induces a thermal phase, which separates the FMBL and MBL phases. We also verify that the time evolution of the system coincides with the static phase diagrams. For suitable fillings, the conservation of the initial-state information and low-growth entanglement dynamics are also observed. These properties depend on the strength of disorder and interactions.

Ergodicity breaking in an incommensurate system observed by OTOCs and loschmidt echoes: From quantum diffusion to sub-diffusion

The metal-insulator transition (MIT), which includes Anderson localization and Mott insulators as extreme regimes, has received renewed interest as the many-body effects often constitute a limitation for the handling of quantum interference. This resulted in the field dubbed many-body localization (MBL), intensively studied theoretically and experimentally as understanding the appearance of equilibration and thermalization becomes relevant in dealing with finite systems. Here, we propose a new observable to study this transition in a spin chain under the ``disorder'' of a Harper-Hofstadter-Aubry-Andr\'e on-site potential. This quantity, which we call zeroth-order gradient entanglement (ZOGE) is extracted from the fundamental Fourier mode of a family of out-of-time-ordered correlators (OTOCs). These are just Loschmidt Echoes, where a field gradient is applied before the time reversal. In the absence of many-body interactions, the ZOGE coincides with the inverse participation ratio of a Fermionic wave function. By adding an Ising interaction to an XY Hamiltonian, one can explore the MBL phase diagram of the system. Close to the critical region, the excitation dynamics is consistent with a diffusion law. However, for weak disorder, quantum diffusion prevails while for strong disorder the excitation dynamics is sub-diffusive.

Flow equations for disordered Floquet systems

In this work, we present a new approach to disordered, periodically driven (Floquet) quantum many-body systems based on flow equations. Specifically, we introduce a continuous unitary flow of Floquet operators in an extended Hilbert space, whose fixed point is both diagonal and time-independent, allowing us to directly obtain the Floquet modes. We first apply this method to a periodically driven Anderson insulator, for which it is exact, and then extend it to driven many-body localized systems within a truncated flow equation ansatz. In particular we compute the emergent Floquet local integrals of motion that characterise a periodically driven many-body localized phase. We demonstrate that the method remains well-controlled in the weakly-interacting regime, and allows us to access larger system sizes than accessible by numerically exact methods, paving the way for studies of two-dimensional driven many-body systems.

Signature of Many-Body Localization of Phonons in Strongly Disordered Superlattices.

Many-body localization (MBL) has attracted significant attention because of its immunity to thermalization, role in logarithmic entanglement entropy growth, and opportunities to reach exotic quantum orders. However, experimental realization of MBL in solid-state systems has remained challenging. Here, we report evidence of a possible phonon MBL phase in disordered GaAs/AlAs superlattices. Through grazing-incidence inelastic X-ray scattering, we observe a strong deviation of the phonon population from equilibrium in samples doped with ErAs nanodots at low temperature, signaling a departure from thermalization. This behavior occurs within finite phonon energy and wavevector windows, suggesting a localization-thermalization crossover. We support our observation by proposing a theoretical model for the effective phonon Hamiltonian in disordered superlattices, and showing that it can be mapped exactly to a disordered 1D Bose-Hubbard model with a known MBL phase. Our work provides momentum-resolved experimental evidence of phonon localization, extending the scope of MBL to disordered solid-state systems.

Scaling properties of a spatial one-particle density-matrix entropy in many-body localized systems

We investigate a spatial subsystem entropy extracted from the one-particle density matrix (OPDM) in one-dimensional disordered interacting fermions that host a many-body localized (MBL) phase. Deep in the putative MBL regime, this OPDM entropy exhibits the salient features of localization, despite not being a proper entanglement measure. We numerically show that the OPDM entropy of the eigenstates obeys an area law. Similar to the von-Neumann entropy, the OPDM entropy grows logarithmically with time after a quantum quench, albeit with a different prefactor. Both these features survive at moderately large interactions and well towards the transition into the ergodic phase. The computational cost to calculate the OPDM entropy scales only polynomially with the system size, suggesting that the OPDM provides a promising starting point for developing diagnostic tools for MBL in simulations and experiments.

Many-Body Delocalization as Symmetry Breaking.

We present a framework in which the transition between a many-body localized (MBL) phase and an ergodic one is symmetry breaking. We consider random Floquet spin chains, expressing their averaged spectral form factor (SFF) as a function of time in terms of a transfer matrix that acts in the space direction. The SFF is determined by the leading eigenvalues of this transfer matrix. In the MBL phase, the leading eigenvalue is unique, as in a symmetry-unbroken phase, while in the ergodic phase and at late times the leading eigenvalues are asymptotically degenerate, as in a system with degenerate symmetry-breaking phases. We identify the broken symmetry of the transfer matrix, introduce a local order parameter for the transition, and show that the associated correlation functions are long-ranged only in the ergodic phase.

Rare thermal bubbles at the many-body localization transition from the Fock space point of view

In this work we study the many-body localization (MBL) transition and relate it to the eigenstate structure in the Fock space. Besides the standard entanglement and multifractal probes, we introduce the radial probability distribution of eigenstate coefficients with respect to the Hamming distance in the Fock space from the wave function maximum and relate the cumulants of this distribution to the properties of the quasi-local integrals of motion in the MBL phase. We demonstrate non-self-averaging property of the many-body fractal dimension $D_q$ and directly relate it to the jump of $D_q$ as well as of the localization length of the integrals of motion at the MBL transition. We provide an example of the continuous many-body transition confirming the above relation via the self-averaging of $D_q$ in the whole range of parameters. Introducing a simple toy-model, which hosts ergodic thermal bubbles, we give analytical evidences both in standard probes and in terms of newly introduced radial probability distribution that the MBL transition in the Fock space is consistent with the avalanche mechanism for delocalization, i.e., the Kosterlitz-Thouless scenario. Thus, we show that the MBL transition can been seen as a transition between ergodic states to non-ergodic extended states and put the upper bound for the disorder scaling for the genuine Anderson localization transition with respect to the non-interacting case.

Multifractality and Fock-space localization in many-body localized states: One-particle density matrix perspective

Many-body localization (MBL) is well characterized in Fock space. To quantify the degree of this Fock-space localization, the multifractal dimension ${D}_{q}$ is employed; it has been claimed that ${D}_{q}$ shows a jump from the delocalized value ${D}_{q}=1$ in the eigenstate thermalization hypothesis (ETH) phase to a smaller value $0<{D}_{q}<1$ at the ETH-MBL transition, yet exhibiting a conspicuous discrepancy from the fully localized value ${D}_{q}=0$, which indicates that multifractality remains inside the MBL phase. Here, to better quantify the situation, we employ, instead of the commonly used computational basis, the one-particle density matrix (OPDM) and use its eigenstates (natural orbitals) as a Fock state basis for representing many-body eigenstates $|\ensuremath{\psi}\ensuremath{\rangle}$ of the system. Using this basis, we compute ${D}_{q}$ and other indices quantifying the Fock-space localization, such as the local purity $S$, which is derived from the occupation spectrum ${{n}_{\ensuremath{\alpha}}}$ (eigenvalues of the OPDM). We highlight the statistical distribution of Hamming distance ${x}_{\ensuremath{\mu}\ensuremath{\nu}}$ occurring in the pairwise coefficients $|{a}_{\ensuremath{\mu}}{|}^{2}{|{a}_{\ensuremath{\nu}}|}^{2}$ in $S$, and compare this with a related quantity considered in the literature.

Local integrals of motion and the quasiperiodic many-body localization transition

We study the many-body localization (MBL) transition for interacting fermions subject to quasiperiodic potentials by constructing the local integrals of motion (LIOMs) in the MBL phase as time-averaged local operators. We study numerically how these time-averaged operators evolve across the MBL transition. We find that the norm of such time-averaged operators drops discontinuously to zero across the transition; as we discuss, this implies that LIOMs abruptly become unstable at some critical localization length of order unity. We analyze the LIOMs using hydrodynamic projections and isolating the part of the operator that is associated with interactions. Equipped with these data we perform a finite-size scaling analysis of the quasiperiodic MBL transition. Our results suggest that the quasiperiodic MBL transition occurs at considerably stronger quasiperiodic modulations and has a larger correlation-length critical exponent than previous studies had found.

Spectral Lyapunov exponents in chaotic and localized many-body quantum systems

We consider the spectral statistics of the Floquet operator for disordered, periodically driven spin chains in their quantum chaotic and many-body localized phases (MBL). The spectral statistics are characterized by the traces of powers $t$ of the Floquet operator, and our approach hinges on the fact that, for integer $t$ in systems with local interactions, these traces can be re-expressed in terms of products of dual transfer matrices, each representing a spatial slice of the system. We focus on properties of the dual transfer matrix products as represented by a spectrum of Lyapunov exponents, which we call \textit{spectral Lyapunov exponents}. In particular, we examine the features of this spectrum that distinguish chaotic and MBL phases. The transfer matrices can be block-diagonalized using time-translation symmetry, and so the spectral Lyapunov exponents are classified according to a momentum in the time direction. For large $t$ we argue that the leading Lyapunov exponents in each momentum sector tend to zero in the chaotic phase, while they remain finite in the MBL phase. These conclusions are based on results from three complementary types of calculation. We find exact results for the chaotic phase by considering a Floquet random quantum circuit with on-site Hilbert space dimension $q$ in the large-$q$ limit. In the MBL phase, we show that the spectral Lyapunov exponents remain finite by systematically analyzing models of non-interacting systems, weakly coupled systems, and local integrals of motion. Numerically, we compute the Lyapunov exponents for a Floquet random quantum circuit and for the kicked Ising model in the two phases. As an additional result, we calculate exactly the higher point spectral form factors (hpSFF) in the large-$q$ limit, and show that the generalized Thouless time scales logarithmically in system size for all hpSFF in the large-$q$ chaotic phase.

Many-body localization: Transitions in spin models

We study the transitions between ergodic and many-body localized phases in spin systems, subject to quenched disorder, including the Heisenberg chain and the central spin model. In both cases systems with common spin lengths $1/2$ and $1$ are investigated via exact numerical diagonalization and random matrix techniques. Particular attention is paid to the sample-to-sample variance $(\Delta_sr)^2$ of the averaged consecutive-gap ratio $\langle r\rangle$ for different disorder realizations. For both types of systems and spin lengths we find a maximum in $\Delta_sr$ as a function of disorder strength, accompanied by an inflection point of $\langle r\rangle$, signaling the transition from ergodicity to many-body localization. The critical disorder strength is found to be somewhat smaller than the values reported in the recent literature. Further information about the transitions can be gained from the probability distribution of expectation values within a given disorder realization.

Unlimited growth of particle fluctuations in many-body localized phases

We study quench dynamics in a t-V chain of spinless fermions (equivalent to the spin-1/2 Heisenberg chain) with strong potential disorder. For this prototypical model of many-body localization we have recently argued that – contrary to the established picture – particles do not become fully localized. Here we summarize and expand on our previous results for various entanglement measures such as the number and the Hartley number entropy. We investigate, in particular, possible alternative interpretations of our numerical data. We find that none of these alternative interpretations appear to hold and, in the process, discover further strong evidence for the absence of localization. Furthermore, we obtain more insights into the entanglement dynamics and the particle fluctuations by comparing with non-interacting systems where we derive several strict bounds. We find that renormalized versions of these bounds also hold in the interacting case where they provide support for numerically discovered scaling relations between number and entanglement entropies.

Quantum supremacy and quantum phase transitions

Demonstrating the ability of existing quantum platforms to perform certain computational tasks intractable to classical computers represents a cornerstone in quantum computing. Despite the growing number of such proposed "quantum supreme" tasks, it remains an important challenge to identify their direct applications. In this work, we describe how the approach proposed in Ref. [arXiv:2002.11946] for demonstrating quantum supremacy in generic driven analog many-body systems, such as those found in cold atom and ion setups, can be extended to explore dynamical quantum phase transitions. We show how key quantum supremacy signatures, such as the distance between the output distribution and the expected Porter Thomas distribution at the supremacy regime, can be used as effective order parameters. We apply this approach to a periodically driven disordered 1D Ising model and show that we can accurately capture the transition between the driven thermalized and many-body localized phases. This approach also captures the transition towards the Floquet prethermalized regime for high-frequency driving. Revisiting quantum phases of matter under the light of the recent discussions about quantum supremacy draws a link between complexity theory and analog many-body systems.