{"id":"9df14aff-24c3-4e6c-b947-f0afe5818ce9","arxiv_id":"2506.18622","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Rapid parametric driving of a dissipative O(N) field theory creates an effective pump that can drive a transition into a rotating limit-cycle time-crystalline phase.","lead":"This paper shows that rapidly shaking a magnetic material's interactions, without breaking any symmetry, can make its ordered spins rotate on their own, forming a slow time-crystal phase. The authors build a general field-theory framework for this effect, giving experimenters a potential recipe for making time crystals in pumped magnets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite-width reservoir may flip the sign of the antidamping; the delta-shell ansatz in the SM is the key unverified assumption.","rationale":"The reader's weakest assumption correctly identifies the sharply peaked reservoir occupation and the leading-order loop expansion as the fragile points of the derivation. My stress test sharpens this into a decisive and concrete check: the sign of the sunset self-energy must be recomputed for a finite-width momentum distribution. The paper's own numerics simulate the full time-dependent equation and do find a limit cycle, which gives indirect support, but they do not measure the effective damping or the sign of δγ, so they cannot rescue the analytic result if the finite-width calculation flips the sign. The prefactor discrepancy between Eq. (4) and the SM is a real inconsistency that further weakens confidence in the quantitative condition for the instability, though it does not by itself change the qualitative picture. Since the reader has already conditioned acceptance on control of this calculation, my recommendation is to keep the verdict unchanged: the paper remains plausible and interesting, but the central sign should be validated by the proposed finite-width computation before full acceptance.","tokens_in":13654,"tokens_out":9114,"duration_ms":99683,"concrete_test":"Recompute the retarded self-energy derivative of SM Eq. (S13) with a normalized Lorentzian momentum occupation n(q) = n_b κ / [π((|q| − qΩ)^2 + κ^2)] in place of the delta shell of SM Eq. (S5), for N = 3, scanning κ/γΩ from 0.1 to 10 and γΩ/Ω from 0.01 to 0.3. If Im ∂ωΣR(0,0) changes sign for any physically relevant width, the antidamping mechanism fails. As a complementary check, simulate Eq. (3) and measure the effective low-frequency damping by the linear response of the slow mode to a weak impulse; if γeff − γ does not become negative in a region of the phase diagram, the central claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result δγ < 0 is derived in the Supplemental Material from the sunset integral (SM Eqs. S6–S16) under the ansatz n(q) = n_b δ(|q|−qΩ), together with the ordering Ω ≫ γΩ ≫ γ and a sharp momentum shell. This calculation is not controlled for a reservoir of finite width: once the delta is replaced by a spectral function of width κ, the momentum convolution q1 + q2 + q3 = 0 can receive contributions where the third mode is far from zero, and the frequency convolution ω1 + ω2 + ω3 = 0 can acquire imaginary parts of opposite sign. The numerical limit cycle of Fig. 2 is indirect evidence, but it does not measure δγ directly. In addition, the main-text formula Eq. (4) and the SM result Eq. (S16) do not agree in their prefactor dependence (Eq. (4) carries λ², while Eq. (S16) carries an explicit γ/γΩ²), so the quantitative condition |δγ| > γ is not reliably established. Because the entire time-crystal mechanism rests on the sign of δγ, the finite-width behavior of the sunset integral is the most load-bearing weak point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a field-theoretic framework for parametrically driven dissipative O(N) models. The authors argue that a rapid parametric drive creates a highly occupied high-momentum reservoir; after integrating out these fast modes, the low-frequency sector acquires an effective antidamping δγ<0 and a heating term ΔD>0. When the bare damping γ is sufficiently small, the effective damping γeff=γ+δγ crosses zero, triggering an instability to a time-crystalline limit-cycle phase with a rotating order parameter. The analytic derivation in the Supplemental Material uses a sunset diagram with a sharply peaked reservoir occupation, and numerical simulations of the corresponding Langevin equation in the ordered phase (r<0) show a limit cycle with the predicted amplitude dip and square-root frequency onset. The paper also discusses universality classes, including critical exceptional points and the noisy Hopf bifurcation, and connects the mechanism to pumped magnetic materials.","tokens_in":13759,"tokens_out":5275,"duration_ms":52707,"significance":"If the central antidamping result is correct, this is a significant theoretical contribution: it proposes a generic route to continuous time-crystalline order in dissipative bosonic field theories, with plausible relevance to laser-pumped magnets and magnonics. The Supplemental Material provides a structured analytic derivation, and the numerics reproduce the qualitative predictions (amplitude dip, square-root frequency onset), which gives the package some coherence. The falsifiable predictions for driven magnets and the connection to prior critical-exceptional-point universality classes are valuable. However, the correctness of the mechanism depends on the sharply peaked reservoir ansatz and on the leading-order loop expansion, and the quantitative inconsistency between the main-text formula and the SM expression must be resolved before the central claim is fully supported.","major_comments":[{"comment":"The sign of the central antidamping result is derived under the sharply peaked reservoir occupation ansatz n(q)≈n_b δ(|q|−qΩ) (SM text above Eq. (S1) and Eq. (S5)). The calculation of the sub-integral I2 in Eq. (S11) then assumes a single relaxation rate γΩ and introduces a shell width δq without a precise definition. If the reservoir spectral function has finite width κ, the momentum convolution q1+q2+q3=0 in Eq. (S6) can receive significant contributions from configurations in which the third mode is far from zero momentum, and the frequency convolution can acquire imaginary parts of opposite sign. Since the entire time-crystal mechanism rests on the sign of δγ, the finite-width behavior of the sunset integral must be shown not to flip the sign; currently no such control is provided.","section":"Supplemental Material, Eqs. (S6)–(S16)"},{"comment":"The main-text result Eq. (4) gives δγ = −Nγ λ² (qΩ n_b)², while the detailed SM expression (S16), with the prefactor Nγ defined immediately after, gives δγ ≈ −Nγ (qΩ n_b)² / γΩ² and contains no explicit λ². These two expressions differ in their parametric dependence on λ and γΩ, so the quantitative condition |δγ| > γ used in the main text is not backed by a single consistent formula. The authors should reconcile the prefactors and specify which expression underlies the phase diagram of Fig. 1b.","section":"Main text Eq. (4) and SM Eq. (S16)"},{"comment":"The analytic elimination of the reservoir is performed in the symmetric phase r>0 (main text, Eq. (4), and SM), whereas the extension to the ordered phase r<0, which is the regime relevant for the antiferromagnetic application and for the simulations of Fig. 2, is argued by analogy and not derived. The numerical observation of a limit cycle in Fig. 2 is indirect evidence: it does not directly measure δγ. A derivation for r<0, or at least a controlled argument that the same loop integral survives and retains its sign in the ordered phase, would materially strengthen the central claim.","section":"Main text and numerics for r<0"}],"minor_comments":[{"comment":"The section heading 'ELIMINA TION OF HIGH FREQUENCY BATH' contains a typo ('ELIMINA TION') and should be corrected to 'ELIMINATION'.","section":"Supplemental Material title"},{"comment":"The quantity δq appears after the delta-function ansatz n(q)≈n_b δ(q−qΩ); its definition as the width of the momentum shell should be stated explicitly and used consistently with Eq. (S5).","section":"SM Eq. (S11)"},{"comment":"The caption contains an errant space in 'Ω = 2 .15'; additionally, specifying that the simulations are performed in three dimensions would improve clarity.","section":"Fig. 2 caption"},{"comment":"The sentence 'It then carries through as a parametric drive into the continuum limit' is grammatically awkward and could be rephrased for clarity.","section":"Main text, paragraph on pumped magnetic materials"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a letter with a substantial Supplemental Material. The mechanism is interesting and timely, but the sign of the antidamping is the load-bearing result, and its dependence on the sharp-peak ansatz and finite reservoir width is not controlled. The quantitative mismatch between Eq. (4) and SM Eq. (S16) should also be resolved. These are fixable within the scope of the manuscript, so I recommend major revision rather than rejection. No concerns about citation practice or overlap beyond the authors' own prior work, which is appropriately referenced."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know upfront. The paper has a real and new idea: rapid parametric pumping creates a high-momentum reservoir which, after integrating out, acts as an incoherent pump for long-wavelength modes, giving antidamping δγ < 0 and driving a continuous time crystal in a dissipative O(N) model. Second, the central sign δγ is computed under a sharply peaked reservoir ansatz, and the derivation is not yet shown to be robust; there is also a prefactor mismatch between the main text and the SM.\n\nThe genuinely new piece is the two-loop sunset calculation of δγ from a parametrically driven reservoir, and the crisp claim that this is a generic route to limit-cycle time-crystalline order. The numerics qualitatively confirm the limit cycle with the predicted amplitude dip and square-root frequency onset. The universality class discussion reuses their prior work, but that is a legitimate extension, not circular.\n\nThe biggest soft spot is the delta-shell occupation n(q) ≈ δ(q−qΩ)n_b. The sunset integral that yields δγ is delicate; once the reservoir has finite width, the momentum and frequency convolutions can pick up contributions with opposite imaginary parts, and the sign of δγ is no longer guaranteed. The ordering Ω ≫ γΩ ≫ γ makes the shell narrow, but the authors should either do a finite-width calculation or give a controlled argument that the sign is stable. As it stands, the sign of δγ is the load-bearing result and its robustness is not established.\n\nSecond, there is an actual inconsistency in the prefactor: main-text Eq. (4) has λ², while SM Eq. (S16) appears to omit it. This matters because the condition |δγ| > γ is the whole point. A reader cannot tell which version is correct.\n\nThird, the mapping from drive amplitude r_D to reservoir occupation n_b is not derived. Treating n_b as a phenomenological parameter is fine, but then the quantitative link to specific experiments stays open.\n\nFourth, the loop expansion is leading order only. That is a standard limitation, but with a highly occupied reservoir, higher orders could matter; a sentence on that would help.\n\nThe numerics are good qualitative evidence, though there are no error bars or data files. For a letter that is acceptable; error bars on the dip and frequency would strengthen it.\n\nThis paper is for people working on time crystals, driven magnonics, and nonequilibrium critical phenomena. It deserves a serious referee. I would send it to peer review, expecting revision to address the prefactor inconsistency and the finite-width robustness. The idea is likely to be influential even if the details need tightening.","headline":"A genuinely new mechanism for continuous time crystals via parametrically pumped antidamping, but the sign of δγ rests on a sharp-peak ansatz and a prefactor mismatch in the SM needs fixing.","tokens_in":14386,"tokens_out":6179,"would_cite":true,"duration_ms":52334,"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":"Rapid parametric driving of a dissipative O(N) field theory can reverse the effective damping of long-wavelength modes, driving a transition into a continuous time-crystalline limit-cycle phase.","keywords":["time crystals","parametric driving","Keldysh field theory","driven-dissipative systems","O(N) models","limit cycles","antiferromagnets","nonequilibrium phase transitions"],"falsifier":"Compute the next-order correction to $\\delta\\gamma$ for a reservoir of finite width $\\delta q$: if the negative sign survives only for $\\delta q/q_\\Omega$ below a few percent, the mechanism fails for realistic pumps. Experimentally, measure the $q=0$ magnon linewidth versus drive power in a low-damping antiferromagnet; the prediction is that the linewidth shrinks to zero and re-emerges with opposite sign at the onset of the rotating phase, while the rotation frequency rises as the square root of the pump strength above threshold.","tokens_in":13329,"feed_emoji":"🕰️","tokens_out":11878,"duration_ms":101175,"temperature":0.7,"pith_summary":"The paper aims to show that a fast, symmetry-preserving parametric drive—implemented as an oscillating mass, the generic way a laser couples to an insulator—acts as an effective pump for the slow modes of a dissipative O(N) order-parameter theory. The drive creates a sharply occupied reservoir of high-momentum modes, and nonlinear scattering makes those modes 'rain down' into the long-wavelength sector, contributing an antidamping $\\delta\\gamma<0$ and a heating term $\\Delta D>0$. When the bare damping is small enough, $\\gamma_{\\rm eff}=\\gamma+\\delta\\gamma$ crosses zero and the steady state becomes a rotating limit cycle: a continuous time crystal whose period is much longer than, and incommensurate with, the drive period. The authors derive the loop corrections analytically, draw the phase diagram, and support the ordered-phase case with lattice simulations. A sympathetic reader would care because this gives a generic, symmetry-preserving route to time-crystalline order in pumped magnetic materials.","feed_headline":"Rapid drive can push a dissipative magnet into a time crystal","feed_subtitle":"Pumping high-frequency modes flips a magnet's effective damping and drives a rotating limit cycle.","key_machinery":"The load-bearing construction is a two-step coarse-graining. Split the field into slow low-momentum modes and fast reservoir modes peaked at momentum $q_\\Omega$ with occupation $n(q)\\approx \\delta(q-q_\\Omega)n_b$, then integrate the reservoir out in a leading-order loop expansion of the response-functional path integral equivalent to the Langevin equation. The damping shift is read off from the sunset (two-loop) self-energy diagram, $\\delta\\gamma=\\operatorname{Im}\\,\\partial_\\omega \\Sigma^R_p(\\omega=0,q=0)$, in which two reservoir correlation functions and one retarded response function meet; the analytic continuation of this sunset integral produces the negative sign. The heating shift $\\Delta D$ comes from the same diagram with three reservoir correlation functions. The generated nonlinear damping vertices $u(\\tilde\\phi\\cdot\\partial_t\\phi)(\\phi\\cdot\\phi)$ and $u'(\\tilde\\phi\\cdot\\phi)\\partial_t(\\phi\\cdot\\phi)$ carry the stability of the limit cycle, with $u'>u>0$ enforcing rotation rather than amplitude oscillation.","core_discovery":"The central claim is that integrating out the fast reservoir created by the drive leaves the low-frequency sector with an effective damping $\\gamma_{\\rm eff}=\\gamma+\\delta\\gamma$, where $\\delta\\gamma=-N_\\gamma \\lambda^2 (q_\\Omega n_b)^2<0$, with $n_b$ the reservoir occupation set by pump power and $q_\\Omega$ the resonant momentum. For sufficiently small bare $\\gamma$, $\\gamma_{\\rm eff}$ crosses zero and the symmetric (or ordered) state becomes unstable toward a limit cycle in which the order parameter rotates at frequency $\\omega_0$ much smaller than the drive frequency; the paper gives the explicit rotating solution for $N=3$. The same reservoir integration generates positive noise $\\Delta D\\propto \\lambda^2 N_D (q_\\Omega n_b)^3/\\Omega^2$ and nonlinear damping terms $u,u'>0$ with $u'>u$, which stabilize the rotating phase and select rotation over amplitude oscillations. Too strong a drive therefore heats the system back into a paramagnet, giving a finite window of time-crystalline order in the pump-power–temperature plane. The transition from the ordered phase is claimed to occur through a critical exceptional point in the O(3)/antiferromagnetic case and through a noisy Hopf bifurcation for the SO(3)/ferromagnetic precession case.","pith_inferences":["A clean experimental discriminator would be the scaling $\\omega_0\\sim\\sqrt{\\text{pump power}-\\text{threshold}}$ together with the amplitude dip at the transition; these features separate the antidamping mechanism from ordinary parametric amplification, which does not produce a self-sustained rotating condensate.","Since $\\delta\\gamma$ grows linearly with $n_b^2$ while heating grows with $n_b^3$, materials with the smallest intrinsic magnetic damping should show the broadest time-crystal window; drive frequency may be most effective slightly off exact resonance where scattering phase space is larger.","The same effective-field-theory reduction may apply to driven exciton-polariton or phonon systems in which an optically hot band pumps a slow condensate, even when the drive does not address the condensate mode directly; the key ingredients are nonlinear scattering between fast and slow bands plus a relaxing bath.","A test of the leading-order treatment is to measure how the effective linewidth and noise scale separately with pump power: the predicted ratio $\\delta\\gamma^2/\\Delta D \\propto n_b$ may acquire corrections at higher occupation, revealing whether the sunset approximation is the whole story."],"forward_implications":["A symmetry-preserving parametric pump is enough to produce continuous time-crystalline order in dissipative bosonic field theories; no explicit symmetry breaking or engineered nonreciprocity is required.","The paramagnet-to-time-crystal transition falls into a nonthermal universality class whose correlations violate fluctuation-dissipation relations, while the ordered antiferromagnet case is governed by a critical exceptional point with a dip in the order-parameter amplitude and $\\omega_0\\propto\\sqrt{r_D-r_{D,c}}$.","In an SO(3) ferromagnet the induced precession term changes the ordered-side transition into a noisy Hopf bifurcation, described by the same universality class as a driven-dissipative complex order-parameter field.","The mechanism transfers to any setting with a high-frequency band nonlinearly coupled to slow modes and to a thermal bath, including the strongly driven Hubbard model, where the numerically seen 'magnon condensation' would be the onset of this time crystal.","Very strong drives heat the low-frequency sector through $\\Delta D\\propto n_b^3$, so time-crystalline order occupies a finite window of pump power and low bath temperature, with a re-entrant paramagnet at high power."],"supporting_citations":[{"why":"Supplies the Keldysh path-integral formalism and the dissipative action used as the starting point in Eq. (2).","marker":"[48]"},{"why":"Provides the driven-open-system effective-action framework and the dissipative O(N) Lagrangian.","marker":"[22]"},{"why":"Established the rotating limit-cycle phase and critical-exceptional-point phenomenology that the paper builds on.","marker":"[41]"},{"why":"Identified the nonthermal O(3) universality class at the onset of time-crystalline order used for the symmetric-phase transition.","marker":"[56]"},{"why":"Provides Floquet-RPA evidence of a magnetic transition in a driven Hubbard model that the paper interprets as the onset of its time crystal.","marker":"[54]"},{"why":"Supplies numerical evidence for a nonequilibrium phase transition in a driven-dissipative antiferromagnet matching the mechanism.","marker":"[55]"},{"why":"Supplies the noisy-coupled-oscillator universality class used for the SO(3) ferromagnet-to-time-crystal transition.","marker":"[58]"}],"fun_headline_variants":["Rapid param drives yield time-crystalline phases in O(N) models","Shaking magnets at high frequency creates a time crystal","Negative effective damping from pump powers time-crystalline order","Fast pump flips dissipation, stabilizes rotating order","Parametric driving induces limit cycle: a time crystal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on the ansatz that the pumped reservoir occupation is sharply peaked at a single momentum, $n(q)\\approx\\delta(q-q_\\Omega)n_b$, with a narrow width and a decay rate $\\gamma_\\Omega$, and on computing only leading-order (one- and two-loop) corrections; if the reservoir is spread in momentum or higher-order loops matter, the claimed antidamping $\\delta\\gamma<0$ could change sign or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Rapid param drives yield time-crystalline phases in O(N) models","Shaking magnets at high frequency creates a time crystal","Negative effective damping from pump powers time-crystalline order","Fast pump flips dissipation, stabilizes rotating order","Parametric driving induces limit cycle: a time crystal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000205,"raw_usage":{"total_tokens":1377,"prompt_tokens":910,"completion_tokens":467,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":387}},"tokens_in":526,"tokens_out":467,"duration_ms":5380,"temperature":1.0,"reasoning_tokens":387,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:45:45.470588+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the next-order correction to $\\delta\\gamma$ for a reservoir of finite width $\\delta q$: if the negative sign survives only for $\\delta q/q_\\Omega$ below a few percent, the mechanism fails for realistic pumps. Experimentally, measure the $q=0$ magnon linewidth versus drive power in a low-damping antiferromagnet; the prediction is that the linewidth shrinks to zero and re-emerges with opposite sign at the onset of the rotating phase, while the rotation frequency rises as the square root of the pump strength above threshold.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Established the rotating limit-cycle phase and critical-exceptional-point phenomenology that the paper builds on."},{"cited_title":"Daviet, C","cited_arxiv_id":null,"evidence_quote":"Identified the nonthermal O(3) universality class at the onset of time-crystalline order used for the symmetric-phase transition."},{"cited_title":"Walldorf, D","cited_arxiv_id":null,"evidence_quote":"Provides Floquet-RPA evidence of a magnetic transition in a driven Hubbard model that the paper interprets as the onset of its time crystal."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies numerical evidence for a nonequilibrium phase transition in a driven-dissipative antiferromagnet matching the mechanism."},{"cited_title":"Risler, J","cited_arxiv_id":null,"evidence_quote":"Supplies the noisy-coupled-oscillator universality class used for the SO(3) ferromagnet-to-time-crystal transition."}],"review_version":2}