Fluctuations and pattern formation in self-propelled particles

Hysteresis, reentrance, and glassy dynamics in systems of self-propelled rods

Nonequilibrium active matter made up of self-driven particles with short-range repulsive interactions is a useful minimal system to study active matter as the system exhibits collective motion and nonequilibrium order-disorder transitions. We studied high-aspect-ratio self-propelled rods over a wide range of packing fractions and driving to determine the nonequilibrium state diagram and dynamic properties. Flocking and nematic-laning states occupy much of the parameter space. In the flocking state, the average internal pressure is high and structural and mechanical relaxation times are long, suggesting that rods in flocks are in a translating glassy state despite overall flock motion. In contrast, the nematic-laning state shows fluidlike behavior. The flocking state occupies regions of the state diagram at both low and high packing fraction separated by nematic-laning at low driving and a history-dependent region at higher driving; the nematic-laning state transitions to the flocking state for both compression and expansion. We propose that the laning-flocking transitions are a type of glass transition that, in contrast to other glass-forming systems, can show fluidization as density increases. The fluid internal dynamics and ballistic transport of the nematic-laning state may promote collective dynamics of rod-shaped micro-organisms.

Hydrodynamics of Turning Flocks

We present a hydrodynamic model of flocking that generalizes the familiar Toner-Tu equations to incorporate turning inertia of well-polarized flocks. The continuum equations controlled by only two dimensionless parameters, orientational inertia and alignment strength, are derived by coarse-graining the inertial spin model recently proposed by Cavagna et al. The interplay between orientational inertia and bend elasticity of the flock yields anisotropic spin waves that mediate the propagation of turning information throughout the flock. The coupling between spin-current density to the local vorticity field through a nonlinear friction gives rise to a hydrodynamic mode with angular-dependent propagation speed at long wavelengths. This mode becomes unstable as a result of the growth of bend and splay deformations augmented by the spin wave, signaling the transition to complex spatiotemporal patterns of continuously turning and swirling flocks.

Pattern formation in flocking models: A hydrodynamic description

We study in detail the hydrodynamic theories describing the transition to collective motion in polar active matter, exemplified by the Vicsek and active Ising models. Using a simple phenomenological theory, we show the existence of an infinity of propagative solutions, describing both phase and microphase separation, that we fully characterize. We also show that the same results hold specifically in the hydrodynamic equations derived in the literature for the active Ising model and for a simplified version of the Vicsek model. We then study numerically the linear stability of these solutions. We show that stable ones constitute only a small fraction of them, which, however, includes all existing types. We further argue that, in practice, a coarsening mechanism leads towards phase-separated solutions. Finally, we construct the phase diagrams of the hydrodynamic equations proposed to qualitatively describe the Vicsek and active Ising models and connect our results to the phenomenology of the corresponding microscopic models.

Flocking with discrete symmetry: The two-dimensional active Ising model

We study in detail the active Ising model, a stochastic lattice gas where collective motion emerges from the spontaneous breaking of a discrete symmetry. On a two-dimensional lattice, active particles undergo a diffusion biased in one of two possible directions (left and right) and align ferromagnetically their direction of motion, hence yielding a minimal flocking model with discrete rotational symmetry. We show that the transition to collective motion amounts in this model to a bona fide liquid-gas phase transition in the canonical ensemble. The phase diagram in the density-velocity parameter plane has a critical point at zero velocity which belongs to the Ising universality class. In the density-temperature "canonical" ensemble, the usual critical point of the equilibrium liquid-gas transition is sent to infinite density because the different symmetries between liquid and gas phases preclude a supercritical region. We build a continuum theory which reproduces qualitatively the behavior of the microscopic model. In particular, we predict analytically the shapes of the phase diagrams in the vicinity of the critical points, the binodal and spinodal densities at coexistence, and the speeds and shapes of the phase-separated profiles.

Kinetic attractor phase diagrams of active nematic suspensions: the dilute regime

