Non-Reciprocal Two-Species Vicsek Model

• The non-reciprocal two-species Vicsek model is a binary active-matter system featuring asymmetric interspecies couplings that fundamentally alter flocking behavior.
• It employs unequal coupling strengths in particle alignment, leading to emergent phenomena such as chiral motion, band formation, demixing, and spatiotemporal chaos.
• Hydrodynamic and kinetic theories uncover distinct instability mechanisms, where density fluctuations and asymmetric alignment trigger both long-wavelength and finite-wavelength order disruptions.

Searching arXiv for papers on non-reciprocal and two-species Vicsek models to ground the article. The non-reciprocal two-species Vicsek model denotes a class of binary active-matter models in which two species of self-propelled particles undergo Vicsek-type orientational updates with asymmetric interspecies couplings, so that the influence of species $A$ on species $B$ differs from the influence of $B$ on $A$. In the strict sense used in recent work, non-reciprocity is implemented through unequal cross-couplings such as $J_{AB}\neq J_{BA}$ or $\chi^{12}\neq \chi^{21}$, while self-alignment is typically retained within each species. This class of models has been used to study chiral motion, demixing, large-scale structure formation, run-and-chase dynamics, and the breakdown of homogeneous flocking states (Dutta et al., 20 Apr 2026). Closely related binary Vicsek models with reciprocal antagonistic couplings or generalized multi-species extensions clarify which phenomena require explicit non-reciprocity and which already arise from heterogeneous but symmetric interactions (Chatterjee et al., 2022).

1. Definition and scope of the model family

A canonical continuous-time non-reciprocal two-species Vicsek model in two dimensions assigns each particle $i$ a position $\mathbf r_i$, a heading $\theta_i$, a species label $s_i\in\{A,B\}$, constant speed $B$0, and metric interactions within a range $B$1. One representative formulation is

$B$2

$B$3

with $B$4 and $B$5 generating non-reciprocity (Woo et al., 7 Apr 2026). A discrete-time metric version uses the standard noisy Vicsek update with a local weighted alignment field

$B$6

followed by

$B$7

with $B$8 in dimensionless units (Dutta et al., 20 Apr 2026).

An alternative microscopic implementation uses discrete-time binary Vicsek mixtures with couplings $B$9: $B$0

$B$1

where $B$2 and the interspecies couplings are decomposed into

$B$3

Here $B$4 is the non-reciprocal part, and the regime of interest is weak nonreciprocity, $B$5 (Myin et al., 29 Oct 2025).

The defining distinction is that a genuine non-reciprocal two-species Vicsek model contains two physical species and an asymmetric $B$6, $B$7 coupling structure. This separates it from reciprocal two-species mixtures with $B$8, from single-species higher-order Vicsek models, and from models whose effective asymmetry arises only from local normalization or many-body effects rather than from an explicitly asymmetric coupling matrix (León et al., 22 Dec 2025).

2. Reciprocal, antagonistic, and non-reciprocal binary Vicsek models

The reciprocal baseline is the two-species Vicsek model (TSVM) of two unfriendly species, where particles of species $B$9 and $A$0 align with their own species and anti-align with the other. In that model the local interaction field is

$A$1

with $A$2, so same-species interactions are aligning and different-species interactions are anti-aligning (Chatterjee et al., 2022). The interaction is antagonistic but not non-reciprocal in the strict sense, because both species treat each other symmetrically through the same $A$3 factor.

A useful rewriting in the reciprocal TSVM introduces

$A$4

for which the update field becomes

$A$5

This rewriting clarifies the origin of parallel flocking (PF) and anti-parallel flocking (APF): PF is favored when the species segregate so that the “wrong” species is absent locally, while APF is favored when both species coexist with compatible $A$6-ordering (Chatterjee et al., 2022).

By contrast, strictly non-reciprocal models impose unequal cross-couplings. In one common antisymmetric case,

$A$7

so $A$8 aligns with $A$9 while $J_{AB}\neq J_{BA}$0 anti-aligns with $J_{AB}\neq J_{BA}$1, and the ratio

$J_{AB}\neq J_{BA}$2

measures the strength of non-reciprocal frustration (Dutta et al., 20 Apr 2026). A related weakly non-reciprocal formulation keeps both species mutually aligning or mutually anti-aligning on average, but with slightly different interspecies couplings, so that $J_{AB}\neq J_{BA}$3 fixes the reciprocal part and $J_{AB}\neq J_{BA}$4 the asymmetry (Myin et al., 29 Oct 2025).

