Swarmalators

Swarmalators ^[1] are generalizations of phase oscillators ^[2] that swarm around in space as they synchronize in time. They were introduced to model the diverse real-world systems which both sync and swarm, such as vinegar eels, ^[3] magnetic domain walls, ^[4] and Japanese tree frogs. ^[5] More formally, they are dynamical units with spatial degrees of freedom and internal degrees of freedom whose dynamics are coupled.

Real world examples

Swarmalation ^[6] occurs in diverse parts of Nature and technology some of which are discussed below. The Figure to the right plots some examples in a (discipline, number of particles) plot.

Examples of real world swarmalators plotted in a (discipline, number of particles) plane

Biological microswimmers. Sperm, vinegar eels and potentially other swimmers such as C elegans swarm through space via the rhythmic beating of their tails. This beating may synchronize with the beating of a neighboring swimmer via hydrodynamic coupling, which in turn causes spatial attraction; sync links to self-assembly. This can lead to vortex arrays, ^[7] trains ^[8] metachronal waves ^[9] and other collective effects.

Magnetic domain walls ^[10] are key features in the field of magnetism and materials science, defined by the boundary between different magnetic domains in ferromagnetic materials. These domains are regions within a material where the magnetic moments of atoms are aligned in the same direction, creating a uniform magnetic field. The hold great promise as memory devices in next generation spintronics. In a simplified ${\displaystyle q,\phi }$ model, ^[11] a domain wall can be described by its center of mass ${\displaystyle q}$ and the in-plane angle ${\displaystyle \phi }$ of its magnetic dipole vector, thereby classifying them as swarmalators. Experiments reveal that the interaction between two such domain walls leads to rich spatiotemporal behaviors some of which is captured by the 1D swarmalator model listed above. ^[12]

Japanese Tree frogs. During courtship rituals, male Japanese Tree frogs attract the attention of females by croaking rhythmically. Neighboring males tend to alternate the croaking (croak ${\displaystyle \pi }$ degree out of phase) so as to avoid "speaking over each other". Evidence ^[13] suggests this (anti)-synchronization influences thee inter-frog spatial dynamics, making them swarmalators.

Janus particles ^[14] are spherical particles with one hemisphere coated in a magnetic substance, the other remaining non-magnetic. They are named after the Roman God Janus who has two faces. This anisotropy gives the particles unusual magnetic properties. When subject to external magnetic fields, their magnetic dipoles vectors begins to oscillate which induces and couples to movements (thus qualifying as a swarmalators). The resultant "sync-selected self-assembly" ^[15] gives rise novel superstructure with potential use in biomedicine contexts such as targeted drug delivery, bio imaging, and bio-sensing. ^[16]

Quincke rollers^{[ citation needed ]} are a class of active particle that exhibits self-propelled motion in a fluid due to an electrohydrodynamic phenomenon known as the Quincke effect. ^[17] This effect occurs when a dielectric (non-conducting) particle is subject to an electric field. The rotation of the particle, combined with frictional interactions with the surrounding fluid and surface, leads to a rolling motion. Thus, the particle has a phase ${\displaystyle \theta }$ and a position ${\displaystyle x}$ which couple, as required of swarmalator. Collections of Quincke rollers produce rich emergent behavior such as activity waves ^[18] and shock waves. ^[19]

Embryonic cells are the foundational building blocks of an embryo, undergoing division and differentiation to form the complex structures of an organism. These cells exhibit remarkable plasticity, allowing them to transform into a wide range of specialized cell types. In the context of swarmalators, embryonic cells display a unique blend of synchronization and swarming behaviors. ^[20] They coordinate their movements and genetic expression patterns in response to various cues, a process essential for proper tissue formation and organ development. This linking of sync and self-assembly make embryonic cells a compelling example of a real-world swarmalators.

Robot swarms. Land based rovers as well as aerial drones programmed with swarmalator models have been created ^[21] and has recreated the five collective states of the swarmalator model (see Mathematical Models section for the plot of these states). The linking of sync and swarming defines a new kind of bio-inspired algorithm which several potential applications. ^[22]