Large-scale simulations by the authors of the kinetic-hydrodynamic equations for active polar nematics revealed a variety of spatio-temporal attractors, including steady and unsteady, banded (1d) and cellular (2d) spatial patterns. These particle scale activation-induced attractors arise at dilute nanorod volume fractions where the passive equilibrium phase is isotropic, whereas all previous model simulations have focused on the semi-dilute, nematic equilibrium regime and mostly on low-moment orientation tensor and polarity vector models. Here we extend our previous results to complete attractor phase diagrams for active nematics, with and without an explicit polar potential, to map out novel spatial and dynamic transitions, and to identify some new attractors, over the parameter space of dilute nanorod volume fraction and nanorod activation strength. The particle-scale activation parameter corresponds experimentally to a tunable force dipole strength (so-called pushers with propulsion from the rod tail) generated by active rod macromolecules, e.g., catalysis with the solvent phase, ATP-induced propulsion, or light-activated propulsion. The simulations allow 2d spatial variations in all flow and orientational variables and full spherical orientational degrees of freedom; the attractors correspond to numerical integration of a coupled system of 125 nonlinear PDEs in 2d plus time. The phase diagrams with and without the polar interaction potential are remarkably similar, implying that polar interactions among the rodlike particles are not essential to long-range spatial and temporal correlations in flow, polarity, and nematic order. As a general rule, above a threshold, low volume fractions induce 1d banded patterns, whereas higher yet still dilute volume fractions yield 2d patterns. Again as a general rule, varying activation strength at fixed volume fraction induces novel dynamic transitions. First, stationary patterns saturate the instability of the isotropic state, consisting of discrete 1d banded or 2d cellular patterns depending on nanorod volume fraction. Increasing activation strength further induces a sequence of attractor bifurcations, including oscillations superimposed on the 1d and 2d stationary patterns, a uniform translational motion of 1d and 2d oscillating patterns, and periodic switching between 1d and 2d patterns. These results imply that active macromolecular suspensions are capable of long-range spatial and dynamic organization at isotropic equilibrium concentrations, provided particle-scale activation is sufficiently strong.

Collective motion of self-propelled particles with memory

We show that memory, in the form of underdamped angular dynamics, is a crucial ingredient for the collective properties of self-propelled particles. Using Vicsek-style models with an Ornstein-Uhlenbeck process acting on angular velocity, we uncover a rich variety of collective phases not observed in usual overdamped systems, including vortex lattices and active foams. In a model with strictly nematic interactions the smectic arrangement of Vicsek waves giving rise to global polar order is observed. We also provide a calculation of the effective interaction between vortices in the case where a telegraphic noise process is at play, explaining thus the emergence and structure of the vortex lattices observed here and in motility assay experiments.

From phase to microphase separation in flocking models: the essential role of nonequilibrium fluctuations

We show that the flocking transition in the Vicsek model is best understood as a liquid-gas transition, rather than an order-disorder one. The full phase separation observed in flocking models with Z(2) rotational symmetry is, however, replaced by a microphase separation leading to a smectic arrangement of traveling ordered bands. Remarkably, continuous deterministic descriptions do not account for this difference, which is only recovered at the fluctuating hydrodynamics level. Scalar and vectorial order parameters indeed produce different types of number fluctuations, which we show to be essential in selecting the inhomogeneous patterns. This highlights an unexpected role of fluctuations in the selection of flock shapes.

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.

Collision dynamics of traveling bands in systems of deformable self-propelled particles

We study the soliton-like character of traveling bands in systems of interacting deformable self-propelled particles in two dimensions. The collision dynamics of the model in which migration velocity increases with increasing local density is investigated numerically by changing the relaxation rate of deformations. The bands are unstable upon head-on collisions for larger relaxation rates. This clearly indicates that deformability plays a crucial role of the soliton-like behavior.

Capillary instability of axisymmetric, active liquid crystal jets

