{"id":"36b5246b-b836-46db-bac6-fba291cae53e","arxiv_id":"2411.16274","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In a random-matrix model of a chaotic many-body system, the late-time OTOC decays as a Gaussian with time scale ℏ/(2Δ); with the conjecture Δ = ℏλ_max, this is governed by λ_max t.","lead":"A theoretical physicist calculated how the 'out-of-time ordered correlator', a quantum probe of chaos, decays at late times in a random-matrix model of a chaotic many-body system. The decay is governed by the energy scale where quantum level statistics become universal; the paper conjectures this scale equals Planck's constant times the classical Lyapunov exponent.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The advertised universal λ_max t dependence rests on the unsupported conjecture Δ=ℏλ_max (Eq. 31), not on the calculation itself; until that link is independently tested, the central physical claim remains unestablished.","rationale":"I read the paper in good faith as an explicit calculation within a postulated random-matrix ensemble, with an honest statement that the connection to λ_max is a conjecture. The reader's weakest_assumption is the Gaussian form of F in Eq. (15), and that is indeed a stated limitation of the model, not derived from microscopic dynamics. However, I judge the most load-bearing concern to be the conjecture Eq. (31), Δ=ℏλ_max. That single equation is what turns a calculation about an auxiliary parameter Δ into the advertised result that late-time OTOC decay is a universal probe of the classical Lyapunov exponent. It has essentially one supporting datum (Ref. [8], 400 Sinai-billiard levels) and is explicitly in tension with Ref. [9]'s exponential decay, which the paper does not resolve. Replacing the Gaussian in Eq. (15) with a Lorentzian would alter the decay shape but would not destroy the λ_max t time-scaling once Eq. (31) is assumed; failing Eq. (31) removes the physical interpretation entirely. Because the paper labels this as conjectural and the reader's CONDITIONAL verdict already accounts for such an unproved link, I do not see a reason to move the verdict. The concrete test above would settle the most important open question. This is a disagreement about which weakness deserves the primary weight, not a rejection of the reader's overall assessment.","tokens_in":10004,"tokens_out":6481,"duration_ms":65222,"concrete_test":"Perform a modern high-statistics test of Eq. (31) on the Sinai billiard, the system underlying Ref. [8]: compute ~10^4 consecutive levels in the semiclassical regime, extract the local-GOE correlation width Δ from a fit to the two-point spectral correlation function (or eigenvector component correlations), independently compute the leading classical Lyapunov exponent λ_max, and test whether Δ=ℏλ_max holds within numerical uncertainty. If the equality fails, the central λ_max t claim is not supported; if it holds, the remaining burden is the Gaussian form of F in Eq. (15).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytic derivation yields a self-contained statement about Δ: Eqs. (27) and (28) give ensemble-averaged OTOC decaying as exp(-2t²Δ²/ℏ²), with Δ the local-GOE correlation width. The abstract's headline conclusion, that late-time OTOC is governed by λ_max t, requires Eq. (31), Δ=ℏλ_max. That equation is only a conjecture, based on Arve's 1991 analysis of 400 Sinai-billiard levels (Ref. [8]); it is not derived from the model of Section III. The paper itself notes the dissonance with Ref. [9], which finds exponential OTOC decay at the Ruelle-Pollicott rate, and does not resolve whether that system satisfies the local-GOE assumptions. This is the most load-bearing link because removing it collapses the paper's advertised physical message: even a perfect calculation of OTOC in terms of Δ would say nothing about Lyapunov exponents. The reader's Gaussian-vs-Lorentzian concern about F in Eq. (15) is real but secondary: a Lorentzian F would change the decay shape (to exponential in t) while preserving the time scale ℏ/Δ, so conditional on Eq. (31) the λ_max t scaling and the physical probe interpretation would survive the shape change.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers a time-reversal-invariant chaotic many-body quantum system described by a Hartree-Fock term plus a residual interaction, and models its Hamiltonian with a random-matrix ensemble in which eigenvectors and eigenvalues obey local GOE statistics over an energy correlation width Δ. Using a non-crossing contraction expansion for moments of the evolution operator, the author computes the ensemble-averaged out-of-time-ordered correlator C(t) and the symmetrized function F(t) to leading order in 1/N, obtaining a Gaussian time decay exp(-2t^2Δ^2/ℏ^2) and a static asymptotic term for C(t). The author also bounds the variance by 1/N and conjectures Δ = ℏλ_max, so that the late-time OTOC is governed by the dimensionless parameter λ_max t.","tokens_in":10308,"tokens_out":17618,"duration_ms":165044,"significance":"If the calculation is correct, the paper provides an explicit analytic formula for the late-time OTOC in a random-matrix model of many-body chaos and identifies a concrete time scale, ℏ/Δ, set by the Bohigas-Giannoni-Schmit correlation width. The moment expansion in Section V and the 1/N bookkeeping are carefully organized, and the variance estimate in Section VIII is a useful self-averaging result. The author is also transparent about two major caveats: the Gaussian form of the correlation function in Eq. (15) is based on numerical evidence, and the connection to Lyapunov exponents through Eq. (31) is a conjecture. However, the internal inconsistency in Eq. (20), the very strong assumption in Eq. (26), and the mismatch between the variance bound and the almost-everywhere claim in Eq. (30) mean that the central results are not yet established as stated.","major_comments":[{"comment":"The printed connected k-th moment has prefactor (√(2π)ρ(E_{m_1})Δ)^{k-1} and a positive exponent +Σ_{j<l}(E_{m_j}-E_{m_l})^2/(2kΔ^2). These two features are internally inconsistent with the equations that are said to follow from Eq. (20): inserting k=2 and χ_1=χ_2=-β gives Eq. (24) only if the prefactor is inverted and the pairwise exponent has a minus sign, and the same negative Gaussian factors appear in Eqs. (27) and (33). As printed, the moment grows with energy separation and has the wrong 1/N scaling, so the derivation of the central results is not self-consistent.","section":"Section V, Eq. (20)"},{"comment":"The 'slightly more stringent' assumption is not slight. For a dense spectrum, requiring Σ_m A_{mm}e^{-βE_m}e^{-(E_m-E_n)^2/(2kΔ^2)}=0 for every E_n, β, and positive integer k forces the measure ρ(E_m)A_{mm}e^{-βE_m} to vanish identically, because the Gaussian factor has a positive Fourier transform. Thus V and W are effectively required to have vanishing diagonal matrix elements in the Hartree-Fock basis. The paper gives no argument that the simple operators of Ref. [4] satisfy this condition, and without Eq. (26) the terms discarded after possibility (ii) in Section VII contribute at an uncontrolled order to Eqs. (27) and (28).","section":"Section VII, Eq. (26)"},{"comment":"The abstract's headline statement that the large-time OTOC is governed by λ_max t is not a consequence of the calculation: the derivation produces a decay with time scale ℏ/Δ, and the identification Δ = ℏλ_max is a conjecture imported from Ref. [8] and not derived within the ensemble of Section III. The paper itself notes the discrepancy with Ref. [9] and does not resolve it. The physical claim should be explicitly marked as conditional on Eq. (31), or Eq. (31) should be supported by independent evidence.","section":"Section IX, Eq. (31)"},{"comment":"Section VIII bounds the variance at fixed t by order 1/N relative to the squared mean. This yields convergence in probability for each fixed t, not the assertion in Eq. (30) that F(t)=⟨F(t)⟩ and C(t)=⟨C(t)⟩ for almost all members of the ensemble. To justify an almost-everywhere statement uniformly in t, the paper needs a continuity or tightness argument, or it should weaken Eq. (30) to a statement about ensemble averages and fixed-time fluctuations.","section":"Section VIII and Eq. (30)"},{"comment":"The explicit Gaussian factor exp(-2t^2Δ^2/ℏ^2) is the Fourier transform of the assumed Gaussian model for F, and the paper concedes that this form is based on numerical evidence and may fail for weak chaos. Since the Gaussian shape is the central quantitative prediction, the paper should either present Eq. (15) as a model assumption whose consequences are conditional, or repeat the calculation for a generic correlation function of width Δ in order to separate the robust time scale ℏ/Δ from the shape-dependent decay.","section":"Section IX and Eq. (15)"}],"minor_comments":[{"comment":"Equation (14) uses F(E_m - E_α) while Eq. (15) defines F(E_m - \\bar E_α); the overline should appear consistently in both places.","section":"Section III, Eqs. (14) and (15)"},{"comment":"The word 'Ljapunov' should be 'Lyapunov' (abstract and Section IX), and there are typographical errors such as 'advantange', 'occuring', and 'diffult' that should be corrected.","section":"Throughout"},{"comment":"Equation (33) has a missing brace in the exponential: it should read exp{-(E_{m_1}-E_{m_3})^2/(4Δ^2)} rather than exp -(E_{m_1}-E_{m_3})^2/(4Δ^2)}.","section":"Appendix, Eq. (33)"},{"comment":"The notation m_{j+1} with j=k implicitly means m_1; this cyclic convention should be stated explicitly when the product of Kronecker deltas is introduced.","section":"Section V, Eq. (20)"},{"comment":"The statement that the factors formally commute but their order must be respected is clear, but the claim that only one connected contraction pattern survives for general k would be easier to verify with an explicit k=3 example illustrating the non-crossing rule.","section":"Section V, Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"This is a compact single-author paper built around a model calculation plus a conjecture. The core calculation is promising but currently contains an internal inconsistency in Eq. (20) and a very strong assumption in Eq. (26), both of which are fixable. The journal's readership may also want the abstract to clearly separate the derived Δ-dependent statement from the conjectured λ_max identification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is the explicit ensemble-averaged OTOC for the parametric random-matrix model of Ref. [3]. That calculation is new: Ref. [4] only estimated the asymptotic value, and the Gaussian factor exp(-2t²Δ²/ℏ²) in Eqs. (27) and (28) is worked out rather than guessed. The contraction-counting in Section V is careful, the 1/N suppression of the variance in Section VIII is plausible, and the paper is candid about what it can and cannot do. Credit where due: this is a legitimate analytic result about a specific many-body random-matrix ensemble, not a hollow claim. The soft spots are real but they are not hidden. The Gaussian time decay is inherited from the assumed Gaussian form of F in Eq. (15). The author says so himself and notes that a Lorentzian would change the decay shape while preserving the ℏ/Δ time scale. That is a fair caveat. The more load-bearing issue is the conjecture Δ = ℏλ_max in Eq. (31). The abstract's physical punchline—late-time OTOC governed by λ_max t—rests entirely on that conjecture, which is supported only by Arve's 400 Sinai-billiard levels and is in tension with Ref. [9]. The stress-test note is right about this. But the paper does not dress the conjecture as a theorem; it is flagged as a conjecture, and the conflict with Ref. [9] is acknowledged rather than ignored. I do not think there is a load-bearing flaw in the calculation itself. If you drop Eq. (31), you still have a well-defined statement about the correlation width Δ, and the shape of the decay follows from the assumed F. The math in the moments section looks consistent to me, and the treatment of Tr(Z) fluctuations is reasonable for N ≫ 1. The main risk is not circularity—it is that the universal Lyapunov claim may be empirically wrong or restricted to a narrower class of systems than the BGS conjecture covers. That is a scientific question worth refereeing, not a reason to desk-reject. Who should read this: anyone working on OTOC in random-matrix or chaotic many-body settings, and anyone interested in the BGS conjecture's quantitative content. It deserves a serious referee—the calculation should be checked, the conjecture should be challenged, and the relation to Ref. [9] should be sharpened. My own verdict would be conditional, not accept, but this is exactly the kind of paper peer review exists for.","headline":"A careful, honest random-matrix calculation of late-time OTOC decay, with the advertised λ_max connection explicitly labeled as a conjecture rather than a derived result.","tokens_in":743,"tokens_out":4358,"would_cite":true,"duration_ms":45585,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q50","15B52"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that for a chaotic many-body quantum system built on the Bohigas–Giannoni–Schmit random-matrix ensemble, the out-of-time ordered correlator decays at large times as a Gaussian with width ℏ/(2Δ), and conjectures that Δ…","keywords":["out-of-time ordered correlator","quantum chaos","random matrix theory","Bohigas-Giannoni-Schmit conjecture","Lyapunov exponent","Gaussian orthogonal ensemble","many-body quantum system","energy correlation width"],"falsifier":"Measure the two-point energy correlation function of a chaotic quantum system (for instance the Sinai billiard or a kicked rotor) in the semiclassical regime and compare it to a Gaussian; if the correlation function is measurably non-Gaussian, or if the measured large-time OTOC decay is exponential rather than Gaussian with the predicted $\\hbar/(2\\Delta)$ width, the central claim fails. A cheaper test: check whether the conjectured relation $\\Delta = \\hbar \\lambda_{\\max}$ holds by computing $\\Delta$ from spectral statistics and $\\lambda_{\\max}$ from the classical dynamics of the same system.","tokens_in":9786,"feed_emoji":"🎲","tokens_out":7291,"duration_ms":100478,"temperature":0.7,"pith_summary":"This paper tries to pin down the large-time behavior of the out-of-time ordered correlator (OTOC) for a chaotic many-body quantum system, using a random-matrix parametrization of the Hamiltonian that implements the Bohigas–Giannoni–Schmit conjecture. The central result is that the ensemble-averaged OTOC decays as exp(-2t²Δ²/ℏ²), so the decay time is ℏ/(2Δ), where Δ is the width of the energy interval in which Gaussian orthogonal ensemble statistics holds. The paper further conjectures that Δ is universally related to the largest Lyapunov exponent of the corresponding classical system by Δ = ℏλ_max. If that conjecture is right, the late-time OTOC decay is governed by the dimensionless parameter λ_max t, making OTOC a direct quantum probe of classical chaos.","feed_headline":"Gaussian OTOC decay ties quantum chaos to Lyapunov time","feed_subtitle":"Late-time OTOC decay runs on ℏ/Δ; if Δ=ℏλ_max, it directly exposes classical chaos.","key_machinery":"The central machinery is the universal parametrization of the chaotic Hamiltonian: $H_{mn} = \\sum_\\alpha O_{m\\alpha} E_\\alpha O_{n\\alpha}$, with zero-centered Gaussian orthogonal-matrix elements satisfying $\\langle O_{m\\alpha} O_{n\\beta} \\rangle = \\delta_{mn} \\delta_{\\alpha\\beta} F(E_m - \\bar{E}_\\alpha)$ and a normalized Gaussian $F$ of width $\\Delta$; within $\\Delta$ the eigenvalues obey Wigner-Dyson statistics. The argument is carried by the correlated-moment formula for products of $Y(\\chi) = \\exp(\\chi H)$, which to leading order in $1/N$ keeps only non-crossing contraction patterns and yields explicit exponentials in the energy differences and in $\\chi^2 \\Delta^2/(2k)$. Applying that formula to the four $Y$-insertions in $C(t)$ and $F(t)$, and using the assumption that the test operators have vanishing thermal averages over a width-$\\Delta$ window, reduces each average to three terms, all sharing the time factor $\\exp(-2t^2\\Delta^2/\\hbar^2)$.","core_discovery":"For a time-reversal-invariant chaotic many-body system whose Hamiltonian H = H_HF + V is represented through random orthogonal matrices O and eigenvalues E with Wigner-Dyson statistics inside a correlation width Δ, the paper shows that the thermal OTOC C(t) and the symmetrized function F(t) are self-averaging: their variances are of order 1/N relative to the squared means, so almost every member of the ensemble exhibits the average behavior. The averages decay at large times as Gaussians: ⟨C(t)⟩ carries the factor exp(-2t²Δ²/ℏ²) and tends to a nonvanishing constant built from thermal averages of V² and W², while ⟨F(t)⟩ vanishes with the same Gaussian factor times exp(-3β²Δ²/8). Since the Gaussian time factor comes by Fourier transformation from the assumed Gaussian energy correlation function, the time scale ℏ/Δ is robust to the precise form of that function. The paper conjectures that Δ = ℏλ_max, which would make the late-time OTOC decay universally controlled by the classical Lyapunov exponent.","pith_inferences":["If the Gaussian decay is confirmed, the Lorentzian alternative for weak chaos suggests a possible order-parameter-like distinction between strong and weak quantum chaos, with the shape of the energy correlation function determining the OTOC decay law.","The approach is restricted to time-reversal-invariant (GOE) systems; an analogous treatment for GUE or GSE ensembles might yield a different numerical coefficient in the Gaussian exponent, worth testing.","The conjectured equality $\\Delta = \\hbar \\lambda_{\\max}$, if combined with the BGS interval, could provide a way to extract Lyapunov exponents from spectral statistics without computing time evolution at all.","The discrepancy with exponential Ruelle–Pollicott decay in finite-size matrix models might be resolved by checking whether those models obey the BGS conjecture within a width $\\Delta$; if they do not, the Gaussian prediction would not apply to them."],"forward_implications":["If $\\Delta = \\hbar \\lambda_{\\max}$ is confirmed, the late-time OTOC decay directly measures the classical Lyapunov exponent in any chaotic many-body system that satisfies the Bohigas–Giannoni–Schmit conjecture.","The decay time scale $\\hbar/(2\\Delta)$ is independent of the detailed form of the energy correlation function, so it should be a universal signature of strong quantum chaos.","The self-averaging property (variances of order $1/N$) means a single realization of a chaotic many-body system, not only an ensemble average, should show the predicted Gaussian decay.","The asymptotic value of $\\langle C(t)\\rangle$ gives a relation between thermal averages of $V^2$ and $W^2$, a prediction checkable in numerical simulations.","The conjecture extends the Bohigas–Giannoni–Schmit conjecture by adding a quantitative link between the spectral correlation width and the classical dynamics."],"supporting_citations":[{"why":"Supplies the universal parametrization of the chaotic Hamiltonian and the Gaussian form of the energy correlation function $F$, which is the input for the whole calculation.","marker":"[3]"},{"why":"Provides the definitions of $C(t)$ and $F(t)$, the $\\beta/4$ insertion trick, and the expected asymptotic behavior that the paper computes explicitly.","marker":"[4]"},{"why":"States the Bohigas–Giannoni–Schmit conjecture that chaotic quantum systems locally obey GOE statistics, the premise the paper builds on.","marker":"[5]"},{"why":"Gives the GOE statistics, Wigner-Dyson energy correlations, and random-matrix properties used in the averages.","marker":"[6]"},{"why":"Establishes the non-crossing contraction rule for Gaussian-distributed orthogonal matrix elements, which the paper uses to derive the correlated-moment formula.","marker":"[7]"},{"why":"Contains the Sinai-billiard data and the conjecture $\\Delta = \\hbar \\lambda_{\\max}$ in Eq. (3), which the paper generalizes to all chaotic systems.","marker":"[8]"},{"why":"Presents a matrix model with exponential large-time OTOC decay via Ruelle-Pollicott resonances, which the paper discusses as an apparent discrepancy to its Gaussian prediction.","marker":"[9]"}],"fun_headline_variants":["OTOC Gaussian decay reveals classical Lyapunov time","Quantum chaos OTOC maps directly to Lyapunov exponent","Late-time OTOC decay controlled by Lyapunov time","Gaussian OTOC decay links quantum to classical chaos","OTOC decay rate set by classical Lyapunov exponent"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation rests on the assumption that the energy correlation function $F(E_m - \\bar{E}_\\alpha)$ is a Gaussian of width $\\Delta$; the paper itself notes this form rests on numerical evidence and would be a Lorentzian for weak chaos, which would turn the Gaussian OTOC decay into exponential decay.","fun_headline_variants_meta":{"raw":{"variants":["OTOC Gaussian decay reveals classical Lyapunov time","Quantum chaos OTOC maps directly to Lyapunov exponent","Late-time OTOC decay controlled by Lyapunov time","Gaussian OTOC decay links quantum to classical chaos","OTOC decay rate set by classical Lyapunov exponent"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000153,"raw_usage":{"total_tokens":1187,"prompt_tokens":902,"completion_tokens":285,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":206}},"tokens_in":518,"tokens_out":285,"duration_ms":3549,"temperature":1.0,"reasoning_tokens":206,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:18:09.208016+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the two-point energy correlation function of a chaotic quantum system (for instance the Sinai billiard or a kicked rotor) in the semiclassical regime and compare it to a Gaussian; if the correlation function is measurably non-Gaussian, or if the measured large-time OTOC decay is exponential rather than Gaussian with the predicted $\\hbar/(2\\Delta)$ width, the central claim fails. A cheaper test: check whether the conjectured relation $\\Delta = \\hbar \\lambda_{\\max}$ holds by computing $\\Delta$ from spectral statistics and $\\lambda_{\\max}$ from the classical dynamics of the same system.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the universal parametrization of the chaotic Hamiltonian and the Gaussian form of the energy correlation function $F$, which is the input for the whole calculation."},{"cited_title":"Maldacena, S","cited_arxiv_id":null,"evidence_quote":"Provides the definitions of $C(t)$ and $F(t)$, the $\\beta/4$ insertion trick, and the expected asymptotic behavior that the paper computes explicitly."},{"cited_title":"Bohigas, M","cited_arxiv_id":null,"evidence_quote":"States the Bohigas–Giannoni–Schmit conjecture that chaotic quantum systems locally obey GOE statistics, the premise the paper builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the GOE statistics, Wigner-Dyson energy correlations, and random-matrix properties used in the averages."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the non-crossing contraction rule for Gaussian-distributed orthogonal matrix elements, which the paper uses to derive the correlated-moment formula."},{"cited_title":"Arve, Phys","cited_arxiv_id":null,"evidence_quote":"Contains the Sinai-billiard data and the conjecture $\\Delta = \\hbar \\lambda_{\\max}$ in Eq. (3), which the paper generalizes to all chaotic systems."},{"cited_title":"Garcia-Mata, M","cited_arxiv_id":null,"evidence_quote":"Presents a matrix model with exponential large-time OTOC decay via Ruelle-Pollicott resonances, which the paper discusses as an apparent discrepancy to its Gaussian prediction."}],"review_version":1}