pith. sign in

arxiv: 2606.18800 · v1 · pith:FILE4Q2Anew · submitted 2026-06-17 · ❄️ cond-mat.str-el · cond-mat.quant-gas

Mimicry of chaos and k-design in higher order OTOCs of Luttinger liquids

Bal\'azs D\'ora , Catalin Pascu Moca , Roderich Moessner This is my paper

Pith reviewed 2026-06-26 19:19 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.quant-gas
keywords out-of-time-order correlatorsLuttinger liquidXXZ chainHarper modelk-designbosonizationquantum chaosinformation scrambling
0
0 comments X

The pith

Late-time saturation values of higher-order OTOCs in Luttinger liquids map exactly onto the partition function of a non-Hermitian Harper model.

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

The paper computes the time evolution of the first three higher-order out-of-time-order correlators in a Luttinger liquid and its lattice realization, the XXZ chain. These correlators rapidly lose memory of the initial state and saturate to steady values. The central result is that the saturation values for the whole sequence of such correlators are given exactly by the partition function of a non-Hermitian Harper model obtained via bosonization. For moderate interaction strengths this mapping yields parametrically small steady-state values through seventh order, which mimics the statistics of higher k-designs. A sympathetic reader would care because the result supplies an exact, non-perturbative handle on the fine structure of apparent scrambling in an integrable model.

Core claim

In a Luttinger liquid the late-time saturation values for the entire sequence of higher-order OTOCs are given exactly by the partition function of a non-Hermitian Harper model. The same mapping holds in the XXZ Heisenberg chain. Through this equivalence the steady-state OTOCs become parametrically small up to seventh order when the interaction is moderately strong, thereby reproducing the diagnostic signature of higher k-designs.

What carries the argument

exact mapping of late-time saturation values of the sequence of higher-order OTOCs onto the partition function of a non-Hermitian Harper model

If this is right

  • The full tower of higher-order OTOC saturation values is obtainable from a single Harper-model partition function.
  • For moderate interactions the OTOCs remain small through at least seventh order.
  • The same Harper-model equivalence applies to both the continuum Luttinger liquid and its lattice XXZ realization.
  • Higher-order OTOCs therefore furnish a finer diagnostic of apparent scrambling than the usual second-order correlator.

Where Pith is reading between the lines

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

  • The Harper-model mapping may extend to other integrable 1D systems whose bosonization yields similar quadratic forms.
  • If the mapping survives weak integrability-breaking perturbations, it could supply an analytic window into the crossover from apparent to genuine higher-order scrambling.
  • Exact diagonalization or tensor-network calculations of higher-order OTOCs in larger XXZ chains would provide an independent test of the saturation values.

Load-bearing premise

The bosonization procedure and the identification of late-time OTOC values with the Harper-model partition function continue to hold in the interaction regime where the correlators become parametrically small up to seventh order.

What would settle it

A numerical evaluation of the seventh-order OTOC saturation value on the XXZ chain at moderate interaction strength that differs substantially from the value predicted by the non-Hermitian Harper partition function would falsify the claimed mapping.

Figures

Figures reproduced from arXiv: 2606.18800 by Bal\'azs D\'ora, Catalin Pascu Moca, Roderich Moessner.

