Search

Measurement has a special role in quantum theory1: by collapsing the wavefunction it can enable phenomena such as teleportation2 and thereby alter the "arrow of time" that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space-time3-10 that go beyond established paradigms for characterizing phases, either in or out of equilibrium11-13. On present-day NISQ processors14, the experimental realization of this physics is challenging due to noise, hardware limitations, and the stochastic nature of quantum measurement. Here we address each of these experimental challenges and investigate measurement-induced quantum information phases on up to 70 superconducting qubits. By leveraging the interchangeability of space and time, we use a duality mapping9,15-17 to avoid mid-circuit measurement and access different manifestations of the underlying phases—from entanglement scaling3,4 to measurement-induced teleportation18—in a unified way. We obtain finite-size signatures of a phase transition with a decoding protocol that correlates the experimental measurement record with classical simulation data. The phases display sharply different sensitivity to noise, which we exploit to turn an inherent hardware limitation into a useful diagnostic. Our work demonstrates an approach to realize measurement-induced physics at scales that are at the limits of current NISQ processors.

Quantum many-body systems display rich phase structure in their low-temperature equilibrium states1. However, much of nature is not in thermal equilibrium. Remarkably, it was recently predicted that out-of-equilibrium systems can exhibit novel dynamical phases2–8 that may otherwise be forbidden by equilibrium thermodynamics, a paradigmatic example being the discrete time crystal (DTC)7,9–15. Concretely, dynamical phases can be defined in periodically driven many-body-localized (MBL) systems via the concept of eigenstate order7,16,17. In eigenstate-ordered MBL phases, the entire many-body spectrum exhibits quantum correlations and long-range order, with characteristic signatures in late-time dynamics from all initial states. It is, however, challenging to experimentally distinguish such stable phases from transient phenomena, or from regimes in which the dynamics of a few select states can mask typical behaviour. Here we implement tunable controlled-phase (CPHASE) gates on an array of superconducting qubits to experimentally observe an MBL-DTC and demonstrate its characteristic spatiotemporal response for generic initial states7,9,10. Our work employs a time-reversal protocol to quantify the impact of external decoherence, and leverages quantum typicality to circumvent the exponential cost of densely sampling the eigenspectrum. Furthermore, we locate the phase transition out of the DTC with an experimental finite-size analysis. These results establish a scalable approach to studying non-equilibrium phases of matter on quantum processors., A study establishes a scalable approach to engineer and characterize a many-body-localized discrete time crystal phase on a superconducting quantum processor.

We address the critical properties of the many-body localization (MBL) phase transition in one-dimensional systems subject to spatially correlated disorder. We consider a general family of disorder models, parameterized by how strong the fluctuations of the disordered couplings are when coarse-grained over a region of size $\ell$. For uncorrelated randomness, the characteristic scale for these fluctuations is $\sqrt{\ell}$; more generally they scale as $\ell^\gamma$. We discuss both positively correlated disorder ($1/2 < \gamma < 1$) and anticorrelated, or "hyperuniform," disorder ($\gamma < 1/2$). We argue that anticorrelations in the disorder are generally irrelevant at the MBL transition. Moreover, assuming the MBL transition is described by the recently developed renormalization-group scheme of Morningstar \emph{et al.} [Phys. Rev. B 102, 125134, (2020)], we argue that even positively correlated disorder leaves the critical theory unchanged, although it modifies certain properties of the many-body localized phase., Comment: 25 pages, including 9 figures, 1 table, and 3 appendices (additional reference added in v2)

We numerically study both the avalanche instability and many-body resonances in strongly-disordered spin chains exhibiting many-body localization (MBL). We distinguish between a finite-size/time MBL regime, and the asymptotic MBL phase, and identify some "landmarks" within the MBL regime. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition, in infinite-length systems, occurs much deeper in the MBL regime than has been suggested by most previous studies. Our other landmarks involve system-wide resonances. We find that the effective matrix elements producing eigenstates with system-wide resonances are enormously broadly distributed. This means that the onset of such resonances in typical samples occurs quite deep in the MBL regime, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regime in to three subregimes: (i) at strongest disorder, typical samples do not have any eigenstates that are involved in system-wide many-body resonances; (ii) there is a substantial intermediate regime where typical samples do have such resonances, but the pair of eigenstates with the minimum spectral gap does not; and (iii) in the weaker randomness regime, the minimum gap is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain non-resonant and the system thus still appears MBL in many respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i)., 23 pages, 17 figures, v2 contains updated avalanche numerics and discussion, v3 contains minor clarifications

