Search

Activity and autonomous motion are fundamental aspects of many living and engineering systems. Here, the scale of biological agents covers a wide range, from nanomotors, cytoskeleton, and cells, to insects, fish, birds, and people. Inspired by biological active systems, various types of autonomous synthetic nano- and micromachines have been designed, which provide the basis for multifunctional, highly responsive, intelligent active materials. A major challenge for understanding and designing active matter is their inherent non-equilibrium nature due to persistent energy consumption, which invalidates equilibrium concepts such as free energy, detailed balance, and time-reversal symmetry. Furthermore, interactions in ensembles of active agents are often non-additive and non-reciprocal. An important aspect of biological agents is their ability to sense the environment, process this information, and adjust their motion accordingly. It is an important goal for the engineering of micro-robotic systems to achieve similar functionality. With many fundamental properties of motile active matter now reasonably well understood and under control, the ground is prepared for the study of physical aspects and mechanisms of motion in complex environments, of the behavior of systems with new physical features like chirality, of the development of novel micromachines and microbots, of the emergent collective behavior and swarming of intelligent self-propelled particles, and of particular features of microbial systems. The vast complexity of phenomena and mechanisms involved in the self-organization and dynamics of motile active matter poses major challenges, which can only be addressed by a truly interdisciplinary effort involving scientists from biology, chemistry, ecology, engineering, mathematics, and physics.

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. In this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region - a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces., Comment: 6 pages, 4 figures

The feedback between mechanical and chemical signals plays a key role in controlling many biological processes and collective cell behavior. Here we focus on the emergence of spatiotemporal density waves in a one-dimensional "cell train." Combining a minimal theoretical model with observations in an in vitro experimental system of MDCK epithelial cells confined to a linear pattern, we examine the spontaneous oscillations driven by the feedback between myosin activation and mechanical deformations and their effect on the response of the tissue to externally applied deformations. We show that the nature and frequency of spontaneous oscillations is controlled by the size of the cell train, with a transition from size-dependent standing waves to intrinsic spontaneous waves at the natural frequency of the tissue. The response to external boundary perturbations exhibit a resonance at this natural frequency, providing a possible venue for inferring the mechanochemical couplings that control the tissue behavior from rheological experiments., Comment: 15 pages, 10 figures

We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid-liquid interfaces., Comment: 26 pages, 22 figures

Active liquid crystals exert nonequilibrium stresses on their surroundings through constant consumption of energy, giving rise to dynamical steady states not present in equilibrium. The paradigmatic example of an active liquid crystal is a suspension of microtubule bundles powered by kinesin motor proteins, which exhibits self-sustained spatiotemporal chaotic flows. This system has been modelled using continuum theories that couple the microtubule orientation to active flows. Recently the focus has shifted to the interfacial properties of mixtures of active liquid crystals and passive fluids. Active/passive interfaces have been shown to support propagating capillary waves in the absence of inertia and offer a promising route for relating experimental parameters to those of the continuum theory. In this paper we report the derivation of a minimal model that captures the linear dynamics of the interface between an active liquid crystal and a passive fluid. We show that the dynamics of the interface, although powered by active flows throughout the bulk, is qualitatively captured by equations that couple non-reciprocally interface height and nematic director at the interface. This minimal model reproduces the dynamical structure factor evaluated from numerical simulations and the qualitative form of the wave dispersion relation seen in experiments., Comment: 14 pages, 7 figures; Updated title, and Fig.5

Multiphase field models have emerged as an important computational tool for understanding biological tissue while resolving single-cell properties. While they have successfully reproduced many experimentally observed behaviors of living tissue, the theoretical underpinnings have not been fully explored. We show that a two-dimensional version of the model, which is commonly employed to study tissue monolayers, can be derived from a three-dimensional version in the presence of a substrate. We also show how viscous forces, which arise from friction between different cells, can be included in the model. Finally, we numerically simulate a tissue monolayer, and find that intercellular friction tends to solidify the tissue., Comment: 16 pages, 4 figures

Spatiotemporal patterns in multicellular systems are important to understanding tissue dynamics, for instance, during embryonic development and disease. Here, we use a multiphase field model to study numerically the behavior of a near-confluent monolayer of deformable cells with intercellular friction. Varying friction and cell motility drives a solid-liquid transition, and near the transition boundary, we find the emergence of nematic order of cell deformation driven by shear-aligning cellular flows. Intercellular friction endows the monolayer with a finite viscosity, which significantly increases the spatial correlation in the flow and, concomitantly, the extent of nematic order. We also show that hexatic and nematic order are tightly coupled and propose a mechanical-geometric model for the colocalization of +1/2 nematic defects and 5-7 disclination pairs, which are the structural defects in the hexatic phase. Such topological defects coincide with regions of high cell-cell overlap, suggesting that they may mediate cellular extrusion from the monolayer, as found experimentally. Our results delineate a mechanical basis for the recent observation of nematic and hexatic order in multicellular collectives in experiments and simulations and pinpoint a generic pathway to couple topological and physical effects in these systems.

Using a multi-phase field model, we examine how particle deformability, which is a proxy for cell stiffness, affects motility induced phase separation (MIPS). We show that purely repulsive deformable, i.e., squishy, cells phase separate more effectively than their rigid counterparts. This can be understood as due to the fact that deformability increases the effective duration of collisions. In addition, the dense regions become increasingly disordered as deformability increases. Our results contextualize the applicability of MIPS to biological systems and have implications for how cells in biological systems may self-organize, Comment: 7 pages, 6 figures

