Contrarian compulsions produce exotic time-dependent flocking of active particles

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.

Pattern formation by bacteria-phage interactions

The interactions between bacteria and phages-viruses that infect bacteria-play critical roles in agriculture, ecology, and medicine; however, how these interactions influence the spatial organization of both bacteria and phages remain largely unexplored. Here, we address this gap in knowledge by developing a theoretical model of motile, proliferating bacteria that aggregate via motility-induced phase separation (MIPS) and encounter phage that infect and lyse the cells. We find that the non-reciprocal predator-prey interactions between phage and bacteria strongly alter spatial organization, in some cases giving rise to a rich array of finite-scale stationary and dynamic patterns in which bacteria and phage coexist. We establish principles describing the onset and characteristics of these diverse behaviors, thereby helping to provide a biophysical basis for understanding pattern formation in bacteria-phage systems, as well as in a broader range of active and living systems with similar predator-prey or other non-reciprocal interactions.

Emergence of chimeras: An impetus by exceptional points

A nonconservative system with nonreciprocal interaction has been found to reveal exotic features where sudden phase transitions can occur. In this paper, the emanation of a chimera in a network of Stuart-Landau oscillators with nonreciprocal interaction is reported. Note that the spins follow the random discrete distribution. In other words, we pick a random oscillator to rotate clockwise or anticlockwise. At the transition points, we find the spectral singularities in the eigenplane, where eigenvalues coalesce, commonly known as exceptional points. We find that the counterrotational symmetry breaking induced by exceptional points antecedents the occurrence of the chimera. We numerically attest to the findings for two cases of initial conditions, namely, bipartite and random. We also extend our study to a two-dimensional array of nonreciprocally interacting, distributed spins. The findings could have pragmatic implications in the areas of active matter, networks, and photonics.

Curvature-controlled geometrical lensing behavior in self-propelled colloidal particle systems

In many biological systems, the curvature of the surfaces cells live on influences their collective properties. Curvature should likewise influence the behavior of active colloidal particles. We show using molecular simulation of self-propelled active particles on surfaces of Gaussian curvature (both positive and negative) how curvature sign and magnitude can alter the system's collective behavior. Curvature acts as a geometrical lens and shifts the critical density of motility-induced phase separation (MIPS) to lower values for positive curvature and higher values for negative curvature, which we explain theoretically by the nature of parallel lines in spherical and hyperbolic space. Curvature also fluidizes dense MIPS clusters due to the emergence of defect patterns disrupting the crystalline order inside the clusters. Using our findings, we engineer three confining surfaces that strategically combine regions of different curvature to produce a host of novel dynamical behaviors, including cyclic MIPS on spherocylinders, directionally biased cyclic MIPS on spherocones, and position dependent cluster fluctuations on metaballs.

Alignment-mediated segregation in an active-passive mixture

We report segregation between the athermal active and passive particles mediated by the local alignment interaction in a confined space. The competition between the alignment interaction and self-propulsion force results in a transition between disordered and ordered phases. We show that as the coordination between the particles increases, they form an ordered mill, which helps the particles to aggregate into isotropic clusters. As a result, particles segregate into active core and passive shells. This segregation phenomenon is adversely affected by the packing fraction and the size dispersion between active and passive particles. We show that this adverse effect can be overcome by incorporating higher coordination in the system. We report that the monodispersed system is more desirable for segregation in a binary mixture than a bidispersed system, as the latter favors the mixed state.

Linear response in the uniformly heated granular gas

We analyze the linear response properties of the uniformly heated granular gas. The intensity of the stochastic driving fixes the value of the granular temperature in the nonequilibrium steady state reached by the system. Here, we investigate two specific situations. First, we look into the "direct" relaxation of the system after a single (small) jump of the driving intensity. This study is carried out by two different methods. Not only do we linearize the evolution equations around the steady state, but we also derive generalized out-of-equilibrium fluctuation-dissipation relations for the relevant response functions. Second, we investigate the behavior of the system in a more complex experiment, specifically a Kovacs-like protocol with two jumps in the driving. The emergence of an anomalous Kovacs response is explained in terms of the properties of the direct relaxation function: it is the second mode changing sign at the critical value of the inelasticity that demarcates anomalous from normal behavior. The analytical results are compared with numerical simulations of the kinetic equation, and a good agreement is found.

Non-reciprocal phase transitions

Out of equilibrium, a lack of reciprocity is the rule rather than the exception. Non-reciprocity occurs, for instance, in active matter^1-6, non-equilibrium systems^7-9, networks of neurons^10,11, social groups with conformist and contrarian members¹², directional interface growth phenomena^13-15 and metamaterials^16-20. Although wave propagation in non-reciprocal media has recently been closely studied^1,16-20, less is known about the consequences of non-reciprocity on the collective behaviour of many-body systems. Here we show that non-reciprocity leads to time-dependent phases in which spontaneously broken continuous symmetries are dynamically restored. We illustrate this mechanism with simple robotic demonstrations. The resulting phase transitions are controlled by spectral singularities called exceptional points²¹. We describe the emergence of these phases using insights from bifurcation theory^22,23 and non-Hermitian quantum mechanics^24,25. Our approach captures non-reciprocal generalizations of three archetypal classes of self-organization out of equilibrium: synchronization, flocking and pattern formation. Collective phenomena in these systems range from active time-(quasi)crystals to exceptional-point-enforced pattern formation and hysteresis. Our work lays the foundation for a general theory of critical phenomena in systems whose dynamics is not governed by an optimization principle.

Tracking collective cell motion by topological data analysis