We introduce a novel class of quantum circuits that are unitary along three distinct "arrows of time". These dynamics share some of the analytical tractability of "dual-unitary" circuits, while exhibiting distinctive and richer phenomenology. We find that two-point correlations in these dynamics are strictly confined to three directions in $(1+1)$-dimensional spacetime -- the two light cone edges, $\delta x=\pm v\delta t$, and the static worldline $\delta x=0$. Along these directions, correlation functions are obtained exactly in terms of quantum channels built from the individual gates that make up the circuit. We prove that, for a class of initial states, entanglement grows at the maximum allowed speed up to an entropy density of at least one half of the thermal value, at which point it becomes model-dependent. Finally, we extend our circuit construction to $2+1$ dimensions, where two-point correlation functions are confined to the one-dimensional edges of a tetrahedral light cone -- a subdimensional propagation of information reminiscent of "fractonic" physics., Comment: 11 pages main text + 2 page appendix, 15 figures. v2: added references, fixed typos

Unitary circuits subject to repeated projective measurements can undergo an entanglement phase transition (EPT) as a function of the measurement rate. This transition is generally understood in terms of a competition between the scrambling effects of unitary dynamics and the disentangling effects of measurements. We find that, surprisingly, EPTs are possible even in the absence of scrambling unitary dynamics, where they are best understood as arising from measurements alone. This motivates us to introduce \emph{measurement-only models}, in which the "scrambling" and "un-scrambling" effects driving the EPT are fundamentally intertwined and cannot be attributed to physically distinct processes. This represents a novel form of an EPT, conceptually distinct from that in hybrid unitary-projective circuits. We explore the entanglement phase diagrams, critical points, and quantum code properties of some of these measurement-only models. We find that the principle driving the EPTs in these models is \emph{frustration}, or mutual incompatibility, of the measurements. Suprisingly, an entangling (volume-law) phase is the generic outcome when measuring sufficiently long but still local ($\gtrsim 3$-body) operators. We identify a class of exceptions to this behavior ("bipartite ensembles") which cannot sustain an entangling phase, but display dual area-law phases, possibly with different kinds of quantum order, separated by self-dual critical points. Finally, we introduce a measure of information spreading in dynamics with measurements and use it to demonstrate the emergence of a statistical light-cone, despite the non-locality inherent to quantum measurements., 14 pages + bibliography and appendices, 10 figures. v2: added section on quantum code properties, added references. v3: substantial changes to the structure of the paper, improved clarity

Tensor Processing Units (TPUs) are specialized hardware accelerators developed by Google to support large-scale machine-learning tasks, but they can also be leveraged to accelerate and scale other linear-algebra-intensive computations. In this paper we demonstrate the usage of TPUs for massively parallel, classical simulations of quantum many-body dynamics on long timescales. We apply our methods to study the phenomenon of Floquet prethermalization, i.e., exponentially slow heating in quantum spin chains subject to high-frequency periodic driving. We simulate the dynamics of L=34 qubits for over $10^5$ Floquet periods, corresponding to circuits with millions of two-qubit gates. The circuits simulated have no additional symmetries and represent a pure-state evolution in the full $2^L$-dimensional Hilbert space. This is achieved by distributing the computation over 128 TPU cores. On that size TPU cluster, we find speedups in wall-clock runtime of 230x and 15x when compared to reference CPU and single-GPU simulations, respectively, for shorter 30-qubit simulations that can be handled by all three platforms. We study the computational cost of the simulations, as a function of both the number of qubits and the number of TPU cores used, up to our maximum capacity of L=40 qubits, which requires a ``full pod" of 2048 TPU cores with tens of terabytes of memory in total. For these simulations, an 8-TPU-core machine is comparable to a single A100 GPU, and thus the full TPU pod is comparable to a machine with hundreds of GPUs. However, the TPU pod is more energy and cost efficient, and readily accessible (via Google Cloud), unlike such large many-GPU configurations. We also study the accumulation of numerical error as a function of circuit depth in very deep circuits. Our work demonstrates that TPUs can offer significant advantages for state-of-the-art simulations of quantum many-body dynamics., Comment: 12 pages, 7 figures; v2 contains substantial improvements including GPU simulations

