Clustering and heterogeneous dynamics in a kinetic Monte Carlo model of self-propelled hard disks

Emergent structures and dynamics of cell colonies by contact inhibition of locomotion

Cells in tissues can organize into a broad spectrum of structures according to their function. Drastic changes of organization, such as epithelial-mesenchymal transitions or the formation of spheroidal aggregates, are often associated either to tissue morphogenesis or to cancer progression. Here, we study the organization of cell colonies by means of simulations of self-propelled particles with generic cell-like interactions. The interplay between cell softness, cell-cell adhesion, and contact inhibition of locomotion (CIL) yields structures and collective dynamics observed in several existing tissue phenotypes. These include regular distributions of cells, dynamic cell clusters, gel-like networks, collectively migrating monolayers, and 3D aggregates. We give analytical predictions for transitions between noncohesive, cohesive, and 3D cell arrangements. We explicitly show how CIL yields an effective repulsion that promotes cell dispersal, thereby hindering the formation of cohesive tissues. Yet, in continuous monolayers, CIL leads to collective cell motion, ensures tensile intercellular stresses, and opposes cell extrusion. Thus, our work highlights the prominent role of CIL in determining the emergent structures and dynamics of cell colonies.

Classical Nucleation Theory Description of Active Colloid Assembly

Nonaligning self-propelled particles with purely repulsive excluded volume interactions undergo athermal motility-induced phase separation into a dilute gas and a dense cluster phase. Here, we use enhanced sampling computational methods and analytic theory to examine the kinetics of formation of the dense phase. Despite the intrinsically nonequilibrium nature of the phase transition, we show that the kinetics can be described using an approach analogous to equilibrium classical nucleation theory, governed by an effective free energy of cluster formation with identifiable bulk and surface terms. The theory captures the location of the binodal, nucleation rates as a function of supersaturation, and the cluster size distributions below the binodal, while discrepancies in the metastable region reveal additional physics about the early stages of active crystal formation. The success of the theory shows that a framework similar to equilibrium thermodynamics can be obtained directly from the microdynamics of an active system, and can be used to describe the kinetics of evolution toward nonequilibrium steady states.

Crystallization of Self-Propelled Hard Discs

We experimentally study the crystallization of a monolayer of vibrated discs with a built-in polar asymmetry, a model system of active liquids, and contrast it with that of vibrated isotropic discs. Increasing the packing fraction ϕ, the quasicontinuous crystallization reported for isotropic discs is replaced by a transition, or a crossover, towards a "self-melting" crystal. Starting from the liquid phase and increasing the packing fraction, clusters of dense hexagonal-ordered packed discs spontaneously form, melt, split, and merge, leading to a highly intermittent and heterogeneous dynamics. For a packing fraction larger than ϕ^{*}, a few large clusters span the system size. The cluster size distribution is monotonically decreasing for ϕ<ϕ^{*}, nonmonotonic for ϕ>ϕ^{*}, and is a power law at the transition. The system is, however, never dynamically arrested. The clusters permanently melt from place to place, forming droplets of an active liquid which rapidly propagate across the system. This self-melting crystalline state subsists up to the highest possible packing fraction, questioning the stability of the crystal for active discs unless it is at ordered close packing.

Coarsening and clustering in run-and-tumble dynamics with short-range exclusion

The emergence of clustering and coarsening in crowded ensembles of self-propelled agents is studied using a lattice model in one dimension. The persistent exclusion process, where particles move at directions that change randomly at a low tumble rate α, is extended allowing sites to be occupied by more than one particle, with a maximum n_{max} per site. Three phases are distinguished. For n_{max}=1 a gas of clusters form, with sizes distributed exponentially and no coarsening takes place. For n_{max}≥3 and small values of α, coarsening takes place and few large clusters appear, with a large fraction of the total number of particles in them. In the same range of n_{max} but for larger values of α, a gas phase where a negligible fraction of particles takes part of clusters. Finally, n_{max}=2 corresponds to a crossover phase. The character of the transitions between phases is studied extending the model to allow n_{max} to take real values and jumps to an occupied site are probabilistic. The transition from the gas of clusters to the coarsening phase is continuous and the mass of the large clusters grows continuously when varying the maximum occupancy, and the crossover found corresponds to values close to the transition. The second transition, from the coarsening to the gaseous phase, can be either continuous or discontinuous depending on the parameters, with a critical point separating both cases.

Active cage model of glassy dynamics

We build up a phenomenological picture in terms of the effective dynamics of a tracer confined in a cage experiencing random hops to capture some characteristics of glassy systems. This minimal description exhibits scale invariance properties for the small-displacement distribution that echo experimental observations. We predict the existence of exponential tails as a crossover between two Gaussian regimes. Moreover, we demonstrate that the onset of glassy behavior is controlled only by two dimensionless numbers: the number of hops occurring during the relaxation of the particle within a local cage and the ratio of the hopping length to the cage size.

Binary Mixtures of Particles with Different Diffusivities Demix

