Pith. sign in

REVIEW 3 major objections 5 minor 5 references

The Yukawa–SYK model, with its boson mass tuned relative to the random coupling, continuously interpolates between SYK2-like single-particle chaos and SYK4-like many-body chaos, and a cavity-QED setup can realize it.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 23:34 UTC pith:ARCVNLXN

load-bearing objection A genuinely useful finite-size study of Yukawa-SYK, but the SYK2-like endpoint leans on the hard-core boson truncation and needs a cutoff-scaling check before the bridge claim is taken literally. the 3 major comments →

arxiv 2511.04762 v3 pith:ARCVNLXN submitted 2025-11-06 quant-ph cond-mat.dis-nncond-mat.quant-gascond-mat.str-elhep-th

From single-particle to many-body chaos in Yukawa--SYK: theory and a cavity-QED proposal

David Pascual Solis , Alex Windey , Soumik Bandyopadhyay , Andrea Legramandi , Philipp Hauke This is my paper
classification quant-ph cond-mat.dis-nncond-mat.quant-gascond-mat.str-elhep-th
keywords Yukawa–SYKquantum chaosSachdev–Ye–Kitaevspectral form factorout-of-time-ordered correlatorsprethermalizationcavity QEDlevel statistics
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper tries to establish that a single physical knob—the ratio of boson mass ω0 to coupling strength g, in the form ω0/g^{2/3}—turns the Yukawa–SYK model from a quadratically integrable system with single-particle chaos into a quartic, many-body chaotic system of the SYK4 type. The authors support this with exact-diagonalization data for small systems (N=8 fermions, M=4 bosons), using gap-ratio statistics, the spectral form factor, and out-of-time-ordered correlators, and introduce a time-rescaling that lets YSYK observables be quantitatively compared with SYK2 and SYK4. In between the two limits they find a crossover region with prethermalization plateaus, two-step scrambling, and partial ergodicity breaking. They also propose an ultracold-atoms-in-an-optical-cavity implementation with estimated interaction rates in the megahertz (Yukawa) and kilohertz (effective SYK4) range, with dissipation suppressed enough for coherent observation. If correct, this makes YSYK a tunable experimental platform for studying the transition from single-particle to many-body chaos.

Core claim

The paper's central claim is that the dimensionless ratio ω0/g^{2/3} acts as a continuous dial between two distinct chaos regimes: for ω0/g^{2/3}≪1 the Yukawa–SYK model behaves like complex SYK2—Poissonian many-body level statistics, a superlinear spectral-form-factor ramp, an OTOC that saturates to a non-zero prethermal plateau and then undergoes a delayed secondary decay on a time scale ~g/ω0^{5/2}—while for ω0/g^{2/3}≫1 it behaves like complex SYK4, with GUE level statistics, a linear ramp, an OTOC that decays to zero, and an effective quartic interaction obtained by Schrieffer–Wolff elimination of the bosons. The authors argue this bridge is quantitative after rescaling time by the varia

What carries the argument

The central object is the Yukawa–SYK Hamiltonian, H = -μΣ c†c + ... + (1/√(MN)) Σ gi j,k c†_i c_j (a_k + a†_k), with GUE-distributed random couplings gi j,k and boson mass ω0. The control parameter is the ratio ω0/g^{2/3}. The argument is carried by two perturbative tools: for small ω0, a Magnus expansion in the boson mass treats the SYK2-like quadratic limit and predicts the second scrambling time ~g/ω0^{5/2}; for large ω0, a Schrieffer–Wolff transformation integrates out the bosons and yields an effective four-fermion SYK4-like coupling J~g^2/ω0^2. The diagnostics—gap-ratio distribution, spectral form factor, and OTOCs—are made quantitatively comparable to the SYK2/SYK4 benchmarks through

Load-bearing premise

The quantitative claims rest on exact diagonalization of very small systems (N=8 fermions, M=4 bosons, hard-core boson cutoff N_b=1), and if the Poisson-to-GUE crossover, the plateau-time jump, or the OTOC secondary decay shift or vanish as N and M grow, the claimed interpolation may not survive in the thermodynamic limit.

