Mesoscale pattern formation of self-propelled rods with velocity reversal

Activity-induced phase transition and coarsening dynamics in dry apolar active nematics

Using the Lebwohl-Lasher interaction for reciprocal local alignment, we present a comprehensive phase diagram for a dry, apolar, active nematic system using its stochastic off-lattice dynamics. The nematic-isotropic transition in this system is first-order and occurs alongside a fluctuation-dominated phase separation. Our phase diagram identifies three distinct regions based on activity and orientational noise relative to alignment strength: a homogeneous isotropic phase, a nematic phase with giant density fluctuations, and a coexistence region. Using mean-field analysis and hydrodynamic theory, we demonstrate that reciprocal interactions lead to a density fluctuation-induced first-order transition and derive a phase boundary consistent with numerical results. Quenching from the isotropic to nematic phase reveals coarsening dynamics where nematic ordering precedes particle clustering. Both the nematic and density fields exhibit similar scaling behaviors, exhibiting dynamic exponents zS ≈ 2.5 and zρ ≈ 2.34, consistently falling within the range of 2 and 3.

Environment-Stored Memory in Active Nematics and Extra-Cellular Matrix Remodeling

Many active systems display nematic order, while interacting with their environment. In this Letter, we show theoretically how environment-stored memory acts an effective external field that aligns active nematics. The coupling to the environment leads to substantial modifications of the known phase diagram and dynamics of active nematics, including nematic order at arbitrarily low densities and arrested domain coarsening. We are motivated mainly by cells that remodel fibers in their extra-cellular matrix (ECM), while being directed by the fibers during migration. Our predictions indicate that remodeling promotes cellular and ECM alignment, and possibly limits the range of ordered ECM domains, in accordance with recent experiments.

Aspect-ratio-dependent void formation in active rhomboidal and elliptical particle systems

We execute a numerical simulation of active nematics with particles interacting by an excluded-volume effect. Systems with rhomboidal particles and with elliptical particles are considered in order to investigate the effect of the direct contact of particles. In our simulation, the void regions, where the local number density is almost zero, appear in both systems when the aspect ratio of the particles is high. We focus on the relationship between the void regions and the particle orientation of the bulk. The particle number density, particle orientation, topological defects, and void regions are analyzed for different aspect ratios in both systems. The systems with rhomboidal particles have characteristic void sizes, which increase with an increase in the aspect ratio. In contrast, the distribution of the void-region size in the systems with elliptical particles is broad. The present results suggest that the void size in the systems with rhomboidal particles is determined by the correlation length of the particle orientational field around the void regions, while the void size might be determined by the system size in the systems with elliptical particles.

Derivation and simulation of a computational model of active cell populations: How overlap avoidance, deformability, cell-cell junctions and cytoskeletal forces affect alignment

Collective alignment of cell populations is a commonly observed phenomena in biology. An important example are aligning fibroblasts in healthy or scar tissue. In this work we derive and simulate a mechanistic agent-based model of the collective behaviour of actively moving and interacting cells, with a focus on understanding collective alignment. The derivation strategy is based on energy minimisation. The model ingredients are motivated by data on the behaviour of different populations of aligning fibroblasts and include: Self-propulsion, overlap avoidance, deformability, cell-cell junctions and cytoskeletal forces. We find that there is an optimal ratio of self-propulsion speed and overlap avoidance that maximises collective alignment. Further we find that deformability aids alignment, and that cell-cell junctions by themselves hinder alignment. However, if cytoskeletal forces are transmitted via cell-cell junctions we observe strong collective alignment over large spatial scales.

Hierarchical defect-induced condensation in active nematics

