Fluctuations and pattern formation in self-propelled particles

Banded phases in topological flocks

Flocking phase transitions arise in many aligning active soft matter systems, and an interesting question concerns the role of "topological" vs. "metric" interactions on these transitions. While recent theoretical work suggests that the order-disorder transition in these polar aligning models is universally first order, numerical studies have suggested that topological models may instead have a continuous transition. Some recent simulations have found that some variations of topologically interacting flocking agents have a discontinuous transition, but unambiguous observations of phase coexistence using common Voronoi-based alignment remains elusive. In this work, we use a custom GPU-accelerated simulation package to perform million-particle-scale simulations of a Voronoi-Vicsek model in which alignment interactions stem from an XY-like Hamiltonian. By accessing such large systems on appropriately long time scales and in the time-continuous limit, we are able to show a regime of stable phase coexistence between the ordered and disordered phases, confirming the discontinuous nature of this transition in the thermodynamic limit.

Visual collective behaviors on spherical robots

The implementation of collective motion, traditionally, disregard the limited sensing capabilities of an individual, to instead assuming an omniscient perception of the environment. This study implements a visual flocking model in a 'robot-in-the-loop' approach to reproduce these behaviors with a flock composed of 10 independent spherical robots. The model achieves robotic collective motion by only using panoramic visual information of each robot, such as retinal position, optical size and optic flow of the neighboring robots. We introduce a virtual anchor to confine the collective robotic movements so to avoid wall interactions. For the first time, a simple visual robot-in-the-loop approach succeed in reproducing several collective motion phases, in particular, swarming, and milling. Another milestone achieved with by this model is bridging the gap between simulation and physical experiments by demonstrating nearly identical behaviors in both environments with the same visual model. To conclude, we show that our minimal visual collective motion model is sufficient to recreate most collective behaviors on a robot-in-the-loop system that be implemented using several individuals, behaves as numerical simulations predict and is easily comparable to traditional models.

Spatiotemporal control of structure and dynamics in a polar active fluid

We apply optimal control theory to a model of a polar active fluid (the Toner-Tu model), with the objective of driving the system into particular emergent dynamical behaviors or programming switching between states on demand. We use the effective self-propulsion speed as the control parameter (i.e. the means of external actuation). We identify control protocols that achieve outcomes such as relocating asters to targeted positions, forcing propagating solitary waves to reorient to a particular direction, and switching between stationary asters and propagating fronts. We analyze the solutions to identify generic principles for controlling polar active fluids. Our findings have implications for achieving spatiotemporal control of active polar systems in experiments, particularly in vitro cytoskeletal systems. Additionally, this research paves the way for leveraging optimal control methods to engineer the structure and dynamics of active fluids more broadly.

Following Your Nose: Autochemotaxis and Other Mechanisms for Spinodal Decomposition in Flocks

We develop the hydrodynamic theory of dry, polar ordered, active matter ("flocking") with autochemotaxis; i.e., self-propelled entities moving in the same direction, each emitting a substance which attracts the others (e.g., ants). We find that sufficiently strong autochemotaxis leads to an instability to phase separation into one high and one low density band. This is very analogous to both equilibrium phase separation, and "motility induced phase separation" and can occur in flocks due to any microscopic mechanism (e.g., sufficiently strong attractive interactions) that makes the entities cohere.

Spinodal decomposition and phase separation in polar active matter

We develop and study the hydrodynamic theory of flocking with autochemotaxis. This describes large collections of self-propelled entities all spontaneously moving in the same direction, each emitting a substance which attracts the others (e.g., ants). The theory combines features of the Keller-Segel model for autochemotaxis with the Toner-Tu theory of flocking. We find that sufficiently strong autochemotaxis leads to an instability of the uniformly moving state (the "flock"), in which bands of different density form moving parallel to the mean flock velocity with different speeds. This instability is, therefore, completely different from the well-known "banding instability," in which bands form perpendicular to the mean flock velocity. The bands we find, which are reminiscent of ant trails, coarsen over time to reach a phase-separated state, in which one high-density and one low-density band fill the entire system. The same instability, described by the same hydrodynamic theory, can occur in flocks phase separating due to any microscopic mechanism (e.g., sufficiently strong attractive interactions). Although in many ways analogous to equilibrium phase separation via spinodal decomposition, the two steady-state densities here are determined not by a common tangent construction, as in equilibrium, but by an uncommon tangent construction very similar to that found for motility-induced phase separation of disordered active particles. Our analytic theory agrees well with our numerical simulations of our equations of motion.

Current reversal in polar flock at order-disorder interface