What would settle it

A calculation or experiment that checked the plateau-time jump and the OTOC collapse at, say, N=12–16 fermions and M=6–8 bosons (keeping the boson cutoff low enough for diagonalization) and compared the crossover position and the secondary-decay time t~g/ω0^{5/2} would settle whether the ω0/g^{2/3} scaling is real or a small-system artifact; alternatively, a cavity-QED experiment measuring the two-step scrambling at a single tunable detuning would test the same prediction directly.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the claim holds, YSYK provides a single Hamiltonian in which one experimental knob (the cavity-laser detuning, which sets ω0) tunes the system from single-particle chaos to many-body chaos, allowing controlled study of the crossover.
  • The paper's time-rescaling framework gives a concrete way to compare any near-integrable or near-SYK4 system to the SYK benchmarks, so the same method can be applied to other boson-mediated models.
  • The predicted two-step scrambling with an intermediate plateau and a late secondary decay at time ~g/ω0^{5/2} is a testable dynamical signature; the collapse of OTOC curves predicted by this scaling can be checked directly.
  • The Schrieffer–Wolff analysis shows the large-ω0 limit is an effective SYK4 with coupling g^2/ω0^2, meaning the model reproduces SYK4 physics without requiring the explicit four-fermion term.
  • The cavity-QED proposal, with a single-atom cooperativity around 10, predicts interaction rates exceeding loss rates, making the YSYK chaos crossover accessible to current experiments.

Where Pith is reading between the lines

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

  • The paper leaves implicit that if the ω0/g^{2/3} scaling survives in the thermodynamic limit, it would identify a universal crossover variable for quantum chaos in boson-mediated models, with implications for models of strange metals and electron-phonon systems.
  • One could test the claimed interpolation by measuring the slope-dip region of the SFF as a function of M at fixed ω0/g^{2/3}: the rank of the effective four-fermion coupling grows with M, so the ramp should sharpen toward SYK4 as M increases.
  • The two-step scrambling mechanism—plateau followed by secondary decay—suggests that weakly broken conservation laws, not the details of the Yukawa coupling, control late-time scrambling; if so, similar prethermal plateaus should appear in any model with an approximate conserved sector.
  • Extending the OTOC analysis to finite temperature would test whether the SYK2-to-SYK4 interpolation also holds at low temperatures, where the SYK4 phase is expected to be maximally chaotic; the paper only treats infinite temperature.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper studies a spinless Yukawa-SYK (YSYK) model in which N complex fermions interact with M bosonic modes through random all-to-all couplings. It introduces a time-rescaling procedure, based on matching the short-time variance of the spectral form factor and OTOC, to compare YSYK observables with those of complex SYK2 and SYK4 benchmark models. Exact-diagonalization results for small systems (mainly N=8, M=4, boson cutoff N_b=1) are used to argue that tuning ω0/g^{2/3} interpolates continuously between a SYK2-like single-particle-chaos regime at small ω0 and a SYK4-like many-body-chaos regime at large ω0. In the strong-coupling (small-ω0) regime the paper reports Poissonian gap statistics, an SFF plateau at the SYK2 time 2N, and a two-step OTOC decay with a prethermal plateau; in the weak-coupling (large-ω0) regime it finds a low-energy-band description governed by a Schrieffer–Wolff effective four-fermion Hamiltonian with SYK4-like spectral rigidity and scrambling. The paper also gives a cavity-QED implementation with concrete parameter estimates and a treatment of dissipation.

Significance. If the central interpolation claim holds, the YSYK model would be a valuable tunable platform for studying the crossover from single-particle to many-body quantum chaos, and the proposed cavity-QED realization would be a timely experimental route. The paper has clear strengths: the perturbative Schrieffer–Wolff and Magnus analysis is presented explicitly, the rescaling protocol is transparent, and several finite-size checks are included. However, the SYK2-like endpoint is established only for hard-core bosons (N_b=1), and the physical model—as well as the cavity-photon implementation—has an unbounded bosonic Hilbert space. Because the time-rescaling factor itself depends on N_b, the quantitative match to SYK2 in the small-ω0 regime is not yet supported for the continuum-boson model. This is a load-bearing issue that tempers the significance of the reported bridge.