This terminological distinction matters because several nearby models are often conflated. The conformity-based higher-order Vicsek model derives effective pairwise and three-body terms from a single-species local-conformity rule,

$J_{AB}\neq J_{BA}$5

leading to

$J_{AB}\neq J_{BA}$6

with $J_{AB}\neq J_{BA}$7 and $J_{AB}\neq J_{BA}$8. Although its bidirectionally ordered state can be interpreted phenomenologically as two oppositely moving groups, it does not define two species or explicit asymmetric interspecies couplings (León et al., 22 Dec 2025). Likewise, generalized reciprocal binary models with $J_{AB}\neq J_{BA}$9 can show rich phase coexistence, but they are not non-reciprocal at the coupling-matrix level (Lardet et al., 22 Mar 2025).

3. Ordered states and collective phases

Reciprocal antagonistic binary Vicsek models exhibit a coexistence phenomenology that already extends the one-species Vicsek picture. In the TSVM, the flocking transition remains liquid-gas–like, with microphase separation into traveling dense bands moving through a gaseous background. Two band families occur: $\chi^{12}\neq \chi^{21}$0-bands composed mainly of $\chi^{12}\neq \chi^{21}$1-particles and $\chi^{12}\neq \chi^{21}$2-bands composed mainly of $\chi^{12}\neq \chi^{21}$3-particles. Within the coexistence region the system supports PF states, in which $\chi^{12}\neq \chi^{21}$4- and $\chi^{12}\neq \chi^{21}$5-bands move in the same direction, and APF states, in which they move in opposite directions (Chatterjee et al., 2022).

These states are characterized by the order parameters

$\chi^{12}\neq \chi^{21}$6

with $\chi^{12}\neq \chi^{21}$7 and $\chi^{12}\neq \chi^{21}$8. In the thermodynamic limit PF has $\chi^{12}\neq \chi^{21}$9 and $i$0, whereas APF has $i$1 and $i$2; APF ordering is reported to be stronger than PF ordering (Chatterjee et al., 2022).

Explicit non-reciprocity broadens the phase repertoire. In weakly non-reciprocal mixtures with $i$3, increasing $i$4 destabilizes the homogeneous flock and produces a single traveling band mostly made of the more strongly aligning species, while the other species forms a more homogeneous liquid background. In the mutually anti-aligning case $i$5, increasing $i$6 destabilizes homogeneous anti-flocking and drives laning, followed at larger asymmetry by polar clusters with chaotic dynamics and little macroscopic order (Myin et al., 29 Oct 2025).

In strongly non-reciprocal models, additional regimes appear. A continuous-time two-species Vicsek model with $i$7 shows a phase diagram containing homogeneous parallel flocking, homogeneous antiparallel flocking, longitudinal antiparallel lane, species-separated parallel-flocking band, and a chaotic phase. Near the antisymmetric line $i$8, the mean-field homogeneous chiral state is replaced by chaos at large scales (Woo et al., 7 Apr 2026). In the discrete-time metric model with $i$9 aligning with $\mathbf r_i$0 and $\mathbf r_i$1 anti-aligning with $\mathbf r_i$2, the observed regimes include chiral, aligned, anti-aligned, independently aligned, and weakly chiral states, with chirality occupying only a restricted window (Dutta et al., 20 Apr 2026).

Closely related but formally distinct non-Vicsek or discrete-symmetry counterparts display analogous behavior. The non-reciprocal two-species active Ising model shows PF and one-species flocking at weak non-reciprocity, a disordered gas at intermediate frustration, and run-and-chase at strong non-reciprocity, with the latter interpreted as a traveling coupled-band state in which $\mathbf r_i$3 chases $\mathbf r_i$4 while $\mathbf r_i$5 avoids $\mathbf r_i$6 (Mangeat et al., 2024).

4. Hydrodynamic and kinetic descriptions

Several complementary coarse-grained descriptions have been developed. In weakly nonreciprocal mixtures, Boltzmann–Ginzburg–Landau coarse-graining introduces for each species a density field $\mathbf r_i$7 and polarization field $\mathbf r_i$8, obeying

$\mathbf r_i$9

and

$\theta_i$0

The coefficient $\theta_i$1 depends on the densities, and this density dependence is central to the instability mechanism (Myin et al., 29 Oct 2025).

Linearization around homogeneous flocking or anti-flocking yields a $\theta_i$2 problem

$\theta_i$3

