Generalized thermodynamics of phase equilibria in scalar active matter

Statistical mechanics of active Ornstein-Uhlenbeck particles

We study the statistical properties of active Ornstein-Uhlenbeck particles (AOUPs). In this simplest of models, the Gaussian white noise of overdamped Brownian colloids is replaced by a Gaussian colored noise. This suffices to grant this system the hallmark properties of active matter, while still allowing for analytical progress. We study in detail the steady-state distribution of AOUPs in the small persistence time limit and for spatially varying activity. At the collective level, we show AOUPs to experience motility-induced phase separation both in the presence of pairwise forces or due to quorum-sensing interactions. We characterize both the instability mechanism leading to phase separation and the resulting phase coexistence. We probe how, in the stationary state, AOUPs depart from their thermal equilibrium limit by investigating the emergence of ratchet currents and entropy production. In the small persistence time limit, we show how fluctuation-dissipation relations are recovered. Finally, we discuss how the emerging properties of AOUPs can be characterized from the dynamics of their collective modes.

Coexistence of active Brownian disks: van der Waals theory and analytical results

At thermal equilibrium, intensive quantities like temperature and pressure have to be uniform throughout the system, restricting inhomogeneous systems composed of different phases. The paradigmatic example is the coexistence of vapor and liquid, a state that can also be observed for active Brownian particles steadily driven away from equilibrium. Recently, a strategy has been proposed that allows to predict phase equilibria of active particles [Solon et al., Phys. Rev. E 97, 020602(R) (2018)2470-004510.1103/PhysRevE.97.020602]. Here we elaborate on this strategy and formulate it in the framework of a van der Waals theory for active disks. For a given equation of state, we derive the effective free energy analytically and show that it yields coexisting densities in very good agreement with numerical results. We discuss the interfacial tension and the relation to Cahn-Hilliard models.

Nonequilibrium grand-canonical ensemble built from a physical particle reservoir

We introduce a nonequilibrium grand-canonical ensemble defined by considering the stationary state of a driven system of particles put in contact with a particle reservoir. When an additivity assumption holds for the large deviation function of density, a chemical potential of the reservoir can be defined. The grand-canonical distribution then takes a form similar to the equilibrium one. At variance with equilibrium, though, the probability weight is "renormalized" by a contribution coming from the contact, with respect to the canonical probability weight of the isolated system. A formal grand-canonical potential can be introduced in terms of a scaled cumulant generating function, defined as the Legendre-Fenchel transform of the large deviation function of density. The role of the formal Legendre parameter can be played, physically, by the chemical potential of the reservoir when the latter can be defined, or by a potential energy difference applied between the system and the reservoir. Static fluctuation-response relations naturally follow from the large deviation structure. Some of the results are illustrated on two different explicit examples, a gas of noninteracting active particles and a lattice model of interacting particles.

Disorder-Induced Long-Ranged Correlations in Scalar Active Matter

We study the impact of quenched random potentials and torques on scalar active matter. Microscopic simulations reveal that motility-induced phase separation is replaced in two dimensions by an asymptotically homogeneous phase with anomalous long-ranged correlations and nonvanishing steady-state currents. Using a combination of phenomenological models and a field-theoretical treatment, we show the existence of a lower-critical dimension d_{c}=4, below which phase separation is only observed for systems smaller than an Imry-Ma length scale. We identify a weak-disorder regime in which the structure factor scales as S(q)∼1/q^{2}, which accounts for our numerics. In d=2, we predict that, at larger scales, the behavior should cross over to a strong-disorder regime. In d>2, these two regimes exist separately, depending on the strength of the potential.

Phase Separation and Multibody Effects in Three-Dimensional Active Brownian Particles

Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behavior of equilibrium fluids with short-ranged pairwise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multibody effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of n-body effective interactions extracted from the configurational structure.

The "isothermal" compressibility of active matter

We demonstrate that the mechanically defined "isothermal" compressibility behaves as a thermodynamic-like response function for suspensions of active Brownian particles. The compressibility computed from the active pressure-a combination of the collision and unique swim pressures-is capable of predicting the critical point for motility induced phase separation, as expected from the mechanical stability criterion. We relate this mechanical definition to the static structure factor via an active form of the thermodynamic compressibility equation and find the two to be equivalent, as would be the case for equilibrium systems. This equivalence indicates that compressibility behaves like a thermodynamic response function, even when activity is large. Finally, we discuss the importance of the phase interface when defining an active chemical potential. Previous definitions of the active chemical potential are shown to be accurate above the critical point but breakdown in the coexistence region. Inclusion of the swim pressure in the mechanical compressibility definition suggests that the interface is essential for determining phase behavior.

