Kinetic theory for systems of self-propelled particles with metric-free interactions

Anti-aligning interaction between active particles induces a finite wavelength instability: The dancing hexagons

By considering a simple model for self-propelled particle interaction, we show that anti-aligning forces induce a finite wavelength instability. Consequently, the system exhibits pattern formation. The formed pattern involves, let us say, a choreographic movement of the active entities. At the level of particle density, the system oscillates between a stripe pattern and a hexagonal one. The underlying dynamics of these density oscillations consists of two counterpropagating and purely hexagonal traveling waves. They are assembling and disassembling a global hexagonal structure and a striped lineup of particles. This self-assembling process becomes quite erratic for long-time simulations, seeming aperiodic.

Flocking in one dimension: Asters and reversals

We study the one-dimensional active Ising model in which aligning particles undergo diffusion biased by the signs of their spins. The phase diagram obtained varying the density of particles, their hopping rate, and the temperature controlling the alignment shows a homogeneous disordered phase but no homogeneous ordered one, as well as two phases with localized dense structures. In the flocking phase, large ordered aggregates move ballistically and stochastically reverse their direction of motion. In what we termed the "aster" phase, dense immobile aggregates of opposite magnetization face each other, exchanging particles, without any net motion of the aggregates. Using a combination of numerical simulations and mean-field theory, we study the evolution of the shapes of the flocks, the statistics of their reversal times, and their coarsening dynamics. Solving exactly for the zero-temperature dynamics of an aster allows us to understand their coarsening, which shows extremal dynamics, while mean-field equations account for their shape.

Selective interaction and its effect on collective motion

Plenty of empirical evidence on biological swarms reveal that interaction between individuals is selective. Each individual’s neighbor is selected based on one or more featured factors. Based on the self-propelled model, we develop a general probability neighbor selection framework to study the effect of four typical featured factors (i.e., distance, bearing, orientation change and bearing change). In this work, two common cases are involved to comprehensively analyze the impact of the four featured factors on the collective motion. One is the flocking, the other is the responsivity to stimulus. The impact of different selection strengths of the featured factors on both cases are investigated. The effect of noise on flocking and different stimulus intensities on responsivity to stimulus are analyzed. This study allows us to get the insight of selective interaction and suggests the potential solution to overcome the trade-off between flocking and responsivity quality.

Systematic Analysis of Emergent Collective Motion Produced by a 3D Hybrid Zonal Model

Emergent patterns of collective motion are thought to arise from local rules of interaction that govern how individuals adjust their velocity in response to the relative locations and velocities of near neighbours. Many models of collective motion apply rules of interaction over a metric scale, based on the distances to neighbouring group members. However, empirical work suggests that some species apply interactions over a topological scale, based on distance determined neighbour rank. Here, we modify an important metric model of collective motion (Couzin et al. in J Theor Biol 218(1):1-11, 2002), so that interactions relating to orienting movements with neighbours and attraction towards more distant neighbours operate over topological scales. We examine the emergent group movement patterns generated by the model as the numbers of neighbours that contribute to orientation- and attraction-based velocity adjustments vary. Like the metric form of the model, simulated groups can fragment (when interactions are influenced by less than 10-15% of the group), swarm and move in parallel, but milling does not occur. The model also generates other cohesive group movements including cases where groups exhibit directed motion without strong overall alignment of individuals. Multiple emergent states are possible for the same set of underlying model parameters in some cases, suggesting sensitivity to initial conditions, and there is evidence that emergent states of the system depend on the history of the system. Groups that do not fragment tend to stay relatively compact in terms of neighbour distances. Even if a group does fragment, individuals remain relatively close to near neighbours, avoiding complete isolation.

Fluctuation-Induced Phase Separation in Metric and Topological Models of Collective Motion

We study the role of noise on the nature of the transition to collective motion in dry active matter. Starting from field theories that predict a continuous transition at the deterministic level, we show that fluctuations induce a density-dependent shift of the onset of order, which in turn changes the nature of the transition into a phase-separation scenario. Our results apply to a range of systems, including models in which particles interact with their "topological" neighbors that have been believed so far to exhibit a continuous onset of order. Our analytical predictions are confirmed by numerical simulations of fluctuating hydrodynamics and microscopic models.

Flocking in complex environments-Attention trade-offs in collective information processing