In recent years, nonreciprocally coupled systems have received growing attention. Previous work has shown that the interplay of nonreciprocal coupling and Goldstone modes can drive the emergence of temporal order such as traveling waves. We show that these phenomena are generically found in a broad class of pattern-forming systems, including mass-conserving reaction--diffusion systems and viscoelastic active gels. All these systems share a characteristic dispersion relation that acquires a non-zero imaginary part at the edge of the band of unstable modes and exhibit a regime of propagating structures (traveling wave bands or droplets). We show that models for these systems can be mapped to a common "normal form" that can be seen as a spatially extended generalization of the FitzHugh--Nagumo model, providing a unifying dynamical-systems perspective. We show that the minimal nonreciprocal Cahn--Hilliard (NRCH) equations exhibit a surprisingly rich set of behaviors, including interrupted coarsening of traveling waves without selection of a preferred wavelength and transversal undulations of wave fronts in two dimensions. We show that the emergence of traveling waves and their speed are precisely predicted from the local dispersion relation at interfaces far away from the homogeneous steady state. Our work thus generalizes previously studied nonreciprocal phase transitions and shows that interfaces are the relevant collective excitations governing the rich dynamical patterns of conserved fields., Comment: 21 pages, 12 figures (including appendices), videos: https://youtube.com/playlist?list=PLdejNN8-Lt6HyU9G9eB8Pp0-dyJvg9AJs; v2: minor changes in introduction and discussion; v3 major revision; v4 minor changes, final accepted version

We use numerical simulations and linear stability analysis to study the dynamics of an active liquid crystal film on a substrate in the regime where the passive system would be isotropic. Extensile activity builds up local orientational order and destabilizes the quiescent isotropic state above a critical activity value, eventually resulting in spatiotemporal chaotic dynamics akin to the one observed ubiquitously in the nematic state. Here we show that tuning substrate friction yields a variety of emergent structures at intermediate activity, including lattices of flow vortices with associated regular arrangements of topological defects and a new state where flow vortices trap pairs of $+1/2$ defect that chase each other tail. These chiral units spontaneously pick the sense of rotation and organize in a hexagonal lattice, surrounded by a diffuse flow of opposite rotation to maintain zero net vorticity. The length scale of these emergent structures is set by the screening length $l_\eta=\sqrt{\eta/\Gamma}$ of the flow, controlled by the shear viscosity $\eta$ and the substrate friction $\Gamma$, and can be captured by simple mode selection of the vortical flows. We demonstrate that the emergence of coherent structures can be interpreted as a phase separation of vorticity, where friction plays a role akin to that of birth/death processes in breaking conservation of the phase separating species and selecting a characteristic scale for the patterns. Our work shows that friction provides an experimentally accessible tuning parameter for designing controlled active flows.

Using a mean field approach and simulation, we study the non-linear mechanical response of the vertex model (VM) of biological tissue under compression and dilation. The VM is known to exhibit a transition between rigid and fluid-like, or floppy, states driven by geometric incompatibility. Target perimeter and area set a target shape which may not be geometrically achievable, thereby engendering frustration. Previously, an asymmetry in the linear elastic response was identified at the rigidity transition between compression and dilation. Here we show and characterize how the asymmetry extends away from the transition point for finite strains. Under finite compression, an initially solid VM can totally relax perimeter tension, and thereby have reduced bulk and shear modulus. Conversely, an initially floppy VM under dilation can rigidify and have a higher bulk and shear modulus. These observations imply that re-scaling of cell area shifts the transition between rigid and floppy states. Based on this insight, we calculate the re-scaling of cell area engendered by intrinsic curvature and write a prediction for the rigidity transition in the presence of curvature. The shift of the rigidity transition in the presence of curvature for the VM provides a new metric for predicting tissue rigidity from image data for curved tissues in a manner analogous to the flat case., Comment: 10 pages, 3 figures

Orientational order, encoded in anisotropic fields, plays an important role during the development of an organism. A striking example of this is the freshwater polyp Hydra, where topological defects in the muscle fiber orientation have been shown to localize to key features of the body plan. This body plan is organized by morphogen concentration gradients, raising the question how muscle fiber orientation, morphogen gradients and body shape interact. Here, we introduce a minimal model that couples nematic orientational order to the gradient of a morphogen field. We show that on a planar surface alignment to a radial concentration gradient can induce unbinding of topological defects, as observed during budding and tentacle formation in Hydra, and stabilize aster/vortex-like defects, as observed at a Hydra's mouth. On curved surfaces mimicking the morphologies of Hydra in various stages of development -- from spheroid to adult -- our model reproduces the experimentally observed reorganization of orientational order. Our results suggest how gradient alignment and curvature effects may work together to control orientational order during development and lays the foundations for future modeling efforts that will include the tissue mechanics that drive shape deformations., Comment: v2: Published version with minor revisions. v1: 12 pages, 7 figures (4 main, 3 suppl.). Movies: https://youtube.com/playlist?list=PLdejNN8-Lt6HJ52YfhyegpI0PGoRq8pSb

