Nonequilibrium glassy dynamics of self-propelled hard disks

Identifying the mechanism for superdiffusivity in mouse fibroblast motility

We seek to characterize the motility of mouse fibroblasts on 2D substrates. Utilizing automated tracking techniques, we find that cell trajectories are super-diffusive, where displacements scale faster than t1/2 in all directions. Two mechanisms have been proposed to explain such statistics in other cell types: run and tumble behavior with Lévy-distributed run times, and ensembles of cells with heterogeneous speed and rotational noise. We develop an automated toolkit that directly compares cell trajectories to the predictions of each model and demonstrate that ensemble-averaged quantities such as the mean-squared displacements and velocity autocorrelation functions are equally well-fit by either model. However, neither model correctly captures the short-timescale behavior quantified by the displacement probability distribution or the turning angle distribution. We develop a hybrid model that includes both run and tumble behavior and heterogeneous noise during the runs, which correctly matches the short-timescale behaviors and indicates that the run times are not Lévy distributed. The analysis tools developed here should be broadly useful for distinguishing between mechanisms for superdiffusivity in other cells types and environments.

A random first-order transition theory for an active glass

How does nonequilibrium activity modify the approach to a glass? This is an important question, since many experiments reveal the near-glassy nature of the cell interior, remodeled by activity. However, different simulations of dense assemblies of active particles, parametrized by a self-propulsion force, [Formula: see text], and persistence time, [Formula: see text], appear to make contradictory predictions about the influence of activity on characteristic features of glass, such as fragility. This calls for a broad conceptual framework to understand active glasses; here, we extend the random first-order transition (RFOT) theory to a dense assembly of self-propelled particles. We compute the active contribution to the configurational entropy through an effective model of a single particle in a caging potential. This simple active extension of RFOT provides excellent quantitative fits to existing simulation results. We find that whereas [Formula: see text] always inhibits glassiness, the effect of [Formula: see text] is more subtle and depends on the microscopic details of activity. In doing so, the theory automatically resolves the apparent contradiction between the simulation models. The theory also makes several testable predictions, which we verify by both existing and new simulation data, and should be viewed as a step toward a more rigorous analytical treatment of active glass.

Active Atoms and Interstitials in Two-Dimensional Colloidal Crystals

We study experimentally and numerically the motion of a self-phoretic active particle in two-dimensional loosely packed colloidal crystals at fluid interfaces. Two scenarios emerge depending on the interactions between the active particle and the lattice: the active particle either navigates throughout the crystal as an interstitial or is part of the lattice and behaves as an active atom. Active interstitials undergo a run-and-tumble-like motion, with the passive colloids of the crystal acting as tumbling sites. Instead, active atoms exhibit an intermittent motion, stemming from the interplay between the periodic potential landscape of the passive crystal and the particle's self-propulsion. Our results constitute the first step towards the realization of non-close-packed crystalline phases with internal activity.

Turbulent Dynamics of Epithelial Cell Cultures

We investigate the large length and long time scales collective flows and structural rearrangements within in vitro human bronchial epithelial cell (HBEC) cultures. Activity-driven collective flows result in ensembles of vortices randomly positioned in space. By analyzing a large population of vortices, we show that their area follows an exponential law with a constant mean value and their rotational frequency is size independent, both being characteristic features of the chaotic dynamics of active nematic suspensions. Indeed, we find that HBECs self-organize in nematic domains of several cell lengths. Nematic defects are found at the interface between domains with a total number that remains constant due to the dynamical balance of nucleation and annihilation events. The mean velocity fields in the vicinity of defects are well described by a hydrodynamic theory of extensile active nematics.

From jamming to collective cell migration through a boundary induced transition

Cell monolayers provide an interesting example of active matter, exhibiting a phase transition from flowing to jammed states as they age. Here we report experiments and numerical simulations illustrating how a jammed cellular layer rapidly reverts to a flowing state after a wound. Quantitative comparison between experiments and simulations shows that cells change their self-propulsion and alignment strength so that the system crosses a phase transition line, which we characterize by finite-size scaling in an active particle model. This wound-induced unjamming transition is found to occur generically in epithelial, endothelial and cancer cells.

Flocking transitions in confluent tissues

Collective cell migration in dense tissues underlies important biological processes, such as embryonic development, wound healing and cancer invasion. While many aspects of single cell movements are now well established, the mechanisms leading to displacements of cohesive cell groups are still poorly understood. To elucidate the emergence of collective migration in mechanosensitive cells, we examine a self-propelled Voronoi (SPV) model of confluent tissues with an orientational feedback that aligns a cell's polarization with its local migration velocity. While shape and motility are known to regulate a density-independent liquid-solid transition in tissues, we find that aligning interactions facilitate collective motion and promote solidification, with transitions that can be predicted by extending statistical physics tools such as effective temperature to this far-from-equilibrium system. In addition to accounting for recent experimental observations obtained with epithelial monolayers, our model predicts structural and dynamical signatures of flocking, which may serve as gateway to a more quantitative characterization of collective motility.

