Understanding collective dynamics of soft active colloids by binary scattering

Breakdown of Boltzmann-type models for the alignment of self-propelled rods

Studies in the collective motility of organisms use a range of analytical approaches to formulate continuous kinetic models of collective dynamics from rules or equations describing agent interactions. However, the derivation of these kinetic models often relies on Boltzmann's "molecular chaos" hypothesis, which assumes that correlations between individuals are short-lived. While this assumption is often the simplest way to derive tractable models, it is often not valid in practice due to the high levels of cooperation and self-organization present in biological systems. In this work, we illustrated this point by considering a general Boltzmann-type kinetic model for the alignment of self-propelled rods where rod reorientation occurs upon binary collisions. We examine the accuracy of the kinetic model by comparing numerical solutions of the continuous equations to an agent-based model that implements the underlying rules governing microscopic alignment. Even for the simplest case considered, our comparison demonstrates that the kinetic model fails to replicate the discrete dynamics due to the formation of rod clusters that violate statistical independence. Additionally, we show that introducing noise to limit cluster formation helps improve the agreement between the analytical model and agent simulations but does not restore the agreement completely. These results highlight the need to both develop and disseminate improved moment-closure methods for modeling biological and active matter systems.

Marginal speed confinement resolves the conflict between correlation and control in collective behaviour

Speed fluctuations of individual birds in natural flocks are moderate, due to the aerodynamic and biomechanical constraints of flight. Yet the spatial correlations of such fluctuations are scale-free, namely they have a range as wide as the entire group, a property linked to the capacity of the system to collectively respond to external perturbations. Scale-free correlations and moderate fluctuations set conflicting constraints on the mechanism controlling the speed of each agent, as the factors boosting correlation amplify fluctuations, and vice versa. Here, using a statistical field theory approach, we suggest that a marginal speed confinement that ignores small deviations from the natural reference value while ferociously suppressing larger speed fluctuations, is able to reconcile scale-free correlations with biologically acceptable group’s speed. We validate our theoretical predictions by comparing them with field experimental data on starling flocks with group sizes spanning an unprecedented interval of over two orders of magnitude.

Bird flocks are known to adjust the orientation and speed of individual birds giving rise to correlations that extend across very large groups. The authors show that marginal control provides an explanation of scale-free correlations of speed fluctuations in natural bird flocks of any sizes.

Dynamic self-organization of migrating cells under constraints by spatial confinement and epithelial integrity

Understanding how migrating cells can establish both dynamic structures and coherent dynamics may provide mechanistic insights to study how living systems acquire complex structures and functions. Recent studies revealed that intercellular contact communication plays a crucial role for establishing cellular dynamic self-organization (DSO) and provided a theoretical model of DSO for migrating solitary cells in a free space. However, to apply those understanding to situations in living organisms, we need to know the role of cell-cell communication for tissue dynamics under spatial confinements and epithelial integrity. Here, we expand the previous numerical studies on DSO to migrating cells subjected spatial confinement and/or epithelial integrity. An epithelial monolayer is simulated by combining the model of cellular DSO and the cellular vertex model in two dimensions for apical integrity. Under confinement to a small space, theoretical models of both solitary and epithelial cells exhibit characteristic coherent dynamics, including apparent swirling. We also find that such coherent dynamics can allow the cells to overcome the strong constraint due to spatial confinement and epithelial integrity. Furthermore, we demonstrate how epithelial cell clusters behave without spatial confinement and find various cluster dynamics, including spinning, migration and elongation.

Contact inhibition of locomotion generates collective cell migration without chemoattractants in an open domain

