Interplay of sync and swarm: Theory and application of swarmalators

TL;DR

This paper surveys swarmalators—agents whose spatial motion and internal phase co-evolve in a bidirectional loop—bridging synchronization and swarming theories. It synthesizes foundational models (Kuramoto, Vicsek, Cucker–Smale, Couzin), introduces the swarmalator framework (notably the 2D model with coupled space and phase dynamics), and reports tractable reductions to 1D ring/line and 2D periodic cases. It highlights key phenomena (static sync/async, static phase wave, splintered and active phase waves), and analyzes extensions including delays, phase lag, higher harmonics, heterogeneous velocities/frequencies, non-Kuramoto dynamics, and higher-order interactions, along with external forcing and predator scenarios. The review emphasizes theoretical barriers to exact analysis (e.g., density Evolution in 2D) while showing progress via solvable lower-dimensional models and OA-type reductions, and discusses predator–prey analogies, local vs nonlocal couplings, and practical applications in robotics, biology, and materials science. Altogether, the work lays a multi-scale, multi-model foundation for understanding and harnessing space–time coherence in complex active matter systems with potential impact on distributed robotics and bio-inspired design.

Abstract

Swarmalators, entities that combine the properties of swarming particles with synchronized oscillations, represent a novel and growing area of research in the study of collective behavior. This review provides a comprehensive overview of the current state of swarmalator research, focusing on the interplay between spatial organization and temporal synchronization. After a brief introduction to synchronization and swarming as separate phenomena, we discuss the various mathematical models that have been developed to describe swarmalator systems, highlighting the key parameters that govern their dynamics. The review also discusses the emergence of complex patterns, such as clustering, phase waves, and synchronized states, and how these patterns are influenced by factors such as interaction range, coupling strength, and frequency distribution. Recently, some minimal models were proposed that are solvable and mimic real-world phenomena. The effect of predators in the swarmalator dynamics is also discussed. Finally, we explore potential applications in fields ranging from robotics to biological systems, where understanding the dual nature of swarming and synchronization could lead to innovative solutions. By synthesizing recent advances and identifying open challenges, this review aims to provide a foundation for future research in this interdisciplinary field.

Figures (29)

• Figure 1: Each panel on the top shows the collection of oscillators situated in the unit circle (when each oscillator $j$ is represented as $e^{i\theta_j}$). The color of each oscillator represents its natural frequency. From left to right we observe how oscillators start to concentrate as the coupling $K$ increases. In the panels below we show the synchronization diagram, i.e., the Kuramoto order parameter $r$ as a function of $K$. It is clear that $K_c=1$ as obtained by using the distribution $g(\omega)$ shown in the right panel. Eq. \ref{['found3']} is used to perform simulations and obtain the results. Source: Reprinted figure with permission from Ref. hermoso2014synchronization.
• Figure 2: Number of oscillators firing as a function of time for $S_0=2$, $\gamma=1$ and $\epsilon=0.3$ and random initial conditions. Time is plotted in multiples of the natural period $T$ of the oscillators. Each period is divided into 10 equal intervals, and the number of oscillators firing during each interval is plotted vertically. For simulation, Eq. \ref{['foundpeskin']} is used with $N=100$ for $S_0=2$, $\gamma=1$ and $\epsilon=0.3$. Source: Reprinted figure with permission from Ref. mirollo1990synchronization.
• Figure 3: Stability diagram for the forced Kuramoto model (Eq. \ref{['found9']}) for $K=5$. (a) Regions A-E correspond to qualitatively different phase portraits. Four types of bifurcations occur: supercritical Hopf bifurcation; homoclinic bifurcation; and two types of saddle-node bifurcations. The filled circle marks a codimension-2 Takens-Bogdanov point, at which the Hopf, homoclinic, and upper saddle-node curve intersect tangentially. (b) Enlargement of the crossover region, where all four bifurcation curves run nearly parallel to one another. (c) Enlargement of the region near the codimension-2 cusp point (filled square), where the upper and lower branches of saddle-node bifurcations meet tangentially. (d) Schematic version of the stability diagram, intended to show how the bifurcation curves connect in the confusing crossover region. Tangential intersections have been opened up for clarity. Source: Reprinted figure with permission from Ref. childs2008stability.
• Figure 4: Simulation of the Vicsek model given by Eqs. \ref{['vicsek1']}-\ref{['vicesk2']}. Velocities of the particles are indicated by small arrows, while their trajectories for the last 20 time steps are shown by a short continuous curve. The number of particles is $N = 300$ and $v=0.03$. (a) initial configuration at $t = 0$. (b) $L = 25$, $\eta = 0.1$; for small densities and noise the particles tend to form groups moving coherently in random directions. (c) $L = 7$, $\eta = 2.0$; at higher densities and noise the particles move randomly with some correlation. (d) $L = 5$, $\eta = 0.1$; for higher density and small noise the motion becomes ordered. Source: Reprinted figure with permission from Ref. vicsek1995novel.
• Figure 5: Patterns formed in a two-dimensional space ($d = 2$) of length $L = 50$, with control parameters $a = 1.8$, $\epsilon = 0.3$, $N = 200$, $R = 15$, and $F = 0.2$. Results are obtained using Eq. \ref{['eq:shibata']}. Panels show the collective states at each four time-steps, from left to right and from top to bottom. Internal states are represented by colors. Source: Reprinted figure with permission from Ref. shibata2003coupled.