Minimal Swarmalator System Fundamentals

• Minimal swarmalator systems are dynamical models defined by bidirectional coupling between spatial movement and phase synchronization on periodic domains.
• They employ a reduced Fourier formulation with two key parameters, J and K, to capture collective states such as synchrony, phase waves, and chaotic regimes.
• These models bridge concepts from Kuramoto and Vicsek dynamics, providing a tractable benchmark for exploring self-organization in active multi-agent systems.

A minimal swarmalator system is a dynamical system of agents (“swarmalators”) that combines collective swarming in physical space with mutual synchronization of internal phases, through explicit bidirectional coupling between positions and phases. “Minimal” refers to the simplest nontrivial construction that retains the irreducible features: (i) explicit, reciprocal interplay between spatial and phase dynamics; (ii) analytically tractable form; and (iii) the ability to support multiple collective states, such as synchrony, asynchronous dispersion, phase waves, and active/chaotic regimes. In minimal models, all agents are typically identical, interactions are mean-field (all-to-all), and all variables are restricted to compact spaces such as circles or tori (e.g., $x, y, \theta \in \mathbb{S}^1$).

1. Mathematical Formulation of Minimal Swarmalator Models

The canonical minimal swarmalator model in two spatial dimensions with periodic boundary conditions—the “solvable 2D swarmalator model”—is defined as follows. Let $N$ identical units carry coordinates $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$ and phase $\theta_i\in\mathbb{S}^1$. The dynamics are

$$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$

Here, $J \in \mathbb{R}$ couples phase similarity to spatial attraction; $K \in \mathbb{R}$ couples spatial proximity to phase synchronization. All agents live on a flat torus, allowing the entire model to be written in first Fourier harmonics (sines and cosines). Repulsive terms and higher harmonics are omitted, maximizing tractability. The model generalizes naturally to one dimension and to variants with noise, external forcing, or parameter disorder.

2. Minimality: Structural and Analytical Considerations

Key minimality features:

• Dimensional Reduction and Periodicity: By confining positions and phases to $[0,2\pi)$, the periodic geometry eliminates boundary effects and makes all couplings Fourier-analytic ($\sin$, $\cos$), enabling explicit linear stability, bifurcation analyses, and closed-form expressions.
• Reduced Coupling Structure: Only two tunable couplings appear: $N$0 (space$N$1phase coupling in spatial motion) and $N$2 (space$N$3phase coupling in phase evolution). No additional parameters or agent heterogeneity exist in the strict minimal form.
• Symmetry and Analytical Tractability: The system possesses permutation symmetry and translation invariance in each variable, ensuring the block-diagonalization of Jacobians and the solvability of spectral problems.
• Fourier-based Representation: Expanding all interaction terms in first harmonics retains only the minimal nontrivial bidirectional space–phase cross-couplings required for qualitative richness.

These constraints yield a system that, while reduced and structurally simple, exhibits the essential phenomenology of swarmalation.

3. Collective States and Bifurcations

Minimal swarmalator systems support a finite set of macroscopic dynamical states, each demarcated by analytically computable boundaries. For the solvable 2D model, varying $N$4 at fixed $N$5 yields the following sequence $N$6

