{"id":"10557bda-8637-4f6e-9381-dfadefe29612","arxiv_id":"2505.09117","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A Rydberg atom chain split into two halves driven at golden-ratio-related frequencies shows stable quasiperiodic time-crystalline order in numerical simulations.","lead":"The authors propose a way to create a discrete time quasi-crystal, a quantum state whose oscillations are ordered but never repeat exactly, in a chain of Rydberg atoms. They show numerically that driving two halves of the chain at incommensurate frequencies produces the predicted quasi-periodic response, and they map where the phase exists.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Missing decoupled control: absolute value in Eq. (3) can generate the claimed (1/2,1/2) and (-1/2,1/2) peaks from two independent subharmonic DTCs, so the interaction-causality claim is not yet established.","rationale":"The reader's weakest assumption identifies exactly the most load-bearing gap: the absence of a decoupled control and the nonlinearity of m. I agree with that assessment. The paper provides useful numerical evidence for a concrete Floquet protocol, including MPS truncation control, size scaling, and phase diagrams, and the observed quasiperiodic signal is plausibly real. However, the central interpretive claim that the aperiodic response is caused by Rydberg blockade interaction is not established by the present data, because all observables used to define DTQC are nonlinear functions of the separate subsystem signals. A straightforward decoupled simulation would settle the question. If the control shows that the sum and difference peaks persist without coupling, the paper's headline claim would be substantially weakened, reducing it to a proposal for generating a quasiperiodic signal by superimposing two DTCs rather than a new interaction-enabled phase. If the peaks vanish, the current interpretation would be supported. The conditional verdict is therefore appropriate and should remain until the control is performed. I do not see a separate, more severe internal inconsistency that would justify rejection on the current evidence; the missing baseline is specific and addressable.","tokens_in":15381,"tokens_out":7128,"duration_ms":73468,"concrete_test":"Run the same Floquet evolution with the cross-boundary coupling removed, i.e., in Eq. (1) set P_{N_L+1}=1 in the term i=N_L and P_{N_L}=1 in the term i=N_L+1, or equivalently evolve two independent chains of lengths N_L and N_R with the same Ω_L, Ω_R, T_L, T_R, θ=π, and Z2 initial product state. Compute m(t) from Eq. (3) and F(t) from Eq. (6) over the same total time (e.g., t=1000 T_L) and compare the Fourier amplitudes at (1/2,1/2) and (-1/2,1/2) with the coupled case. If the peaks survive with comparable amplitude, the paper's causal claim in the abstract and Sec. VII is falsified; if they vanish or are strongly suppressed, the interaction attribution is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in the abstract and Sec. VII is that the aperiodic response is 'indeed caused by interaction between systems via Rydberg blockade effect.' This is load-bearing because the paper's novelty is precisely that coupling two DTCs creates a new DTQC phase, not just a trivial combination of independent DTCs. The evidence does not yet rule out the trivial combination. The order parameter m(t)=|Σ_i (-1)^i n_i|/N (Eq. 3) is a non-additive function of the separate left and right staggered magnetizations M_L and M_R. In the decoupled limit each subsystem is expected to oscillate subharmonically at f_L/2 and f_R/2. Then |M_L(t)+M_R(t)| contains intermodulation products: for M_L≈A cos(π f_L t) and M_R≈B cos(π f_R t), the rectified sum has Fourier components at (f_L+f_R)/2 and |f_L-f_R|/2, exactly the claimed (1/2,1/2) and (-1/2,1/2) peaks. The fidelity F=|⟨Z2|ψ(t)⟩| (Eq. 6) has the same weakness: for a product state it factorizes into the product of the two subsystem fidelities, and a product of two incommensurate subharmonic signals generates sum and difference frequencies without any coupling. No decoupled control is reported in which the boundary PXP projectors between sites N_L and N_L+1 are switched off. The statement in Sec. VI that nonzero entanglement entropy demonstrates coupling is insufficient: it shows that some boundary entanglement exists, but does not show that the dominant Fourier peaks of m and F are caused by that coupling. The high-frequency decoupling limit in Fig. 2(d) is not a controlled decoupling, since it changes the drive parameters rather than isolating the interaction term.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a scheme to generate a discrete time quasi-crystal (DTQC) in a one-dimensional Rydberg atom array by dividing the chain into two halves and driving them with two Floquet frequencies whose ratio is the golden ratio. The authors argue that the Rydberg blockade interaction at the boundary couples the two halves, each of which individually supports a discrete time crystal (DTC), and that this coupling produces a stable aperiodic response. They study the antiferromagnetic order parameter m defined in Eq. (3), the fidelity F defined in Eq. (6), compute Fourier spectra that show peaks at half-integer combinations of the two drive frequencies, map phase diagrams using a lifetime criterion tau > 30 T_L, and compute the bipartite entanglement entropy between the two halves. The central claim, stated in the abstract and in Section VII, is that the aperiodic response is 'indeed caused by interaction between systems via Rydberg blockade effect.'","tokens_in":15824,"tokens_out":14338,"duration_ms":129192,"significance":"If the interaction-causality claim were established, the work would be a notable proposal that extends DTC physics to quasi-periodic time order on a Rydberg array, which is a leading quantum simulation platform. The numerics are carefully specified (exact diagonalization for N<14, MPS time evolution with truncation error below 5e-9 and bond dimension up to 500), the drive ratio is fixed a priori to the golden ratio rather than fitted, and the DTC benchmark is grounded in Refs. 31 and 32. The paper also makes concrete falsifiable predictions about the location of Fourier peaks. However, as written, the central claim is not yet supported because the chosen order parameter and fidelity are nonlinear or factorizing functions of the two subsystem observables in a way that produces the claimed sum- and difference-frequency peaks even when the two subsystems are completely decoupled. The missing decoupled control is a load-bearing gap, but it is addressable within the scope of the manuscript.","major_comments":[{"comment":"The central claim that the aperiodic response is 'indeed caused by interaction between systems via Rydberg blockade effect' is not supported by the evidence, because the order parameter m = |Σ_i (-1)^i n_i|/N is a nonlinear function of the separate staggered magnetizations M_L and M_R. In the decoupled limit, each subsystem is expected to oscillate subharmonically at ω_L/2 and ω_R/2. For A=B, |M_L+M_R| = 2|A cos((ω_L+ω_R)t/4) cos((ω_L-ω_R)t/4)|; each |cos| factor has a DC component and harmonics, so the product contains spectral lines at (ω_L+ω_R)/2 and |ω_L-ω_R|/2, exactly the (1/2,1/2) and (-1/2,1/2) peaks shown in Fig. 1(c). The manuscript does not report a control simulation in which the boundary PXP projectors between sites N_L and N_L+1 are switched off, so the observed peaks could arise from two independent DTCs even with zero coupling. Please add such a decoupled control (e.g., remove the terms at i=N_L and i=N_L+1 in Eq. (1)) and compare the resulting Fourier spectra with the coupled case; if the peaks persist, the causal attribution in the abstract and Section VII must be revised.","section":"Section III, Eq. (3); Abstract and Section VII"},{"comment":"The fidelity F = |⟨Z2|ψ(t)⟩| has the same weakness as m. In a decoupled evolution with |ψ(t)⟩ = |ψ_L(t)⟩ ⊗ |ψ_R(t)⟩, F factorizes as F_L(t) F_R(t) (up to a constant), and the product of two subharmonic signals at angular frequencies ω_L/2 and ω_R/2 contains components at (ω_L+ω_R)/2 and |ω_L-ω_R|/2. The fidelity peaks labeled (1/2,1/2) and (-1/2,1/2) in Fig. 3(b) are therefore fully consistent with a zero-coupling product state. The paper should either provide a decoupled-control fidelity spectrum or state explicitly that the fidelity data do not distinguish coupled from uncoupled dynamics.","section":"Section V, Eq. (6), Fig. 3"},{"comment":"The argument that a non-zero entanglement entropy between the left and right halves 'demonstrates that there exists the coupling between two DTCs' only shows that some boundary coupling exists; it does not establish that the dominant Fourier peaks of m and F are caused by that coupling. To support the causal claim, the paper should connect coupling strength to the spectral response, for example by tuning the boundary PXP coupling (or comparing with the decoupled control suggested above) and showing that the (1/2,1/2) and (-1/2,1/2) peak amplitudes vanish or change qualitatively as the coupling is removed. The manuscript's own description of the large-frequency limit in Fig. 2(d) as 'simple additivity of two independent DTCs' further illustrates that the global observables can produce the same spectral structure without interaction.","section":"Section VI, first two paragraphs"}],"minor_comments":[{"comment":"The sentence beginning 'While we analysis its robustness' contains a grammatical error and should read 'While we analyze its robustness'; the phrase 'We significantly calculate the entanglement entropy' is awkward and should be rephrased.","section":"Abstract"},{"comment":"The text refers to Fig. 7 for the initial-state dependence and to Fig. 6 for the fidelity phase diagrams, but the figure captions describe Fig. 6 as initial-state dependence and Fig. 7 as fidelity phase diagrams; the references or the captions should be swapped so that they agree.","section":"Appendices B and C"},{"comment":"The caption parameters do not match the panel assignment in the text: the text assigns (a) uniform quench, (b) non-uniform quench, (c) DTC, and (d) DTQC, while the caption lists parameter sets that place the DTC in (b) and the non-uniform quench in (c). Please reconcile the caption with the text.","section":"Appendix A, Fig. 5 caption"},{"comment":"Equation (5) uses f both as a summation index and as the argument of cos(ft), which is confusing; using a dummy index such as j would avoid the conflict.","section":"Section IV, around Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The missing decoupled control is the key issue and is fixable: the authors can simulate the Hamiltonian with the boundary PXP terms removed and show whether the (1/2,1/2) and (-1/2,1/2) peaks persist. If they persist, the manuscript should be reframed as a study of a driven two-part PXP system whose global observables exhibit quasiperiodic peaks, without claiming that interaction causes those peaks. If they disappear, the authors would need to explain why the existing entanglement argument is insufficient and provide the new evidence. The paper is within scope for a quantum-information journal and the numerical methods are solid, so with the control and careful revision of the causal language the work could become publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a decent numerical proposal for realizing a discrete time quasicrystal in a Rydberg PXP chain by driving two halves at incommensurate rates. The numerics are clearly described and the observed spectral features are real, but the paper's central causal claim is not established. The order parameter m and fidelity F are nonlinear enough that two independent DTCs can produce the exact sum and difference peaks without any coupling. The authors need to run the decoupled control.\n\nWhat is genuinely new is the specific two-coupled-DTC construction on a Rydberg chain — DTQC theory exists (Refs. 84–86) and the NV experiment exists (Ref. 25), but not this platform or the coupled-subchain construction. The phase diagrams, size dependence, and the explicit TeNPy/MPS parameters are useful. The paper is also honest about parameter regimes and includes initial-state dependence in the appendices. That part is solid.\n\nThe main flaw: there is no simulation with the boundary blockade coupling switched off. If the PXP projectors between sites N_L and N_L+1 are removed, each half should evolve independently and show subharmonic response at f_L/2 and f_R/2. For m = |M_L + M_R|/N, with M_L and M_R oscillating at those half-frequencies, the absolute value produces intermodulation at (f_L+f_R)/2 and |f_L−f_R|/2 — exactly the claimed (1/2,1/2) and (−1/2,1/2) peaks. Fidelity F = |⟨Z2|ψ⟩| has the same weakness: for a product state it factorizes into subsystem fidelities, and the product of two incommensurate subharmonic signals also generates sum and difference frequencies. Nonzero entanglement entropy only shows there is some boundary coupling; it does not show that the dominant Fourier peaks of m and F arise from that coupling. Fig. 2(d)'s high-frequency decoupling limit changes the drive parameters rather than isolating the interaction, so it is not a controlled control. The Eq. (4) criterion is partially self-referential, but since the drive parameters are not fitted to the peaks, I'd call that a minor issue.\n\nThis warrants peer review, because the missing control is specific and addressable. If the authors add the decoupled baseline and the peaks persist — or differ from the independent-subsystem prediction — the result becomes solid and interesting. If not, the abstract's \"indeed caused by interaction\" claim should be withdrawn. I'd send it out with a clear request for that baseline. I'd cite it as a proposed construction, not yet as a demonstration.","headline":"A concrete Rydberg proposal for DTQC with clean numerics, but the central claim that the observed intermodulation peaks require Rydberg coupling is not yet supported: the nonlinear order parameter can fake them.","tokens_in":16326,"tokens_out":3209,"would_cite":true,"duration_ms":33507,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two golden-ratio-driven time crystals can form a discrete time quasi-crystal in a Rydberg chain.","keywords":["discrete time quasi-crystal","Rydberg atom chain","PXP model","quantum many-body scars","Floquet driving","golden ratio","time-translation symmetry","entanglement entropy"],"falsifier":"Remove the Rydberg blockade constraint on the single bond connecting the two halves (equivalently, set the boundary PXP projector to zero) and repeat the golden-ratio double-Floquet protocol. If the Fourier peaks at $(1/2,1/2)$ and $(-1/2,1/2)$ persist with comparable amplitude while the bipartite entanglement entropy drops to zero, the aperiodic response would be shown to come from two independent time crystals rather than from their interaction.","tokens_in":15189,"feed_emoji":"⏳","tokens_out":6474,"duration_ms":65627,"temperature":0.7,"pith_summary":"Discrete time quasi-crystals extend discrete time crystals by replacing periodic time order with quasi-periodic order. The paper proposes a concrete way to create one in a Rydberg atomic chain: split the chain into two halves, drive each half as its own discrete time crystal, and choose the two drive periods in the golden ratio, the most incommensurate possible ratio. It claims the two subharmonic responses mix through the Rydberg blockade interaction at the boundary, producing stable aperiodic oscillations in the antiferromagnetic order parameter and fidelity, with sharp Fourier peaks at the sum and difference of half the two drive frequencies. If correct, this gives Rydberg arrays a protocol for a new non-equilibrium phase and connects time-crystal physics to quasiperiodic order.","feed_headline":"Two golden-ratio time crystals form one time quasi-crystal","feed_subtitle":"Rydberg-chain simulation shows stable aperiodic oscillations whose Fourier peaks mix the two drive frequencies.","key_machinery":"The machinery is the bipartite Floquet-driven PXP model: the Rydberg blockade is encoded by projectors $\\hat P_{i-1}\\hat X_i\\hat P_{i+1}$, so no two neighbouring atoms can be simultaneously excited, and the drive is a sequence of two $\\delta$-function kicks with periods $T_L$ and $T_R$ satisfying $T_L/T_R=(\\sqrt5+1)/2$. The paper couples this with the observable $m(t)=|\\sum_i(-1)^i\\hat n_i|/N$, whose Fourier peaks are indexed by half-integer combinations $f=(k_1/2)f_L+(k_2/2)f_R$. The golden-ratio choice makes the two drive frequencies maximally incommensurate, so the sum-frequency peaks are genuine quasiperiodic responses rather than harmonics of a common period; the PXP constraint at the boundary is what lets the two halves talk to each other. The antiferromagnetic initial state and the $\\theta=\\pi$ Floquet strength, where the drive anti-commutes with the subsystem Hamiltonian, stabilize the scarred revivals on each half.","core_discovery":"The central claim is that coupling two discrete time crystals with maximally incommensurate driving frequencies produces a discrete time quasi-crystal phase. Concretely, the paper studies a one-dimensional Rydberg chain described by the PXP model, divided into left and right subsystems with Rabi frequencies and Floquet pulse periods whose ratio is the golden ratio $(\\sqrt5+1)/2$. Starting from the experimentally natural $Z_2$ antiferromagnetic state, the antiferromagnetic order parameter $m=|\\sum_i(-1)^i\\hat n_i|/N$ keeps oscillating without an identifiable period; its Fourier spectrum shows stable subharmonic peaks at $(1/2,1/2)$ and $(-1/2,1/2)$ in units of the two driving frequencies, exactly the sums and differences of the two half-frequencies. The paper interprets these peaks as the signature of the quasi-crystal phase and argues that the aperiodic response is caused by the Rydberg blockade coupling between the two subsystems, not merely by adding two independent time crystals. Fidelity and bipartite entanglement entropy simulations are used to support the phase's existence, robustness, and boundary-localized entanglement.","pith_inferences":["A decoupled-control simulation, with the boundary blockade term switched off, would settle whether the $(1/2,1/2)$ peaks require interaction; the absolute value in $m$ means intermodulation of two independent half-frequency oscillators could in principle mimic them.","The golden-ratio construction suggests a general recipe: any pair of irrational drive frequencies with strong incommensurability might produce analogous quasiperiodic time order, so scanning the frequency ratio would test how special the golden ratio is.","Because the entanglement entropy is boundary-localized and small, the quasi-crystal phase may be fragile to strengthening the boundary coupling or adding bulk interactions; varying the boundary spacing or Rabi frequency would probe this fragility.","The observed exponential decay of fidelity peaks with system size implies that direct revival diagnostics become impractical in large chains, so the phase is best certified through the order-parameter spectrum or boundary entanglement rather than wavefunction overlap."],"forward_implications":["At intermediate driving frequencies, the Fourier spectrum of the antiferromagnetic order parameter is dominated by the peaks $(1/2,1/2)$ and $(-1/2,1/2)$, the sum and difference of half the two drive frequencies.","Low driving frequency pushes the response into a chaotic regime, while high driving frequency decouples the two halves and leaves only the simple addition of two independent discrete time crystals.","The quasi-crystal order parameter is robust to system size, with the typical Fourier amplitudes saturating as the chain grows, whereas fidelity amplitudes decay exponentially with system size.","The bipartite entanglement entropy stays low and appears concentrated at the boundary, indicating that the coupling between the two time crystals comes from the Rydberg blockade of the nearest boundary sites.","The phase region widens near modulation strength $\\theta=\\pi$ and a driving frequency about twice the inherent scar oscillation frequency, giving concrete parameter targets for an experiment."],"supporting_citations":[{"why":"Experimental Floquet control of many-body dynamics in Rydberg atom arrays, the platform and scar-control setting the proposal builds on.","marker":"[31]"},{"why":"Supplies the discrete time crystal protocol and the delta-function Floquet modulation scheme that the two subsystems adapt.","marker":"[32]"},{"why":"Prior experimental realization of discrete time quasicrystals in a different platform, defining the target phenomenon this paper moves to Rydberg arrays.","marker":"[25]"},{"why":"Observation of quantum many-body scars and $Z_2$ revivals in a Rydberg chain, the dynamical basis for the scarred time crystal.","marker":"[56]"},{"why":"Introduces the PXP model used throughout to encode the Rydberg blockade constraint.","marker":"[87]"},{"why":"Provides the scarred eigenstate picture and the linear relation between the inherent oscillation frequency and the Rabi frequency.","marker":"[67]"},{"why":"Establishes weak ergodicity breaking from quantum many-body scars, justifying the choice of drive periods proportional to Rabi frequencies.","marker":"[68]"}],"fun_headline_variants":["Golden-ratio drives yield time quasi-crystal in Rydberg chain","Rydberg chain couples two time crystals to form quasi-crystal","Aperiodic order from gold-ratio coupling of time crystals","Two incommensurate drives create Rydberg time quasi-crystal","Time quasi-crystal emerges from golden-ratio Rydberg coupling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The causal claim that the aperiodic response arises from the Rydberg-blockade interaction between the two halves rests on the assumption that the absolute-value order parameter $m(t)$ reports genuine coupled dynamics; if two decoupled subharmonic oscillators already produce the same Fourier peaks, the claimed mechanism is not established.","fun_headline_variants_meta":{"raw":{"variants":["Golden-ratio drives yield time quasi-crystal in Rydberg chain","Rydberg chain couples two time crystals to form quasi-crystal","Aperiodic order from gold-ratio coupling of time crystals","Two incommensurate drives create Rydberg time quasi-crystal","Time quasi-crystal emerges from golden-ratio Rydberg coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000244,"raw_usage":{"total_tokens":1539,"prompt_tokens":961,"completion_tokens":578,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":486}},"tokens_in":577,"tokens_out":578,"duration_ms":5156,"temperature":1.0,"reasoning_tokens":486,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:39:18.636619+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Remove the Rydberg blockade constraint on the single bond connecting the two halves (equivalently, set the boundary PXP projector to zero) and repeat the golden-ratio double-Floquet protocol. If the Fourier peaks at $(1/2,1/2)$ and $(-1/2,1/2)$ persist with comparable amplitude while the bipartite entanglement entropy drops to zero, the aperiodic response would be shown to come from two independent time crystals rather than from their interaction.","supporting_citations":[{"cited_title":"Quantum Lotka-Volterra dynamics","cited_arxiv_id":"2408.01726","evidence_quote":"Experimental Floquet control of many-body dynamics in Rydberg atom arrays, the platform and scar-control setting the proposal builds on."},{"cited_title":"Bluvstein, A","cited_arxiv_id":null,"evidence_quote":"Supplies the discrete time crystal protocol and the delta-function Floquet modulation scheme that the two subsystems adapt."},{"cited_title":"Randall, C","cited_arxiv_id":null,"evidence_quote":"Prior experimental realization of discrete time quasicrystals in a different platform, defining the target phenomenon this paper moves to Rydberg arrays."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the PXP model used throughout to encode the Rydberg blockade constraint."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes weak ergodicity breaking from quantum many-body scars, justifying the choice of drive periods proportional to Rabi frequencies."}],"review_version":1}