Lamellar to Micellar Phases and Beyond: When Tactic Active Systems Admit Free Energy Functionals

We consider microscopic models of active particles whose velocities, rotational diffusivities, and tumbling rates depend on the gradient of a local field that is either externally imposed or depends on all particle positions. Despite the fundamental differences between active and passive dynamics at the microscopic scale, we show that a large class of such tactic active systems admit fluctuating hydrodynamics equivalent to those of interacting Brownian colloids in equilibrium. We exploit this mapping to show how taxis may lead to the lamellar and micellar phases observed for soft repulsive colloids. In the context of chemotaxis, we show how the competition between chemoattractant and chemorepellent may lead to a bona fide equilibrium liquid-gas phase separation in which a loss of thermodynamic stability of the fluid signals the onset of a chemotactic collapse.

Self-Organized Critical Coexistence Phase in Repulsive Active Particles

We revisit motility-induced phase separation in two models of active particles interacting by pairwise repulsion and uncover new qualitative features: the resulting dense phase contains gas bubbles distributed algebraically up to a typically extremely large cutoff scale. At large enough system size and/or global density, all the gas may be contained inside the bubbles, at which point the system is microphase separated with a finite cutoff bubble scale. We further observe that the ordering is clearly anomalous, with different dynamics for the coarsening of the dense phase and of the gas bubbles. This self-organized critical phenomenology is reproduced by a "reduced bubble model" that implements the basic idea of reverse Ostwald ripening put forward in Tjhung et al. [Phys. Rev. X 8, 031080 (2018)PRXHAE2160-330810.1103/PhysRevX.8.031080].

Surface Tensions between Active Fluids and Solid Interfaces: Bare vs Dressed

We analyze the surface tension exerted at the interface between an active fluid and a solid boundary in terms of tangential forces. Focusing on active systems known to possess an equation of state for the pressure, we show that interfacial forces are of a more complex nature. Using a number of macroscopic setups, we show that the surface tension is a combination of an equation-of-state abiding part and of setup-dependent contributions. The latter arise from generic setup-dependent steady currents which "dress" the measurement of the "bare" surface tension. The former shares interesting properties with its equilibrium counterpart, and can be used to generalize the Young-Laplace law to active systems. Finally, we show how a suitably designed probe can directly access this bare surface tension, which can also be computed using a generalized virial formula.

Stealth Entropy Production in Active Field Theories near Ising Critical Points

We address the steady-state entropy production rate (EPR) of active scalar ϕ^{4} theories, which lack time-reversal symmetry, close to a phase-separation critical point. We consider both nonconserved (model A) and conserved (model B) dynamics at Gaussian level, and also address the former at leading order in ε=4-d. In each case, activity is irrelevant in the RG sense: the active model lies in the same (dynamic Ising) universality class as its time-reversible counterpart. Hence one might expect that activity brings no new critical behavior. Here we show instead that, on approach to criticality in these models, the singular part of the EPR per (diverging) spacetime correlation volume either remains finite or itself diverges. A nontrivial critical scaling for entropy production thus ranks among universal dynamic Ising-class properties.

Different-shaped micro-objects driven by active particle aggregations

The dynamics of passive micro-objects in an active bath has been receiving much attention. However, the influence of the shapes of micro-objects remains unclear. Here, we use 2D simulation to investigate the interaction between active Brownian particles and different-shaped passive micro-objects. We show that active particles accumulate around micro-objects and self-assemble into living aggregations at a high active velocity and high volume fraction. The shapes of micro-objects affect the distributions of the aggregations. In turn, the different distribution of aggregations influences the motion of micro-objects and induces abnormal diffusive behaviors. We further demonstrate that polar distributed aggregations at a high active velocity and the inhibition of the active bath at a low active velocity induce the counterintuitive anisotropic enhanced diffusion of rods, and the steric interaction between active particles induces the reverse translation-rotation coupled diffusion of chevrons.

The 2020 motile active matter roadmap