Neural crest cells are embryonic stem cells that migrate throughout embryos and, at different target locations, give rise to the formation of a variety of tissues and organs. The directional migration of the neural crest cells is experimentally described using a process referred to as contact inhibition of locomotion, by which cells redirect their movement upon the cell-cell contacts. However, it is unclear how the migration alignment is affected by the motility properties of the cells. Here, we theoretically model the migration alignment as a function of the motility dynamics and interaction of the cells in an open domain with a channel geometry. The results indicate that by increasing the influx rate of the cells into the domain a transition takes place from random movement to an organized collective migration, where the migration alignment is maximized and the migration time is minimized. This phase transition demonstrates that the cells can migrate efficiently over long distances without any external chemoattractant information about the direction of migration just based on local interactions with each other. The analysis of the dependence of this transition on the characteristic properties of cellular motility shows that the cell density determines the coordination of collective migration whether the migration domain is open or closed. In the open domain, this density is determined by a feedback mechanism between the flux and order parameter, which characterises the alignment of collective migration. The model also demonstrates that the coattraction mechanism proposed earlier is not necessary for collective migration and a constant flux of cells moving into the channel is sufficient to produce directed movement over arbitrary long distances.

Persistence length regulates emergent dynamics in active roller ensembles

Active colloidal fluids, biological and synthetic, often demonstrate complex self-organization and the emergence of collective behavior. Spontaneous formation of multiple vortices has been recently observed in a variety of active matter systems, however, the generation and tunability of the active vortices not controlled by geometrical confinement remain challenging. Here, we exploit the persistence length of individual particles in ensembles of active rollers to tune the formation of vortices and to orchestrate their characteristic sizes. We use two systems and employ two different approaches exploiting shape anisotropy or polarization memory of individual units for control of the persistence length. We characterize the dynamics of emergent multi-vortex states and reveal a direct link between the behavior of the persistence length and properties of the emergent vortices. We further demonstrate common features between the two systems including anti-ferromagnetic ordering of the neighboring vortices and active turbulent behavior with a characteristic energy cascade in the particles velocity field energy spectra. Our findings provide insights into the onset of spatiotemporal coherence in active roller systems and suggest a control knob for manipulation of dynamic self-assembly in active colloidal ensembles.

Dynamic Self-Organization of Idealized Migrating Cells by Contact Communication

This Letter investigates what forms of cellular dynamic self-organization are caused through intercellular contact communication based on a theoretical model in which migrating cells perform contact following and contact inhibition and attraction of locomotion. Tuning those strengths causes varieties of dynamic patterns. This further includes a novel form of collective migration, snakelike dynamic assembly. Scrutinizing this pattern reveals that cells in this state can accurately respond to an external directional cue but have no spontaneous global polar order.

Reconfigurable emergent patterns in active chiral fluids

Active fluids comprised of autonomous spinning units injecting energy and angular momentum at the microscopic level represent a promising platform for active materials design. The complexity of the accessible dynamic states is expected to dramatically increase in the case of chiral active units. Here, we use shape anisotropy of colloidal particles to introduce chiral rollers with activity-controlled curvatures of their trajectories and spontaneous handedness of their motion. By controlling activity through variations of the energizing electric field, we reveal emergent dynamic phases, ranging from a gas of spinners to aster-like vortices and rotating flocks, with either polar or nematic alignment of the particles. We demonstrate control and reversibility of these dynamic states by activity. Our findings provide insights into the onset of spatial and temporal coherence in a broad class of active chiral systems, both living and synthetic, and hint at design pathways for active materials based on self-organization and reconfigurability.

Adhesion-stabilizing long-distance transport of cells on tissue surface

The stable transport of migrating eukaryotic cells is essential in organ development and repair processes. However, the mechanism that preserves transport stability over long distances in organs is not fully understood. As the driving mechanism of cell migration, the expressions of heterophilic cell-cell adhesion between moving cells and scaffolding tissue have been observed in such transport. In this paper, we theoretically investigate this heterophilic adhesion, which is persistently polarized in the migrating cell, as a possible transport stabilization mechanism. The adhesion was examined on the basis of the cellular Potts model, and our results confirm the stabilization of the transport to be an effect of the persistence.

Emergence of self-organized multivortex states in flocks of active rollers

