It is widely expected that the von Neumann entanglement of a local operator grows at most logarithmically in integrable many-body systems in infinite volume. We give a counterexample in a four-state brickwork circuit whose gate is a permutation matrix solving the constant Yang--Baxter equation. The von Neumann operator entropy grows as $(\log2)\sqrt{t/π}+O(\log t)$, whereas fixed-index Rényi entropies grow linearly below index one and logarithmically above it. At fixed relative Hilbert--Schmidt error below one, matrix product operator simulations require a bond dimension growing at least as $\exp(c\sqrt t)$ for some $c>0$.
Comments
4.5 pages and 17 pages of Supplemental Material, v2: references added
It is widely expected that the von Neumann entanglement of a local operator grows at most logarithmically in integrable many-body systems in infinite volume. We give a counterexample in a four-state brickwork circuit whose gate is a permutation matrix solving the constant Yang--Baxter equation. The von Neumann operator entropy grows as (log2)t/π+O(logt), whereas fixed-index Rényi entropies grow linearly below index one and logarithmically above it. At fixed relative Hilbert--Schmidt error below one, matrix product operator simulations require a bond dimension growing at least as exp(ct) for some c>0.
We study the operator entanglement of local operators in one-dimensional brickwork circuits whose two-site gate satisfies the braid relation; throughout this work, we call such a gate a Yang--Baxter gate. We establish upper bounds for several structured, overlapping classes of Yang--Baxter gates. We show that the operator Schmidt rank remains uniformly bounded in time for all qubit Yang--Baxter gates and, in arbitrary local dimension, for permutation gates obtained from non-degenerate Yang--Baxter maps. We also show that it grows at most polynomially for involutive dual-unitary Yang--Baxter gates and for arbitrary phase dressings of permutation gates obtained from non-degenerate Yang--Baxter maps. These results imply, respectively, constant and logarithmic upper bounds on the operator entanglement. Conversely, we construct a seven-state involutive Yang--Baxter gate without dual unitarity and a one-site operator whose exact operator Schmidt rank grows exponentially, although the corresponding operator entropies remain undetermined. Entanglement growth in the general Yang--Baxter case remains open. All proofs and selected examples were constructed by ChatGPT 5.6 Sol.
We study a so-called semi-ergodic brickwork dual-unitary circuits where, in the infinite volume limit, the two-point correlation functions of single-site operators exhibit ergodic behavior along one light ray and non-ergodic behavior along the other light ray. Here, however, we study intermediate and long-time dynamics of a system in a finite, large volume. Under such dynamics, the Heisenberg evolution of a single traceless single-site operator lies within a restricted subspace, and this time evolution can be mapped to a simpler problem of a single qutrit scattering with a bunch of qubits sequentially. Despite the model being non-integrable and free from any quenched disorder, the operator entanglement grows at most logarithmic in time, contrary to prior expectations. The auto-correlation function can be written in terms of a sum of products of SO(3) matrices, allowing for a random matrix prediction for the auto-correlation function at late times. The operator size distribution also becomes bimodal at certain times, displaying intermediate behavior between chaotic and free systems.
We aim to rigorously address how local noise affects classical simulability of quantum dynamics benchmarked by tensor-network methods. Using operator entanglement entropy (OEE), we prove the following: (1) For single-qubit depolarizing noise on arbitrary circuits, tensor networks with poly(n) bond dimension suffice for fixed absolute Hilbert-Schmidt error after $\order{1}$ depth, while relative error demands $\order{\log n}$ depth; and this bound is optimal. (2) For single-qubit depolarizing noise on 1D local circuits, the existence of whole-trajectory error-bounded matrix product operator (MPO) of poly(n) bond dimension at all depths. (3) For general single-qubit noise in 1D brickwall circuits, random two-design gates with contraction coefficient c<1/3 yield an $\order{1}$ OEE plateau with probability 1−Te−Ω(n), while arbitrary gates with c<1/48 give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth poly(n) average boundary-bond dimensions for projected entangled pair operators~(PEPO) across every cut -- under depolarizing noise at either absolute or relative accuracy, and under general noise with strong contraction at absolute accuracy. Our results establish a rigorous connection between certain noise models, circuit types, and their classical simulability.
We compute the dynamics of entanglement in the minimal setup producing ergodic and mixing quantum many-body dynamics, which we previously dubbed boundary chaos. This consists of a free, non-interacting brickwork quantum circuit, in which chaos and ergodicity is induced by an impurity interaction, i.e., an entangling two-qudit gate, placed at the system’s boundary. We compute both the conventional bipartite entanglement entropy with respect to a connected subsystem including the impurity interaction for initial product states as well as the so-called operator entanglement entropy of initial local operators. Thereby we provide exact results in a particular scaling limit of both time and system size going to infinity for either very small or very large subsystems. We show that different classes of impurity interactions lead to very distinct entanglement dynamics. For impurity gates preserving a local product state forming the bulk of the initial state, entanglement entropies of states show persistent spikes with period set by the system size and suppressed entanglement in between, contrary to the expected linear growth in ergodic systems. We observe similar dynamics of operator entanglement for generic impurities. In contrast, for T-dual impurities, which remain unitary under partial transposition, we find entanglement entropies of both states and operators to grow linearly in time with the maximum possible speed allowed by the geometry of the system. The intensive nature of interactions in all cases causes entanglement to grow on extensive time scales proportional to system size.
Local-operator entanglement (LOE) quantifies the nonlocal structure of Heisenberg operators and serves as a diagnostic of many-body chaos. We provide rigorous bounds showing when an operator can be well-approximated by a matrix-product operator (MPO), given asymptotic scaling of its LOE α-Rényi entropies. Specifically, we prove that a volume law scaling for α≥1 implies that the operator cannot be approximated efficiently as an MPO while faithfully reproducing all expectation values. On the other hand, if we restrict to correlations over a relevant sub-class of (ensembles of) states, then logarithmic scaling of the α<1 entropies implies MPO simulability. This result covers a range of relevant quantities, including infinite temperature autocorrelation functions, out-of-time-ordered correlators, and average-case expectation values over ensembles of computational basis states. Beyond this regime, we provide numerical evidence together with a random matrix model to argue that this simulability result also typically holds for arbitrary states. Our results put on firm footing the heuristic expectation that a low operator entanglement implies efficient tensor network representability, extending celebrated foundational results from the theory of matrix-product states and providing a formal link between quantum chaos and classical simulability.
Free fermion models with an operator entanglement growing as O(t)
We study the time evolution of the Operator Space Entanglement Entropy (OSEE) of one-dimensional number-conserving free fermions subject to a random space- and time-dependent evolution. It is known that, in free-fermion systems without randomness, the OSEE of \emph{local} operators grows at most \emph{logarithmically} in time. Here we show that, in the presence of randomness, the OSEE of local spin operators involving a Jordan-Wigner string grows as O(t), for all values of the Rényi index. We establish this in two models: a Brownian Hamiltonian and a brick-wall matchgate circuit. We also consider the OSEE of the evolution operator, and we show that it also grows as O(t), in contrast with the linear growth generically expected without randomness. We analytically compute the average entropies of the two operators by mapping the quantum models to the symmetric simple exclusion process, relying on an assumption of self-averaging which we support numerically. We also consider two other types of randomness: random evolution in time that keeps spatial translation invariance, and quenched disorder. In the former case, we find that the OSEE also grows as t. In the latter case, we find that the entropy grows as O(loglogt) for Hamiltonian evolution both for local operators with a Jordan-Wigner string and for the evolution operator, whereas it remains bounded at all times in the brick-wall matchgate circuit.