No unjamming transition in a Voronoi model of biological tissue

Vertex models are a popular approach to simulating the mechanical and dynamical properties of dense biological tissues, describing the tissue as a network of polygons. Recently a class of two-dimensional vertex models was shown to exhibit a disordered rigidity transition controlled by the preferred cellular geometry, which was subsequently echoed by experimental findings. An attractive variant of these models uses a Voronoi tessellation to describe the cells, which reduces the number of degrees of freedom as compared the original vertex model. The Voronoi model was also endowed with a non-equilibrium model of cellular motility, leading to rich, glassy behavior. This glassy behavior was suggested to be inextricably linked to an underlying jamming transition. We test this conjecture, exploring the low-effective-temperature limit of the 2D Voronoi model by studying cell trajectories from detailed dynamical simulations in combination with rigidity measurements of energy-minimized disordered cell configurations. We find that the zero-temperature limit of this model has no unjamming transition. We show that this absence of an unjamming transition is intimately linked to the marginality of the model, i.e. the fact that the constraints imposed on cell areas and perimeters precisely balance the number of degrees of freedom in the model. Our work suggests that constraint counting arguments are useful to understand rigidity in a broad class of models of dense biological tissues.

Static structure of active Brownian hard disks

We explore the changes in static structure of a two-dimensional system of active Brownian particles (ABP) with hard-disk interactions, using event-driven Brownian dynamics simulations. In particular, the effect of the self-propulsion velocity and the rotational diffusivity on the orientationally-averaged fluid structure factor is discussed. Typically activity increases structural ordering and generates a structure factor peak at zero wave vector which is a precursor of motility-induced phase separation. Our results provide reference data to test future statistical theories for the fluid structure of active Brownian systems. This manuscript was submitted for the special issue of the Journal of Physics: Condensed Matter associated with the Liquid Matter Conference 2017.

Criticality of the zero-temperature jamming transition probed by self-propelled particles

We perform simulations of athermal systems of self-propelled particles (SPPs) interacting via harmonic repulsion in the vicinity of the zero-temperature jamming transition at point J. Every particle is propelled by a constant force f pointing to a randomly assigned and fixed direction. When f is smaller than the yield force f_y, the system is statically jammed. We find that f_y increases with packing fraction and exhibits finite size scaling, implying the criticality of point J. When f > f_y, SPPs flow forever and their velocities satisfy the k-Gamma distribution. Velocity distributions at various packing fractions and f collapse when the particle velocity is scaled by the average velocity v[combining macron], suggesting that v[combining macron] is a reasonable quantity to characterize the response to f. We thus define a response function R(ϕ,f) = v[combining macron](ϕ,f)/f. The function exhibits critical scaling nicely, implying again the criticality of point J. Our analysis and results indicate that systems of SPPs behave analogically to sheared systems, although their driving mechanisms are apparently distinct.

Mode-coupling theory for active Brownian particles

We present a mode-coupling theory (MCT) for the slow dynamics of two-dimensional spherical active Brownian particles (ABPs). The ABPs are characterized by a self-propulsion velocity v_{0} and by their translational and rotational diffusion coefficients D_{t} and D_{r}, respectively. Based on the integration-through-transients formalism, the theory requires as input only the equilibrium static structure factors of the passive system (where v_{0}=0). It predicts a nontrivial idealized-glass-transition diagram in the three-dimensional parameter space of density, self-propulsion velocity, and rotational diffusivity that arise because at high densities, the persistence length of active swimming ℓ_{p}=v_{0}/D_{r} interferes with the interaction length ℓ_{c} set by the caging of particles. While the low-density dynamics of ABPs is characterized by a single Péclet number Pe=v_{0}^{2}/D_{r}D_{t}, close to the glass transition the dynamics is found to depend on Pe and ℓ_{p} separately. At fixed density, increasing the self-propulsion velocity causes structural relaxation to speed up, while decreasing the persistence length slows down the relaxation. The active-MCT glass is a nonergodic state that is qualitatively different from the passive glass. In it, correlations of initial density fluctuations never fully decay, but also an infinite memory of initial orientational fluctuations is retained in the positions.

Universal glass-forming behavior of in vitro and living cytoplasm