Active matter, both synthetic and biological, demonstrates complex spatiotemporal self-organization and the emergence of collective behavior. A coherent rotational motion, the vortex phase, is of great interest because of its ability to orchestrate well-organized motion of self-propelled particles over large distances. However, its generation without geometrical confinement has been a challenge. Here, we show by experiments and computational modeling that concentrated magnetic rollers self-organize into multivortex states in an unconfined environment. We find that the neighboring vortices more likely occur with the opposite sense of rotation. Our studies provide insights into the mechanism for the emergence of coherent collective motion on the macroscale from the coupling between microscale rotation and translation of individual active elements. These results may stimulate design strategies for self-assembled dynamic materials and microrobotics.

Phase transition and emergence of active temperature in an active Brownian system in underdamped background

We explore the role of inertia in the properties of active Brownian particles (ABPs) immersed in an underdamped background in two dimensions using Langevin dynamics computer simulation. Similar to an equilibrium two-dimensional passive interacting particle system, the system of ABPs transits from a liquid phase to a solid phase with the change in the coupling parameter, which is the ratio of interaction potential energy and thermal energy of the background solvent. Important qualitative and quantitative differences are found in the liquid-solid phase transition with increasing strength of activity as compared to those found in the conventional overdamped background limit. In the underdamped background, inherent activity is found to lead to a temperature, called the active temperature and defined by average velocity fluctuations of the ABPs, that is different from the fixed background solvent temperature. A new scaling law for active temperature as a function of activity strength is found near the liquid-solid boundary. Active temperature, which behaves similar to the thermodynamic equilibrium temperature, is also found to depend upon the interaction strength between the active particles and the strength of the background dissipation. With an increase in background dissipation, the difference between active temperature and the background solvent temperature decreases and the difference is found to eventually vanish in the overdamped limit, demonstrating the correctness of the calculation.

Three-body interactions drive the transition to polar order in a simple flocking model

A large class of mesoscopic or macroscopic flocking theories are coarse-grained from microscopic models that feature binary interactions as the chief aligning mechanism. However, while such theories seemingly predict the existence of polar order with just binary interactions, actomyosin motility assay experiments show that binary interactions are insufficient to obtain polar order, especially at high densities. To resolve this paradox, here we introduce a solvable one-dimensional flocking model and derive its stochastic hydrodynamics. We show that two-body interactions are insufficient to generate polar order unless the noise is non-Gaussian. We show that noisy three-body interactions in the microscopic theory allow us to capture all essential dynamical features of the flocking transition, in systems that achieve orientational order above a critical density.

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

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

Two types of exclusion interactions for self-propelled objects and collective motion induced by their combination

Exclusive interactions between self-driven objects may play crucial roles in their collective behavior, e.g., in collective migration of living cells. Here, such collective behavior is studied based on a simple but sufficient model taking account the exclusion effects, which incorporate the following two distinct kinds of exclusion interactions in two dimensions: The first is the mechanical exclusion wherein two objects mechanically repel each other when they overlap. The second is the scattering exclusion, wherein the directions along which each object tries to move are modulated to avoid overlapping. We propose a theoretical model based on two principles: (1) Each object maintains its own polarity with a fixed strength and attempts to move into the polarity direction and (2) objects interact with each other through the abovementioned exclusions. Based on this model, we look at the difference of consequences and combinatory effects of these two kinds of exclusions. Furthermore, we calculate the polar order of polarity directions without an external directional bias. Our results suggest that the combination of these two kinds of exclusions leads to effectively inelastic scattering of two objects, which eventually gives rise to global polar ordering. We also find that the traveling band can arise by this mechanism of alignment at the intermediate density, as generally seen in collective motion with polar alignment and investigated in various earlier works. Characteristics of transitions among disordered, traveling band, and homogeneously ordered states of the presented model are investigated, and their similarities and differences with those given by the explicit alignment interaction are discussed.

Three-body correlations and conditional forces in suspensions of active hard disks