The extension of many-body quantum dynamics to the non-unitary domain has led to a series of exciting developments, including new out-of-equilibrium entanglement phases and phase transitions. We show how a duality transformation between space and time on one hand, and unitarity and non-unitarity on the other, can be used to realize steady state phases of non-unitary dynamics that exhibit a rich variety of behavior in their entanglement scaling with subsystem size -- from logarithmic to extensive to \emph{fractal}. We show how these outcomes in non-unitary circuits (that are "spacetime-dual" to unitary circuits) relate to the growth of entanglement in time in the corresponding unitary circuits, and how they differ, through an exact mapping to a problem of unitary evolution with boundary decoherence, in which information gets "radiated away" from one edge of the system. In spacetime-duals of chaotic unitary circuits, this mapping allows us to uncover a non-thermal volume-law entangled phase with a logarithmic correction to the entropy distinct from other known examples. Most notably, we also find novel steady state phases with \emph{fractal} entanglement scaling, $S(\ell) \sim \ell^{\alpha}$ with tunable $0 < \alpha < 1$ for subsystems of size $\ell$ in one dimension. These fractally entangled states add a qualitatively new entry to the families of many-body quantum states that have been studied as energy eigenstates or dynamical steady states, whose entropy almost always displays either area-law, volume-law or logarithmic scaling. We also present an experimental protocol for preparing these novel steady states with only a very limited amount of postselection via a type of "teleportation" between spacelike and timelike slices of quantum circuits., Comment: v2: updated interpretation of volume-law phase, added discussion on breaking unitarity. v3: published version

The dynamics of entanglement in `hybrid' non-unitary circuits (for example, involving both unitary gates and quantum measurements) has recently become an object of intense study. A major hurdle toward experimentally realizing this physics is the need to apply \emph{postselection} on random measurement outcomes in order to repeatedly prepare a given output state, resulting in an exponential overhead. We propose a method to sidestep this issue in a wide class of non-unitary circuits by taking advantage of \emph{spacetime duality}. This method maps the purification dynamics of a mixed state under non-unitary evolution onto a particular correlation function in an associated unitary circuit. This translates to an operational protocol which could be straightforwardly implemented on a digital quantum simulator. We discuss the signatures of different entanglement phases, and demonstrate examples via numerical simulations. With minor modifications, the proposed protocol allows measurement of the purity of arbitrary subsystems, which could shed light on the properties of the quantum error correcting code formed by the mixed phase in this class of hybrid dynamics., 4 pages main text + 9 page supplemental material. v2: added references, expanded supplement, minor modifications to text

The speed of information propagation is finite in quantum systems with local interactions. In many such systems, local operators spread ballistically in time and can be characterized by a "butterfly velocity", which can be measured via out-of-time-ordered correlation functions. In general, the butterfly velocity can depend asymmetrically on the direction of information propagation. In this work, we construct a family of simple 2-local Hamiltonians for understanding the asymmetric hydrodynamics of operator spreading. Our models live on a one dimensional lattice and exhibit asymmetric butterfly velocities between the left and right spatial directions. This asymmetry is transparently understood in a free (non-interacting) limit of our model Hamiltonians, where the butterfly speed can be understood in terms of quasiparticle velocities., 16 pages, 8 figures

We present a general symmetry-based framework for obtaining many-body Hamiltonians with scarred eigenstates that do not obey the eigenstate thermalization hypothesis. Our models are derived from parent Hamiltonians with a non-Abelian (or q-deformed) symmetry, whose eigenspectra are organized as degenerate multiplets that transform as irreducible representations of the symmetry (`tunnels'). We show that large classes of perturbations break the symmetry, but in a manner that preserves a particular low-entanglement multiplet of states -- thereby giving generic, thermal spectra with a `shadow' of the broken symmetry in the form of scars. The generators of the Lie algebra furnish operators with `spectrum generating algebras' that can be used to lift the degeneracy of the scar states and promote them to equally spaced `towers'. Our framework applies to several known models with scars, but we also introduce new models with scars that transform as irreducible representations of symmetries such as SU(3) and $q$-deformed SU(2), significantly generalizing the types of systems known to harbor this phenomenon. Additionally, we present new examples of generalized AKLT models with scar states that do not transform in an irreducible representation of the relevant symmetry. These are derived from parent Hamiltonians with enhanced symmetries, and bring AKLT-like models into our framework., v2- published version