2D swarmalator model

A mathematical model for swarmalators moving in 2D has been proposed. This 2D swarmalator model in generic form is

${\displaystyle {\begin{array}{rcl}{\dot {x}}_{i}&=&{\frac {1}{N}}\sum _{j}I_{att}(x_{j}-x_{i})F(\theta _{j}-\theta _{i})-I_{rep}(x_{j}-x_{i})\\{\dot {\theta }}_{i}&=&{\frac {1}{N}}\sum _{j}G_{sync}(\theta _{j}-\theta _{i})H(x_{j}-x_{i})\end{array}}}$

The spatial dynamics combine pairwise interaction ${\displaystyle I_{att}}$ with pairwise ${\displaystyle I_{rep}}$, which produces swarming / aggregation. The novelty is the attraction is modified by a phase term ${\displaystyle F}$; thus the aggregation becomes phase-dependent. Likewise, the phase dynamics contain a sync term ${\displaystyle G_{sync}}$ modified by a spatial term so the synchronization becomes position dependent. In short, the swarmalators model the interaction between self-synchronization and self-assembly in space.

While in general the position could be in 2D or 3D, the instance of the swarmalator model originally introduced is a 2D model ${\displaystyle x\in \mathbb {R} ^{2}}$ and the choices for ${\displaystyle I_{att}(x),I_{rep}(x)}$ etc. were

${\displaystyle {\begin{array}{rcl}{\dot {x}}_{i}&=&{\frac {1}{N}}\sum _{j}{\frac {(x_{j}-x_{i})}{|x_{j}-x_{i}|}}\left(1+J\cos(\theta _{j}-\theta _{i})\right)-{\frac {x_{j}-x_{i}}{|x_{j}-x_{i}|^{2}}}\\{\dot {\theta }}_{i}&=&{\frac {K}{N}}\sum _{j}{\frac {\sin(\theta _{j}-\theta _{i})}{|x_{j}-x_{i}|}}\end{array}}}$

There are two parameters ${\displaystyle J}$ and ${\displaystyle K}$ are parameters: ${\displaystyle J}$ controls the strength of phase-space attraction/repulsion, while ${\displaystyle K}$ describes the phase coupling strength. The above can be considered a blending of the aggregation model introduced from biological swarming ^[23] (the spatial part) and the Kuramoto model of phase oscillators (the phase part).

Phenomena

Figure 1: A depiction of the five collective states of the swarmalator model.

The model above produces 5 collective states depicted in the Figure 1 below.

• Static sync: Swarmalators form a disk in space and are fully synchronized in phase.
• Static async: Swarmalators form a disk in space and are fully asynchronous in phase.
• Static phase wave: Swarmalators form an annulus in space with a full phase wave (e.g. a full color-wheel or rainbow.)
• Splintered phase wave: The phase wave splinters into clusters of synchronous swarmalators. Within each clusters, the swarmalators execute periodic motion in space and phase.
• Active phase wave: Swarmalators run in a space-phase vortex, with half running clockwise and the remaining half running counter-clockwise.

To demarcate where each state arises and disappears as a parameters are changed, the rainbow order parameters,

${\displaystyle {\begin{array}{rcl}W_{\pm }&=&S_{\pm }e^{i\phi _{\pm }}=&{\frac {1}{N}}\sum _{j}e^{i(\phi _{j}-\theta _{i})}\end{array}}}$

Figure 2: The rainbow order parameters and ${\displaystyle \gamma }$ can discern the various collective states. See main text for definitions of these quantities.

where ${\displaystyle \phi _{i}:=\arctan(y_{i}/x_{i})}$ are used. Figure 2 plots ${\displaystyle S_{\pm }}$ versus ${\displaystyle K}$ for fixed ${\displaystyle J}$. As can be seen, ${\displaystyle S_{+}=1}$ in the rainbow-like static phased wave state (at ${\displaystyle K}$ = 0), and then declines as ${\displaystyle K}$ decreases. A second order parameter ${\displaystyle \gamma }$, defined as the fraction of swarmalators that have completed at least one cycle in space and phase after transients in also plotted, which can distinguish between the active phase wave and splintered phase wave states.