whose spectrum reveals two distinct instabilities: the usual Vicsek-type instability near onset with $\theta_i$4 as $\theta_i$5, and a new non-reciprocal instability deeper in the ordered phase with $\theta_i$6 as $\theta_i$7. The latter is the signature of a distinct long-wavelength mechanism (Myin et al., 29 Oct 2025).

In the $\theta_i$8 limit the dynamics reduces to three slow fields $\theta_i$9, with reduced matrix

$s_i\in\{A,B\}$0

and longitudinal eigenvalues

$s_i\in\{A,B\}$1

where

$s_i\in\{A,B\}$2

The instability criterion

$s_i\in\{A,B\}$3

requires the density cross-couplings to have opposite signs or sufficiently strong asymmetry, and is absent in reciprocal mixtures (Myin et al., 29 Oct 2025).

For the strongly non-reciprocal chiral problem, a Boltzmann kinetic theory expands the one-particle distributions $s_i\in\{A,B\}$4 in angular Fourier modes $s_i\in\{A,B\}$5, with density $s_i\in\{A,B\}$6 and complex polarization $s_i\in\{A,B\}$7. The hierarchy is truncated at $s_i\in\{A,B\}$8, and stability of the homogeneous chiral limit cycle is analyzed via Floquet exponents $s_i\in\{A,B\}$9. In all studied parameters, $B$00 peaks at a finite wave number $B$01, so the instability of the homogeneous chiral state is finite-wavelength rather than long-wavelength (Woo et al., 7 Apr 2026).

A different but related kinetic framework has been derived for generalized multi-species Vicsek systems with additive sine couplings,

$B$02

using Fourier modes $B$03. In the reciprocal binary case, this theory predicts a finite-wavelength Turing-Hopf instability with $B$04 and $B$05, matching a measured stripe wavelength $B$06. The main binary treatment is reciprocal, but the general framework extends to directed cyclic couplings for $B$07 (Lardet et al., 22 Dec 2025).

5. Instability mechanisms, demixing, and chirality

The principal mechanism identified in weakly non-reciprocal binary flocks is a feedback loop between local density fluctuations, local polarization, and species-dependent alignment strength. Because the local ordering coefficient $B$08 depends on density, a density fluctuation of one species modifies local order; if the couplings are non-reciprocal, the species respond differently, so one can polarize more strongly, advect differently, and amplify the original fluctuation. This mechanism is absent in a single-species flock and also absent in theories that ignore density fluctuations (Myin et al., 29 Oct 2025).

For mutually aligning species, this instability yields a single traveling band containing mostly one species, moving through a more homogeneous liquid of the other species. For mutually anti-aligning species, the same non-reciprocal mechanism destabilizes homogeneous anti-flocking and drives lanes, then polar clusters with large-scale chaotic dynamics. In both cases, species demixing appears without repulsive interactions (Myin et al., 29 Oct 2025).

Chirality constitutes a different non-reciprocal route. In the discrete-time metric model with antisymmetric interspecies couplings, the mean-field non-motile limit defines a phase-difference map

$B$09

with fixed points

$B$10

The branches $B$11 are chiral, with mean-field angular velocity

$B$12

and they are stable for any $B$13 in that non-motile mean-field picture (Dutta et al., 20 Apr 2026).

In the motile finite-range model, however, chirality is not generic. It requires high density, very low motility, small system size relative to the interaction range, and sufficiently mixed local neighborhoods. The local composition indicator

$B$14

is close to $B$15 in the chiral regime, indicating strong local interspecies mixing (Dutta et al., 20 Apr 2026). This suggests that sustained mixing, rather than non-reciprocity alone, is necessary to synchronize the rotational frustration into a coherent chiral state.

A stronger statement is obtained in the continuous-time model: the homogeneous chiral state survives only below a characteristic length scale. The relevant intrinsic scale is the rotation radius

$B$16

and the time-averaged chirality obeys a finite-size scaling form in the variable $B$17. For small $B$18 the chiral state is stable, while for large $B$19 the global chirality vanishes, indicating that the homogeneous chiral state is not asymptotically stable in large two-dimensional systems (Woo et al., 7 Apr 2026).

6. Finite-size effects, switching, and spatiotemporal chaos

Finite-size effects are central throughout the binary Vicsek literature. In the reciprocal TSVM, PF and APF can both exist in the low-density coexistence region and perform stochastic transitions from one to the other. The steady-state distribution $B$20 becomes bimodal in this regime, and the dwell-time distribution has an exponential tail. At fixed aspect ratio $B$21, the PF dwell time scales as

