{"id":"95d3e9b3-024c-4a27-8d44-1b0890c296b7","arxiv_id":"2506.12414","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A Floquet analysis of the atom-only master equation gives the complex excitation spectrum of a dissipative time crystal, including a gapless mode at the continuous transition and non-crossing branches at the discontinuous transition.","lead":"This paper computes the low-frequency excitation spectrum of atoms forming a discrete dissipative time crystal inside an optical cavity. It predicts how the excitation frequency and damping behave across the phase transition and shows how to measure them with a weak probe laser.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main-text V1 in Eqs. (4)-(6) omits the factor Δ that the SM derivation (Eq. S30) contains; as printed, the central atom-only equations are internally inconsistent and the Floquet spectra are not reproducible.","rationale":"The paper's central claim is that the atom-only mean-field Floquet spectrum quantitatively reproduces the full cavity-atom response. That requires the printed equations (4)-(6) to be the equations actually solved, and those equations to be a valid reduction of the full model. The main text defines V1 without Δ, while the SM derivation (S30) includes Δ; this is not a stylistic difference. Using the printed definition, the dissipative rate is an order of magnitude larger at the stated parameters (Δ = 0.1κ), so the phase boundaries and linewidths in Fig. 2 would shift if the printed formula were used. This is the most concrete, checkable threat to the central claim: it is an internal inconsistency, not a question of regime validity. The truncated-Wigner benchmark is helpful evidence for the bad-cavity reduction, but it cannot resolve the printed-formula ambiguity because the reader cannot know which V1 entered the numerics. The reader's weakest assumption was the bad-cavity timescale separation; I regard that as partially mitigated by the comparison in Fig. 3 and SM Fig. S2, though those comparisons lack error bars. The V1 discrepancy is more immediately load-bearing because it is a definite inconsistency in the defining equations of the theory. A single recomputation of Fig. 2(c,d) with both definitions settles the issue. If the spectra are unchanged, the main text has a typo and the conditional verdict stands; if they change, the main-text presentation is not reproducible and the claim would need to be re-evaluated. I therefore keep the verdict unchanged at CONDITIONAL.","tokens_in":17524,"tokens_out":21361,"duration_ms":280239,"concrete_test":"Recompute the g1/g0 = 0.6 cut of Fig. 2(c,d) using Eq. (7) with (a) V1 = 4 g^2 δc κ/(δc^2 + κ^2)^2 as printed in the main text and (b) V1 = 4 δc Δ κ g^2/(δc^2 + κ^2)^2 from Eq. (S30), keeping all other parameters fixed. Since Δ = 0.1κ, any change in the phase boundary γFl = 0, the cusp position, or the DTC linewidths by more than a few percent identifies the printed equations as not the ones used; no change would confirm a benign typographical omission.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the atom-only mean-field equations (4)-(6) faithfully encode the cavity-mediated dynamics, since the Floquet matrix Eq. (7), the spectra in Fig. 2, and the Lorentzians in Fig. 3 are all built from them. The main text defines V1 = 4 g^2 δc κ/(δc^2 + κ^2)^2, but the SM derivation defines V1 = 4 δc Δ κ g^2/(δc^2 + κ^2)^2 (Eq. S30). The extra factor Δ is not cosmetic: it enters through the difference |c_+|^2 - |c_-|^2, where c± have denominators δc ± Δ - iκ. For the stated parameters Δ = 0.1κ, the two definitions differ by an order of magnitude, so a reader following the main text cannot reproduce the phase boundaries, linewidths, or transition behavior in Fig. 2. This is the load-bearing inconsistency: if the numerics used the SM definition, the main text is wrong as written; if they used the main-text definition, the benchmark in Fig. 3 is not testing the model described.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the collective excitation spectrum of a dissipative Dicke time crystal. Starting from a cavity-atom master equation, the authors adiabatically eliminate the cavity in the bad-cavity limit to obtain a nonlinear, dissipative atom-only mean-field model (Eqs. (4)-(6)). Linearizing this model around the normal-phase fixed point or the 2T-periodic DTC limit cycle leads to the fluctuation matrix Sigma in Eq. (7); its Floquet eigenvalues over two drive periods define a complex excitation frequency lambda_Fl = gamma_Fl - i nu_Fl. The paper maps gamma_Fl and nu_Fl across the continuous and discontinuous transitions (Fig. 2) and, using truncated Wigner simulations of the full cavity-atom dynamics, shows that probe-induced photon-number changes are well described by Lorentzians centered at nu_Fl with width gamma_Fl (Fig. 3).","tokens_in":17713,"tokens_out":12479,"duration_ms":125849,"significance":"If the results hold, the paper provides a computationally light and experimentally feasible way to determine the low-frequency collective excitations of dissipative time crystals, complementing Liouvillian spectral methods. The approach is explicit: the atom-only equations, the fluctuation matrix, and the Floquet extraction are spelled out; the benchmark against truncated Wigner is direct, and the only fitted quantity in Fig. 3(b)-(c) is the Lorentzian amplitude. The predicted non-analytic vanishing of the frequency at the continuous transition and the non-crossing branches at the discontinuous transition are falsifiable signatures. The main obstacle to accepting the manuscript as written is the inconsistency between the printed master equations and the Supplemental derivation of the dissipative coupling V1, which affects every subsequent numerical result.","major_comments":[{"comment":"The main text defines V1 = 4 g^2(t) delta_c kappa/(delta_c^2+kappa^2)^2, whereas the SM derivation, Eq. (S30), gives V1 = 4 delta_c Delta kappa g^2(t)/(delta_c^2+kappa^2)^2. The factor Delta is not optional: it originates from the difference |c_+|^2 - |c_-|^2 in Eqs. (S27)-(S28) and controls the asymmetry that stabilizes the limit cycle. For the parameters used in Figs. 1-3 (Delta = 0.1 kappa), the two definitions differ by an order of magnitude, so a reader who implements the main-text equations cannot reproduce the phase diagram in Fig. 2 or the Lorentzian positions and widths in Fig. 3. The main text's own definition of gamma_0 below Eq. (7) contains Delta, which underscores that the printed V1 is inconsistent with the rest of the paper. Please correct the main-text definition (or explicitly state that the numerics used the SM definition) and confirm that the Floquet spectra and all comparisons are obtained with the SM V1.","section":"Equations (4)-(6) and the definition of V1 in the main text; SM Eq. (S30)"}],"minor_comments":[{"comment":"The extraction of lambda_Fl from the logarithm of the monodromy eigenvalues requires a choice of branch; since the system is 2T-periodic, nu_Fl is defined modulo omega/2. Please specify how the branch is chosen so that the plotted nu_Fl corresponds to the resonance probed in Fig. 3, especially when eigenvalues cross the branch cut.","section":"SM, 'Numerical simulation of the fluctuation matrix'"},{"comment":"The value of the probe strength eta_0 used in the truncated Wigner simulations of Fig. 3 is not stated in the main text; please provide it (or refer explicitly to the SM section) for reproducibility.","section":"Main text, 'Probing Excitations'"},{"comment":"The color maps in Fig. 2(a)-(b) and Fig. 3(a) have no color bars, so the quantitative values of gamma_Fl/gamma_0, nu_Fl/(omega/2), and the intensity difference cannot be read from the figures. Please add color bars or describe the scales in the captions.","section":"Figures 2 and 3"},{"comment":"The matrix in Eq. (7) is typeset in a way that makes the third-row entries run together; please use a proper matrix environment with clear column separation so that each entry is unambiguous.","section":"Equation (7)"},{"comment":"The expansion of c_+/- in Eq. (S20) is garbled in the typeset version; the intermediate algebra is hard to follow. Please rewrite this expansion with explicit steps.","section":"SM, Eq. (S20)"},{"comment":"The symbol Sigma is used both for the fluctuation matrix and for the time-evolution/monodromy operator in Eq. (S47); this notational conflict should be resolved, for example by using Phi for the monodromy operator.","section":"SM, Eq. (S47)"}],"recommendation":"major_revision","confidential_remarks":"The paper is otherwise solid, and the V1 discrepancy in the main text looks like a typographical omission of Delta rather than a fundamental error: the SM derivation, the definition of gamma_0, and the numerics appear to use the factor Delta. If the authors confirm that all results were obtained with the SM expression and correct the main-text equations accordingly, I would be happy to accept after revision. The formal recommendation of major revision is driven by reproducibility: the central equations as printed do not define the model that is actually solved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: this paper gives a usable, concrete way to compute the low-frequency collective excitation spectrum of a dissipative time crystal from an atom-only mean-field model, and it benchmarks the result against truncated Wigner simulations of the full cavity-atom system. That is a real step beyond the existing bifurcation-theory and Liouvillian-gap treatments. Second: there is a typo-level inconsistency in the definition of V1 that must be fixed before a reader can reproduce anything. The main text defines V1 without the factor Δ; the SM derivation, Eq. (S30), includes it. For the stated parameters Δ = 0.1κ, that is a factor of ten difference. The SM is internally consistent, so I read it as a typo in the main text rather than a wrong derivation, but as printed the central equations are not reproducible.\n\nThe Floquet extraction itself is well laid out: linearize around the 2T-periodic limit cycle, integrate the time-ordered exponential, identify the single zero eigenvalue with conserved spin length, and read off the complex pair as frequency and linewidth. The spectral features across the continuous and discontinuous transitions are concrete and novel—the gapless mode at the continuous transition and the non-crossing branches at the discontinuous one. The probe protocol is practical, and the comparison to truncated Wigner is a genuine benchmark, not a fit; only the Lorentzian amplitude is adjusted.\n\nNow the soft spots, in proportion. The Fig. 3 benchmark is a single parameter region, shown without error bars or scatter statistics, and the claim of \"excellent agreement\" rests on visual overlap; the (c) panel matches reasonably but not as tightly as the text suggests. The bad-cavity assumption is stated clearly but not stress-tested at its validity boundary, and the truncated Wigner simulation is itself semiclassical, so it does not independently validate the atom-only reduction against a fully quantum calculation. No code or data artifacts are provided, which is a real omission for a paper whose headline is that this spectrum is computable.\n\nThis paper is for people working on dissipative time crystals and driven-dissipative cavity QED. It deserves a serious referee: the method is a genuine contribution, the derivation is coherent apart from the V1 typo, and the probe protocol is likely useful for experiments. I would send it to review with a request to correct V1, add error bars or scatter characterization to the Wigner data, and state the simulation parameters fully. With those fixes, I would be comfortable with publication.","headline":"Useful Floquet toolkit for DTC excitations, but the main text has a factor-Δ typo in V1 that must be fixed before any of it is reproducible.","tokens_in":18284,"tokens_out":2962,"would_cite":true,"duration_ms":35539,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.50.Pq","05.70.Fh"],"model":"deepseek-v4-flash","headline":"The low-frequency collective excitations of a dissipative time crystal reduce to a single complex Floquet frequency that quantitatively predicts the measured probe response.","keywords":["dissipative time crystal","Floquet theory","collective excitations","dissipative Dicke model","bad-cavity regime","mean-field approximation","truncated Wigner simulation","phase transition"],"falsifier":"A decisive check is to compare the atom-only Floquet spectrum with the probe response of the full dissipative Dicke model when the cavity is not in the bad-cavity regime (for example, with cavity linewidth $\\kappa$ comparable to $\\omega_\\mathrm{res}$, so that $\\omega\\tau_c \\sim 1$); a mismatch between the predicted and measured Lorentzian centers and widths in that regime would show that the adiabatic elimination is the load-bearing assumption.","tokens_in":17276,"feed_emoji":"🕰️","tokens_out":13942,"duration_ms":148523,"temperature":0.7,"pith_summary":"The paper aims to show that the low-frequency collective excitations of a dissipative time crystal — a driven, lossy atom-cavity system whose atoms spontaneously oscillate at half the drive frequency — are described by a single complex Floquet frequency $\\lambda_\\mathrm{Fl} = \\gamma_\\mathrm{Fl} - i\\nu_\\mathrm{Fl}$. In the bad-cavity limit, where the cavity is fast enough to be eliminated, the atoms obey three nonlinear mean-field equations, and linearizing fluctuations around the time-crystalline limit cycle yields this spectrum via Floquet theory. The spectrum behaves differently at the two kinds of transition: at the continuous transition the excitation frequency vanishes with a cusp-like non-analyticity, while at the discontinuous transition the two phases keep non-crossing frequency branches and coexist in a bistable region. When a weak probe drive is added, the full atom-cavity simulations show Lorentzian response lines whose centers and widths match the Floquet predictions, indicating that the spectrum is measurable and the theory quantitative. This gives an efficient tool for characterizing out-of-equilibrium phases without diagonalizing the full quantum master equation.","feed_headline":"One complex number captures a time crystal's collective excitations","feed_subtitle":"Atom-only Floquet theory predicts the probe line shapes in both phases across the transition.","key_machinery":"The machinery is self-consistent Floquet analysis — linear theory for equations with time-periodic coefficients — applied to the fluctuation matrix $\\Sigma(t)$ obtained by linearizing the atom-only mean-field equations around the time-periodic steady state. Because $\\Sigma(t)$ inherits the $2T$-periodicity of the limit cycle, excitations are classified by the Floquet eigenvalues $\\lambda_j$ obtained from the monodromy matrix $\\Phi = \\mathcal{T}\\exp\\!\\int_0^{2T} d\\tau\\, \\Sigma(\\tau)$, with $\\phi_j = e^{2\\lambda_j T}$. Conservation of spin length forces one eigenvalue to zero, and the remaining pair appears as complex conjugates, so the whole spectrum reduces to one complex number $\\lambda_\\mathrm{Fl} = \\gamma_\\mathrm{Fl} - i\\nu_\\mathrm{Fl}$; the coherent and dissipative cavity-mediated interactions enter through the time-periodic coefficients $V_0(t)$ and $V_1(t)$ derived from the bad-cavity elimination. The same linearization, supplemented by a probe term $\\Sigma_\\mathrm{pr}$, gives the atom-only probe response used to interpret the full simulations.","core_discovery":"The paper's central claim is that, in the bad-cavity regime where the cavity adiabatically follows the atoms, the dissipative time crystal and its low-frequency excitations are fully captured by the atom-only mean-field equations (4)–(6). Linearizing fluctuations $\\delta\\vec{v}$ around the $2T$-periodic limit-cycle solution yields the time-periodic fluctuation matrix $\\Sigma(t)$, and the Floquet eigenvalues of the monodromy matrix define a single complex excitation frequency $\\lambda_\\mathrm{Fl} = \\gamma_\\mathrm{Fl} - i\\nu_\\mathrm{Fl}$ for the atomic degrees of freedom. The paper shows that this spectrum behaves differently at the two phase boundaries: a cusp-like vanishing of $\\nu_\\mathrm{Fl}$ with a rapidly dropping $\\gamma_\\mathrm{Fl}$ at the continuous transition, and non-crossing frequency branches with a bistable region at the discontinuous transition. The same Floquet prediction is then compared with truncated Wigner simulations of the full cavity-atom system driven by a weak probe; the simulated probe-induced photon-number changes reproduce Lorentzian lines centered at $\\nu_\\mathrm{Fl}$ with width $\\gamma_\\mathrm{Fl}$, matching the atom-only theory.","pith_inferences":["A natural extension, left open by the paper, is that the same Floquet linearization applied to a continuous time crystal will show a gapless critical mode, because the broken symmetry is continuous rather than discrete.","Because $\\nu_\\mathrm{Fl}$ vanishes at the continuous transition, the probe response near that boundary could act as a sensitive amplifier or sensor; the paper does not explore metrological consequences.","The truncated-Wigner benchmark is itself semiclassical, so genuine quantum fluctuations are not tested here; a fully quantum treatment in the critical region could shift linewidths even when the bad-cavity condition holds.","The non-crossing branches at the discontinuous transition provide a spectral fingerprint of bistability that could be used to identify first-order character in other driven-dissipative systems."],"forward_implications":["At the continuous transition the excitation frequency $\\nu_\\mathrm{Fl}$ vanishes with a cusp-like non-analyticity, giving a direct spectral signature of the transition that can be looked for in cavity transmission.","At the discontinuous transition the normal-phase and time-crystal frequency branches do not cross, and the bistable region is visible in the spectrum, providing a clear way to distinguish first-order from continuous transitions.","A weak probe drive creates Lorentzian features in the transmitted photon number whose center and width equal $\\nu_\\mathrm{Fl}$ and $\\gamma_\\mathrm{Fl}$, so the complex spectrum is experimentally readable through the cavity output.","The Floquet method maps the excitation spectrum over a wide parameter range at low numerical cost, avoiding exact diagonalization of the full master equation."],"supporting_citations":[{"why":"Supplies the atom-only Lindblad master equation and the mean-field equations (4)–(6), including the coherent $V_0$ and dissipative $V_1$ terms that stabilize the time-crystal limit cycle.","marker":"[52]"},{"why":"Provides the experimental realization of the dissipative time crystal and the modulated-coupling scheme $g(t)=g_0+g_1\\cos(\\omega t)$ that the paper models.","marker":"[29]"},{"why":"Establishes the parametrically driven dissipative Dicke model and the parametric-resonance physics at $\\omega\\approx 2\\omega_\\mathrm{res}$ that underlies the transition.","marker":"[27]"},{"why":"Establishes discrete time-crystalline order in cavity QED systems, defining the spontaneous $\\mathbb{Z}_2$ and discrete time-translation symmetry breaking studied here.","marker":"[28]"},{"why":"Underlies the displaced-frame/Schrieffer-Wolff cavity elimination used to derive the atom-only description in the bad-cavity regime.","marker":"[51]"},{"why":"Supplies the Liouvillian-gap logic for classifying dissipative phase transitions, which the Floquet spectrum extends to periodically driven systems.","marker":"[47]"},{"why":"Provides the mean-field Floquet numerical method that makes the spectrum calculation efficient enough to map over a wide parameter range.","marker":"[48]"}],"fun_headline_variants":["A single collective frequency defines dissipative time crystal","Floquet spectrum reveals time crystal's complex transition signatures","Probe lines reveal time crystal's two-phase spectra","Floquet eigenvalues map time crystal's excitation spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assumption that the cavity responds much faster than the atoms, so that it can be eliminated from the dynamics entirely; if that time-scale separation fails, the atom-only equations and the predicted excitation spectrum are not guaranteed to match the full atom-cavity system.","fun_headline_variants_meta":{"raw":{"variants":["A single collective frequency defines dissipative time crystal","Floquet spectrum reveals time crystal's complex transition signatures","Probe lines reveal time crystal's two-phase spectra","Floquet eigenvalues map time crystal's excitation spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000572,"raw_usage":{"total_tokens":2674,"prompt_tokens":889,"completion_tokens":1785,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":1722}},"tokens_in":505,"tokens_out":1785,"duration_ms":16825,"temperature":1.0,"reasoning_tokens":1722,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:53:30.281982+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to compare the atom-only Floquet spectrum with the probe response of the full dissipative Dicke model when the cavity is not in the bad-cavity regime (for example, with cavity linewidth $\\kappa$ comparable to $\\omega_\\mathrm{res}$, so that $\\omega\\tau_c \\sim 1$); a mismatch between the predicted and measured Lorentzian centers and widths in that regime would show that the adiabatic elimination is the load-bearing assumption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the atom-only Lindblad master equation and the mean-field equations (4)–(6), including the coherent $V_0$ and dissipative $V_1$ terms that stabilize the time-crystal limit cycle."},{"cited_title":"Chitra and O","cited_arxiv_id":null,"evidence_quote":"Establishes the parametrically driven dissipative Dicke model and the parametric-resonance physics at $\\omega\\approx 2\\omega_\\mathrm{res}$ that underlies the transition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the mean-field Floquet numerical method that makes the spectrum calculation efficient enough to map over a wide parameter range."}],"review_version":1}