Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth.We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries.We study both 1 þ 1D and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems.In 1 þ 1D, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed v B .We find that in 1 þ 1D, the "front" of the OTOC broadens diffusively, with a width scaling in time as t 1=2 .We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization.Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem.This implies that the width of the front of the averaged OTOC scales as t 1=3 in 2 þ 1D and as t 0.240 in 3 þ 1D (exponents of the Kardar-Parisi-Zhang universality class).We support our analytic argument with simulations in 2 þ 1D.We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical.However, when full spatial rotational symmetry is present in 2 þ 1D, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution.For an alternative perspective on the OTOC in 1 þ 1D, we map it to the partition function of an Ising-like statistical mechanics model.As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly.We also use this mapping to give exact results for entanglement growth in 1 þ 1D circuits.
Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth.We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries.We study both 1 þ 1D and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems.In 1 þ 1D, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed v B .We find that in 1 þ 1D, the "front" of the OTOC broadens diffusively, with a width scaling in time as t 1=2 .We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization.Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem.This implies that the width of the front of the averaged OTOC scales as t 1=3 in 2 þ 1D and as t 0.240 in 3 þ 1D (exponents of the Kardar-Parisi-Zhang universality class).We support our analytic argument with simulations in 2 þ 1D.We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical.However, when full spatial rotational symmetry is present in 2 þ 1D, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution.For an alternative perspective on the OTOC in 1 þ 1D, we map it to the partition function of an Ising-like statistical mechanics model.As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly.We also use this mapping to give exact results for entanglement growth in 1 þ 1D circuits.
Characterizing how entanglement grows with time in a many-body system, for example, after a quantum quench, is a key problem in nonequilibrium quantum physics. We study this problem for the case of random unitary dynamics, representing either Hamiltonian evolution with time-dependent noise or evolution by a random quantum circuit. Our results reveal a universal structure behind noisy entanglement growth, and also provide simple new heuristics for the “entanglement tsunami” in Hamiltonian systems without noise. In 1D, we show that noise causes the entanglement entropy across a cut to grow according to the celebrated Kardar-Parisi-Zhang (KPZ) equation. The mean entanglement grows linearly in time, while fluctuations grow like (time)[superscript 1/3] and are spatially correlated over a distance ∝(time)[superscript 2/3]. We derive KPZ universal behavior in three complementary ways, by mapping random entanglement growth to (i) a stochastic model of a growing surface, (ii) a “minimal cut” picture, reminiscent of the Ryu-Takayanagi formula in holography, and (iii) a hydrodynamic problem involving the dynamical spreading of operators. We demonstrate KPZ universality in 1D numerically using simulations of random unitary circuits. Importantly, the leading-order time dependence of the entropy is deterministic even in the presence of noise, allowing us to propose a simple coarse grained minimal cut picture for the entanglement growth of generic Hamiltonians, even without noise, in arbitrary dimensionality. We clarify the meaning of the “velocity” of entanglement growth in the 1D entanglement tsunami. We show that in higher dimensions, noisy entanglement evolution maps to the well-studied problem of pinning of a membrane or domain wall by disorder. Subject Areas: Condensed Matter Physics, Quantum Information, Statistical Physics
The physics of information scrambling in quantum many-body systems is intimately related to thermalisation and emergence of chaos. However, its standard characterisation through bipartite entanglement or operator growth remains inherently coarse-grained, obscuring the spatiotemporal anatomy of how quantum information flows between different regions of the system and across various length scales. In this work, we develop a theory for the dynamics of local information in local random unitary circuits which resolves the fine-grained microscopic structures of local information flow in such systems. Using the framework of information lattice which systematises the information content at each length scale within each subsystem, we identify a distribution of length scales at which information resides and obtain the dynamics of the distribution. For Haar-random circuits with infinite local Hilbert-space dimensions, we map this dynamics onto an exact classical stochastic process which reveals that this distribution is described by a Tracy-Widom form which moves ballistically in time ($ \propto t$) toward larger scales, accompanied by a ∼t1/3 broadening. We also find the same qualitative behaviour for random Clifford circuits acting on qubits. The similar scaling behaviour in two rather different settings hints strongly towards the universality of our results. In the case of Clifford circuits, we develop a phenomenological Markov process for the dynamics of the stabiliser generators in a specific gauge which confirms the Tracy-Widom distribution and the t1/3 scaling of the fluctuations. Ultimately, our results establish the universal properties of the dynamics of local quantum information in a length scale-resolved fashion and provide a possible route towards bridging exact microscopic theories for information dynamics with emergent hydrodynamic descriptions of information flow.
Thermalization of chaotic quantum many-body systems under unitary time evolution is related to the growth in complexity of initially simple Heisenberg operators.Operator growth is a manifestation of information scrambling and can be diagnosed by out-of-time-order correlators (OTOCs).However, the behavior of OTOCs of local operators in generic chaotic local Hamiltonians remains poorly understood, with some semiclassical and large-N models exhibiting exponential growth of OTOCs and a sharp chaos wave front and other random circuit models showing a diffusively broadened wave front.In this paper, we propose a unified physical picture for scrambling in chaotic local Hamiltonians.We construct a random time-dependent Hamiltonian model featuring a large-N limit where the OTOC obeys a Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) type equation and exhibits exponential growth and a sharp wave front.We show that quantum fluctuations manifest as noise (distinct from the randomness of the couplings in the underlying Hamiltonian) in the FKPP equation and that the noise-averaged OTOC exhibits a crossover to a diffusively broadened wave front.At small N, we demonstrate that operator growth dynamics, averaged over the random couplings, can be efficiently simulated for all time using matrix product state techniques.To show that time-dependent randomness is not essential to our conclusions, we push our previous matrix product operator methods to very large size and show that data for a timeindependent Hamiltonian model are also consistent with a diffusively broadened wave front.
Dual-unitary circuits are paradigmatic examples of exactly solvable yet chaotic quantum many-body systems, but solvability naturally goes along with a degree of nongeneric behavior. By investigating the effect of weakly broken dual unitarity on the spreading of local operators, we study whether, and how, small deviations from dual unitarity recover fully generic many-body dynamics. We present a discrete path-integral formula for the out-of-time-order correlator and recover a butterfly velocity smaller than the light-cone velocity, v_{B}<v_{LC}, and a diffusively broadening operator front, two generic features of ergodic quantum spin chains absent in dual-unitary circuit dynamics. The butterfly velocity and diffusion constant are determined by a small set of microscopic quantities, and the operator entanglement of the gates has a crucial role.
The scrambling of quantum information in closed many-body systems, as measured by out-of-timeordered correlation functions (OTOCs), has received considerable attention lately.Recently, a hydrodynamical description of OTOCs has emerged from considering random local circuits.Numerical work suggests that aspects of this description are universal to ergodic many-body systems, even without randomness; a conjectured explanation for this is that while the random circuits have noise built into them, deterministic quantum systems, much like classically chaotic ones, "generate their own noise" and look effectively random on sufficient length scales and timescales.In this paper, we extend this approach to systems with locally conserved quantities (e.g., energy).We do this by considering local random unitary circuits with a conserved U(1) charge and argue, with numerical and analytical evidence, that the presence of a conservation law slows relaxation in both time-ordered and out-of-time-ordered correlation functions; both can have a diffusively relaxing component or "hydrodynamic tail" at late times.We verify the presence of such tails also in a deterministic, periodically driven system.We show that for OTOCs, the combination of diffusive and ballistic components leads to a wave front with a specific asymmetric shape, decaying as a power law behind the front.These results also explain existing numerical investigations in non-noisy ergodic systems with energy conservation.Moreover, we consider OTOCs in Gibbs states, parametrized by a chemical potential μ, and apply perturbative arguments to show that for μ ≫ 1 the ballistic front of information spreading can only develop at times exponentially large in μ-with the information traveling diffusively at earlier times.We also develop a new formalism for describing OTOCs and operator spreading, which allows us to interpret the saturation of OTOCs as a form of thermalization on the Hilbert space of operators.
Random unitary circuits have become a model system to investigate information scrambling in quantum systems. In the literature, mostly random circuits with Haar-distributed gate operations have been considered. In this work, we investigate operator spreading in random unitary circuits in which the elementary gate operations are drawn from general unitary-invariant ensembles, which include the well-studied Haar-distributed random unitary circuits as a special case. Similar to the Haar-distributed case, the long-time behavior of operator spreading with the more general unitary-invariant gate distribution is governed by drift-diffusion equations characterized by the butterfly velocity v B and a diffusion constant D . Differences with the Haar-random case are (i) that it takes a finite time τ b until ensemble-averaged Pauli-string weights take a “binary” form, in which they depend only on whether Pauli operators inside the support of the Pauli strong are equal to the identity matrix, and (ii) that the operator spreading is characterized by a finite “domain-wall width” n DW separating regions with a random-matrix-like Pauli-string distribution. To illustrate these findings, we perform explicit calculations for random unitary circuits distributed according to the Poisson kernel, which interpolates between the trivial and Haar-distributed circuits.