• Abstract and Figures
• Public Full-text

On-lattice Vicsek model in confined geometries Andreas Kuhn∗ and Sabine C. Fischery Center for Computationaland Theoretical Biology, Fakult¨atf¨urBiologie, Universit¨atW¨urzburg, Klara-Oppenheimer-Weg 32, 97074 W¨urzburg, Germany (Dated: June 22, 2021) The Vicsek model (Vicsek et al. 1995) is a very popular minimalist model to study active matter with a number of applications to biological systems at different length scales. With its off-lattice implementation and the periodic boundary conditions, it aims at the analysis of bulk behaviour of a limited number of particles. To expand the applicability of the model to further biological systems, we introduce an on-lattice implementation and analyse its behaviour for three different geometries with reflective boundary conditions. For sufficiently fine lattices, the model behaviour does not differ between off-lattice and on-lattice implementation. The reflective boundary conditions introduce an alignment of the particles with the boundary for low levels of noise. Numerical sensitivity analysis of the swarming behaviour results in a detailed characterisation of the Vicsek model for confined geometries with reflective boundary conditions. In a channel geometry, the boundary alignment causes swarms to move along the channel. In a box, the edges act as swarm traps and the trapping shows a discontinuous noise dependence. In a disk geometry, an ordered rotational state arises. This state is well described by a novel order parameter. These results provide the basis for applications of the Vicsek model to biological questions involving large particle numbers in confined environments. I. INTRODUCTION on a square simulation domain with periodic boundary conditions. Such a setup is useful for approximating bulk Swarms of fish, flocks of birds or herds of big mammals properties from finite size simulations. However, collec- are spectacular macroscopic phenomena that are famil- tive behaviour in many biological systems relies on a large iar to everybody. At first glance, they appear to be the number of individuals interacting in a confined environ- result of a complex interplay between highly developed ment. In this case, the effects of the geometry and the animals. But a deeper investigation reveals that the un- boundaries of the environment on the behaviour of a large derlying mechanisms could in fact be very simple. Bio- number of particles have to be examined by direct simu- logical systems with similar properties can also be found lations. at much smaller length scales including microorganisms Including local repulsive forces in the basic VM limits [1, 2] and even subcelullar components [3]. However, col- the maximal local density in the system. In this adap- lective behaviour is not limited to living organisms. Mi- tation of the VM, confining the particles by a circular cromotors as well as macroscopic and granular rods also reflective boundary results in rotations of the particles exhibit collective behaviour [4{6]. Surprisingly, systems both clockwise and anti-clockwise in the high density and with collective behaviour show effects that have previ- low noise regime [11]. A similar behaviour has been ob- ously only been known for thermodynamic equilibrium served for the continuous version of the VM [12]. Ex- systems [7]. For example, migrating tissue cells can un- tending the analysis to a channel with two parallel peri- dergo a phase transition from a disordered to an ordered odic and two reflective boundaries and a box with four state [8]. reflective boundaries reveals shape independent features A complete picture of many of these so called active of the collective behaviour. For all three geometries, the matter systems does not exist yet, but the results of boundaries introduce a spatial coherence of swarms in the ongoing research underline the universal properties and continuous VM that has not been observed in the basic the close connection of collective behaviour to statisti- version. Hence, in finite domains of reasonable size, the cal physics. The Vicsek model (VM) [9] is one of the boundaries may act as attractors. most prominent theoretical models to study active mat- Analysis of the collective behaviour of repulsive dis- ter. It is able to create complex macroscopic behaviour crete particles in the same three shapes shows a different like swarming or phase transitions [10] through a very picture [13]. Spatially incoherent transition states are ob- simple velocity alignment rule for the interaction of parti- served, where particles are ordered and densely packed in arXiv:2105.08792v2 [cond-mat.stat-mech] 21 Jun 2021 cles with their neighbours. It has been applied to a wide one part of the confined environment and the remaining range of systems and different variants including alter- system shows disorder. After a slow transition, phase ations to the interaction rule and an extension to a con- diagrams for order parameters adapted to the different tinuous model have been implemented [11]. Numerical geometries show that these states are eventually resolved simulations for the basic VM are implemented off-lattice into coherent behaviour as for periodic boundary condi- tions. This leads to the hypothesis that the boundaries hinder the spreading of the correlations. ∗ [email protected] These observations point to the importance of study- y sabine.fi[email protected] ing the effects of boundary conditions in the discrete VM. 2 Having the large particle numbers in biological systems counteracts the alignment and hence the order in the sys- in mind, we employ an on-lattice implementation. As tem. geometries, we consider a channel, a box and a disk with 2. Update of the position. Each particle is assigned a reflective boundary conditions and analyse their influence new position according to on the collective behaviour of the particles. We find that for a sufficiently small lattice, the behaviour of the VM cos Θi(t) xi(t + ∆t) = xi(t) + v0∆t : (4) is independent of the type of implementation. The re- sin Θi(t) flective boundary conditions result in particle alignment parallel to the tangential vector of the closest boundary point for sufficiently low noise levels in all geometries. For high levels of noise, the particle movement is uncor- related similar to the basic VM. The boundary alignment results in behaviour particular to each geometry. For the channel, separate swarms moving in opposite directions at the top and the bottom wall can occur. In the box ge- ometry, the boundary alignment can yield particle trap- ping in the corners. For the disk geometry, the boundary effect results in an ordered rotational state. Character- isation of the VM for these different geometries extends its applicability to a wider range of biological systems. II. MODEL AND ORDER PARAMETER A system in the basic VM [9] consists of N particles with two properties: position ~xi(t) and velocity ~vi(t). All FIG. 1: Alignment mechanism in the VM. The black particles move with the same absolute velocity v0. The dots mark the current position of the particles and the direction of movement is expressed by the angle Θi(t). Hence, the velocity is described by arrows point towards their current direction of movement. The i-th particle aligns its direction with all cos Θi(t) particles inside its interaction radius R (indicated by ~vi(t) = v0 (1) sin Θi(t) the blue circle and the blue arrows). All particles exist in a square shaped system (edge length L) with periodic boundary conditions. In the ini- In summary, the VM has five free parameters, particle tialisation phase, the system gets populated with N par- number N, system size L, velocity v0, interaction range ticles with randomly assigned positions and directions R and the noise strength η. These parameters are not of movement. Afterwards, the model evolves in discrete independent of each other and can be reduced to four time steps ∆t. During each time step, two updates are effective parameters. The behaviour of the VM only de- conducted for each particle: pends on the system size L, the noise strength η, the 1. Update of the particle direction. Each particle is NR2 density of interaction spheres ρ = L2 and the ratio of assigned a new direction v0 velocity and interaction range vr = R . To quantify the behaviour of the model, Vicsek et al. Θ (t + ∆t) = hΘ (t)i + ζ (t); (2) i j jxi−xj j<R i [9] have introduced the polar order parameter va, such that with N 1 X hsin Θj(t)i v = ~v : (5) hΘj(t)i = arctan ; (3) a Nv i hcos Θj(t)i 0 i=1 where hΘj(t)ijxi−xj j<R is the average direction of all The order parameter va is the average normalised ve- particles j surrounding the i-th particle (including itself) locity of the system. That means, if va = 1, all particles within a radius R (Fig. 1). The parameter R is the in- move in the same direction, and if va = 0 all particles teraction range of each particle. This alignment term move uncorrelated in the system. The VM shows a phase creates order in the system. transition, from an ordered motion state to a disordered The second term ζi(t) is a random angle drawn from motion state, with increasing η [14]. In large systems, a uniform probability distribution over (−η=2; η=2). The the state of the system changes from coherently moving range of the distribution η can be interpreted as a tem- swarms, to a band phase to a "cross sea" phase and fi- perature of the system. It is the main mechanism that nally to a completely uncorrelated state [15]. 3 To provide a good intuition on how the system is dis- cretized and how precise it is, v0, R and L are always given in units of lattice sites.

Load more

  12

Copy Link