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).
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).
Understanding the spreading of the operator space entanglement entropy (OSEE) is key in order to explore out-of-equilibrium quantum many-body systems. Here we argue that for integrable models the dynamics of the OSEE is related to the diffusion of the operator front. We derive the logarithmic bound 1/2ln(t) for the OSEE of some simple, i.e., low-rank, diagonal local operators. We numerically check that the bound is saturated in the rule 54 chain, which is representative of interacting integrable systems. Remarkably, the same bound is saturated in the spin-1/2 Heisenberg XXZ chain. Away from the isotropic point and from the free-fermion point, the OSEE grows as 1/2ln(t), irrespective of the chain anisotropy, suggesting universality. Finally, we discuss the effect of integrability breaking. We show that strong finite-time effects are present, which prevent us from probing the asymptotic behavior of the OSEE.
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.
In a many-body quantum system, local operators in the Heisenberg picture O(t)=e^{iHt}Oe^{-iHt} spread as time increases. Recent studies have attempted to find features of that spreading which could distinguish between chaotic and integrable dynamics. The operator entanglement-the entanglement entropy in operator space-is a natural candidate to provide such a distinction. Indeed, while it is believed that the operator entanglement grows linearly with time t in chaotic systems, we present evidence that it grows only logarithmically in generic interacting integrable systems. Although this logarithmic growth has been previously established for noninteracting fermions, there has been no progress on interacting integrable systems to date. In this Letter we provide an analytical upper bound on operator entanglement for all local operators in the "Rule 54" qubit chain, a cellular automaton model introduced in the 1990s [Bobenko et al., CMP 158, 127 (1993)CMPHAY0010-361610.1007/BF02097234], and recently advertised as the simplest representative of interacting integrable systems. Physically, the logarithmic bound originates from the fact that the dynamics of the models is mapped onto the one of stable quasiparticles that scatter elastically. The possibility of generalizing this scenario to other interacting integrable systems is briefly discussed.
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.
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.
Scrambling is a key concept in the analysis of nonequilibrium properties of quantum many-body systems. Most studies focus on its characterization via out-of-time-ordered correlation (OTOC) functions, particularly through the early-time decay of the OTOC. However, scrambling is a complex process which involves operator spreading and operator entanglement, and a full characterization requires one to access more refined information on the operator dynamics at several timescales. In this work we analyze operator scrambling by expanding the target operator in a complete basis and studying the structure of the expansion coefficients treated as a coarse-grained probability distribution in the space of operators. We study different features of this distribution, such as its mean, variance, and participation ratio, for the Ising model with longitudinal and transverse fields, kicked collective spin models, and random circuit models. We show that the long-time properties of the operator distribution display common features across these cases and discuss how these properties can be used as a proxy for the onset of quantum chaos. Finally, we discuss the connection with OTOCs and analyze the cost of probing the operator distribution experimentally using these correlation functions.