Self-propelled Brownian particles show rich out-of-equilibrium physics, for instance, the motility-induced phase separation (MIPS). While decades of studying the structure of liquids have established a deep understanding of passive systems, not much is known about correlations in active suspensions. In this work we derive an approximate analytic theory for three-body correlations and forces in systems of active Brownian disks starting from the many-body Smoluchowski equation. We use our theory to predict the conditional forces that act on a tagged particle and their dependence on the propulsion speed of self-propelled disks. We identify preferred directions of these forces in relation to the direction of propulsion and the positions of the surrounding particles. We further relate our theory to the effective swimming speed of the active disks, which is relevant for the physics of MIPS. To test and validate our theory, we additionally run particle-resolved computer simulations, for which we explicitly calculate the three-body forces. In this context, we discuss the modeling of active Brownian swimmers with nearly hard interaction potentials. We find very good agreement between our simulations and numerical solutions of our theory, especially for the nonequilibrium pair-distribution function. For our analytical results, we carefully discuss their range of validity in the context of the different levels of approximation we applied. This discussion allows us to study the individual contribution of particles to three-body forces and to the emerging structure. Thus, our work sheds light on the collective behavior, provides the basis for further studies of correlations in active suspensions, and makes a step towards an emerging liquid state theory.

Do hydrodynamically assisted binary collisions lead to orientational ordering of microswimmers?

We have investigated the onset of collective motion in systems of model microswimmers, by performing a comprehensive analysis of the binary collision dynamics using three-dimensional direct numerical simulations (DNS) with hydrodynamic interactions. From this data, we have constructed a simplified binary collision model (BCM) which accurately reproduces the collective behavior obtained from the DNS for most cases. Thus, we show that global alignment can mostly arise solely from binary collisions. Although the agreement between both models (DNS and BCM) is not perfect, the parameter range in which notable differences appear is also that for which strong density fluctuations are present in the system (where pseudo-sound mound can be observed (N. Oyama et al., Phys. Rev. E 93, 043114 (2016))).

The emergence and transient behaviour of collective motion in active filament systems

Most living systems, ranging from animal flocks, self-motile microorganisms to the cytoskeleton rely on self-organization processes to perform their own specific function. Despite its importance, the general understanding of how individual active constituents initiate the intriguing pattern formation phenomena on all these different length scales still remains elusive. Here, using a high density actomyosin motility assay system, we show that the observed collective motion arises from a seeding process driven by enhanced acute angle collisions. Once a critical size is reached, the clusters coarsen into high and low density phases each with fixed filament concentrations. The steady state is defined by a balance of collision induced randomization and alignment effects of the filaments by multi-filament collisions within ordered clusters.

Self-organization is observed in cytoskeletal systems but emergence of order from disorder is poorly understood. Using a high density actomyosin system, the authors capture the transition from disorder to order, which is driven by enhanced alignment effects caused by increase in multi-filament collisions.

Cell-alignment patterns in the collective migration of cells with polarized adhesion

Dictyostelium discoideum (Dd) utilizes inhomogeneities in the distribution of cell-cell adhesion molecules on cell membranes for collective cell migration. A simple example of an inhomogeneity is a front-side (leading-edge) polarization in the distribution at the early streaming stage. Experiments have shown that the polarized cell-cell adhesion induces side-by-side contact between cells [Beug et al., Nature (London) 274, 445 (1978)NATUAS0028-083610.1038/274445a0]. This result is counterintuitive, as one would expect cells to align front to front in contact with each other on the basis of front-side polarization. In this work, we theoretically examine whether front-side polarization induces side-by-side contact in collective cell migration. We construct a model for expressing cells with this polarization based on the two-dimensional cellular Potts model. By a numerical simulation with this model, we find cell-cell alignment wherein cells form lateral arrays with side-by-side contacts as observed in the experiments.

Order-disorder transition in repulsive self-propelled particle systems

We study the collective dynamics of repulsive self-propelled particles. The particles are governed by coupled equations of motion that include polar self-propulsion, damping of velocity and of polarity, repulsive particle-particle interaction, and deterministic dynamics. Particle dynamics simulations show that the collective coherent motion with large density fluctuations spontaneously emerges from a disordered, isotropic state. In the parameter region where the rotational damping of polarity is strong, the system undergoes an abrupt shift to the absorbing ordered state after a waiting period in the metastable disordered state. In order to obtain a simple understanding of the mechanism underlying the collective behavior, we analyze the binary particle scattering process. We show that this approach correctly predicts the order-disorder transition at a dilute limit. The same approach is expanded for finite densities, although it disagrees with the result from many-particle simulations due to many-body correlations and density fluctuations.