Topological defects play a central role in the physics of many materials, including magnets, superconductors and liquid crystals. In active fluids, defects become autonomous particles that spontaneously propel from internal active stresses and drive chaotic flows stirring the fluid. The intimate connection between defect textures and active flow suggests that properties of active materials can be engineered by controlling defects, but design principles for their spatiotemporal control remain elusive. Here we propose a symmetry-based additive strategy for using elementary activity patterns, as active topological tweezers, to create, move and braid such defects. By combining theory and simulations, we demonstrate how, at the collective level, spatial activity gradients act like electric fields which, when strong enough, induce an inverted topological polarization of defects, akin to a negative susceptibility dielectric. We harness this feature in a dynamic setting to collectively pattern and transport interacting active defects. Our work establishes an additive framework to sculpt flows and manipulate active defects in both space and time, paving the way to design programmable active and living materials for transport, memory and logic., Comment: 11 pages (including Methods), 5 figures. Changed title and format, final version

The vertex model of epithelia describes the apical surface of a tissue as a tiling of polygonal cells, with a mechanical energy governed by deviations in cell shape from preferred, or target, area, $A_0$, and perimeter, $P_0$. The model exhibits a rigidity transition driven by geometric incompatibility as tuned by the target shape index, $p_0 = P_0 / \sqrt{A_0}$. For $p_0 > p_*(6) = \sqrt{8 \sqrt{3}} \approx 3.72$, with $p_*(6)$ the perimeter of a regular hexagon of unit area, a cell can simultaneously attain both the preferred area and preferred perimeter. As a result, the tissue is in a mechanically soft compatible state, with zero shear and Young's moduli. For $p_0 < p_*(6)$, it is geometrically impossible for any cell to realize the preferred area and perimeter simultaneously, and the tissue is in an incompatible rigid solid state. Using a mean-field approach, we present a complete analytical calculation of the linear elastic moduli of an ordered vertex model. We analyze a relaxation step that includes non-affine deformations, leading to a softer response than previously reported. The origin of the vanishing shear and Young's moduli in the compatible state is the presence of zero-energy deformations of cell shape. The bulk modulus exhibits a jump discontinuity at the transition and can be lower in the rigid state than in the fluid-like state. The Poisson's ratio can become negative which lowers the bulk and Young's moduli. Our work provides a unified treatment of linear elasticity for the vertex model and demonstrates that this linear response is protocol-dependent.

The ability of cells to reorganize in response to external stimuli is important in areas ranging from morphogenesis to tissue engineering. Elongated cells can co-align due to steric effects, forming states with local order. We show that molecular-scale substrate anisotropy can direct cell organization, resulting in the emergence of nematic order on tissue scales. To quantitatively examine the disorder-order transition, we developed a high-throughput imaging platform to analyze velocity and orientational correlations for several thousand cells over days. The establishment of global, seemingly long-ranged order is facilitated by enhanced cell division along the substrate's nematic axis, and associated extensile stresses that restructure the cells' actomyosin networks. Our work, which connects to a class of systems known as active dry nematics, provides a new understanding of the dynamics of cellular remodeling and organization in weakly interacting cell collectives. This enables data-driven discovery of cell-cell interactions and points to strategies for tissue engineering., Comment: 29 pages, 7 figures

The rheology of biological tissue is key to processes such as embryo development, wound healing and cancer metastasis. Vertex models of confluent tissue monolayers have uncovered a spontaneous liquid-solid transition tuned by cell shape; and a shear-induced solidification transition of an initially liquid-like tissue. Alongside this jamming/unjamming behaviour, biological tissue also displays an inherent viscoelasticity, with a slow time and rate dependent mechanics. With this motivation, we combine simulations and continuum theory to examine the rheology of the vertex model in nonlinear shear across a full range of shear rates from quastistatic to fast, elucidating its nonlinear stress-strain curves after the inception of shear of finite rate, and its steady state flow curves of stress as a function of strain rate. We formulate a rheological constitutive model that couples cell shape to flow and captures both the tissue solid-liquid transition and its rich linear and nonlinear rheology., Comment: 5 pages, 3 figures; plus supplemental material

We study the dynamics of topological defects in active nematic films with spatially-varying activity and consider two setups: i) a constant activity gradient, and ii) a sharp jump in activity. A constant gradient of extensile (contractile) activity endows the comet-like $+1/2$ defect with a finite vorticity that drives the defect to align its nose in the direction of decreasing (increasing) gradient. A constant gradient does not, however, affect the known self-propulsion of the $+1/2$ defect and has no effect on the $-1/2$ that remains a non-motile particle. A sharp jump in activity acts like a wall that traps the defects, affecting the translational and rotational motion of both charges. The $+1/2$ defect slows down as it approaches the interface and the net vorticity tends to reorient the defect polarization so that it becomes perpendicular to the interface. The $-1/2$ defect acquires a self-propulsion towards the activity interface, while the vorticity-induced active torque tends to align the defect to a preferred orientation. This effective attraction of the negative defects to the wall is consistent with the observation of an accumulation of negative topological charge at both active/passive interfaces and physical boundaries., Comment: 14 pages, 8 figures

