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

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

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