Activity and autonomous motion are fundamental in living and engineering systems. This has stimulated the new field of 'active matter' in recent years, which focuses on the physical aspects of propulsion mechanisms, and on motility-induced emergent collective behavior of a larger number of identical agents. The scale of agents ranges from nanomotors and microswimmers, to cells, fish, birds, and people. Inspired by biological microswimmers, various designs of autonomous synthetic nano- and micromachines have been proposed. Such machines provide the basis for multifunctional, highly responsive, intelligent (artificial) active materials, which exhibit emergent behavior and the ability to perform tasks in response to external stimuli. A major challenge for understanding and designing active matter is their inherent nonequilibrium nature due to persistent energy consumption, which invalidates equilibrium concepts such as free energy, detailed balance, and time-reversal symmetry. Unraveling, predicting, and controlling the behavior of active matter is a truly interdisciplinary endeavor at the interface of biology, chemistry, ecology, engineering, mathematics, and physics. The vast complexity of phenomena and mechanisms involved in the self-organization and dynamics of motile active matter comprises a major challenge. Hence, to advance, and eventually reach a comprehensive understanding, this important research area requires a concerted, synergetic approach of the various disciplines. The 2020 motile active matter roadmap of Journal of Physics: Condensed Matter addresses the current state of the art of the field and provides guidance for both students as well as established scientists in their efforts to advance this fascinating area.

Clustering and phase separation in mixtures of dipolar and active particles

The self-assembly of colloidal particles in dynamic environments has become an important field of study because of potential applications in fabricating out-of-equilibrium materials. We investigate the phase behavior of mixtures of passive dipolar colloids and active soft spheres using Brownian dynamics simulations in two dimensions. The phase behaviors exhibited include dipolar percolated network, dipolar string-fluid, isotropic fluid, and a phase-separated state. We find that the clustering of dipolar colloids is enhanced in the presence of slow-moving active particles compared to the clustering of dipolar particles mixed with passive particles. When the active particle motility is high, the chains of dipolar particles are either broken into short chains or pushed into dense clusters. Motility-induced phase separation into dense and dilute phases is also present. The area fraction of particles in the dilute phase increases as the fraction of active particles in the system decreases, while the area fraction of particles in the dense phase remains constant. Our findings are relevant to the development of reconfigurable self-assembled materials.

Collective forces in scalar active matter

Large-scale collective behavior in suspensions of active particles can be understood from the balance of statistical forces emerging beyond the direct microscopic particle interactions. Here we review some aspects of the collective forces that can arise in suspensions of self-propelled active Brownian particles: wall forces under confinement, interfacial forces, and forces on immersed bodies mediated by the suspension. Even for non-aligning active particles, these forces are intimately related to a non-uniform polarization of particle orientations induced by walls and bodies, or inhomogeneous density profiles. We conclude by pointing out future directions and promising areas for the application of collective forces in synthetic active matter, as well as their role in living active matter.

Resonant osmosis across active switchable membranes

To overcome the traditional paradigm of filtration, where separation is essentially performed upon steric sieving principles, we explore the concept of dynamic osmosis through active membranes. A partially permeable membrane presents a time-tuneable feature that changes the effective pore interaction with the solute and thus actively changes permeability with time. In general, we find that slow flickering frequencies effectively decrease the osmotic pressure and large flickering frequencies do not change it. In the presence of an asymmetric membrane, we find a resonant frequency where pumping of the solute is performed and can be analyzed in terms of ratchet transport. We discuss and highlight the properties of this resonant osmotic transport. Furthermore, we show that dynamic osmosis allows us to pump the solute at the nanoscale using less energy than reverse osmosis. This opens new possibilities to build advanced filtration devices and design artificial ionic machinery.

Microscopic origins of the swim pressure and the anomalous surface tension of active matter

The unique pressure exerted by active particles-the "swim" pressure-has proven to be a useful quantity in explaining many of the seemingly confounding behaviors of active particles. However, its use has also resulted in some puzzling findings including an extremely negative surface tension between phase separated active particles. Here, we demonstrate that this contradiction stems from the fact that the swim pressure is not a true pressure. At a boundary or interface, the reduction in particle swimming generates a net active force density-an entirely self-generated body force. The pressure at the boundary, which was previously identified as the swim pressure, is in fact an elevated (relative to the bulk) value of the traditional particle pressure that is generated by this interfacial force density. Recognizing this unique mechanism for stress generation allows us to define a much more physically plausible surface tension. We clarify the utility of the swim pressure as an "equivalent pressure" (analogous to those defined from electrostatic and gravitational body forces) and the conditions in which this concept can be appropriately applied.

Linear Response Theory and Green-Kubo Relations for Active Matter