Demixing of binary liquids is a ubiquitous transition, which is explained using a well-established thermodynamic formalism that requires equality of intensive thermodynamics parameters across the phase boundaries. Demixing transitions also occur when binary fluid mixtures are driven away from equilibrium, for example, by external shear flow. Predicting demixing transition under non-equilibrium non-potential conditions remains, however, a challenge. We drive liquid-liquid phase separation (LLPS) of attractive DNA nanostar molecules away from equilibrium using an internally driven microtubule-based active fluid. Activity lowers the critical temperature and narrows the coexistence concentrations, but only when there are mechanical bonds between the liquid droplets and the reconfiguring active fluid. Similar behaviors are observed in numerical simulations, suggesting that activity suppression of liquid-liquid phase separation is a generic feature of active LLPS. Our work describes a platform for building soft active materials with feedback control while also providing insight into cell biology, where phase separation emerged as a ubiquitous self-organizational principle., Comment: 9 pages, 6 figures, supplementary information

Controlling interfaces of phase separating fluid mixtures is key to creating diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we study how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and non-inertial propagating active waves. At high activities, stresses disrupt interface continuity and drive droplet generation, producing an emulsion-like active state comprised of finite-sized droplets. When in contact with a solid boundary, active interfaces exhibit non-equilibrium wetting transitions, wherein the fluid climbs the wall against gravity. These results demonstrate the promise of mechanically driven interfaces for creating a new class of soft active matter., Comment: 62 pages, 5 figures, 13 supplementary figures

We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint.Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, whilst they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smectic-like pattern, which induce an associated global polar ordering of $+1/2$ nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the $+1/2$ defects, consisting of anti-aligned pairs of neighboring defects, as established in previous studies., Comment: 10 pages, 4 figures

We use a continuum model to examine the effect of activity on a phase separating mixture of an extensile active nematic and a passive fluid. We highlight the distinct role of previously considered interfacial active stresses and bulk active stresses that couple to liquid crystalline degrees of freedom. Interfacial active stresses can arrest phase separation, as previously demonstrated. Bulk extensile active stresses can additionally strongly suppress phase separation by sustained self-stirring of the fluid, substantially reducing the size of the coexistence region in the temperature/concentration plane relative to that of the passive system. The phase separated state is a dynamical emulsion-like steady state of continuously splitting and merging droplets, as suggested by recent experiments. Using scaling analysis and simulations, we identify various regimes for the dependence of droplet size on activity. These results can provide a criterion for identifying the mechanisms responsible for arresting phase separation in experiments.

Boundary conditions dictate how fluids, including liquid crystals, flow when pumped through a channel. Can boundary conditions also be used to control internally driven active fluids that generate flows spontaneously? By using numerical simulations and stability analysis we explore how surface anchoring of active agents at the boundaries and substrate drag can be used to rectify coherent flow of an active polar fluid in a 2D channel. Upon increasing activity, a succession of dynamical states is obtained, from laminar flow to vortex arrays to eventual turbulence, that are controlled by the interplay between the hydrodynamic screening length and the extrapolation length quantifying the anchoring strength of the orientational order parameter. We highlight the key role of symmetry in both flow and order and show that coherent laminar flow with net throughput is only possible for weak anchoring and intermediate activity. Our work demonstrates the possibility of controlling the nature and properties of active flows in a channel simply by patterning the confining boundaries., Comment: 12 pages, 9 figures

The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multi-phase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solid-like state and a fluid-like state tuned by cell motility and deformability - the ratio of the energetic costs of steric cell-cell repulsion and cell surface tension. When perturbed, solid tissues exhibit yield-stress behavior, with a threshold force for the onset of motion of a probe particle that vanishes upon approaching the solid-to-liquid transition. This onset of motion is qualitatively different in the low and high deformability regimes. At high deformability, the tissue is amorphous when solid, it responds compliantly to deformations, and the probe transition to motion is smooth. At low deformability, the monolayer is more ordered translationally and stiffer, and the onset of motion appears discontinuous. Our results suggest that cellular or nanoparticle transport in different types of tissues can be fundamentally different, and point to ways in which it can be controlled., Comment: Main text (6 pages, 4 figures) plus SM (8 pages, 9 figures)

We study the active flow around isolated defects and the self-propulsion velocity of $+1/2$ defects in an active nematic film with both viscous dissipation (with viscosity $\eta$) and frictional damping $\Gamma$ with a substrate. The interplay between these two dissipation mechanisms is controlled by the hydrodynamic dissipation length $\ell_d=\sqrt{\eta/\Gamma}$ that screens the flows. For an isolated defect, in the absence of screening from other defects, the size of the vortical flows around the defect is controlled by the system size $R$. In the presence of friction that leads to a finite value of $\ell_d$, the vorticity field decays to zero on the lengthscales larger than $\ell_d$. We show that the self-propulsion velocity of $+1/2$ defects grows with $R$ in small systems where $R<\ell_d$, while in the infinite system limit or when $R\gg \ell_d$, it approaches a constant value determined by $\ell_d$., Comment: 21 pages, 8 figures