We study linear stability of an infinitely long, axisymmetric, cylindrical active liquid crystal (ALC) jet in a passive isotropic fluid matrix using a polar active liquid crystal (ALC) model. We identify three possible unstable modes (or mechanisms) as the result of the interaction between the flow and the active (or self-propelled) molecular motion. The first unstable mode is related to the polarity vector instability when coupled to the flow field in the presence of the molecular activity. It can be traced back to the inherent polarity vector instability in a bulk active liquid crystal flow. However, it can be grossly amplified in the ALC jet to encompass up to infinitely many unstable growth rates when the long range distortional elastic interaction is weak in certain parameter regimes; it can also be suppressed in other parameter regimes completely. The second unstable mode is related to the classical capillary or Rayleigh instability, which exists in a finite wave interval [0, k(cutoff)]. The new feature for this instability lies in the dependence of the cutoff wave number (k(cutoff)) on the activity of the active matter system. For ALC jets with sufficiently strong contractile activity, the instability can be completely suppressed though. The third unstable mode is due to the active viscous stress. This unstable mode can emerge in the intermediate wave number regime at a sufficiently strong active viscosity and even expand all the way to the zero wave number limit when the Rayleigh unstable mode is absent. It can also be suppressed in the regime of weak active viscous stress. At any given values of the model parameters, the three types of instabilities can show up either individually or in a certain combination, or be completely suppressed altogether. In this paper, we discuss the positive growth rates associated with the instabilities, windows of instability and their dependence on model parameters through extensive numerical computations aided by asymptotic analyses.

Large-scale chaos and fluctuations in active nematics

We show that dry active nematics, e.g., collections of shaken elongated granular particles, exhibit large-scale spatiotemporal chaos made of interacting dense, ordered, bandlike structures in a parameter region including the linear onset of nematic order. These results are obtained from the study of both the well-known (deterministic) hydrodynamic equations describing these systems and of the self-propelled particle model they were derived from. We prove, in particular, that the chaos stems from the generic instability of the band solution of the hydrodynamic equations. Revisiting the status of the strong fluctuations and long-range correlations in the particle model, we show that the giant number fluctuations observed in the chaotic phase are a trivial consequence of density segregation. However anomalous, curvature-driven number fluctuations are present in the homogeneous quasiordered nematic phase and characterized by a nontrivial scaling exponent.

Activating membranes

We present a general dynamical theory of a membrane coupled to an actin cortex containing polymerizing filaments with active stresses and currents, and demonstrate that active membrane dynamics [S. Ramaswamy et al., Phys. Rev. Lett. 84, 3494 (2000)] and spontaneous shape oscillations emerge from this description. We also consider membrane instabilities and patterns induced by the presence of filaments with polar orientational correlations in the tangent plane of the membrane. The dynamical features we predict should be seen in a variety of cellular contexts involving the dynamics of the membrane-cytoskeleton composite and cytoskeletal extracts coupled to synthetic vesicles.

Macroscopic traveling packet and soliton states of quasi-one-dimensional flocks

Using a continuum model for inhomogeneous flocks, we show that a finite but arbitrarily large moving "packet" of active particles (e.g., moving creatures) can form in a background of a lower density disordered phase of these particles, like a liquid drop surrounded by vapor. The "vapor density" of the disordered background can be made arbitrarily low. We find three basic types of quasi-one-dimensional states: "longitudinal", "transverse", and "oblique" states, with their internal velocity fields, respectively, parallel, perpendicular, and oblique to the interface. The transitions between these states are also studied.

Emergent spatial structures in flocking models: a dynamical system insight

We show that hydrodynamic theories of polar active matter generically possess inhomogeneous traveling solutions. We introduce a unifying dynamical-system framework to establish the shape of these intrinsically nonlinear patterns, and show that they correspond to those hitherto observed in experiments and numerical simulation: periodic density waves, and solitonic bands, or polar-liquid droplets both cruising in isotropic phases. We elucidate their respective multiplicity and mutual relations, as well as their existence domain.

Dynamics of self-organized rotating spiral-coils in bacterial swarms

