From phase to microphase separation in flocking models: the essential role of nonequilibrium fluctuations

Hydrodynamically Enhanced Brownian Motion in Flowing Polymer Solutions

Diffusion of a molecule in solution typically occurs via thermally driven Brownian motion, with solvent collisions leading to a random-walk trajectory for the solute. This physical principle guides our understanding of molecular transport in a wide variety of situations, ranging from protein diffusion in biological systems and mixing in solution processes to charge transport in polyelectrolyte solutions. Thermal diffusion represents a "speed limit" for molecular transport, which can typically only be surpassed by imposing a directional, external field. The other way particles expedite diffusion is via self-propulsion. This "active Brownian motion" is famously seen in some single-cell organisms and can also be shown in some colloidal systems, but because it requires self-propulsion, this is not seen at the molecular level. We show that it is possible to dramatically increase the diffusion of small molecules in a way that mimics active Brownian motion, instead driven by the disturbance flows of highly stretched polymers in strong flows to induce propulsion at a distance. We use molecular simulations that account for these hydrodynamic disturbances to demonstrate that it is possible to increase the effective diffusion constant by more than an order of magnitude, and we provide a mechanistic model for how the interplay of polymer concentration, flow-induced polymer stretching, and chain length gives rise to hydrodynamically enhanced Brownian motion. This effect has important implications for molecular transport, and we show that strong flows and a low concentration of stretched polymers can be used to promote rapid diffusion.

Banded phases in topological flocks

Flocking phase transitions arise in many aligning active soft matter systems, and an interesting question concerns the role of "topological" vs. "metric" interactions on these transitions. While recent theoretical work suggests that the order-disorder transition in these polar aligning models is universally first order, numerical studies have suggested that topological models may instead have a continuous transition. Some recent simulations have found that some variations of topologically interacting flocking agents have a discontinuous transition, but unambiguous observations of phase coexistence using common Voronoi-based alignment remains elusive. In this work, we use a custom GPU-accelerated simulation package to perform million-particle-scale simulations of a Voronoi-Vicsek model in which alignment interactions stem from an XY-like Hamiltonian. By accessing such large systems on appropriately long time scales and in the time-continuous limit, we are able to show a regime of stable phase coexistence between the ordered and disordered phases, confirming the discontinuous nature of this transition in the thermodynamic limit.

Collective dynamics of intelligent active Brownian particles with visual perception and velocity alignment in 3D: spheres, rods, and worms

Many living systems, such as birds and fish, exhibit collective behaviors like flocking and swarming. Recently, an experimental system of active colloidal particles has been developed, where the motility of each particle is adjusted based on its visual detection of surrounding particles. These particles with visual-perception-dependent motility exhibit group formation and cohesion. Inspired by these behaviors, we investigate intelligent active Brownian particles (iABPs) equipped with visual perception and velocity alignment in three dimensions using computer simulations. The visual-perception-based self-steering describes the tendency of iABPs to move toward the center of mass of particles within their visual cones, while velocity alignment encourages alignment with neighboring particles. We examine how the behavior varies with the visual cone angle θ, self-propulsion speed (Péclet number Pe), and the interaction strengths of velocity alignment (Ωa) and visual-based self-steering (Ωv). Our findings show that spherical iABPs form dense clusters, worm-like clusters, milling behaviors, and dilute-gas phases, consistent with 2D studies. By reducing the simulation box size, we observe additional structures like band-like clusters and dense baitball formations. Additionally, rod-like iABPs form band-like, worm-like, radiating, and helical structures, while iABP worms exhibit band-like, streamlined, micellar-like and entangled structures. Many of these patterns resemble collective behaviors in nature, such as ant milling, fish baitballs, and worm clusters. Advances in synthetic techniques could enable nanorobots with similar capabilities, offering insights into multicellular systems through active matter.

Collective motion of energy depot active disks

In the present work we have studied collectives of active disks with an energy depot, moving in the two-dimensional plane and interacting via an excluded volume. The energy depot accounts for the extraction of energy taking place at the level of each particle in order to perform self-propulsion, included in an underdamped Langevin dynamics. We show that this model undergoes a flocking transition, exhibiting some of the key features of the Vicsek model, namely, band formation and giant number fluctuations. These bands, either single or multiple, are dense and very strongly polarised propagating structures. Large density bands disappear as the activity is further increased, eventually reaching a homogeneous polar state. We unravel an effective alignment interaction at the level of two-particle collisions that can be controlled by activity and gives rise to flocking at large scales.

Emergent flocking in mixtures of antialigning self-propelled particles

We observe a flocking mechanism, the emergence of a state with global polar order, in mixed systems of self-propelled particles with purely antialigning interactions, i.e., the ground state for any pair of particles is to be opposedly oriented. In binary mixtures, we find that flocking can be realized by cross-species antialigning that is dominant compared to intraspecies antialignment. While the key mechanism can be understood within a mean-field description, beyond mean-field we develop an asymptotically exact Boltzmann-scattering theory from first principles. This theory yields analytical predictions for the flocking transition and shows excellent quantitative agreement with simulations of dilute systems. For large systems, we find either microphase separation or static patterns with patches or stripes that carry different polarization orientations.

Sedimentation and structure of squirmer suspensions under gravity

