{"id":"8e0e8c0f-7e80-4ae3-9617-972a48ee28e7","arxiv_id":"2512.02755","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Kinetic driving eliminates nearest-neighbor hopping in a Bose-Hubbard lattice and produces an effective four-boson Hamiltonian whose spectral statistics and OTOC dynamics match the bosonic Sachdev-Ye-Kitaev model.","lead":"This paper shows that periodically modulating the hopping amplitude of a Bose-Hubbard lattice—'kinetic driving'—can produce an effective all-to-all four-particle interaction that mimics the Sachdev-Ye-Kitaev (SYK) model. The authors match the spectral form factor and out-of-time-ordered correlator (OTOC) signatures of this driven model, known as the KDBH model, to the bosonic SYK model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Untested Magnus truncation: all SYK diagnostics rely on the lowest-order effective Hamiltonian, never on the full driven lattice","rationale":"The reader's weakest assumption correctly identifies the gap between validating H_KDBH and validating the actual driven system. My independent reading confirms that the paper's derivation of H_KDBH is internally consistent: Eq. (EM9) follows from the Fourier transform and Jacobi-Anger expansion, and the SFF collapse for H_KDBH is a nontrivial indicator of chaos. However, the step from 'H_KDBH has SYK-like SFF/OTOC' to 'the driven system realizes SYK physics' requires the Floquet-Magnus expansion to be controlled. The paper's own End Matter acknowledges the neglect of higher-order terms, so this is an explicit, testable limitation rather than a hidden flaw. The proposed concrete test—a direct Floquet simulation for small systems—would resolve whether the neglected terms are consequential. If the test passes, the practicality claim is much stronger; if it fails, the effective-model results are misleading. Therefore the reader's CONDITIONAL verdict remains appropriate, and I recommend no change.","tokens_in":13460,"tokens_out":9935,"duration_ms":98426,"concrete_test":"Diagonalize the full one-period Floquet operator U(T) for the driven Bose-Hubbard Hamiltonian H(t) = -J0 cos(ωt) Σ_j (b†_{j+1} b_j + h.c.) + U Σ_j b†_j b†_j b_j b_j, for N=6 bosons on L=6 sites with periodic boundary conditions, at κ=J0/ω=0.8 and ω/U = 5, 10, 20. Restrict to the zero-momentum, reflection-symmetric sector used in Fig. 2. Compute the spectral form factor from the quasienergy spectrum and the OTOC of Eq. (5) at stroboscopic times t=nT, and compare with the corresponding H_KDBH predictions. If the scaled SFF agrees in its drop-ramp-plateau structure (ramp time within a factor of 2) and the OTOC mean deviation is below 0.1 for all three frequencies, the concern is resolved. If deviations exceed this, the effective-Hamiltonian validation is insufficient and the practical claim is unsupported for those parameters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that a kinetically driven Hubbard model 'indeed reproduces SYK physics' (abstract). All numerical evidence—SFF in Fig. 2 and OTOC in Fig. 3—is computed from the static effective Hamiltonian H_KDBH, Eqs. (EM8)-(EM9), obtained as the leading (time-averaged) term of a Magnus expansion in the interaction picture. The authors explicitly acknowledge in the End Matter after Eq. (EM9): 'higher-order terms in the effective Hamiltonian would correspond to many-body interactions... These are neglected in the present study, which assumes a high frequency regime ω≫U.' No estimate or numerical test of the magnitude of these corrections is provided. For an experimental realization, the drive frequency ω must be large enough that the neglected terms (including six-operator interactions) are small compared to the U-scale couplings, but the paper never checks whether e.g. ω=10U and κ=0.8 preserves the SFF or OTOC. If the next-order terms are not negligible, the actual shaken lattice is not described by H_KDBH, and the claimed 'practical and accurate platform' fails. This is load-bearing because every SYK diagnostic is computed on the approximate model, not on the actual time-dependent system; the paper explicitly flags this as a limitation, making it an acknowledged gap rather than a hidden assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a Floquet-engineering scheme, termed kinetic driving, in which the hopping amplitude of a Hubbard model is modulated as J(t)=J0 cos ωt so that the time-averaged single-particle hopping vanishes. In the high-frequency limit, a Magnus expansion yields an effective static Hamiltonian, the KDBH model, consisting of all-to-all quartic interactions with deterministic amplitudes Q_{ijkl} (Eqs. EM8-EM9). The authors argue that this structure resembles the bosonic SYK model, and they support this by exact-diagonalization studies of the spectral form factor and infinite-temperature OTOCs, comparing against a Gaussian random bosonic SYK model. The central claim is that a shaken optical lattice provides a practical and accurate platform for quantum simulation of SYK physics. All numerical diagnostics are computed on the lowest-order effective Hamiltonian H_KDBH, not on the full time-dependent Floquet system.","tokens_in":13767,"tokens_out":8295,"duration_ms":90624,"significance":"If the leading-order effective Hamiltonian faithfully represents the driven lattice, the proposal is significant: it offers a relatively simple cold-atom route toward a bosonic SYK-like model, with analytic control over the generated couplings. The derivation of H_KDBH is transparent and the SFF collapse in Fig. 2 is a genuinely suggestive signature of nonlocal chaotic dynamics. The paper also provides a useful contrast with the undriven Bose-Hubbard model, showing that its SFF lacks the universal ramp. The use of exact diagonalization on finite systems is appropriate, and the comparison with an independently defined Gaussian SYK model for OTOCs is a reasonable benchmark. However, the evidence is currently incomplete in one load-bearing respect: the validation never tests whether the truncated Magnus expansion describes the actual time-dependent driven system, despite the paper explicitly acknowledging that higher-order terms are neglected.","major_comments":[{"comment":"The central validation is performed exclusively on the lowest-order Magnus Hamiltonian H_KDBH. The End Matter states: “higher-order terms in the effective Hamiltonian would correspond to many-body interactions … These are neglected in the present study, which assumes a high frequency regime ω≫U.” No estimate or numerical test of these corrections is provided. Since the claim is that the driven system “reproduces SYK physics” and that the lattice setup is a “practical and accurate platform,” the truncation is load-bearing. The next-order Magnus correction in the rotating frame is roughly of order U^2/ω (with Bessel-function factors depending on κ); for plausible parameters such as ω=10U and κ=0.8 this is not obviously negligible. I ask the authors to either simulate the full time-dependent Hamiltonian H(t) and compute the same SFF and OTOC diagnostics, or to estimate the norm/effect of th","section":"§Model and End Matter after Eq. (EM9)"},{"comment":"The claimed quantitative agreement between the KDBH and SYK OTOCs relies on a single fitting parameter: the time coordinate of the SYK data is rescaled by a factor of 9. Because a global time rescale can always align two decay curves of similar shape, this weakens the claim that the SYK model and KDBH model are quantitatively equivalent. The factor 9 should be derived from the model parameters, e.g., from the root-mean-square of the Q_{ijkl} amplitudes relative to the Gaussian J_{ij;kl} variance, or at least shown to be consistent with such an estimate. Without this, the OTOC comparison demonstrates a common functional form but not a quantitative match.","section":"Fig. 3(c) and Results (OTOC comparison)"},{"comment":"The SFF is compared to the universal RMT drop-ramp-plateau form, not directly to the SFF of the bosonic SYK Hamiltonian of Eq. (3). Since the universal SFF form is shared by many chaotic nonlocal models, this provides only indirect evidence for SYK specifically. A direct overlay of the KDBH SFF and the SYK SFF for matched Hilbert-space dimension and the same symmetry sector would make the claimed equivalence much more compelling. This is not a fatal flaw, but it is a gap between the abstract's wording (“direct comparison of the spectral form factor”) and what is actually plotted.","section":"Fig. 2 and Table I"},{"comment":"The paper describes H_KDBH as having “quasi-random all-to-all interactions,” but the amplitudes Q_{ijkl} are deterministic, translationally invariant, and non-Gaussian with gaps (Fig. 1(b)). The justification relies on sparse-SYK results from Refs. [63,64], which are formulated for random sparsification of SYK couplings. It is not self-evident that those criteria transfer to a deterministic structured coupling matrix. The diagnostics shown are reassuring, but the paper should either state more precisely what is meant by “quasi-random” or provide a more direct test, e.g., comparing the KDBH diagnostics against a random model with the same sparsity and coupling distribution. This would clarify the resemblance to SYK beyond the two specific observables shown.","section":"Fig. 1(b) and §Model (sparsity and determinism)"}],"minor_comments":[{"comment":"The closed-form expression for f(λ) is quoted from Ref. [58]. A one-line commutator derivation would make the End Matter self-contained and avoid requiring the reader to consult the prior work.","section":"End Matter around Eq. (EM5)"},{"comment":"The SFF is averaged over four values of κ, with the statement that the form does not depend on κ. Showing individual curves or a measure of spread would substantiate this claim and help the reader judge the robustness of the collapse.","section":"Fig. 2 caption and text"},{"comment":"The SYK OTOC is averaged over 25 random samples, but no sample-to-sample fluctuations or error bars are shown. Since the comparison is central to the paper, a statement of the statistical uncertainty would be useful.","section":"Fig. 3(c) caption and text"},{"comment":"The main text uses the BH model at U=J for OTOCs, while the End Matter SFF for the BH model uses U=0.2J and the text cites GOE behavior for U/J<0.1. The parameter regime relevant to Table I should be specified consistently, so that the reader knows which BH model is being compared.","section":"Results and End Matter (SFF of BH model)"},{"comment":"The two-frequency implementation is described only qualitatively, with “a proper separation of energy scales” assumed. A rough estimate of the required ω/U, the achievable κ range, and the expected heating rate would strengthen the claim that the proposal is practical.","section":"Concluding remarks (experimental realization)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's central idea is appealing and the derivation is clean, but the gap between the effective Hamiltonian and the actual driven system is acknowledged yet not addressed. I would recommend major revision rather than rejection because the missing test (full Floquet simulation or a quantitative estimate of higher-order corrections) is within the scope of the manuscript. The authors should also consider whether the title and abstract overstate the directness of the comparison, given the fitted time rescale and the indirect SFF comparison. The paper draws heavily on the authors' previous work (Refs. [58,65]); the novelty relative to those papers should be made explicit in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nYou should know two things before reading this paper. First, it does something new: it takes the kinetically-driven Bose-Hubbard effective Hamiltonian (KDBH) and shows numerically that its spectral form factor and OTOC decay match the bosonic SYK model. That specific comparison is not in the earlier kinetic-driving or spectral-statistics papers. Second, the whole SYK claim is tested only on the lowest-order Magnus effective Hamiltonian, never on the actual time-dependent shaken lattice. The authors admit this in the End Matter but do not quantify the error.\n\nWhat it does well: the derivation of H_KDBH is coherent and self-contained, and the exact diagonalization numerics are clean. The contrast with the undriven Bose-Hubbard model is useful: BH has GOE statistics in some parameter range but does not show the universal drop-ramp-plateau SFF, so they correctly argue that spectral statistics alone is not enough. The paper is honest about the non-Gaussian, sparse structure of the Q amplitudes and cites the sparse-SYK literature. That is good scholarship.\n\nThe soft spots are real but not fatal to the core theoretical observation. The main one: every diagnostic—SFF in Fig. 2, OTOC in Fig. 3—comes from H_KDBH, the time-averaged first term of a Magnus expansion. The condition ω≫U is stated, but there is no estimate or numerical test of the neglected higher-order many-body terms. For the claimed 'practical and accurate platform,' this is load-bearing. It is an acknowledged gap, not a hidden assumption, but it is still a gap. A simple check on a small system with the full time-dependent Hamiltonian (e.g., ω=10U, κ=0.8) would settle a large part of it.\n\nThe OTOC agreement also rests on a fitted time-rescale factor of 9. That is a free parameter, and the factor is large enough to make the reader wonder whether the agreement is a forced shape match. It would help if they showed the rescale is consistent with the bandwidth difference between the two Hamiltonians, rather than just chosen for the fit. The phrase 'quasi-random all-to-all interactions' overstates the deterministic, spatially structured Q. And the SFF is compared to RMT universality rather than directly to the SYK SFF; that is a minor quibble because SYK itself follows RMT, but a direct overlay would strengthen the claim.\n\nWho is it for? People working on Floquet engineering of many-body chaos, cold-atom quantum simulation, and sparse SYK variants. It deserves a serious referee: the idea is plausible, the methods are standard, and the missing checks are addressable. My recommendation is to send it to peer review, with the request that the authors either simulate the full driven system at finite ω or estimate the next-order Magnus terms, and that they justify the OTOC time rescale rather than present it as a single fitted number. Conditional on those, it is a solid contribution.","headline":"Useful Floquet-SYK proposal with an honest but untested Magnus truncation; worth refereeing, but the practical claim needs finite-frequency checks.","tokens_in":14250,"tokens_out":2896,"would_cite":true,"duration_ms":29652,"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":"A periodically driven Bose-Hubbard chain can reproduce the defining signatures of the Sachdev-Ye-Kitaev model, including its universal spectral form factor and non-local scrambling, offering a practical optical-lattice route to SYK physics.","keywords":["Sachdev-Ye-Kitaev model","Floquet engineering","Bose-Hubbard model","spectral form factor","out-of-time-ordered correlator","kinetic driving","quantum chaos","optical lattices"],"falsifier":"Simulate the full periodically driven Bose-Hubbard model at finite driving frequency (e.g., ω ~ 3–10 U) and compute the stroboscopic SFF and OTOC; if these deviate from the effective KDBH predictions or require ω so large that heating dominates, the claimed practical platform would fail. An experiment could measure OTOCs in a shaken optical lattice and compare directly with the effective-model prediction.","tokens_in":13356,"feed_emoji":"⚛️","tokens_out":3675,"duration_ms":39393,"temperature":0.7,"pith_summary":"This paper claims that a simple periodic modulation of hopping in the Bose-Hubbard model—kinetic driving—produces an effective Hamiltonian with the same quartic all-to-all interaction structure as the bosonic Sachdev-Ye-Kitaev (SYK) model. The drive cancels single-particle hopping and generates long-range, quasi-random couplings. Exact diagonalization shows that the spectral form factor of this effective model collapses onto the universal drop-ramp-plateau curve, and the out-of-time-ordered correlators match the SYK result after a single time rescaling. The authors conclude that shaking an optical lattice is a practical way to simulate SYK physics in cold atoms.","feed_headline":"Shaking a Hubbard lattice yields SYK quantum chaos","feed_subtitle":"Driven Bose-Hubbard mimics SYK spectrum and scrambling, offering a cold-atom route to the model.","key_machinery":"The central object is the effective Hamiltonian H_KDBH (Eq. EM8) with interaction amplitudes Q_ijkl given by Eq. EM9, obtained from the lowest-order Magnus expansion of the driven Bose-Hubbard model. These amplitudes are Bessel-function-weighted momentum sums that eliminate single-particle hopping and generate quartic, number-conserving, long-range interactions. The comparison machinery is the spectral form factor (SFF) and the normalized out-of-time-order correlator (OTOC), computed by exact diagonalization in a symmetry-restricted sector.","core_discovery":"The paper establishes that the kinetically driven Bose-Hubbard model, with hopping J(t)=J0 cos(ωt), has a high-frequency effective Hamiltonian H_KDBH consisting solely of quartic terms U∑Q_ijkl b†_i b†_j b_k b_l, where the Q amplitudes are determined by Bessel functions and decay with distance. Although these amplitudes are not random and the Hamiltonian is about 50% sparse, the spectral form factor shows the universal SYK drop-ramp-plateau structure with ramp time scaling as 1/D, and the averaged OTOC decays in the same form as the bosonic SYK model. The authors verify this by exact diagonalization for several system sizes and driving parameters, concluding that kinetic driving of Hubbard-t","pith_inferences":["The paper validates only the static effective Hamiltonian, not the full time-dependent Floquet dynamics; a natural next step is to check whether finite-frequency corrections preserve the SFF and OTOC signatures.","Since the Q amplitudes fall off with distance, the model is effectively long-range but not fully random; one might expect a size or energy-scale crossover beyond which SYK signatures degrade—this crossover is not analyzed.","The authors show that breaking translational symmetry improves the Gaussian character of the amplitude distribution, suggesting that engineered disorder could make the effective model even closer to the SYK ensemble.","If heating can be controlled, the same kinetic-driving mechanism could be adapted to simulate other all-to-all random-interaction models beyond SYK, such as generalized sparse random Hamiltonians or higher-dimensional driven lattices."],"forward_implications":["If correct, the Bose-Hubbard model with periodically modulated hopping provides a feasible experimental path to SYK physics in existing cold-atom optical lattices.","The effective model exhibits the universal drop-ramp-plateau SFF, indicating spectral rigidity and quantum chaos at all timescales.","OTOCs decay immediately regardless of operator separation, confirming the absence of spatial locality and the presence of fast scrambling.","The method extends to fermionic Hubbard models with nearest-neighbor interactions, giving the same Q amplitude distribution and thus the potential to simulate fermionic SYK models.","Sparse and non-Gaussian interaction distributions (about 50% non-zero amplitudes) do not break the SYK universality, supporting the use of sparse SYK models in experiments."],"fun_headline_variants":["Kinetic driving makes Bose-Hubbard a SYK simulator","Hubbard under drive: SYK scrambling from cold atoms","Floquet Hubbard: a practical route to SYK physics","Shaking hopping in Hubbard yields SYK chaos"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire validation is performed on the static effective Hamiltonian from the lowest-order Magnus expansion, relying on ω≫U to make higher-order many-body terms negligible, but the paper does not test this by simulating the actual time-dependent driven system.","fun_headline_variants_meta":{"raw":{"variants":["Kinetic driving makes Bose-Hubbard a SYK simulator","Hubbard under drive: SYK scrambling from cold atoms","Floquet Hubbard: a practical route to SYK physics","Shaking hopping in Hubbard yields SYK chaos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000583,"raw_usage":{"total_tokens":2559,"prompt_tokens":702,"completion_tokens":1857,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":446,"completion_tokens_details":{"reasoning_tokens":1789}},"tokens_in":446,"tokens_out":1857,"duration_ms":16610,"temperature":1.0,"reasoning_tokens":1789,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T18:54:35.137825+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the full periodically driven Bose-Hubbard model at finite driving frequency (e.g., ω ~ 3–10 U) and compute the stroboscopic SFF and OTOC; if these deviate from the effective KDBH predictions or require ω so large that heating dominates, the claimed practical platform would fail. An experiment could measure OTOCs in a shaken optical lattice and compare directly with the effective-model prediction.","supporting_citations":[],"review_version":1}