The ability of biological and artificial collectives to outperform solitary individuals in a wide variety of tasks depends crucially on the efficient processing of social and environmental information at the level of the collective. Here, we model collective behavior in complex environments with many potentially distracting cues. Counter-intuitively, large-scale coordination in such environments can be maximized by strongly limiting the cognitive capacity of individuals, where due to self-organized dynamics the collective self-isolates from disrupting information. We observe a fundamental trade-off between coordination and collective responsiveness to environmental cues. Our results offer important insights into possible evolutionary trade-offs in collective behavior in biology and suggests novel principles for design of artificial swarms exploiting attentional bottlenecks.

Understanding how consensus is reached and information is processed within a collective is fundamental to many aspects of social dynamics in animals and humans. It is widely accepted that high connectivity among individuals facilitates group consensus, and being in a group provides benefits to individuals through social information about the environment provided by other group members. We show that this does not hold for collectives in complex environments: Limited attention capacity, that severely reduces connectivity among individuals, is highly beneficial for global coordination. However, this comes at a price: Collectives outperform isolated individuals in responding to the environment only at sufficiently high attention capacities, where global coordination breaks down. Thus, we demonstrate a fundamental trade-off in collective behavior between social coordination and responsiveness to environmental cues. Our work demonstrates the importance of sensory and cognitive limitations for the emergence and function of animal collectives, and poses fundamental questions about co-evolution of social behavior and individual attention capacity. The observed trade-off in collective information processing has implications for human social systems and for the design of robotic swarms operating in complex environments.

Transport coefficients of self-propelled particles: Reverse perturbations and transverse current correlations

The reverse perturbation method [Phys. Rev. E 59, 4894 (1999)1063-651X10.1103/PhysRevE.59.4894] for shearing simple liquids and measuring their viscosity is extended to the Vicsek model (VM) of active particles [Phys. Rev. Lett. 75, 1226 (1995)PRLTAO0031-900710.1103/PhysRevLett.75.1226] and its metric-free version. The sheared systems exhibit a phenomenon that is similar to the skin effect of an alternating electric current: Momentum that is fed into the boundaries of a layer decays mostly exponentially toward the center of the layer. It is shown how two transport coefficients, i.e., the shear viscosity ν and the momentum amplification coefficient λ, can be obtained by fitting this decay with an analytical solution of the hydrodynamic equations for the VM. The viscosity of the VM consists of two parts, a kinetic and a collisional contribution. While analytical predictions already exist for the former, a novel expression for the collisional part is derived by an Enskog-like kinetic theory. To verify the predictions for the transport coefficients, Green-Kubo relations were evaluated and transverse current correlations were measured in independent simulations. Not too far to the transition to collective motion, we find excellent agreement between the different measurements of the transport coefficients. However, the measured values of ν and 1-λ are always slightly higher than the mean-field predictions, even at large mean free paths and at state points quite far from the threshold to collective motion, that is, far in the disordered phase. These findings seem to indicate that the mean-field assumption of molecular chaos is much less reliable in systems with velocity-alignment rules such as the VM, compared to models obeying detailed balance such as multiparticle collision dynamics.

Contrarian compulsions produce exotic time-dependent flocking of active particles

Animals having a tendency to align their velocities to an average of those of their neighbors may flock as illustrated by the Vicsek model and its variants. If, in addition, they feel a systematic contrarian trend, the result may be a time periodic adjustment of the flock or period doubling in time. These exotic phases are predicted from kinetic theory and numerically found in a modified two-dimensional Vicsek model of self-propelled particles. Numerical simulations demonstrate striking effects of alignment noise on the polarization order parameter measuring particle flocking: maximum polarization length is achieved at an optimal nonzero noise level. When contrarian compulsions are more likely than conformist ones, nonuniform polarized phases appear as the noise surpasses threshold.

Giant Kovacs-Like Memory Effect for Active Particles

Dynamical properties of self-propelled particles obeying a bounded confidence rule are investigated by means of kinetic theory and agent-based simulations. While memory effects are observed in disordered systems, we show that they also occur in active matter systems. In particular, we find that the system exhibits a giant Kovacs-like memory effect that is much larger than predicted by a generic linear theory. Based on a separation of time scales we develop a nonlinear theory to explain this effect. We apply this theory to driven granular gases and propose further applications to spin glasses.

Gaussian theory for spatially distributed self-propelled particles