The effect of gravity on the collective motion of living microswimmers, such as bacteria and micro-algae, is pivotal to unravel not only bio-convection patterns but also the settling of bacterial biofilms on solid surfaces. In this work, we investigate suspensions of microswimmers under the influence of a gravitational field and hydrodynamics, simulated via the dissipative particle dynamics (DPD) coarse-grained model. We first study the collective sedimentation of passive colloids and microswimmers of the puller and pusher types upon increasing the imposed gravitational field and compare them with previous results. Once sedimentation occurs, we observe that, as the gravitational field increases, the bottom layer undergoes a transition to an ordered state compatible with a hexagonal crystal. In comparison with passive colloids, both pullers and pushers easily rearrange at the bottom layer to anneal defects. Specifically, pullers are better than pushers in preserving the hexagonal order of the bottom mono-layer at high gravitational fields.

Motility-Induced Pinning in Flocking System with Discrete Symmetry

We report a motility-induced pinning transition in the active Ising model for a self-propelled particle system with discrete symmetry. This model was known to exhibit a liquid-gas type flocking phase transition, but a recent study reveals that the polar order is metastable due to droplet excitation. Using extensive Monte Carlo simulations, we demonstrate that, for an intermediate alignment interaction strength, the steady state is characterized by traveling local domains, which renders the polar order short-ranged in both space and time. We further demonstrate that interfaces between colliding domains become pinned as the alignment interaction strength increases. A resonating back-and-forth motion of individual self-propelled particles across interfaces is identified as a mechanism for the pinning. We present a numerical phase diagram for the motility-induced pinning transition, and an approximate analytic theory for the growth and shrink dynamics of pinned interfaces. Our results show that pinned interfaces grow to a macroscopic size preventing the polar order in the regime where the particle diffusion rate is sufficiently smaller than the self-propulsion rate. The growth behavior in the opposite regime and its implications on the polar order remain unresolved and require further investigation.

Strong Casimir-like Forces in Flocking Active Matter

Confining in space the equilibrium fluctuations of statistical systems with long-range correlations is known to result into effective forces on the boundaries. Here we demonstrate the occurrence of Casimir-like forces in the nonequilibrium context provided by flocking active matter. In particular, we consider a system of aligning self-propelled particles in two spatial dimensions that are transversally confined by reflecting or partially reflecting walls. We show that in the ordered flocking phase this confined active vectorial fluid is characterized by extensive boundary layers, as opposed to the finite ones usually observed in confined scalar active matter. Moreover, a finite-size, fluctuation-induced contribution to the pressure on the wall emerges, which decays slowly and algebraically upon increasing the distance between the walls. We explain our findings-which display a certain degree of universality-within a hydrodynamic description of the density and velocity fields.

Tuning collective behaviour in zebrafish with genetic modification

Zebrafish collective behaviour is widely used to assess their physical and mental state, serving as a valuable tool to assess the impact of ageing, disease genetics, and the effect of drugs. The essence of these macroscopic phenomena can be represented by active matter models, where the individuals are abstracted as interactive self-propelling agents. The behaviour of these agents depends on a set of parameters in a manner reminiscent of those between the constituents of physical systems. In a few cases, the system may be controlled at the level of the individual constituents such as the interactions between colloidal particles, or the enzymatic behaviour of de novo proteins. Usually, however, while the collective behaviour may be influenced by environmental factors, it typically cannot be changed at will. Here, we challenge this scenario in a biological context by genetically modifying zebrafish. We thus demonstrate the potential of genetic modification in the context of controlling the collective behaviour of biological active matter systems at the level of the constituents, rather than externally. In particular, we probe the effect of the lack of col11a2 gene in zebrafish, which causes the early onset of osteoarthritis. The resulting col11a2 -/- zebrafish exhibited compromised vertebral column properties, bent their body less while swimming, and took longer to change their orientations. Surprisingly, a group of 25 mutant fish exhibited more orderly collective motion than the wildtype. We show that the collective behaviour of wildtype and col11a2 -/- zebrafish are captured with a simple active matter model, in which the mutant fish are modelled by self-propelling agents with a higher orientational noise on average. In this way, we demonstrate the possibility of tuning a biological system, changing the state space it occupies when interpreted with a simple active matter model.

Self-organization of anti-aligning active particles: Waving pattern formation and chaos

Recently, it has been shown that purely anti-aligning interaction between active particles may induce a finite wavelength instability. The formed patterns display intricate spatiotemporal dynamics, suggesting the presence of chaos. Here, we propose a quasi-one-dimensional simplification of the particle interaction model. This simplified model allows us to deduce amplitude equations that describe the collective motion of the active entities. We show that these equations exhibit chaotic orbits. Furthermore, via direct numerical simulations of the particle's system, we discuss the pertinence of these amplitude equations approach for describing the particle's self-coordinated motions.

Collective Hall current in chiral active fluids: Coupling of phase and mass transport through traveling bands

Active fluids composed of constituents that are constantly driven away from thermal equilibrium can support spontaneous currents and can be engineered to have unconventional transport properties. Here, we report the emergence of (meta)stable traveling bands in computer simulations of aligning circle swimmers. These bands are different from polar flocks and, through coupling phase with mass transport, induce a bulk particle current with a component perpendicular to the propagation direction, thus giving rise to a collective Hall (or Magnus) effect. Traveling bands require sufficiently small orbits and undergo a discontinuous transition into a synchronized state with transient polar clusters for large orbital radii. Within a minimal hydrodynamic theory, we show that the bands can be understood as nondispersive soliton solutions fully accounting for the numerically observed properties.