Vertex Models, as used to describe cellular tissue, have an energy controlled by deviations of each cell area and perimeter from target values. The constrained nonlinear relation between area and perimeter leads to new mechanical response. Here we provide a mean-field treatment of a highly simplified model: a uniform network of regular polygons with no topological rearrangements. Since all polygons deform in the same way, we only need to analyze the ground states and the response to deformations of a single polygon (cell). The model exhibits the known transition between a fluid/compatible state, where the cell can accommodate both target area and perimeter, and a rigid/incompatible state. %The rigid solid-like state has a single gapped ground state. We calculate and measure the mechanical resistance to various deformation protocols and discover that at the onset of rigidity, where a single zero-energy ground-state exists, %We show that in the incompatible state, where a single frustrated ground-state exists, linear elasticity fails to describe the mechanical response to even infinitesimal deformations. In particular we identify a breakdown of reciprocity expressed via different moduli for compressive and tensile loads, implying non-analyticity of the energy functional. We give a pictorial representation in configuration space that reveals that the complex elastic response of the Vertex Model arises from the presence of two distinct sets of reference states (associated with target area and target perimeter).

Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confluent tissues under simple shear deformation. We show that an initially undeformed fluidlike tissue acquires finite rigidity above a critical applied strain. This is akin to the shear-driven rigidity observed in other soft matter systems. Interestingly, shear-driven rigidity can be understood by a critical scaling analysis in the vicinity of the second order critical point that governs the liquid-solid transition of the undeformed system. We further show that a solidlike tissue responds linearly only to small strains and but then switches to a nonlinear response at larger stains, with substantial stiffening. Finally, we propose a mean-field formulation for cells under shear that offers a simple physical explanation of shear-driven rigidity and nonlinear response in a tissue.

Experimental evidence shows that there is a feedback between cell shape and cell motion. How this feedback impacts the collective behavior of dense cell monolayers remains an open question. We investigate the effect of a feedback that tends to align the cell crawling direction with cell elongation in a biological tissue model. We find that the alignment interaction promotes nematic patterns in the fluid phase that eventually undergo a non-equilibrium phase transition into a quasi-hexagonal solid. Meanwhile, highly asymmetric cells do not undergo the liquid-to-solid transition for any value of the alignment coupling. In this regime, the dynamics of cell centers and shape fluctuation show features typical of glassy systems.

This article summarizes some of the open questions in the field of active matter that have emerged during Active20, a nine-week program held at the Kavli Institute for Theoretical Physics (KITP) in Spring 2020. The article does not provide a review of the field, but rather a personal view of the authors, informed by contributions of all participants, on new directions in active matter research. The topics highlighted include: the ubiquitous occurrence of spontaneous flows and active turbulence and the theoretical and experimental challenges associated with controlling and harnessing such flows; the role of motile topological defects in ordered states of active matter and their possible biological relevance; the emergence of non-reciprocal effective interactions and the role of chirality in active systems and their intriguing connections to non-Hermitian quantum mechanics; the progress towards a formulation of the thermodynamics of active systems thanks to the feedback between theory and experiments; the impact of the active matter framework on our understanding of the emergent mechanics of biological tissue. These seemingly diverse phenomena all stem from the defining property of active matter - assemblies of self-driven entities that individually break time-reversal symmetry and collectively organize in a rich variety of nonequilibrium states., Comment: 23 pages, 10 figures

We adapt the Halperin-Mazenko formalism to analyze two-dimensional active nematics coupled to a generic fluid flow. The governing hydrodynamic equations lead to evolution laws for nematic topological defects and their corresponding density fields. We find that $\pm 1/2$ defects are propelled by the local fluid flow and by the nematic orientation coupled with the flow shear rate. In the overdamped and compressible limit, we recover the previously obtained active self-propulsion of the +1/2 defects. Non-local hydrodynamic effects are primarily significant for incompressible flows, for which it is not possible to eliminate the fluid velocity in favor of the local defect polarization alone. For the case of two defects with opposite charge, the non-local hydrodynamic interaction is mediated by non-reciprocal pressure-gradient forces. Finally, we derive continuum equations for a defect gas coupled to an arbitrary (compressible or incompressible) fluid flow.

Using the Poisson bracket method, we derive continuum equations for a fluid of deformable particles in two dimensions. Particle shape is quantified in terms of two continuum fields: an anisotropy density field that captures the deformations of individual particles from regular shapes and a shape tensor density field that quantifies both particle elongation and nematic alignment of elongated shapes. We explicitly consider the example of a dense biological tissue as described by the Vertex model energy, where cell shape has been proposed as a structural order parameter for a liquid-solid transition. The hydrodynamic model of biological tissue proposed here captures the coupling of cell shape to flow, and provides a starting point for modeling the rheology of dense tissue., Comment: 8 pages, 1 figure

Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix., Comment: 6 pages, 1 figure. v2: minor edits and clarifying comments

Active matter encompasses different nonequilibrium systems in which individual constituents convert energy into non-conservative forces or motion at the microscale. This review provides an elementary introduction to the role of topology in active matter through experimentally relevant examples. Here, the focus lies on topological defects and topologically protected edge modes with an emphasis on the distinctive properties they acquire in active media. These paradigmatic examples represent two physically distinct classes of phenomena whose robustness can be traced to a common mathematical origin: the presence of topological invariants. These invariants are typically integer numbers that cannot be changed by continuous deformations of the relevant order parameters or physical parameters of the underlying medium. We first explain the mechanisms whereby topological defects self propel and proliferate in active nematics, leading to collective states which can be manipulated by geometry and patterning. Possible implications for active microfluidics and biological tissues are presented. We then illustrate how the propagation of waves in active fluids and solids is affected by the presence of topological invariants characterizing their dispersion relations. We discuss the relevance of these ideas for the design of robotic metamaterials and the properties of active granular and colloidal systems. Open theoretical and experimental challenges are presented as future research prospects., Comment: 23 pages, 3 figures, 4 boxes