$B$22

for small $B$23, then crosses over to

$B$24

with an exponential scaling not ruled out. The crossover is controlled by the band width $B$25 relative to $B$26, and occurs roughly when $B$27 (Chatterjee et al., 2022).

In non-reciprocal chiral flocks, finite-size control is even more stringent. Small systems can sustain long-lived coherent chirality, but larger systems show decay of the global chirality $B$28 to zero, accompanied by species segregation and inhomogeneity. The decay-time distribution is inconsistent with simple Poisson nucleation of defects, which argues against a metastability picture based on independent droplet nucleation (Woo et al., 7 Apr 2026).

Beyond the finite-wavelength instability, the large-scale state is identified as extensive spatiotemporal chaos. The evidence includes a positive largest Lyapunov exponent, short-ranged velocity-velocity correlations, a broad energy spectrum, and an extensive number of unstable Floquet modes,

$B$29

A master-slave synchronization protocol reveals a finite chaotic length scale: when the uncoupled core radius is smaller than this scale, the slave is entrained by the master, while for larger cores the slave develops independent chaotic dynamics. The crossover is found around

$B$30

(Woo et al., 7 Apr 2026).

The complementary discrete-time metric study reaches a related conclusion from another angle. There, global chiral order weakens as $B$31 increases; for large systems the phase-difference correlation $B$32 decays rapidly, indicating only local chirality. Population imbalance drives transitions from chiral states to porous PF or APF liquids, and motility imbalance induces asynchronous oscillations and, for sufficiently large mismatch,

$B$33

segregation into moving clusters of the faster species within a more dispersed background of slower particles (Dutta et al., 20 Apr 2026).

7. Conceptual boundaries and related generalizations

A recurrent source of confusion is the tendency to classify any binary or bidirectional Vicsek-like ordering as non-reciprocal. Several recent models show why that identification is too broad. The local-conformity higher-order Vicsek model yields pairwise and three-body terms and supports a bidirectionally ordered phase with $B$34 and $B$35, but it contains one species, symmetric metric interactions, no explicit two-population decomposition, and no asymmetric $B$36, $B$37 coupling matrix (León et al., 22 Dec 2025).

Likewise, generalized reciprocal two-species Vicsek models with

$B$38

already display parallel flocking, antiparallel flocking, flocking stripes, antiparallel flocking stripes, nematic stripes, and disordered hyperuniform phases. In that setting the couplings are reciprocal at the matrix level, even though local normalization by neighbor number can make the effective torques non-reciprocal in practice (Lardet et al., 22 Mar 2025). This suggests that some phenomena often associated with non-reciprocal active matter can also emerge from reciprocal competition between alignment and anti-alignment.

Other extensions broaden the field beyond the strict binary case. Metric-free topological mixtures with weak non-reciprocity, where species $B$39 aligns with both $B$40 and $B$41 while species $B$42 does not align with anyone, exhibit reentrant band–polar liquid–band behavior as a function of noise, and require higher angular harmonics for accurate coarse-grained stability analysis (Tang et al., 2024). Non-reciprocal active Brownian mixtures with Vicsek-like torques and steric repulsion show asymmetric clustering in which single-species clusters chase more dilute accumulations of the other species; this is not a pure Vicsek point-particle model, but it demonstrates that broken action-reaction symmetry in orientational couplings alone can rotate the unstable density mode away from pure clustering or pure demixing (Kreienkamp et al., 2024).

At the multi-species level, cyclic directed couplings provide a fully non-reciprocal generalization. One framework uses a coupling graph

$B$43

while another introduces a constant inter-species phase shift $B$44 and full $B$45 Potts symmetry. These models generate chasing stripes, parity-dependent ordering, chirality, species separation, and coexistence phases, and clarify how the binary phenomenology extends to directed interaction networks (Lardet et al., 22 Dec 2025).

Taken together, these results define the non-reciprocal two-species Vicsek model as a minimal but nontrivial active-matter setting in which asymmetry of interspecies alignment reorganizes flocking through density–polarization feedback, finite-wavelength or long-wavelength instabilities, and frustration-induced rotation. The strongest current conclusion is not that non-reciprocity generically produces a single universal ordered state, but that it opens several competing routes—demixing, chiral motion, bands, lanes, clusters, and chaos—whose selection depends sensitively on whether the asymmetry is weak or strong, reciprocal on average or antisymmetric, and compatible or incompatible with sustained local mixing (Myin et al., 29 Oct 2025).