We studied a system of polar self-propelled particles (SPPs) on a thin rectangular channel designed into three regions of order-disorder-order. The division of the three regions is made on the basis of the noise SPPs experience in the respective regions. The noise in the two wide regions is chosen lower than the critical noise of order-disorder transition and noise in the middle region or interface is higher than the critical noise. This makes the geometry of the system analogous to the Josephson junction (JJ) in solid-state physics. Keeping all other parameters fixed, we study the properties of the moving SPPs in the bulk as well as along the interface for different widths of the junction. On increasing interface width, the system shows an order-to-disorder transition from coherent moving SPPs in the whole system to the interrupted current for large interface width. Surprisingly, inside the interface, we observed the current reversal for intermediate widths of the interface. Such current reversal is due to the strong randomness present inside the interface, which makes the wall of the interface reflecting. Hence, our study gives new interesting collective properties of SPPs at the interface which can be useful to design switching devices using active agents.

Symbiotic dynamics in living liquid crystals

An amalgam of nematic liquid crystals and active matter, referred to as living liquid crystals, is a promising self-healing material with futuristic applications for targeted delivery of information and microcargo. We provide a phenomenological model to study the symbiotic pattern dynamics in this contemporary system using the Toner-Tu model for active matter (AM), the Landau-de Gennes free energy for liquid crystals (LCs), and an experimentally motivated coupling term that favours coalignment of the active and nematic components. Our extensive theoretical studies unfold two novel steady states, chimeras and solitons, with sharp regions of distinct orientational order that sweep through the coupled system in synchrony. The induced dynamics in the passive nematic is unprecedented. We show that the symbiotic dynamics of the AM and LC components can be exploited to induce and manipulate order in an otherwise disordered system.

Phase ordering, topological defects, and turbulence in the three-dimensional incompressible Toner-Tu equation

We investigate the phase-ordering dynamics of the incompressible Toner-Tu equation in three dimensions. We show that the phase ordering proceeds via defect merger events and the dynamics is controlled by the Reynolds number Re. At low Re, the dynamics is similar to that of the Ginzburg-Landau equation. At high Re, turbulence controls phase ordering. In particular, we observe a forward energy cascade from the coarsening length scale to the dissipation scale, clustering of defects, and multiscaling in velocity correlations.

Swarming Bottom Feeders: Flocking at Solid-Liquid Interfaces

We present the hydrodynamic theory of coherent collective motion ("flocking") at a solid-liquid interface, and many of its predictions for experiment. We find that such systems are stable, and have long-range orientational order, over a wide range of parameters. When stable, these systems exhibit "giant number fluctuations," which grow as the 3/4th power of the mean number. Stable systems also exhibit anomalous rapid diffusion of tagged particles suspended in the passive fluid along any directions in a plane parallel to the solid-liquid interface, whereas the diffusivity along the direction perpendicular to the plane is not anomalous. In the remaining parameter space, the system becomes unstable.

Density fluctuations and energy spectra of 3D bacterial suspensions

Giant number fluctuations are often considered as a hallmark of the emergent nonequilibrium dynamics of active fluids. However, these anomalous density fluctuations have only been reported experimentally in two-dimensional dry active systems heretofore. Here, we investigate density fluctuations of bulk Escherichia coli suspensions, a paradigm of three-dimensional (3D) wet active fluids. Our experiments demonstrate the existence and quantify the scaling relation of giant number fluctuations in 3D bacterial suspensions. Surprisingly, the anomalous scaling persists at small scales in low-concentration suspensions well before the transition to active turbulence, reflecting the long-range nature of hydrodynamic interactions of 3D wet active fluids. To illustrate the origin of the density fluctuations, we measure the energy spectra of suspension flows and explore the density-energy coupling in both the steady and transient states of active turbulence. A scale-invariant density-independent correlation between density fluctuations and energy spectra is uncovered across a wide range of length scales. In addition, our experiments show that the energy spectra of bacterial turbulence exhibit the scaling of 3D active nematic fluids, challenging the common view of dense bacterial suspensions as active polar fluids.

Hydrodynamic theory of flocking at a solid-liquid interface: Long-range order and giant number fluctuations

We construct the hydrodynamic theory of coherent collective motion ("flocking") at a solid-liquid interface. The polar order parameter and concentration of a collection of "active" (self-propelled) particles at a planar interface between a passive, isotropic bulk fluid and a solid surface are dynamically coupled to the bulk fluid. We find that such systems are stable, and have long-range orientational order, over a wide range of parameters. When stable, these systems exhibit "giant number fluctuations," i.e., large fluctuations of the number of active particles in a fixed large area. Specifically, these number fluctuations grow as the 3/4th power of the mean number within the area. Stable systems also exhibit anomalously rapid diffusion of tagged particles suspended in the passive fluid along any directions in a plane parallel to the solid-liquid interface, whereas the diffusivity along the direction perpendicular to the plane is nonanomalous. In other parameter regimes, the system becomes unstable.

Non-reciprocal phase transitions