major comments (3)
  1. [Sec. 4.2, App. B.1, App. C.1, Fig. 10] The SYK2-like regime is analyzed exclusively with the hard-core boson cutoff N_b=1. The rescaling factor in Eq. (B.7), α_SFF = g sqrt(N_b/(2ω0)), depends explicitly on N_b, and Eq. (B.6) shows that σ_H^2 grows linearly with N_b. For the model defined in Eqs. (1)–(2), and for the cavity-QED realization of Sec. 5 where photons are unbounded, the physically relevant limit is N_b→∞; in that limit the variance entering the short-time matching diverges and the rescaling procedure is ill-defined. The Magnus expansion of App. C.1 is also performed explicitly for hard-core bosons. The only N_b>1 test is the OTOC curve in Fig. 10b for N=6, M=3, N_b=3; no spectral statistics, SFF plateau time, or secondary-scrambling collapse is shown for N_b=2,3. The claim that YSYK 'quantitatively reproduces' SYK2-like behavior at small ω0 therefore needs to be either restricted to the truncated model or supporte
  2. [Sec. 3.3, Eqs. (14) and (16)] The time rescaling is defined by imposing that the quadratic short-time expansions of the SFF and OTOC match the target model. Thus the early-time 'quantitative agreement' in the slope/dip region is imposed by construction, not discovered. The paper should state this more prominently and identify which later-time features are genuinely emergent, e.g., the plateau at 2N, the secondary-decay collapse at t∼g/ω0^{5/2} shown in Fig. 5, and the SYK4 ramp–plateau in Fig. 7. Without this separation, the headline 'quantitative reproduction' overstates the evidence.
  3. [Sec. 4.2 and App. A.2] The claimed sharp crossover in the plateau time (Fig. 4, top-middle inset) is obtained with extraction tolerances δ_pl=0.1 and α Δt_pl=1.9 that are chosen ad hoc. No sensitivity analysis of the extracted t_plateau or t_ramp to these tolerances is given, and no error bars are shown despite disorder averaging. Because the 'jump' of the plateau time from 2N to t_H is a central qualitative result, the authors should show that this feature is robust to the choice of detection thresholds and to sampling noise.
minor comments (5)
  1. [App. B.1, Eq. (B.5)] The expression for the bosonic trace average is confusing: it reads as (N_b+1)^M N_b/(N_b+1)^M = N_b. For a standard harmonic-oscillator truncation one would expect an average of (2n+1), so please clarify the hard-core-boson convention and its relation to the spin representation used in App. C.1.
  2. [Abstract and Sec. 4] The abstract says 'the interaction strength acts as a tunable control parameter,' but in the body the control parameter is the ratio ω0/g^{2/3}. Please rephrase for consistency.
  3. [Sec. 3.1.3] Typo: 'chaothic' should be 'chaotic' in the sentence beginning 'As the system becomes less chaothic'.
  4. [Fig. 2 caption] Minor grammar: 'across coupling range' should be 'across the coupling range'.
  5. [Sec. 5.1, Eq. (30)] The identification in Eq. (30) is ambiguous: it compares g_{ij,k}/√(2ω0) to the experimental scale, but the YSYK coupling enters as g_{ij,k}/√(2ω0 M N) in Eq. (2). Please clarify the normalization used in the mapping.

Circularity Check

0 steps flagged

No significant circularity: the rescaling is explicit calibration and the central interpolation claims rest on independent spectral and dynamical evidence.

full rationale

