{"id":"d0d4b467-a682-4847-bfd8-a5efd6ccfead","arxiv_id":"2509.06028","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A quasiclassical Green function path integral reproduces and extends known spectral form factor results for quantum thermalization, including non-perturbative and diffusive regimes.","lead":"This paper develops a path-integral method for describing quantum thermalization, where a system acts as its own environment, and tests it on two model systems. It connects short-time scattering to the Heisenberg-time plateau in the spectral form factor, a key observable of quantum chaos.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8) fails the decoupling limit: at Γ=0 it gives (1+K_1)^L instead of the exact K_1^L, because the nonperturbative replacement was applied to the large-t form Eq. (6) rather than the exact Potts result Eq. (50).","rationale":"The reader correctly flags the non-perturbative extension as the least secure part of the paper, specifically the unproven equal damping of Altshuler-Andreev saddles and the symmetry-inferred late-time action. My stress test found a more concrete and more easily falsified problem in the very same extension: Eq. (8) does not reduce to the exact noninteracting SFF when the coupling is turned off. The root cause is that the replacement t→K_1(t) was applied to the large-t approximate formula Eq. (6), not to the exact t-state Potts expression Eq. (50) that Eq. (6) approximates. Consequently Eq. (8) misses the '-1' that is essential for the Γ→0 limit. This is not a matter of an unproven expectation; it is an internal algebraic inconsistency in the central result. The error is numerically significant at the small D values used for validation (D=8,12), where t∼D is not large enough to justify the t≈t-1 replacement. The numerical agreement in Figs. 2–3 may therefore be achieved partly by fitting Γ away from its golden-rule value. The concrete test is decisive and inexpensive: set Λ=0 and compare Eq. (8) to the exact factorized limit; the formula fails. Because the paper is otherwise careful and the issue is likely repairable by using the exact Potts form with K_1-1, the reader's CONDITIONAL verdict remains appropriate rather than outright rejection. I disagree with the reader only in the precise location: the weakest assumption is not just the AA-saddle damping expectation, but the algebraic form of the interpolation used to build Eq. (8).","tokens_in":35850,"tokens_out":17509,"duration_ms":181206,"concrete_test":"Analytical check: evaluate Eq. (8) at Λ=0 (Γ=0), t<D. It yields (1+K_1(t))^L=(1+t)^L; the exact factorized result is K_1(t)^L=t^L. These differ for every L>1 and t>0. Equivalently, run the brickwork simulation with Λ=0 and confirm K(t)=K_1(t)^L; any agreement with Eq. (8) would indicate the formula is not the SFF of the stated model. If the authors intend Eq. (8) only for t≫L, they should state the correction terms and re-fit Fig. 2 with the exact Potts form (1+(K_1-1)e^{-Γt})^L + ... to show the fitted Γ then matches Λ^2/4.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central non-perturbative formula Eq. (8) is internally inconsistent in a limit it must satisfy. The exact semiclassical Potts result, Eq. (50), is K=(1+(t-1)e^{-Γt})^L + (t-1)(1-e^{-Γt})^L, which reduces to K=t^L at Γ=0 and to K=1 at t=1. Eq. (6) is only the t≫1 simplification of Eq. (50), replacing both t-1 factors by t. The non-perturbative extension in §V.D replaces every t by K_1(t) (or K_L(t)) in Eq. (6), producing Eq. (8) with (1+K_1 e^{-Γt})^L. Setting Γ=0 in Eq. (8) gives (1+K_1(t))^L, whereas the model with Λ=0 is L decoupled Haar qudits and must give K_1(t)^L. The discrepancy is O(L/K_1) and is not a D^{-1} correction; for the D=8,12 numerics used to validate Eq. (8) it is tens of percent. The same issue affects the second term (factor t vs t-1). Thus Eq. (8) is at best an asymptotic large-t formula, not the claimed non-perturbative generalization, and the numerical fits of Γ may be absorbing the error. This is distinct from, but compounds, the acknowledged expectation about AA-saddle damping in §V.D.1.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a real-time path-integral formalism for many-body quantum thermalization, built from quasiclassical Green functions in a tensor-product Hilbert space. The central objects are collective G-variables, whose fluctuations are controlled by a nonlinear sigma model; the GΣ/Luttinger-Ward construction is used to couple qudits. For a brickwork circuit model the authors derive the spectral form factor, first semiclassically as a t-state Potts model (Eqs. (49)–(50)) and then propose a nonperturbative generalization Eq. (8) designed to incorporate single-qudit and many-body Heisenberg-time cutoffs. For an array of capacitively coupled quantum dots they derive a diffusive late-time form factor Eq. (9). The results are compared with exact diagonalization data for L=2–5 and D=8–64, with one fit parameter (Γ) in the discrete case.","tokens_in":36325,"tokens_out":6864,"duration_ms":74151,"significance":"Should the central claims hold, this is a substantial step: an effective field theory for the crossover from short-time semiclassical dynamics to the many-body Heisenberg time, with a modular construction that could apply to other models. The paper is largely self-contained, recovers the known Potts-model result, and provides a clear semiclassical-exactness argument for the cancellation of higher-order D^{-1} corrections (Appendix D). Reproducible numerical data are deposited on Zenodo. However, the headline nonperturbative result Eq. (8) fails a basic consistency test in the decoupling limit, and the Heisenberg-time extension relies on unproved expectations about AA saddles and the universal ergodic sigma model. Because these are load-bearing, the paper in its present form does not fully establish the announced beyond-perturbation-theory framework.","major_comments":[{"comment":"Equation (8) fails the decoupling limit Γ→0. The exact semiclassical Potts result, Eq. (50), is K=(1+(t-1)e^{-Γt})^L+(t-1)(1-e^{-Γt})^L and reduces to K=t^L at Γ=0. Equation (8) is obtained by replacing t→K_1(t) and t→K_L(t) in Eq. (6), the large-t simplification of Eq. (50) in which every (t-1) is replaced by t. At Γ=0, Eq. (8) gives (1+K_1(t))^L, but the model with Λ=0 consists of L decoupled Haar qudits and must give K_1(t)^L. The discrepancy is O(L/K_1), not O(D^{-1}); for the D=8,12 systems in Figs. 2–3 it is tens of percent at t∼D where Γt≪1. Thus Eq. (8) is not the claimed non-perturbative generalization; it is at best an asymptotic large-t formula, and the fitted Γ values may be absorbing this systematic error. The extension should be applied to Eq. (50) (e.g. by replacing (t-1) with K_1(t)-1 and K_L(t)-1), or the precise domain of validity of Eq. (8) must be stated.","section":"§V.D, Eq. (8)"},{"comment":"The non-perturbative extension rests on two unproved inputs: (i) the AA-saddle contributions are damped at the same rate Γ as the standard fluctuations ('we thus expect', §V.D.1), and (ii) the late-time effective theory is the universal ergodic sigma model with D^L levels, Eq. (54), inferred from symmetry rather than derived from the microscopic action. Both are load-bearing for Eq. (8) and for the claim of describing dynamics up to the many-body Heisenberg time. The first is not a routine technicality: AA saddles are stationary points distinct from the standard saddle, and their coupling to interactions need not coincide with the Golden-rule damping of near-standard fluctuations. Since the paper defers the full technical execution to a forthcoming publication, the current manuscript does not yet provide a first-principles derivation of its central nonperturbative result. The authors sho","section":"§V.D.1–D.2"}],"minor_comments":[{"comment":"The sentence 'we obtain a find the form factor' appears to have a missing word; it should likely read 'we find the form factor'.","section":"§II.B"},{"comment":"The sentence fragment 'Here, the unit operators on the diagonal' is incomplete; the definition of G^{-1}(B) should be finished.","section":"Appendix B.1"},{"comment":"The displayed action contains an unreadable block of symbols; please check the typesetting of the fermion bilinear term.","section":"Eq. (21)"},{"comment":"The role of the two fit parameters (C and Γ) and the softening of the max-function over λ-scales should be stated explicitly in the main text; currently the lack of numerical control is mentioned only in the caption.","section":"§VII.D and Eq. (9)"},{"comment":"Reference [20] lists only an arXiv identifier; for a journal submission, please provide the published version if available.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The paper is ambitious but relies heavily on the authors' own prior work (Refs. [22,25]) and defers key nonperturbative derivations to a forthcoming publication. The decoupling-limit failure of Eq. (8) is a concrete correctness issue that should be resolved before the paper can be accepted; it may be fixable in revision by re-deriving Eq. (8) from Eq. (50). The numerical comparison at large D supports the semiclassical part of the paper, but the beyond-perturbation-theory claim is not yet fully proven."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nTwo things to know before reading this paper. First, it is an ambitious, largely self-contained attempt to build a quasiclassical Green function path integral for quantum thermalization, covering brickwork circuits and Hamiltonian arrays. Second, the flagship nonperturbative result, Eq. (8), is internally inconsistent: in the decoupling limit Γ→0 it does not reduce to the exact form factor of L independent qudits.\n\nWhat is genuinely new and good: the formalism combines the color-flavor transform, GΣ functional, and nonlinear sigma model into a coherent framework. The continuous-time treatment of capacitively coupled dots, yielding diffusive relaxation through time-translation modes, is new, and the symmetry-enriched result Eq. (58) extends earlier work on conserved charges. The paper is careful to state where it is relying on expectations, and the semiclassical derivation of Eq. (6) is detailed and matches Ref. [5].\n\nThe soft spots, in proportion: the nonperturbative extension in Section V.D is the weak underbelly. The stress-test note is correct: the exact Potts sum, Eq. (50), contains factors (t−1), not t. Replacing t by K_1 in the large-t form Eq. (6) to obtain Eq. (8) gives (1+K_1 e^{−Γt})^L, which at Γ=0 is (1+K_1)^L rather than K_1^L. That is not a small D−1 correction; for D=8, L=2 it is an error of tens of percent. The authors' fitted Γ may be absorbing this error. The correct extension would need to replace t−1 by K_1−1, at least in the first term, and the physical justification for that replacement is still missing. The late-time effective theory in Eq. (54) is likewise inferred from symmetry, not derived. These gaps are acknowledged, but the Γ=0 inconsistency is a concrete mathematical contradiction that should have been caught.\n\nWho this is for: people working on many-body quantum chaos and non-equilibrium field theory will want to see the framework, and the paper deserves a serious referee because the formalism is promising and the issue is likely fixable. But the central claim of a first-principles description up to the Heisenberg time is not currently supported.\n\nRecommendation: send it to peer review, but ask the referee to check the decoupling limit and the t→K_1 replacement. I would be skeptical until that is resolved.","headline":"The paper builds a serious quasiclassical framework, but Eq. (8) fails the decoupling limit because the nonperturbative replacement was applied to the large-t form instead of the exact Potts result.","tokens_in":36734,"tokens_out":3548,"would_cite":false,"duration_ms":37270,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper introduces a quasiclassical Green function path integral that describes unitary yet irreversible thermalization dynamics, and derives a closed spectral form factor for chaotic brickwork circuits that extends the semiclassical resu","keywords":["quantum thermalization","spectral form factor","quasiclassical Green function","GΣ (Luttinger-Ward) formalism","Altshuler-Andreev saddle","color-flavor transform","many-body Heisenberg time","brickwork quantum circuit"],"falsifier":"Run a brickwork-circuit simulation at L=3, D=64 with a coupling set so Γ≈1/D, resolving the spectral form factor around t=D and t=D^3. Equation (8) predicts a kink of height D at t=D and a plateau onset at t=D^3; if the Altshuler-Andreev contributions decay with a rate different from Γ, the kink height or the plateau onset will shift visibly without any fit parameter.","tokens_in":35753,"feed_emoji":"🔁","tokens_out":9300,"duration_ms":97744,"temperature":0.7,"pith_summary":"The paper introduces a real-time path integral over slow 'quasiclassical Green functions' as a first-principles description of quantum thermalization in systems that act as their own environment. Its central result is a closed formula for the spectral form factor of a brickwork circuit, K(t) = (1 + K_1(t)e^{-Γt})^L + K_L(t)(1-e^{-Γt})^L, which interpolates from the short-time t^L ramp of independent qudits, through an exponential synchronization at rate Γ, to the plateau of a single ergodic system at the many-body Heisenberg time t ~ D^L. The same construction, adapted to energy-conserving quantum-dot arrays, yields a soft, diffusive approach to ergodicity on times ~L^2. If correct, this provides a transferable effective field theory for strongly entangled many-body chaos, with parameter-free agreement against numerics when the local Hilbert-space dimension is large.","feed_headline":"One formula now tracks thermalization to the Heisenberg time","feed_subtitle":"A quasiclassical Green function field theory derives the spectral form factor through the full thermalization crossover.","key_machinery":"The central object is the quasiclassical Green function G^ab_tu, a collective variable formed by summing products of retarded and advanced wave amplitudes over a qudit's D internal states. Three pieces of machinery carry the argument: (1) a GΣ (Luttinger–Ward) construction that elevates G to an integration variable, turning the low-energy sector into a nonlinear sigma-model of Goldstone modes associated with a spontaneously broken causal symmetry; (2) the color–flavor transform, which converts Haar-averaged unitary evolution into a matrix B-field living on the symmetric space U(2t)/U(t)×U(t); and (3) Altshuler–Andreev saddle points, non-perturbative stationary points whose inclusion replaces","core_discovery":"The central discovery is that the unitary dynamics of an ergodic many-body system can be represented by a real-time path integral whose slow variables are quasiclassical Green functions G^ab_tu = -(ia/D)∑_μ ψ^a_{μt}ψ^b*_{μu}. Interactions between subsystems act as an environment: they damp all fluctuation modes except one globally synchronized 'Hadamard' mode X_mn = ∏_j B_{j,mn}, which carries the late-time ergodic phase. Inclusion of Altshuler-Andreev saddle points—the non-perturbative configurations that terminate the ramp at time D—yields the closed form factor formula Eq. (8), and a symmetry-restoring effective theory with D^L levels extends it to the many-body Heisenberg time. The paper","pith_inferences":["A natural next test is to push the brickwork formula to L=3 or 4 with D significantly larger than 12 (e.g., via tensor-network methods), where the kink at t=D and the plateau at t=D^L become sharp enough to distinguish Eq. (8) from the semiclassical Eq. (6).","If the uniform damping of Altshuler-Andreev saddles holds generally, the same construction should extend to spatially non-uniform couplings, producing site-dependent rates Γ_j and a generalized synchronization mode that interpolates between local and global ergodicity.","The equality of the L-qudit propagator with a single D^L-level ergodic propagator suggests an attractor: any chaotic arrangement of pairwise interactions should flow to the same late-time universal spectral form factor, making the framework a practical probe for device thermalization without solving the full Schrödinger equation."],"forward_implications":["For the brickwork model, the spectral form factor crosses over from t^L to t exponentially in Γt, and the non-perturbative replacements t→K_1 and t→K_L keep the formula valid at and beyond the single-qudit Heisenberg time t~D.","In symmetry-enriched brickwork circuits, the late-time form factor factorizes into the singlet form factor times an SU(N)-Heisenberg evolution, so conserved-charge diffusion sets a minimum thermalization (Thouless) time ~L^2.","For continuous Hamiltonian systems, energy diffusion replaces exponential relaxation, giving Eq. (9) with its characteristic max(...) factors and the same L^2 scaling of the slowest mode.","Because the construction requires only fast-thermalizing qudits with local pairwise correlations, it is modular and transferable to other system classes whose randomness breaks integrability.","The same path integral provides a diagnostic for quantum thermalization in engineered devices, distinguishing the ergodic, diffusive, and plateau regimes from a single analytical expression rather than from exact diagonalization."],"supporting_citations":[{"why":"supplies the semiclassical t-state Potts result (Eq. 6) that this paper extends beyond perturbation theory.","marker":"[5]"},{"why":"provides the conserved-charge brickwork analysis and the diffusive L^2 time scale that the symmetry section generalizes.","marker":"[6]"},{"why":"supplies the coherent-state Schwinger-Keldysh path integral from which the GΣ action is built.","marker":"[8]"},{"why":"is the Luttinger-Ward functional whose GΣ structure allows interactions to be traded for quasiclassical Green functions.","marker":"[12]"},{"why":"gives the ergodic single-system form factor K_n(t) used to replace semiclassical t factors at Heisenberg time.","marker":"[15]"},{"why":"introduces the Altshuler-Andreev saddle points whose inclusion produces the non-perturbative contributions to Eq. (8).","marker":"[16]"},{"why":"provides the color-flavor transform that converts Haar-averaged unitary dynamics into the B-field representation.","marker":"[21]"},{"why":"gives the quantitative energy-domain treatment of the form factor with Altshuler-Andreev saddles, extended here to the time domain.","marker":"[22]"}],"fun_headline_variants":["Path integral tracks thermalization to Heisenberg time","One framework now maps full thermalization crossover","Quasiclassical Green functions reach Heisenberg time","From scattering to ergodicity: a unified path integral","Spectral form factor derived across all time scales"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that interaction damps the rare Altshuler-Andreev saddle fluctuations at the same golden-rule rate Γ as ordinary fluctuations, a step the paper states as an expectation rather than a derivation; the late-time symmetry action is likewise inferred from symmetry, not derived from the microscopic action.","fun_headline_variants_meta":{"raw":{"variants":["Path integral tracks thermalization to Heisenberg time","One framework now maps full thermalization crossover","Quasiclassical Green functions reach Heisenberg time","From scattering to ergodicity: a unified path integral","Spectral form factor derived across all time scales"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000358,"raw_usage":{"total_tokens":1781,"prompt_tokens":751,"completion_tokens":1030,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":957}},"tokens_in":495,"tokens_out":1030,"duration_ms":11171,"temperature":1.0,"reasoning_tokens":957,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T04:36:15.654515+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a brickwork-circuit simulation at L=3, D=64 with a coupling set so Γ≈1/D, resolving the spectral form factor around t=D and t=D^3. Equation (8) predicts a kink of height D at t=D and a plateau onset at t=D^3; if the Altshuler-Andreev contributions decay with a rate different from Γ, the kink height or the plateau onset will shift visibly without any fit parameter.","supporting_citations":[{"cited_title":"(41) reads as 0=δ Σj S[ ¯G, ¯Σ, B]=D ( ¯Gj+iτ 3(G−1(Bj)−i ¯Σjτ3)−1) , 0=δ Gj S[ ¯G, ¯Σ, B]=D ¯Σj+δ Gj Sint[ ¯G]","cited_arxiv_id":null,"evidence_quote":"supplies the semiclassical t-state Potts result (Eq. 6) that this paper extends beyond perturbation theory."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the conserved-charge brickwork analysis and the diffusive L^2 time scale that the symmetry section generalizes."},{"cited_title":"However, it reappears in opera- tors representing correlation functions in the form of gauge-phase factors containing discrete time integrals overΘ n","cited_arxiv_id":null,"evidence_quote":"supplies the coherent-state Schwinger-Keldysh path integral from which the GΣ action is built."},{"cited_title":"It implies the vector integral generalization δµρ = ∫ Dψ e− ¯ψψ ψµ ¯ψρ (A1) whereψ={ψ µ}now is aD-component Grassmann vector, Dψ= ∏µ d ¯ψµdψµ and ¯ψψ= ∑µ ¯ψµψµ","cited_arxiv_id":null,"evidence_quote":"is the Luttinger-Ward functional whose GΣ structure allows interactions to be traded for quasiclassical Green functions."},{"cited_title":"Our starting point is the configurations Eq","cited_arxiv_id":null,"evidence_quote":"gives the ergodic single-system form factor K_n(t) used to replace semiclassical t factors at Heisenberg time."},{"cited_title":"(3) as de- scribed by an integral overG=−iT τ 3T−1 with the action Eq","cited_arxiv_id":null,"evidence_quote":"introduces the Altshuler-Andreev saddle points whose inclusion produces the non-perturbative contributions to Eq. (8)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the color-flavor transform that converts Haar-averaged unitary dynamics into the B-field representation."},{"cited_title":"Fran¸ ca, Making quantum dynamics exact, Physics14, 60 (2021), publisher: American Physical Society","cited_arxiv_id":null,"evidence_quote":"gives the quantitative energy-domain treatment of the form factor with Altshuler-Andreev saddles, extended here to the time domain."}],"review_version":1}