Three-state active lattice gas: A discrete Vicsek-like model with excluded volume

We study a discrete-space model of active matter with excluded volume. Particles are restricted to the sites of a triangular lattice and can assume one of three orientations. Varying the density and noise intensity, Monte Carlo simulations reveal a variety of spatial patterns. Ordered states occur in the form of condensed structures, which (away from the full occupancy limit) coexist with a low-density vapor. The condensed structures feature low particle mobility, particularly those that wrap the system via the periodic boundaries. As the noise intensity is increased, dense structures give way to a disordered phase. We characterize the parameter values associated with the condensed phases and perform a detailed study of the order-disorder transition at (i) full occupation and (ii) a density of 0.1. In the former case, the model possesses the same symmetry as the three-state Potts model and exhibits a continuous phase transition, as expected, with critical exponents consistent with those of the associated Potts model. In the low-density case, the transition is clearly discontinuous, with a strong dependence of the final state upon the initial configuration, hysteresis, and nonmonotonic dependence of the Binder cumulant upon noise intensity.

Ordering kinetics in the active Ising model

We undertake a numerical study of the ordering kinetics in the two-dimensional (2D) active Ising model (AIM), a discrete flocking model with a conserved density field coupled to a nonconserved magnetization field. We find that for a quench into the liquid-gas coexistence region and in the ordered liquid region, the characteristic length scale of both the density and magnetization domains follows the Lifshitz-Cahn-Allen growth law, R(t)∼t^{1/2}, consistent with the growth law of passive systems with scalar order parameter and nonconserved dynamics. The system morphology is analyzed with the two-point correlation function and its Fourier transform, the structure factor, which conforms to the well-known Porod's law, a manifestation of the coarsening of compact domains with smooth boundaries. We also find the domain growth exponent unaffected by different noise strengths and self-propulsion velocities of the active particles. However, transverse diffusion is found to play the most significant role in the growth kinetics of the AIM. We extract the same growth exponent by solving the hydrodynamic equations of the AIM.

Finite-size scaling in kinetics of phase separation in certain models of aligning active particles

To study the kinetics of phase separation in active matter systems, we consider models that impose a Vicsek-type self-propulsion rule on otherwise passive particles interacting via the Lennard-Jones potential. Two types of kinetics are of interest: one conserves the total momentum of all the constituents and the other does not. We carry out molecular dynamics simulations to obtain results on structural, growth, and aging properties. Results from our studies, with various finite boxes, show that there exist scalings with respect to the system sizes, in both the latter quantities, as in the standard passive cases. We have exploited this scaling picture to accurately estimate the corresponding exponents, in the thermodynamically large system size limit, for power-law time dependences. It is shown that certain analytical functions describe the behavior of these quantities quite accurately, including the finite-size limits. Our results demonstrate that even though the conservation of velocity has at best weak effects on the dynamics of evolution in the thermodynamic limit, the finite-size behavior is strongly influenced by the presence (or the absence) of it.

Nematic Torques in Scalar Active Matter: When Fluctuations Favor Polar Order and Persistence

We study the impact of nematic alignment on scalar active matter in the disordered phase. We show that nematic torques control the emergent physics of particles interacting via pairwise forces and can either induce or prevent phase separation. The underlying mechanism is a fluctuation-induced renormalization of the mass of the polar field that generically arises from nematic torques. The correlations between the fluctuations of the polar and nematic fields indeed conspire to increase the particle persistence length, contrary to what phenomenological computations predict. This effect is generic and our theory also quantitatively accounts for how nematic torques enhance particle accumulation along confining boundaries and opposes demixing in mixtures of active and passive particles.

Collective motion of chiral particles in complex noise environments

Collective motion of chiral particles in complex noise environments is investigated based on the Vicsek model. In the model, we added chirality, along with complex noise, affecting particles clustering motion. Particles can only avoid noise interference in a specific channel, and this consideration is more realistic due to the complexity of the environment. Via simulations, we find that the channel proportion, p, critically influences chiral particle synchronization. Specifically, we observe a disorder-order transition at critical [Formula: see text], only when [Formula: see text], the system can achieve global synchronization. Combined with our definition of spatial distribution parameter and observation of the model, the reason is that particles begin to escape from the noise region under the influence of complex noise. In addition, the value of [Formula: see text] increases linearly with velocity, while it decreases monotonically with the increase in chirality and interaction radius. Interestingly, an appropriate noise amplitude minimizes [Formula: see text]. Our findings may inspire novel strategies to manipulate self-propelled particles of distinct chirality to achieve desired spatial migration and global synchronization.

Anti-aligning interaction between active particles induces a finite wavelength instability: The dancing hexagons