Out of equilibrium, a lack of reciprocity is the rule rather than the exception. Non-reciprocity occurs, for instance, in active matter^1-6, non-equilibrium systems^7-9, networks of neurons^10,11, social groups with conformist and contrarian members¹², directional interface growth phenomena^13-15 and metamaterials^16-20. Although wave propagation in non-reciprocal media has recently been closely studied^1,16-20, less is known about the consequences of non-reciprocity on the collective behaviour of many-body systems. Here we show that non-reciprocity leads to time-dependent phases in which spontaneously broken continuous symmetries are dynamically restored. We illustrate this mechanism with simple robotic demonstrations. The resulting phase transitions are controlled by spectral singularities called exceptional points²¹. We describe the emergence of these phases using insights from bifurcation theory^22,23 and non-Hermitian quantum mechanics^24,25. Our approach captures non-reciprocal generalizations of three archetypal classes of self-organization out of equilibrium: synchronization, flocking and pattern formation. Collective phenomena in these systems range from active time-(quasi)crystals to exceptional-point-enforced pattern formation and hysteresis. Our work lays the foundation for a general theory of critical phenomena in systems whose dynamics is not governed by an optimization principle.

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.

How do clusters in phase-separating active matter systems grow? A study for Vicsek activity in systems undergoing vapor-solid transition

Via molecular dynamics simulations, we have studied the kinetics of vapor-"solid" phase transition in an active matter model in which self-propulsion is introduced via the well-known Vicsek rule. The overall density of the particles is chosen in such a way that the evolution morphology consists of disconnected clusters that are defined as regions of high density of particles. Our focus has been on understanding the influence of the above-mentioned self-propulsion on structure and growth of these clusters by comparing the results with those for the passive limit of the model that also exhibits vapor-"solid" transition. While in the passive case growth occurs due to a standard diffusive mechanism, the Vicsek activity leads to very rapid growth, via a process that is practically equivalent to the ballistic aggregation mechanism. The emerging growth law in the latter case has been accurately estimated and explained by invoking information on velocity and structural aspects of the clusters into a relevant theory. Some of these results are also discussed with reference to a model for active Brownian particles.

Coarsening in the two-dimensional incompressible Toner-Tu equation: Signatures of turbulence

We investigate coarsening dynamics in the two-dimensional, incompressible Toner-Tu equation. We show that coarsening proceeds via vortex merger events, and the dynamics crucially depend on the Reynolds number Re. For low Re, the coarsening process has similarities to Ginzburg-Landau dynamics. On the other hand, for high Re, coarsening shows signatures of turbulence. In particular, we show the presence of an enstrophy cascade from the intervortex separation scale to the dissipation scale.

Flocking with a q-fold discrete symmetry: Band-to-lane transition in the active Potts model

We study the q-state active Potts model (APM) on a two-dimensional lattice in which self-propelled particles have q internal states corresponding to the q directions of motion. A local alignment rule inspired by the ferromagnetic q-state Potts model and self-propulsion via biased diffusion according to the internal particle states elicits collective motion at high densities and low noise. We formulate a coarse-grained hydrodynamic theory with which we compute the phase diagrams of the APM for q=4 and q=6 and analyze the flocking dynamics in the coexistence region, where the high-density (polar liquid) phase forms a fluctuating stripe of coherently moving particles on the background of the low-density (gas) phase. A reorientation transition of the phase-separated profiles from transversal band motion to longitudinal lane formation is found, which is absent in the Vicsek model and the active Ising model. The origin of this reorientation transition is revealed by a stability analysis: for large velocities the transverse diffusivity approaches zero and stabilizes lanes. Computer simulations corroborate the analytical predictions of the flocking and reorientation transitions and validate the phase diagrams of the APM.

Towards an analytical description of active microswimmers in clean and in surfactant-covered drops

Geometric confinements are frequently encountered in the biological world and strongly affect the stability, topology, and transport properties of active suspensions in viscous flow. Based on a far-field analytical model, the low-Reynolds-number locomotion of a self-propelled microswimmer moving inside a clean viscous drop or a drop covered with a homogeneously distributed surfactant, is theoretically examined. The interfacial viscous stresses induced by the surfactant are described by the well-established Boussinesq-Scriven constitutive rheological model. Moreover, the active agent is represented by a force dipole and the resulting fluid-mediated hydrodynamic couplings between the swimmer and the confining drop are investigated. We find that the presence of the surfactant significantly alters the dynamics of the encapsulated swimmer by enhancing its reorientation. Exact solutions for the velocity images for the Stokeslet and dipolar flow singularities inside the drop are introduced and expressed in terms of infinite series of harmonic components. Our results offer useful insights into guiding principles for the control of confined active matter systems and support the objective of utilizing synthetic microswimmers to drive drops for targeted drug delivery applications.

Traveling bands, clouds, and vortices of chiral active matter