1. Static Asynchronous State ($N$7): All agents evenly fill phase and spatial space, $N$8, $N$9.
2. Thick Phase Wave ($(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$0): A saddle–node bifurcation at $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$1 creates a state with nonzero $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$2 (meaning $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$3 correlation), $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$4.
3. Thin Phase Wave ($(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$5): For $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$6 the phase wave narrows further; both spatial coordinates and phases wind around $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$7 in locked synchrony modulo $(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$8.
4. Full Synchrony (Clustered State) ($(x_i, y_i)\in\mathbb{S}^1\times\mathbb{S}^1$9): All particles collapse to a point in $\theta_i\in\mathbb{S}^1$0-space; this state is stable iff $\theta_i\in\mathbb{S}^1$1.

The explicit stability thresholds are:

State	Existence Condition	Stability Boundaries
Synchronized (clustered)	$\theta_i\in\mathbb{S}^1$2, $\theta_i\in\mathbb{S}^1$3	Jacobian eigenvalues: $\theta_i\in\mathbb{S}^1$4
Thin phase wave	$\theta_i\in\mathbb{S}^1$5	$\theta_i\in\mathbb{S}^1$6 from spectrum
Thick phase wave	$\theta_i\in\mathbb{S}^1$7	Saddle-node at $\theta_i\in\mathbb{S}^1$8
Static async	$\theta_i\in\mathbb{S}^1$9	Destabilizes at $$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$0

These transitions match continuum and finite-$$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$1 calculations.

4. Order Parameters and Macroscopic Diagnostics

The principal observables are the “rainbow” order parameters, capturing correlations among space–phase “sum” and “difference” coordinates. Let

$$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$2

Define

$$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$3

$$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$4, $$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$5 are amplitudes measuring spatial–phase correlations along $$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$6, $$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$7.

In practice, $$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$8 suffices to summarize macroscopic coherence. $$\begin{aligned} \dot x_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(x_j - x_i)\,\cos(\theta_j-\theta_i), \ \dot y_i &=\; \frac{J}{N}\sum_{j=1}^N \sin(y_j - y_i)\,\cos(\theta_j-\theta_i), \ \dot \theta_i &=\; \frac{K}{N}\sum_{j=1}^N \sin(\theta_j-\theta_i)\ [\cos(x_j - x_i)+\cos(y_j - y_i)] . \end{aligned}$$9 implies perfect space–phase locking; $J \in \mathbb{R}$0 represents fully incoherent (uniform) states.

5. Connections to Kuramoto and Vicsek Models

The minimal swarmalator system sits at the intersection of synchronization and active matter:

• Kuramoto Limit: If $J \in \mathbb{R}$1, the model reduces to coupled Kuramoto oscillators on a torus, with phases synchronizing under interaction terms weighted by spatial proximity. The phase equation generalizes Kuramoto,

$J \in \mathbb{R}$2

• Vicsek Limit: In contrast to the Vicsek model ($J \in \mathbb{R}$3, $J \in \mathbb{R}$4 plus noise), spatial motion in the minimal swarmalator system arises from pairwise cross-couplings, and heading (phase) affects motion only via $J \in \mathbb{R}$5. The interplay is genuinely bidirectional.

These correspondences elucidate why the minimal model is analytically tractable yet phenomenologically rich: it contains both the essential structure of phase oscillators and of collectively moving agents, but with only two degrees of freedom per agent (space and phase) and no extraneous parameters.

6. Minimality, Variants, and Future Extensions

The “minimal” moniker is justified by several criteria:

• Parameter Reduction: Only two free couplings ($J \in \mathbb{R}$6), no agent-level disorder or noise.
• Periodic Geometry: By taking all variables on tori/circles, the model enables spectral analysis and avoids boundary artifacts present in Euclidean space.
• Truncation to Lowest Fourier Modes: Omitting higher harmonics and non-Fourier couplings preserves only the simplest nontrivial forms necessary for the space–phase feedback loop.
• Absence of Repulsion: Exclusion of explicit repulsive kernels further simplifies analysis, yet the model still supports all principal collective states.

Future directions proposed include:

• Introduction of natural frequency heterogeneity or spatial drift, analyzable via generalized Ott–Antonsen ansätze.
• Additional Fourier modes, higher-dimensional tori ($J \in \mathbb{R}$7), or reintroduction of weak repulsion, which, according to weak-kernel analysis, does not qualitatively alter bifurcation structure.
• Experimental realization in motile agents with internal oscillators, subject to periodic geometries and tunable couplings.

7. Broader Significance and Open Problems

Minimal swarmalator models serve as a “solvable playground” for the study of bidirectional coupling between swarming and synchronization, occupying a role analogous to that of the Kuramoto model for internal synchronization phenomena. They furnish explicit order parameters and phase diagrams for the emergence and stability of macroscopic patterns. The analytic solvability provides benchmarks for more complex models (including real biological or robotic swarms) and a foundation for further work on fluctuations, heterogeneity, finite-size effects, higher-order couplings, and externally driven or spatially disordered systems.

Outstanding open challenges include:

• Analytical derivation of the full structure of active (non-stationary) states, including transient “active phase waves” and transitions to chaos in high-dimensional phase space.
• Systematic characterization of multistability and bifurcation mechanisms under parameter variation.
• Extension to networks with finite-range or complex topology, or coupling disorder.
• Experimental design principles for engineering swarmalator-like behavior in synthetic active matter or distributed robotic systems.

As the first analytically tractable, bidirectionally coupled space–phase oscillator model, the minimal swarmalator system delineates the essential ingredients for spatiotemporal self-organization in coupled active systems and provides a reference platform for theoretical and applied investigations (O'Keeffe et al., 2023, Sar et al., 10 Oct 2025, Sar et al., 2022, Lizarraga et al., 2023, O'Keeffe et al., 2019).