Operator Spreading in Random Unitary Circuits

• Published in: Physical review. X 2018-04, Vol.8 (2), p.021014, Article 021014
• Main Authors: Nahum, Adam, Vijay, Sagar, Haah, Jeongwan
• Format: Article
• Language: English
• Subjects: Circuits; CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; Hydrodynamic equations; Information retrieval; Ising model; Mapping; Partitions (mathematics); Quantum entanglement; Quantum phenomena; Quantum theory; Random walk; Statistical mechanics; System effectiveness
• ISSN: 2160-3308; 2160-3308
• DOI: 10.1103/PhysRevX.8.021014

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 both1+1Dand higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems. In1+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 speedvB. We find that in1+1D, the “front” of the OTOC broadens diffusively, with a width scaling in time ast1/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 ast1/3in2+1Dand ast0.240in3+1D(exponents of the Kardar-Parisi-Zhang universality class). We support our analytic argument with simulations in2+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 in2+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 in1+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 in1+1Dcircuits.