Puzzles

There are several unresolved puzzles and open questions related to swarmalators:

• Melting point ${\displaystyle K_{m}(J)}$: What is the value ${\displaystyle K_{m}(J)}$ at which the static async state melts into the active phase wave state?
• Splitting point ${\displaystyle K_{s}(J)}$: What is the splitting point ${\displaystyle K_{s}}$ at which the active phase wave splits into the splintered phase wave?
• Rainbow order parameters: Can you derive an expression for the supercritical branch of the ${\displaystyle S_{\pm }(K)}$ for fixed ${\displaystyle J}$ in the active phase wave and splintered phase wave states?
• ${\displaystyle N_{c}}$: What determines the number of clusters formed in the splintered phase wave?

1D swarmalator model

A simpler swarmalator model where the spatial motion is confined to a 1D ring ${\displaystyle x_{i}\in \mathbb {S} ^{1}}$ has also been proposed ^{[24] [25]}

${\displaystyle {\begin{array}{rcl}{\dot {x}}_{i}&=\nu _{i}+{\frac {J}{N}}\sum _{j}\sin(x_{j}-x_{i})\cos(\theta _{j}-\theta _{i})\\{\dot {\theta }}_{i}&=\omega _{i}+{\frac {K}{N}}\sum _{j}\sin(\theta _{j}-\theta _{i})\cos(x_{j}-x_{i})\end{array}}}$

where ${\displaystyle \nu _{i},\omega _{i}}$ are the (random) natural frequencies of the i-th swarmalator and are drawn from certain distributions ${\displaystyle g(.)}$. This 1D model corresponds to the angular component of the 2D swarmalator model. The restriction to this simpler topology allows for a greater analysis. For instance, the model with natural frequencies can be solved by defining the sum/difference coordinates ${\displaystyle \xi ,\eta =x_{i}\pm \theta _{i}}$ the model simplifies into a pair of linearly coupled Kuramoto models

${\displaystyle {\begin{array}{rcl}{\dot {\xi }}_{i}&=\omega _{+,i}+J_{+}S_{+}\sin(\phi _{+}-\xi )+J_{-}S_{-}\sin(\phi _{-}-\eta )\\{\dot {\eta }}_{i}&=\omega _{-,i}+J_{-}S_{+}\sin(\phi _{+}-\xi )+J_{+}S_{-}\sin(\phi _{-}-\eta )\end{array}}}$

where ${\displaystyle J_{\pm }=(J\pm K)/2}$, ${\displaystyle \omega _{\pm }=\nu \pm \omega }$ and the rainbow order parameters are the equivalent of the 2D model

${\displaystyle {\begin{array}{rcl}W_{\pm }=S_{\pm }e^{i\phi _{\pm }}:=N^{-1}\sum _{j}e^{i(x_{j}\pm \theta _{j})}\\\end{array}}}$

Collective states of 1D swarmalator model, where swarmalators are depicted as dots in the ${\displaystyle (x,\theta )}$ plane.

Top row, sync order parameter ${\displaystyle R}$ versus phase coupling strength ${\displaystyle K}$ for the Kuramoto model with Cauchy natural frequencies with zero center and unit width. Bottom row, rainbow order parameters ${\displaystyle S_{\pm }}$ for the 1D swarmalator model also with unit Cauchy natural frequencies and ${\displaystyle J=9}$. The black dotted lines mark the transitions between the async state, phase wave, and mixed state in that order. Data were collecting using an RK4 solver with timestep ${\displaystyle dt=0.1}$ for ${\displaystyle T=200}$ time units with ${\displaystyle N=10^{6}}$ swarmalators.