By considering a simple model for self-propelled particle interaction, we show that anti-aligning forces induce a finite wavelength instability. Consequently, the system exhibits pattern formation. The formed pattern involves, let us say, a choreographic movement of the active entities. At the level of particle density, the system oscillates between a stripe pattern and a hexagonal one. The underlying dynamics of these density oscillations consists of two counterpropagating and purely hexagonal traveling waves. They are assembling and disassembling a global hexagonal structure and a striped lineup of particles. This self-assembling process becomes quite erratic for long-time simulations, seeming aperiodic.

Nonequilibrium phase transitions in a Brownian p-state clock model

We introduce a Brownian p-state clock model in two dimensions and investigate the nature of phase transitions numerically. As a nonequilibrium extension of the equilibrium lattice model, the Brownian p-state clock model allows spins to diffuse randomly in the two-dimensional space of area L^{2} under periodic boundary conditions. We find three distinct phases for p>4: a disordered paramagnetic phase, a quasi-long-range-ordered critical phase, and an ordered ferromagnetic phase. In the intermediate critical phase, the magnetization order parameter follows a power-law scaling m∼L^{-β[over ̃]}, where the finite-size scaling exponent β[over ̃] varies continuously. These critical behaviors are reminiscent of the double Berezinskii-Kosterlitz-Thouless (BKT) transition picture of the equilibrium system. At the transition to the disordered phase, the exponent takes the universal value β[over ̃]=1/8, which coincides with that of the equilibrium system. This result indicates that the BKT transition driven by the unbinding of topological excitations is robust against the particle diffusion. On the contrary, the exponent at the symmetry-breaking transition to the ordered phase deviates from the universal value β[over ̃]=2/p^{2} of the equilibrium system. The deviation is attributed to a nonequilibrium effect from the particle diffusion.

A Mutant of Dictyostelium discoideum, KI-Cell, as a Model of Collective Cell Migration Independent of Chemotaxis

Collective cell migration occurs when the orientation of cell polarity is aligned with each other in a group of cells. Such collective polarization depends on a reciprocal process between cell intrinsic mechanisms such as cell-cell adhesion and extracellular guidance mechanism such as wound healing and chemotaxis. As part of its development life cycle, individual single cells of Dictyostelium discoideum exhibit chemotaxis toward cAMP, which is secreted from a certain population of cells. During the formation of multicellular body by chemotaxis-dependent cell aggregation, D. discoideum is also known to relay on multiple cell-cell adhesion mechanisms. In particular, tail-following behavior at the contact site, called contact following of locomotion (CFL), plays a pivotal role on the formation of the multicellular body. However, whether and how CFL alone can lead to a formation of collective behavior was not well understood. KI cell is a mutant of D. discoideum that lacks all chemotactic activity. Yet, it can exhibit the CFL activity and show nontrivial collective cell migration. This mutant provides an excellent model system to analyze the mechanism of the CFL and the macroscopic phenomena brought by the CFL. This chapter describes protocols for using KI cell to understand the biophysics and cell biology behind the collective cell migration induced by CFL.

Coarsening and wavelength selection far from equilibrium: A unifying framework based on singular perturbation theory

Intracellular protein patterns are described by (nearly) mass-conserving reaction-diffusion systems. While these patterns initially form out of a homogeneous steady state due to the well-understood Turing instability, no general theory exists for the dynamics of fully nonlinear patterns. We develop a unifying theory for nonlinear wavelength-selection dynamics in (nearly) mass-conserving two-component reaction-diffusion systems independent of the specific mathematical model chosen. Previous work has shown that these systems support an extremely broad band of stable wavelengths, but the mechanism by which a specific wavelength is selected has remained unclear. We show that an interrupted coarsening process selects the wavelength at the threshold to stability. Based on the physical intuition that coarsening is driven by competition for mass and interrupted by weak source terms that break strict mass conservation, we develop a singular perturbation theory for the stability of stationary patterns. The resulting closed-form analytical expressions enable us to quantitatively predict the coarsening dynamics and the final pattern wavelength. We find excellent agreement with numerical results throughout the diffusion- and reaction-limited regimes of the dynamics, including the crossover region. Further, we show how, in these limits, the two-component reaction-diffusion systems map to generalized Cahn-Hilliard and conserved Allen-Cahn dynamics, therefore providing a link to these two fundamental scalar field theories. The systematic understanding of the length-scale dynamics of fully nonlinear patterns in two-component systems provided here builds the basis to reveal the mechanisms underlying wavelength selection in multicomponent systems with potentially several conservation laws.

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.

Alignment interaction and band formation in assemblies of autochemorepulsive walkers

Chemotaxis refers to the motion of an organism induced by chemical stimuli and is a motility mode shared by many living species that has been developed by evolution to optimize certain biological processes such as foraging or immune response. In particular, autochemotaxis refers to chemotaxis mediated by a cue produced by the chemotactic particle itself. Here, we investigate the collective behavior of autochemotactic particles that are repelled by the cue and therefore migrate preferentially towards low-concentration regions. To this end, we introduce a lattice model inspired by the true self-avoiding walk which reduces to the Keller-Segel model in the continuous limit, for which we describe the rich phase behavior. We first rationalize the chemically mediated alignment interaction between walkers in the limit of stationary concentration fields, and then describe the various large-scale structures that can spontaneously form and the conditions for them to emerge, among which we find stable bands traveling at constant speed in the direction transverse to the band.

Jamming and flocking in the restricted active Potts model