Figure 1
Figure 1. Figure 1: FIG. 1. The time evolution of the [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The long time limit of the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
read the original abstract

Out-of-time-order correlators (OTOCs) provide a fundamental metric for quantum chaos, but capturing the fine structure of information scrambling requires exploring their higher-order generalizations. Here, we systematically investigate the sequence of higher-order OTOCs in a Luttinger liquid and its lattice realization, the XXZ Heisenberg chain. Using bosonization and numerics, we extract the full temporal dynamics of the first three OTOCs, revealing that they rapidly erase memory of the initial state, and quickly saturate to their steady state values. Strikingly, we show that calculating the late time saturation values for the entire sequence of higher-order OTOCs maps exactly onto determining the partition function of a non-Hermitian Harper model. Through this mapping, we demonstrate that for moderately strong interactions, the steady-state OTOCs become parametrically small up to the seventh order, mimicking higher $k$-design. Our results reveal that Luttinger liquids exhibit an unexpectedly profound degree of apparent scrambling when viewed through the lens of higher-order OTOCs.

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

2 major / 2 minor

Summary. The paper investigates higher-order out-of-time-order correlators (OTOCs) in Luttinger liquids and their lattice realization (XXZ chain). Using bosonization and numerics, it computes the dynamics of the first three OTOCs, finds rapid saturation to steady-state values, and asserts that the late-time saturation values for the full sequence of higher-order OTOCs map exactly onto the partition function of a non-Hermitian Harper model. For moderately strong interactions, these values become parametrically small up to seventh order, which the authors interpret as mimicry of higher k-designs.

Significance. If the exact mapping holds without uncontrolled corrections from the bosonization approximations, the result would provide a concrete link between higher-order scrambling diagnostics in an integrable 1D model and the partition function of a non-Hermitian Harper model, offering a new perspective on apparent chaos in Luttinger liquids. The use of both analytic bosonization and numerics for low-order cases is a positive feature.

major comments (2)
  1. [mapping section] § on the mapping (around the paragraph asserting exact equivalence): the identification of late-time OTOC saturation values with the Harper-model partition function is presented as exact, yet the derivation relies on the standard bosonization dictionary (linear dispersion, vertex-operator mapping). No explicit error bound or regime of validity is given for the interaction strengths where the OTOCs become parametrically small up to order 7; any O(1) correction from nonlinear dispersion or higher-gradient terms would invalidate the exact mapping.
  2. [numerical results section] Abstract and § on numerical results: the claim that OTOCs saturate rapidly and become parametrically small up to seventh order for moderate interactions rests on the bosonization procedure remaining controlled in that regime. The manuscript provides no quantitative estimate (e.g., comparison of linear vs. nonlinear dispersion contributions) showing that the suppression survives beyond the low-energy approximation.
minor comments (2)
  1. [figures] Figure captions for the OTOC time series should explicitly state the interaction parameter values (K or Δ) used and the system sizes for the numerics.
  2. [definitions] Notation for the higher-order OTOCs (e.g., the precise operator ordering and normalization) should be defined once in a dedicated subsection before the mapping is introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments, which help clarify the scope of our results. We address the two major comments point by point below. We will revise the manuscript to better delineate the regime of validity of the bosonization mapping while preserving the exactness of the mapping inside that effective theory.

read point-by-point responses

  1. Referee: [mapping section] § on the mapping (around the paragraph asserting exact equivalence): the identification of late-time OTOC saturation values with the Harper-model partition function is presented as exact, yet the derivation relies on the standard bosonization dictionary (linear dispersion, vertex-operator mapping). No explicit error bound or regime of validity is given for the interaction strengths where the OTOCs become parametrically small up to order 7; any O(1) correction from nonlinear dispersion or higher-gradient terms would invalidate the exact mapping.

    Authors: The mapping between the late-time saturation values of the higher-order OTOCs and the partition function of the non-Hermitian Harper model is exact within the Luttinger-liquid effective theory (linear dispersion, standard vertex-operator correspondence). This is the same controlled approximation routinely used for the low-energy physics of the XXZ chain. We agree, however, that the manuscript should explicitly state the regime of validity and note that corrections from nonlinear dispersion or umklapp terms are parametrically small at low energies. We will add a dedicated paragraph in the mapping section that (i) recalls the standard validity conditions of bosonization for the XXZ model and (ii) argues why, for the moderate interaction strengths considered, higher-gradient corrections remain sub-leading for the saturation values up to order 7. revision: yes

  2. Referee: [numerical results section] Abstract and § on numerical results: the claim that OTOCs saturate rapidly and become parametrically small up to seventh order for moderate interactions rests on the bosonization procedure remaining controlled in that regime. The manuscript provides no quantitative estimate (e.g., comparison of linear vs. nonlinear dispersion contributions) showing that the suppression survives beyond the low-energy approximation.

    Authors: The rapid saturation and the numerical values for the first three OTOCs are obtained from lattice numerics (exact diagonalization / DMRG) on the XXZ chain and do not rely on bosonization. The parametric smallness up to seventh order follows from the exact mapping inside the Luttinger-liquid theory. We will revise the abstract and the numerical-results section to (i) clearly separate the lattice results from the effective-theory predictions and (ii) include a brief discussion, supported by existing literature on the XXZ model, of the energy scales at which nonlinear corrections become appreciable. A fully quantitative error estimate for every order would require new calculations outside the present scope; we therefore provide a qualitative but explicit statement of the regime rather than a numerical bound. revision: partial

Circularity Check

0 steps flagged

No circularity: bosonization mapping to Harper partition function is independently derived

full rationale

The central claim equates late-time OTOC saturation values to the partition function of a non-Hermitian Harper model via standard bosonization applied to the Luttinger liquid/XXZ chain. This is presented as an exact mapping obtained from the effective theory and numerics, without any indication that the Harper partition function was fitted to OTOC data or that the result reduces to a self-definition, renamed known result, or self-citation chain. Bosonization is an external low-energy technique whose validity is an assumption (addressed by the skeptic as a correctness issue, not circularity). No load-bearing step reduces by construction to the target OTOC sequence; the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of bosonization for higher-order OTOCs and on the exact identification of late-time values with the Harper-model partition function; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Bosonization correctly captures the higher-order OTOC dynamics in the Luttinger liquid for the interaction strengths considered
    Invoked to extract the full temporal dynamics and the saturation values

pith-pipeline@v0.9.1-grok · 5730 in / 1353 out tokens · 25447 ms · 2026-06-26T19:19:35.438590+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

43 extracted references · 2 canonical work pages

  1. [1]

    We have checked numerically that by using the matrix representation[37] XandPfor the Harper model, its moments reproduce Eq

    All of its moments are real due to the self duality of the non-hermitian Harper Hamiltonian. We have checked numerically that by using the matrix representation[37] XandPfor the Harper model, its moments reproduce Eq. (10), (11) and (13). Numerics.These analytical results forλ= 2 are com- pared with many-body exact diagonalization of the XXZ Heisenberg ch...

    2023

  2. [2]

    Kitaev,A simple model of quantum holog- raphy, online.kitp.ucsb.edu/online/entangled15/kitaev/, online.kitp.ucsb.edu/online/entangled15/kitaev2/

    A. Kitaev,A simple model of quantum holog- raphy, online.kitp.ucsb.edu/online/entangled15/kitaev/, online.kitp.ucsb.edu/online/entangled15/kitaev2/. Talks at KITP, April 7, 2015 and May 27, 2015

    2015

  3. [3]

    Kitaev and S

    A. Kitaev and S. J. Suh,The soft mode in the sachdev- ye-kitaev model and its gravity dual, J. High Energ. Phys. 2018, 183 (2018)

    2018

  4. [4]

    Garc´ ıa-Mata, R

    I. Garc´ ıa-Mata, R. A. Jalabert, and D. A. Wisni- acki,Out-of-time-order correlations and quantum chaos, Scholarpedia18(4), 55237 (2023), revision #204529

    2023

  5. [5]

    Xu and B

    S. Xu and B. Swingle,Scrambling dynamics and out-of- time-ordered correlators in quantum many-body systems, PRX Quantum5, 010201 (2024)

    2024

  6. [6]

    X. Mi, P. Roushan, C. Quintana, S. Mandr` a, J. Marshall, C. Neill, F. Arute, K. Arya, J. Atalaya, R. Babbush, J. C. Bardin, R. Barends,et al.,Information scrambling in quantum circuits, Science374(6574), 1479 (2021), https://www.science.org/doi/pdf/10.1126/science.abg5029

    work page doi:10.1126/science.abg5029 2021

  7. [7]

    K. A. Landsman, C. Figgatt, T. Schuster, N. M. Linke, B. Yoshida, N. Y. Yao, and C. Monroe,Verified quantum information scrambling, Nature567, 61 (2019)

    2019

  8. [8]

    A. I. Larkin and Y. N. Ovchinnikov,Quasiclassical method in the theory of superconductivity, Sov. Phys. JETP28, 1200 (1969)

    1969

  9. [9]

    D. A. Roberts and B. Swingle,Lieb-robinson bound and the butterfly effect in quantum field theories, Phys. Rev. Lett.117, 091602 (2016)

    2016

  10. [10]

    Stanford,Many-body chaos at weak coupling, Journal of High Energy Physics2016, 9 (2016)

    D. Stanford,Many-body chaos at weak coupling, Journal of High Energy Physics2016, 9 (2016)

    2016

  11. [11]

    I. L. Aleiner, L. Faoro, and L. B. Ioffe,Microscopic model of quantum butterfly effect: Out-of-time-order correlators and traveling combustion waves, Annals of Physics375, 378 (2016)

    2016

  12. [12]

    G¨ arttner, J

    M. G¨ arttner, J. G. Bohnet, A. Safavi-Naini, M. L. Wall, J. J. Bollinger, and A. M. Rey,Measuring out- of-time-order correlations and multiple quantum spectra in a trapped ion quantum magnet, Nature Phys13, 781 (2017)

    2017

  13. [13]

    J. Li, R. Fan, H. Wang, B. Ye, B. Zeng, H. Zhai, X. Peng, and J. Du,Measuring out-of-time-order correlators on a nuclear magnetic resonance quantum simulator, Phys. Rev. X7, 031011 (2017)

    2017

  14. [14]

    Maldacena and D

    J. Maldacena and D. Stanford,Remarks on the sachdev- ye-kitaev model, Phys. Rev. D94, 106002 (2016)

    2016

  15. [15]

    G. Zhu, M. Hafezi, and T. Grover,Measurement of many- body chaos using a quantum clock, Phys. Rev. A94, 062329 (2016)

    2016

  16. [16]

    C. W. von Keyserlingk, T. Rakovszky, F. Pollmann, and S. L. Sondhi,Operator hydrodynamics, otocs, and en- tanglement growth in systems without conservation laws, Phys. Rev. X8, 021013 (2018)

    2018

  17. [17]

    Maldacena, S

    J. Maldacena, S. H. Shenker, and D. Stanford,A bound on chaos, Journal of High Energy Physics2016, 106 (2016)

    2016

  18. [18]

    Bohrdt, C

    A. Bohrdt, C. B. Mendl, M. Endres, and M. Knap, Scrambling and thermalization in a diffusive quantum many-body system, New Journal of Physics19(6), 063001 (2017)

    2017

  19. [19]

    D. A. Roberts and D. Stanford,Diagnosing chaos using four-point functions in two-dimensional conformal field theory, Phys. Rev. Lett.115, 131603 (2015)

    2015

  20. [20]

    D. A. Roberts and B. Yoshida,Chaos and complexity by design, J. High Energ. Phys.2017, 121 (2017)

    2017

  21. [21]

    Google Quantum AI and Collaborators,Observation of constructive interference at the edge of quantum ergodic- 6 ity, Nature646, 825 (2025)

    2025

  22. [22]

    Iaconis,Quantum state complexity in computation- ally tractable quantum circuits, PRX Quantum2, 010329 (2021)

    J. Iaconis,Quantum state complexity in computation- ally tractable quantum circuits, PRX Quantum2, 010329 (2021)

    2021

  23. [23]

    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

    work page arXiv

  24. [24]

    Giamarchi,Quantum Physics in One Dimension(Ox- ford University Press, Oxford, 2004)

    T. Giamarchi,Quantum Physics in One Dimension(Ox- ford University Press, Oxford, 2004)

    2004

  25. [25]

    A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, Bosonization and Strongly Correlated Systems(Cam- bridge University Press, Cambridge, 1998)

    1998

  26. [26]

    M. A. Cazalilla,Bosonizing one-dimensional cold atomic gases, J. Phys. B: At. Mol. Opt. Phys.37, S1 (2004)

    2004

  27. [27]

    Nielsen and I

    M. Nielsen and I. Chuang,Quantum Computation and Quantum Information(Cambridge University Press, Cambridge, 2000)

    2000

  28. [28]

    M. Fava, J. Kurchan, and S. Pappalardi,Designs via free probability, Phys. Rev. X15, 011031 (2025)

    2025

  29. [29]

    von Delft and H

    J. von Delft and H. Schoeller,Bosonization for beginners - refermionization for experts, Ann. Phys. (Leipzig)7, 225 (1998)

    1998

  30. [30]

    Luther and I

    A. Luther and I. Peschel,Calculation of critical expo- nents in two dimensions from quantum field theory in one dimension, Phys. Rev. B12, 3908 (1975)

    1975

  31. [31]

    Iucci, G

    A. Iucci, G. A. Fiete, and T. Giamarchi,Fourier trans- form of the2k F luttinger liquid density correlation func- tion with different spin and charge velocities, Phys. Rev. B75, 205116 (2007)

    2007

  32. [32]

    10c) in terms of cos(2c) when comparing to numerics

    We use cos(6c) = 4 cos3(2c)−3 cos(2c) to rewrite all har- monics but the highest (i.e. 10c) in terms of cos(2c) when comparing to numerics

  33. [33]

    P. G. Harper,Single band motion of conduction electrons in a uniform magnetic field, Proceedings of the Physical Society. Section A68(10), 874 (1955)

    1955

  34. [34]

    Aubry and G

    A. Aubry and G. Andre,Analyticity breaking and ander- son localization in incommensurate lattices, Ann. Israel. Phys. Soc.3, 133 (1980)

    1980

  35. [35]

    A. C. Quillen, N. Skerrett, D. R. Sowinski, and A. S. Mi- akhel,Notions of adiabatic drift in the quantized harper model, Phys. Rev. A112, 042226 (2025)

    2025

  36. [36]

    A. G. Moghaddam, V. K¨ onye, L. Mertens, J. van Wezel, and J. van den Brink,Synthetic horizons in an atomic chain: Horizon-induced effects and connections to quan- tum hall systems, Phys. Rev. B113, 195430 (2026)

    2026

  37. [37]

    Q.-R. Xu, E. Cobanera, and G. Ortiz,Bloch and bethe ans¨ atze for the harper model: A butterfly with a bound- ary, Phys. Rev. B104, 165140 (2021)

    2021

  38. [38]

    H. J. Korsch and K. Rapedius,Computations in quan- tum mechanics made easy, European Journal of Physics 37(5), 055410 (2016)

    2016

  39. [39]

    Lukyanov,Correlation amplitude for the XXZ spin chain in the disordered regime, Phys

    S. Lukyanov,Correlation amplitude for the XXZ spin chain in the disordered regime, Phys. Rev. B59, 11163 (1999)

    1999

  40. [40]

    D´ ora and R

    B. D´ ora and R. Moessner,Out-of-time-ordered density correlators in luttinger liquids, Phys. Rev. Lett.119, 026802 (2017)

    2017

  41. [41]

    Lin and O

    C.-J. Lin and O. I. Motrunich,Out-of-time-ordered cor- relators in short-range and long-range hard-core boson models and in the luttinger-liquid model, Phys. Rev. B 98, 134305 (2018)

    2018

  42. [42]

    Gisti, D

    M. Gisti, D. J. Luitz, and M. Debertolis,Symmetry re- solved out-of-time-order correlators of Heisenberg spin chains using projected matrix product operators, Quan- tum9, 1871 (2025), ISSN 2521-327X

    2025

  43. [43]

    Hosur, X.-L

    P. Hosur, X.-L. Qi, D. A. Roberts, and B. Yoshida, Chaos in quantum channels, J. High Energ. Phys2016, 4 (2016)

    2016