For unimodal distribution of ${\displaystyle \omega ,\nu }$ such as the Cauchy distribution, the model exhibits four collective states depicted in the figure on the right.

• Async or ${\displaystyle (S_{+},S_{-})=(0,0)}$ state. Swarmalators do not exhibit any coherence either in space or phase, being distributed uniformly in position and phase. This state is characterized by the absence of any synchronization or spatial clustering among the swarmalators, as reflected by the zero values of both order parameters ${\displaystyle S_{+}=S_{-}=0}$.
• Phase wave or ${\displaystyle (S_{+},S_{-})=(S,0)}$ state. In the phase wave state, swarmalators form a band or wave pattern with the position ${\displaystyle x_{i}}$ and phase ${\displaystyle \theta _{i}}$ are correlated. The wave can either run clockwise or counterclockwise.
• Mixed or ${\displaystyle (S_{+},S_{-})=(S_{1},S_{2})}$ state. Swaramlators again form a phase wave, but now the wave is distorted, forming two rough clusters; thus it is the a mixture of the phase wave and the sync state (described next).
• Sync or ${\displaystyle (S_{+},S_{-})=(S,S)}$ state. Swarmalators form two synchronous cluster in both space and phase. Single cluster states are formed for some initial conditions.

Note in each state, the swarmalators split into a locked/drifting sub-populations, just like the Kuramoto model. The locked population are the denser regions in the Figure, the drifters the light grey regions.

The figure to the right compares the bifurcations of the Kuramoto model to those of the 1D swarmalator model. For the Kuramoto model (top row), the sync order parameter ${\displaystyle R:=|\langle e^{i\theta }\rangle |}$ bifurcates from the async state (${\displaystyle R=0}$) and then increases monontonically in the sync state (${\displaystyle R>0}$). For the 1D swarmalator model, the bifurcations are richer. Starting with the phase coupling ${\displaystyle K<K_{c}}$ and increasing, ${\displaystyle S_{+},S_{-}}$ bifurcate from the async state (${\displaystyle S_{+}=S_{-}=0}$) to the phase wave (${\displaystyle S_{+}=S_{pw},S_{-}=0}$) then to the mixed state (${\displaystyle S_{+},S_{-}=S_{1},S_{2}}$) before finally ending up in the sync state (${\displaystyle S_{+}=S_{-}=S_{sync}}$). Note we have taken ${\displaystyle S_{+}>S_{-}}$ without loss of generality and ${\displaystyle S,S_{1},S_{2}}$ are constants that depend on ${\displaystyle J,K\Delta }$. Expressions for ${\displaystyle S_{pw},S_{sync}}$ have been worked out, those for ${\displaystyle S_{1},S_{2}}$ in the mixed state are unknown (see ref [25]).

References