The influence of size differences, shape, mass, and persistent motion on phase separation in binary mixtures has been intensively studied. Here we focus on the exclusive role of diffusivity differences in binary mixtures of equal-sized particles. We find an effective attraction between the less diffusive particles, which are essentially caged in the surrounding species with the higher diffusion constant. This effect leads to phase separation for systems above a critical size: A single close-packed cluster made up of the less diffusive species emerges. Experiments for testing our predictions are outlined.

Clustering and phase behaviour of attractive active particles with hydrodynamics

We simulate clustering, phase separation and hexatic ordering in a monolayered suspension of active squirming disks subject to an attractive Lennard-Jones-like pairwise interaction potential, taking hydrodynamic interactions between the particles fully into account. By comparing the hydrodynamic case with counterpart simulations for passive and active Brownian particles, we elucidate the relative roles of self-propulsion, interparticle attraction, and hydrodynamic interactions in determining clustering and phase behaviour. Even in the presence of an attractive potential, we find that hydrodynamic interactions strongly suppress the motility induced phase separation that might a priori have been expected in a highly active suspension. Instead, we find only a weak tendency for the particles to form stringlike clusters in this regime. At lower activities we demonstrate phase behaviour that is broadly equivalent to that of the counterpart passive system at low temperatures, characterized by regimes of gas-liquid, gas-solid and liquid-solid phase coexistence. In this way, we suggest that a dimensionless quantity representing the level of activity relative to the strength of attraction plays the role of something like an effective non-equilibrium temperature, counterpart to the (dimensionless) true thermodynamic temperature in the passive system. However there are also some important differences from the equilibrium case, most notably with regards the degree of hexatic ordering, which we discuss carefully.

Formation, compression and surface melting of colloidal clusters by active particles

We demonstrate with experiments and numerical simulations that the structure and dynamics of a suspension of passive particles is strongly altered by adding a very small (<1%) number of active particles. With increasing passive particle density, we observe first the formation of dynamic clusters comprised of passive particles being surrounded by active particles, then the merging and compression of these clusters, and eventually the local melting of crystalline regions by enclosed active particles.

Self-phoretic active particles interacting by diffusiophoresis: A numerical study of the collapsed state and dynamic clustering

Self-phoretic active colloids move and orient along self-generated chemical gradients by diffusiophoresis, a mechanism reminiscent of bacterial chemotaxis. In combination with the activity of the colloids, this creates effective repulsive and attractive interactions between particles depending on the sign of the translational and rotational diffusiophoretic parameters. A delicate balance of these interactions causes dynamic clustering and for overall strong effective attraction the particles collapse to one single cluster. Using Langevin dynamics simulations, we extend the state diagram of our earlier work (Phys. Rev. Lett. 112, 238303 (2014)) to regions with translational phoretic repulsion. With increasing repulsive strength, the collapsed cluster first starts to fluctuate strongly, then oscillates between a compact form and a colloidal cloud, and ultimately the colloidal cloud becomes static. The oscillations disappear if the phoretic interactions within compact clusters are not screened. We also study dynamic clustering at larger area fractions by exploiting cluster size distributions and mean cluster sizes. In particular, we identify the dynamic clustering 2 state as a signature of phoretic interactions. We analyze fusion and fission rate functions to quantify the kinetics of cluster formation and identify them as local signatures of phoretic interactions, since they can be measured on single clusters.

Rotational and translational diffusion in an interacting active dumbbell system

We study the dynamical properties of a two-dimensional ensemble of self-propelled dumbbells with only repulsive interactions. This model undergoes a phase transition between a homogeneous and a segregated phase and we focus on the former. We analyze the translational and rotational mean-square displacements in terms of the Péclet number, describing the relative role of active forces and thermal fluctuations, and of particle density. We find that the four distinct regimes of the translational mean-square displacement of the single active dumbbell survive at finite density for parameters that lead to a separation of time scales. We establish the Péclet number and density dependence of the diffusion constant in the last diffusive regime. We prove that the ratio between the diffusion constant and its value for the single dumbbell depends on temperature and active force only through the Péclet number at all densities explored. We also study the rotational mean-square displacement proving the existence of a rich behavior with intermediate regimes only appearing at finite density. The ratio of the rotational late-time diffusion constant and its vanishing density limit depends on the Péclet number and density only. At low Péclet number it is a monotonically decreasing function of density. At high Péclet number it first increases to reach a maximum and then decreases as a function of density. We interpret the latter result advocating the presence of large-scale fluctuations close to the transition, at large-enough density, that favor coherent rotation inhibiting, however, rotational motion for even larger packing fractions.

Multidimensional stationary probability distribution for interacting active particles

We derive the stationary probability distribution for a non-equilibrium system composed by an arbitrary number of degrees of freedom that are subject to Gaussian colored noise and a conservative potential. This is based on a multidimensional version of the Unified Colored Noise Approximation. By comparing theory with numerical simulations we demonstrate that the theoretical probability density quantitatively describes the accumulation of active particles around repulsive obstacles. In particular, for two particles with repulsive interactions, the probability of close contact decreases when one of the two particle is pinned. Moreover, in the case of isotropic confining potentials, the radial density profile shows a non trivial scaling with radius. Finally we show that the theory well approximates the “pressure” generated by the active particles allowing to derive an equation of state for a system of non-interacting colored noise-driven particles.

