Pith. sign in

REVIEW 2 major objections 4 minor 30 references

The standard 1D swarmalator model is a normal form of chemotactic oscillators, with all four collective states meeting at one cusp.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-14 15:23 UTC pith:YWSXA77P

load-bearing objection Solid normal-form grounding plus the first real stability map for the 1D swarmalator, with an exact phase-wave boundary and a needed correction to the prior sync order parameter. the 2 major comments →

arxiv 2607.09810 v1 pith:YWSXA77P submitted 2026-07-10 nlin.AO math.DS

A solvable normal form for coupled swarmalators

Kevin P. O'Keeffe This is my paper
classification nlin.AO math.DS
keywords swarmalatorsnormal formchemotactic oscillatorsphase wavessynchronizationbifurcation cusporder parametersKuramoto model
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Swarmalators combine synchronization in phase with self-assembly in space, but their simplest models have been treated as ad hoc toys whose stability boundaries were largely unsolved. This paper shows that the standard one-dimensional ring model is the first-harmonic, zero-lag reduction of adapted chemotactic limit-cycle oscillators, so its collective behavior is generic rather than special. It then derives the stability boundaries of the four known states—async, phase wave, mixed, and sync—showing that they meet at a single cusp in coupling space. Along the way it corrects an earlier order-parameter formula for the fully synchronized state and finds a non-monotonic rise-and-dip of coherence that has no counterpart in the ordinary Kuramoto model. The result supplies the same three pillars that made Kuramoto central: a normal-form origin, solvable mean-field theory, and a concrete experimental candidate in mobile Belousov–Zhabotinsky droplets on a ring.

Core claim

The canonical one-dimensional swarmalator equations are recovered exactly as the first-harmonic, zero-lag normal form of adapted chemotactic Hopf oscillators, and the stability boundaries of their four collective states meet at a single organizing cusp at (K, J) = (4, 2). The phase-wave boundary is given by an exact elementary self-consistency F(K, J) = 0 that is a finite sum of arctangents; the sync order parameter must be decomposed into locked plus tongue contributions that correct the earlier Ott–Antonsen formula for all nonzero J.

What carries the argument

The continuum linearization of the continuity equation about each partially locked density, reduced to a finite-rank susceptibility self-consistency. For the phase wave this collapses to the exact elementary condition F(K, J) = χ_locked + χ_drift − 1 = 0; for sync it becomes an exchange-odd scalar condition H(K, J) = 0 whose asymptotics and locking diamond organize the remaining boundaries.

Load-bearing premise

The chemical field is assumed both to adapt to contrast rather than absolute level and to be dominated by its first spatial harmonic with zero lag; if either fails, the claim that the toy model is a generic normal form collapses.

What would settle it

Direct N ≫ 1 simulations (or a microfluidic ring of mobile BZ droplets) that track whether the four-state cascade and the measured phase-wave and sync boundaries coincide with the predicted curves F(K, J) = 0 and the large-J asymptote K3 = J + (2/π) log(2J) + 4/π.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The manuscript argues that the canonical 1D swarmalator model (Eqs. 1–2) is the first-harmonic, zero-lag normal form of adapted chemotactic Hopf oscillators (building on Tanaka), so its collective behavior is generic rather than ad hoc. In sum/difference coordinates it then derives the stability boundaries of the four known states (async, phase wave, mixed, sync), shows they meet at a single cusp (K,J)=(4,2), supplies an exact elementary self-consistency F(K,J)=0 (finite arctangent sums) for the phase-wave boundary, corrects the prior Ott–Antonsen sync order parameter by a locked-plus-tongue decomposition that matches N=10^6 simulations, and reports a non-monotonic sync response for J>Jc together with a parameter-free large-J sync boundary.

Significance. If the results hold, the paper supplies the missing normal-form foundation and the first nontrivial stability/bifurcation theory for the main theoretical testbed of swarmalator dynamics. The phase-wave boundary is reduced to a closed elementary condition with matched asymptotics; the sync order parameter is rebuilt from locked and tongue contributions and demonstrably corrects an ~40% OA error; and the large-J sync boundary K3=J+(2/π)log(2J)+4/π is parameter-free and matches simulation. These are concrete, falsifiable advances that place the 1D model on a footing closer to Kuramoto and open a clear experimental route via annular BZ droplets.