Topological defects play a central role in the formation and organization of various biological systems. Historically, such nonequilibrium defects have been mainly studied in the context of homogeneous active nematics. Phase-separated systems, in turn, are known to form dense and dynamic nematic bands, but typically lack topological defects. In this paper, we use agent-based simulations of weakly aligning, self-propelled polymers and demonstrate that contrary to the existing paradigm phase-separated active nematics form -1/2 defects. Moreover, these defects, emerging due to interactions among dense nematic bands, constitute a novel second-order collective state. We investigate the morphology of defects in detail and find that their cores correspond to a strong increase in density, associated with a condensation of nematic fluxes. Unlike their analogs in homogeneous systems, such condensed defects form and decay in a different way and do not involve positively charged partners. We additionally observe and characterize lateral arc-like structures that separate from a band's bulk and move in transverse direction. We show that the key control parameters defining the route from stable bands to the coexistence of dynamic lanes and defects are the total density of particles and their path persistence length. We introduce a hydrodynamic theory that qualitatively recapitulates all the main features of the agent-based model, and use it to show that the emergence of both defects and arcs can be attributed to the same anisotropic active fluxes. Finally, we present a way to artificially engineer and position defects, and speculate about experimental verification of the provided model.

Quadrupolar active stress induces exotic patterns of defect motion in compressible active nematics

A wide range of living and artificial active matter exists in close contact with substrates and under strong confinement, where in addition to dipolar active stresses, quadrupolar active stresses can become important. Here, we numerically investigate the impact of quadrupolar non-equilibrium stresses on the emergent patterns of self-organisation in non-momentum conserving active nematics. Our results reveal that beyond having stabilising effects, the quadrupolar active forces can induce various modes of topological defect motion in active nematics. In particular, we find the emergence of both polar and nematic ordering of the defects, as well as new patterns of self-organisation that comprise topological defect chains and transient topological defect asters. The results contribute to further understanding of emergent patterns of collective motion and non-equilibrium self-organisation in active matter.

Susceptibility of Orientationally Ordered Active Matter to Chirality Disorder

We investigate the susceptibility of long-range ordered phases of two-dimensional dry aligning active matter to population disorder, taken in the form of a distribution of intrinsic individual chiralities. Using a combination of particle-level models and hydrodynamic theories derived from them, we show that while in finite systems all ordered phases resist a finite amount of such chirality disorder, the homogeneous ones (polar flocks and active nematics) are unstable to any amount of disorder in the infinite-size limit. On the other hand, we find that the inhomogeneous solutions of the coexistence phase (bands) may resist a finite amount of chirality disorder even asymptotically.

Pattern-induced local symmetry breaking in active-matter systems

The emergence of macroscopic order and patterns is a central paradigm in systems of (self-)propelled agents and a key component in the structuring of many biological systems. The relationships between the ordering process and the underlying microscopic interactions have been extensively explored both experimentally and theoretically. While emerging patterns often show one specific symmetry (e.g., nematic lane patterns or polarized traveling flocks), depending on the symmetry of the alignment interactions patterns with different symmetries can apparently coexist. Indeed, recent experiments with an actomysin motility assay suggest that polar and nematic patterns of actin filaments can interact and dynamically transform into each other. However, theoretical understanding of the mechanism responsible remains elusive. Here, we present a kinetic approach complemented by a hydrodynamic theory for agents with mixed alignment symmetries, which captures the experimentally observed phenomenology and provides a theoretical explanation for the coexistence and interaction of patterns with different symmetries. We show that local, pattern-induced symmetry breaking can account for dynamically coexisting patterns with different symmetries. Specifically, in a regime with moderate densities and a weak polar bias in the alignment interaction, nematic bands show a local symmetry-breaking instability within their high-density core region, which induces the formation of polar waves along the bands. These instabilities eventually result in a self-organized system of nematic bands and polar waves that dynamically transform into each other. Our study reveals a mutual feedback mechanism between pattern formation and local symmetry breaking in active matter that has interesting consequences for structure formation in biological systems.

A particle-field approach bridges phase separation and collective motion in active matter