Recent progress in the realm of noisy, intermediate scale quantum (NISQ) devices represents an exciting opportunity for many-body physics, by introducing new laboratory platforms with unprecedented control and measurement capabilities. We explore the implications of NISQ platforms for many-body physics in a practical sense: we ask which {\it physical phenomena}, in the domain of quantum statistical mechanics, they may realize more readily than traditional experimental platforms. As a particularly well-suited target, we identify discrete time crystals (DTCs), novel non-equilibrium states of matter that break time translation symmetry. These can only be realized in the intrinsically out-of-equilibrium setting of periodically driven quantum systems stabilized by disorder induced many-body localization. While precursors of the DTC have been observed across a variety of experimental platforms - ranging from trapped ions to nitrogen vacancy centers to NMR crystals - none have \emph{all} the necessary ingredients for realizing a fully-fledged incarnation of this phase, and for detecting its signature long-range \emph{spatiotemporal order}. We show that a new generation of quantum simulators can be programmed to realize the DTC phase and to experimentally detect its dynamical properties, a task requiring extensive capabilities for programmability, initialization and read-out. Specifically, the architecture of Google's Sycamore processor is a remarkably close match for the task at hand. We also discuss the effects of environmental decoherence, and how they can be distinguished from `internal' decoherence coming from closed-system thermalization dynamics. Already with existing technology and noise levels, we find that DTC spatiotemporal order would be observable over hundreds of periods, with parametric improvements to come as the hardware advances., v2: added appendices B and C, added Fig.1, expanded discussion. v3: added Fig. 8

We study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fracton-like conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents., 12 pages, 7 figures, 1 table; added new Appendix and additional results in v2; added new Appendix and clarified explanations in v3; published in Physical Review B

A recent experiment on a 51-atom Rydberg blockaded chain observed anomalously long-lived temporal oscillations of local observables after quenching from an antiferromagnetic initial state. This coherence is surprising as the initial state should have thermalized rapidly to infinite temperature. In this article, we show that the experimental Hamiltonian exhibits non-thermal behavior across its entire many-body spectrum, with similar finite-size scaling properties as models proximate to integrable points. Moreover, we construct an explicit small local deformation of the Hamiltonian which enhances both the signatures of integrability and the coherent oscillations observed after the quench. Our results suggest that a parent proximate integrable point controls the early-to-intermediate time dynamics of the experimental system. The unconventional quench dynamics in the parent model could signal a novel class of integrable system., Published version

We show how a finite number of conservation laws can globally `shatter' Hilbert space into exponentially many dynamically disconnected subsectors, leading to an unexpected dynamics with features reminiscent of both many body localization and quantum scars. A crisp example of this phenomenon is provided by a `fractonic' model of quantum dynamics constrained to conserve both charge and dipole moment. We show how the Hilbert space of the fractonic model dynamically fractures into disconnected emergent subsectors within a particular charge and dipole symmetry sector. This shattering can occur in arbitrary spatial dimensions. A large number of the emergent subsectors, exponentially many in system volume, have dimension one and exhibit strictly localized quantum dynamics---even in the absence of spatial disorder and in the presence of temporal noise. Other emergent subsectors display non-trivial dynamics and may be constructed by embedding finite sized non-trivial blocks into the localized subspace. While `fractonic' models provide a particularly clean realization, the shattering phenomenon is more general, as we discuss. We also discuss how the key phenomena may be readily observed in near term ultracold atom experiments. In experimental realizations, the conservation laws are approximate rather than exact, so the localization only survives up to a prethermal timescale that we estimate. We comment on the implications of these results for recent predictions of Bloch/Stark many-body localization., Comment: v2 is the published version, which combines arXiv:1904.04815 and arXiv:1910.01137 into a single publication. Text overlap with arXiv:1904.04815 expected