Physiological processes in cells are performed efficiently without getting jammed although cytoplasm is highly crowded with various macromolecules. Elucidating the physical machinery is challenging because the interior of a cell is so complex and driven far from equilibrium by metabolic activities. Here, we studied the mechanics of in vitro and living cytoplasm using the particle-tracking and manipulation technique. The molecular crowding effect on cytoplasmic mechanics was selectively studied by preparing simple in vitro models of cytoplasm from which both the metabolism and cytoskeletons were removed. We obtained direct evidence of the cytoplasmic glass transition; a dramatic increase in viscosity upon crowding quantitatively conformed to the super-Arrhenius formula, which is typical for fragile colloidal suspensions close to jamming. Furthermore, the glass-forming behaviors were found to be universally conserved in all the cytoplasm samples that originated from different species and developmental stages; they showed the same tendency for diverging at the macromolecule concentrations relevant for living cells. Notably, such fragile behavior disappeared in metabolically active living cells whose viscosity showed a genuine Arrhenius increase as in typical strong glass formers. Being actively driven by metabolism, the living cytoplasm forms glass that is fundamentally different from that of its non-living counterpart.

Discontinuous fluidization transition in time-correlated assemblies of actively deforming particles

Tracking experiments in dense biological tissues reveal a diversity of sources for local energy injection at the cell scale. The effect of cell motility has been largely studied, but much less is known about the effect of the observed volume fluctuations of individual cells. We consider a microscopic model of "actively deforming" particles where local fluctuations of the particle size constitute a unique source of motion. We demonstrate that collective motion can emerge under the sole influence of such active volume fluctuations. We interpret the onset of diffusive motion as a nonequilibrium first-order phase transition, which arises at a well-defined amplitude of self-deformation. This behavior contrasts with the glassy dynamics produced by self-propulsion, but resembles the mechanical response of soft solids under mechanical deformation. It thus constitutes an example of an active yielding transition.

Jamming transitions in cancer

The traditional picture of tissues, where they are treated as liquids defined by properties such as surface tension or viscosity has been redefined during the last few decades by the more fundamental question: under which conditions do tissues display liquid-like or solid-like behaviour? As a result, basic concepts arising from the treatment of tissues as solid matter, such as cellular jamming and glassy tissues, have shifted into the current focus of biophysical research. Here, we review recent works examining the phase states of tissue with an emphasis on jamming transitions in cancer. When metastasis occurs, cells gain the ability to leave the primary tumour and infiltrate other parts of the body. Recent studies have shown that a linkage between an unjamming transition and tumour progression indeed exists, which could be of importance when designing surgery and treatment approaches for cancer patients.

Nonequilibrium mode-coupling theory for dense active systems of self-propelled particles

The physics of active systems of self-propelled particles, in the regime of a dense liquid state, is an open puzzle of great current interest, both for statistical physics and because such systems appear in many biological contexts. We develop a nonequilibrium mode-coupling theory (MCT) for such systems, where activity is included as a colored noise with the particles having a self-propulsion force f₀ and a persistence time τ_p. Using the extended MCT and a generalized fluctuation-dissipation theorem, we calculate the effective temperature T_eff of the active fluid. The nonequilibrium nature of the systems is manifested through a time-dependent T_eff that approaches a constant in the long-time limit, which depends on the activity parameters f₀ and τ_p. We find, phenomenologically, that this long-time limit is captured by the potential energy of a single, trapped active particle (STAP). Through a scaling analysis close to the MCT glass transition point, we show that τ_α, the α-relaxation time, behaves as τ_α ∼ f₀^-2γ, where γ = 1.74 is the MCT exponent for the passive system. τ_α may increase or decrease as a function of τ_p depending on the type of active force correlations, but the behavior is always governed by the same value of the exponent γ. Comparison with the numerical solution of the nonequilibrium MCT and simulation results give excellent agreement with scaling analysis.

Glassy swirls of active dumbbells

Is an active glass different from a conventional passive glass? To address this, we study the dynamics of a dense binary mixture of soft dumbbells, each subject to an active propulsion force and thermal fluctuations. This dense assembly shows dynamical arrest, first to a translational and then to a rotational glass, as one reduces temperature T or the self-propulsion force f. We monitor the dynamics along an iso-relaxation-time contour in the (T-f) plane. We find dramatic differences both in the fragility and in the nature of dynamical heterogeneity, which characterize the onset of glass formation-the activity-induced glass exhibits large swirls or vortices, whose scale is set by activity, and it appears to diverge as one approaches the glass transition. This large collective swirling movement should have implications for collective cell migration in epithelial layers. We construct continuum hydrodynamic equations for the simulated system, and we show that the observed behavior of this growing dynamic length scale can be understood from these equations.

