Flocking with a q-fold discrete symmetry: Band-to-lane transition in the active Potts model

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.

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.

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.

Disentangling growth and decay of domains during phase ordering

Using Monte Carlo simulations we study phase-ordering dynamics of a multispecies system modeled via the prototype q-state Potts model. In such a multispecies system, we identify a spin state or species as the winner if it has survived as the majority in the final state, otherwise, we mark them as loser. We disentangle the time (t) dependence of the domain length of the winner from losers, rather than monitoring the average domain length obtained by treating all spin states or species alike. The kinetics of domain growth of the winner at a finite temperature in space dimension d=2 reveal that the expected Lifshitz-Cahn-Allen scaling law ∼t^{1/2} can be realized with no early-time corrections, even for system sizes much smaller than what is traditionally used. Up to a certain period, all others species, i.e., the losers, also show a growth that, however, is dependent on the total number of species, and slower than the expected ∼t^{1/2} growth. Afterwards, the domains of the losers start decaying with time, for which our numerical data appear to be consistent with a ∼t^{-2} behavior. We also demonstrate that this approach of looking into the kinetics even provides new insights for the special case of phase ordering at zero temperature, both in d=2 and d=3.

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.

Suppression of discontinuous phase transitions by particle diffusion

We investigate the phase transitions of the q-state Brownian Potts model in two dimensions (2D) comprising Potts spins that diffuse like Brownian particles and interact ferromagnetically with other spins within a fixed distance. With extensive Monte Carlo simulations we find a continuous phase transition from a paramagnetic to a ferromagnetic phase even for q>4. This is in sharp contrast to the existence of a discontinuous phase transition in the equilibrium q-state Potts model in 2D with q>4. We present detailed numerical evidence for a continuous phase transition and argue that diffusion generated dynamical positional disorder suppresses phase coexistence leading to a continuous transition.