We explore quantum dynamics in Floquet many-body systems with local conservation laws in one spatial dimension, focusing on sectors of the Hilbert space which are highly polarized. We numerically compare the predicted charge diffusion constants and quantum butterfly velocity of operator growth between models of chaotic Floquet dynamics (with discrete spacetime translation invariance) and random unitary circuits which vary both in space and time. We find that for small but nonzero density of charge (in the thermodynamic limit), the random unitary circuit correctly predicts the scaling of the butterfly velocity but incorrectly predicts the scaling of the diffusion constant. We argue that this is a consequence of quantum coherence on short time scales. Our work clarifies the settings in which random unitary circuits provide correct physical predictions for nonrandom chaotic systems and sheds light into the origin of the slow down of the butterfly effect in highly polarized systems or at low temperature.
We explore quantum dynamics in Floquet many-body systems with local conservation laws in one spatial dimension, focusing on sectors of the Hilbert space which are highly polarized. We numerically compare the predicted charge diffusion constants and quantum butterfly velocity of operator growth between models of chaotic Floquet dynamics (with discrete spacetime translation invariance) and random unitary circuits which vary both in space and time. We find that for small but nonzero density of charge (in the thermodynamic limit), the random unitary circuit correctly predicts the scaling of the butterfly velocity but incorrectly predicts the scaling of the diffusion constant. We argue that this is a consequence of quantum coherence on short time scales. Our work clarifies the settings in which random unitary circuits provide correct physical predictions for nonrandom chaotic systems and sheds light into the origin of the slow down of the butterfly effect in highly polarized systems or at low temperature.
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 slow 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) nonconserved operators.This emission leads to dissipation in the operator hydrodynamics, with the dissipative process being the slow 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 nonobservable, 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 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.
In two-dimensional Floquet systems, many-body localized dynamics in the bulk may give rise to a chaotic evolution at the one-dimensional edges that is characterized by a nonzero chiral topological index. Such anomalous dynamics is qualitatively different from local-Hamiltonian evolution. Here we show how the presence of a nonzero index affects entanglement generation and the spreading of local operators, focusing on the coarse-grained description of generic systems. We tackle this problem by analyzing exactly solvable models of random quantum cellular automata (QCA) that generalize random circuits. We find that a nonzero index leads to asymmetric butterfly velocities with different diffusive broadening of the light cones and to a modification of the order relations between the butterfly and entanglement velocities. We propose that these results can be understood via a generalization of the recently introduced entanglement membrane theory, by allowing for a spacetime entropy current, which in the case of a generic QCA is fixed by the index. We work out the implications of this current on the entanglement "membrane tension" and show that the results for random QCA are recovered by identifying the topological index with a background velocity for the coarse-grained entanglement dynamics.
Random quantum circuits are paradigmatic models of minimally structured and analytically tractable chaotic dynamics. We study a family of Floquet unitary circuits with Haar random U(1) charge conserving dynamics; the minimal such model has nearest-neighbor gates acting on spin-1/2 qubits, and a single layer of even/odd gates repeated periodically in time. We find that this minimal model is not robustly thermalizing at numerically accessible system sizes, and displays slow subdiffusive dynamics for long times. We map out the thermalization dynamics in a broader parameter space of charge conserving circuits, and understand the origin of the slow dynamics in terms of proximate localized and integrable regimes in parameter space. In contrast, we find that small extensions to the minimal model are sufficient to achieve robust thermalization; these include (i) increasing the interaction range to three-site gates, (ii) increasing the local Hilbert space dimension by appending an additional unconstrained qubit to the conserved charge on each site, or (iii) using a larger Floquet period comprised of two independent layers of gates. Our results should inform future studies of charge conserving circuits, which are relevant for a wide range of topical theoretical questions.
Understanding the spread of quantum entanglement and scrambling of information across a quantum many-body system is a fundamental problem in quantum dynamics. A necessary part of the spread of the entanglement is the growth in time of Heisenberg operators of initially localized operators. This operator growth can be diagnosed by studying the space-time structure of commutators of well-separated operators. Previous studies in this area focused on translation invariant models, but the effects of disorder, which can radically alter the motion of heat and charge, had not been studied. Here, the authors study the growth of operators in interacting systems with quenched disorder. In the localized phase, they find that operator sizes grow logarithmically in time similar to other measures of entanglement spreading. In the disordered metal, they find that operator sizes grow linearly with time, i.e., ballistically, in contrast to the diffusive motion of charge and heat. The ballistic growth of operators is quantified by the butterfly velocity which the authors relate to the charge diffusion constant and the interaction-induced inelastic scattering rate. When the diffusion of charge is slow, the resulting butterfly velocity is much smaller than the maximum speed allowed by microscopic causality constraints.
Phenomena in quantum chaotic systems such as spectral fluctuations are known to be described remarkably well by the random matrix theory (RMT) in small energy intervals and over large time scales. However, in exchange for analytical tractability, RMT has been constructed to leave out an important element of realistic physical systems { the notion of locality. Since quantum information propagates at a finite speed in a spatially-extended quantum chaotic system, there exists time scales before which RMT is effective. How do we describe the non-equilibrium quantum dynamics at intermediate timescales which are much longer than the local relaxation time but too short for RMT to apply? What are the signatures of many-body quantum chaotic systems beyond the description of RMT? In this thesis, we address these questions by analysing minimal models that incorporate a notion of locality and at the same time retain some analytical tractability of RMT. More specifically, we study Floquet quantum circuits that are random in space but periodic in time. These circuits consist of local and random quantum gates that act on a q-state `spin' at each site. After a brief review of classical chaos, quantum chaos, and quantum circuits in Chapter 1, we present three main results. In Chapter 2, we formulate a diagrammatical approach for computing ensemble averages of observables over random unitary circuits, and derive exact results on spectral statistics and dynamical observables for a Floquet circuit in the large-q limit. In Chapter 3, we further investigate the spectral statistics in spatiallyextended many-body quantum chaotic systems by computing the spectral form factor K(t) analytically and numerically. We show that K(t) follows the RMT description at times longer than the Thouless time tTh and identify the scaling of tTh at the large-q limit of a Floquet circuit. In Chapter 4, we discuss eigenstate correlations for spatially-extended many-body quantum chaotic systems, in terms of the statistical properties of matrix elements of local observables. While the eigenstate thermalization hypothesis (ETH) is known to give an excellent description of these quantities, the phenomenon of scrambling and the butterfly effect imply structure beyond ETH. We determine analytically and verify numerically the universal form of this structure at long distances and small eigenvalue separations for Floquet systems.
We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading.