pith. machine review for the scientific record. sign in

arxiv: 2604.16720 · v1 · submitted 2026-04-17 · ❄️ cond-mat.stat-mech · nlin.CD· quant-ph

Recognition: unknown

Quantum many-body operator cascade as a route to chaos

Authors on Pith no claims yet

Pith reviewed 2026-05-10 06:44 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech nlin.CDquant-ph
keywords quantum many-body chaosoperator cascadefractal dimensionRuelle-Pollicott eigenvectorsKolmogorov cascadetruncated propagatorunitarity constraintnon-local operators
0
0 comments X

The pith

In quantum many-body chaos, local operators develop fractal non-locality as they decay, linking time and space scales through unitarity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that the slowest-decaying operators in quantum chaotic systems, identified as leading eigenvectors of the truncated operator propagator, exhibit a fractal dimension that quantifies their non-locality. This structure is also seen in the divergence of their condition numbers. Unitarity imposes an approximate equality between the rate at which local correlations decay in time and this spatial fractal dimension. This leads to a picture where operators cascade from local to non-local fractal forms, causing effective relaxation in the local operator space, similar to a Kolmogorov cascade but in quantum operator space. The findings are shown in models like the kicked Ising model and dual-unitary circuits where they are exact.

Core claim

By examining the spectral properties of the truncated operator propagator, the slowest-decaying operators are shown to be the leading Ruelle-Pollicott eigenvectors with nontrivial fractal dimension quantifying non-locality. Unitarity constrains the temporal decay rate of local correlations to equal this spatial operator fractal dimension approximately. This reveals a many-body quantum chaos scenario where local operators evolve towards increasingly non-local fractal structures, naturally inducing effective non-unitary relaxation on the local operators subspace, akin to a Kolmogorov cascade in operator space.

What carries the argument

The truncated operator propagator and its leading Ruelle-Pollicott eigenvectors, which carry the fractal dimension of operator non-locality.

If this is right

  • Local operators evolve over time into non-local ones with quantifiable fractal structure.
  • Effective non-unitary relaxation emerges on the subspace of local operators.
  • The fractal dimension appears through diverging condition numbers of the operators.
  • Predictions are exact in dual-unitary circuits and hold in kicked Ising and random brickwall circuits.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The fractal measure of operator non-locality could be used to test analogies between quantum many-body chaos and classical turbulence cascades.
  • The method offers a route to estimate relaxation timescales from spatial operator properties without simulating full dynamics.
  • Increasing system size in simulations would test whether the fractal dimension stabilizes as assumed.

Load-bearing premise

The truncated operator propagator's leading eigenvectors faithfully capture the long-time non-local structure of the full unitary evolution, and the fractal dimension remains well-defined and stable in the thermodynamic limit.

What would settle it

Numerical checks in larger systems showing that condition numbers of slowest-decaying operators stay bounded or that decay rates fail to match the measured fractal dimension would falsify the proposed link.

Figures

