Vicsek model

The Vicsek model is a mathematical model used to describe active matter. One motivation of the study of active matter by physicists is the rich phenomenology associated to this field. Collective motion and swarming are among the most studied phenomena. Within the huge number of models that have been developed to catch such behavior from a microscopic description, the most famous is the model introduced by Tamás Vicsek et al. in 1995. ^[1]

Physicists have a great interest in this model as it is minimal and describes a kind of universality. It consists in point-like self-propelled particles that evolve at constant speed and align their velocity with their neighbours' one in presence of noise. Such a model shows collective motion at high density of particles or low noise on the alignment.

Model (mathematical description)

As this model aims at being minimal, it assumes that flocking is due to the combination of any kind of self propulsion and of effective alignment.Since velocities of each particle is a constant, the net momentum of the system is not conserved during collisions.

An individual ${\displaystyle i}$ is described by its position ${\displaystyle \mathbf {r} _{i}(t)}$ and the angle defining the direction of its velocity ${\displaystyle \Theta _{i}(t)}$ at time ${\displaystyle t}$. The discrete time evolution of one particle is set by two equations:

1. At each time step ${\displaystyle \Delta t}$, each agent aligns with its neighbours within a given distance ${\displaystyle r}$ with an uncertainty due to a noise ${\displaystyle \eta _{i}(t)}$:
   • ${\displaystyle \Theta _{i}(t+\Delta t)=\langle \Theta _{j}\rangle _{|r_{i}-r_{j}|<r}+\eta _{i}(t)}$
2. The particle then moves at constant speed ${\displaystyle v}$ in the new direction:
   • ${\displaystyle \mathbf {r} _{i}(t+\Delta t)=\mathbf {r} _{i}(t)+v\Delta t{\begin{pmatrix}\cos \Theta _{i}(t)\\\sin \Theta _{i}(t)\end{pmatrix}}}$

In these equations, ${\displaystyle \langle \Theta _{j}\rangle _{|r_{i}-r_{j}|<r}}$ denotes the average direction of the velocities of particles (including particle ${\displaystyle i}$) within a circle of radius ${\displaystyle r}$ surrounding particle ${\displaystyle i}$. The average normalized velocity acts as the order parameter for this system, and is given by ${\displaystyle v_{a}={\frac {1}{Nv}}|\sum _{i=1}^{N}v_{i}|}$.

The whole model is controlled by three parameters: the density of particles, the amplitude of the noise on the alignment and the ratio of the travel distance ${\displaystyle v\Delta t}$ to the interaction range ${\displaystyle r}$ . From these two simple iteration rules, various continuous theories ^[2] have been elaborated such as the Toner Tu theory ^[3] which describes the system at the hydrodynamic level. An Enskog-like kinetic theory, which is valid at arbitrary particle density, has been developed. ^[4] This theory quantitatively describes the formation of steep density waves, also called invasion waves, near the transition to collective motion. ^[5]

Phenomenology

This model shows a phase transition ^[6] from a disordered motion to large-scale ordered motion. At large noise or low density, particles are on average not aligned, and they can be described as a disordered gas. At low noise and large density, particles are globally aligned and move in the same direction (collective motion). This state is interpreted as an ordered liquid. The transition between these two phases is not continuous, indeed the phase diagram of the system exhibits a first order phase transition with a microphase separation. In the co-existence region, finite-size liquid bands ^[7] emerge in a gas environment and move along their transverse direction. Recently, a new phase has been discovered: a polar ordered Cross sea phase of density waves with inherently selected crossing angle. ^[8] This spontaneous organization of particles epitomizes collective motion.

Extensions

Since its appearance in 1995 this model has been very popular within the physics community; many scientists have worked on and extended it. For example, one can extract several universality classes from simple symmetry arguments concerning the motion of the particles and their alignment. ^[9]

Moreover, in real systems, many parameters can be included in order to give a more realistic description, for example attraction and repulsion between agents (finite-size particles), chemotaxis (biological systems), memory, non-identical particles, the surrounding liquid.