We study the active Potts model with either site occupancy restriction or on-site repulsion to explore jamming and kinetic arrest in a flocking model. The incorporation of such volume exclusion features leads to a surprisingly rich variety of self-organized spatial patterns. While bands and lanes of moving particles commonly occur without or under weak volume exclusion, strong volume exclusion along with low temperature, high activity, and large particle density facilitates jams due to motility-induced phase separation. Through several phase diagrams, we identify the phase boundaries separating the jammed and free-flowing phases and study the transition between these phases which provide us with both qualitative and quantitative predictions of how jamming might be delayed or dissolved. We further formulate and analyze a hydrodynamic theory for the restricted APM which predicts various features of the microscopic model.

Motional consensus of self-propelled particles

The motional consensus of self-propelled particles is studied in both noise-free cases and cases with noise by the standard Vicsek model. In the absence of noise, we propose a simple method, using grid-based technique and defining the normalized variance of the ratio of the number of particles locally to globally, to quantitatively study the movement pattern of the system by the spatial distribution of the particles and the degree of aggregation of particles. It is found that the weaker correlation of velocity leads to larger degree of aggregation of the particles. In the cases with noise, we quantify the competition between velocity alignment and noise by considering the difference of the variety of order parameter result from the velocity alignment and noise. The variation of the effect of noise on motional consensus is non-monotonic for the change of the probability distribution of noise from uniform to non-uniform. Our results may be useful and encourage further efforts in exploring the basic principles of collective motion.

Phase behavior and dynamics in a colloid-polymer mixture under spherical confinement

From studies via molecular dynamics simulations, we report results on structure and dynamics in mixtures of active colloids and passive polymers that are confined inside a spherical container with a repulsive boundary. All interactions in the fully passive limit are chosen in such a way that in equilibrium coexistence between colloid-rich and polymer-rich phases occurs. For most part of the studies the chosen compositions give rise to Janus-like structure: nearly one side of the sphere is occupied by the colloids and the rest by the polymers. This partially wet situation mimics approximately a neutral wall in the fully passive scenario. Following the introduction of a velocity-aligning activity to the colloids, the shape of the polymer-rich domain changes to that of an ellipsoid, around the long axis of which the colloid-rich domain attains a macroscopic angular momentum. In the steady state, the orientation of this axis evolves via diffusion, magnitude of which depends upon the strength of activity, but only weakly.

Flocking without Alignment Interactions in Attractive Active Brownian Particles

Within a simple model of attractive active Brownian particles, we predict flocking behavior and challenge the widespread idea that alignment interactions are necessary to observe this collective phenomenon. Here, we show that even nonaligning attractive interactions can lead to a flocking state. Monitoring the velocity polarization as the order parameter, we reveal the onset of a first-order transition from a disordered phase, characterized by several small clusters, to a flocking phase, where a single flocking cluster is emerging. The scenario is confirmed by studying the spatial connected correlation function of particle velocities, which reveals scale-free behavior in flocking states and exponential-like decay for nonflocking configurations. Our predictions can be tested in microscopic and macroscopic experiments showing flocking, such as animals, migrating cells, and active colloids.

Flocking of two unfriendly species: The two-species Vicsek model

We consider the two-species Vicsek model (TSVM) consisting of two kinds of self-propelled particles, A and B, that tend to align with particles from the same species and to antialign with the other. The model shows a flocking transition that is reminiscent of the original Vicsek model: it has a liquid-gas phase transition and displays micro-phase-separation in the coexistence region where multiple dense liquid bands propagate in a gaseous background. The interesting features of the TSVM are the existence of two kinds of bands, one composed of mainly A particles and one mainly of B particles, the appearance of two dynamical states in the coexistence region: the PF (parallel flocking) state in which all bands of the two species propagate in the same direction, and the APF (antiparallel flocking) state in which the bands of species A and species B move in opposite directions. When PF and APF states exist in the low-density part of the coexistence region they perform stochastic transitions from one to the other. The system size dependence of the transition frequency and dwell times show a pronounced crossover that is determined by the ratio of the band width and the longitudinal system size. Our work paves the way for studying multispecies flocking models with heterogeneous alignment interactions.

Ant waves-Spontaneous activity waves in fire-ant columns

Active matter, which includes crowds of organisms, is composed of constituents that independently consume and dissipate energy. Some active matter systems have been shown to sustain the propagation of various types of waves, resulting from the interplay between density and alignment. Here, we examine a type of solitary wave in dense two-dimensional columns of Solenopsis invicta, fire ants, in which the local activity, density and alignment all play a key role. We demonstrate that these waves are nonlinear and that they are composed of aligned ants that are constrained at the top by the time it takes disordered ants to activate and align and at the bottom by a density minimum enforced by gravity. Our results suggest that intrinsically switchable activity can be a productive framework to understand and trigger a broad range of wave-like behaviors, including stampedes in crowds and herds.

Scale-free chaos in the confined Vicsek flocking model