The paper's quantitative comparison between YSYK and SYK2/SYK4 uses time rescalings defined in Sec. 3.3 (Eqs. 14 and 16) by matching short-time expansions of the SFF and OTOC. This is an openly declared calibration procedure, not a hidden prediction: the paper states 'as an operational procedure we impose that the time evolution of two models coincide at early times.' The subsequently claimed agreements at later times—plateau time 2N, linear ramp emergence, prethermal plateau, secondary scrambling with the derived g/omega_0^{5/2} timescale, and the Schrieffer-Wolff-derived SYK4 effective coupling g^2/omega_0^2—are not fixed by those rescaling factors and are verified by independent exact-diagonalization data and perturbation theory. The alpha factors are computed from Hamiltonian moments or commutators (Appendices B.1, B.2), not fitted to the target curves. The Poisson-to-GUE crossover in the gap ratio, the plateau-time jump to the Heisenberg time, and the approach to SYK4 with increasing omega_0 are all nontrivial numerical observations. Self-citations (e.g., Refs. [75,76]) are used for the cavity-QED implementation context, not to justify the chaos interpolation. The acknowledged finite-size and boson-cutoff limitations (App. D.1, Fig. 10b) are correctness risks rather than circularity. I therefore find no circular step that reduces the paper's central claim to its own inputs.

Axiom & Free-Parameter Ledger

3 free parameters · 6 axioms · 0 invented entities

The ledger shows that the paper introduces no new particles or forces; the bosons and cavity photons are physical or already present in the YSYK model. The main loaded assumptions are computational truncation (N_b=1), the by-construction early-time matching through time rescaling, and extrapolation from N≤8 numerics. The free parameters are mostly diagnostic cutoffs and analysis thresholds, not fitted physics constants.

free parameters (3)
  • Boson occupation cutoff N_b = 1
    All main-text numerics truncate each bosonic mode to N_b=1 (hard-core bosons). This makes the Hilbert space finite and shapes the DOS/SFF cluster structure; App. D.2 checks N_b=3 only for OTOCs, with no systematic N_b→∞ extrapolation.
  • SFF plateau/ramp extraction tolerances = δ_pl=0.1, αΔt_pl=1.9, δ_r=5e-3
    App. A.2 defines plateau and ramp onset times using hand-chosen tolerances; the quantitative statements about t_plateau and t_ramp in Fig. 4 depend on these choices.
  • Savitzky-Golay low-pass filter parameters = unspecified
    App. A.3/Fig. 8 uses a numerical low-pass filter to isolate the oscillating part of F(t); the filter parameters are not given, affecting the collapse shown in Fig. 6b.