By modifying and calibrating an active vertex model to experiments, we have simulated numerically a confluent cellular monolayer spreading on an empty space and the collision of two monolayers of different cells in an antagonistic migration assay. Cells are subject to inertial forces and to active forces that try to align their velocities with those of neighboring ones. In agreement with experiments in the literature, the spreading test exhibits formation of fingers in the moving interfaces, there appear swirls in the velocity field, and the polar order parameter and the correlation and swirl lengths increase with time. Numerical simulations show that cells inside the tissue have smaller area than those at the interface, which has been observed in recent experiments. In the antagonistic migration assay, a population of fluidlike Ras cells invades a population of wild type solidlike cells having shape parameters above and below the geometric critical value, respectively. Cell mixing or segregation depends on the junction tensions between different cells. We reproduce the experimentally observed antagonistic migration assays by assuming that a fraction of cells favor mixing, the others segregation, and that these cells are randomly distributed in space. To characterize and compare the structure of interfaces between cell types or of interfaces of spreading cellular monolayers in an automatic manner, we apply topological data analysis to experimental data and to results of our numerical simulations. We use time series of data generated by numerical simulations to automatically group, track and classify the advancing interfaces of cellular aggregates by means of bottleneck or Wasserstein distances of persistent homologies. These techniques of topological data analysis are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant for tissue engineering, biological effects of materials, tissue and organ regeneration.

Spatiotemporal dynamics of a self-propelled system with opposing alignment and repulsive forces

Effect of concurrent alignment and repulsion is studied in the purview of a confined active matter system using a modified force-based Vicsek model. On alteration of the alignment and the repulsive force parameters, a low alignment random phase, a midrange alignment milling phase, and a high alignment oscillatory phase are identified. Based on the particle aggregations, the milling phase is further classified into three subphases, two of which are spatial patterns: one consisting of compact ring-shaped mills and the other incorporating both rings and clusters. A correlation function based on the inner product of spatial velocity fluctuations of the particles shows a high correlation length for the ringed milling and the rings-clusters hybrid milling state. On analyzing temporal velocity fluctuations of particles through chaos detection techniques, low alignment and high alignment states are indicative of chaos, while the middle order alignment is symbolic of periodicity. The extent of synchronization of the particles' motion is analyzed through a Hilbert transform-based mean frequency approach, leading to the detection of a weak chimera state in the case of the spatial structures. The ringed milling state shows a unique category of weak chimera consisting of multiple oscillator groups showcasing different synchronization frequencies coexisting with desynchronized oscillators.

Nonreciprocity as a generic route to traveling states

We examine a nonreciprocally 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 nonreciprocity, 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 nonreciprocity, the static demixed pattern acquires a finite velocity, resulting in a state that breaks both spatial and time-reversal 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 nonreciprocal 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 nonconservative dynamics, from reaction-diffusion and prey-predators models to multispecies mixtures of microorganisms with antagonistic interactions.

Analysis of aligning active local searchers orbiting around their common home position

We discuss effects of pairwise aligning interactions in an ensemble of central place foragers or of searchers that are connected to a common home. In a wider sense, we also consider self-moving entities that are attracted to a central place such as, for instance, the zooplankton Daphnia being attracted to a beam of light. Single foragers move with constant speed due to some propulsive mechanism. They explore at random loops the space around and return rhytmically to their home. In the ensemble, the direction of the velocity of a searcher is aligned to the motion of its neighbors. At first, we perform simulations of this ensemble and find a cooperative behavior of the entities. Above an overcritical interaction strength the trajectories of the searcher qualitatively changes and searchers start to move along circles around the home position. Thereby, all searchers rotate either clockwise or anticlockwise around the central home position as it was reported for the zooplankton Daphnia. At second, the computational findings are analytically explained by the formulation of transport equations outgoing from the nonlinear mean field Fokker-Planck equation of the considered situation. In the asymptotic stationary limit, we find expressions for the critical interaction strength, the mean radial and orbital velocities of the searchers and their velocity variances. We also obtain the marginal spatial and angular densities in the undercritical regime where the foragers behave like individuals as well as in the overcritical regime where they rotate collectively around the considered home. We additionally elaborate the overdamped Smoluchowski-limit for the ensemble.

Active Ornstein-Uhlenbeck particles

Active Ornstein-Uhlenbeck particles (AOUPs) are overdamped particles in an interaction potential subject to external Ornstein-Uhlenbeck noises. They can be transformed into a system of underdamped particles under additional velocity dependent forces and subject to white noise forces. There has been some discussion in the literature on whether AOUPs can be in equilibrium for particular interaction potentials and how far from equilibrium they are in the limit of small persistence time. By using a theorem on the time reversed form of the AOUP Langevin-Ito equations, I prove that they have an equilibrium probability density invariant under time reversal if and only if their smooth interaction potential has zero third derivatives. In the limit of small persistence Ornstein-Uhlenbeck time τ, a Chapman-Enskog expansion of the Fokker-Planck equation shows that the probability density has a local equilibrium solution in the particle momenta modulated by a reduced probability density that varies slowly with the position. The reduced probability density satisfies a continuity equation in which the probability current has an asymptotic expansion in powers of τ. Keeping up to O(τ) terms, this equation is a diffusion equation, which has an equilibrium stationary solution with zero current. However, O(τ^{2}) terms contain fifth- and sixth-order spatial derivatives and the continuity equation no longer has a zero current stationary solution. The expansion of the overall stationary solution now contains odd terms in the momenta, which clearly shows that it is not an equilibrium.