We address the question of how interacting active systems in a nonequilibrium steady state respond to an external perturbation. We establish an extended fluctuation-dissipation theorem for active Brownian particles (ABP), which highlights the role played by the local violation of detailed balance due to activity. By making use of a Markovian approximation we derive closed Green-Kubo expressions for the diffusivity and mobility of ABP and quantify the deviations from the Stokes-Einstein relation. We compute the linear response function to an external force using unperturbed simulations of ABP and compare the results with the analytical predictions of the transport coefficients. Our results show the importance of the interplay between activity and interactions in the departure from equilibrium linear response.

Motility-Induced Temperature Difference in Coexisting Phases

Unlike in thermodynamic equilibrium where coexisting phases always have the same temperature, here we show that systems comprising "active" self-propelled particles can self-organize into two coexisting phases at different kinetic temperatures, which are separated from each other by a sharp and persistent temperature gradient. Contrasting previous studies that have focused on overdamped descriptions of active particles, we show that a "hot-cold coexistence" occurs if and only if accounting for inertia, which is significant, e.g., in activated dusty plasmas, microflyers, whirling fruits, or beetles at interfaces. Our results exemplify a route to use active particles to create a self-sustained temperature gradient across coexisting phases. This phenomenon is fundamentally beyond equilibrium physics and is accompanied by a slow coarsening law with an exponent significantly smaller than the universal 1/3 exponent seen in both equilibrium systems and overdamped active Brownian particles.

Phase coexistence of active Brownian particles

We investigate motility-induced phase separation of active Brownian particles, which are modeled as purely repulsive spheres that move due to a constant swim force with freely diffusing orientation. We develop on the basis of power functional concepts an analytical theory for nonequilibrium phase coexistence and interfacial structure. Theoretical predictions are validated against Brownian dynamics computer simulations. We show that the internal one-body force field has four nonequilibrium contributions: (i) isotropic drag and (ii) interfacial drag forces against the forward motion, (iii) a superadiabatic spherical pressure gradient, and (iv) the quiet life gradient force. The intrinsic spherical pressure is balanced by the swim pressure, which arises from the polarization of the free interface. The quiet life force opposes the adiabatic force, which is due to the inhomogeneous density distribution. The balance of quiet life and adiabatic forces determines bulk coexistence via equality of two bulk state functions, which are independent of interfacial contributions. The internal force fields are kinematic functionals which depend on density and current but are independent of external and swim forces, consistent with power functional theory. The phase transition originates from nonequilibrium repulsion, with the agile gas being more repulsive than the quiet liquid.

Nonequilibrium chemical potentials of steady-state lattice gas models in contact: A large-deviation approach

We introduce a general framework to describe the stationary state of two driven systems exchanging particles or mass through a contact, in a slow exchange limit. The definition of chemical potentials for the systems in contact requires that the large-deviation function describing the repartition of mass between the two systems is additive, in the sense of being a sum of contributions from each system. We show that this additivity property does not hold for an arbitrary contact dynamics, but is satisfied on condition that a macroscopic detailed balance condition holds at contact, and that the coarse-grained contact dynamics satisfies a factorization property. However, the nonequilibrium chemical potentials of the systems in contact keep track of the contact dynamics, and thus do not obey an equation of state. These nonequilibrium chemical potentials can be related either to the equilibrium chemical potential, or to the nonequilibrium chemical potential of the isolated systems. Results are applied both to an exactly solvable driven lattice gas model and to the Katz-Lebowitz-Spohn model using a numerical procedure to evaluate the chemical potential. The breaking of the additivity property is also illustrated on the exactly solvable model.

Reduced Thermodynamic Description of Phase Separation in a Quasi-One-Dimensional Granular Gas

We describe simulations of a quasi-one-dimensional, vibrated granular gas which exhibits an apparent phase separation into a liquidlike phase and a gaslike phase. In thermal equilibrium, such a phase separation in one dimension is prohibited by entropic considerations. We propose that the granular gas minimizes a function of the conserved mechanical variables alone: the particle number and volume. Simulations in small cells can be used to extract the equation of state and predict the coexisting pressure and densities, as confirmation of the minimization principle. Fluctuations in the system manifest themselves as persistent density waves but they do not destroy the phase-separated state.

Active nematic-isotropic interfaces in channels