Active matter consists of units that generate mechanical work by consuming energy. Examples include living systems, such as assemblies of bacteria and biological tissues, biopolymers driven by molecular motors, and suspensions of synthetic self-propelled particles. A central question in the field is to understand and control the self-organization of active assemblies in space and time. Most active systems exhibit either spatial order mediated by interactions that coordinate the spatial structure and the motion of active agents or the temporal synchronization of individual oscillatory dynamics. The simultaneous control of spatial and temporal organization is more challenging and generally requires complex interactions, such as reaction-diffusion hierarchies or genetically engineered cellular circuits. Here, we report a novel and simple means to simultaneously control the spatial and temporal self-organization of bacterial active matter. By confining an active bacterial suspension and manipulating a single macroscopic parameter, namely the viscoelasticity of the suspending fluid, we have found that the bacterial fluid first self-organizes in space into a millimeter-scale rotating vortex; then displays temporal organization as the giant vortex switches its global chirality periodically with tunable frequency, reminiscent of a torsional pendulum - a self-driven one. Combining experiments with an active matter model, we explain this striking behavior in terms of the interplay between active forcing and viscoelastic stress relaxation. Our findings advance the understanding of bacterial behavior in complex fluids, and demonstrate experimentally for the first time that rheological properties can be harnessed to control active matter flows. Coupled with actuation, our tunable self-oscillating bacterial vortex may be used as a "clock" for locomotion of soft robots and microfluidic pumping., Comment: 29 pages with 4 figure including Methods and Extended Data (with 10 figures); SI available upon request (accepted version)

Recent experiments and numerical studies have drawn attention to the dynamics of active nematics. Two-dimensional active nematics flow spontaneously and exhibit spatiotemporal chaotic flows with proliferation of topological defects in the nematic texture. It has been proposed that the dynamics of active nematics can be understood in terms of the dynamics of interacting defects, propelled by active stress. Previous work has derived effective equations of motion for individual defects as quasi-particles moving in the mean field generated by other defects, but an effective theory governing multi-defect dynamics has remained out of reach. In this paper, we examine the dynamics of 2D active nematics in the limit of strong order and overdamped compressible flow. The activity-induced defect dynamics is formulated as a perturbation of the manifold of quasi-static nematic textures explicitly parameterized by defect positions. This makes it possible to derive a set of coupled ordinary differential equations governing defect (and therefore texture) dynamics. Interestingly, because of the non-orthogonality of textures associated with individual defects, their motion is coupled through a position dependent ``collective mobility" matrix. In addition to the familiar active self-propulsion of the $+1/2$ defect, we obtain new collective effects of activity that can be interpreted in terms of non-central and non-reciprocal interactions between defects., Comment: 14 pages, 3 figures

Cellularized tissue and polymer networks can both transition from floppy to rigid as a function of their control parameters, and, yet, the two systems often mechanically interact, which may affect their respective rigidities. To study this interaction, we consider a vertex model with interfacial tension (a spheroid) embedded in a spring network in two dimensions. We identify two regimes with different global spheroid shapes, governed by the pressure resulting from competition between interfacial tension and tension in the network. In the first regime, the tissue remains compact, while in the second, a cavitation-like instability leads to the emergence of gaps at the tissue-network interface. Intriguingly, compression of the tissue promotes fluidization, while tension promotes cellular alignment and rigidification, with the mechanisms driving rigidification differing on either side of the instability., Comment: Main text: 5 pages, 3 figures. Supplementary Material: 13 pages, 14 figures

We examine a non-reciprocally coupled dynamical model of a mixture of two diffusing species. We demonstrate that nonreciprocity, which is encoded in the model via antagonistic cross diffusivities, provides a generic mechanism for the emergence of traveling patterns in purely diffusive systems with conservative dynamics. In the absence of non-reciprocity, the binary fluid mixture undergoes a phase transition from a homogeneous mixed state to a demixed state with spatially separated regions rich in one of the two components. Above a critical value of the parameter tuning non-reciprocity, the static demixed pattern acquires a finite velocity, resulting in a state that breaks both spatial and time translational symmetry, as well as the reflection parity of the static pattern. We elucidate the generic nature of the transition to traveling patterns using a minimal model that can be studied analytically. Our work has direct relevance to nonequilibrium assembly in mixtures of chemically interacting colloids that are known to exhibit non-reciprocal effective interactions, as well as to mixtures of active and passive agents where traveling states of the type predicted here have been observed in simulations. It also provides insight on transitions to traveling and oscillatory states seen in a broad range of nonreciprocal systems with non-conservative dynamics, from reaction-diffusion and prey-predators models to multispecies mixtures of microorganisms with antagonistic interactions., Comment: 8 pages, 3 figures