major comments (2)
  1. [Sec. III (adaptation and assumptions (i)–(ii))] Sec. III: the normal-form claim rests on rewriting Tanaka’s decay as (S−S̄)/τ (explicitly “to ensure the reduction maps onto” Eqs. 1–2) and then truncating to first-harmonic zero lag so G(x)≈g cos x. Both steps are transparent but load-bearing for genericity. The manuscript should state more sharply which physical regimes (ring geometry, source structure, τ window) make higher harmonics and lag negligible, and what qualitative changes are expected when they are not, so the “generic” claim is not overstated.
  2. [Sec. IV.B.3 (near-cusp asymptotics)] Sec. IV.B.3, Eqs. (38)–(39): the near-cusp slope j*=7/20 is obtained by setting the exchange-odd inner-flow constant ΔH=6/5, which is described as “consistent with” a direct evaluation but left without a closed-form proof. Because the organizing cusp is a central claim of the phase diagram, either a short derivation of ΔH or an explicit statement that j* remains conjectural (with the numerical evidence that supports it) is needed before the slope can be treated as analytic.
minor comments (4)
  1. [Fig. 4] Fig. 4 caption and surrounding text: the dashed sync boundary is labeled “semi-analytic” and “anchored to simulation markers.” Clarify in the caption which parts of J3(K) are analytic (large-J asymptotics, cusp slope) versus numerically solved from H=0.
  2. [Sec. IV.C] Sec. IV.C: mixed-state order parameters are obtained by self-consistency with a small numerical both-drift remainder. A one-sentence statement of the residual size (already given as ≲0.3% of r) and the grid size used would help readers assess the numerical error bar.
  3. [Sec. IV (opening paragraphs)] Notation: the same symbols (K,J) are used for the original couplings and for the rotated (K,J)=((J′+K′)/2,(J′−K′)/2). A brief reminder at first use in Sec. IV would reduce momentary confusion.
  4. [Abstract / Sec. IV intro] The non-monotonic r(K) for J>Jc is highlighted in the abstract and introduction; a short sentence quantifying the depth of the dip (or pointing to the mixed-branch self-consistency) would make the claim easier to verify from the figures alone.

Circularity Check

1 steps flagged

No load-bearing circularity: normal-form recovery is conditional on stated assumptions; stability and order-parameter results are independent continuum calculations checked against simulation.

specific steps
  1. other [Sec. III, adaptation paragraph and assumptions (i)–(ii)]
    "Tanaka used S/τ instead; we replace it with (S−S̄)/τ to ensure the reduction maps onto the 1D swarmalator model (as we soon show). ... Note that G(x)=cos x recovers the 1D swarmalator model. We now show this is precisely what emerges under two independent assumptions. (i) First harmonic. ... (ii) Zero lag."

    The adaptation form and first-harmonic/zero-lag truncations are chosen so the reduction recovers the pre-existing toy model whose genericity is then claimed. This is a deliberate, disclosed modeling step rather than a hidden self-definition of a numerical prediction; it does not force the subsequent stability boundaries or order-parameter corrections, which are derived independently from the continuum equations of the recovered model.

full rationale

The paper’s derivation chain does not reduce its central predictions to their inputs by construction. The normal-form claim is explicitly conditional: after adapting Tanaka’s chemotactic reduction, the authors replace S/τ by (S−S̄)/τ and then impose first-harmonic plus zero-lag so that G(x)≈g cos x recovers Eqs. (1)–(2). That is a transparent modeling choice (biologically motivated by contrast sensing), not a self-definitional prediction of a quantity already fitted or defined as the answer. Once the 1D model is accepted, the phase-wave boundary F(K,J)=0 is obtained from continuum linearization of the continuity equation around the two-piece density, collapsing to an elementary arctangent self-consistency; the sync order parameter is recomputed as r_lock+r_tongue (correcting the prior OA formula of Yoon et al.); the cusp at (4,2) and the large-J asymptotic K3=J+(2/π)log(2J)+4/π follow from those same expansions. All are checked against independent N=10^6 simulations. Self-citation of the authors’ earlier swarmalator papers is present but not load-bearing for the new stability or corrected order-parameter results; Tanaka’s reduction is external. Mixed-state stability and a fully closed cusp slope remain open and are acknowledged. No fitted-input-called-prediction, uniqueness-import, or ansatz-smuggling step forces the main claims. Score 1 only for the mild, disclosed engineering of the adaptation term so the reduction lands on the pre-existing toy model—proportionate and non-central.