Whereas self-propelled hard discs undergo motility-induced phase separation, self-propelled rods exhibit a variety of nonequilibrium phenomena, including clustering, collective motion, and spatio-temporal chaos. In this work, we present a theoretical framework representing active particles by continuum fields. This concept combines the simplicity of alignment-based models, enabling analytical studies, and realistic models that incorporate the shape of self-propelled objects explicitly. By varying particle shape from circular to ellipsoidal, we show how nonequilibrium stresses acting among self-propelled rods destabilize motility-induced phase separation and facilitate orientational ordering, thereby connecting the realms of scalar and vectorial active matter. Though the interaction potential is strictly apolar, both, polar and nematic order may emerge and even coexist. Accordingly, the symmetry of ordered states is a dynamical property in active matter. The presented framework may represent various systems including bacterial colonies, cytoskeletal extracts, or shaken granular media.

Interacting self-propelled particles exhibit phase separation or collective motion depending on particle shape. A unified theory connecting these paradigms represents a major challenge in active matter, which the authors address here by modeling active particles as continuum fields.

Rules of collective migration: from the wildebeest to the neural crest

Collective migration, the movement of groups in which individuals affect the behaviour of one another, occurs at practically every scale, from bacteria up to whole species' populations. Universal principles of collective movement can be applied at all levels. In this review, we will describe the rules governing collective motility, with a specific focus on the neural crest, an embryonic stem cell population that undergoes extensive collective migration during development. We will discuss how the underlying principles of individual cell behaviour, and those that emerge from a supracellular scale, can explain collective migration. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Modelling collective cell motion: are on- and off-lattice models equivalent?

Biological processes, such as embryonic development, wound repair and cancer invasion, or bacterial swarming and fruiting body formation, involve collective motion of cells as a coordinated group. Collective cell motion of eukaryotic cells often includes interactions that result in polar alignment of cell velocities, while bacterial patterns typically show features of apolar velocity alignment. For analysing the population-scale effects of these different alignment mechanisms, various on- and off-lattice agent-based models have been introduced. However, discriminating model-specific artefacts from general features of collective cell motion is challenging. In this work, we focus on equivalence criteria at the population level to compare on- and off-lattice models. In particular, we define prototypic off- and on-lattice models of polar and apolar alignment, and show how to obtain an on-lattice from an off-lattice model of velocity alignment. By characterizing the behaviour and dynamical description of collective migration models at the macroscopic level, we suggest the type of phase transitions and possible patterns in the approximative macroscopic partial differential equation descriptions as informative equivalence criteria between on- and off-lattice models. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Collective Dynamics of Model Pili-Based Twitcher-Mode Bacilliforms

Pseudomonas aeruginosa, like many bacilliforms, are not limited only to swimming motility but rather possess many motility strategies. In particular, twitching-mode motility employs hair-like pili to transverse moist surfaces with a jittery irregular crawl. Twitching motility plays a critical role in redistributing cells on surfaces prior to and during colony formation. We combine molecular dynamics and rule-based simulations to study twitching-mode motility of model bacilliforms and show that there is a critical surface coverage fraction at which collective effects arise. Our simulations demonstrate dynamic clustering of twitcher-type bacteria with polydomains of local alignment that exhibit spontaneous correlated motions, similar to rafts in many bacterial communities.

The 2020 motile active matter roadmap

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

Extracellular matrix anisotropy is determined by TFAP2C-dependent regulation of cell collisions

The isotropic or anisotropic organization of biological extracellular matrices has important consequences for tissue function. We study emergent anisotropy using fibroblasts that generate varying degrees of matrix alignment from uniform starting conditions. This reveals that the early migratory paths of fibroblasts are correlated with subsequent matrix organization. Combined experimentation and adaptation of Vicsek modelling demonstrates that the reorientation of cells relative to each other following collision plays a role in generating matrix anisotropy. We term this behaviour 'cell collision guidance'. The transcription factor TFAP2C regulates cell collision guidance in part by controlling the expression of RND3. RND3 localizes to cell-cell collision zones where it downregulates actomyosin activity. Cell collision guidance fails without this mechanism in place, leading to isotropic matrix generation. The cross-referencing of alignment and TFAP2C gene expression signatures against existing datasets enables the identification and validation of several classes of pharmacological agents that disrupt matrix anisotropy.