The Vicsek model encompasses the paradigm of active dry matter. Motivated by collective behavior of insects in swarms, we have studied finite-size effects and criticality in the three-dimensional, harmonically confined Vicsek model. We have discovered a phase transition that exists for appropriate noise and small confinement strength. On the critical line of confinement versus noise, swarms are in a state of scale-free chaos characterized by minimal correlation time, correlation length proportional to swarm size and topological data analysis. The critical line separates dispersed single clusters from confined multicluster swarms. Scale-free chaotic swarms occupy a compact region of space and comprise a recognizable "condensed" nucleus and particles leaving and entering it. Susceptibility, correlation length, dynamic correlation function, and largest Lyapunov exponent obey power laws. The critical line and a narrow criticality region close to it move simultaneously to zero confinement strength for infinitely many particles. At the end of the first chaotic window of confinement, there is another phase transition to infinitely dense clusters of finite size that may be termed flocking black holes.

Flocking in one dimension: Asters and reversals

We study the one-dimensional active Ising model in which aligning particles undergo diffusion biased by the signs of their spins. The phase diagram obtained varying the density of particles, their hopping rate, and the temperature controlling the alignment shows a homogeneous disordered phase but no homogeneous ordered one, as well as two phases with localized dense structures. In the flocking phase, large ordered aggregates move ballistically and stochastically reverse their direction of motion. In what we termed the "aster" phase, dense immobile aggregates of opposite magnetization face each other, exchanging particles, without any net motion of the aggregates. Using a combination of numerical simulations and mean-field theory, we study the evolution of the shapes of the flocks, the statistics of their reversal times, and their coarsening dynamics. Solving exactly for the zero-temperature dynamics of an aster allows us to understand their coarsening, which shows extremal dynamics, while mean-field equations account for their shape.

Intermittent attractive interactions lead to microphase separation in nonmotile active matter

Nonmotile active matter exhibits a wide range of nonequilibrium collective phenomena yet examples are crucially lacking in the literature. We present a microscopic model inspired by the bacteria Neisseria meningitidis in which diffusive agents feel intermittent attractive forces. Through a formal coarse-graining procedure, we show that this truly scalar model of active matter exhibits the time-reversal-symmetry breaking terms defining the Active Model B+ class. In particular, we confirm the presence of microphase separation by solving the kinetic equations numerically. We show that the switching rate controlling the interactions provides a regulation mechanism tuning the typical cluster size, e.g., in populations of bacteria interacting via type IV pili.

Lane and band formation of oppositely driven colloidal particles in two-dimensional ring geometries

We study the segregation phenomena for oppositely driven colloidal particles in two-dimensional ring geometries by means of Brownian dynamics simulations without hydrodynamic interactions. The particles interact via a repulsive Yukawa potential and are confined to a two-dimensional circular channel by hard walls, in which half of the particles are driven clockwise and the other half are driven counterclockwise. In addition to lane formation, which is commonly found in oppositely driven systems, we found band formation along the angular direction in channels with a very large radius. This indicates that a formation of lanes is prevented in the limit of channels with an infinitely large inner radius. The dependency of this segregation has been examined for the two control parameters, the interaction strength between the particles and the width of the circular channel.

Signatures of directed and spontaneous flocking

Collective motion-or flocking-is an emergent phenomena that underlies many biological processes of relevance, from cellular migrations to animal group movements. In this work, we derive scaling relations for the fluctuations of the mean direction of motion and for the static density structure factor (which encodes static density fluctuations) in the presence of a homogeneous, small external field. This allows us to formulate two different and complementary criteria capable of detecting instances of directed motion exclusively from easily measurable dynamical and static signatures of the collective dynamics, without the need to detect correlations with environmental cues. The static one is informative in large enough systems, while the dynamical one requires large observation times to be effective. We believe these criteria may prove useful to detect or confirm the directed nature of collective motion in in vivo experimental observations, which are typically conducted in complex and not fully controlled environments.

Helical flow states in active nematics

We show that confining extensile nematics in three-dimensional (3D) channels leads to the emergence of two self-organized flow states with nonzero helicity. The first is a pair of braided antiparallel streams-this double helix occurs when the activity is moderate, anchoring negligible, and reduced temperature high. The second consists of axially aligned counter-rotating vortices-this grinder train arises between spontaneous axial streaming and the vortex lattice. These two unanticipated helical flow states illustrate the potential of active fluids to break symmetries and form complex but organized spatiotemporal structures in 3D fluidic devices.

Small Obstacle in a Large Polar Flock

We show that arbitrarily large polar flocks are susceptible to the presence of a single small obstacle. In a wide region of parameter space, the obstacle triggers counterpropagating dense bands leading to reversals of the flow. In very large systems, these bands interact, yielding a never-ending chaotic dynamics that constitutes a new disordered phase of the system. While most of these results were obtained using simulations of aligning self-propelled particles, we find similar phenomena at the continuous level, not when considering the basic Toner-Tu hydrodynamic theory, but in simulations of truncations of the relevant Boltzmann equation.

Susceptibility of Polar Flocks to Spatial Anisotropy

We study the effect of spatial anisotropy on polar flocks by investigating active q-state clock models in two dimensions. In contrast to the equilibrium case, we find that any amount of anisotropy is asymptotically relevant, drastically altering the phenomenology from that of the rotationally invariant case. All of the well-known physics of the Vicsek model, from giant density fluctuations to microphase separation, is replaced by that of the active Ising model, with short-range correlations and complete phase separation. These changes appear beyond a length scale that diverges in the q→∞ limit, so that the Vicsek-model phenomenology is observed in finite systems for weak enough anisotropy, i.e., sufficiently high q. We provide a scaling argument which explains why anisotropy has such different effects in the passive and active cases.