Obtaining a reduced description with particle and momentum flux densities outgoing from the microscopic equations of motion of the particles requires approximations. The usual method, we refer to as truncation method, is to zero Fourier modes of the orientation distribution starting from a given number. Here we propose another method to derive continuum equations for interacting self-propelled particles. The derivation is based on a Gaussian approximation (GA) of the distribution of the direction of particles. First, by means of simulation of the microscopic model, we justify that the distribution of individual directions fits well to a wrapped Gaussian distribution. Second, we numerically integrate the continuum equations derived in the GA in order to compare with results of simulations. We obtain that the global polarization in the GA exhibits a hysteresis in dependence on the noise intensity. It shows qualitatively the same behavior as we find in particles simulations. Moreover, both global polarizations agree perfectly for low noise intensities. The spatiotemporal structures of the GA are also in agreement with simulations. We conclude that the GA shows qualitative agreement for a wide range of noise intensities. In particular, for low noise intensities the agreement with simulations is better as other approximations, making the GA to an acceptable candidates of describing spatially distributed self-propelled particles.

On effective temperature in network models of collective behavior

Collective behavior of self-propelled units is studied analytically within the Vectorial Network Model (VNM), a mean-field approximation of the well-known Vicsek model. We propose a dynamical systems framework to study the stochastic dynamics of the VNM in the presence of general additive noise. We establish that a single parameter, which is a linear function of the circular mean of the noise, controls the macroscopic phase of the system-ordered or disordered. By establishing a fluctuation-dissipation relation, we posit that this parameter can be regarded as an effective temperature of collective behavior. The exact critical temperature is obtained analytically for systems with small connectivity, equivalent to low-density ensembles of self-propelled units. Numerical simulations are conducted to demonstrate the applicability of this new notion of effective temperature to the Vicsek model. The identification of an effective temperature of collective behavior is an important step toward understanding order-disorder phase transitions, informing consistent coarse-graining techniques and explaining the physics underlying the emergence of collective phenomena.

The evolution of distributed sensing and collective computation in animal populations

Many animal groups exhibit rapid, coordinated collective motion. Yet, the evolutionary forces that cause such collective responses to evolve are poorly understood. Here, we develop analytical methods and evolutionary simulations based on experimental data from schooling fish. We use these methods to investigate how populations evolve within unpredictable, time-varying resource environments. We show that populations evolve toward a distinctive regime in behavioral phenotype space, where small responses of individuals to local environmental cues cause spontaneous changes in the collective state of groups. These changes resemble phase transitions in physical systems. Through these transitions, individuals evolve the emergent capacity to sense and respond to resource gradients (i.e. individuals perceive gradients via social interactions, rather than sensing gradients directly), and to allocate themselves among distinct, distant resource patches. Our results yield new insight into how natural selection, acting on selfish individuals, results in the highly effective collective responses evident in nature.

DOI:http://dx.doi.org/10.7554/eLife.10955.001

In nature, we see many examples of highly coordinated movements of groups of individuals; think of a flock of birds turning swiftly in unison or a crowd of people filing through the exit of a building. A common feature of these behaviors is that they occur without any centralized control, and that they involve sudden and often dramatic changes in the 'collective state' of the group (i.e. speed, or the distances between individuals). In the past, researchers have likened these transitions in collective behavior to phase transitions in physical systems, for example, the transition between liquid water and water vapor. However, it is not clear how such collective responses could have evolved.

Natural selection is an evolutionary process whereby individuals with particularly 'fit' traits produce more offspring than others. Over many generations, these beneficial traits tend to become more common in the population. Hein, Rosenthal, Hagstrom et al. developed a mathematical model to investigate whether the capacity of a population to perform collective motions could evolve through natural selection.

The model shows that over many generations, populations consistently evolve a unique collective trait whereby small responses of individuals to an environmental cue can cause spontaneous changes in the collective state of the local population. These transitions in collective state greatly enhance the ability of individuals to locate and exploit resources. Hein, Rosenthal, Hagstrom et al.’s findings suggest that natural selection acting on the behavior of individuals can cause a population to evolve a distinctive, collective behavior.

The next challenge will be to identify a biological system in which the evolution of collective motion can be studied experimentally to test these predictions.

DOI:http://dx.doi.org/10.7554/eLife.10955.002

Active matter beyond mean-field: ring-kinetic theory for self-propelled particles