Self-propelled particles (SPP) exhibit complex collective motions, mimicking autonomous behaviors that are often seen in the natural world, but essentially are generated by simple mutual interactions. Previous research on SPP systems focuses on collective behaviors of a uniform population. However, very little is known about the evolution of individual particles under the same global influence. Here we show self-organized rotating spiral coils in a two-dimensional (2D) active system. By using swarming bacteria Vibrio alginolyticus as an ideal experimental realization of a well-controlled 2D self-propelled system, we study the interaction between ultra-long cells and short background active cells. The self-propulsion of long cells and their interactions with neighboring short cells leads to a self-organized, stable spiral rotational state in 2D. We find four types of spiral coils with two main features: the rotating direction (clockwise or counter-clockwise) and the central structure (single or double spiral). The body length of the spiral coils falls between 32 and 296 μm and their rotational speed is within a range from 2.22 to 22.96 rad s(-1). The dynamics of these spiral coils involves folding and unfolding processes, which require local velocity changes of the long bacterium. This phenomenon can be qualitatively replicated by a Brownian dynamics simulation using a simple rule of the propulsion thrust, imitating the reorientation of bacterial flagella. Apart from the physical and biological interests in swarming cells, the formation of self-organized spiral coils could be useful for the next generation of microfabrication.

Active crystals and their stability

A recently introduced active phase field crystal model describes the formation of ordered resting and traveling crystals in systems of self-propelled particles. Increasing the active drive, a resting crystal can be forced to perform collectively ordered migration as a single traveling object. We demonstrate here that these ordered migrating structures are linearly stable. In other words, during migration, the single-crystalline texture together with the globally ordered collective motion is preserved even on large length scales. Furthermore, we consider self-propelled particles on a substrate that are surrounded by a thin fluid film. We find that in this case the resulting hydrodynamic interactions can destabilize the order.

Self-organized vortices of circling self-propelled particles and curved active flagella

Self-propelled pointlike particles move along circular trajectories when their translocation velocity is constant and the angular velocity related to their orientation vector is also constant. We investigate the collective behavior of ensembles of such circle swimmers by Brownian dynamics simulations. If the particles interact via a "velocity-trajectory coordination" rule within neighboring particles, a self-organized vortex pattern emerges. This vortex pattern is characterized by its particle-density correlation function Gρ, the density correlation function Gc of trajectory centers, and an order parameter S representing the degree of the aggregation of the particles. Here we systematically vary the system parameters, such as the particle density and the interaction range, in order to reveal the transition of the system from a light-vortex-dominated to heavy-vortex-dominated state, where vortices contain mainly a single and many self-propelled particles, respectively. We also study a semidilute solution of curved, sinusoidal-beating flagella, as an example of circling self-propelled particles with explicit propulsion mechanism and excluded-volume interactions. Our simulation results are compared with previous experimental results for the vortices in sea-urchin sperm solutions near a wall. The properties of the vortices in simulations and experiments are found to agree quantitatively.

Spiral and never-settling patterns in active systems

We present a combined numerical and analytical study of pattern formation in an active system where particles align, possess a density-dependent motility, and are subject to a logistic reaction. The model can describe suspensions of reproducing bacteria, as well as polymerizing actomyosin gels in vitro or in vivo. In the disordered phase, we find that motility suppression and growth compete to yield stable or blinking patterns, which, when dense enough, acquire internal orientational ordering to give asters or spirals. We predict these may be observed within chemotactic aggregates in bacterial fluids. In the ordered phase, the reaction term leads to previously unobserved never-settling patterns which can provide a simple framework to understand the formation of motile and spiral patterns in intracellular actin systems.

Formation and collision of traveling bands in interacting deformable self-propelled particles

We study the collective dynamics of interacting deformable self-propelled particles whose migration velocity increases with increasing local density. In two-dimensional numerical simulations of this system, the local density dependence on migration velocity leads to traveling bands similar to those previously reported for Vicsek-type models. We show that a pair of straight bands moving in opposite directions survives a head-on collision. Although traveling bands also appear in systems of constant migration velocity subjected to random noise, they are found to be unstable in a head-on collision.

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.

Topological defects and density fluctuations in collectively moving systems