Aging and rejuvenation of active matter under topological constraints

The coupling of active, self-motile particles to topological constraints can give rise to novel non-equilibrium dynamical patterns that lack any passive counterpart. Here we study the behavior of self-propelled rods confined to a compact spherical manifold by means of Brownian dynamics simulations. We establish the state diagram and find that short active rods at sufficiently high density exhibit a glass transition toward a disordered state characterized by persistent self-spinning motion. By periodically melting and revitrifying the spherical spinning glass, we observe clear signatures of time-dependent aging and rejuvenation physics. We quantify the crucial role of activity in these non-equilibrium processes, and rationalize the aging dynamics in terms of an absorbing-state transition toward a more stable active glassy state. Our results demonstrate both how concepts of passive glass phenomenology can carry over into the realm of active matter, and how topology can enrich the collective spatiotemporal dynamics in inherently non-equilibrium systems.

Using cell deformation and motion to predict forces and collective behavior in morphogenesis

In multi-cellular organisms, morphogenesis translates processes at the cellular scale into tissue deformation at the scale of organs and organisms. To understand how biochemical signaling regulates tissue form and function, we must understand the mechanical forces that shape cells and tissues. Recent progress in developing mechanical models for tissues has led to quantitative predictions for how cell shape changes and polarized cell motility generate forces and collective behavior on the tissue scale. In particular, much insight has been gained by thinking about biological tissues as physical materials composed of cells. Here we review these advances and discuss how they might help shape future experiments in developmental biology.

Mode coupling theory for nonequilibrium glassy dynamics of thermal self-propelled particles

We present a mode coupling theory study for the relaxation and glassy dynamics of a system of strongly interacting self-propelled particles, wherein the self-propulsion force is described by Ornstein-Uhlenbeck colored noise and thermal noises are included. Our starting point is an effective Smoluchowski equation governing the distribution function of particle positions, from which we derive a memory function equation for the time dependence of density fluctuations in nonequilibrium steady states. With the basic assumption of the absence of macroscopic currents and standard mode coupling approximation, we can obtain expressions for the irreducible memory function and other relevant dynamic terms, wherein the nonequilibrium character of the active system is manifested through an averaged diffusion coefficient D[combining macron] and a nontrivial structural function S₂(q) with q being the magnitude of wave vector q. D[combining macron] and S₂(q) enter the frequency term and the vertex term for the memory function, and thus influence both the short time and the long time dynamics of the system. With these equations obtained, we study the glassy dynamics of this thermal self-propelled particle system by investigating the Debye-Waller factor f_q and relaxation time τ_α as functions of the persistence time τ_p of self-propulsion, the single particle effective temperature T_eff as well as the number density ρ. Consequently, we find the critical density ρ_c for given τ_p shifts to larger values with increasing magnitude of propulsion force or effective temperature, in good accordance with previously reported simulation work. In addition, the theory facilitates us to study the critical effective temperature T for fixed ρ as well as its dependence on τ_p. We find that T increases with τ_p and in the limit τ_p → 0, it approaches the value for a simple passive Brownian system as expected. Our theory also well recovers the results for passive systems and can be easily extended to more complex systems such as active-passive mixtures.

Transitions between multiple dynamical states in a confined dense active-particle system

We study the collective motion of a dense suspension of active swimmers in a viscous fluid medium. The swimmers are modeled as soft spheres moving in a highly viscous fluid medium. The magnitude of the propelling thrust exerted by each particle is taken to be a constant and the direction is aligned to its velocity. Depending on the magnitude of the exerted thrust, several nonequilibrium steady states are observed. The transitions between the steady states are characterized using the total dissipation as a function of the magnitude of the thrust. The transitions between the nonequilibrium states are characterized by changes in exponent at low thrust values. At high thrust values, hysteretic transitions between ordered and disordered states are observed.

Cell division and death inhibit glassy behaviour of confluent tissues

We investigate the effects of cell division and apoptosis on collective dynamics in two-dimensional epithelial tissues. Our model includes three key ingredients observed across many epithelia, namely cell-cell adhesion, cell death and a cell division process that depends on the surrounding environment. We show a rich non-equilibrium phase diagram depending on the ratio of cell death to cell division and on the adhesion strength. For large apoptosis rates, cells die out and the tissue disintegrates. As the death rate decreases, however, we show, consecutively, the existence of a gas-like phase, a gel-like phase, and a dense confluent (tissue) phase. Most striking is the observation that the tissue is self-melting through its own internal activity, ruling out the existence of any glassy phase.