1. O'Keeffe, Kevin P., Hyunsuk Hong, and Steven H. Strogatz. "Oscillators that sync and swarm." Nature communications 8.1 (2017): 1504.
2. Strogatz, Steven H. "From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators." Physica D: Nonlinear Phenomena 143.1-4 (2000): 1-20.
3. Peshkov, Anton, Sonia McGaffigan, and Alice C. Quillen. "Synchronized oscillations in swarms of nematode Turbatrix aceti." Soft Matter 18.6 (2022): 1174-1182.
4. Hrabec, Aleš, et al. "Velocity enhancement by synchronization of magnetic domain walls." Physical review letters 120.22 (2018): 227204.
5. Aihara, Ikkyu, et al. "Spatio-temporal dynamics in collective frog choruses examined by mathematical modeling and field observations." Scientific reports 4.1 (2014): 3891.
6. Verberck, Bart (2022). "Wavy worms". Nature Physics. 18 (2): 131–131. doi:10.1038/s41567-022-01516-1.
7. Riedel, Ingmar H.; Kruse, Karsten; Howard, Jonathon (2005). "A self-organized vortex array of hydrodynamically entrained sperm cells". Science. 309 (5732): 300–303. doi:10.1126/science.1110329.
8. Schoeller, Simon F.; Holt, William V.; Keaveny, Eric E. (2020). "Collective dynamics of sperm cells". Philosophical Transactions of the Royal Society B. 375 (1807): 20190384. doi: 10.1098/rstb.2019.0384 . PMC   7423380 .
9. Quillen, A. C. (2021). "Metachronal waves in concentrations of swimming Turbatrix aceti nematodes and an oscillator chain model for their coordinated motions". Physical Review E. 104 (1): 014412. arXiv: 2101.06809 . doi:10.1103/PhysRevE.104.014412.
10. Hrabec, Aleš (2018). "Velocity enhancement by synchronization of magnetic domain walls". Physical Review Letters. 120 (22): 227204. arXiv: 1804.01385 . doi:10.1103/PhysRevLett.120.227204.
11. Slonczewski, J. C. (1972). "Dynamics of magnetic domain walls". AIP Conference Proceedings. Vol. 5. American Institute of Physics.
12. Sar, Gourab Kumar; Ghosh, Dibakar; O'Keeffe, Kevin (2023). "Solvable model of driven matter with pinning". arXiv: 2306.09589 .
13. Aihara, Ikkyu (2014). "Spatio-temporal dynamics in collective frog choruses examined by mathematical modeling and field observations". Scientific Reports. 4 (1): 3891. doi: 10.1038/srep03891 . PMC   5384602 .
14. Zhang, Jie; Grzybowski, Bartosz A.; Granick, Steve (2017). "Janus particle synthesis, assembly, and application". Langmuir. 33 (28): 6964–6977. doi:10.1021/acs.langmuir.7b01088.
15. Yan, Jing (2012). "Linking synchronization to self-assembly using magnetic Janus colloids". Nature. 491 (7425): 578–581. doi:10.1038/nature11619.
16. "Viscosity decrease induced by a DC electric field in a suspension". ScienceDirect. Here, the viscosity decrease is obtained by making use of Quincke rotation: the spontaneous rotation of insulating particles suspended in a weakly conducting liquid when the system is submitted to a DC electric field. In such a case, particles rotate around any axis perpendicular to the applied field. ... T.B. Jones Quincke rotation of spheres ...
17. Jones, Thomas B. (1984). "Quincke rotation of spheres". IEEE Transactions on Industry Applications. 4: 845–849.
18. Liu, Zeng Tao (2021). "Activity waves and freestanding vortices in populations of subcritical Quincke rollers". Proceedings of the National Academy of Sciences. 118 (40): e2104724118. doi:10.1073/pnas.2104724118. PMC   8501844 .
19. Zhang, Bo (2023). "Spontaneous shock waves in pulse-stimulated flocks of Quincke rollers". Nature Communications. 14 (1): 7050. doi:10.1038/s41467-023-7050-0.
20. Tsiairis, Charisios D.; Aulehla, Alexander (2016). "Self-organization of embryonic genetic oscillators into spatiotemporal wave patterns". Cell. 164 (4): 656–667. doi: 10.1016/j.cell.2015.12.053 . PMC   4752873 .
21. Barciś, Agata; Bettstetter, Christian (2020). "Sandsbots: Robots that sync and swarm". IEEE Access. 8: 218752–218764.
22. O'Keeffe, Kevin; Bettstetter, Christian (2019). "A review of swarmalators and their potential in bio-inspired computing". Micro-and Nanotechnology Sensors, Systems, and Applications XI. 10982: 383–394.
23. Bernoff, Andrew J., and Chad M. Topaz. "Nonlocal aggregation models: A primer of swarm equilibria." siam REVIEW 55.4 (2013): 709-747.
24. O'Keeffe, Kevin, Steven Ceron, and Kirstin Petersen. "Collective behavior of swarmalators on a ring." Physical Review E 105.1 (2022): 014211.
25. Yoon, S., et al. "Sync and Swarm: Solvable Model of Nonidentical Swarmalators." Physical Review Letters 129.20 (2022): 208002.