Ensembles of collectively moving particles like flocks of birds, bacteria, or filamentous polymers show a broad range of intriguing phenomena, yet seem to obey very similar physical principles. These generic principles have been predicted to lead to characteristic density fluctuations, which are in sharp contrast to normal fluctuations determining the properties of ordered systems in thermal equilibrium. Using high-density motility assays of driven filaments, we characterize here the origin and nature of giant fluctuations that emerge in this class of systems. By showing that these unique statistical properties result from the coupling between particle density and the topology of the velocity field of the particles, we provide insight in the physics of collective motion.

Traveling and resting crystals in active systems

A microscopic field theory for crystallization in active systems is proposed which unifies the phase-field-crystal model of freezing with the Toner-Tu theory for self-propelled particles. A wealth of different active crystalline states are predicted and characterized. In particular, for increasing strength of self-propulsion, a transition from a resting crystal to a traveling crystalline state is found where the particles migrate collectively while keeping their crystalline order. Our predictions, which are verifiable in experiments and in particle-resolved computer simulations, provide a starting point for the design of new active materials.

Nonlinear field equations for aligning self-propelled rods

We derive a set of minimal and well-behaved nonlinear field equations describing the collective properties of self-propelled rods from a simple microscopic starting point, the Vicsek model with nematic alignment. Analysis of their linear and nonlinear dynamics shows good agreement with the original microscopic model. In particular, we derive an explicit expression for density-segregated, banded solutions, allowing us to develop a more complete analytic picture of the problem at the nonlinear level.

Emergent states in dense systems of active rods: from swarming to turbulence

Dense suspensions of self-propelled rod-like particles exhibit a fascinating variety of non-equilibrium phenomena. By means of computer simulations of a minimal model for rigid self-propelled colloidal rods with variable shape we explore the generic diagram of emerging states over a large range of rod densities and aspect ratios. The dynamics is studied using a simple numerical scheme for the overdamped noiseless frictional dynamics of a many-body system in which steric forces are dominant over hydrodynamic ones. The different emergent states are identified by various characteristic correlation functions and suitable order parameter fields. At low density and aspect ratio, a disordered phase with no coherent motion precedes a highly cooperative swarming state with giant number fluctuations at large aspect ratio. Conversely, at high densities weakly anisometric particles show a distinct jamming transition whereas slender particles form dynamic laning patterns. In between there is a large window corresponding to strongly vortical, turbulent flow. The different dynamical states should be verifiable in systems of swimming bacteria and artificial rod-like micro-swimmers.

Nucleation-induced transition to collective motion in active systems

While the existence of polar ordered states in active systems is well established, the dynamics of the self-assembly processes are still elusive. We study a lattice gas model of self-propelled elongated particles interacting through excluded volume and alignment interactions, which shows a phase transition from an isotropic to a polar ordered state. By analyzing the ordering process we find that the transition is driven by the formation of a critical nucleation cluster and a subsequent coarsening process. Moreover, the time to establish a polar ordered state shows a power-law divergence.

Self-regulation in self-propelled nematic fluids

We consider the hydrodynamic theory of an active fluid of self-propelled particles with nematic aligning interactions. This class of materials has polar symmetry at the microscopic level, but forms macrostates of nematic symmetry. We highlight three key features of the dynamics. First, as in polar active fluids, the control parameter for the order-disorder transition, namely the density, is dynamically convected by the order parameter via active currents. The resulting dynamical self-regulation of the order parameter is a generic property of active fluids and destabilizes the uniform nematic state near the mean-field transition. Secondly, curvature-driven currents render the system unstable deep in the nematic state, as found previously. Finally, and unique to self-propelled nematics, nematic order induces local polar order that in turn leads to the growth of density fluctuations. We propose this as a possible mechanism for the smectic order of polar clusters seen in numerical simulations.

Reanalysis of the hydrodynamic theory of fluid, polar-ordered flocks

