{"id":"8e477ff5-6218-41ba-b4b3-55f53bdb2f0c","arxiv_id":"2507.18950","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Quenches into dissipative discrete time crystals show Kibble-Zurek defect scaling, with a universality class set by the dissipative linear parametric oscillator's vz=1.","lead":"This paper shows that defects formed when a system is rapidly driven into a discrete time crystal follow the same universal power laws as ordinary phase transitions. The result extends the Kibble-Zurek mechanism to time crystals and identifies a damped oscillator model that predicts when this universality appears.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The universality claim rests on the asserted DLPO mapping for the SGM and DLM, which is cited to a same-group paper rather than derived or numerically verified here; if the linearized relaxation rate near the transition does not vanish as |A-Ac|, the vz=1 prediction and the inferred KZM scaling…","rationale":"The paper has real strengths: the DLPO analysis gives an explicit mechanism for vz=1, the KZM scaling of t_hat, xi, and n_d is tested in two very different models, averages over 100 noise realizations are used, and finite-size effects are checked. I do not see an internal inconsistency in the DLPO calculation itself; the multi-scale derivation is standard, and the threshold Ac=2gamma/Omega follows when damping is balanced against the parametric drive. The weakest point is the leap from a single parametric oscillator to the many-body lattice models. That leap is load-bearing because the title and abstract claim a universality class, not merely a curve collapse for two examples. A direct linear-stability test would settle it: if the Floquet exponent is linear in Ac-A in both models, the concern is resolved and the central claim is materially supported. The additional observations about hand-set thresholds, truncated Wigner approximation, and the differing v,z values are worth reporting, but they are secondary and would lower confidence rather than break the argument. Because the reader already marked the paper CONDITIONAL on essentially this issue, my stress-test does not move the verdict.","tokens_in":18859,"tokens_out":11110,"duration_ms":123687,"concrete_test":"Linearize the SGM equation (8) around theta=0 and the truncated-Wigner DLM equations (B2) around a=0, S^z=-N/2, S^x=S^y=0, using the same parameters and k_r values (0 or pi) as in Fig. 4. For fixed drive amplitude A below Ac, compute the dominant Floquet multiplier of the linearized system over one drive period and extract the relaxation rate lambda(A). Repeat for several A approaching Ac and fit lambda proportional to (Ac-A)^theta. If theta=1 for both models and both DTC configurations, the DLPO mapping and the predicted vz=1 are confirmed; if theta differs from 1, the analytic basis for the universality claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central statement is that any DTC-forming system mappable onto a DLPO obeys the adiabatic-impulse approximation with vz=1, and that the SGM and DLM belong to this class. The DLPO calculation in Sec. II and Appendix A is explicit for a single oscillator, but the required extension to the finite-lattice SGM and DLM is not derived in this paper. The text delegates the mapping to Ref. [78] (same group) and, for the DLM, to Refs. [43,77]; no effective DLPO parameters are written down, and no check is given that the soft spatial mode at k_r is decoupled from the other modes during the ramps used in Fig. 4. The numerical evidence that beta_t is close to 1/2 is a consistency check with the assumed mapping, not independent support: the same measured exponent could in principle arise even if the mapping were inaccurate. Concretely, if the linearized relaxation rate of either model near the transition were to vanish as |Ac-A|^theta with theta different from 1, or with a logarithmic correction, the analytic prediction vz=1 would be inapplicable and the universality-class claim would lose its foundation. The paper's own extracted v and z differ between the SGM and the DLM and between FM and AFM configurations; this does not falsify the DLPO mapping, but it underscores that the claimed shared universality class rests entirely on the single product vz, which is exactly the quantity imported from the unverified mapping.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the Kibble-Zurek mechanism (KZM) applies to quenches from a disordered phase into a discrete time crystal (DTC), and that systems mappable onto a dissipative linear parametric oscillator (DLPO) form a universality class with critical exponent product vz = 1. An analytic multi-scale calculation for the DLPO shows the relaxation time diverges as |A - A_c|^{-1}, and this is used to predict transition-delay, correlation-length, and defect-number power laws. The authors test the prediction numerically in two one-dimensional models: the classical Sine-Gordon model (SGM) and the open Dicke lattice model (DLM) simulated with a truncated Wigner approximation. For both models and for both ferromagnetic and antiferromagnetic DTC configurations, they report scaling exponents consistent with the predicted KZM forms, and they extract v and z separately, finding values near the mean-field Ising class for the SGM and somewhat different values for the DLM.","tokens_in":19173,"tokens_out":4245,"duration_ms":49016,"significance":"If the central claim is correct, the work would meaningfully extend the KZM to dynamical (spatiotemporal) order and offer a concrete, falsifiable criterion: any DTC-forming system that can be mapped to a DLPO should exhibit vz = 1 scaling. The analytic DLPO derivation in Sec. II and Appendix A is clean and self-contained, and the numerical study of two very different models (classical coupled pendula and a dissipative spin-boson lattice) is a valuable consistency test. The paper also explicitly identifies the regime of fast-quench breakdown and checks the relation between correlation length and defect number, which strengthens the case for KZM behavior. However, the universality-class claim is only as solid as the asserted DLPO mapping for the SGM and DLM, and that mapping is not derived here; moreover, the numerical extraction of the critical point partly assumes the scaling it is used to verify. These issues are load-bearing for the paper's main conclusion, so the present version is not yet fully convincing.","major_comments":[{"comment":"The central universality claim requires that the finite-lattice SGM and DLM near their DTC transitions be described by a DLPO whose slow mode has relaxation time diverging as |A - A_c|^{-1}. This mapping is not derived in the manuscript: the general statement is delegated to Ref. [78] (same group) and the DLM mapping to Refs. [43,77], and no effective DLPO parameters (effective gamma, Omega, and the coupling of the soft spatial mode to other modes) are provided for either model. If the linearized relaxation rate of the soft mode vanishes as |A - A_c|^theta with theta different from 1, or with a logarithmic correction, the prediction vz = 1 and the shared-universality-class conclusion would not follow. I ask the authors to either derive the mapping for the two lattice models or provide a direct numerical verification, for example by measuring the exponential decay rate of small perturbations toward the disordered state as a function of A - A_c.","section":"Sec. II, Eq. (6), and Sec. III"},{"comment":"The critical point A_c used to compute t_c is obtained by fitting <A(t_p)> to the KZM-derived functional form A(tau_q) = A_0 tau_q^{-1/(1+vz)} + A_c. The same KZM scaling is then used to measure beta_t, beta_xi, and beta_nd, so the reported consistency with vz = 1 is partially built into the analysis: the extrapolated A_c can shift t_c in a way that favors the assumed exponent. Please determine A_c independently (for instance, from the static phase boundary or from the divergence of the relaxation time), or at least report the fitted exponent b from Eq. (C2) and show that the extracted A_c is stable when the assumed fitting form is changed.","section":"Appendix C, Eq. (C2), and Sec. IV"},{"comment":"The DLM exponents are only weakly constraining for the claimed universality class: the paper reports vz_FM = 0.8(5) and vz_AFM = 0.8(2), so consistency with vz = 1 is established only within large error bars, and the extracted v and z differ between models and between FM and AFM configurations. The manuscript should state the confidence intervals on beta_t for each model and configuration and quantify whether the differences in v and z between the SGM and DLM are statistically significant. As written, the separate v and z values do not independently support a common universality class; only the product vz is used, and that product is the quantity imported from the asserted DLPO mapping.","section":"Sec. V, Figs. 5(b) and 6"}],"minor_comments":[{"comment":"The caption says '(a)-(b)' but the panels are labeled (a)-(d), and the text refers to '(c) the SGM' and '(d) the DLM'; the caption should be corrected to match the panel labels.","section":"Fig. 6 caption"},{"comment":"Equation (A1) ends with a stray prime after 'theta = 0'; this appears to be a typo and should be removed.","section":"Eq. (A1)"},{"comment":"The threshold delta = 0.15 is arbitrary; please report a sensitivity check (e.g., delta = 0.1 and 0.2) to show that the extracted scaling exponents do not depend on the chosen threshold.","section":"Sec. IV, Eq. (14)"},{"comment":"The dashed lines in Fig. 4 correspond to v = 1/2, z = 2, i.e., the mean-field Ising universality class. Since the DLPO prediction fixes only the product vz = 1, the dashed lines are not predictions of the DLPO universality class alone but involve an additional assumption about v; this should be clarified in the text.","section":"Fig. 4 and Sec. IV"},{"comment":"Reference [47] contains a typo: 'time crsytals' should be 'time crystals'.","section":"References"},{"comment":"The phrase 'Given that our results here are in the regime in which the noise is sufficiently weak to preserve the coherence of the DTCs' is grammatically awkward; consider revising for clarity.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"This is an interesting and well-written manuscript with a clean analytic core, but the main claim is currently supported by a mapping cited to earlier same-group work rather than derived or verified here, and the numerical analysis contains a partial circularity in the extraction of the critical point. Both issues are fixable within the scope of a revision, so I do not recommend rejection. The editor may also wish to check that the journal's standards permit the use of a same-author paper as the sole basis for the key mapping step."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives the first clear evidence that the Kibble-Zurek mechanism (KZM) applies when a driven-dissipative system is quenched into a discrete time crystal phase. That is new and worth taking seriously. The analytic piece — the dissipative linear parametric oscillator (DLPO) relaxation time diverging as |A-Ac|^{-1}, so vz=1 — is clean, and the numerics for the sine-Gordon and Dicke lattice models show the expected power laws with a sensible saturation at fast quenches. The authors deserve credit for testing both a classical and a semiclassical quantum model.\n\nThe soft spot is proportionate to the strength of the claim. The paper wants to establish a universality class for DTC-forming systems, but the evidence for the universality class reduces to the single product vz. That product is not fitted freely; it is imported from the DLPO mapping, and the mapping from the lattice models to the DLPO is delegated to a same-group reference rather than derived or numerically verified here. So the observed beta_t ~ 1/2 is a consistency check with the assumed mapping, not independent confirmation. The Dicke lattice exponents also come with large error bars (vz=0.8(5) and 0.8(2)), compatible with 1 but not a tight test. There are also smaller issues: the transition-time threshold delta=0.15 is arbitrary, the critical point A_c is extracted using a fit to the KZM form, and no code or data are released. None of these look fatal. The central observation — that quenches into a DTC show KZM scaling — appears solid. What is not solid is the stronger universality class claim.\n\nWho is this for? People working on time crystals, parametric oscillator arrays, and KZM in driven-dissipative systems. They will get a clear, useful numerical study plus a testable prediction (vz=1) for other DLPO-mappable systems. I would engage with the paper in review, but I would ask the authors to either derive or directly compute the linearized relaxation rate near the transition for the lattice models, or to present the universality claim as a scaling observation rather than a proven class. Three months of additional numerics could tighten the DLM exponents.\n\nRecommendation: send it to peer review. It deserves referee time; the structural weakness is addressable and the core scaling result is likely right.","headline":"First credible KZM scaling evidence for DTC formation, but the universality class claim outruns the evidence.","tokens_in":19754,"tokens_out":2178,"would_cite":true,"duration_ms":22079,"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":"This paper claims that quenches into discrete time crystals obey Kibble-Zurek scaling set by a single damped oscillator class, with both classical and quantum lattice models falling into the same universality class.","keywords":["discrete time crystals","Kibble-Zurek mechanism","dissipative linear parametric oscillator","adiabatic-impulse approximation","universality class","Sine-Gordon model","Dicke lattice model","dynamical phase transition"],"falsifier":"Measure the relaxation time directly by holding the system just below the period-doubling threshold, perturbing it, and observing the decay rate across a range of drive amplitudes near $A_c$; if the divergence exponent is not 1, the predicted $\\beta_{\\hat t}=1/2$ should fail. Alternatively, in any DTC-forming array, extract the transition-delay exponent from a slow linear ramp of drive amplitude: a value measurably different from $1/2$ in the scaling regime would falsify membership in the DLPO universality class.","tokens_in":1963,"feed_emoji":"🔁","tokens_out":2249,"duration_ms":86422,"temperature":0.7,"pith_summary":"The paper argues that forming a discrete time crystal by driving a dissipative many-body system through its critical point is a genuine phase transition governed by the Kibble-Zurek mechanism, the standard scaling theory for how defects appear when a system is quenched across a transition. Its central claim is that any system that can be mapped onto a single damped, parametrically driven harmonic oscillator should show the same universal scaling: the delay before time-crystal order appears and the number of spatial defects both follow power laws in the quench time, with the product of the static and dynamical critical exponents equal to one. The authors verify this in a classical chain of driven pendula and in a quantum array of lossy spin-cavity systems. If correct, the result extends the notion of universality from static phases to dynamical, periodically ordered phases and gives concrete predictions for any experiment that quenches into a time crystal.","feed_headline":"Time crystal formation follows Kibble-Zurek scaling","feed_subtitle":"Classical pendulum chains and quantum cavity lattices share the same transition delay and defect powers.","key_machinery":"The central object is the dissipative linear parametric oscillator (DLPO), a damped harmonic oscillator with periodically modulated drive amplitude. In the resonance limit, a multi-scale analysis yields relaxation times $\\tau_\\pm = 8(A_c \\pm A)^{-1}$, and the slower one diverges at the critical drive amplitude $A_c = 2\\gamma/\\Omega$. That divergence is exactly what the adiabatic-impulse approximation requires, and it sets the exponent product $vz=1$; through the Kibble-Zurek scaling formulas, it determines the power laws for the transition delay, correlation length, and defect density. The DLPO serves as a template: any lattice or cavity model mappable onto it is predicted to belong to the same universality class.","core_discovery":"On the paper's own terms, the discovery is that the adiabatic-impulse approximation—the assumption that a system freezes as its relaxation time diverges while a control parameter is ramped—holds for the dissipative linear parametric oscillator (DLPO), the damped harmonic oscillator whose restoring force is modulated at twice its natural frequency. Near resonance, the oscillator's relaxation time diverges as $(A_c-A)^{-1}$ as the drive amplitude approaches the period-doubling threshold, giving $vz=1$. Since many DTC-forming models reduce to a DLPO, any such system quenched from disorder into a discrete time crystal should exhibit Kibble-Zurek scaling: transition delay proportional to $\\tau_q^{1/2}$, correlation length proportional to $\\tau_q^{v/(1+vz)}$, and defect number proportional to $\\tau_q^{-v/(1+vz)}$, which in one dimension with point defects gives exponents $1/4$ and $-1/4$. The paper shows numerically that both the classical Sine-Gordon model and the open Dicke lattice follow these laws, and that their measured $\\beta_{\\hat t}=1/2$ and $\\beta_{n_d}/\\beta_\\xi=1$ place them in the same universality class even though their individual static and dynamical exponents differ.","pith_inferences":["If the DLPO mapping is as general as the paper claims, existing experimental arrays of coupled nanomechanical resonators, superconducting parametric oscillators, and atom-cavity lattices should exhibit the same $\\tau_q^{1/2}$ transition-delay law when driven across their period-doubling threshold; this is a direct test that does not require measuring full correlation functions.","Because the universality class is fixed by the product $vz$, measuring the transition-delay exponent alone may suffice to classify a DTC-forming system, simplifying experimental protocols.","In the strong-noise regime, the closed-loop spatiotemporal defects the paper observes arise from fluctuations, not from Kibble-Zurek correlation build-up, so their density should not follow quench-time scaling; this suggests a separate theory is needed for fluctuation-seeded defects.","For effectively zero-dimensional all-to-all coupled systems, the same reasoning predicts only the temporal part of the Kibble-Zurek mechanism (transition delay, with no spatial defect network), so single-mode cavity experiments could test the universality class by measuring delay scaling alone."],"forward_implications":["Any open DTC-forming system mappable onto a DLPO, classical or quantum, should show a transition delay scaling as $\\tau_q^{1/2}$ in the slow-quench regime, independent of microscopic details.","Defect production and correlation growth during DTC formation obey the Kibble-Zurek mechanism: in one dimension with point defects, defect number scales as $\\tau_q^{-1/4}$ while correlation length grows as $\\tau_q^{1/4}$, with $\\beta_{n_d}/\\beta_\\xi = 1$.","The classical Sine-Gordon model and the open Dicke lattice belong to the same universality class, characterized by $vz=1$, even though their individual $v$ and $z$ values differ; universality therefore extends to spatiotemporally ordered dynamical phases.","Fast quenches should show a universal breakdown of scaling, with delay, correlation length, and defect number saturating at finite values, consistent with the reference prediction.","The work implies that DTC formation is a genuine phase transition with a many-body character, not merely a single-oscillator nonlinear effect, because the observed scaling is set by collective critical exponents."],"supporting_citations":[{"why":"Supplies the mapping of a broad class of DTC-forming systems onto the dissipative linear parametric oscillator, the basis for the proposed universality class.","marker":"[78]"},{"why":"Provides the multi-scale analytic solution of the parametric oscillator used to obtain the relaxation-time divergence in Eq. (6).","marker":"[79]"},{"why":"States the adiabatic-impulse approximation and the Kibble-Zurek scaling relations that the paper verifies.","marker":"[5]"},{"why":"Defines the classical Sine-Gordon model of discrete time crystals and its thermal noise properties used in the numerics.","marker":"[48]"},{"why":"Introduces the open Dicke lattice model and its mapping to a parametric oscillator below threshold.","marker":"[77]"},{"why":"Prior derivation that quenched open systems show an apparent Kibble-Zurek delay; supplies the DLM-to-DLPO connection and the static-limit baseline.","marker":"[43]"},{"why":"Establishes the universal breakdown of Kibble-Zurek scaling in fast quenches, used to interpret the saturation at small quench times.","marker":"[42]"},{"why":"Provides the truncated Wigner approximation used to simulate the quantum Dicke lattice dynamics.","marker":"[80]"}],"fun_headline_variants":["Dissipative time crystals follow Kibble-Zurek scaling","Universal defect scaling in dissipative time crystal formation","Open-system time crystals enter Kibble-Zurek universality class","Parametric oscillator explains universal time crystal quench laws"],"cache_read_input_tokens":21760,"weakest_assumption_plain":"The whole argument rests on the assumption that near the transition the many-body Sine-Gordon and Dicke-lattice systems really are captured by the single-oscillator dissipative linear parametric oscillator description, so that their relaxation time diverges as $|A-A_c|^{-1}$; if that mapping is inaccurate, the predicted $vz=1$ scaling would not apply.","fun_headline_variants_meta":{"raw":{"variants":["Dissipative time crystals follow Kibble-Zurek scaling","Universal defect scaling in dissipative time crystal formation","Open-system time crystals enter Kibble-Zurek universality class","Parametric oscillator explains universal time crystal quench laws"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1521,"prompt_tokens":989,"completion_tokens":532,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":465}},"tokens_in":605,"tokens_out":532,"duration_ms":5424,"temperature":1.0,"reasoning_tokens":465,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:05:20.021754+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the relaxation time directly by holding the system just below the period-doubling threshold, perturbing it, and observing the decay rate across a range of drive amplitudes near $A_c$; if the divergence exponent is not 1, the predicted $\\beta_{\\hat t}=1/2$ should fail. Alternatively, in any DTC-forming array, extract the transition-delay exponent from a slow linear ramp of drive amplitude: a value measurably different from $1/2$ in the scaling regime would falsify membership in the DLPO universality class.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the mapping of a broad class of DTC-forming systems onto the dissipative linear parametric oscillator, the basis for the proposed universality class."},{"cited_title":"Kovacic, R","cited_arxiv_id":null,"evidence_quote":"Provides the multi-scale analytic solution of the parametric oscillator used to obtain the relaxation-time divergence in Eq. (6)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the classical Sine-Gordon model of discrete time crystals and its thermal noise properties used in the numerics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the open Dicke lattice model and its mapping to a parametric oscillator below threshold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior derivation that quenched open systems show an apparent Kibble-Zurek delay; supplies the DLM-to-DLPO connection and the static-limit baseline."},{"cited_title":"Zeng, C.-Y","cited_arxiv_id":null,"evidence_quote":"Establishes the universal breakdown of Kibble-Zurek scaling in fast quenches, used to interpret the saturation at small quench times."}],"review_version":2}