{"id":"32f69654-05d9-4818-919e-fb9bb4cc12d6","arxiv_id":"1909.00266","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A phase-modulated transverse pump in a cavity-BEC system is predicted to produce a rigid, tunable incommensurate time crystal, with the cavity photon number pulsing at a period set by the detuning from a parametric resonance.","lead":"This paper predicts that shaking the optical lattice created by a pump beam in a cavity-BEC system produces a long-lived time crystal whose oscillations repeat at a non-integer multiple of the drive period. The signature is a slow pulsing of the photons leaking from the cavity, measurable with existing atom-cavity setups.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermodynamic-limit persistence rests on an unvalidated finite-N TWA extrapolation; no lifetime scaling or exact benchmark is provided.","rationale":"The reader's conditionality is well placed. The paper's strongest claim is a long-lived, thermodynamically persistent incommensurate time crystal, and the evidence for that persistence is the finite-N TWA decay in Fig. 12 together with the assertion that the mean-field limit is exact as Na→∞ at fixed NaΔ0. The mean-field exactness argument is standard for the cavity-mediated all-to-all coupling, but the finite-N lifetime extrapolation is not quantified: no scaling exponent, no fit, no error bar, and no comparison with an exact small-system calculation. The TWA adds leading-order quantum fluctuations to the initial state and cavity noise, but both the dephasing of the BDW1 order and the cavity noise contribution are controlled only by the approximate phase-space method. These are exactly the places where a driven-dissipative many-body system could heat or dephase faster than predicted, which would shorten the lifetime at any finite Na and weaken the N→∞ inference. My reading identifies the same weakest assumption as the reader's, so I do not change the verdict. The proposed lifetime-scaling test is a direct, relatively inexpensive check: if the decay time does not grow without bound with Na, or if the TWA envelope decay disagrees with a small exact master-equation benchmark, the central persistence claim would need to be lowered; if the lifetime does diverge and the benchmark agrees, the current CONDITIONAL verdict can be upgraded. I considered whether the conceptual ambiguity of an incommensurate period (Eq. 16 vs. Eq. 17) should be the primary concern, but it is secondary to the numerical persistence bridge, because even a quasiperiodic response would still be interesting if it persists in the thermodynamic limit. The paper gives clear credit where due: the parameter-free combination in Eq. 17 and the robustness checks against drive perturbations are meaningful evidence, and the authors explicitly acknowledge finite-N decay and strong-interaction metastability. Those acknowledgements were weighed and do not resolve the missing scaling analysis, which is why the concern remains load-bearing rather than fatal.","tokens_in":17281,"tokens_out":16198,"duration_ms":181755,"concrete_test":"Run TWA simulations for Na = 6×10^4, 1.2×10^5, 2.4×10^5, 4.8×10^5, and 9.6×10^5, rescaling Δ0 so that NaΔ0 is fixed and using the same filtered |α|^2 signal as in Fig. 12; extract the envelope decay time τ(Na) and fit its scaling. The persistence claim is supported only if τ(Na) diverges (at least linearly) in this range; as a control, compare one TWA decay curve against an exact Lindblad master-equation calculation for a small truncated few-mode system to validate the TWA dephasing rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the subharmonic response to persist in the thermodynamic limit, but the only quantitative bridge from finite-N simulations to N→∞ is Fig. 12, where a few TWA decay curves are shown and the text asserts that the oscillation lifetime increases with Na. No lifetime scaling law is fitted, no data or code are shipped, and the TWA dephasing rates are not benchmarked against an exact solution of the quantum master equation for a truncated few-mode model. If the true dissipative dynamics produce dephasing beyond TWA, or if τ(Na) grows sublinearly or saturates for larger Na, the infinite-lifetime extrapolation fails, and the mean-field limit cycle of Eq. (12) could be an artifact of the semiclassical approximation rather than a genuine time-crystalline phase. Because Eq. (16) also asserts exact periodicity with TB while Eq. (17) permits an irrational TB/T, the incommensurate claim itself needs a quantitative spectral characterization; the finite-N lifetime analysis is the most load-bearing weakness in establishing persistence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a cavity-BEC system in which a transverse pump is phase-modulated, producing a shaken lattice. It reports a dynamical phase in which the cavity photon number and a bond-density-wave order parameter oscillate with period TB = 2π/[ω_d(1 − ω_res/ω_d)], where ω_res = 2ω_rec√ϵ_p, so that the response can be incommensurate with the drive. The results are obtained from semiclassical mean-field equations of motion and truncated Wigner approximation (TWA) simulations using parameters motivated by the Hamburg experiment. The authors map a dynamical phase diagram and test robustness to stochastic noise, drive imperfections, contact interactions, and finite atom number, concluding that the time crystal persists in the thermodynamic limit.","tokens_in":17425,"tokens_out":4454,"duration_ms":109275,"significance":"If the central claim holds, the paper offers an experimentally accessible platform for an incommensurate time crystal, with a parameter-free prediction for the subharmonic period and a directly measurable cavity output. The strengths of the paper are its explicit use of experimentally motivated parameters, a defined momentum-mode truncation, TWA simulations that go beyond mean field, and several independent robustness checks. The main weakness is that the thermodynamic-limit persistence rests on a small set of finite-size TWA curves without a fitted scaling law, uncertainty estimates, or an exact benchmark, so the conclusion is plausible but not yet fully established.","major_comments":[{"comment":"The statement that the finite-size oscillations decay but 'their lifetime increases with Na, which leads to an infinite decay time in the thermodynamic limit' is an extrapolation without quantitative support. Figure 12 shows a small number of TWA curves with no fitted lifetime τ(Na), no error bars, and no comparison against an exact solution of the quantum master equation for a truncated few-mode model. If τ(Na) grows sublinearly or saturates, the persistence claim fails. Please add a scaling analysis of the lifetime versus Na and an exact small-system benchmark (or a quantitative TWA-validity argument) to make this load-bearing step reproducible.","section":"Sec. IV, Fig. 12"},{"comment":"The incommensurate character of the time crystal is asserted from Eq. (17), which uses ω_BDW1 ≈ ω_res. The paper does not provide a quantitative spectral characterization showing that the numerically observed peak is at the predicted frequency to within the spectral resolution and is distinct from any rational multiple of ω_d over the accessible time window. Figure 5 reports the ratio ω_BDW1/ω_d but gives no linewidths or fit residuals. Please add a quantitative comparison of the extracted frequency with Eq. (17), including uncertainties, to support the 'incommensurate' claim.","section":"Sec. III B, Eq. (17)"},{"comment":"The construction of the dynamical phase diagram is not fully specified. The text refers to long-time averaged Φ_DW1 and |Φ_BDW1| and to the 'leading order parameter and the long-time dynamics,' but the exact criteria for separating BEC, DW1, TC, and DW2 (including how the metastable DW2 state is excluded from the TC region) are not stated. Please define the classification rules, the time windows, and any thresholds used so that the phase boundaries in Fig. 1(d) are reproducible.","section":"Sec. III C, Fig. 1(d)"}],"minor_comments":[{"comment":"The number of momentum modes is inconsistent: Sec. II states that modes spanning {n,m} ∈ [−6,6]ℏk are used, which gives 169 modes, while Sec. IV states M = 81 momentum modes. Please reconcile these numbers.","section":"Sec. II vs. Sec. IV"},{"comment":"The inset captions in Fig. 3 refer to a 'DW4' phase, while the main text and Fig. 3(l) describe the DW2 phase. Please correct this notation.","section":"Fig. 3 caption"},{"comment":"The Gaussian filter used to remove fast oscillations in Figs. 11 and 12, and also in Fig. 3(l), is not specified. Please give the filter width and comment on how the filtering affects the reported oscillation lifetimes.","section":"Sec. IV, Figs. 11 and 12"},{"comment":"The TWA calculations are based on 10^3 trajectories, but no statistical uncertainties are reported. Adding error bars to Fig. 12 would strengthen the finite-size comparison.","section":"Sec. IV, Fig. 12"},{"comment":"Equation (17) introduces m without defining it immediately. Please define m = ω_res/ω_d and clarify that the 'noninteger multiple' of the driving period corresponds to TB/T = 1/(1−m).","section":"Eq. (17)"},{"comment":"No data or code availability statement is included; for a numerical study, shipping the simulation code or depositing the data would materially aid verification by readers.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"This is a solid but not fully verified numerical study. The main risk is the finite-N to thermodynamic-limit extrapolation; the authors should be encouraged to provide a lifetime scaling law and an exact benchmark. The paper also leans heavily on earlier parametric-resonance results from the same group (Refs. [56,62,63]); the editor may wish to ensure the novelty relative to those works is clearly delineated in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a theory proposal for an incommensurate time crystal in a shaken cavity-BEC system. The headline result is that a phase-modulated transverse pump excites a bond-ordered density wave that switches between two checkerboard states at a period TB = 2π/(ωd − ωres), where ωres = 2ωrec√εp. The period is not fitted; it is a parameter-free combination of known frequencies, so the circularity burden is low. The claim is specific: the cavity photon number should oscillate at 2/TB and this should be visible in the emitted light of an existing apparatus. That is a genuinely new prediction—previous work on Dicke time crystals gave a commensurate period, and the shaken-cavity work of Ref. [51] did not identify TTSB.\n\nWhat the paper does well: it uses experimentally motivated parameters, gives a clear physical picture of the BDW1 order, maps a phase diagram, and tests robustness with truncated Wigner simulations against drive imperfections and contact interactions. The comparison with the commensurate amplitude-modulated case is useful. The numerics look internally consistent, and the mode truncation is plausible.\n\nThe soft spots are real but not fatal. The biggest one is the thermodynamic-limit extrapolation. The paper shows finite-N TWA decay curves in Fig. 12, states that the lifetime grows with N, and concludes it becomes infinite in the N → ∞ limit. But it fits no lifetime scaling law and offers no benchmark of TWA against an exact small-system master equation. This is the load-bearing step for calling it a true time crystal rather than a long-lived metastable state. The authors acknowledge the metastability for strong interactions, but they don't quantify the crossover. A referee should ask for either a scaling law or an exact calculation in a truncated few-mode model.\n\nSecondary issues: no code or data are shipped, so the 81-mode convergence assertion is hard to check; the phase boundaries in Fig. 1(d) come from long-time averages without error bars; and Eq. (16) asserts exact periodicity with TB, which is stronger than what the simulations actually show—a direct spectral test of incommensurability would be cleaner.\n\nWho it's for: theorists and experimentalists working on time crystals in driven-dissipative atom-cavity systems. The prediction is falsifiable and specific. I'd send it to peer review, and I'd ask for the lifetime analysis and data/code before accepting. The central idea is sound enough to warrant the effort.","headline":"Specific, falsifiable prediction of an incommensurate time crystal in a shaken cavity-BEC system, with the main weakness being the unquantified finite-N to thermodynamic-limit extrapolation.","tokens_in":18036,"tokens_out":3004,"would_cite":true,"duration_ms":27691,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A shaken atom-cavity system produces a long-lived incommensurate time crystal, with cavity photons oscillating at a tunable non-integer multiple of the drive period.","keywords":["time crystal","incommensurate time crystal","atom-cavity system","bond density wave","parametric resonance","shaken lattice","driven-dissipative dynamics","truncated Wigner approximation"],"falsifier":"An exact master-equation simulation (for example via quantum trajectories) for a small but finite atom number in the same parameter regime would show whether the cavity photon peak at $2\\omega_{\\mathrm{DW1}}$ survives beyond the truncated-Wigner prediction; alternatively, an experiment could measure the period $T_B = 2\\pi/(\\omega_d(1 - \\omega_{\\mathrm{res}}/\\omega_d))$ and check whether the oscillation amplitude decays over thousands of cycles at fixed $N_a$, with a lifetime that scales with $N_a$ as plotted in Fig. 12.","tokens_in":1766,"feed_emoji":"⚛️","tokens_out":2115,"duration_ms":58752,"temperature":0.7,"pith_summary":"This paper claims that phase-modulating the transverse pump of an atom-cavity system, which shakes the optical lattice, creates an incommensurate time crystal: the cavity photon number and atomic order parameters oscillate at a non-integer multiple of the driving period, spontaneously breaking the discrete time-translation symmetry of the Hamiltonian. The subharmonic response corresponds to dynamical switching between bond-ordered density wave states, where atoms self-organize at the nodes of the pump-cavity potential, a state that does not exist in equilibrium. The authors map out a wide time-crystal region in the pump-intensity versus shaking-amplitude phase diagram and show, using semiclassical simulations, that the oscillations are rigid against quantum fluctuations, small drive imperfections, and weak collisional interactions, with lifetimes that grow with atom number. If the claim holds, the cavity's emitted light provides a direct in situ readout of time-translation symmetry breaking in a driven-dissipative many-body system.","feed_headline":"Shaken cavity atoms form a tunable incommensurate time crystal","feed_subtitle":"A phase-modulated pump makes cavity photons oscillate at a non-multiple of the drive period, persisting over thousands of cycles.","key_machinery":"The load-bearing object is the bond-density-wave order parameter $\\Phi_{\\mathrm{BDW1}}$, a nonequilibrium density order in which atoms sit at the nodes ('bonds') of the combined pump-cavity potential. The driving phase $\\phi(t) = f_0\\sin(\\omega_d t)$ couples parametrically to the product of $\\Phi_{\\mathrm{BDW1}}$ and the real part of the cavity field, and when the drive frequency exceeds the superradiant resonance frequency $\\omega_{\\mathrm{res}}$, this coupling excites a large-amplitude oscillation of $\\Phi_{\\mathrm{BDW1}}$ near $\\omega_{\\mathrm{res}}$. The beating between the drive frequency and the BDW1 frequency produces the slow period $T_B$. The simulations solve coupled c-number equations of motion for the cavity mode and up to 81 momentum modes, with cavity noise included in the truncated Wigner approximation; the mean-field limit is argued to become exact as the atom number $N_a \\to \\infty$ at fixed $N_a\\Delta_0$.","core_discovery":"The paper establishes that resonant excitation of the dynamical bond-density-wave order parameter $\\Phi_{\\mathrm{BDW1}} = \\langle\\sin(ky)\\cos(kz)\\rangle$ at frequency $\\omega_{\\mathrm{BDW1}} \\approx \\omega_{\\mathrm{res}} = 2\\omega_{\\mathrm{rec}}\\sqrt{\\epsilon_p}$ produces a time-translation-symmetry-breaking response of the cavity mode with period $T_B = 2\\pi/\\omega_{\\mathrm{DW1}} = 2\\pi/(\\omega_d(1 - \\omega_{\\mathrm{res}}/\\omega_d))$. During each such period the system switches between two checkerboard-like BDW1 configurations, and the cavity photon occupation $|\\alpha|^2$ oscillates at $2\\omega_{\\mathrm{DW1}}$, a non-integer multiple of the drive period when $\\omega_{\\mathrm{res}}/\\omega_d$ is not rational. The authors demonstrate through mean-field and truncated Wigner simulations that this incommensurate time crystal persists for thousands of drive cycles, is robust to cavity-loss noise, to drive-amplitude perturbations up to about 20%, and to weak contact interactions, and that its finite-size decay time increases with atom number, extrapolating to infinite lifetime in the thermodynamic limit.","pith_inferences":["The beating mechanism should generalize beyond cavity QED: any driven-dissipative system with a parametric resonance whose natural frequency is slightly detuned from the drive can exhibit an incommensurate subharmonic period, making this a generic route to time quasicrystals.","Because the time crystal is read out through the emitted cavity field, one could test beyond-TWA effects by measuring photon correlations or heterodyne spectra and comparing them with the TWA predictions; disagreement would signal quantum corrections not captured by the semiclassical treatment.","Tuning $\\omega_{\\mathrm{res}}/\\omega_d$ from a rational value (commensurate Dicke time crystal) to an irrational value should produce a continuous transition from period doubling to a genuinely incommensurate oscillation, which may clarify how discrete time-translation breaking evolves into quasiperiodic behavior."],"forward_implications":["The incommensurate time crystal can be observed in situ simply by monitoring the cavity photons, since $|\\alpha|^2$ oscillates at $2\\omega_{\\mathrm{DW1}}$, a frequency much slower than the drive, so a photodetector suffices.","The subharmonic frequency is tunable by adjusting either the pump intensity, which sets $\\omega_{\\mathrm{res}}$, or the drive frequency, allowing both incommensurate and commensurate responses (for example $T_B = 80T$) in the same platform.","The phase is stable against drive-amplitude imperfections up to $\\delta/f_0 \\approx 20\\%$ and against weak contact interactions up to $E_{\\mathrm{int}}/E_{\\mathrm{rec}} = 0.05$, indicating that existing atom-cavity experiments can realize it without perfect control.","For finite atom numbers the oscillations decay, but their lifetime grows with $N_a$, supporting a genuine time crystal in the thermodynamic limit rather than a transient prethermal effect.","The BDW1 state is a new nonequilibrium density order that does not exist in equilibrium and is stabilized only by the periodic driving, expanding the set of attainable ordered phases in cavity-BEC systems."],"supporting_citations":[{"why":"Supplies the experimental parameters (Na = 60x10^3, omega_rec, kappa, Delta_0) used in all simulations.","marker":"[49]"},{"why":"Provides the density-wave order parameter and the mean-field treatment underlying the superradiant phase.","marker":"[47]"},{"why":"Establishes dissipative discrete time crystals in cavity QED and the Na to infinity mean-field exactness used in the thermodynamic-limit argument.","marker":"[25]"},{"why":"Introduces the parametric resonance and the dynamical normal phase that the BDW1 resonance generalizes.","marker":"[62]"},{"why":"Supplies the amplitude-modulation scheme and DW2 excitation against which the shaken case is compared.","marker":"[56]"},{"why":"Provides the truncated Wigner approximation used for beyond-mean-field rigidity and persistence tests.","marker":"[64]"},{"why":"Gives the Dicke time crystal as the commensurate comparison and motivates the expectation of stability against many-body correlations.","marker":"[28]"}],"fun_headline_variants":["Incommensurate time crystal emerges in shaken atom-cavity","Cavity atoms create incommensurate time crystal via shaking","Shaken cavity atoms form persistent incommensurate time crystal","Incommensurate time crystal from shaken cavity persists","Cavity atoms tick at noninteger times in shaken lattice"],"cache_read_input_tokens":20224,"weakest_assumption_plain":"The predicted long-lived time crystal rests on the semiclassical equations of motion, with cavity noise added only in the truncated Wigner approximation, being faithful to the full quantum master equation; if genuine quantum fluctuations or heating not captured by this treatment scramble the BDW1 oscillations, the subharmonic response could decay instead of persisting.","fun_headline_variants_meta":{"raw":{"variants":["Incommensurate time crystal emerges in shaken atom-cavity","Cavity atoms create incommensurate time crystal via shaking","Shaken cavity atoms form persistent incommensurate time crystal","Incommensurate time crystal from shaken cavity persists","Cavity atoms tick at noninteger times in shaken lattice"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001479,"raw_usage":{"total_tokens":5951,"prompt_tokens":960,"completion_tokens":4991,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":4905}},"tokens_in":576,"tokens_out":4991,"duration_ms":32731,"temperature":1.0,"reasoning_tokens":4905,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:56:57.105160+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An exact master-equation simulation (for example via quantum trajectories) for a small but finite atom number in the same parameter regime would show whether the cavity photon peak at $2\\omega_{\\mathrm{DW1}}$ survives beyond the truncated-Wigner prediction; alternatively, an experiment could measure the period $T_B = 2\\pi/(\\omega_d(1 - \\omega_{\\mathrm{res}}/\\omega_d))$ and check whether the oscillation amplitude decays over thousands of cycles at fixed $N_a$, with a lifetime that scales with $N_a$ as plotted in Fig. 12.","supporting_citations":[{"cited_title":"Dynamical phase transition in the open Dicke model,","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental parameters (Na = 60x10^3, omega_rec, kappa, Delta_0) used in all simulations."},{"cited_title":"Self-organization of a Bose-Einstein condensate in an optical cavity,","cited_arxiv_id":null,"evidence_quote":"Provides the density-wave order parameter and the mean-field treatment underlying the superradiant phase."},{"cited_title":"Discrete Time- Crystalline Order in Cavity and Circuit QED Systems,","cited_arxiv_id":null,"evidence_quote":"Establishes dissipative discrete time crystals in cavity QED and the Na to infinity mean-field exactness used in the thermodynamic-limit argument."},{"cited_title":"Dynamical many-body phases of the parametrically driven, dissipative Dicke model,","cited_arxiv_id":null,"evidence_quote":"Introduces the parametric resonance and the dynamical normal phase that the BDW1 resonance generalizes."},{"cited_title":"Dy- namical Control of Order in a Cavity-BEC System,","cited_arxiv_id":null,"evidence_quote":"Supplies the amplitude-modulation scheme and DW2 excitation against which the shaken case is compared."},{"cited_title":"Phase space representation of quantum dy- namics,","cited_arxiv_id":null,"evidence_quote":"Provides the truncated Wigner approximation used for beyond-mean-field rigidity and persistence tests."},{"cited_title":"Dicke time crystals in driven-dissipative quantum many-body systems,","cited_arxiv_id":null,"evidence_quote":"Gives the Dicke time crystal as the commensurate comparison and motivates the expectation of stability against many-body correlations."}],"review_version":1}