A simpler theory, the Active Ising model, ^[10] has been developed to facilitate the analysis of the Vicsek model.

References

1. Vicsek, Tamás; Czirók, András; Ben-Jacob, Eshel; Cohen, Inon; Shochet, Ofer (1995-08-07). "Novel Type of Phase Transition in a System of Self-Driven Particles". Physical Review Letters. 75 (6): 1226–1229. arXiv: cond-mat/0611743 . Bibcode:1995PhRvL..75.1226V. doi:10.1103/PhysRevLett.75.1226. PMID   10060237. S2CID   15918052.
2. Bertin, Eric; Droz, Michel; Grégoire, Guillaume (2006-08-02). "Boltzmann and hydrodynamic description for self-propelled particles". Physical Review E. 74 (2): 022101. arXiv: cond-mat/0601038 . Bibcode:2006PhRvE..74b2101B. doi:10.1103/PhysRevE.74.022101. PMID   17025488. S2CID   19658705.
3. Toner, John; Tu, Yuhai (1995-12-04). "Long-Range Order in a Two-Dimensional Dynamical $\mathrm{XY}$ Model: How Birds Fly Together". Physical Review Letters. 75 (23): 4326–4329. Bibcode:1995PhRvL..75.4326T. doi:10.1103/PhysRevLett.75.4326. PMID   10059876.
4. Ihle, Thomas (2011-03-16). "Kinetic theory of flocking: Derivation of hydrodynamic equations". Physical Review E. 83 (3): 030901. arXiv: 1006.1825 . Bibcode:2011PhRvE..83c0901I. doi: 10.1103/PhysRevE.83.030901 . PMID   21517447.
5. Ihle, Thomas (2013-10-18). "Invasion-wave-induced first-order phase transition in systems of active particles". Physical Review E. 88 (4): 040303. arXiv: 1304.0149 . Bibcode:2013PhRvE..88d0303I. doi:10.1103/PhysRevE.88.040303. PMID   24229097. S2CID   14951536.
6. Grégoire, Guillaume; Chaté, Hugues (2004-01-15). "Onset of Collective and Cohesive Motion". Physical Review Letters. 92 (2): 025702. arXiv: cond-mat/0401208 . Bibcode:2004PhRvL..92b5702G. doi:10.1103/PhysRevLett.92.025702. PMID   14753946. S2CID   37159324.
7. Solon, Alexandre P.; Chaté, Hugues; Tailleur, Julien (2015-02-12). "From Phase to Microphase Separation in Flocking Models: The Essential Role of Nonequilibrium Fluctuations". Physical Review Letters. 114 (6): 068101. arXiv: 1406.6088 . Bibcode:2015PhRvL.114f8101S. doi:10.1103/PhysRevLett.114.068101. PMID   25723246. S2CID   43186543.
8. Kürsten, Rüdiger; Ihle, Thomas (2020-10-30). "Dry Active Matter Exhibits a Self-Organized Cross Sea Phase". Physical Review Letters. 125 (18): 188003. arXiv: 2002.03198 . Bibcode:2020PhRvL.125r8003K. doi:10.1103/PhysRevLett.125.188003. PMID   33196272. S2CID   211069694.
9. Chaté, H.; Ginelli, F.; Grégoire, G.; Peruani, F.; Raynaud, F. (2008-07-11). "Modeling collective motion: variations on the Vicsek model" (PDF). The European Physical Journal B. 64 (3–4): 451–456. Bibcode:2008EPJB...64..451C. doi:10.1140/epjb/e2008-00275-9. ISSN   1434-6028. S2CID   49363896.
10. Solon, A. P.; Tailleur, J. (2013-08-13). "Revisiting the Flocking Transition Using Active Spins". Physical Review Letters. 111 (7): 078101. arXiv: 1303.4427 . Bibcode:2013PhRvL.111g8101S. doi:10.1103/PhysRevLett.111.078101. PMID   23992085.