While a clean driven system generically absorbs energy until it reaches `infinite temperature', it may do so very slowly exhibiting what is known as a prethermal regime. Here, we show that the emergence of an additional approximately conserved quantity in a periodically driven (Floquet) system can give rise to an analogous long-lived regime. This can allow for non-trivial dynamics, even from initial states that are at a high or infinite temperature with respect to an effective Hamiltonian governing the prethermal dynamics. We present concrete settings with such a prethermal regime, one with a period-doubled (time-crystalline) response. We also present a direct diagnostic to distinguish this prethermal phenomenon from its infinitely long-lived many-body localised cousin. We apply these insights to a model of the recent NMR experiments by Rovny et al., [Phys. Rev. Lett. 120, 180603 (2018)] which, intriguingly, detected signatures of a Floquet time crystal in a clean three-dimensional material. We show that a mild but subtle variation of their driving protocol can increase the lifetime of the time-crystalline signal by orders of magnitude., Comment: v2- published version; discussions expanded and clarified

We address the hydrodynamics of operator spreading in interacting integrable lattice models. In these models, operators spread through the ballistic propagation of quasiparticles, with an operator front whose velocity is locally set by the fastest quasiparticle velocity. In interacting integrable systems, this velocity depends on the density of the other quasiparticles, so equilibrium density fluctuations cause the front to follow a biased random walk, and therefore to broaden diffusively. Ballistic front propagation and diffusive front broadening are also generically present in non-integrable systems in one dimension; thus, although the mechanisms for operator spreading are distinct in the two cases, these coarse grained measures of the operator front do not distinguish between the two cases. We present an expression for the front-broadening rate; we explicitly derive this for a particular integrable model (the "Floquet-Fredrickson-Andersen" model), and argue on kinetic grounds that it should apply generally. Our results elucidate the microscopic mechanism for diffusive corrections to ballistic transport in interacting integrable models., Comment: Published version

We consider periodically driven Anderson insulators. The short time behavior for weak, monochromatic, uniform electric fields is given by linear response theory and was famously derived by Mott. We go beyond this to consider both long times---which is the physics of Floquet late time states---and strong electric fields. This results in a `phase diagram' in the frequency-field strength plane, in which we identify four distinct regimes. These are: a linear response regime dominated by pre-existing Mott resonances, which exists provided Floquet saturation is not reached within a period; a non-linear perturbative regime, which exhibits multiphoton-absorption in response to the field; a near-adiabatic regime, which exhibits a primarily reactive response spread over the entire sample and is insensitive to pre-existing resonances; and finally an enhanced dissipative regime., Comment: 15 pages, 9 figures

We re-examine attempts to study the many-body localization transition using measures that are physically natural on the ergodic/quantum chaotic regime of the phase diagram. Using simple scaling arguments and an analysis of various models for which rigorous results are available, we find that these measures can be particularly adversely affected by the strong finite-size effects observed in nearly all numerical studies of many-body localization. This severely impacts their utility in probing the transition and the localized phase. In light of this analysis, we argue that a recent study [\v{S}untajs et al., arXiv:1905.06345] of the behavior of the Thouless energy and level repulsion in disordered spin chains likely reaches misleading conclusions, in particular as to the absence of MBL as a true phase of matter., Comment: 11 pages, 5 figures

A logarithmic signature Some one-dimensional disordered interacting quantum systems have been theoretically predicted to display a property termed many-body localization (MBL), where the system retains the memory of its initial state and fails to thermalize. However, proving experimentally that something does not occur is tricky. Instead, physicists have proposed monitoring the entanglement entropy of the system, which should grow logarithmically with evolution time in an MBL system. Lukin et al. observed this characteristic logarithmic trend in a disordered chain of interacting atoms of rubidium-87. This method should be generalizable to other experimental platforms and higher dimensions. Science , this issue p. 256