I reanalyze the hydrodynamic theory of fluid, polar-ordered flocks. I find new linear terms in the hydrodynamic equations which slightly modify the anisotropy, but not the scaling, of the damping of sound modes. I also find that the nonlinearities allowed in equilibrium do not stabilize long-ranged order in spatial dimensions d=2, in accord with the Mermin-Wagner theorem. Nonequilibrium nonlinearities do stabilize long-ranged order in d=2, as argued by earlier work. Some of these were missed by earlier work; it is unclear whether or not they change the scaling exponents in d=2.

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

A model of self-driven particles similar to the Vicsek model [Phys. Rev. Lett. 75, 1226 (1995)] but with metric-free interactions is studied by means of a novel Enskog-type kinetic theory. In this model, N particles of constant speed v(0) try to align their travel directions with the average direction of a fixed number of closest neighbors. At strong alignment a global flocking state forms. The alignment is defined by a stochastic rule, not by a Hamiltonian. The corresponding interactions are of genuine multibody nature. The theory is based on a Master equation in 3N-dimensional phase space, which is made tractable by means of the molecular chaos approximation. The phase diagram for the transition to collective motion is calculated and compared to direct numerical simulations. A linear stability analysis of a homogeneous ordered state is performed using the kinetic but not the hydrodynamic equations in order to achieve high accuracy. In contrast to the regular metric Vicsek-model no instabilities occur. This confirms previous direct simulations that, for Vicsek-like models with metric-free interactions, there is no formation of density bands and that the flocking transition is continuous.

Thermal and athermal three-dimensional swarms of self-propelled particles

Swarms of self-propelled particles exhibit complex behavior that can arise from simple models, with large changes in swarm behavior resulting from small changes in model parameters. We investigate the steady-state swarms formed by self-propelled Morse particles in three dimensions using molecular dynamics simulations optimized for graphics processing units. We find a variety of swarms of different overall shape assemble spontaneously and that for certain Morse potential parameters at most two competing structures are observed. We report a rich "phase diagram" of athermal swarm structures observed across a broad range of interaction parameters. Unlike the structures formed in equilibrium self-assembly, we find that the probability of forming a self-propelled swarm can be biased by the choice of initial conditions. We investigate how thermal noise influences swarm formation and demonstrate ways it can be exploited to reconfigure one swarm into another. Our findings validate and extend previous observations of self-propelled Morse swarms and highlight open questions for predictive theories of nonequilibrium self-assembly.

Collective dynamics of self-propelled particles with variable speed

Understanding the organization of collective motion in biological systems is an ongoing challenge. In this paper we consider a minimal model of self-propelled particles with variable speed. Inspired by experimental data from schooling fish, we introduce a power-law dependency of the speed of each particle on the degree of polarization order in its neighborhood. We derive analytically a coarse-grained continuous approximation for this model and find that, while the specific variable speed rule used does not change the details of the ordering transition leading to collective motion, it induces an inverse power-law correlation between the speed or the local polarization order and the local density. Using numerical simulations, we verify the range of validity of this continuous description and explore regimes beyond it. We discover, in disordered states close to the transition, a phase-segregated regime where most particles cluster into almost static groups surrounded by isolated high-speed particles. We argue that the mechanism responsible for this regime could be present in a wide range of collective motion dynamics.

Pattern formation in self-propelled particles with density-dependent motility

We study the behavior of interacting self-propelled particles, whose self-propulsion speed decreases with their local density. By combining direct simulations of the microscopic model with an analysis of the hydrodynamic equations obtained by explicitly coarse graining the model, we show that interactions lead generically to the formation of a host of patterns, including moving clumps, active lanes, and asters. This general mechanism could explain many of the patterns seen in recent experiments and simulations.

Spatial structures and giant number fluctuations in models of active matter

The large scale fluctuations of the ordered state in active matter systems are usually characterized by studying the "giant number fluctuations" of particles in any finite volume, as compared to the expectations from the central limit theorem. However, in ordering systems, the fluctuations in density ordering are often captured through their structure functions deviating from Porod's law. In this Letter we study the relationship between giant number fluctuations and structure functions for different models of active matter as well as other nonequilibrium systems. A unified picture emerges, with different models falling in four distinct classes depending on the nature of their structure functions. For one class, we show that experimentalists may find Porod's law violation, by measuring subleading corrections to the number fluctuations.