A broad perspective to particle-laden fluid interfaces systems: from chemically homogeneous particles to active colloids

Particles adsorbed to fluid interfaces are ubiquitous in industry, nature or life. The wide range of properties arising from the assembly of particles at fluid interface has stimulated an intense research activity on shed light to the most fundamental physico-chemical aspects of these systems. These include the mechanisms driving the equilibration of the interfacial layers, trapping energy, specific inter-particle interactions and the response of the particle-laden interface to mechanical perturbations and flows. The understanding of the physico-chemistry of particle-laden interfaces becomes essential for taking advantage of the particle capacity to stabilize interfaces for the preparation of different dispersed systems (emulsions, foams or colloidosomes) and the fabrication of new reconfigurable interface-dominated devices. This review presents a detailed overview of the physico-chemical aspects that determine the behavior of particles trapped at fluid interfaces. This has been combined with some examples of real and potential applications of these systems in technological and industrial fields. It is expected that this information can provide a general perspective of the topic that can be exploited for researchers and technologist non-specialized in the study of particle-laden interfaces, or for experienced researcher seeking new questions to solve.

Understanding topological defects in fluidized dry active nematics

Dense assemblies of self-propelling rods (SPRs) may exhibit fascinating collective behaviors and anomalous physical properties that are far away from equilibrium. Using large-scale Brownian dynamics simulations, we investigate the dynamics of disclination defects in 2D fluidized swarming motions of dense dry SPRs (i.e., without hydrodynamic effects) that form notable local positional topological structures that are reminiscent of smectic order. We find the deformations of smectic-like rod layers can create unique polar structures that lead to slow translations and rotations of ±1/2-order defects, which are fundamentally different from the fast streaming defect motions observed in wet active matter. We measure and characterize the statistical properties of topological defects and reveal their connections with the coherent structures. Furthermore, we construct a bottom-up active-liquid-crystal model to analyze the instability of polar lanes, which effectively leads to defect formation between interlocked polar lanes and serves as the origin of the large-scale swarming motions.

Quantitative kinetic theory of flocking with three-particle closure

We consider aligning self-propelled particles in two dimensions. Their motion is given by Langevin equations and includes nonadditive N-particle interactions. The qualitative behavior is as for the famous Vicsek model. We develop a kinetic theory of flocking beyond mean field. In particular, we self-consistently take into account the full pair correlation function. We find excellent quantitative agreement of the pair correlations with direct agent-based simulations within the disordered regime. Furthermore we use a closure relation to incorporate spatial correlations of three particles. In that way we achieve good quantitative agreement of the onset of flocking with direct simulations. Compared to mean-field theory, the flocking transition is shifted significantly toward lower noise because directional correlations favor disorder. We compare our theory with a recently developed Landau-kinetic theory.

Kinetic Monte Carlo Algorithms for Active Matter Systems

We study kinetic Monte Carlo (KMC) descriptions of active particles. We show that, when they rely on purely persistent, active steps, their continuous-time limit is ill-defined, leading to the vanishing of trademark behaviors of active matter such as the motility-induced phase separation, ratchet effects, as well as to a diverging mechanical pressure. We then show how, under an appropriate scaling, mixing passive steps with active ones leads to a well-defined continuous-time limit that however differs from standard active dynamics. Finally, we propose new KMC algorithms whose continuous-time limits lead to the dynamics of active Ornstein-Uhlenbeck, active Brownian, and run-and-tumble particles.

Vicsek model by time-interlaced compression: A dynamical computable information density

Collective behavior, both in real biological systems and in theoretical models, often displays a rich combination of different kinds of order. A clear-cut and unique definition of "phase" based on the standard concept of the order parameter may therefore be complicated, and made even trickier by the lack of thermodynamic equilibrium. Compression-based entropies have been proved useful in recent years in describing the different phases of out-of-equilibrium systems. Here, we investigate the performance of a compression-based entropy, namely, the computable information density, within the Vicsek model of collective motion. Our measure is defined through a coarse graining of the particle positions, in which the key role of velocities in the model only enters indirectly through the velocity-density coupling. We discover that such entropy is a valid tool in distinguishing the various noise regimes, including the crossover between an aligned and misaligned phase of the velocities, despite the fact that velocities are not explicitly used. Furthermore, we unveil the role of the time coordinate, through an encoding recipe, where space and time localities are both preserved on the same ground, and find that it enhances the signal, which may be particularly significant when working with partial and/or corrupted data, as is often the case in real biological experiments.

Breakdown of Ergodicity and Self-Averaging in Polar Flocks with Quenched Disorder

We show that spatial quenched disorder affects polar active matter in ways more complex and far reaching than heretofore believed. Using simulations of the 2D Vicsek model subjected to random couplings or a disordered scattering field, we find in particular that ergodicity is lost in the ordered phase, the nature of which we show to depend qualitatively on the type of quenched disorder: for random couplings, it remains long-range ordered, but qualitatively different from the pure (disorderless) case. For random scatterers, polar order varies with system size but we find strong non-self-averaging, with sample-to-sample fluctuations dominating asymptotically, which prevents us from elucidating the asymptotic status of order.