Axiom & Free-Parameter Ledger

2 free parameters · 6 axioms · 0 invented entities

The central claims rest on Tanaka’s chemotactic reduction plus two modeling limits (first harmonic, zero lag), a biologically motivated but target-selected adaptation form, Cauchy disorder (rescaled to width 1), and the infinite-N continuum limit. Stability analysis uses standard continuity-equation linearization and finite-rank order-parameter feedback. No new physical entities are postulated. A few asymptotic constants (ε² coefficient, ΔH) are not fully closed-form and are the main soft spots; they affect the cusp slope detail, not the existence of the four states or the phase-wave boundary.

free parameters (2)
  • cusp inner-flow constant ΔH (taken as 6/5) = 6/5 (conjectural)
    Sets the near-cusp slope j*=1/2−ΔH/8 of the sync boundary. Authors report a direct evaluation consistent with 6/5 but leave a closed-form proof to future work; the slope is therefore not fully derived.
  • phase-wave near-cusp ε² coefficient −31/64 = −31/64 (numerical estimate)
    Appears in the small-ε expansion of J2(K); stated as a high-precision numerical estimate rather than an exact coefficient, so it is not load-bearing for the leading cusp structure.
axioms (6)
  • domain assumption Chemical response after adaptation is dominated by the first spatial harmonic: G(x)≈g1 cos x with higher modes suppressed.
    Sec. III assumption (i); required to recover the sin/cos kernels of the 1D swarmalator model.
  • domain assumption Zero chemical lag: ωτ/(1+Dτ)≪1 so g1 is real.
    Sec. III assumption (ii); without it the kernels acquire a phase lag α and the exact model (1)–(2) is not recovered.
  • ad hoc to paper Adaptation is modeled by (S−S̄)/τ rather than Tanaka’s S/τ, so the uniform chemical background is removed.
    Sec. III explicitly replaces Tanaka’s term “to ensure the reduction maps onto the 1D swarmalator model”; biologically motivated by contrast sensing but selected for the target equations.
  • domain assumption Natural frequencies (v′i,ω′i) are i.i.d. Cauchy with common width Δ, set to 1 by time rescaling; transformed (ν,μ) are Lorentzian of half-width 2.
    Sec. II and footnote 1; enables closed-form self-consistency. Robustness to other full-support densities is claimed in SM but not proved for compact support.
  • standard math Infinite-N continuum limit: density ρ obeys the continuity equation with mean-field velocities depending only on the two complex order parameters.
    Used throughout Sec. IV for linearization around async, phase wave, and sync; standard mean-field assumption.
  • domain assumption Tanaka’s weakly coupled chemotactic Hopf reduction (near-onset amplitude elimination then phase reduction) is valid on a 1D ring.
    Sec. III builds the normal form on this framework; validity near Hopf onset is inherited from that literature.

pith-pipeline@v1.1.0-grok45 · 17062 in / 3971 out tokens · 45180 ms · 2026-07-14T15:23:04.609222+00:00 · methodology

0 comments
read the original abstract

Swarmalators are mobile generalizations of phase oscillators. Introduced to model systems in which sync and self-assembly interact, they remain poorly understood theoretically. Unlike the Kuramoto model for coupled oscillators, existing swarmalator models lack a normal-form foundation, and their basic stabilities and bifurcations remain largely unsolved. Here we address both problems. Building on Tanaka's reduction of chemotactic oscillators, we show that the canonical one-dimensional swarmalator model -- previously introduced as an ad hoc toy model -- is recovered in the first-harmonic, zero-lag limit, implying its behavior is generic. We then derive the stability boundaries organizing its four collective states, show they meet at a single cusp, correct a previously published order-parameter formula, and uncover a non-monotonic sync response absent in the Kuramoto model.

