Operator spreading under unitary time evolution has attracted a lot of attention recently as a way to probe many-body quantum chaos. While quantities such as out-of-time-ordered correlators (OTOCs) do distinguish interacting from noninteracting systems, it has remained unclear to what extent they can truly diagnose chaotic vs integrable dynamics in many-body quantum systems. Here, we analyze operator spreading in generic one-dimensional many-body quantum systems using a combination of matrix product operator (MPO) and analytical techniques, focusing on the operator right weight. First, we show that while small bond dimension MPOs allow one to capture the exponentially decaying tail of the operator front, in agreement with earlier results, they lead to significant quantitative and qualitative errors for the actual front---defined by the maximum of the right weight. We find that while the operator front broadens diffusively in both integrable and chaotic interacting spin chains, the precise shape and scaling of the height of the front in integrable systems is anomalous for all accessible times. We interpret these results using a quasiparticle picture. This provides a sharp, although rather subtle, signature of many-body quantum chaos in the operator front.
Operator spreading under unitary time evolution has attracted a lot of attention recently as a way to probe many-body quantum chaos. While quantities such as out-of-time-ordered correlators (OTOCs) do distinguish interacting from noninteracting systems, it has remained unclear to what extent they can truly diagnose chaotic vs integrable dynamics in many-body quantum systems. Here, we analyze operator spreading in generic one-dimensional many-body quantum systems using a combination of matrix product operator (MPO) and analytical techniques, focusing on the operator right weight. First, we show that while small bond dimension MPOs allow one to capture the exponentially decaying tail of the operator front, in agreement with earlier results, they lead to significant quantitative and qualitative errors for the actual front---defined by the maximum of the right weight. We find that while the operator front broadens diffusively in both integrable and chaotic interacting spin chains, the precise shape and scaling of the height of the front in integrable systems is anomalous for all accessible times. We interpret these results using a quasiparticle picture. This provides a sharp, although rather subtle, signature of many-body quantum chaos in the operator front.
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.
Abstract In this article we study a set of integrable quantum cellular automata, the quantum hardcore gases (QHCG), with an arbitrary local Hilbert space dimension, and discuss the matrix product ansatz based approach for solving the dynamics of local operators analytically. Subsequently, we focus on the dynamics of operator spreading, in particular on the out-of-time ordered correlation functions (OTOCs), operator weight spreading and operators space entanglement entropy (OSEE). All of the quantities were conjectured to provide signifying features of integrable systems and quantum chaos. We show that in QHCG OTOCs spread diffusively and that in the limit of the large local Hilbert space dimension they increase linearly with time, despite their integrability. On the other hand, it was recently conjectured that operator weight front, which is associated with the extent of operators, spreads diffusively in both, integrable and generic systems, but its decay seems to differ in these two cases (Lopez-Piqueres et al 2021 Phys. Rev. B 104 104307). We observe that the spreading of the operator weight front in QHCG is markedly different from chaotic, generic integrable and free systems, as the front freezes in the long time limit. Finally, we discuss the OSEE in QHCG and show that it grows at most logarithmically with time in accordance with the conjectured behaviour for interacting integrable systems (Alba et al 2019 Phys. Rev. Lett. 122 250603).
The decay of out-of-time-ordered correlators (OTOCs) to a small residual value signals the scrambling of quantum information. This behavior is generic: it occurs in both chaotic and interacting integrable many-body systems, and therefore cannot by itself distinguish between them. In this paper we analytically study the late-time saturation value of a closely related observable, operator purity, and show that its scaling with system size robustly distinguishes chaotic from integrable systems. Our approach is based on the Mazur bound, which provides concrete scaling estimates and makes explicit why the locality of integrable charges matters. We support our findings with numerical simulations of one-dimensional spin chains.
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.
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.
Measuring out-of-time-ordered correlators (OTOCs) is experimentally challenging, often requiring effective reversal of many-body dynamics. Here, we propose an alternative approach to probing operator spreading through local quantum fluctuations, using only forward evolution and local measurements. Motivated by the Heisenberg uncertainty relation, we introduce a differential local variance and establish a rigorous bound connecting it to operator noncommutativity. For a ferromagnetic Heisenberg chain in the single-magnon regime, we derive an exact relation between the local variance and a squared-commutator OTOC. In integrable and nonintegrable Ising chains, the differential variance and OTOCs reveal qualitatively consistent propagation fronts. Our results establish local quantum fluctuations as an experimentally accessible, state-dependent, and complementary diagnostic of operator spreading.