Neural network based machine learning is emerging as a powerful tool for obtaining phase diagrams when traditional regression schemes using local equilibrium order parameters are not available, as in many-body localized or topological phases. Nevertheless, instances of machine learning offering new insights have been rare up to now. Here we show that a single feed-forward neural network can decode the defining structures of two distinct MBL phases and a thermalizing phase, using entanglement spectra obtained from individual eigenstates. For this, we introduce a simplicial geometry based method for extracting multi-partite phase boundaries. We find that this method outperforms conventional metrics (like the entanglement entropy) for identifying MBL phase transitions, revealing a sharper phase boundary and shedding new insight into the topology of the phase diagram. Furthermore, the phase diagram we acquire from a single disorder configuration confirms that the machine-learning based approach we establish here can enable speedy exploration of large phase spaces that can assist with the discovery of new MBL phases. To our knowledge this work represents the first example of a machine learning approach revealing new information beyond conventional knowledge., Comment: 5 pages, 4 figures

Time crystals are proposed states of matter that spontaneously break time translation symmetry. Despite much recent interest, there is no settled definition for such states and existing definitions only tangentially refer to the well-established foundations of spontaneous symmetry breaking. We offer a definition of time crystals, which treats time translation much like any other symmetry and follows the traditional recipe for defining Wigner symmetries and order parameters. Using our definition, we find: (i) systems with time independent Hamiltonians should not exhibit time translation symmetry breaking, and (ii) the many-body localized $\ensuremath{\pi}$ spin glass/Floquet time crystal can be viewed as breaking both a global internal symmetry and the time translation symmetry, as befits the two aspects of its name.

We provide a systematic comparison of the many-body localization transition in spin chains with nonrandom quasiperiodic vs. random fields. We find evidence that these belong to two separate universality classes: one dominated by "intrinsic" intra-sample randomness, and the second dominated by external inter-sample quenched randomness. We show that the effects of inter-sample variations are strongly growing, but not yet dominant, at the system sizes probed by exact-diagonalization studies on random models. Thus, the observed finite-size critical scaling collapses in such studies appear to be in a preasymptotic regime near the nonrandom universality class, but showing signs of the initial crossover towards the quenched-disorder dominated universality class. We also show that the MBL phase is more stable for the quasiperiodic model as compared to the random one, and the transition in the quasiperiodic model suffers less from certain finite-size effects. Altogether, our results motivate the need for a greater focus on quasiperiodic models in theoretical studies of the MBL transition., v2- published version; discussions expanded

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 transition from a many-body localized phase to a thermalizing one is a dynamical quantum phase transition which lies outside the framework of equilibrium statistical mechanics. We provide a detailed study of the critical properties of this transition at finite sizes in one dimension. We find that the entanglement entropy of small subsystems looks strongly subthermal in the quantum critical regime, which indicates that it varies discontinuously across the transition as the system-size is taken to infinity, even though many other aspects of the transition look continuous. We also study the variance of the half-chain entanglement entropy which shows a peak near the transition, and find substantial variation in the entropy across eigenstates of the same sample. Further, the sample-to-sample variations in this quantity are strongly growing, and larger than the intra-sample variations. We posit that these results are consistent with a picture in which the transition to the thermal phase is driven by an eigenstate-dependent sparse resonant "backbone" of long-range entanglement, which just barely gains enough strength to thermalize the system on the thermal side of the transition as the system size is taken to infinity. This discontinuity in a global quantity --- the presence of a fully functional bath --- in turn implies a discontinuity even for local properties. We discuss how this picture compares with existing renormalization group treatments of the transition., v2 - published version, discussions clarified in various places

We introduce a unitary matrix-product operator (UMPO) based variational method that approximately finds all the eigenstates of fully many-body localized (fMBL) one-dimensional Hamiltonians. The computational cost of the variational optimization scales linearly with system size for a fixed bond dimension of the UMPO ansatz. We demonstrate the usefulness of our approach by considering the Heisenberg chain in a strongly disordered magnetic field for which we compare the approximation to exact diagonalization results., Comment: 6 pages, 4 figures

Clean and interacting periodically driven systems are believed to exhibit a single, trivial "infinite-temperature" Floquet-ergodic phase. In contrast, here we show that their disordered Floquet many-body localized counterparts can exhibit distinct ordered phases delineated by sharp transitions. Some of these are analogs of equilibrium states with broken symmetries and topological order, while others-genuinely new to the Floquet problem-are characterized by order and nontrivial periodic dynamics. We illustrate these ideas in driven spin chains with Ising symmetry.