Dynamical self-regulation in self-propelled particle flows

We study a continuum model of overdamped self-propelled particles with aligning interactions in two dimensions. Combining analytical theory and computations, we map out the phase diagram for the parameter space covered by the model. We find that the system self-organizes into two robust structures in different regions of parameter space: solitary waves composed of ordered swarms moving through a low density disordered background, and stationary radially symmetric asters. The self-regulating nature of the flow yields phase separation, ubiquitous in this class of systems, and controls the formation of solitary waves. Self-propulsion and the associated active convection play a crucial role in aster formation.

Effect of dynamic and static friction on an asymmetric granular piston

We investigate the influence of dry friction on an asymmetric, granular piston of mass M, composed of two materials, undergoing inelastic collisions with bath particles of mass m. Numerical simulations of the Boltzmann-Lorentz equation reveal the existence of two scaling regimes depending on the friction strength. In the large friction limit, we introduce an exact model giving the asymptotic behavior of the Boltzmann-Lorentz equation. For small friction and for large mass ratio M/m, we derive a Fokker-Planck equation for which the exact solution is also obtained. Static friction attenuates the motor effect and results in a discontinuous velocity distribution.

Collective motion of binary self-propelled particle mixtures

In this study, we investigate the phenomenon of collective motion in binary mixtures of self-propelled particles. More precisely, we consider two particle species, each of which consisting of pointlike objects that propel with a velocity of constant magnitude. Within each species, the particles try to achieve polar alignment of their velocity vectors, whereas we analyze the cases of preferred polar, antiparallel, as well as perpendicular alignment between particles of different species. Our focus is on the effect that the interplay between the two species has on the threshold densities for the onset of collective motion and on the nature of the solutions above onset. For this purpose, we start from suitable Langevin equations in the particle picture, from which we derive mean field equations of the Fokker-Planck type and finally macroscopic continuum field equations. We perform particle simulations of the Langevin equations and linear stability analyses of the Fokker-Planck and macroscopic continuum equations, and we numerically solve the Fokker-Planck equations. Both spatially homogeneous and inhomogeneous solutions are investigated, where the latter correspond to stripelike flocks of collectively moving particles. In general, the interaction between the two species reduces the threshold density for the onset of collective motion of each species. However, this interaction also reduces the spatial organization in the stripelike flocks. The case that shows the most interesting behavior is the one of preferred perpendicular alignment between different species. There a competition between polar and truly nematic orientational ordering of the velocity vectors takes place within each particle species. Finally, depending on the alignment rule for particles of different species and within certain ranges of particle densities, identical and inverted spatial density profiles can be found for the two particle species. The system under investigation is confined to two spatial dimensions.

Fast wandering of slow birds

I study a single slow bird moving with a flock of birds of a different and faster (or slower) species. I find that every species of flocker has a characteristic speed γ ≠ v(0), where v(0) is the mean speed of the flock such that if the speed v(s) of the slow bird equals γ, it will randomly wander transverse to the mean direction of flock motion far faster than the other birds will: Its mean-squared transverse displacement will grow in d = 2 with time t like t(5/3), in contrast to t(4/3) for the other birds. In d = 3, the slow bird's mean-squared transverse displacement grows like t(5/4), in contrast to t for the other birds. If v(s) ≠ γ, the mean-squared displacement of the slow bird crosses over from t(5/3) to t(4/3) scaling in d = 2 and from t(5/4) to t scaling in d = 3 at a time t(c) that scales according to t(c) proportionally |v(s) - γ|(-2).

Shearing self-propelled suspensions: arrest of coarsening and suppression of giant density fluctuations

We study the effect of a linear shear flow on a collection of interacting active, self-propelled particles modeled via the Vicsek model. The imposed flow has a dramatic effect on the behavior of the model. We find that in the presence of shear there is no order-disorder transition, and that coarsening of the domains is arrested. Shear also suppresses the so-called giant density fluctuations that are observed in the quiescent limit.