Study of dynamic heterogeneity of an active particle system

We have studied spatial and temporal dynamic heterogeneity (DH) in a system of hard-sphere particles, subjected to active forces with constant amplitude and random direction determined by rotational diffusion with correlation time τ. We have used a variety of observables to characterize the DH behavior, including the deviation from standard Stokes-Einstein (SE) relation, a non-Gaussian parameter α_{2}(Δt) for the distribution of particle displacement within a certain time interval Δt, a four-point susceptibility χ_{4}(Δt,ΔL) for the correlation in dynamics between any two points in space separated by distance ΔL within some time window Δt, and a vector spatial-temporal correlation function S_{vec}(R,Δt) for vector displacements within time interval Δt of particle pairs originally separated by R. By mapping the particle motion into a continuous-time random walk with constant jump length, we can obtain the average waiting time 〈t_{x}〉∝D_{s}^{-1} and persistence time 〈t_{p}〉∝η, with D_{s} the self-diffusion coefficient and η the shear viscosity, such that the observable λ=〈t_{p}〉/〈t_{x}〉∝D_{s}η can be calculated as a function of the control parameter τ to show how it deviates from its SE value λ_{0}. Interestingly, we find λ/λ_{0} shows a nonmonotonic behavior for large volume fraction φ_{a}, wherein λ/λ_{0} undergoes a minimum at a certain intermediate value of τ, indicating that both small and large particle activity may lead to strong DH. Such a reentrance phenomenon is further demonstrated in terms of the non-Gaussian parameters α_{2}, four-point susceptibility χ_{4}, and vector spatiotemporal correlation functions S_{vec}, respectively. Detail analysis shows that it is the competition between the dual roles of particle activity, namely, activity-induced higher effective temperature and activity-induced clustering, that leads to such nontrivial nonmonotonic behaviors. In addition, we find that DH may also show a maximum level at an intermediate value of φ_{a} if τ is large enough, implying that a more crowded system may be less heterogeneous than a less crowded one for a system with high particle activity.

Nonlinear Rheology in a Model Biological Tissue

The rheological response of dense active matter is a topic of fundamental importance for many processes in nature such as the mechanics of biological tissues. One prominent way to probe mechanical properties of tissues is to study their response to externally applied forces. Using a particle-based model featuring random apoptosis and environment-dependent division rates, we evidence a crossover from linear flow to a shear-thinning regime with an increasing shear rate. To rationalize this nonlinear flow we derive a theoretical mean-field scenario that accounts for the interplay of mechanical and active noise in local stresses. These noises are, respectively, generated by the elastic response of the cell matrix to cell rearrangements and by the internal activity.

Excitation of vibrational soft modes in disordered systems using active oscillation

We propose a new method to characterize the spatial distribution of particles' vibrations in solids with much lower computational costs compared to the usual normal mode analysis. We excite the specific vibrational mode in a two dimensional athermal jammed system by giving a small amplitude of active oscillation to each particle's size with an identical driving frequency. The response is then obtained as the real time displacements of the particles. We show remarkable correlations between the real time displacements and the eigen vectors obtained from conventional normal mode analysis. More importantly, from these real time displacements, we can measure the participation ratio and spatial polarization of particles' vibrations. From these measurements, we find three distinct frequency regimes which characterize the spatial distribution and correlation of particles' vibrations in jammed amorphous solids. Furthermore, we can possibly apply this method to a much larger system to examine the low frequency behaviour of amorphous solids with a much higher resolution of the frequency space.

Shape and Displacement Fluctuations in Soft Vesicles Filled by Active Particles

We investigate numerically the dynamics of shape and displacement fluctuations of two-dimensional flexible vesicles filled with active particles. At low concentration most of the active particles accumulate at the boundary of the vesicle where positive particle number fluctuations are amplified by trapping, leading to the formation of pinched spots of high density, curvature and pressure. At high concentration the active particles cover the vesicle boundary almost uniformly, resulting in fairly homogeneous pressure and curvature, and nearly circular vesicle shape. The change between polarized and spherical shapes is driven by the number of active particles. The center-of-mass of the vesicle performs a persistent random walk with a long time diffusivity that is strongly enhanced for elongated active particles due to orientational correlations in their direction of propulsive motion. In our model shape-shifting induces directional sensing and the cell spontaneously migrate along the polarization direction.

Structural entropy of glassy systems from graph isomorphism