Understanding Dense Active Nematics from Microscopic Models

We study dry, dense active nematics at both particle and continuous levels. Specifically, extending the Boltzmann-Ginzburg-Landau approach, we derive well-behaved hydrodynamic equations from a Vicsek-style model with nematic alignment and pairwise repulsion. An extensive study of the phase diagram shows qualitative agreement between the two levels of description. We find in particular that the dynamics of topological defects strongly depends on parameters and can lead to "arch" solutions forming a globally polar, smecticlike arrangement of Néel walls. We show how these configurations are at the origin of the defect ordered states reported previously. This work offers a detailed understanding of the theoretical description of dense active nematics directly rooted in their microscopic dynamics.

Dry active turbulence in a model for microtubule-motor mixtures

We study the dynamics and phase behaviour of a dry suspension of microtubules and molecular motors. We obtain a set of continuum equations by rigorously coarse graining a microscopic model where motor-induced interactions lead to parallel or antiparallel ordering. Through numerical simulations, we show that this model generically creates either stable stripes, or a never-settling pattern where stripes periodically form, rotate and then split up. We derive a minimal model which displays the same instability as the full model, and clarifies the underlying physical mechanism. The necessary ingredients are an extensile flux arising from microtubule sliding and an interfacial torque favouring ordering along density gradients. We argue that our minimal model unifies various previous observations of chaotic behaviour in dry active matter into a general universality class.

Self-Driven Phase Transitions Drive Myxococcus xanthus Fruiting Body Formation

Combining high-resolution single cell tracking experiments with numerical simulations, we show that starvation-induced fruiting body formation in Myxococcus xanthus is a phase separation driven by cells that tune their motility over time. The phase separation can be understood in terms of cell density and a dimensionless Péclet number that captures cell motility through speed and reversal frequency. Our work suggests that M. xanthus takes advantage of a self-driven nonequilibrium phase transition that can be controlled at the single cell level.

Shared behavioral mechanisms underlie C. elegans aggregation and swarming

In complex biological systems, simple individual-level behavioral rules can give rise to emergent group-level behavior. While collective behavior has been well studied in cells and larger organisms, the mesoscopic scale is less understood, as it is unclear which sensory inputs and physical processes matter a priori. Here, we investigate collective feeding in the roundworm C. elegans at this intermediate scale, using quantitative phenotyping and agent-based modeling to identify behavioral rules underlying both aggregation and swarming-a dynamic phenotype only observed at longer timescales. Using fluorescence multi-worm tracking, we quantify aggregation in terms of individual dynamics and population-level statistics. Then we use agent-based simulations and approximate Bayesian inference to identify three key behavioral rules for aggregation: cluster-edge reversals, a density-dependent switch between crawling speeds, and taxis towards neighboring worms. Our simulations suggest that swarming is simply driven by local food depletion but otherwise employs the same behavioral mechanisms as the initial aggregation.

A statistical physics view of swarming bacteria

Bacterial swarming is a collective mode of motion in which cells migrate rapidly over surfaces, forming dynamic patterns of whirls and jets. This review presents a physical point of view of swarming bacteria, with an emphasis on the statistical properties of the swarm dynamics as observed in experiments. The basic physical principles underlying the swarm and their relation to contemporary theories of collective motion and active matter are reviewed and discussed in the context of the biological properties of swarming cells. We suggest a paradigm according to which bacteria have optimized some of their physical properties as a strategy for rapid surface translocation. In other words, cells take advantage of favorable physics, enabling efficient expansion that enhances survival under harsh conditions.

Dynamical subclasses of dry active nematics