Emergence of coherent structures and large-scale flows in motile suspensions

The emergence of coherent structures, large-scale flows and correlated dynamics in suspensions of motile particles such as swimming micro-organisms or artificial microswimmers is studied using direct particle simulations. A detailed model is proposed for a slender rod-like particle that propels itself in a viscous fluid by exerting a prescribed tangential stress on its surface, and a method is devised for the efficient calculation of hydrodynamic interactions in large-scale suspensions of such particles using slender-body theory and a smooth particle-mesh Ewald algorithm. Simulations are performed with periodic boundary conditions for various system sizes and suspension volume fractions, and demonstrate a transition to large-scale correlated motions in suspensions of rear-actuated swimmers, or Pushers, above a critical volume fraction or system size. This transition, which is not observed in suspensions of head-actuated swimmers, or Pullers, is seen most clearly in particle velocity and passive tracer statistics. These observations are consistent with predictions from our previous mean-field kinetic theory, one of which states that instabilities will arise in uniform isotropic suspensions of Pushers when the product of the linear system size with the suspension volume fraction exceeds a given threshold. We also find that the collective dynamics of Pushers result in giant number fluctuations, local alignment of swimmers and strongly mixing flows. Suspensions of Pullers, which evince no large-scale dynamics, nonetheless display interesting deviations from the random isotropic state.

Excitable patterns in active nematics

We analyze a model of mutually propelled filaments suspended in a two-dimensional solvent. The system undergoes a mean-field isotropic-nematic transition for large enough filament concentrations, and the nematic order parameter is allowed to vary in space and time. We show that the interplay between nonuniform nematic order, activity, and flow results in spatially modulated relaxation oscillations, similar to those seen in excitable media. In this regime the dynamics consists of nearly stationary periods separated by "bursts" of activity in which the system is elastically distorted and solvent is pumped throughout. At even higher activity, the dynamics becomes chaotic.

Kinetic theory of flocking: derivation of hydrodynamic equations

It is shown how to explicitly coarse-grain the microscopic dynamics of the rule-based Vicsek model for self-propelled agents. The hydrodynamic equations are derived by means of an Enskog-type kinetic theory. Expressions for all transport coefficients are given. The transition from a disordered to a flocking state, which at large particle speeds appears to be a fluctuation-induced first-order phase transition, is studied numerically and analytically.

Polar patterns of driven filaments

The emergence of collective motion exhibited by systems ranging from flocks of animals to self-propelled microorganisms to the cytoskeleton is a ubiquitous and fascinating self-organization phenomenon. Similarities between these systems, such as the inherent polarity of the constituents, a density-dependent transition to ordered phases or the existence of very large density fluctuations, suggest universal principles underlying pattern formation. This idea is followed by theoretical models at all levels of description: micro- or mesoscopic models directly map local forces and interactions using only a few, preferably simple, interaction rules, and more macroscopic approaches in the hydrodynamic limit rely on the systems' generic symmetries. All these models characteristically have a broad parameter space with a manifold of possible patterns, most of which have not yet been experimentally verified. The complexity of interactions and the limited parameter control of existing experimental systems are major obstacles to our understanding of the underlying ordering principles. Here we demonstrate the emergence of collective motion in a high-density motility assay that consists of highly concentrated actin filaments propelled by immobilized molecular motors in a planar geometry. Above a critical density, the filaments self-organize to form coherently moving structures with persistent density modulations, such as clusters, swirls and interconnected bands. These polar nematic structures are long lived and can span length scales orders of magnitudes larger than their constituents. Our experimental approach, which offers control of all relevant system parameters, complemented by agent-based simulations, allows backtracking of the assembly and disassembly pathways to the underlying local interactions. We identify weak and local alignment interactions to be essential for the observed formation of patterns and their dynamics. The presented minimal polar-pattern-forming system may thus provide new insight into emerging order in the broad class of active fluids and self-propelled particles.