Configurational entropy plays a central role in thermodynamic scenarios of the glass transition. As a measure of the number of basins in the potential energy landscape, configurational entropy for a glass-forming liquid can be evaluated by explicitly counting distinct inherent structures. In this work, we propose a graph-theory based method to examine local structure and obtain the corresponding entropy of hard-particle systems. Voronoi diagrams of associated clusters are classified using a graph isomorphism algorithm. The statistics of these clusters reveal structural motifs such as icosahedron-like order, and also allow us to calculate the structural entropy SG. We find the structural entropy of an n-particle subsystem grows linearly with n. Thus the structural entropy per particle can be obtained from the slope dSG/dn. Our results are consistent with previous values for configurational entropy obtained via thermodynamic integration. Structural entropies per particle are measured for hard-disk and hard-sphere polydisperse systems, and extrapolated for monodisperse hard disks, all of which are nonzero at the dynamic glass transition.

The nonequilibrium glassy dynamics of self-propelled particles

We study the glassy dynamics taking place in dense assemblies of athermal active particles that are driven solely by a nonequilibrium self-propulsion mechanism. Active forces are modeled as an Ornstein-Uhlenbeck stochastic process, characterized by a persistence time and an effective temperature, and particles interact via a Lennard-Jones potential that yields well-studied glassy behavior in the Brownian limit, which is obtained as the persistence time vanishes. By increasing the persistence time, the system departs more strongly from thermal equilibrium and we provide a comprehensive numerical analysis of the structure and dynamics of the resulting active fluid. Finite persistence times profoundly affect the static structure of the fluid and give rise to nonequilibrium velocity correlations that are absent in thermal systems. Despite these nonequilibrium features, for any value of the persistence time we observe a nonequilibrium glass transition as the effective temperature is decreased. Surprisingly, increasing departure from thermal equilibrium is found to promote (rather than suppress) the glassy dynamics. Overall, our results suggest that with increasing persistence time, microscopic properties of the active fluid change quantitatively, but the general features of the nonequilibrium glassy dynamics observed with decreasing the effective temperature remain qualitatively similar to those of thermal glass-formers.

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.

Active fluidization in dense glassy systems

Dense soft glasses show strong collective caging behavior at sufficiently low temperatures. Using molecular dynamics simulations of a model glass former, we show that the incorporation of activity or self-propulsion, f0, can induce cage breaking and fluidization, resulting in the disappearance of the glassy phase beyond a critical f0. The diffusion coefficient crosses over from being strongly to weakly temperature dependent as f0 is increased. In addition, we demonstrate that activity induces a crossover from a fragile to a strong glass and a tendency of active particles to cluster. Our results are of direct relevance to the collective dynamics of dense active colloidal glasses and to recent experiments on tagged particle diffusion in living cells.

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.

Phase-field model for collective cell migration

We construct a phase-field model for collective cell migration based on a Ginzburg-Landau free-energy formulation. We model adhesion, surface tension, repulsion, coattraction, and polarization, enabling us to follow the cells' morphologies and the effect of their membranes fluctuations on collective motion. We were able to measure the tissue surface tension as a function of the individual cell cortical tension and adhesion and identify a density threshold for cell-sheet formation.

Motility-driven glass and jamming transitions in biological tissues

Cell motion inside dense tissues governs many biological processes, including embryonic development and cancer metastasis, and recent experiments suggest that these tissues exhibit collective glassy behavior. To make quantitative predictions about glass transitions in tissues, we study a self-propelled Voronoi (SPV) model that simultaneously captures polarized cell motility and multi-body cell-cell interactions in a confluent tissue, where there are no gaps between cells. We demonstrate that the model exhibits a jamming transition from a solid-like state to a fluid-like state that is controlled by three parameters: the single-cell motile speed, the persistence time of single-cell tracks, and a target shape index that characterizes the competition between cell-cell adhesion and cortical tension. In contrast to traditional particulate glasses, we are able to identify an experimentally accessible structural order parameter that specifies the entire jamming surface as a function of model parameters. We demonstrate that a continuum Soft Glassy Rheology model precisely captures this transition in the limit of small persistence times, and explain how it fails in the limit of large persistence times. These results provide a framework for understanding the collective solid-to-liquid transitions that have been observed in embryonic development and cancer progression, which may be associated with Epithelial-to-Mesenchymal transition in these tissues.

Emergent self-organization in active materials

Biological systems exhibit large-scale self-organized dynamics and structures which enable organisms to perform the functions of life. The field of active matter strives to develop and understand microscopically driven nonequilibrium materials, with emergent properties comparable to those of living systems. This review will describe two recently developed classes of active matter systems, in which simple building blocks-self-propelled colloidal particles or extensile rod-like particles-self-organize to form macroscopic structures with features not possible in equilibrium systems. We summarize the recent experimental and theoretical progress on each of these systems, and we present simple descriptions of the physics underlying their emergent behaviors.

