Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2603.05656 · Physical Review A

Classical Simulability from Operator Entanglement Scaling

Neil Dowling

quant-phcond-mat.stat-mechhep-thSubmitted Mar 5, 2026• Announced Mar 9, 2026• 0 citations• Revised Oct 19, 2026 (v4)
v1 · Mar 5, 2026v2 · Mar 23, 2026v3 · May 26, 2026v4 · Oct 19, 2026
arXivHTMLPDF

Abstract

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 $α\geq 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.

Comments
6+11 pages, 2+2 figures. v4: fixed minor errors

Related papers

Page Curve for Local-Operator Entanglement from Free Probability

Neil Dowling, Silvia Pappalardi

The local-operator entanglement (LOE) measures the classical simulability of a Heisenberg operator and is conjectured to witness many-body chaos in locally interacting systems. Using tools from free probability, we analytically compute its asymptotic value for Haar random dynamics for all Rényi indices. We find that to the first two orders in large dimension, for traceless operators it exactly reproduces the Page curve for random states. In contrast to higher-order out-of-time-ordered correlators (OTOCs), which depend on more complex operator correlations via free cumulants, the leading-order LOE is independent of the initial operator. \add{We detail exactly how these higher-order correlations are encoded in the coefficients of the finite-size corrections. Therefore, guided by our Haar result, we argue that the long-time value of the LOE entropies in chaotic systems is robust: they depend only on a simple sub-sector of the hierarchy of correlations encoded by the full Eigenstate Thermalization Hypothesis. Our results therefore provide a conceptual boundary between operator simulability and higher-order scrambling diagnostics.

quant-pharXiv:2605.02995•May 6, 2026•0 citations

Bridging Entanglement and Magic Resources within Operator Space

Neil Dowling, Kavan Kishore Modi, Gregory A. L. White

Local-operator entanglement (LOE) dictates the complexity of simulating Heisenberg evolution using tensor network methods, and it bears witness to many-body chaos for local dynamics. We show that LOE is also sensitive to how non-Clifford a unitary is-its magic resources. In particular, we prove that LOE is always upper-bound by three distinct magic monotones: T-count, unitary nullity, and operator-stabilizer Rényi entropy. Moreover, in the average case for large, random circuits, LOE and magic monotones approximately coincide. Our results imply that an operator evolution that is expensive to simulate using tensor network methods must also be inefficient using both stabilizer and Pauli truncation methods. In terms of a previous conjecture on the characteristic scaling of LOE, our results imply that nonintegrable spin chains cannot be simulated classically. Entanglement in operator space therefore measures a unified picture of nonclassical resources, in stark contrast to the Schrödinger picture.

arXiv:2501.18679•Physical Review Letters•Sep 18, 2025•12 citations

A derivation of the late-time volume law for local operator entanglement

Guilherme Ilário Correr, John Goold, Marco Cattaneo

Local Operator Entanglement (LOE) has emerged an indicator of quantum chaos in many-body systems. Numerical studies have shown that, in chaotic systems, LOE grows linearly in time and displays a volume-law behavior at late times, scaling proportionally with the number of local degrees of freedom. Despite extensive numerical evidence, complemented by analytical studies in integrable systems, a fully analytical understanding of the emergence of the volume law remains incomplete. In this paper, we contribute toward this goal by deriving a late-time expression for LOE in chaotic systems that exhibits volume-law scaling. Our derivation proceeds by expressing the late-time LOE in the Liouville eigenstate basis and relies on three main assumptions: a higher-order non-resonance condition for the Hamiltonian eigenenergies, the Eigenstate Thermalization Hypothesis (ETH) ansatz for the matrix elements of the initial local operator, and the replacement of Hamiltonian eigenstates with random states in the final expression for LOE. Under these assumptions, we obtain an explicit formula displaying volume-law scaling. Finally, we complement our analytical derivation with numerical simulations of the 1D mixed-field Ising model, testing the resulting formula and exploring the regime of validity of our assumptions.

quant-pharXiv:2603.25387•Journal of Physics A Mathematical and Theoretical•Mar 27, 2026•0 citations

Complexity of tensor network simulation for noisy quantum circuits

Yuguo Shao, Zishuo Zhao, Song Cheng +1 more

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 $\mathrm{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 $\mathrm{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 $\mathrm{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.

cond-mat.str-elCross-listedarXiv:2606.00474•Jun 2, 2026•0 citations

Long-time limits of local operator entanglement in interacting integrable models

J. Alexander Jacoby, Sarang Gopalakrishnan

We explore the long-time behavior of Local Operator Entanglement entropy (LOE) in finite-size interacting integrable systems. For certain operators in the Rule 54 automaton, we prove that the LOE saturates to a value that is at most logarithmic in system size. This bound extends previous work [Phys. Rev. Lett. 122, 250603 (2019); Commun. Math. Phys. 371, 651 (2019)] showing LOE grows logarithmically in the early time regime, t\ll L t ≪ L , to the late time regime, $tL $. However, the late-time logarithmic bound relies on a feature of Rule 54 that does not generalize to other interacting integrable systems: namely, that there are only two types of quasiparticles, and therefore only two possible values of the phase shift between quasiparticles. We present a heuristic argument, supported by numerical evidence, that for generic interacting integrable systems (such as the Heisenberg spin chain) the saturated value of the LOE is volume-law in system size.

arXiv:2503.09578•SciPost Physics•Jun 11, 2026•0 citations

Mapping twist fields to local operators via tensor networks

Andrea Bulgarelli, Marco Panero, Paolo Stornati +1 more

Twist fields are a powerful formal tool to compute Rényi entropies in quantum many-body systems, but their conventional formulation in tensor network states involves operations acting on virtual degrees of freedom, which are not directly accessible in experiments. In this work, we construct explicit local operators acting on the physical Hilbert space whose expectation values reproduce the action of twist fields in matrix product states. Our construction is exact in the injectivity limit and when the tensor is chosen at the center of orthogonality, and provides a direct operational method to evaluate Rényi entropies without accessing auxiliary tensor indices. We test our formulation numerically in the transverse-field Ising model, demonstrating rapid convergence to the exact entanglement entropy as the injectivity scale is reached. Furthermore, we show that twist operators determined from relatively small reference systems can be reliably transferred to larger systems, once the reference size exceeds a characteristic scale set by the correlation length. Since the resulting operators admit a decomposition in terms of a finite number of local observables, our results provide a scalable and experimentally accessible framework to probe entanglement in quantum simulators.

quant-pharXiv:2605.26208•May 27, 2026•0 citations
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α\geq 1α≥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α< 1α<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.
The local-operator entanglement (LOE) measures the classical simulability of a Heisenberg operator and is conjectured to witness many-body chaos in locally interacting systems. Using tools from free probability, we analytically compute its asymptotic value for Haar random dynamics for all Rényi indices. We find that to the first two orders in large dimension, for traceless operators it exactly reproduces the Page curve for random states. In contrast to higher-order out-of-time-ordered correlators (OTOCs), which depend on more complex operator correlations via free cumulants, the leading-order LOE is independent of the initial operator. \add{We detail exactly how these higher-order correlations are encoded in the coefficients of the finite-size corrections. Therefore, guided by our Haar result, we argue that the long-time value of the LOE entropies in chaotic systems is robust: they depend only on a simple sub-sector of the hierarchy of correlations encoded by the full Eigenstate Thermalization Hypothesis. Our results therefore provide a conceptual boundary between operator simulability and higher-order scrambling diagnostics.
Local-operator entanglement (LOE) dictates the complexity of simulating Heisenberg evolution using tensor network methods, and it bears witness to many-body chaos for local dynamics. We show that LOE is also sensitive to how non-Clifford a unitary is-its magic resources. In particular, we prove that LOE is always upper-bound by three distinct magic monotones: T-count, unitary nullity, and operator-stabilizer Rényi entropy. Moreover, in the average case for large, random circuits, LOE and magic monotones approximately coincide. Our results imply that an operator evolution that is expensive to simulate using tensor network methods must also be inefficient using both stabilizer and Pauli truncation methods. In terms of a previous conjecture on the characteristic scaling of LOE, our results imply that nonintegrable spin chains cannot be simulated classically. Entanglement in operator space therefore measures a unified picture of nonclassical resources, in stark contrast to the Schrödinger picture.
Local Operator Entanglement (LOE) has emerged an indicator of quantum chaos in many-body systems. Numerical studies have shown that, in chaotic systems, LOE grows linearly in time and displays a volume-law behavior at late times, scaling proportionally with the number of local degrees of freedom. Despite extensive numerical evidence, complemented by analytical studies in integrable systems, a fully analytical understanding of the emergence of the volume law remains incomplete. In this paper, we contribute toward this goal by deriving a late-time expression for LOE in chaotic systems that exhibits volume-law scaling. Our derivation proceeds by expressing the late-time LOE in the Liouville eigenstate basis and relies on three main assumptions: a higher-order non-resonance condition for the Hamiltonian eigenenergies, the Eigenstate Thermalization Hypothesis (ETH) ansatz for the matrix elements of the initial local operator, and the replacement of Hamiltonian eigenstates with random states in the final expression for LOE. Under these assumptions, we obtain an explicit formula displaying volume-law scaling. Finally, we complement our analytical derivation with numerical simulations of the 1D mixed-field Ising model, testing the resulting formula and exploring the regime of validity of our assumptions.
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)\mathrm{poly}(n)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)\mathrm{poly}(n)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/3c<1/3c<1/3 yield an $\order{1}$ OEE plateau with probability 1−Te−Ω(n)1-Te^{-Ω(n)}1−Te−Ω(n), while arbitrary gates with c<1/48c<1/48c<1/48 give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth poly(n)\mathrm{poly}(n)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 explore the long-time behavior of Local Operator Entanglement entropy (LOE) in finite-size interacting integrable systems. For certain operators in the Rule 54 automaton, we prove that the LOE saturates to a value that is at most logarithmic in system size. This bound extends previous work [Phys. Rev. Lett. 122, 250603 (2019); Commun. Math. Phys. 371, 651 (2019)] showing LOE grows logarithmically in the early time regime, t\ll L t ≪ L , to the late time regime, $tL $. However, the late-time logarithmic bound relies on a feature of Rule 54 that does not generalize to other interacting integrable systems: namely, that there are only two types of quasiparticles, and therefore only two possible values of the phase shift between quasiparticles. We present a heuristic argument, supported by numerical evidence, that for generic interacting integrable systems (such as the Heisenberg spin chain) the saturated value of the LOE is volume-law in system size.
Twist fields are a powerful formal tool to compute Rényi entropies in quantum many-body systems, but their conventional formulation in tensor network states involves operations acting on virtual degrees of freedom, which are not directly accessible in experiments. In this work, we construct explicit local operators acting on the physical Hilbert space whose expectation values reproduce the action of twist fields in matrix product states. Our construction is exact in the injectivity limit and when the tensor is chosen at the center of orthogonality, and provides a direct operational method to evaluate Rényi entropies without accessing auxiliary tensor indices. We test our formulation numerically in the transverse-field Ising model, demonstrating rapid convergence to the exact entanglement entropy as the injectivity scale is reached. Furthermore, we show that twist operators determined from relatively small reference systems can be reliably transferred to larger systems, once the reference size exceeds a characteristic scale set by the correlation length. Since the resulting operators admit a decomposition in terms of a finite number of local observables, our results provide a scalable and experimentally accessible framework to probe entanglement in quantum simulators.