Polar Pattern Formation in Driven Filament Systems Require Non-Binary Particle Collisions

Living matter has the extraordinary ability to behave in a concerted manner, which is exemplified throughout nature ranging from the self-organisation of the cytoskeleton to flocks of animals [1-4]. The microscopic dynamics of constituents have been linked to the system's meso- or macroscopic behaviour in silico via the Boltzmann equation for propelled particles [5-10]. Thereby, simplified binary collision rules between the constituents had to be assumed due to the lack of experimental data. We report here experimentally determined binary collision statistics by studying the recently introduced molecular system, the high density actomyosin motility assay [11-13]. We demonstrate that the alignment effect of the binary collision statistics is too weak to account for the observed ordering transition. The transition density for polar pattern formation decreases quadratically with filament length, which indicates that multi-filament collisions drive the observed ordering phenomenon and that a gas-like picture cannot explain the transition of the system to polar order. The presented findings demonstrate that the unique properties of biological active matter systems require a description that goes well beyond a gas-like picture developed in the framework of kinetic theories.

Formation and dissolution of bacterial colonies

Many organisms form colonies for a transient period of time to withstand environmental pressure. Bacterial biofilms are a prototypical example of such behavior. Despite significant interest across disciplines, physical mechanisms governing the formation and dissolution of bacterial colonies are still poorly understood. Starting from a kinetic description of motile and interacting cells we derive a hydrodynamic equation for their density on a surface, where most of the kinetic coefficients are estimated from experimental data for N. gonorrhoeae bacteria. We use it to describe the formation of multiple colonies with sizes consistent with experimental observations. Finally, we show how the changes in the cell-to-cell interactions lead to the dissolution of the bacterial colonies. The successful application of kinetic theory to a complex far from equilibrium system such as formation and dissolution of living bacterial colonies potentially paves the way for the physical quantification of the initial stages of biofilm formation.

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.

Tailoring the interactions between self-propelled bodies

We classify the interactions between self-propelled particles moving at a constant speed from symmetry considerations. We establish a systematic expansion for the two-body forces in the spirit of a multipolar expansion. This formulation makes it possible to rationalize most of the models introduced so far within a common framework. We distinguish between three classes of physical interactions: i) potential forces, ii) inelastic collisions and iii) non-reciprocal interactions involving polar or nematic alignment with an induced field. This framework provides simple design rules for the modeling and the fabrication of self-propelled bodies interacting via physical interactions. A class of possible interactions that should yield new phases of active matter is highlighted.

Defect-mediated phase transitions in active soft matter

How do topological defects affect the degree of order in active matter? To answer this question we investigate an agent-based model of self-propelled particles, which accounts for polar alignment and short-ranged repulsive interactions. For strong alignment forces we find collectively moving polycrystalline states with fluctuating networks of grain boundaries. In the regime where repulsive forces dominate, the fluctuations generated by the active system give rise to quasi-long-range transitional order, but-unlike the thermal system-without creating topological defects.

Controlling active self-assembly through broken particle-shape symmetry

Many structural properties of conventional passive materials are known to arise from the symmetries of their microscopic constituents. By contrast, it is largely unclear how the interplay between particle shape and self-propulsion controls the meso- and macroscale behavior of active matter. Here we use large-scale simulations of homo- and heterogeneous self-propelled particle systems to identify generic effects of broken particle-shape symmetry on collective motion. We find that even small violations of fore-aft symmetry lead to fundamentally different collective behaviors, which may facilitate demixing of differently shaped species as well as the spontaneous formation of stable microrotors. These results suggest that variation of particle shape yields robust physical mechanisms to control self-assembly of active matter, with possibly profound implications for biology and materials design.