Effective interactions in active Brownian suspensions

Active colloids exhibit persistent motion, which can lead to motility-induced phase separation (MIPS). However, there currently exists no microscopic theory to account for this phenomenon. We report a first-principles theory, free of fit parameters, for active spherical colloids, which shows explicitly how an effective many-body interaction potential is generated by activity and how this can rationalize MIPS. For a passively repulsive system the theory predicts phase separation and pair correlations in quantitative agreement with simulation. For an attractive system the theory shows that phase separation becomes suppressed by moderate activity, consistent with recent experiments and simulations, and suggests a mechanism for reentrant cluster formation at high activity.

Active microrheology in active matter systems: Mobility, intermittency, and avalanches

We examine the mobility and velocity fluctuations of a driven particle moving through an active matter bath of self-mobile disks for varied density or area coverage and varied activity. We show that the driven particle mobility can exhibit nonmonotonic behavior that is correlated with distinct changes in the spatiotemporal structures that arise in the active media. We demonstrate that the probe particle velocity distributions exhibit specific features in the different dynamic regimes and identify an activity-induced uniform crystallization that occurs for moderate activity levels and is distinct from the previously observed higher activity cluster phase. The velocity distribution in the cluster phase has telegraph noise characteristics produced when the probe particle moves alternately through high-mobility areas that are in the gas state and low-mobility areas that are in the dense phase. For higher densities and large activities, the system enters what we characterize as an active jamming regime. Here the probe particle moves in intermittent jumps or avalanches that have power-law-distributed sizes that are similar to the avalanche distributions observed for nonactive disk systems near the jamming transition.

Dynamics of a homogeneous active dumbbell system

We analyze the dynamics of a two-dimensional system of interacting active dumbbells. We characterize the mean-square displacement, linear response function, and deviation from the equilibrium fluctuation-dissipation theorem as a function of activity strength, packing fraction, and temperature for parameters such that the system is in its homogeneous phase. While the diffusion constant in the last diffusive regime naturally increases with activity and decreases with packing fraction, we exhibit an intriguing nonmonotonic dependence on the activity of the ratio between the finite-density and the single-particle diffusion constants. At fixed packing fraction, the time-integrated linear response function depends nonmonotonically on activity strength. The effective temperature extracted from the ratio between the integrated linear response and the mean-square displacement in the last diffusive regime is always higher than the ambient temperature, increases with increasing activity, and, for small active force, monotonically increases with density while for sufficiently high activity it first increases and next decreases with the packing fraction. We ascribe this peculiar effect to the existence of finite-size clusters for sufficiently high activity and density at the fixed (low) temperatures at which we worked. The crossover occurs at lower activity or density the lower the external temperature. The finite-density effective temperature is higher (lower) than the single dumbbell one below (above) a crossover value of the Péclet number.

Absorbing phase transitions and dynamic freezing in running active matter systems

We examine a two-dimensional system of sterically repulsive interacting disks where each particle runs in a random direction. This system is equivalent to a run-and-tumble dynamics system in the limit where the run time is infinite. At low densities, we find a strongly fluctuating state composed of transient clusters. Above a critical density that is well below the density at which non-active particles would crystallize, the system can organize into a drifting quiescent or frozen state where the fluctuations are lost and large crystallites form surrounded by a small density of individual particles. Although all the particles are still moving, their paths form closed orbits. The average transient time to organize into the quiescent state diverges as a power law upon approaching the critical density from above. We compare our results to the random organization observed for periodically sheared systems that can undergo an absorbing transition from a fluctuating state to a dynamical non-fluctuating state. In the random organization studies, the system organizes to a state in which the particles no longer interact; in contrast, we find that the randomly running active matter organizes to a strongly interacting dynamically jammed state. We show that the transition to the frozen state is robust against a certain range of stochastic fluctuations. We also examine the effects of adding a small number of pinned particles to the system and find that the transition to the frozen state shifts to significantly lower densities and arises via the nucleation of faceted crystals centered at the obstacles.

Self-propelled particle in an external potential: existence of an effective temperature

We study a stationary state of a single self-propelled, athermal particle in linear and quadratic external potentials. The self-propulsion is modeled as a fluctuating internal driving force evolving according to the Ornstein-Uhlenbeck process, independently of the state of the particle. Without an external potential, in the long time limit, the self-propelled particle moving in a viscous medium performs diffusive motion, which allows one to identify an effective temperature. We show that in the presence of a linear external potential the stationary state distribution has an exponential form with the sedimentation length determined by the effective temperature of the free self-propelled particle. In the presence of a quadratic external potential the stationary state distribution has a Gaussian form. However, in general, this distribution is not determined by the effective temperature of the free self-propelled particle.