Recently, Hanke et al. [Phys. Rev. E 88, 052309 (2013)] showed that mean-field kinetic theory fails to describe collective motion in soft active colloids and that correlations must not be neglected. Correlation effects are also expected to be essential in systems of biofilaments driven by molecular motors and in swarms of midges. To obtain correlations in an active matter system from first principles, we derive a ring-kinetic theory for Vicsek-style models of self-propelled agents from the exact N-particle evolution equation in phase space. The theory goes beyond mean-field and does not rely on Boltzmann's approximation of molecular chaos. It can handle precollisional correlations and cluster formation, which are both important to understand the phase transition to collective motion. We propose a diagrammatic technique to perform a small-density expansion of the collision operator and derive the first two equations of the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. An algorithm is presented that numerically solves the evolution equation for the two-particle correlations on a lattice. Agent-based simulations are performed and informative quantities such as orientational and density correlation functions are compared with those obtained by ring-kinetic theory. Excellent quantitative agreement between simulations and theory is found at not-too-small noises and mean free paths. This shows that there are parameter ranges in Vicsek-like models where the correlated closure of the BBGKY hierarchy gives correct and nontrivial results. We calculate the dependence of the orientational correlations on distance in the disordered phase and find that it seems to be consistent with a power law with an exponent around -1.8, followed by an exponential decay. General limitations of the kinetic theory and its numerical solution are discussed.

Tricritical points in a Vicsek model of self-propelled particles with bounded confidence

We study the orientational ordering in systems of self-propelled particles with selective interactions. To introduce the selectivity we augment the standard Vicsek model with a bounded-confidence collision rule: a given particle only aligns to neighbors who have directions quite similar to its own. Neighbors whose directions deviate more than a fixed restriction angle α are ignored. The collective dynamics of this system is studied by agent-based simulations and kinetic mean-field theory. We demonstrate that the reduction of the restriction angle leads to a critical noise amplitude decreasing monotonically with that angle, turning into a power law with exponent 3/2 for small angles. Moreover, for small system sizes we show that upon decreasing the restriction angle, the kind of the transition to polar collective motion changes from continuous to discontinuous. Thus, an apparent tricritical point with different scaling laws is identified and calculated analytically. We investigate the shifting and vanishing of this point due to the formation of density bands as the system size is increased. Agent-based simulations in small systems with large particle velocities show excellent agreement with the kinetic theory predictions. We also find that at very small interaction angles, the polar ordered phase becomes unstable with respect to the apolar phase. We derive analytical expressions for the dependence of the threshold noise on the restriction angle. We show that the mean-field kinetic theory also permits stationary nematic states below a restriction angle of 0.681π. We calculate the critical noise, at which the disordered state bifurcates to a nematic state, and find that it is always smaller than the threshold noise for the transition from disorder to polar order. The disordered-nematic transition features two tricritical points: At low and high restriction angle, the transition is discontinuous but continuous at intermediate α. We generalize our results to systems that show fragmentation into more than two groups and obtain scaling laws for the transition lines and the corresponding tricritical points. A numerical method to evaluate the nonlinear Fredholm integral equation for the stationary distribution function is also presented. This method is shown to give excellent agreement with agent-based simulations, even in strongly ordered systems at noise values close to zero.

Invasion-wave-induced first-order phase transition in systems of active particles

An instability near the transition to collective motion of self-propelled particles is studied numerically by Enskog-like kinetic theory. While hydrodynamics breaks down, the kinetic approach leads to steep solitonlike waves. These supersonic waves show hysteresis and lead to an abrupt jump of the global order parameter if the noise level is changed. Thus they provide a mean-field mechanism to change the second-order character of the phase transition to first order. The shape of the wave is shown to follow a scaling law and to quantitatively agree with agent-based simulations.

Revisiting the flocking transition using active spins

We consider an active Ising model in which spins both diffuse and align on lattice in one and two dimensions. The diffusion is biased so that plus or minus spins hop preferably to the left or to the right, which generates a flocking transition at low temperature and high density. We construct a coarse-grained description of the model that predicts this transition to be a first-order liquid-gas transition in the temperature-density ensemble, with a critical density sent to infinity. In this first-order phase transition, the magnetization is proportional to the liquid fraction and thus varies continuously throughout the phase diagram. Using microscopic simulations, we show that this theoretical prediction holds in 2D whereas the fluctuations alter the transition in 1D, preventing, for instance, any spontaneous symmetry breaking.