We consider stochastic dynamics of self-propelled particles with nonlocal normalized alignment interactions subject to phase lag. The role of the lag is to indirectly generate chirality into particle motion. To understand large-scale behavior, we derive a continuum description of an active Brownian particle flow with macroscopic scaling in the form of a partial differential equation for a one-particle probability density function. Due to indirect chirality, we find a spatially homogeneous nonstationary analytic solution for this class of equations. Our development of kinetic and hydrodynamic theories towards such a solution reveals the existence of a wide variety of spatially nonhomogeneous patterns reminiscent of traveling bands, clouds, and vortical structures of linear active matter. Our model may thereby serve as the basis for understanding the nature of chiral active media and designing multiagent swarms with designated behavior.

Flocking transition within the framework of Kuramoto paradigm for synchronization: Clustering and the role of the range of interaction

A Kuramoto-type approach to address flocking phenomena is presented. First, we analyze a simple generalization of the Kuramoto model for interacting active particles, which is able to show the flocking transition (the emergence of coordinated movements in a group of interacting self-propelled agents). In the case of all-to-all interaction, the proposed model reduces to the Kuramoto model for phase synchronization of identical motionless noisy oscillators. In general, the nature of this non-equilibrium phase transition depends on the range of interaction between the particles. Namely, for a small range of interaction, the transition is first order, while for a larger range of interaction, it is a second order transition. Moreover, for larger interaction ranges, the system exhibits the same features as in the case of all-to-all interaction, showing a spatially homogeneous flux when flocking phenomenon has emerged, while for lower interaction ranges, the flocking transition is characterized by cluster formation. We compute the phase diagram of the model, where we distinguish three phases as a function of the range of interaction and the effective coupling strength: a disordered phase, a spatially homogeneous flocking phase, and a cluster-flocking phase. Then, we present a general discussion about the applicability of this way of modeling to more realistic and general situations, ending with a brief presentation of a second example (a second model with a conservative interaction) where the flocking transition may be studied within the framework that we are proposing.

Modelling collective cell motion: are on- and off-lattice models equivalent?

Biological processes, such as embryonic development, wound repair and cancer invasion, or bacterial swarming and fruiting body formation, involve collective motion of cells as a coordinated group. Collective cell motion of eukaryotic cells often includes interactions that result in polar alignment of cell velocities, while bacterial patterns typically show features of apolar velocity alignment. For analysing the population-scale effects of these different alignment mechanisms, various on- and off-lattice agent-based models have been introduced. However, discriminating model-specific artefacts from general features of collective cell motion is challenging. In this work, we focus on equivalence criteria at the population level to compare on- and off-lattice models. In particular, we define prototypic off- and on-lattice models of polar and apolar alignment, and show how to obtain an on-lattice from an off-lattice model of velocity alignment. By characterizing the behaviour and dynamical description of collective migration models at the macroscopic level, we suggest the type of phase transitions and possible patterns in the approximative macroscopic partial differential equation descriptions as informative equivalence criteria between on- and off-lattice models. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Decentralized Control for Swarm Robots That Can Effectively Execute Spatially Distributed Tasks

A swarm robotic system is a system in which multiple robots cooperate to fulfill a macroscopic function. Many swarm robots have been developed for various purposes. This study aims to design swarm robots capable of executing spatially distributed tasks effectively, which can be potentially used for tasks such as search-and-rescue operation and gathering scattered garbage in rooms. We propose a simple decentralized control scheme for swarm robots by extending our previously proposed non-reciprocal-interaction-based model. Each robot has an internal state, called its workload. Each robot first moves randomly to find a task, and when it does, its workload increases, and then it attracts its neighboring robots to ask for their help. We demonstrate, via simulations, that the proposed control scheme enables the robots to effectively execute multiple tasks in parallel under various environments. Fault tolerance of the proposed system is also demonstrated.

Phase separation and state oscillation of active inertial particles

Active matter systems are of great interest for their novel out-of-equilibrium collective behavior. Active Brownian particles (ABPs) are known to exhibit clustering and motility-induced phase separation, and there have been many studies revealing this rich behavior in the overdamped limit of Brownian motion, where inertial effects are insignificant. Here we simulate an Active Inertial Particle (AIP) model where we focus on the underdamped, rather than overdamped limit, to study the interplay between particle inertia and collective behavior, such as phase separation. We show that inertia reduces particle motility due to collisions and a longer time delay for particles to regain speed, thereby suppressing phase separation relative to that observed in the overdamped limit. Additionally, we observe interesting oscillatory behavior between a phase separated steady-state and a homogeneous fluid state that results from inertia-induced collective motion within active clusters due to momentum transfer. Such oscillatory behavior has been reported for ABP systems with particle shape anisotropy, where collective motion is mediated by particle shape anisotropy. Furthermore, we confirm that there is no single characteristic frequency for the oscillatory behavior. The power spectral density is a power law in the high frequency domain, with an exponent close to -2.5.

Coarsening dynamics in the Vicsek model of active matter