Understanding quantum dynamics away from equilibrium is an outstanding challenge in the modern physical sciences. It is well known that out-of-equilibrium systems can display a rich array of phenomena, ranging from self-organized synchronization to dynamical phase transitions. More recently, advances in the controlled manipulation of isolated many-body systems have enabled detailed studies of non-equilibrium phases in strongly interacting quantum matter. As a particularly striking example, the interplay of periodic driving, disorder, and strong interactions has recently been predicted to result in exotic "time-crystalline" phases, which spontaneously break the discrete time-translation symmetry of the underlying drive. Here, we report the experimental observation of such discrete time-crystalline order in a driven, disordered ensemble of $\sim 10^6$ dipolar spin impurities in diamond at room-temperature. We observe long-lived temporal correlations at integer multiples of the fundamental driving period, experimentally identify the phase boundary and find that the temporal order is protected by strong interactions; this order is remarkably stable against perturbations, even in the presence of slow thermalization. Our work opens the door to exploring dynamical phases of matter and controlling interacting, disordered many-body systems., Comment: 6 + 3 pages, 4 figures

One fundamental assumption in statistical physics is that generic closed quantum many-body systems thermalize under their own dynamics. Recently, the emergence of many-body localized systems has questioned this concept, challenging our understanding of the connection between statistical physics and quantum mechanics. Here we report on the observation of a many-body localization transition between thermal and localized phases for bosons in a two-dimensional disordered optical lattice. With our single site resolved measurements we track the relaxation dynamics of an initially prepared out-of-equilibrium density pattern and find strong evidence for a diverging length scale when approaching the localization transition. Our experiments mark the first demonstration and in-depth characterization of many-body localization in a regime not accessible with state-of-the-art simulations on classical computers., Comment: 10 pages and 9 figures

The eigenstates of many-body localized (MBL) Hamiltonians exhibit low entanglement. We adapt the highly successful density-matrix renormalization group method, which is usually used to find modestly entangled ground states of local Hamiltonians, to find individual highly excited eigenstates of many body localized Hamiltonians. The adaptation builds on the distinctive spatial structure of such eigenstates. We benchmark our method against the well studied random field Heisenberg model in one dimension. At moderate to large disorder, we find that the method successfully obtains excited eigenstates with high accuracy, thereby enabling a study of MBL systems at much larger system sizes than those accessible to exact-diagonalization methods., Comment: Published version. Slightly expanded discussion; supplement added

We examine the response of a system localized by disorder to a time dependent local perturbation which varies smoothly with a characteristic timescale $\tau$. We find that such a perturbation induces a non-local response, involving a rearrangement of conserved quantities over a length scale $\sim \ln \tau$. This effect lies beyond linear response, is absent in undisordered insulators and highlights the remarkable subtlety of localized phases. The effect is common to both single particle and many body localized phases. Our results have implications for numerous fields, including topological quantum computation in quantum Hall systems, quantum control in disordered environments, and time dependent localized systems. For example, they indicate that attempts to braid quasiparticles in quantum Hall systems or Majorana nanowires will surely fail if the manipulations are performed asymptotically slowly, and thus using such platforms for topological quantum computation will require considerable engineering. They also establish that disorder localized insulators suffer from a statistical orthogonality catastrophe.

We argue that the a.c. conductivity $\sigma(\omega)$ in the many-body localized phase is a power law of frequency $\omega$ at low frequency: specifically, $\sigma(\omega) \sim \omega^\alpha$ with the exponent $\alpha$ approaching 1 at the phase transition to the thermal phase, and asymptoting to 2 deep in the localized phase. We identify two separate mechanisms giving rise to this power law: deep in the localized phase, the conductivity is dominated by rare resonant pairs of configurations; close to the transition, the dominant contributions are rare regions that are locally critical or in the thermal phase. We present numerical evidence supporting these claims, and discuss how these power laws can also be seen through polarization-decay measurements in ultracold atomic systems., Comment: 10 pages, 6 figures