Mechanics of fire ant aggregations

Fire ants link their bodies to form aggregations; these can adopt a variety of structures, they can drip and spread, or withstand applied loads. Here, by using oscillatory rheology, we show that fire ant aggregations are viscoelastic. We find that, at the lowest ant densities probed and in the linear regime, the elastic and viscous moduli are essentially identical over the spanned frequency range, which highlights the absence of a dominant mode of structural relaxation. As ant density increases, the elastic modulus rises, which we interpret by alluding to ant crowding and subsequent jamming. When deformed beyond the linear regime, the aggregation flows, exhibiting shear-thinning behaviour with a stress load that is comparable to the maximum load the aggregation can withstand before individual ants are torn apart. Our findings illustrate the rich, collective mechanical behaviour that can arise in aggregations of active, interacting building blocks.

Theory for the dynamics of dense systems of athermal self-propelled particles

We present a derivation of a recently proposed theory for the time dependence of density fluctuations in stationary states of strongly interacting, athermal, self-propelled particles. The derivation consists of two steps. First, we start from the equation of motion for the joint distribution of particles' positions and self-propulsions and we integrate out the self-propulsions. In this way we derive an approximate, many-particle evolution equation for the probability distribution of the particles' positions. Second, we use this evolution equation to describe the time dependence of steady-state density correlations. We derive a memory function representation of the density correlation function and then we use a factorization approximation to obtain an approximate expression for the memory function. In the final equation of motion for the density correlation function the nonequilibrium character of the active system manifests itself through the presence of a new steady-state correlation function that quantifies spatial correlations of the velocities of the particles. This correlation function enters into the frequency term, and thus it describes the dependence of the short-time dynamics on the properties of the self-propulsions. More importantly, the correlation function of particles' velocities enters into the vertex of the memory function and through the vertex it modifies the long-time glassy dynamics.

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.

Unjamming and cell shape in the asthmatic airway epithelium

From coffee beans flowing in a chute to cells remodelling in a living tissue, a wide variety of close-packed collective systems-both inert and living-have the potential to jam. The collective can sometimes flow like a fluid or jam and rigidify like a solid. The unjammed-to-jammed transition remains poorly understood, however, and structural properties characterizing these phases remain unknown. Using primary human bronchial epithelial cells, we show that the jamming transition in asthma is linked to cell shape, thus establishing in that system a structural criterion for cell jamming. Surprisingly, the collapse of critical scaling predicts a counter-intuitive relationship between jamming, cell shape and cell-cell adhesive stresses that is borne out by direct experimental observations. Cell shape thus provides a rigorous structural signature for classification and investigation of bronchial epithelial layer jamming in asthma, and potentially in any process in disease or development in which epithelial dynamics play a prominent role.

Glass-like dynamics in confined and congested ant traffic

The collective movement of animal groups often occurs in confined spaces. As animal groups are challenged to move at high density, their mobility dynamics may resemble the flow of densely packed non-living soft materials such as colloids, grains, or polymers. However, unlike inert soft-materials, self-propelled collective living systems often display social interactions whose influence on collective mobility are only now being explored. In this paper, we study the mobility of bi-directional traffic flow in a social insect (the fire ant Solenopsis invicta) as we vary the diameter of confining foraging tunnels. In all tunnel diameters, we observe the emergence of spatially heterogeneous regions of fast and slow traffic that are induced through two phenomena: physical obstruction, arising from the inability of individual ants to interpenetrate, and time-delay resulting from social interaction in which ants stop to briefly antennate. Density correlation functions reveal that the relaxation dynamics of high density traffic fluctuations scale linearly with fluctuation size and are sensitive to tunnel diameter. We separate the roles of physical obstruction and social interactions in traffic flow using cellular automata based simulation. Social interaction between ants is modeled as a dwell time (Tint) over which interacting ants remain stationary in the tunnel. Investigation over a range of densities and Tint reveals that the slowing dynamics of collective motion in social living systems are consistent with dynamics near a fragile glass transition in inert soft-matter systems. In particular, flow is relatively insensitive to density until a critical density is reached. As social interaction affinity is increased (increasing Tint) traffic dynamics change and resemble a strong glass transition. Thus, social interactions play an important role in the mobility of collective living systems at high density. Our experiments and model demonstrate that the concepts of soft-matter physics aid understanding of the mobility of collective living systems, and motivate further inquiry into the dynamics of densely confined social living systems.

Glassy dynamics of athermal self-propelled particles: Computer simulations and a nonequilibrium microscopic theory