Figures

Figures reproduced from arXiv: 2607.09810 by Kevin P. O'Keeffe.

Figure 1
Figure 1. Figure 1: FIG. 1. Kuramoto and the swarmalator model have parallel [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Scatter plots of the four states in the ( [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Order parameters versus coupling [PITH_FULL_IMAGE:figures/full_fig_p002_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Phase diagram in the ( [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Sync order parameter vs [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

30 extracted references

  1. [1]

    [7] derived the sync order parame- ter via the Ott–Antonsen ansatz, obtainingS OA =p 1−4K/(K 2 −J 2)

    Order parameter Yoon et al. [7] derived the sync order parame- ter via the Ott–Antonsen ansatz, obtainingS OA =p 1−4K/(K 2 −J 2). Here we find this result is incor- rect for generalJ: it holds only on the decoupled line J= 0. ForJ̸= 0 it undercounts the true order pa- rameter, discarding a partially-locked population (the “tongue”—oscillators withξlocked ...

  2. [2]

    In rotated coordinatesp= (ξ+η)/2, q= (ξ−η)/2 with Ω = (ν+µ)/2, Λ = (ν−µ)/2, the frozen flow decouples into ˙p= Ω−Asinpcosq and ˙q= Λ−Bcospsinq

    Stability and the locking diamond The sync state’s locking region is not a square but a diamond. In rotated coordinatesp= (ξ+η)/2, q= (ξ−η)/2 with Ω = (ν+µ)/2, Λ = (ν−µ)/2, the frozen flow decouples into ˙p= Ω−Asinpcosq and ˙q= Λ−Bcospsinq. A locked point requires u≡Ω/A= sinpcosqandv≡Λ/B= cospsinq, which satisfy|u|+|v| ≤1, i.e. |Ω| S(K+J) + |Λ| S(K−J) ≤1....

  3. [3]

    Asymptotics We solveH= 0 numerically in the interior and ana- lytically in two limits. LargeJ.Theη-locking edge becomes a thin barrier near cosq= 0; rescaling the barrier coordinatey, the first-harmonic response obeys (2 +iy)Z 1 = 1, so the in- quadrature kernel isK(y) =y/(y 2 + 4), whose finite part FP Z ∞ 0 h y y2 + 4 − 1y>1 y i dy=−log 2 sets theO(1) o...

  4. [4]

    Kuramoto, inInternational Symposium on Mathemat- ical Problems in Theoretical Physics, edited by H

    Y. Kuramoto, inInternational Symposium on Mathemat- ical Problems in Theoretical Physics, edited by H. Araki, Lecture Notes in Physics Vol. 39 (Springer, Berlin, 1975), pp. 420–422

    1975

  5. [5]

    S. H. Strogatz, From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators, Physica D143, 1 (2000)

    2000

  6. [6]

    J. A. Acebr´ on, L. L. Bonilla, C. J. P´ erez Vicente, F. Ritort, and R. Spigler, The Kuramoto model: A sim- ple paradigm for synchronization phenomena, Rev. Mod. Phys.77, 137 (2005)

    2005

  7. [7]

    Wiesenfeld, P

    K. Wiesenfeld, P. Colet, and S. H. Strogatz, Synchroniza- tion transitions in a disordered Josephson series array, Phys. Rev. Lett.76, 404 (1996)

    1996

  8. [8]

    K. P. O’Keeffe, H. Hong, and S. H. Strogatz, Oscillators that sync and swarm, Nat. Commun.8, 1504 (2017)

    2017

  9. [9]

    O’Keeffe, S

    K. O’Keeffe, S. Ceron, and K. Petersen, Collective behav- ior of swarmalators on a ring, Phys. Rev. E105, 014211 (2022)

    2022

  10. [10]

    S. Yoon, K. P. O’Keeffe, J. F. F. Mendes, and A. V. Goltsev, Sync and swarm: Solvable model of nonidentical swarmalators, Phys. Rev. Lett.129, 208002 (2022). 8

    2022

  11. [11]

    Tanaka, General chemotactic model of oscillators, Phys

    D. Tanaka, General chemotactic model of oscillators, Phys. Rev. Lett.99, 134103 (2007)

    2007

  12. [12]

    S. H. Strogatz and R. E. Mirollo, Stability of incoherence in a population of coupled oscillators, J. Stat. Phys.63, 613 (1991)

    1991

  13. [13]

    Mirollo and S

    R. Mirollo and S. H. Strogatz, The spectrum of the par- tially locked state for the Kuramoto model, J. Nonlinear Sci.17, 309 (2007)

    2007

  14. [14]

    Barkai and S

    N. Barkai and S. Leibler, Robustness in simple biochem- ical networks, Nature387, 913 (1997)

    1997

  15. [15]

    Delgado, N

    J. Delgado, N. Li, M. Leda, H. O. Gonz´ alez-Ochoa, S. Fraden, and I. R. Epstein, Coupled oscillations in a 1D emulsion of Belousov–Zhabotinsky droplets, Soft Matter 7, 3155 (2011)

    2011

  16. [16]

    Torbensen, F

    K. Torbensen, F. Rossi, S. Ristori, and A. Abou- Hassan, Chemical communication and dynamics of droplet emulsions in networks of Belousov–Zhabotinsky micro-oscillators produced by microfluidics, Lab Chip17, 1179 (2017)

    2017

  17. [17]

    Thutupalli and S

    S. Thutupalli and S. Herminghaus, Tuning active emul- sion dynamics via surfactants and topology, Eur. Phys. J. E36, 91 (2013)

    2013

  18. [18]

    N. J. Suematsu, Y. Mori, T. Amemiya, and S. Nakata, Oscillation of speed of a self-propelled Belousov– Zhabotinsky droplet, J. Phys. Chem. Lett.7, 3424 (2016)

    2016

  19. [19]

    Heuthe, P

    V.-L. Heuthe, P. Iyer, G. Gompper, and C. Bechinger, Tunable colloidal swarmalators with hydrodynamic cou- pling, Nat. Commun.16, 10984 (2025)

    2025

  20. [20]

    S. G. Leyva, Z. Zhang, M. Olvera de la Cruz, and K. J. M. Bishop, Self-oscillating synchronematic colloids, Nat. Commun.17, 1841 (2026)

    2026

  21. [21]

    H. Hong, K. P. O’Keeffe, J. S. Lee, and H. Park, Swar- malators with thermal noise, Phys. Rev. Research5, 023105 (2023)

    2023

  22. [22]

    M. S. Anwar, D. Ghosh, and K. O’Keeffe, Forced one- dimensional swarmalator model, Phys. Rev. E110, 054205 (2024)

    2024

  23. [23]

    N. Blum, A. Li, K. O’Keeffe, and O. Kogan, Swarmala- tors with delayed interactions, Phys. Rev. E109, 014205 (2024)

    2024

  24. [24]

    K. P. O’Keeffe, R. T. Djeudjo, and J. Hindes, Time delay in the 1D swarmalator model, arXiv:2602.08156 (2026)

    arXiv 2026

  25. [25]

    G. K. Sar, D. Ghosh, and K. O’Keeffe, Pinning in a sys- tem of swarmalators, Phys. Rev. E107, 024215 (2023)

    2023

  26. [26]

    G. K. Sar, D. Ghosh, and K. O’Keeffe, Solvable model of driven matter with pinning, Phys. Rev. E109, 044603 (2024)

    2024

  27. [27]

    B. Hao, M. Zhong, and K. O’Keeffe, Attractive and re- pulsive interactions in the one-dimensional swarmalator model, Phys. Rev. E108, 064214 (2023)

    2023

  28. [28]

    M. S. Anwar, G. K. Sar, M. Perc, and D. Ghosh, Collec- tive dynamics of swarmalators with higher-order interac- tions, Commun. Phys.7, 59 (2024)

    2024

  29. [29]

    Ghosh, K

    S. Ghosh, K. O’Keeffe, G. K. Sar, and D. Ghosh, Dy- namics of pulsating swarmalators on a ring, Phys. Rev. E112, 054217 (2025)

    2025

  30. [30]

    G. K. Sar, K. O’Keeffe, and D. Ghosh, Dynamics of swar- malators in the presence of a contrarian, Phys. Rev. E 111, 014209 (2025)

    2025