Search

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

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

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

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 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

We study the scrambling of local quantum information in chaotic many-body systems in the presence of a locally conserved quantity like charge or energy that moves diffusively. The interplay between conservation laws and scrambling sheds light on the mechanism by which unitary quantum dynamics, which is reversible, gives rise to diffusive hydrodynamics, which is a dissipative process. We obtain our results in a random quantum circuit model that is constrained to have a conservation law. We find that a generic spreading operator consists of two parts: (i) a conserved part which comprises the weight of the spreading operator on the local conserved densities, whose dynamics is described by diffusive charge spreading. This conserved part also acts as a source that steadily emits a flux of (ii) non-conserved operators. This emission leads to dissipation in the operator hydrodynamics, with the dissipative process being the conversion of operator weight from local conserved operators to nonconserved, at a rate set by the local diffusion current. The emitted nonconserved parts then spread ballistically at a butterfly speed, thus becoming highly nonlocal and hence essentially non-observable, thereby acting as the "reservoir" that facilitates the dissipation. In addition, we find that the nonconserved component develops a power law tail behind its leading ballistic front due to the slow dynamics of the conserved components. This implies that the out-of-time-order commutator (OTOC) between two initially separated operators grows sharply upon the arrival of the ballistic front but, in contrast to systems with no conservation laws, it develops a diffusive tail and approaches its asymptotic late-time value only as a power of time instead of exponentially. We also derive these results within an effective hydrodynamic description which contains multiple coupled diffusion equations., Published version; Figure added, title modified, discussion on "observable entropy" added

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

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.

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

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.

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.

Periodic driving has emerged as a powerful tool in the quest to engineer new and exotic quantum phases. While driven many-body systems are generically expected to absorb energy indefinitely and reach an infinite-temperature state, the rate of heating can be exponentially suppressed when the drive frequency is large compared to the local energy scales of the system -- leading to long-lived 'prethermal' regimes. In this work, we experimentally study a bosonic cloud of ultracold atoms in a driven optical lattice and identify such a prethermal regime in the Bose-Hubbard model. By measuring the energy absorption of the cloud as the driving frequency is increased, we observe an exponential-in-frequency reduction of the heating rate persisting over more than 2 orders of magnitude. The tunability of the lattice potentials allows us to explore one- and two-dimensional systems in a range of different interacting regimes. Alongside the exponential decrease, the dependence of the heating rate on the frequency displays features characteristic of the phase diagram of the Bose-Hubbard model, whose understanding is additionally supported by numerical simulations in one dimension. Our results show experimental evidence of the phenomenon of Floquet prethermalization, and provide insight into the characterization of heating for driven bosonic systems.