We use numerical simulations to investigate the hydrodynamic behavior of the interface between nematic (N) and isotropic (I) phases of a confined active liquid crystal. At low activities, a stable interface with constant shape and velocity is observed separating the two phases. For nematics in homeotropic channels, the velocity of the interface at the NI transition increases from zero (i) linearly with the activity for contractile systems and (ii) quadratically for extensile ones. Interestingly, the nematic phase expands for contractile systems while it contracts for extensile ones, as a result of the active forces at the interface. Since both activity and temperature affect the stability of the nematic, for active nematics in the stable regime the temperature can be tuned to observe static interfaces, providing an operational definition for the coexistence of active nematic and isotropic phases. At higher activities, beyond the stable regime, an interfacial instability is observed for extensile nematics. In this regime defects are nucleated at the interface and move away from it. The dynamics of these defects is regular and persists asymptotically for a finite range of activities. We used an improved hybrid model of finite differences and the lattice Boltzmann method with a multi-relaxation-time collision operator, the accuracy of which allowed us to characterize the dynamics of the distinct interfacial regimes.

Critical Motility-Induced Phase Separation Belongs to the Ising Universality Class

A collection of self-propelled particles with volume exclusion interactions can exhibit the phenomenology of a gas-liquid phase separation, known as motility-induced phase separation (MIPS). The nonequilibrium nature of the system is fundamental to the phase transition; however, it is unclear whether MIPS at criticality contributes a novel universality class to nonequilibrium physics. We demonstrate here that this is not the case by showing that a generic critical MIPS belongs to the Ising universality class with conservative dynamics.

Systematic extension of the Cahn-Hilliard model for motility-induced phase separation

We consider a continuum model for motility-induced phase separation (MIPS) of active Brownian particles (ABP) (J. Chem. Phys. 142, 224149 (2015)). Using a recently introduced perturbative analysis (Phys. Rev. E 98, 020604(R) (2018)), we show that this continuum model reduces to the classic Cahn-Hilliard (CH) model near the onset of MIPS. This makes MIPS another example of the so-called active phase separation. We further introduce a generalization of the perturbative analysis to the next higher order. This results in a generic higher-order extension of the CH model for active phase separation. Our analysis establishes the mathematical link between the basic mean-field ABP model on the one hand, and the leading order and extended CH models on the other hand. Comparing numerical simulations of the three models, we find that the leading-order CH model agrees nearly perfectly with the full continuum model near the onset of MIPS. We also give estimates of the control parameter beyond which the higher-order corrections become relevant and compare the extended CH model to recent phenomenological models.

Local stress and pressure in an inhomogeneous system of spherical active Brownian particles

The stress of a fluid on a confining wall is given by the mechanical wall forces, independent of the nature of the fluid being passive or active. At thermal equilibrium, an equation of state exists and stress is likewise obtained from intrinsic bulk properties; even more, stress can be calculated locally. Comparable local descriptions for active systems require a particular consideration of active forces. Here, we derive expressions for the stress exerted on a local volume of a systems of spherical active Brownian particles (ABPs). Using the virial theorem, we obtain two identical stress expressions, a stress due to momentum flux across a hypothetical plane, and a bulk stress inside of the local volume. In the first case, we obtain an active contribution to momentum transport in analogy to momentum transport in an underdamped passive system, and we introduce an active momentum. In the second case, a generally valid expression for the swim stress is derived. By simulations, we demonstrate that the local bulk stress is identical to the wall stress of a confined system for both, non-interacting ABPs as well as ABPs with excluded-volume interactions. This underlines the existence of an equation of state for a system of spherical ABPs. Most importantly, our calculations demonstrated that active stress is not a wall (boundary) effect, but is caused by momentum transport. We demonstrate that the derived stress expression permits the calculation of the local stress in inhomogeneous systems of ABPs.

Optimizing active work: Dynamical phase transitions, collective motion, and jamming

Active work measures how far the local self-forcing of active particles translates into real motion. Using population Monte Carlo methods, we investigate large deviations in the active work for repulsive active Brownian disks. Minimizing the active work generically results in dynamical arrest; in contrast, despite the lack of aligning interactions, trajectories of high active work correspond to a collectively moving, aligned state. We use heuristic and analytic arguments to explain the origin of dynamical phase transitions separating the arrested, typical, and aligned regimes.

Swimming to Stability: Structural and Dynamical Control via Active Doping