We show that the dominant mode of alignment plays an important role in dry active nematics, leading to two dynamical subclasses defined by the nature of the instability of the nematic bands that characterize, in these systems, the coexistence phase separating the isotropic and fluctuating nematic states. In addition to the well-known instability inducing long undulations along the band, another stronger instability leading to the breakup of the band in many transversal segments may arise. We elucidate the origin of this strong instability for a realistic model of self-propelled rods and determine the high-order nonlinear terms responsible for it at the hydrodynamic level.

Active nematics

Active matter extracts energy from its surroundings at the single particle level and transforms it into mechanical work. Examples include cytoskeleton biopolymers and bacterial suspensions. Here, we review experimental, theoretical and numerical studies of active nematics - a type of active system that is characterised by self-driven units with elongated shape. We focus primarily on microtubule-kinesin mixtures and the hydrodynamic theories that describe their properties. An important theme is active turbulence and the associated motile topological defects. We discuss ways in which active turbulence may be controlled, a pre-requisite to harvesting energy from active materials, and we consider the appearance, and possible implications, of active nematics and topological defects to cellular systems and biological processes.

Twist-induced crossover from two-dimensional to three-dimensional turbulence in active nematics

While studies of active nematics in two dimensions have shed light on various aspects of the flow regimes and topology of active matter, three-dimensional properties of topological defects and chaotic flows remain unexplored. By confining a film of active nematics between two parallel plates, we use continuum simulations and analytical arguments to demonstrate that the crossover from quasi-two-dimensional (quasi-2D) to three-dimensional (3D) chaotic flows is controlled by the morphology of the disclination lines. For small plate separations, the active nematic behaves as a quasi-2D material, with straight topological disclination lines spanning the height of the channel and exhibiting effectively 2D active turbulence. Upon increasing channel height, we find a crossover to 3D chaotic flows due to the contortion of disclinations above a critical activity. Above this critical activity highly contorted disclination lines and disclination loops are formed. We further show that these contortions are engendered by twist perturbations producing a sharp change in the curvature of disclinations.

Self-Propelled Particles with Velocity Reversals and Ferromagnetic Alignment: Active Matter Class with Second-Order Transition to Quasi-Long-Range Polar Order

We introduce and study in two dimensions a new class of dry, aligning active matter that exhibits a direct transition to orientational order, without the phase-separation phenomenology usually observed in this context. Characterized by self-propelled particles with velocity reversals and a ferromagnetic alignment of polarities, systems in this class display quasi-long-range polar order with continuously varying scaling exponents, yet a numerical study of the transition leads to conclude that it does not belong to the Berezinskii-Kosterlitz-Thouless universality class but is best described as a standard critical point with an algebraic divergence of correlations. We rationalize these findings by showing that the interplay between order and density changes the role of defects.

Giant Kovacs-Like Memory Effect for Active Particles

Dynamical properties of self-propelled particles obeying a bounded confidence rule are investigated by means of kinetic theory and agent-based simulations. While memory effects are observed in disordered systems, we show that they also occur in active matter systems. In particular, we find that the system exhibits a giant Kovacs-like memory effect that is much larger than predicted by a generic linear theory. Based on a separation of time scales we develop a nonlinear theory to explain this effect. We apply this theory to driven granular gases and propose further applications to spin glasses.

On the mechanism of long-range orientational order of fibroblasts

Long-range alignment ordering of fibroblasts have been observed in the vicinity of cancerous tumors and can be recapitulated with in vitro experiments. However, the mechanisms driving their ordering are not understood. Here, we show that local collision-driven nematic alignment interactions among fibroblasts are insufficient to explain observed long-range alignment. One possibility is that there exists another orientation field coevolving with the cells and reinforcing their alignment. We propose that this field reflects the mechanical cross-talk between the fibroblasts and the underlying fibrous material on which they move. We show that this long-range interaction can give rise to high nematic order and to the observed patterning of the cancer microenvironment.