Figures reproduced from arXiv: 2604.16720 by Marko \v{Z}nidari\v{c}, Urban Duh.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: that the ring radii are indeed roughly given by Eq. (32). Note that this noisy ring-like spectrum is very [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: (c). Two important observations are in place: (i) to achieve an approximate collaps of data for different r and s one has to use a large scaling factor ≈ 1.38 on the verti￾cal axis, and (ii) growth does not seem to be exponential at large s (the whole “collapse” is much worse and it is hard to judge whether the growth is algebraic or expo￾nential in s). The fact that the scaling in ws has a larger factor 1… view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
read the original abstract

Dynamical properties of classical chaotic systems, for instance relaxation, can be understood as emerging from the time evolution of initially smooth long-wavelength densities to ever finer short-wavelength densities with fractal structure. Whether there is any analogous fractality by which one could characterize quantum many-body chaos is not known. By studying the spectral properties of the truncated operator propagator, we provide such structures. Namely, we show that the slowest-decaying operators, i.e., the leading Ruelle-Pollicott eigenvectors, have a nontrivial fractal dimension quantifying their non-locality, visible also in the divergence of their condition numbers. Furthermore, we find that unitarity imposes a constraint, i.e., an (approximate) equality, between the temporal decay rate of local correlations and this spatial operator fractal dimension. With this insight, a scenario for many-body quantum chaos becomes clear: over time, local operators evolve towards increasingly non-local ones with a quantifiable fractal structure, thereby naturally leading to effective non-unitary relaxation on the subspace of local operators - a kind of many-body Kolmogorov cascade in the space of operators. Our predictions are demonstrated in various quantum circuits: the kicked Ising model, brickwall circuits with a random 2-qubit gate, and dual-unitary circuits, where our results are exact.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript studies the spectral properties of the truncated operator propagator in quantum circuits to characterize many-body quantum chaos. It claims that the leading (slowest-decaying) Ruelle-Pollicott eigenvectors possess a nontrivial fractal dimension that quantifies operator non-locality, visible also via diverging condition numbers, and that unitarity enforces an approximate equality between the temporal decay rate of local correlations and this spatial fractal dimension. The results are presented as exact in dual-unitary circuits and supported by numerics in the kicked Ising model and random brickwall circuits, framing quantum chaos as an operator-space Kolmogorov cascade that produces effective non-unitary relaxation on local operators.

Significance. If the central claims hold, the work supplies a concrete, quantifiable notion of fractality in the operator space of quantum many-body systems and a direct link between spatial operator structure and temporal correlation decay. The exact results in dual-unitary circuits constitute a clear strength, as does the reproducible numerical evidence for the approximate equality in other models. This perspective could influence studies of operator growth, scrambling, and thermalization by providing a geometric diagnostic of chaos beyond standard out-of-time-order correlators.

major comments (3)
  1. [numerical results and dual-unitary analysis] The central claim that the fractal dimension remains nontrivial and stable in the thermodynamic limit rests on finite truncations of the operator propagator. No explicit scaling of the extracted dimension with truncation depth or system size L is provided to demonstrate convergence as both limits are taken simultaneously; this is load-bearing for the assertion that the leading eigenvectors faithfully capture long-time non-local structure (abstract and numerical sections).
  2. [abstract and § on unitarity constraint] The reported (approximate) equality between the temporal decay rate of local correlations and the spatial fractal dimension is presented as imposed by unitarity. It is unclear whether this relation is derived parameter-free from the unitary evolution or whether the dimension is extracted from the same truncated correlation data, which would introduce circularity; an independent derivation or cross-validation would be required to substantiate the constraint (abstract).
  3. [dual-unitary circuits] In the dual-unitary circuits the exactness is stated to hold only within the truncated operator subspace. The manuscript should clarify whether the fractal dimension and the equality survive removal of the truncation or remain artifacts of the finite subspace (dual-unitary section).
minor comments (2)
  1. [methods] Notation for the truncated propagator and the precise definition of the fractal dimension (e.g., via box-counting or correlation integral) should be stated explicitly in the main text rather than only in appendices.
  2. [figures] Figure captions should include the truncation depth and system sizes used for each panel to allow direct assessment of finite-size effects.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their positive evaluation of the significance of our work and for the constructive major comments. We have revised the manuscript to strengthen the evidence for convergence, clarify the derivation of the unitarity constraint, and address the status of the exact results in dual-unitary circuits.

read point-by-point responses
  1. Referee: The central claim that the fractal dimension remains nontrivial and stable in the thermodynamic limit rests on finite truncations of the operator propagator. No explicit scaling of the extracted dimension with truncation depth or system size L is provided to demonstrate convergence as both limits are taken simultaneously; this is load-bearing for the assertion that the leading eigenvectors faithfully capture long-time non-local structure (abstract and numerical sections).

    Authors: We agree that explicit scaling with both truncation depth and system size is necessary to support the thermodynamic-limit claims. In the revised manuscript we have added new panels to the numerical sections (for the kicked Ising model and random brickwall circuits) that display the extracted fractal dimension versus truncation depth for several values of L. These plots show rapid stabilization of the dimension with increasing truncation depth and consistency of the limiting value as L grows, thereby demonstrating the required convergence. revision: yes

  2. Referee: The reported (approximate) equality between the temporal decay rate of local correlations and the spatial fractal dimension is presented as imposed by unitarity. It is unclear whether this relation is derived parameter-free from the unitary evolution or whether the dimension is extracted from the same truncated correlation data, which would introduce circularity; an independent derivation or cross-validation would be required to substantiate the constraint (abstract).

    Authors: The approximate equality follows directly from the unitarity of the underlying evolution, which constrains the spectrum and eigenvectors of the truncated propagator. We have revised the abstract and the dedicated unitarity-constraint section to present an explicit, parameter-free derivation of the relation that does not rely on numerically extracting the dimension from the same correlation data. The derivation is cross-validated by the exact analytic results obtained in dual-unitary circuits, where the equality holds independently of any numerical fitting procedure. revision: yes

  3. Referee: In the dual-unitary circuits the exactness is stated to hold only within the truncated operator subspace. The manuscript should clarify whether the fractal dimension and the equality survive removal of the truncation or remain artifacts of the finite subspace (dual-unitary section).

    Authors: In dual-unitary circuits the dual-unitary property renders the relevant operator subspace invariant, so that the leading spectral features—including the fractal dimension of the Ruelle-Pollicott eigenvectors and the equality with the correlation decay rate—are exact within that subspace and remain unchanged upon increasing the truncation depth. We have added an analytic argument and supporting plots in the dual-unitary section showing that both quantities converge immediately to their exact values once the minimal truncation depth required by the dual-unitary structure is reached; hence they are not artifacts of the finite subspace. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivations are numerical and self-contained

full rationale

The paper computes the truncated operator propagator explicitly in concrete models (kicked Ising, brickwall, dual-unitary circuits) and extracts the leading Ruelle-Pollicott eigenvectors and their fractal dimensions directly from those finite matrices. The reported approximate equality between local correlation decay rate and operator fractal dimension is presented as a numerical observation under unitarity within the same truncated spaces, with exactness claimed only inside dual-unitary truncations. No equation reduces a claimed prediction to a fitted input by construction, no load-bearing premise rests on self-citation, and no ansatz is smuggled via prior work. The central chain (truncation → eigenvector spectrum → fractal dimension → observed equality) remains independent of the target claims and is supported by explicit model-specific calculations rather than tautological redefinition.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the spectral properties of a truncated propagator and on the assumption that unitarity directly equates a temporal rate to a spatial fractal dimension; no explicit free parameters or invented entities are named in the abstract.

axioms (2)
  • domain assumption The truncated operator propagator preserves the leading long-time spectral features of the full unitary evolution.
    Invoked to extract Ruelle-Pollicott eigenvectors and their fractal properties.
  • domain assumption A well-defined fractal dimension can be assigned to the slowest-decaying operators in the many-body limit.
    Required for the non-locality quantification and the unitarity constraint.
invented entities (1)
  • Operator fractal dimension no independent evidence
    purpose: Quantify non-locality of leading Ruelle-Pollicott eigenvectors
    New diagnostic introduced to characterize the cascade; no independent falsifiable prediction supplied in abstract.

pith-pipeline@v0.9.0 · 5530 in / 1515 out tokens · 30293 ms · 2026-05-10T06:44:44.283370+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

102 extracted references · 11 canonical work pages · 4 internal anchors

  1. [1]

    attractor

    Diverging leading eigenvector Because we have|λ 1|<1, i.e.d T >0, we know that in r→ ∞limit,|R 1⟩(and|L 1⟩) cannot be fromℓ 2. Remem- ber that|R 1⟩is an operator to which any local operator Awith a finite overlap⟨L 1|α⟩will approach forward in time,t→ ∞, while|L 1⟩is a similar asymptotic operator for evolution backward in time,t→ −∞. We are going to show ...

  2. [2]

    information

    Fractal dimension Let us recall the characterization of fractals by one of the simplest measures, the box-counting dimension. Taking a box of linear sizeϵ, one measures the number of boxesN(ϵ) covered by a set (a fractal). The frac- tal dimensiondis then defined by the asymptotic scal- ingN(ϵ)∼(1/ϵ) d, formallyd= lim ϵ→0 lnN(ϵ)/ln (1/ϵ). Structures in cla...

  3. [3]

    Multifractality We have checked for the presence of multifractality in the kicked Ising model by calculating partial norms of |R1⟩for few exponentsm, ws(m) = X α,sup(α)=s |bα|m.(A1) Due to fractality one asymptotically expectsw s(m)∼ µs(m), and correspondingly anm-dependent dimension dx(m) := lnµ(m). In the main text we have focused on dx =d x(2), here we...

  4. [4]

    collapse

    Non-scaling form ofw s Among various parameters of the kicked Ising model that we have checked we found one case forτ= 0.45 where the partial normsw s do not seem to have a nice scaling form like for other parameters. Smallerτ= 0.45 is in fact a situation where one has strong prethermal- ization with a very long timescale on which the almost- conserved ef...

  5. [5]

    III A we heuristically argued for the exponential convergence of RP resonances

    Perturbation theory in increasing support In Sec. III A we heuristically argued for the exponential convergence of RP resonances. In this appendix, we work towards making the argument more precise by employ- ing perturbation theory. Our approximation, however, is mathematically not well-controlled because, as men- tioned, the perturbation is small only on...

  6. [6]

    |v1⟩ |v2⟩ # . In this notation, the eigenvalue equation is

    0B i+1 Ai+1   (C2) with knowing the solution to the eigenproblem ofU r, i.e., its leading eigenvalueλ (r) 1 and corresponding left and right eigenvectors⟨L (r)|and|R (r)⟩. We denoted i=r/(2v) and assumed thatris a multiple of 2vfor simplicity. First, we relabel the whole top right and bot- tom left block and introduce a constantεthat will help in ...

  7. [7]

    V. I. Arnold and A. Avez,Ergodic problems of classical mechanics(W. A. Benjamin, New York, 1968)

  8. [8]

    Ott,Chaos in dynamical systems(Cambridge Univer- sity Press, 1993)

    E. Ott,Chaos in dynamical systems(Cambridge Univer- sity Press, 1993)

  9. [9]

    A. J. Lichtenberg and M. A. Lieberman,Regular and chaotic dynamics(Springer, New York, 1992)

  10. [10]

    A stable (unstable) manifold is a set of points that flows to (away) from the fixed point

  11. [11]

    Cvitanovi´ c, R

    P. Cvitanovi´ c, R. Artuso, R. Mainieri, G. Tan- ner, and G. Vattay,Chaos: Classical and Quantum chaosbook.org(Niels Bohr Institute, Copenhagen 2020)

  12. [12]

    M. A. Nielsen and I. L. Chuang,Quantum Computation and Quantum Information(Cambridge University Press, Cambridge, 2010)

  13. [13]

    Haake,Quantum signatures of chaos(Springer, 2010)

    F. Haake,Quantum signatures of chaos(Springer, 2010)

  14. [14]

    M. C. Gutzwiller,Chaos in classical and quantum me- chanics, (Springer, New York, 1990)

  15. [15]

    D’Alessio, Y

    L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics, Adv. Phys. 65, 239 (2016)

  16. [16]

    Gaspard,Chaos, scattering and statistical mechanics (Cambridge University Press, Cambridge, 1998)

    P. Gaspard,Chaos, scattering and statistical mechanics (Cambridge University Press, Cambridge, 1998)

  17. [17]

    J. B. Conway,A Course in Functional Analysis (Springer, New York, 1990)

  18. [18]

    I. M. Gel’fand, and N. Ya. Vilenkin,Generalized Func- tions: Applications of Harmonic Analysis(Academic Press, New York, 1964)

  19. [19]

    Pollicott,On the rate of mixing in axiom A flows, Invent

    M. Pollicott,On the rate of mixing in axiom A flows, Invent. Math.81, 413 (1985)

  20. [20]

    Ruelle,Resonances in chaotic dynamical systems, Phys

    D. Ruelle,Resonances in chaotic dynamical systems, Phys. Rev. Lett.56, 405 (1986)

  21. [21]

    Prosen,Ruelle resonances in quantum many-body sys- tems, J

    T. Prosen,Ruelle resonances in quantum many-body sys- tems, J. Phys. A35, L737 (2002)

  22. [22]

    Braun,Dissipative Quantum Chaos and Decoherence, (Springer, 2003)

    D. Braun,Dissipative Quantum Chaos and Decoherence, (Springer, 2003)

  23. [23]

    A. N. Kolmogorov,The local structure of turbulence in incompressible viscous fluid for very large Reynolds num- bers, Dokl. Akad. Nauk SSSR30(1941); english transla- tion in Proc. R. Soc. Lond. A434, 9 (1991)

  24. [24]

    L. F. Richardson,Weather prediction by numerical pro- cesses(The University Press, Cambridge, 1922)

  25. [25]

    L. N. Trefethen and M. Embree,Spectra and pseudospec- tra(Princeton University press, 2005)

  26. [26]

    Mi et.al.,Information scrambling in quantum circuits, Science374, 1479 (2021)

    X. Mi et.al.,Information scrambling in quantum circuits, Science374, 1479 (2021)

  27. [27]

    Kim et al.,Evidence for the utility of quantum com- puting before fault tolerance, Nature618, 500 (2023)

    Y. Kim et al.,Evidence for the utility of quantum com- puting before fault tolerance, Nature618, 500 (2023)

  28. [28]

    Shtanko et

    O. Shtanko et. al.,Uncovering local integrability in quantum many-body dynamics, Nat. Commun.16, 2552 (2025)

  29. [29]

    Digital quantum magnetism on a trapped-ion quantum computer

    R. Haghshenas et al.,Digital quantum magnetism at the frontier of classical simulations,arXiv:2503.20870 (2025)

  30. [30]

    L. E. Fischer et al.,Dynamical simulations of many-body quantum chaos on a quantum computer, Nat. Phys.22, 302 (2026)

  31. [31]

    ˇZnidariˇ c,Momentum-dependent quantum Ruelle- Pollicott resonances in translationally invariant many- body systems, Phys

    M. ˇZnidariˇ c,Momentum-dependent quantum Ruelle- Pollicott resonances in translationally invariant many- body systems, Phys. Rev. E110, 054204 (2024)

  32. [32]

    There is of course in principle no contradiction between the unitarity ofUand exponential relaxation, keeping in mind (i) that one will have relaxation only for appropri- ateA, e.g. local observables, and (ii) one has to carefully take the correct TDL of first sendingL→ ∞and only thent→ ∞, and in this case exponential decay will indeed emerge out of an ex...

  33. [33]

    L. E. Ballentine,Quantum mechanics: A modern devel- opment(World Scientific, Singapore, 1998)

  34. [34]

    Bohm, and M

    A. Bohm, and M. Gadella,Direc Kets, Gamow Vectors and Gel’fand Triplets: The Rigged Hilbert Space Formu- lation of Quantum Mechanics(Springer-Verlga, Berlin, 1969)

  35. [35]

    de la Madrid,The role of the rigged Hilbert space in quantum mechanics, Eur

    R. de la Madrid,The role of the rigged Hilbert space in quantum mechanics, Eur. J. Phys.26287 (2005)

  36. [36]

    de la Madrid, and M

    R. de la Madrid, and M. Gadella,A pedestrian introduc- tion to Gamow vectors, Am. J. Phys.70, 626 (2002)

  37. [37]

    H. H. Hasegawa and W. C. Saphir,Unitarity and ir- reversibility in chaotic systems, Phys. Rev. A46, 7401 (1992)

  38. [38]

    Antoniou and S

    I. Antoniou and S. Tasaki,Spectral decomposition of the Reny maps, J. Phys. A26, 73 (1993)

  39. [39]

    Gaspard, G

    P. Gaspard, G. Nicolis, A. Provata and S. Tasaki,Spectral signature of the pitchfork bifurcation: Liouville equation approach,, Phys. Rev. E51, 74 (1995)

  40. [40]

    Antoniou, B

    I. Antoniou, B. Qiao, and Z. Suchanecki,General- ized spectral decomposition and intrinsic irreversibility of Arnold cat map, Chaos, Solitons & Fractals8, 77 (1997)

  41. [41]

    Prethermalization, shadowing breakdown, and the absence of Trotterization transition in quantum circuits

    M. ˇZnidariˇ c,Prethermalization, shadowing breakdown, and the absence of Trotterization transition in quantum circuits, Phys. Rev. Xin press(arXiv:2505.15521)

  42. [42]

    Duh and M

    U. Duh and M. ˇZnidariˇ c,Ruelle-Pollicott resonances of diffusive U(1)-invariant qubit circuits, SciPost Phys.20, 061 (2026)

  43. [43]

    Klobas,Exact time-dependent solutions of interacting systems, PhD Thesis, (University of Ljubljana, 2020)

    K. Klobas,Exact time-dependent solutions of interacting systems, PhD Thesis, (University of Ljubljana, 2020)

  44. [44]

    Sharipov, M

    R. Sharipov, M. Koterle, S. Grozdanov, and T. Prosen, Ergodic behaviors in reversible 3-state cellular automata, arXiv:2503.16593(2025)

  45. [45]

    ˇZnidariˇ c, U

    M. ˇZnidariˇ c, U. Duh, and L. Zadnik,Integrability is generic in homogeneous U(1)-invariant nearest-neighbor qubit circuits, Phys. Rev. B112, L020302 (2025)

  46. [46]

    ˇZnidariˇ c,Inhomogeneous SU(2) symmetries in homo- geneous integrable U(1) circuits and transport, Nat

    M. ˇZnidariˇ c,Inhomogeneous SU(2) symmetries in homo- geneous integrable U(1) circuits and transport, Nat. Com- mun.16, 4336 (2025)

  47. [47]

    Mori,Liouvillian-gap analysis of open quantum many- body systems in the weak dissipation limit, Phys

    T. Mori,Liouvillian-gap analysis of open quantum many- body systems in the weak dissipation limit, Phys. Rev. B 109, 064311 (2024)

  48. [48]

    J. A. Jacoby, D. A. Huse, and S. Gopalakrishnan,Spectral gaps of local quantum channels in the weak-dissipation limit, Phys. Rev. B111, 104303 (2025)

  49. [49]

    Zhang, L

    C. Zhang, L. Nie, and C. von Keyserlingk,Thermaliza- tion rates and quantum Ruelle-Pollicott resonances: in- sights from operator hydrodynamics,arXiv:2409.17251 (2024)

  50. [50]

    Yoshimura and L

    T. Yoshimura and L. S´ a,Theory of irreversibility in quantum many-body systems, Phys. Rev. E111, 064135 (2025)

  51. [51]

    Duarte, I

    J. Duarte, I. Garc´ ıa-Mata, D. A. Wisniacki,Ruelle- Pollicott decay of out-of-time-order correlators in many- body systems, Phys. Rev. E113, 024209 (2026). 19

  52. [52]

    Teretenkov, F

    A. Teretenkov, F. Uskov, and O. LychkovskiyPseu- domode expansion of many-body correlation functions, Phys. Rev. B111, 174308 (2025)

  53. [53]

    Rakovszky, C

    T. Rakovszky, C. W. von Keyserlingk, and F. Pollmann, Dissipation-assisted operator evolution method for cap- turing hydrodynamic transport, Phys. Rev. B105, 075131 (2022)

  54. [54]

    Beguˇ si´ c and G

    T. Beguˇ si´ c and G. K-L. Chan,Real-time operator evo- lution in two and three dimensions via sparse Pauli dy- namics, PRX Quantum6, 020302 (2025)

  55. [55]

    Schuster, C

    T. Schuster, C. Yin, X. Gao, and N. Y. Yao,A polynomial-time classical algorithm for noisy quantum circuits, Phys. Rev. X15, 041018 (2025)

  56. [56]

    Fontana, M

    E. Fontana, M. S. Rudolph, R. Duncan, I. Rungger, and C. Cˆ ırstoiu,Classical simulations of noisy variational quantum circuits, npj Quantum Inf.11, 84 (2025)

  57. [57]

    Angrisani, A

    A. Angrisani, A. Schmidhuber, M. S. Rudolph, M. Cerezo, Z. Holmes, and H.-Y. Huang,Classically Estimating Observables of Noiseless Quantum Circuits, Phys. Rev. Lett.135, 170602 (2025)

  58. [58]

    M. S. Rudolph, T. Jones, Y. Teng, A. Angrisani, and Z. HolmesPauli Propagation: A Computa- tional Framework for Simulating Quantum Systems, arXiv:2505.21606(2025)

  59. [59]

    Fritzsch, G

    F. Fritzsch, G. O. Alves, M. A. Rampp, and P. W. Claeys, Free Cumulants and Full Eigenstate Thermalization from Boundary Scrambling,arXiv:2509.08060(2025)

  60. [60]

    For instance, an operator with a densitya α 1 =σ x 1 σy 2 has s= 2, whileσ z 11 2σx 3 hass= 3

  61. [61]

    D. S. Watkins,Fundamentals of matrix computations (John Wiley & Sons, New York, 2002)

  62. [62]

    G. W. Stewart and J. Sun,Matrix perturbation theory (Academic Press, 1990)

  63. [63]

    The inequality holds for diagonalizable matrices and eigenvalues with an algebraic multiplicity 1. One can also define a condition number of a diagonalizable ma- trixM=W DW −1 (with diagonalD) byκ(M) := ∥W∥ p · ∥W −1∥p, and the Bauer-Fike theorem puts an upper bound on the spectral variation forVof any size, i.e., difference in any two eigenvalues ofMandM...

  64. [64]

    Petermann,Calculated spontaneous emission fac- tor for double-heterostructure injection lasers with gain- induced waveguiding, IEEE J

    K. Petermann,Calculated spontaneous emission fac- tor for double-heterostructure injection lasers with gain- induced waveguiding, IEEE J. Quant. Electron.15, 566 (1979)

  65. [65]

    Schomerus, K

    H. Schomerus, K. M. Frahm, M. Patra, and C. W. J. Beenakker,Quantum limit of the laser line width in chaotic cavities and statistics of residues of scattering matrix poles, Physica A278, 469 (2000)

  66. [66]

    D. V. Savin and V. V. Sokolov,Quantum versus classical decay laws in open chaotic systems, Phys. Rev. E56, R4911(R) (1998)

  67. [67]

    Evers and A

    F. Evers and A. D. Mirlin,Anderson transitions, Rev. Mod. Phys.80, 1355 (2008)

  68. [68]

    J. L. Kaplan and J. A. Yorke,Chaotic behavior of multidimensional difference equations, in H. O. Peit- gen, H. O. Walther (eds)Functional Differential Equa- tions and Approximation of Fixed Points, pp. 204-227 (Springer, Berlin, 1979)

  69. [69]

    J. D. Farmer, E. Ott, and J. A. Yorke,The dimension of chaotic attractors, Physica D7, 153 (1983)

  70. [70]

    Bertini, P

    B. Bertini, P. Kos, and T. Prosen,Exact Correlation Functions for Dual-Unitary Lattice Models in1 + 1Di- mensions, Phys. Rev. Lett.123, 210601 (2019)

  71. [71]

    Exactly solvable many-body dynamics from space-time duality

    B. Bertini, P. W. Claeys, and T. Prosen,Exactly solvable many-body dynamics from space-time duality, arXiv:2505.11489 (2025)

  72. [72]

    The truncated propagatorU r ink= 0 for DU circuits is actually equal to the non-translationally invariant prop- agator aligning left-moving operators to the left of the causal cone and right-moving operators to the right of the causal cone. On the level of the graph, the choice of the direction of the shiftSdoes not matter, since we are only interested in...

  73. [73]

    Holden-Dye, L

    T. Holden-Dye, L. Masanes, and A. Pal,Fundamen- tal charges for dual-unitary circuits, Quantum9, 1615 (2025)

  74. [74]

    3 is valid forr >4, but the structure of fixed points is the same also for smallerr

    The graph in Fig. 3 is valid forr >4, but the structure of fixed points is the same also for smallerr. All arguments are valid in the general case, largeris considered merely for simplicity. For the complete graph see Ref. [67]

  75. [75]

    The structure atr 0 must be calculated numerically, either by diagonalizingU r or the transfer matrices defined in Ref. [64]

  76. [76]

    For|λ 1|= 1, |R1⟩describes a conserved quantity localized inr 0

    At this point we assume that|λ 1|<1. For|λ 1|= 1, |R1⟩describes a conserved quantity localized inr 0. Triv- ially, there is no cascade in this case. For a more general discussion on the structure of conserved quantities see Ref. [67]

  77. [77]

    von Keyserlingk, F

    C. von Keyserlingk, F. Pollmann, and T. Rakovszky,Op- erator backflow and the classical simulation of quantum transport, Phys. Rev. B105, 245101 (2022)

  78. [78]

    C. W. von Keyserlingk, T. Rakovszky, F. Pollmann, and S. L. Sondhi,Operator Hydrodynamics, OTOCs, and Entanglement Growth in Systems without Conservation Laws, Phys. Rev. X8, 021013 (2018)

  79. [79]

    Single ring theorem

    J. Feinberg, R. Scalettar, and A. Zee,“Single ring theorem” and the disk-annulus phase transition, J. Math. Phys.42, 5718 (2001)

  80. [80]

    Guionnet, M

    A. Guionnet, M. Krishnapur, and O. Zeitouni,The single ring theorem,, Annals of Mathematics174, 1189 (2011)

Showing first 80 references.