We study the flocking model introduced by Vicsek et al. (Phys. Rev. Lett. 75, 1226 (1995)) in the "coarsening" regime. At standard self-propulsion speeds, we find two distinct growth laws for the coupled density and velocity fields. The characteristic length scale of the density domains grows as [Formula: see text] (with [Formula: see text] , while the velocity length scale grows much faster, viz., [Formula: see text] (with [Formula: see text] . The spatial fluctuations in the density and velocity fields are studied by calculating the two-point correlation function and the structure factor, which show deviations from the well-known Porod's law. This is a natural consequence of scattering from irregular morphologies that dynamically arise in the system. At large values of the scaled wave vector, the scaled structure factors for the density and velocity fields decay with powers -2.6 and -1.52 , respectively.

Flip motion of solitary wave in an Ising-type Vicsek model

An Ising-type Vicsek model is proposed for collective motion and sudden direction change in a population of self-propelled particles. Particles move on a linear lattice with velocity +1 or -1 in the one-dimensional model. The probability of the velocity of a particle at the next step is determined by the number difference of the right- and left-moving particles at the present lattice site and its nearest-neighboring sites. A solitary wave appears also in our model similarly to previous models. In some parameter range, the moving direction of the solitary wave sometimes changes rather suddenly, which is like the sudden change of direction of a flock of birds. We study the average reversal time of traveling direction numerically and compare the results with a mean-field theory. The one-dimensional model is generalized to a two-dimensional model. Flip motion of a bandlike soliton is observed in the two-dimensional model.

Frequency-dependent higher-order Stokes singularities near a planar elastic boundary: Implications for the hydrodynamics of an active microswimmer near an elastic interface

The emerging field of self-driven active particles in fluid environments has recently created significant interest in the biophysics and bioengineering communities owing to their promising future for biomedical and technological applications. These microswimmers move autonomously through aqueous media, where under realistic situations they encounter a plethora of external stimuli and confining surfaces with peculiar elastic properties. Based on a far-field hydrodynamic model, we present an analytical theory to describe the physical interaction and hydrodynamic couplings between a self-propelled active microswimmer and an elastic interface that features resistance toward shear and bending. We model the active agent as a superposition of higher-order Stokes singularities and elucidate the associated translational and rotational velocities induced by the nearby elastic boundary. Our results show that the velocities can be decomposed in shear and bending related contributions which approach the velocities of active agents close to a no-slip rigid wall in the steady limit. The transient dynamics predict that contributions to the velocities of the microswimmer due to bending resistance are generally more pronounced than those due to shear resistance. Bending can enhance (suppress) the velocities resulting from higher-order singularities whereas the shear related contribution decreases (increases) the velocities. Most prominently, we find that near an elastic interface of only energetic resistance toward shear deformation, such as that of an elastic capsule designed for drug delivery, a swimming bacterium undergoes rotation of the same sense as observed near a no-slip wall. In contrast to that, near an interface of only energetic resistance toward bending, such as that of a fluid vesicle or liposome, we find a reversed sense of rotation. Our results provide insight into the control and guidance of artificial and synthetic self-propelling active microswimmers near elastic confinements.

A statistical physics view of swarming bacteria

Bacterial swarming is a collective mode of motion in which cells migrate rapidly over surfaces, forming dynamic patterns of whirls and jets. This review presents a physical point of view of swarming bacteria, with an emphasis on the statistical properties of the swarm dynamics as observed in experiments. The basic physical principles underlying the swarm and their relation to contemporary theories of collective motion and active matter are reviewed and discussed in the context of the biological properties of swarming cells. We suggest a paradigm according to which bacteria have optimized some of their physical properties as a strategy for rapid surface translocation. In other words, cells take advantage of favorable physics, enabling efficient expansion that enhances survival under harsh conditions.

Sounds and hydrodynamics of polar active fluids

Spontaneously flowing liquids have been successfully engineered from a variety of biological and synthetic self-propelled units^1-11. Together with their orientational order, wave propagation in such active fluids has remained a subject of intense theoretical studies^12-17. However, the experimental observation of this phenomenon has remained elusive. Here, we establish and exploit the propagation of sound waves in colloidal active materials with broken rotational symmetry. We demonstrate that two mixed modes, coupling density and velocity fluctuations, propagate along all directions in colloidal-roller fluids. We then show how the six material constants defining the linear hydrodynamics of these active liquids can be measured from their spontaneous fluctuation spectrum, while being out of reach of conventional rheological methods. This active-sound spectroscopy is not specific to synthetic active materials and could provide a quantitative hydrodynamic description of herds, flocks and swarms from inspection of their large-scale fluctuations^18-21.

Collective behavior of self-propelled rods with quorum sensing

Active agents-like phoretic particles, bacteria, sperm, and cytoskeletal filaments in motility assays-show a large variety of motility-induced collective behaviors, such as aggregation, clustering, and phase separation. The behavior of dense suspensions of engineered phoretic particles and of bacteria during biofilm formation is determined by two qualitatively different physical mechanisms: (i) volume exclusion (short-range steric repulsion) and (ii) quorum sensing (longer-range reduced propulsion due to alteration of the local chemical environment). To systematically characterize such systems, we study semi-penetrable self-propelled rods in two dimensions, with a propulsion force that decreases with increasing local rod density, by employing Brownian dynamics simulations. Volume exclusion and quorum sensing both lead to phase separation; however, the structure of the systems and the rod dynamics vastly differ. Quorum sensing enhances the polarity of the clusters, induces perpendicularity of rods at the cluster borders, and enhances cluster formation. For systems where the rods essentially become passive at high densities, formation of asters and stripes is observed. Systems of rods with larger aspect ratios show more ordered structures compared to those with smaller aspect ratios, due to their stronger alignment, with almost circular asters for strongly density-dependent propulsion force. With increasing range of the quorum-sensing interaction, the local density decreases, asters become less stable, and polar hedgehog clusters and clusters with domains appear.

Hydrodynamics of shape-driven rigidity transitions in motile tissues

In biological tissues, it is now well-understood that mechanical cues are a powerful mechanism for pattern regulation. While much work has focused on interactions between cells and external substrates, recent experiments suggest that cell polarization and motility might be governed by the internal shear stiffness of nearby tissue, deemed "plithotaxis". Meanwhile, other work has demonstrated that there is a direct relationship between cell shapes and tissue shear modulus in confluent tissues. Joining these two ideas, we develop a hydrodynamic model that couples cell shape, and therefore tissue stiffness, to cell motility and polarization. Using linear stability analysis and numerical simulations, we find that tissue behavior can be tuned between largely homogeneous states and patterned states such as asters, controlled by a composite "morphotaxis" parameter that encapsulates the nature of the coupling between shape and polarization. The control parameter is in principle experimentally accessible, and depends both on whether a cell tends to move in the direction of lower or higher shear modulus, and whether sinks or sources of polarization tend to fluidize the system.

A continuous-time persistent random walk model for flocking

A classical random walker is characterized by a random position and velocity. This sort of random walk was originally proposed by Einstein to model Brownian motion and to demonstrate the existence of atoms and molecules. Such a walker represents an inanimate particle driven by environmental fluctuations. On the other hand, there are many examples of so-called "persistent random walkers," including self-propelled particles that are able to move with almost constant speed while randomly changing their direction of motion. Examples include living entities (ranging from flagellated unicellular organisms to complex animals such as birds and fish), as well as synthetic materials. Here we discuss such persistent non-interacting random walkers as a model for active particles. We also present a model that includes interactions among particles, leading to a transition to flocking, that is, to a net flux where the majority of the particles move in the same direction. Moreover, the model exhibits secondary transitions that lead to clustering and more complex spatially structured states of flocking. We analyze all these transitions in terms of bifurcations using a number of mean field strategies (all to all interaction and advection-reaction equations for the spatially structured states), and compare these results with direct numerical simulations of ensembles of these interacting active particles.

Hydrodynamic theory of active matter

We review the general hydrodynamic theory of active soft materials that is motivated in particular by biological matter. We present basic concepts of irreversible thermodynamics of spatially extended multicomponent active systems. Starting from the rate of entropy production, we identify conjugate thermodynamic fluxes and forces and present generic constitutive equations of polar active fluids and active gels. We also discuss angular momentum conservation which plays a role in the the physics of active chiral gels. The irreversible thermodynamics of active gels provides a general framework to discuss the physics that underlies a wide variety of biological processes in cells and in multicellular tissues.

Emergence of coexisting ordered states in active matter systems

Active systems can produce a far greater variety of ordered patterns than conventional equilibrium systems. In particular, transitions between disorder and either polar- or nematically ordered phases have been predicted and observed in two-dimensional active systems. However, coexistence between phases of different types of order has not been reported. We demonstrate the emergence of dynamic coexistence of ordered states with fluctuating nematic and polar symmetry in an actomyosin motility assay. Combining experiments with agent-based simulations, we identify sufficiently weak interactions that lack a clear alignment symmetry as a prerequisite for coexistence. Thus, the symmetry of macroscopic order becomes an emergent and dynamic property of the active system. These results provide a pathway by which living systems can express different types of order by using identical building blocks.

State diagram of a three-sphere microswimmer in a channel

Geometric confinements are frequently encountered in soft matter systems and in particular significantly alter the dynamics of swimming microorganisms in viscous media. Surface-related effects on the motility of microswimmers can lead to important consequences in a large number of biological systems, such as biofilm formation, bacterial adhesion and microbial activity. On the basis of low-Reynolds-number hydrodynamics, we explore the state diagram of a three-sphere microswimmer under channel confinement in a slit geometry and fully characterize the swimming behavior and trajectories for neutral swimmers, puller- and pusher-type swimmers. While pushers always end up trapped at the channel walls, neutral swimmers and pullers may further perform a gliding motion and maintain a stable navigation along the channel. We find that the resulting dynamical system exhibits a supercritical pitchfork bifurcation in which swimming in the mid-plane becomes unstable beyond a transition channel height while two new stable limit cycles or fixed points that are symmetrically disposed with respect to the channel mid-height emerge. Additionally, we show that an accurate description of the averaged swimming velocity and rotation rate in a channel can be captured analytically using the method of hydrodynamic images, provided that the swimmer size is much smaller than the channel height.

Pattern, growth, and aging in aggregation kinetics of a Vicsek-like active matter model

Via molecular dynamics simulations, we study kinetics in a Vicsek-like phase-separating active matter model. Quantitative results, for isotropic bicontinuous pattern, are presented on the structure, growth, and aging. These are obtained via the two-point equal-time density-density correlation function, the average domain length, and the two-time density autocorrelation function. Both the correlation functions exhibit basic scaling properties, implying self-similarity in the pattern dynamics, for which the average domain size exhibits a power-law growth in time. The equal-time correlation has a short distance behavior that provides reasonable agreement between the corresponding structure factor tail and the Porod law. The autocorrelation decay is a power-law in the average domain size. Apart from these basic similarities, the overall quantitative behavior of the above-mentioned observables is found to be vastly different from those of the corresponding passive limit of the model which also undergoes phase separation. The functional forms of these have been quantified. An exceptionally rapid growth in the active system occurs due to fast coherent motion of the particles, mean-squared-displacements of which exhibit multiple scaling regimes, including a long time ballistic one.

Collective motion of active Brownian particles with polar alignment

We present a comprehensive computational study of the collective behavior emerging from the competition between self-propulsion, excluded volume interactions and velocity-alignment in a two-dimensional model of active particles. We consider an extension of the active brownian particles model where the self-propulsion direction of the particles aligns with the one of their neighbors. We analyze the onset of collective motion (flocking) in a low-density regime (10% surface area) and show that it is mainly controlled by the strength of velocity-alignment interactions: the competition between self-propulsion and crowding effects plays a minor role in the emergence of flocking. However, above the flocking threshold, the system presents a richer pattern formation scenario than analogous models without alignment interactions (active brownian particles) or excluded volume effects (Vicsek-like models). Depending on the parameter regime, the structure of the system is characterized by either a broad distribution of finite-sized polar clusters or the presence of an amorphous, highly fluctuating, large-scale traveling structure which can take a lane-like or band-like form (and usually a hybrid structure which is halfway in between both). We establish a phase diagram that summarizes collective behavior of polar active brownian particles and propose a generic mechanism to describe the complexity of the large-scale structures observed in systems of repulsive self-propelled particles.

Collective Behavior of Chiral Active Matter: Pattern Formation and Enhanced Flocking

We generalize the Vicsek model to describe the collective behavior of polar circle swimmers with local alignment interactions. While the phase transition leading to collective motion in 2D (flocking) occurs at the same interaction to noise ratio as for linear swimmers, as we show, circular motion enhances the polarization in the ordered phase (enhanced flocking) and induces secondary instabilities leading to structure formation. Slow rotations promote macroscopic droplets with late time sizes proportional to the system size (indicating phase separation) whereas fast rotations generate patterns consisting of phase synchronized microflocks with a controllable characteristic size proportional to the average single-particle swimming radius. Our results defy the viewpoint that monofrequent rotations form a vapid extension of the Vicsek model and establish a generic route to pattern formation in chiral active matter with possible applications for understanding and designing rotating microflocks.

Curvature-controlled defect dynamics in active systems

We have studied the collective motion of polar active particles confined to ellipsoidal surfaces. The geometric constraints lead to the formation of vortices that encircle surface points of constant curvature (umbilics). We have found that collective motion patterns are particularly rich on ellipsoids with four umbilics where vortices tend to be located near pairs of umbilical points to minimize their interaction energy. Our results provide a perspective on the migration of living cells, which most likely use the information provided from the curved substrate geometry to guide their collective motion.

Pattern Formation in Polymerizing Actin Flocks: Spirals, Spots, and Waves without Nonlinear Chemistry

We propose a model solely based on actin treadmilling and polymerization which describes many characteristic states of actin-wave formation: spots, spirals, and traveling waves. In our model, as in experiments on cells recovering motility following actin depolymerization, we choose an isotropic low-density initial condition; polymerization of actin filaments then raises the density towards the Onsager threshold where they align. We show that this alignment, in turn, destabilizes the isotropic phase and generically induces transient actin spots or spirals as part of the dynamical pathway towards a polarized phase which can either be uniform or consist of a series of actin-wave trains (flocks). Our results uncover a universal route to actin-wave formation in the absence of any system-specific nonlinear biochemistry, and it may help to understand the mechanism underlying the observation of actin spots and waves in vivo. They also suggest a minimal setup to design similar patterns in vitro.

Front structure and dynamics in dense colonies of motile bacteria: Role of active turbulence

We study the spreading of a bacterial colony undergoing turbulentlike collective motion. We present two minimalistic models to investigate the interplay between population growth and coherent structures arising from turbulence. Using direct numerical simulation of the proposed models we find that turbulence has two prominent effects on the spatial growth of the colony: (a) the front speed is enhanced, and (b) the front gets crumpled. Both these effects, which we highlight by using statistical tools, are markedly different in our two models. We also show that the crumpled front structure and the passive scalar fronts in random flows are related in certain regimes.

Mode instabilities and dynamic patterns in a colony of self-propelled surfactant particles covering a thin liquid layer

We consider a colony of point-like self-propelled surfactant particles (swimmers) without direct interactions that cover a thin liquid layer on a solid support. The particles predominantly swim normal to the free film surface with only a small component parallel to the film surface. The coupled dynamics of the swimmer density and film height profile is captured in a long-wave model allowing for diffusive and convective transport of the swimmers (including rotational diffusion). The dynamics of the film height profile is determined by i) the upward pushing force of the swimmers onto the liquid-gas interface, ii) the solutal Marangoni force due to gradients in the swimmer concentration, and iii) the rotational diffusion of the swimmers together with the in-plane active motion. After reviewing and extending the analysis of the linear stability of the uniform state, we analyse the fully nonlinear dynamic equations and show that point-like swimmers, which only interact via long-wave deformations of the liquid film, self-organise in highly regular (standing, travelling, and modulated waves) and various irregular patterns.