Fluctuation-Induced Phase Separation in Metric and Topological Models of Collective Motion

We study the role of noise on the nature of the transition to collective motion in dry active matter. Starting from field theories that predict a continuous transition at the deterministic level, we show that fluctuations induce a density-dependent shift of the onset of order, which in turn changes the nature of the transition into a phase-separation scenario. Our results apply to a range of systems, including models in which particles interact with their "topological" neighbors that have been believed so far to exhibit a continuous onset of order. Our analytical predictions are confirmed by numerical simulations of fluctuating hydrodynamics and microscopic models.

Phase separation of an active colloidal suspension via quorum-sensing

We present the Brownian dynamics simulation of an active colloidal suspension in two dimensions, where the self-propulsion speed of a colloid is regulated according to the local density sensed by it. The role of concentration-dependent motility in the phase-separation of colloids and their dynamics is investigated in detail. Interestingly, the system phase separates at a very low packing fraction (Φ≈ 0.125) at higher self-propulsion speeds (Pe), into a dense phase coexisting with a homogeneous phase and attains a long-range crystalline order beyond the transition point. The transition point is quantified here from the local density profiles and local and global-bond order parameters. We have shown that the characteristics of the phase diagram are qualitatively akin to the active Brownian particle (ABP) model. Moreover, our investigation reveals that the density-dependent motility amplifies the slow-down of the directed speed, which facilitates phase-separation even at low packing fractions. The effective diffusivity shows a crossover from quadratic rise to a power-law behavior of exponent 3/2 with Pe in the phase-separated regime. Furthermore, we have shown that the effective diffusion decreases exponentially with packing fraction in the phase-separated regime, while it shows a linear decrease in the single phase regime.

Wavelength Selection by Interrupted Coarsening in Reaction-Diffusion Systems

Wavelength selection in reaction-diffusion systems can be understood as a coarsening process that is interrupted by counteracting processes at certain wavelengths. We first show that coarsening in mass-conserving systems is driven by self-amplifying mass transport between neighboring high-density domains. We derive a general coarsening criterion and show that coarsening is generically uninterrupted in two-component systems that conserve mass. The theory is then generalized to study interrupted coarsening and anticoarsening due to weakly broken mass conservation, providing a general path to analyze wavelength selection in pattern formation far from equilibrium.

Collective motion in large deviations of active particles

We analyze collective motion that occurs during rare (large deviation) events in systems of active particles, both numerically and analytically. We discuss the associated dynamical phase transition to collective motion, which occurs when the active work is biased towards larger values, and is associated with alignment of particles' orientations. A finite biasing field is needed to induce spontaneous symmetry breaking, even in large systems. Particle alignment is computed exactly for a system of two particles. For many-particle systems, we analyze the symmetry breaking by an optimal-control representation of the biased dynamics, and we propose a fluctuating hydrodynamic theory that captures the emergence of polar order in the biased state.

Scale-Free Dynamics in Animal Groups and Brain Networks

Collective phenomena fascinate by the emergence of order in systems composed of a myriad of small entities. They are ubiquitous in nature and can be found over a vast range of scales in physical and biological systems. Their key feature is the seemingly effortless emergence of adaptive collective behavior that cannot be trivially explained by the properties of the system's individual components. This perspective focuses on recent insights into the similarities of correlations for two apparently disparate phenomena: flocking in animal groups and neuronal ensemble activity in the brain. We first will summarize findings on the spontaneous organization in bird flocks and macro-scale human brain activity utilizing correlation functions and insights from critical dynamics. We then will discuss recent experimental findings that apply these approaches to the collective response of neurons to visual and motor processing, i.e., to local perturbations of neuronal networks at the meso- and microscale. We show how scale-free correlation functions capture the collective organization of neuronal avalanches in evoked neuronal populations in nonhuman primates and between neurons during visual processing in rodents. These experimental findings suggest that the coherent collective neural activity observed at scales much larger than the length of the direct neuronal interactions is demonstrative of a phase transition and we discuss the experimental support for either discontinuous or continuous phase transitions. We conclude that at or near a phase-transition neuronal information can propagate in the brain with similar efficiency as proposed to occur in the collective adaptive response observed in some animal groups.

Dynamic Self-Organization of Idealized Migrating Cells by Contact Communication

This Letter investigates what forms of cellular dynamic self-organization are caused through intercellular contact communication based on a theoretical model in which migrating cells perform contact following and contact inhibition and attraction of locomotion. Tuning those strengths causes varieties of dynamic patterns. This further includes a novel form of collective migration, snakelike dynamic assembly. Scrutinizing this pattern reveals that cells in this state can accurately respond to an external directional cue but have no spontaneous global polar order.

Dry Active Matter Exhibits a Self-Organized Cross Sea Phase

The Vicsek model of self-propelled particles is known in three different phases: a polar ordered homogeneous phase, also called the Toner-Tu phase, a phase of polar ordered regularly arranged high density bands with surrounding low density regions without polar order, and a homogeneous phase without polar order. Here, we show that the standard Vicsek model has a fourth phase for large system sizes: a polar ordered cross sea phase. We demonstrate that the cross sea phase is not just a superposition of two waves, but it is an independent complex pattern with an inherently selected crossing angle.