Multi-Species Nonreciprocal Active Matter Model

• The paper introduces a minimal framework of multi-species nonreciprocal active matter that uses a constant inter-species phase shift to break action–reaction symmetry.
• It employs Vicsek-type dynamics and hydrodynamic continuum modeling to derive order parameters, phase boundaries, and scaling laws, providing clear analytic tractability.
• The work reveals emergent phenomena such as chiral flocking and vortex-cell species separation, emphasizing the role of permutation symmetry in non-equilibrium collective behavior.

A multi-species nonreciprocal active matter model describes an ensemble of $N$ self-propelled agents, each belonging to one of $Q$ distinguishable species, with interaction rules breaking action–reaction symmetry while retaining full permutation symmetry over species (Potts symmetry). The fundamental microscopic symmetry is that action–reaction is violated solely by a constant phase shift in inter-species velocity alignment, yet all species are otherwise dynamically equivalent, giving rise to a unique set of emergent phenomena including chiral flocking, vortex-cell species separation, and rich phase coexistence (Woo et al., 21 Dec 2025). This paradigm establishes a minimal, analytically tractable framework for nonreciprocal active mixtures with maximal symmetry.

1. Vicsek-Type Microscopic Dynamics with Permutational Nonreciprocity

The system consists of self-propelled particles in two spatial dimensions indexed by species label $s_n \in \{1, ..., Q\}$, with position $\mathbf{r}_n$ and orientation $\theta_n$. The evolution is governed by discrete or continuous-time generalizations of the Vicsek alignment protocol: $\begin{aligned} \theta_n(t+\Delta t) &= \Arg\bigg[\sum_{m \in \mathcal{N}_n} \exp\big(i(\theta_m(t)+\alpha_{nm})\big)\bigg] + \zeta_n(t),\ \mathbf{r}_n(t+\Delta t) &= \mathbf{r}_n(t) + v_0\,\hat{\mathbf{e}}(\theta_n(t))\,\Delta t, \end{aligned}$ where $\mathcal{N}_n$ is the set of neighbors within a radius $r_0$, $v_0$ is the self-propulsion speed, and $\zeta_n$ is angular noise of magnitude $Q$0.

Nonreciprocity is encoded via the phase shift

$Q$1

so intra-species alignment is standard ($Q$2), but all inter-species alignments carry the same constant shift $Q$3. This ensures \emph{permutation invariance (Potts symmetry)}: the dynamics is invariant under species relabeling, $Q$4.

The continuous-time Langevin representation is

$Q$5

with $Q$6 the alignment strength and $Q$7 angular white noise.

2. Hydrodynamic and Boltzmann Continuum Description

At the continuum level, fields are resolved by species. The local angular Fourier modes are introduced: $Q$8 with $Q$9 (density) and $s_n \in \{1, ..., Q\}$0 (polarization field). The Boltzmann moment expansion and Ginzburg–Landau truncation yield, near onset,

$s_n \in \{1, ..., Q\}$1

with the complex coefficients $s_n \in \{1, ..., Q\}$2 depending on $s_n \in \{1, ..., Q\}$3.

A further reduction in the mixed-chiral regime gives an effective $s_n \in \{1, ..., Q\}$4-Langevin model in a co-rotating frame,

$s_n \in \{1, ..., Q\}$5

where $s_n \in \{1, ..., Q\}$6.

3. Spontaneous Symmetry Breaking and Collective Phases

Key order parameters characterize macroscopic states:

• Global polarization: $s_n \in \{1, ..., Q\}$7
• Net chirality: $s_n \in \{1, ..., Q\}$8
• Species (Potts) order: $s_n \in \{1, ..., Q\}$9

Depending on $\mathbf{r}_n$0 and density, the model displays (Woo et al., 21 Dec 2025):

• Chiral–mixed phase: $\mathbf{r}_n$1, $\mathbf{r}_n$2, $\mathbf{r}_n$3, quasi-long-range order (QLRO): all species participate in a synchronized, rotating flock.
• Species–separated (“vortex cell”) phase: $\mathbf{r}_n$4, $\mathbf{r}_n$5, $\mathbf{r}_n$6, Potts symmetry is spontaneously broken, and each species occupies rotating vortex domains.
• Disordered phase: $\mathbf{r}_n$7, $\mathbf{r}_n$8, no global order.
• Coexistence: spatially heterogeneous coexistence of chiral clusters and vortex cells.

Phase boundaries are set by:

• Hopf bifurcation condition: $\mathbf{r}_n$9.
• Species-separation instability: for large $\theta_n$0, the antisymmetric mode destabilizes.

4. Linear Stability and Bifurcation Analysis

Linearizing the hydrodynamic equations near homogeneous states allows the identification of instabilities:

• The Stuart–Landau amplitude equation for the collective order parameter $\theta_n$1,

$\theta_n$2

with Hopf bifurcation at $\theta_n$3.

• Near $\theta_n$4, a two-species antisymmetric fluctuation $\theta_n$5 acquires a positive Lyapunov exponent:

$\theta_n$6

corresponding to lateral repulsion and rotation of antiparallel flocks.

Finite-size scaling of QLRO is quantified by

$\theta_n$7

indicating a BKT-type transition at the edge of the QLRO phase.

5. Spatio-Temporal Patterns and Scaling Laws

Direct simulations reveal:

• Mixed-chiral state: all particles (all species) traverse large-scale (usually counterclockwise) orbits; local clusters exhibit coherent rotation.
• Vortex-cell state: sharply demixed, each cell predominantly a single species, with clockwise circulation.
• Coexistence regime: spontaneous nucleation, growth, and dissolution of both cluster types, exhibiting lane-like and hybrid “bubble” patterns.

Polarization correlations in the QLRO regime decay algebraically: $\theta_n$8 with $\theta_n$9 the BKT threshold.

6. Physical Mechanism: Interaction of Nonreciprocity and Permutational Symmetry

The model features a uniform inter-species phase shift $\begin{aligned} \theta_n(t+\Delta t) &= \Arg\bigg[\sum_{m \in \mathcal{N}_n} \exp\big(i(\theta_m(t)+\alpha_{nm})\big)\bigg] + \zeta_n(t),\ \mathbf{r}_n(t+\Delta t) &= \mathbf{r}_n(t) + v_0\,\hat{\mathbf{e}}(\theta_n(t))\,\Delta t, \end{aligned}$0 as the only source of nonreciprocity, yet treats all species equivalently (full $\begin{aligned} \theta_n(t+\Delta t) &= \Arg\bigg[\sum_{m \in \mathcal{N}_n} \exp\big(i(\theta_m(t)+\alpha_{nm})\big)\bigg] + \zeta_n(t),\ \mathbf{r}_n(t+\Delta t) &= \mathbf{r}_n(t) + v_0\,\hat{\mathbf{e}}(\theta_n(t))\,\Delta t, \end{aligned}$1 symmetry). For small $\begin{aligned} \theta_n(t+\Delta t) &= \Arg\bigg[\sum_{m \in \mathcal{N}_n} \exp\big(i(\theta_m(t)+\alpha_{nm})\big)\bigg] + \zeta_n(t),\ \mathbf{r}_n(t+\Delta t) &= \mathbf{r}_n(t) + v_0\,\hat{\mathbf{e}}(\theta_n(t))\,\Delta t, \end{aligned}$2, global chiral flocking emerges through a Hopf bifurcation, with QLRO enforced by the underlying two-dimensional symmetry.

At large $\begin{aligned} \theta_n(t+\Delta t) &= \Arg\bigg[\sum_{m \in \mathcal{N}_n} \exp\big(i(\theta_m(t)+\alpha_{nm})\big)\bigg] + \zeta_n(t),\ \mathbf{r}_n(t+\Delta t) &= \mathbf{r}_n(t) + v_0\,\hat{\mathbf{e}}(\theta_n(t))\,\Delta t, \end{aligned}$3, inter-species alignment becomes antagonistic (anti-alignment), yielding mutual repulsion between species and breaking Potts symmetry: the minimal mechanism for spontaneous species separation into vortex-cell mosaics.

Unlike generic nonreciprocal models, where interaction matrices are asymmetric and potentially hierarchy-forming or explicit chase–run–evade, this construction yields nonreciprocal but permutation-symmetric macrodynamics, with rich collective outcomes not observed in either symmetric-reciprocal or strongly asymmetric (e.g., predator–prey) settings.

7. General Significance and Links with Broader Active Matter Models

The multi-species nonreciprocal active matter model with Potts symmetry (Woo et al., 21 Dec 2025) is distinct from previously studied classes of nonreciprocal mixtures, such as Cahn-Hilliard-type (scalar) models (Saha et al., 2020), quorum-sensing mixtures (Duan et al., 2023, Duan et al., 2024), and asymmetric Vicsek-type models (Lardet et al., 22 Dec 2025), in which the interaction matrices lack full species-exchange symmetry.

The chiral phases, vortex-cell tiling, and dynamic bubble/lane coexistence observed here depend crucially on the coexistence of uniform nonreciprocity and permutation invariance, and provide a minimal route to symmetry-breaking phenomena with analytically tractable order parameters, phase diagrams, and scaling laws. These results underscore the importance of symmetry constraints, even in driven, non-equilibrium active mixtures, for determining the repertoire of emergent collective states.