We combine computer simulations and analytical theory to investigate the glassy dynamics in dense assemblies of athermal particles evolving under the sole influence of self-propulsion. Our simulations reveal that when the persistence time of the self-propulsion is increased, the local structure becomes more pronounced, whereas the long-time dynamics first accelerates and then slows down. We explain these surprising findings by constructing a nonequilibrium microscopic theory that yields nontrivial predictions for the glassy dynamics. These predictions are in qualitative agreement with the simulations and reveal the importance of steady-state correlations of the local velocities to the nonequilibrium dynamics of dense self-propelled particles.

Understanding the ideal glass transition: lessons from an equilibrium study of hard disks in a channel

We use an exact transfer-matrix approach to compute the equilibrium properties of a system of hard disks of diameter σ confined to a two-dimensional channel of width 1.95σ at constant longitudinal applied force. At this channel width, which is sufficient for next-nearest-neighbor disks to interact, the system is known to have a great many jammed states. Our calculations show that the longitudinal force (pressure) extrapolates to infinity at a well-defined packing fraction ϕ(K) that is less than the maximum possible ϕ(max), the latter corresponding to a buckled crystal. In this quasi-one-dimensional problem there is no question of there being any real divergence of the pressure at ϕ(K). We give arguments that this avoided phase transition is a structural feature, the remnant in our narrow channel system of the hexatic to crystal transition, but that it has the phenomenology of the (avoided) ideal glass transition. We identify a length scale ξ̃(3) as our equivalent of the penetration length for amorphous order: In the channel system, it reaches a maximum value of around 15σ at ϕ(K), which is larger than the penetration lengths that have been reported for three-dimensional systems. It is argued that the α-relaxation time would appear on extrapolation to diverge in a Vogel-Fulcher manner as the packing fraction approaches ϕ(K).

Reentrance in an active glass mixture

Active matter, whose motion is driven, and glasses, whose dynamics are arrested, seem to lie at opposite ends of the spectrum in nonequilibrium systems. In spite of this, both classes of systems exhibit a multitude of stable states that are dynamically isolated from one another. While this defining characteristic is held in common, its origin is different in each case: for active systems, the irreversible driving forces can produce dynamically frozen states, while glassy systems vitrify when they get kinetically trapped on a rugged free energy landscape. In a mixture of active and glassy particles, the interplay between these two tendencies leads to novel phenomenology. We demonstrate this with a spin glass model that we generalize to include an active component. In the absence of a ferromagnetic bias, we find that the spin glass transition temperature depresses with the active fraction, consistent with what has been observed for fully active glassy systems. When a bias does exist, however, a new type of transition becomes possible: the system can be cooled out of the glassy phase. This unusual phenomenon, known as reentrance, has been observed before in a limited number of colloidal and micellar systems, but it has not yet been observed in active glass mixtures. Using low order perturbation theory, we study the origin of this reentrance and, based on the physical picture that results, suggest how our predictions might be measured experimentally.

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.

Active matter transport and jamming on disordered landscapes

We numerically examine the transport of active run-and-tumble particles with steric particle-particle interactions driven with a drift force over random disordered landscapes composed of fixed obstacles. For increasing run lengths, the net particle transport initially increases before reaching a maximum and decreasing at larger run lengths. The transport reduction is associated with the formation of cluster or living crystal states that become locally jammed or clogged by the obstacles. We also find that the system dynamically jams at lower particle densities when the run length is increased. Our results indicate that there is an optimal activity level for transport of run-and-tumble type active matter through quenched disorder and could be important for understanding biological transport in complex environments or for applications of active matter particles in random media.

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

We introduce a kinetic Monte Carlo model for self-propelled hard disks to capture with minimal ingredients the interplay between thermal fluctuations, excluded volume, and self-propulsion in large assemblies of active particles. We analyze in detail the resulting (density, self-propulsion) nonequilibrium phase diagram over a broad range of parameters. We find that purely repulsive hard disks spontaneously aggregate into fractal clusters as self-propulsion is increased and rationalize the evolution of the average cluster size by developing a kinetic model of reversible aggregation. As density is increased, the nonequilibrium clusters percolate to form a ramified structure reminiscent of a physical gel. We show that the addition of a finite amount of noise is needed to trigger a nonequilibrium phase separation, showing that demixing in active Brownian particles results from a delicate balance between noise, interparticle interactions, and self-propulsion. We show that self-propulsion has a profound influence on the dynamics of the active fluid. We find that the diffusion constant has a nonmonotonic behavior as self-propulsion is increased at finite density and that activity produces strong deviations from Fickian diffusion that persist over large time scales and length scales, suggesting that systems of active particles generically behave as dynamically heterogeneous systems.