Additivity, density fluctuations, and nonequilibrium thermodynamics for active Brownian particles

Using an additivity property, we study particle-number fluctuations in a system of interacting self-propelled particles, called active Brownian particles (ABPs), which consists of repulsive disks with random self-propulsion velocities. From a fluctuation-response relation, a direct consequence of additivity, we formulate a thermodynamic theory which captures the previously observed features of nonequilibrium phase transition in the ABPs from a homogeneous fluid phase to an inhomogeneous phase of coexisting gas and liquid. We substantiate the predictions of additivity by analytically calculating the subsystem particle-number distributions in the homogeneous fluid phase away from criticality where analytically obtained distributions are compatible with simulations in the ABPs.

Active Curved Polymers Form Vortex Patterns on Membranes

Recent in vitro experiments with FtsZ polymers show self-organization into different dynamic patterns, including structures reminiscent of the bacterial Z ring. We model FtsZ polymers as active particles moving along chiral, circular paths by Brownian dynamics simulations and a Boltzmann approach. Our two conceptually different methods point to a generic phase behavior. At intermediate particle densities, we find self-organization into vortex structures including closed rings. Moreover, we show that the dynamics at the onset of pattern formation is described by a generalized complex Ginzburg-Landau equation.

Nonequilibrium phase transitions, fluctuations and correlations in an active contractile polar fluid