In equilibrium, disorder conspires with topological defects to redefine the ordered states of matter in systems as diverse as crystals, superconductors and liquid crystals. Far from equilibrium, however, the consequences of quenched disorder on active condensed matter remain virtually uncharted. Here, we reveal a state of strongly disordered active matter with no counterparts in equilibrium: a dynamical vortex glass. Combining high-content microfluidic experiments and theory, we show how colloidal flocks collectively cruise through disorder without relaxing the topological singularities of their flows. The resulting state is highly dynamical but the flow patterns, shaped by a finite density of frozen vortices, are stationary and exponentially degenerated. Quenched isotropic disorder acts as a random gauge field turning active liquids into dynamical vortex glasses. We argue that this robust mechanism should shape the collective dynamics of a broad class of disordered active matter, from synthetic active nematics to collections of living cells exploring heterogeneous media., Comment: 9 pages, 3 figures (corrected a Grant number)

Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways - the presence of an apical-basal polarity, spatial variability of cellular thickness, and their nonequilibrium active nature. Here, we develop a minimal continuum model of a planar epithelial tissue as an active elastic material that incorporates all these features. We start from a full three-dimensional (3D)description of the tissue and derive an effective 2D model that captures, through the curvature of the apical surface, both the apical-basal asymmetry and the spatial geometry of the tissue. Crucially, variations of active stresses across the apical-basal axis lead to active torques that can drive curvature transitions. By identifying four distinct sources of activity, we find that bulk active stresses arising from actomyosin contractility and growth compete with boundary active tensions due to localized actomyosin cables and lamellipodial activity to generate the various states spanning the morphospace of a planar epithelium. Our treatment hence unifies 3D shape deformations through the coupled mechanics of apical curvature change and in-plane expansion/contraction of substrate-adhered tissues. Finally, we discuss the implications of our results for some biologically relevant processes such as tissue folding at the onset of lumen formation., Comment: 12 pages, 8 figures (including Appendix), published version including more calculations without changing results

We examine the interplay of motility and information exchange in a model of run-and-tumble active particles where the particle's motility is encoded as a bit of information that can be exchanged upon contact according to the rules of AND and OR logic gates in a circuit. Motile AND particles become non-motile upon contact with a non-motile particle. Conversely, motile OR particles remain motile upon collision with their non-motile counterparts. AND particles that have become non-motile additionally "reawaken", i.e., recover their motility, at a fixed rate $\mu$, as in the SIS (Susceptible, Infected, Susceptible) model of epidemic spreading, where an infected agent can become healthy again, but keeps no memory of the recent infection, hence it is susceptible to a renewed infection. For $\mu=0$, both AND and OR particles relax irreversibly to absorbing states of all non-motile or all motile particles, respectively. The relaxation kinetics is, however, faster for OR particles that remain active throughout the process. At finite $\mu$, the AND dynamics is controlled by the interplay between reawakening and collision rates. The system evolves to a state of all motile particles (an absorbing state in the language of absorbing phase transitions) for $\mu>\mu_c$ and to a mixed state with coexisting motile and non-motile particles (an active state in the language of absorbing phase transitions) for $\mu<\mu_c$. The final state exhibits a rich structure controlled by motility-induced aggregation. Our work can be relevant to biochemical signaling in motile bacteria, the spreading of epidemics and of social consensus, as well as light-controlled organization of active colloids.

Experiments and theory have shown that cell monolayers and epithelial tissues exhibit solid-liquid and glass-liquid transitions. These transitions are biologically relevant to our understanding of embryonic development, wound healing, and cancer. Current models typically consider purely two-dimensional monolayers with no overlaps between neighboring cells. In reality, overlaps are important, and they may be precursors of cell extrusion -- a key biophysical process to maintain homeostasis in epithelial tissues. Here, we use a multi-phase field model to study the solid-liquid transition in a confluent monolayer of deformable cells which can overlap. When cells overlap rather than deform, we find that the melting transition changes from continuous to discontinuous, and that there is an intermittent regime close to the transition, where solid and liquid states alternate over time. By studying the dynamics of $5$- and $7$-fold disclinations in the hexagonal lattice formed by the cell centers, we observe that these correlate with spatial fluctuations in the cellular overlap, and that cell extrusion tends to initiate near $5$-fold disclinations., Comment: Figures have been updated and the text expanded. Main text (6 pages, 4 figures) plus SM (12 pages, 16 figures). Supplementary movies can be found at https://drive.google.com/open?id=1yfMf0aFqFlgjsGGjE0Zr3zo_eiIcREB9

Topological defects play a prominent role in the physics of two-dimensional materials. When driven out of equilibrium in active nematics, disclinations can acquire spontaneous self-propulsion and drive self-sustained flows upon proliferation. Here we construct a general hydrodynamic theory for a two-dimensional active nematic interrupted by a large number of such defects. Our equations describe the flows and spatio-temporal defect chaos characterizing active turbulence, even close to the defect unbinding transition. At high activity, nonequilibrium torques combined with many-body screening cause the active disclinations to spontaneously break rotational symmetry forming a collectively moving defect ordered polar liquid. By recognizing defects as the relevant quasiparticle excitations, we construct a comprehensive phase diagram for two-dimensional active nematics. Using our hydrodynamic approach, we additionally show that activity gradients can act like "electric fields", driving the sorting of topological charge. This demonstrates the versatility of our continuum model and its relevance for quantifying the use of spatially inhomogeneous activity for controlling active flows and for the fabrication of active devices with targeted transport capabilities., Comment: 18 pages, 7 figures, additional explanation provided with results unchanged

