{"id":"b426ec11-17aa-4ff3-86f6-d5734f8ebd90","arxiv_id":"2512.16193","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Delayed pairwise Kuramoto coupling expands into an effective delay-free model with three-body (higher-order) interactions that reproduces delay-induced bistability and synchronization transitions.","lead":"Time-delayed pairwise coupling in Kuramoto oscillators is approximately equivalent, to second order in the coupling, to a delay-free model with three-body interactions. This gives an analytically tractable explanation of delay-induced bistability and a stability diagram that matches the original delayed system.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed regime ε≪1 and τ<1/ε is too broad: the derivation requires ετ≪1, and several numerical comparisons use ετ close to or above 1, so the central 'quantitatively reproduces' assertion is not established.","rationale":"The reader's weakest assumption identifies exactly the small-ετ requirement as load-bearing. My analysis sharpens this: the paper's stated condition τ<1/ε is not sufficient, and the numerical figures themselves include ετ values that are not small. This reinforces the reader's CONDITIONAL verdict rather than overturning it: the central derivation is plausible, but the claimed quantitative agreement is not established over the stated parameter regime. The proposed numerical test would settle whether the approximation degrades as ετ→1, which is the key unresolved point.","tokens_in":12527,"tokens_out":7243,"duration_ms":72304,"concrete_test":"Fix ε=0.05 (or 0.02) and simulate both Eq. (1) and Eq. (6) (or solve Eq. (11)) for τ chosen so that ετ = 0.1, 0.3, 0.5, 0.7, 0.9, and 1.1. Compute the order parameter R in the bistable region and the stability boundaries in the (ε,τ) plane, and measure the maximum deviation |R_delay − R_reduced| over the same initial conditions. If the error grows sharply as ετ approaches 1, the claimed τ<1/ε regime is invalid and the quantitative-reproduction claim must be restricted to ετ≪1; if the error stays small for all ετ<1, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reduction's validity condition is not simply τ<1/ε. In Appendix A1, the delayed phase is replaced via θ_k(t−τ) ≈ θ_k(t)−ω_kτ + O(ε), where the omitted term is actually O(ετ) when integrated over the delay. Substituting this into the O(ε) coupling produces a first neglected contribution of order ε·ετ = ε^2τ·ε, i.e. O(ε^3 τ). For the retained O(ε^2 τ) three-body term to be the leading correction, one needs ετ≪1, not merely ετ<1. The paper's stated regime τ<1/ε allows ετ arbitrarily close to 1; at ε=0.05, τ=19 gives ετ=0.95, where the neglected term is comparable to the retained three-body term. The numerical comparisons also enter this regime: Fig. 2(a) uses τ=2.8 with ε up to 0.27, so ετ reaches about 0.76, and Fig. 2(b) uses ε=0.3 with τ extending beyond 3.3, giving ετ>1. Thus the stability-diagram comparison does not cleanly test the small-ετ regime. Consequently, the statement that the higher-order model 'quantitatively reproduces' the delayed model within τ<1/ε is not supported by the present evidence; the claim can be maintained only for ετ sufficiently small, and the paper's own regime statement is internally inconsistent with the requirement it states one sentence earlier.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies globally coupled Kuramoto oscillators with a uniform time delay τ in the pairwise coupling, Eq. (1). Under a weak-coupling assumption, it Taylor-expands the delayed phase and substitutes the equations of motion repeatedly to obtain a delay-free model, Eq. (6), consisting of a Kuramoto–Sakaguchi pairwise term and a (2,−1,−1) three-body harmonic whose strength is ε²τ. The paper demonstrates bistability for N=2, compares finite-N simulations of the delayed and reduced models, applies the Ott–Antonsen ansatz to derive the order-parameter dynamics, Eqs. (10)–(12), and derives stability conditions, Eqs. (14)–(15), which it benchmarks against the Yeung–Strogatz result. The central assertion is that, within the regime ε≪1 and τ<1/ε, the higher-order Kuramoto model quantitatively reproduces the synchronization transitions and bistability of the time-delayed model.","tokens_in":12899,"tokens_out":21253,"duration_ms":198507,"significance":"The conceptual message — that delayed pairwise interactions can be recast, through a controlled weak-coupling expansion, as higher-order interactions — is interesting and potentially useful, connecting two active fields. The derivation is explicit and self-contained, with no fitted parameters; the reduction to a low-dimensional Ott–Antonsen equation makes the reduced model analytically tractable; and the stability boundaries are checked against the independent Yeung–Strogatz diagram. These are genuine strengths. The main weakness is that the stated domain of validity is broader than the expansion actually supports: the accuracy is controlled by ετ, not by ε alone, and the numerical comparisons include parameter values with ετ close to or above 1.","major_comments":[{"comment":"The derivation is an expansion in ετ, not in ε alone. In Appendix A1, θ_k(t−τ) is replaced by θ_k(t)−ω_kτ plus a remainder of order ετ. When this remainder is substituted into the O(ε) coupling, the first omitted correction to the retained O(ε²τ) three-body term is smaller by a factor O(ετ). The paper itself states that the accuracy is characterized by the smallness of ετ, but then asserts the regime τ<1/ε, which allows ετ arbitrarily close to 1. At ετ≈1 the neglected corrections are comparable to the retained three-body term, so the statement that Eq. (6) approximates Eq. (1) to O(ε²) is not valid on the claimed domain. The validity condition should be ετ≪1, and the phrase τ<1/ε in the main text and in Fig. 3 should be replaced or supplemented.","section":"Appendix A1 and Eq. (6)"},{"comment":"The numerical evidence for the central 'quantitatively reproduces' claim lies largely outside the small-ετ regime. In Fig. 2(a), τ=2.8 and the reported bifurcations occur at ε=0.12 and ε=0.27, giving ετ≈0.34 and 0.76; if the plotted ε-range extends beyond 0.3, ετ exceeds 1. In Fig. 2(b), ε=0.3, so ετ>1 for any τ>3.3, which covers most of the displayed periodic sequence. Fig. 3 compares stability boundaries over the entire region below the line τ=1/ε, including ετ close to 1. These comparisons therefore do not cleanly test the asymptotic regime in which the reduction is derived. The authors should either restrict the numerical and analytical comparisons to ετ≪1, or provide an explicit error estimate showing that the agreement persists when ετ is not small.","section":"Fig. 2 and Fig. 3"},{"comment":"The reduction to the identical-oscillator limit takes γ→0 after applying the Ott–Antonsen ansatz. This is a singular limit, and the replacement of f1 and f2 by τ and 0 uses (ω_l−ω_k)τ=O(ετ), which again requires ετ≪1. The paper does not justify that the N→∞ and γ→0 limits commute, nor does it discuss the fact that the Lorentzian becomes a δ-distribution. Since the stability diagram is a central quantitative output, the authors should either derive the identical-oscillator case directly or verify the γ→0 limit against finite-N simulations of Eq. (6) in the ετ≪1 regime.","section":"Appendix A3 and Eqs. (10)–(15)"}],"minor_comments":[{"comment":"The abstract says the reduced model exhibits 'qualitatively consistent' transitions, while the main text and conclusion say it 'quantitatively reproduces' the phenomena. These statements should be aligned to the same strength, especially after the validity regime is corrected.","section":"Abstract and §5"},{"comment":"The sums over k and l allow l=j, so the 'three-body' term includes O(1/N) pairwise and constant contributions. In the thermodynamic limit they vanish, but the paper should either exclude l=j explicitly or state that these contributions are negligible for large N.","section":"Eq. (6)"},{"comment":"The phrase 'low-dimensional mode' should read 'low-dimensional model'. Also, the redefinition of z1 in the γ→0 limit is stated but the notation could be cleaner.","section":"Eqs. (10) and (11)"},{"comment":"Please specify the full parameter ranges, the number of oscillators, initial conditions, and how the stable and unstable branches were obtained (e.g., forward/backward sweeps). It would also be helpful to mark the ετ≪1 validity boundary or indicate the ετ values at the reported transition points.","section":"Fig. 2 caption"},{"comment":"The sentence 'The case with ω_l=ω_k is trivial' should be expanded: the limiting values f1=τ and f2=0 are used in the central Eq. (6), so they deserve an explicit derivation.","section":"Appendix A1"}],"recommendation":"major_revision","confidential_remarks":"The core reduction is likely correct in the ετ≪1 regime, and the conceptual finding is worth publishing after revision. The main task for the authors is to narrow the claimed validity domain, redo or re-express the numerical and stability comparisons so that they actually test the small-ετ regime, and justify the γ→0 limiting procedure. I do not see a need for rejection, but the current overly broad 'quantitative' claim is not supportable as written."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: the paper genuinely shows that, to O(ε²), the delayed Kuramoto model reduces to a delay-free model with a Sakaguchi pairwise term plus a (2,−1,−1) three-body interaction whose strength scales as ε²τ. The derivation in the appendix is explicit and the OA order-parameter equation yields stability conditions that reproduce the independent Yeung–Strogatz diagram. That is a real step, and the paper deserves a serious referee.\n\nThe main soft spot is the claimed validity regime. The text says the approximation is controlled by the smallness of ετ, but then states ε≪1 and τ<1/ε. These are not the same: τ<1/ε allows ετ close to 1. At ετ~1, the omitted O(ε³τ²) terms are comparable to the retained O(ε²τ) term (their ratio is ετ). So the \"quantitatively reproduces\" statement holds only for ετ≪1. The numerics in Fig. 2 do not respect that: with τ=2.8 and ε up to about 0.27 you get ετ≈0.76, and with ε=0.3 and τ beyond 3.3 you get ετ>1. The stability-diagram comparison in Fig. 3 also enters that range. This is an overstatement rather than an algebraic error; the derivation itself is fine.\n\nTwo smaller points. The paper does not clearly discuss the relation to the higher-order phase reduction for delay-coupled oscillators in [52,53], one of which has a co-author in common; that needs explicit positioning. And the \"quantitative\" claim in the abstract is stronger than the visual comparisons support.\n\nWho gets value: people working on higher-order Kuramoto models, delay-coupled oscillator networks, or phase reduction. The OA stability analysis is directly usable. I would cite it if I worked in that intersection. It deserves peer review; the main request to the authors should be to state the validity condition as ετ≪1, re-run the numerics in a genuinely small-ετ window (or explain why agreement persists beyond it), and tighten the discussion of [52,53].","headline":"A clean O(ε²) reduction from delayed pairwise to three-body Kuramoto coupling, but the claimed τ<1/ε regime is too broad; the numerics run where ετ is not small.","tokens_in":13394,"tokens_out":9006,"would_cite":true,"duration_ms":80377,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34C15","34K18","34K20","37N25"],"pacs":["05.45.Xt","89.75.-k"],"model":"deepseek-v4-flash","headline":"Delayed pairwise coupling in Kuramoto oscillators is equivalent to three-body coupling to second order, and the three-body term explains delay-induced bistability.","keywords":["higher-order interactions","Kuramoto model","time delay","three-body coupling","order-parameter reduction","bistability","synchronization transition","phase oscillators"],"falsifier":"Simulate the original delayed model and the reduced higher-order model for a parameter point where epsilon*tau is not small, e.g., epsilon = 0.5 and tau = 4 (product 2, violating tau < 1/epsilon), and compare synchronization transition curves: if the reduced model fails to reproduce bistable regions or transition thresholds, the claimed equivalence is limited to the small-parameter regime.","tokens_in":12391,"feed_emoji":"🕰️","tokens_out":4952,"duration_ms":51577,"temperature":0.7,"pith_summary":"The paper argues that time-delayed pairwise interactions in Kuramoto oscillators are not merely a phase lag: when expanded to second order in the coupling strength, the delayed system becomes a delay-free higher-order Kuramoto model with both pairwise phase-lagged coupling and a three-body (2,-1,-1) interaction whose strength grows as the product of coupling and delay. This three-body term is what produces the bistability between incoherence and full synchrony seen in the original delayed model. If true, delays provide a concrete mechanistic origin for higher-order interactions, and delay-induced phenomena become tractable with standard oscillator-analysis tools.","feed_headline":"Delayed pairwise coupling is really three-body coupling","feed_subtitle":"Second-order expansion shows the delay term creates a three-body interaction that explains bistability of synchronization.","key_machinery":"The central device is the second-order Taylor expansion of the delayed sine coupling in powers of the coupling strength, equivalently in powers of the product epsilon*tau. Iterated substitution of the equations of motion converts delayed derivatives into combinations of phase differences, and reorganizing terms by powers of epsilon yields the delay-free higher-order model. The Ott-Antonsen ansatz, which assumes the phase distribution's Fourier coefficients decay geometrically, then collapses the infinite-dimensional system to a closed equation for the synchronization degree R, from which the stability conditions for incoherent and fully synchronized states are read off.","core_discovery":"For identical oscillators with small coupling epsilon and delay tau, the time-delayed Kuramoto model reduces to an ordinary differential equation: each pair experiences sine coupling with phase lag omega0 tau, plus a three-body (2,-1,-1) harmonic with amplitude proportional to epsilon^2 tau. The derivation expands the delayed argument in a Taylor series and repeatedly substitutes the equations of motion to eliminate delay derivatives, keeping terms through O(epsilon^2). Numerical simulations and the standard Ott-Antonsen order-parameter reduction show that this reduced model reproduces the synchronization transitions and, crucially, the bistable region of the original delayed system, provide","pith_inferences":["Because the three-body strength scales as epsilon^2 tau, delay-induced effects grow linearly with delay but quadratically with coupling; experiments or numerics that vary epsilon and tau separately could test this scaling directly.","The same Taylor-expansion route should generate four-body interactions at O(epsilon^3); these may become visible for delays near tau ~ 1/epsilon and could produce cluster states or other phenomena absent from the truncated model.","For non-identical frequencies the effective couplings acquire frequency-dependent phase lags, suggesting a frequency-filtered interaction that could matter for broad frequency distributions.","The equivalence between delay and higher-order interactions may extend to non-global network topologies, where the three-body term would couple triangles of oscillators and make delayed network dynamics amenable to hypergraph analyses."],"forward_implications":["Bistability between incoherence and full synchrony in the delayed Kuramoto model is caused by the three-body term; the pairwise phase-lagged part alone cannot produce hysteresis.","The effective higher-order model is a system of ordinary differential equations, so all standard analysis and control methods for higher-order Kuramoto networks become applicable to time-delayed systems.","Stability boundaries of the incoherent and synchronized states reduce to simple trigonometric inequalities in epsilon and tau, which match the independently known diagram for the delayed model in the small-parameter regime.","The reduction offers a systematic route to derive many-body effective couplings for other delay-coupled oscillator systems beyond Kuramoto dynamics."],"fun_headline_variants":["Time delay hides three-body interactions in sync","Delay-induced three-body terms explain bistability","Pairwise delay reduces to triple coupling","How time delay creates three-body coupling"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The approximation is valid only when the product of coupling strength and delay, epsilon*tau, is small enough that truncating the Taylor expansion at second order is safe; the paper assumes tau < 1/epsilon, and if this fails the residual delay and four-body terms can change the dynamics.","fun_headline_variants_meta":{"raw":{"variants":["Time delay hides three-body interactions in sync","Delay-induced three-body terms explain bistability","Pairwise delay reduces to triple coupling","How time delay creates three-body coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000505,"raw_usage":{"total_tokens":2241,"prompt_tokens":622,"completion_tokens":1619,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":366,"completion_tokens_details":{"reasoning_tokens":1566}},"tokens_in":366,"tokens_out":1619,"duration_ms":12197,"temperature":1.0,"reasoning_tokens":1566,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T15:36:32.648004+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the original delayed model and the reduced higher-order model for a parameter point where epsilon*tau is not small, e.g., epsilon = 0.5 and tau = 4 (product 2, violating tau < 1/epsilon), and compare synchronization transition curves: if the reduced model fails to reproduce bistable regions or transition thresholds, the claimed equivalence is limited to the small-parameter regime.","supporting_citations":[],"review_version":1}