We study the patterning, fluctuations and correlations of an active polar fluid consisting of contractile polar filaments on a two-dimensional substrate, using a hydrodynamic description. The steady states generically consist of arrays of inward pointing asters and show a continuous transition from a moving lamellar phase, a moving aster street, to a stationary aster lattice with no net polar order. We next study the effect of spatio-temporal athermal noise, parametrized by an active temperature TA, on the stability of the ordered phases. In contrast to its equilibrium counterpart, we find that the active crystal shows true long range order at low TA. On increasing TA, the asters dynamically remodel, concomitantly we find novel phase transitions characterized by bond-orientational and polar order upon "heating".

Role of the active viscosity and self-propelling speed in channel flows of active polar liquid crystals

We study channel flows of active polar liquid crystals (APLCs) focusing on the role played by the active viscosity (β) and the self-propelling speed (ω) on the formation and long time evolution of spontaneous flows using a continuum model. First, we study the onset of spontaneous flows by carrying out a linear stability analysis on two special steady states subject to various physical boundary conditions. We identify a single parameter b1, proportional to a linear combination of the active viscosity and the self-propelling speed, and inversely proportional to a Frank elastic constant, the solvent viscosity, and the liquid crystal relaxation time. We show that the active viscosity and the self-propelling speed influence the onset of spontaneous flows through b1 in that for any fixed value of the bulk activity parameter ζ, large enough |b1| can suppress the spontaneous flow. We then follow spontaneous flows in long time to further investigate the role of β and ω on spatial-temporal structures in the nonlinear regime numerically. The numerical study demonstrates a strong correlation between the most unstable eigenfunction obtained from the linear analysis and the terminal steady state or the persistent, traveling wave structure, revealing the genesis of flow and orientational structures in the active matter system. In the nonlinear regime, a nonzero b1 facilitates the formation of traveling waves in the case of boundary anchoring (the Dirichlet boundary condition) so long as the linear stability analysis predicts an onset of spontaneous flows; in the case of the free boundary condition (the Neumann boundary condition), a stable, spatially homogeneous tilted state always emerges in the presence of two active effects. Finally, we note that various fully out-of-plane spatio-temporal structures can emerge in long time dynamics depending on the boundary condition as well as the initial state of the polarity vector field.