External fields can decidedly alter the free energy landscape of soft materials and can be exploited as a powerful tool for the assembly of targeted nanostructures and colloidal materials. Here, we use computer simulations to demonstrate that nonequilibrium internal fields or forces-forces that are generated by driven components within a system-in the form of active particles can precisely modulate the dynamical free energy landscape of a model soft material, a colloidal gel. Embedding a small fraction of active particles within a gel can provide a unique pathway for the dynamically frustrated network to circumvent the kinetic barriers associated with reaching a lower free energy state through thermal fluctuations alone. Moreover, by carefully tuning the active particle properties (the propulsive swim force and persistence length) in comparison to those of the gel, the active particles may induce depletion-like forces between the constituent particles of the gel despite there being no geometric size asymmetry between the particles. These resulting forces can rapidly push the system toward disparate regions of phase space. Intriguingly, the state of the material can be altered by tuning macroscopic transport properties such as the solvent viscosity. Our findings highlight the potential wide-ranging structural and kinetic control facilitated by varying the dynamical properties of a remarkably small fraction of driven particles embedded in a host material.

Full Phase Diagram of Active Brownian Disks: From Melting to Motility-Induced Phase Separation

We establish the complete phase diagram of self-propelled hard disks in two spatial dimensions from the analysis of the equation of state and the statistics of local order parameters. The equilibrium melting scenario is maintained at small activities, with coexistence between active liquid and hexatic order, followed by a proper hexatic phase, and a further transition to an active solid. As activity increases, the emergence of hexatic and solid order is shifted towards higher densities. Above a critical activity and for a certain range of packing fractions, the system undergoes motility-induced phase separation and demixes into low and high density phases; the latter can be either disordered (liquid) or ordered (hexatic or solid) depending on the activity.

Isotropic-nematic transition of self-propelled rods in three dimensions

Using overdamped Brownian dynamics simulations we investigate the isotropic-nematic (IN) transition of self-propelled rods in three spatial dimensions. For two well-known model systems (Gay-Berne potential and hard spherocylinders) we find that turning on activity moves to higher densities the phase boundary separating an isotropic phase from a (nonpolar) nematic phase. This active IN phase boundary is distinct from the boundary between isotropic and polar-cluster states previously reported in two-dimensional simulation studies and, unlike the latter, is not sensitive to the system size. We thus identify a generic feature of anisotropic active particles in three dimensions.

Exact Hydrodynamic Description of Active Lattice Gases

We introduce lattice gas models of active matter systems whose coarse-grained "hydrodynamic" description can be derived exactly. We illustrate our approach by considering two systems exhibiting two of the most studied collective behaviors in active matter: the motility-induced phase separation and the transition to collective motion. In both cases, we derive coupled partial differential equations describing the dynamics of the local density and polarization fields and show how they quantitatively predict the emerging properties of the macroscopic lattice gases.

Collective cell migration without proliferation: density determines cell velocity and wave velocity

Collective cell migration contributes to embryogenesis, wound healing and tumour metastasis. Cell monolayer migration experiments help in understanding what determines the movement of cells far from the leading edge. Inhibiting cell proliferation limits cell density increase and prevents jamming; we observe long-duration migration and quantify space-time characteristics of the velocity profile over large length scales and time scales. Velocity waves propagate backwards and their frequency depends only on cell density at the moving front. Both cell average velocity and wave velocity increase linearly with the cell effective radius regardless of the distance to the front. Inhibiting lamellipodia decreases cell velocity while waves either disappear or have a lower frequency. Our model combines conservation laws, monolayer mechanical properties and a phenomenological coupling between strain and polarity: advancing cells pull on their followers, which then become polarized. With reasonable values of parameters, this model agrees with several of our experimental observations. Together, our experiments and model disantangle the respective contributions of active velocity and of proliferation in monolayer migration, explain how cells maintain their polarity far from the moving front, and highlight the importance of strain-polarity coupling and density in long-range information propagation.

Phase separation of self-propelled ballistic particles

Self-propelled particles phase-separate into coexisting dense and dilute regions above a critical density. The statistical nature of their stochastic motion lends itself to various theories that predict the onset of phase separation. However, these theories are ill-equipped to describe such behavior when noise becomes negligible. To overcome this limitation, we present a predictive model that relies on two density-dependent timescales: τ_{F}, the mean time particles spend between collisions; and τ_{C}, the mean lifetime of a collision. We show that only when τ_{F}<τ_{C} do collisions last long enough to develop a growing cluster and initiate phase separation. Using both analytical calculations and active particle simulations, we measure these timescales and determine the critical density for phase separation in both two and three dimensions.