The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been proposed. Here, we provide a justification of such models for the case of dense suspensions in two dimensions (2d). We subsequently carry out an in-depth numerical study of the properties of one-fluid models as a function of the active driving in view of possible transition scenarios from active turbulence to large-scale pattern, referred to as condensate, formation induced by the classical inverse energy cascade in Newtonian 2d turbulence. Using a one-fluid model it was recently shown (Linkmann et al., Phys. Rev. Lett. (in press)) that two-dimensional active suspensions support two non-equilibrium steady states, one with a condensate and one without, which are separated by a subcritical transition. Here, we report further details on this transition such as hysteresis and discuss a low-dimensional model that describes the main features of the transition through nonlocal-in-scale coupling between the small-scale driving and the condensate.

Recent work on particle-based models of tissues has suggested that any finite rate of cell division and cell death is sufficient to fluidize an epithelial tissue. At the same time, experimental evidence has indicated the existence of glassy dynamics in some epithelial layers despite continued cell cycling. To address this discrepancy, we quantify the role of cell birth and death on glassy states in confluent tissues using simulations of an active vertex model that includes cell motility, cell division, and cell death. Our simulation data is consistent with a simple ansatz in which the rate of cell-life cycling and the rate of relaxation of the tissue in the absence of cell cycling contribute independently and additively to the overall rate of cell motion. Specifically, we find that a glass-like regime with caging behavior indicated by subdiffusive cell displacements can be achieved in systems with sufficiently low rates of cell cycling., Comment: Main text (10 pages, 5 figures) and Appendices (7 pages, 7 figures)

Surface tension governed by differential adhesion can drive fluid particle mixtures to sort into separate regions, i.e., demix. Does the same phenomenon occur in confluent biological tissues? We begin to answer this question for epithelial monolayers with a combination of theory via a vertex model and experiments on keratinocyte monolayers. Vertex models are distinct from particle models in that the interactions between the cells are shape-based, as opposed to distance-dependent. We investigate whether a disparity in cell shape or size alone is sufficient to drive demixing in bidisperse vertex model fluid mixtures. Surprisingly, we observe that both types of bidisperse systems robustly mix on large lengthscales. On the other hand, shape disparity generates slight demixing over a few cell diameters, a phenomenon we term micro-demixing. This result can be understood by examining the differential energy barriers for neighbor exchanges (T1 transitions). Experiments with mixtures of wild-type and E-cadherin-deficient keratinocytes on a substrate are consistent with the predicted phenomenon of micro-demixing, which biology may exploit to create subtle patterning. The robustness of mixing at large scales, however, suggests that despite some differences in cell shape and size, progenitor cells can readily mix throughout a developing tissue until acquiring means of recognizing cells of different types., Comment: 32 pages, 17 figures

We present a comprehensive theory of the dynamics and fluctuations of a two-dimensional suspension of polar active particles in an incompressible fluid confined to a substrate. We show that, depending on the sign of a single parameter, a state with polar orientational order is anomalously stable (or anomalously unstable), with a nonzero relaxation (or growth) rate for angular fluctuations at zero wavenumber. This screening of the broken-symmetry mode in the stable state does lead to conventional rather than giant number fluctuations as argued by Bricard et al., Nature ${\bf 503}$, 95 (2013), but their bend instability in a splay-stable flock does not exist and the polar phase has long-range order in two dimensions. Our theory also describes confined three-dimensional thin-film suspensions of active polar particles as well as dense compressible active polar rods, and predicts a flocking transition without a banding instability

We study active run-and-tumble particles with an additional two-state internal variable characterizing their motile or non-motile state. Motile particles change irreversibly into non-motile ones upon collision with a non-motile particle. The system evolves towards an absorbing state where all particles are non-motile. We initialize the system with one non-motile particles in a bath of motile ones and study numerically the kinetics of relaxation to absorbing state and its structure as function of the density of the initial bath of motile particles and of their tumbling rate. We find a crossover from fractal aggregates at low density to homogeneous ones at high density. The persistence of single-particle dynamics as quantified by the tumbling rate pushes this crossover to higher density and can be used to tune the porosity of the aggregate. At the lowest density the fractal dimension of the aggregate approaches that obtained in single-particle diffusion limited aggregation. Our results could be exploited for the design of structures of desired porosity. The model is a first step towards the study of the collective dynamics of active particles that can exchange biological information., Comment: 6 pages, 4 figures

The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small scale forcing of hydrodynamic flow by microswimmers, we study the properties of a dense supensions of microswimmers in two dimensions, where the conservation of enstrophy can drive an inverse cascade through which energy is accumulated on the largest scales. We find that the dynamical and statistical properties of the flow show a sharp transition to the formation of vortices at the largest length scale. The results show that 2d bacterial and hydrodynamic turbulence are separated by a subcritical phase transition., Comment: postprint version

Collective cell migration plays a central role in tissue development, morphogenesis, wound repair and cancer progression. With the growing realization that physical forces mediate cell motility in development and physiology, a key biological question is how cells integrate molecular activities for force generation on multicellular scales. In this review we discuss recent advances in modeling collective cell migration using quantitative tools and approaches rooted in soft matter physics. We focus on theoretical models of cell aggregates as continuous active media, where the feedback between mechanical forces and regulatory biochemistry gives rise to rich collective dynamical behavior. This class of models provides a powerful predictive framework for the physiological dynamics that underlies many developmental processes, where cells need to collectively migrate like a viscous fluid to reach a target region, and then stiffen to support mechanical stresses and maintain tissue cohesion., Comment: Review article; Figures reproduced with permission from relevant author/publisher