axioms (6)
  • domain assumption Random all-to-all GUE couplings g_{ij,k} with variance g^2 and disorder-averaged observables represent the YSYK ensemble.
    Eqs. (1)-(2); the entire numerical and analytical treatment relies on this ensemble and on self-averaging at finite N,M.
  • domain assumption Truncating each bosonic mode at N_b=1 (or 3) captures the relevant physics of the infinite-dimensional bosonic Hilbert space for the probes studied.
    Eq. (2) with cutoff; App. D.2 gives limited checks, but no exact or analytic proof for all regimes.
  • ad hoc to paper Short-time variance matching (σ_H/σ_H') makes YSYK dynamics quantitatively comparable to SYK2/SYK4 at later times.
    Sec. 3.3 Eqs. (13)-(16) introduce the rescaling; matching early-time quadratic decay is imposed, and subsequent agreement at later times is interpreted as evidence for the interpolation.
  • standard math Schrieffer-Wolff and Magnus perturbative expansions converge in the large-/small-ω0 limits at finite N,M.
    App. C uses standard degenerate perturbation theory and Magnus expansion; no rigorous error bounds are given for N=8, M=4, N_b=1.
  • domain assumption Adiabatic elimination of the excited atomic state and cavity modes is valid under the stated hierarchy |Δ_da|≫|Ω_d|≫|Ω_m| and |Δ_cd|≫|Ω_dΩ_m/Δ_da|.
    Sec. 5.1, Eqs. (26)-(34); the claimed mapping to YSYK and SYK4 depends on this hierarchy and on speckle disorder producing sufficiently random couplings.
  • domain assumption Finite-size results at N≤8, M≤8 are representative of thermodynamic-limit behavior.
    App. D.1 explicitly states that finite-size scaling to N→∞ is not feasible; the claims of a robust crossover and two-stage scrambling rely on extrapolation from small sizes.

pith-pipeline@v1.3.0-alltime-deepseek · 39477 in / 13754 out tokens · 122717 ms · 2026-08-03T23:34:33.369663+00:00 · methodology

0 comments
read the original abstract

Understanding how quantum systems transition from integrable to fully chaotic behavior remains a central open problem in physics. The Sachdev--Ye--Kitaev (SYK) model provides a paradigmatic framework for studying many-body chaos and holography, yet it captures only the strongly correlated limit, leaving intermediate regimes unexplored. Here, we investigate the Yukawa--SYK (YSYK) model, where bosonic fields mediate random fermionic interactions, and demonstrate that it naturally bridges single-particle and many-body chaos. Using spectral and dynamical chaos markers, we perform a comprehensive finite-size characterization of the YSYK model. We show that the interaction strength acts as a tunable control parameter interpolating between the SYK$_2$ and SYK$_4$ limits, and introduce a framework enabling direct and quantitative comparison with these benchmark models. In the intermediate regimes, we uncover distinct dynamical regimes marked by partial ergodicity breaking, prethermalization plateaus, and incomplete scrambling. Finally, we propose a feasible optical-cavity implementation of the YSYK model using ultra-cold atoms. Our results establish the YSYK model as a unifying platform connecting single-particle and many-body chaos, paving the way for experimental observation of these phenomena.

Figures

Figures reproduced from arXiv: 2511.04762 by Alex Windey, Andrea Legramandi, David Pascual Solis, Philipp Hauke, Soumik Bandyopadhyay.

Figure 1
Figure 1. Figure 1: (a) Schematic illustration of the interaction term in the spinless Yukawa– SYK model. A set of N spinless fermionic modes (blue dots) interact pairwise via M independent bosons (wiggly lines), through random couplings gi j,k . Each bosonic mode is shown in a different color, and the line shading represents the coupling strengths. The couplings are drawn independently from a Gaussian Unitary Ensem￾ble (GUE)… view at source ↗
Figure 2
Figure 2. Figure 2: Spectral diagnostics of the YSYK model with [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Averaged gap ratio 〈r〉 versus the interaction variable ω0/g 2/3 for boson–fermion mode ratios M/N ∈ {2.0, 1.0, 0.5} (obtained for (M,N) = (8, 4),(6, 6),(4, 8), with bosonic cutoff Nb = 1). Data points show averages over 500 disorder realizations for each (M,N). At small ω0/g 2/3 , 〈r〉 is Poissonian (dotted line). It reaches the GUE value (dashed) in a broad crossover region ω0/g 2/3 ∼ 0.3 − 1 (grey band). … view at source ↗
Figure 4
Figure 4. Figure 4: Main panel: Log-log plot of the SFF K(t) versus the rescaled time t g/ p 2ω0 for N = 8 fermionic modes and M = 4 bosonic modes with cutoff Nb = 1. Curves are color–coded by ω0 (color bar at top) in the range ω0 ∈ [0.01, 0.5]. Dashed vertical lines mark the SYK2 plateau time 2N and the Heisenberg time tH. Solid circles indicates the numerically extracted plateau–onset times tplateau(ω0 ), while squares deno… view at source ↗
Figure 5
Figure 5. Figure 5: Disorder averaged OTOCs of the YSYK model for small values of [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Disorder averaged YSYK OTOCs F(t) for N = 8, M = 4 and Nb = 1 in the large ω0 limit averaged over 1000 samples. (a) F(t) for ω0 = 3 and 10 versus the rescaled time in Eq. (18). The full OTOC (solid lines) shows oscillations that decrease with ω0 . These are absent when restricting the OTOCs to the lowest energy peak in the DOS (dashed lines). The difference between full and restricted OTOC decreases as ω0 … view at source ↗
Figure 7
Figure 7. Figure 7: Disorder averaged SFF of the YSYK model in the weak fermion-boson cou [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Early-time disorder-averaged YSYK OTOC for [PITH_FULL_IMAGE:figures/full_fig_p026_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Disorder-averaged mean gap ratio 〈r〉 versus ω0/g 2/3 . Data for Nf =4, 6, 8 (panels) and M =2−8 (colors). As ω0/g 2/3 increases, 〈r〉 crosses from Poisson statistics (dashed) to a GUE-like plateau (dotted); the gray band indicates the chaotic window as in [PITH_FULL_IMAGE:figures/full_fig_p032_9.png] view at source ↗
Figure 10
Figure 10. Figure 10 [PITH_FULL_IMAGE:figures/full_fig_p032_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Spectral form factor K(t) for the normal-ordered effective model. Solid curves from lighter to darker shades are for increasing M, i.e., increasing rank of the antisymmetric couplings Ji j;kl defined in Eq. (D.2). As M increases, the curves closely approach the SYK4 benchmark (dashed). Data for N = 8 at half filling, averaged over 500 realizations, and normalized by D 2 with D = [PITH_FULL_IMAGE:figures/… view at source ↗

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

5 extracted references · 1 linked inside Pith

  1. [1]

    [1]M. V . Berry and M. Tabor,Level clustering in the regular spectrum, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences356, 375 (1977). [2]O. Bohigas, M. J. Giannoni and C. Schmit,Characterization of chaotic quantum spectra and universality of level fluctuation laws, Phys. Rev. Lett.52, 1 (1984), doi:10.1103/PhysRevLett.52...

    arXiv 1977

  2. [2020]

    [127]G. C. Santra, A. Windey , S. Bandyopadhyay , A. Legramandi and P . Hauke,Complex- ity transitions in chaotic quantum systems: Nonstabilizerness, entanglement, and fractal dimension in syk and random matrix models, doi:10.48550/arXiv.2505.09707 (2025), 2505.09707. [128]S. Denisov, T . Laptyeva, W . Tarnowski, D. Chru ´sci´nski and K. ˙Zyczkowski,Uni- ...

    doi:10.48550/arxiv.2505.09707 2025

  3. [2405]

    19260. [70]I. Kukuljan, S. c. v. Grozdanov and T . c. v. Prosen,Weak quantum chaos, Phys. Rev. B 96, 060301 (2017), doi:10.1103/PhysRevB.96.060301. [71]R.-Q. He and Z.-Y. Lu,Characterizing many-body localization by out-of-time-ordered cor- relation, Phys. Rev. B95, 054201 (2017), doi:10.1103/PhysRevB.95.054201. [72]B. Swingle and D. Chowdhury ,Slow scramb...

    arXiv 2017

  4. [2411]

    17802. [77]M. Baghdad, P .-A. Bourdel, S. Schwartz, F . Ferri, J. Reichel and R. Long,Spectral en- gineering of cavity-protected polaritons in an atomic ensemble, Nature Physics19(8), 1104–1109 (2023), doi:10.1038/s41567-023-02035-1. [78]N. Sauerwein, F . Orsi, P . Uhrich, S. Bandyopadhyay , F . Mattiotti, T . Cantat-Moltrecht, G. Pupillo, P . Hauke and J...

    arXiv 2023

  5. [2412]

    14280. [64]A. M. García-García, B. Loureiro, A. Romero-Bermúdez and M. Tezuka,Chaotic- integrable transition in the sachdev-ye-kitaev model, Phys. Rev. Lett.120, 241603 (2018), doi:10.1103/PhysRevLett.120.241603. [65]Y. Huang, F . G. S. L. Brandão and Y.-L. Zhang,Finite-size scaling of out-of- time-ordered correlators at late times, Phys. Rev. Lett.123, 0...

    Pith/arXiv arXiv doi:10.48550/arxiv.2405.19260 2018