We benchmark a variant of the recently introduced DMRG-X algorithm against exact results for the localized random field XX chain. We find that the eigenstates obtained via DMRG-X exhibit a highly accurate l-bit description for system sizes much bigger than the direct, many body, exact diagonalization in the spin variables is able to access. We take advantage of the underlying free fermion description of the XX model to accurately test the strengths and limitations of this algorithm for large system sizes. We discuss the theoretical constraints on the performance of the algorithm from the entanglement properties of the eigenstates, and its actual performance at different values of disorder. A small but significant improvement to the algorithm is also presented, which helps significantly with convergence. We find that at high entanglement, DMRG-X shows a bias towards eigenstates with low entanglement, but can be improved with increased bond dimension. This result suggests that one must be careful when applying the algorithm for interacting many body localized spin models near a transition., Updated ver

We consider a quantum system $A\ensuremath{\cup}B$ made up of degrees of freedom that can be partitioned into spatially disjoint regions $A$ and $B$. When the full system is in a pure state in which regions $A$ and $B$ are entangled, the quantum mechanics of region $A$ described without reference to its complement is traditionally assumed to require a reduced density matrix on $A$. While this is certainly true as an exact matter, we argue that under many interesting circumstances expectation values of typical operators anywhere inside $A$ can be computed from a suitable pure state on $A$ alone, with a controlled error. We use insights from quantum statistical mechanics---specifically the eigenstate thermalization hypothesis (ETH)---to argue for the existence of such ``representative states.''

Recent work shows that highly excited many-body localized eigenstates can exhibit broken symmetries and topological order, including in dimensions where such order would be forbidden in equilibrium. In this paper we extend this analysis to discrete symmetry protected order via the explicit examples of the Haldane phase of one dimensional spin chains and the topological Ising paramagnet in two dimensions. We comment on the challenge of extending these results to cases where the protecting symmetry is continuous.

It is now commonly believed that the ground state entanglement spectrum (ES) exhibits universal features characteristic of a given phase. In this letter, we show that this belief is false in general. Most significantly, we show that the entanglement Hamiltonian can undergo quantum phase transitions in which its ground state and low energy spectrum exhibit singular changes, even when the physical system remains in the same phase. For broken symmetry problems, this implies that the ES and the Renyi entropies can mislead entirely, while for quantum Hall systems the ES has much less universal content than assumed to date., Comment: Figured updated

A class of three-dimensional classical lattice systems with macroscopic ground state degeneracies, most famously the spin ice system, are known to exhibit ``Coulomb'' phases wherein long wavelength correlations within the ground state manifold are described by an emergent Maxwell electrodynamics. We discuss a more elaborate example of this phenomenon---the four-state Potts model on the pyrochlore lattice---where the long wavelength description now involves three independent gauge fields, as we confirm via simulation. The excitations above the ground state manifold are bions, defects that are simultaneously charged under two of the three gauge fields, and they exhibit an entropic interaction dictated by these charges. We also show that the distribution of flux loops exhibits a scaling with loop length and system size previously identified as characteristic of Coulomb phases.

The decoupling limit of a certain configuration of D3 branes in a Melvin universe defines a sector of string theory known as Puff Field Theory (PFT) - a theory with non-local dynamics but without gravity. In this work, we present a systematic analysis of the non-local states of strongly coupled PFT using gravitational holography. And we are led to a remarkable new holographic dictionary. We show that the theory admits states that may be viewed as brane protrusions from the D3 brane worldvolume. The footprint of a protrusion has finite size - the scale of non-locality in the PFT - and corresponds to an operator insertion in the PFT. We compute correlators of these states, and we demonstrate that only part of the holographic bulk is explored by this computation. We then show that the remaining space holographically encodes the dynamics of the D3 brane tentacles. The two sectors are coupled: in this holographic description, this is realized via quantum entanglement across a holographic screen - a throat in the geometry - that splits the bulk into the two regions in question. We then propose a description of PFT through a direct product of two Fock spaces - akin to other non-local settings that employ quantum group structures., 44 pages, 13 figures; v2: minor corrections, citations added; v3: typos corrected in section on local operators, some asymptotic expansions improved and made more consistent with rest of paper in section on non-local operators