Pith. sign in

REVIEW 3 major objections 5 minor 2 cited by

Emergent complex phases in a discrete flocking model with reciprocal and non-reciprocal interactions

T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read A two-species active Ising model stabilizes a high-density parallel flocking state that continuous-symmetry flocking models do not, and its non-reciprocal version exhibits run-and-chase dynamics.

desk verdict Genuinely new phases in a two-species flocking model, with a hydrodynamic theory that partially overreaches—still worth serious refereeing. read the letter →

arxiv 2412.02501 v2 pith:C2DUZX7J submitted 2024-12-03 cond-mat.stat-mech

classification cond-mat.stat-mech MSC 82B2682C2282C31 PACS 05.65.+b05.20.-y
keywords two-speciesactiveIsingmodelflockingnon-reciprocalinteractionsAshkin-Tellerrun-and-chasedynamicsmicrophaseseparationmotility-inducedinterfacepinninghydrodynamictheory
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper introduces the two-species active Ising model (TSAIM), a discrete-symmetry lattice flocking model with two species that can interact reciprocally (same-species alignment, cross-species anti-alignment), via species interconversion, or non-reciprocally. Its central claim is that the reciprocal model has a stable high-density parallel flocking (HDPF) state, two dense bands moving in the same direction, a phase that the continuous-symmetry two-species Vicsek model and all prior flocking models do not exhibit. Adding species interconversion maps the model to an active Ashkin-Teller system with four spin/species regions and microphase-separated traveling bands; adding non-reciprocal coupling produces a run-and-chase state and, in the non-motile limit, an oscillatory swap state. The paper also shows that the ordered states are metastable to spontaneous droplet nucleation at low diffusivity and that a motility-induced interface pinning transition occurs at low temperature. A refined mean-field hydrodynamic theory reproduces the simulated phase diagrams.

What carries the argument

The central objects are the two species' local magnetizations: the total polarization vs = mA + mB, the difference va = mA − mB, and the species magnetization m = ρA − ρB, together with the total density ρ. The hydrodynamic equations are closed by a refined mean-field Gaussian ansatz in which m and va are independent Gaussian variables with variance linear in the local density (σ²_m = αmρ, σ²_a = αaρ), introducing two free parameters r1 and r2 that are fitted to the simulated non-motile transition densities; this closure is the step that allows phase-separated and HDPF profiles to emerge in the PDE theory. For the non-reciprocal model the same closure is applied to the per-species magnetizations mA and mB, giving rise to the oscillatory instability that underlies the swap state.

What would settle it

In the microscopic TSAIM, measure the conditional distribution of local va and m given local density inside the coexistence and HDPF regimes; if the distribution is non-Gaussian or the variance is not linear in density, or if the r1 and r2 values required to match the binodals change with system size or activity, then the hydrodynamic validation collapses.

Watch

Extended reading notes

Core claim

The reciprocal TSAIM with conserved species has a high-density parallel flocking (HDPF) state in which the two species form dense bands that propagate in the same direction, occupying half the domain with no gas phase; this state is bistable with the liquid anti-parallel flocking state at high density and becomes the preferred attractor in the infinite-size limit starting from disordered initial conditions. The same model with species interconversion is an active generalization of the Ashkin-Teller model, whose order parameters ⟨vs⟩, ⟨m⟩, and ⟨va⟩ give four distinct spin/species phases and, when spin coupling dominates species coupling, microphase-separated parallel bands. The non-reciprocal TSAIM, with JAB = -JBA, exhibits a run-and-chase state at strong non-reciprocity, where A-bands chase B-bands that flee, and in the non-motile limit an oscillatory swap state with limit-cycle magnetizations that survives in two dimensions. All liquid states are metastable at low diffusivity because spontaneously nucleated counter-propagating droplets destroy the high-density bands, producing stripe-like cluster morphologies, and at sufficiently low temperature the interfaces jam via motility-induced interface pinning.

Load-bearing premise

The hydrodynamic theory's phase diagrams rest on treating local magnetizations as Gaussian variables with variance proportional to density, with two free parameters r1 and r2 fitted to the non-motile transition densities of the same simulations; if that closure fails, the theory's confirmation of the phase diagrams does not follow.

Editorial extensions

If this is right

  • If the central claim is correct, the reciprocal TSAIM is the first flocking model to display a stable high-density parallel flocking phase, distinct from the anti-parallel liquid and the phase-separated coexistence region.
  • Species interconversion turns the model into an active Ashkin-Teller system, so binary mixtures with conversion should show four spin/species states and, when spin coupling dominates, microphase-separated traveling bands.
  • Non-reciprocal coupling generically produces run-and-chase bands at strong coupling and an oscillatory swap state in the non-motile limit, with the two-dimensional ordered state destroyed by any nonzero non-reciprocity.
  • At low diffusivity, the ordered states are metastable to spontaneous droplet nucleation, so experimental or numerical systems in that regime should show striped or clustered morphologies instead of homogeneous polar bands.
  • The hydrodynamic equations with the fitted r1 and r2 reproduce the simulated phase diagrams, providing a closed description that can predict phase boundaries at untested parameters.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable extension not in the paper: if HDPF stability relies on the discrete Ising symmetry of the spin, then softening the spin to Q discrete states (as in the active clock model) should destroy the HDPF phase at finite Q; measuring the parallel-flock lifetime as a function of Q would isolate that mechanism.
  • The Gaussian closure parameters r1 and r2 are the only free inputs to the hydrodynamic theory, so directly measuring the conditional distribution of local va and m in simulations would show whether the theory's phase boundaries are genuine predictions or a two-parameter fit.
  • The run-and-chase state is a minimal predator-prey analogue, suggesting that binary bacterial populations engineered with opposing chemotactic responses could exhibit band chasing without explicit attractive forces.
  • The droplet-nucleation metastability implies that any low-diffusivity realization of two-species flocks will appear as a transient or striped pattern rather than a homogeneous polar band, which could be tested in existing active colloid experiments.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper studies a two-species active Ising model (TSAIM) on a lattice, with reciprocal antiferromagnetic interspecies interactions, species interconversion, and non-reciprocal interactions. Through Monte Carlo simulations and a refined mean-field hydrodynamic theory, it reports a high-density parallel flocking (HDPF) state, Ashkin-Teller-like phase structure with species flip, run-and-chase dynamics under non-reciprocal coupling, and metastability due to spontaneous droplet nucleation and motility-induced interface pinning. The central empirical contribution is the identification of these phases in particle simulations, with supporting hydrodynamic calculations.

Significance. If the direct simulation observations hold, the HDPF state and the run-and-chase dynamics are genuinely new phases in discrete-symmetry multi-species active matter, and the paper provides a useful map of the parameter space. Strengths include a transparent model definition, systematic finite-size and Binder-cumulant analysis, and released C++ codes (Ref. 72). The hydrodynamic arguments are less convincing as an independent validation: the parameters r1 and r2 are fitted to the same simulations, and in the non-motile non-reciprocal sector the theory predicts an ordered state that the microscopic simulations do not exhibit. These issues weaken the abstract's headline claim but do not by themselves negate the simulation evidence.

major comments (3)
  1. [II.B and II.C] Section II.B, Eq. (15): the hydrodynamic order-disorder density is rho* = r1/(2*beta1 - 1), and the authors fit r ~ 2.37 to the simulation Binder crossing in Fig. 2(c), with r2 ~ 1.85 fitted similarly from Fig. 5(c) in Section II.D. The phase-separated profiles and binodals in Figs. 3 and 6 are then generated with these same fitted parameters. The match therefore cannot be presented as an independent confirmation of the phase diagrams; it shows only that the Gaussian closure with fitted parameters reproduces the selected transition lines. I recommend replacing 'validates' with language that distinguishes the direct simulation evidence from the parameterized hydrodynamic reproduction.
  2. [II.D] Section II.D: after deriving the species-flip hydrodynamics, the text states 'we will restrain to the special case r1 = r2 = r' but does not give the numerical value of r used in the calculations shown in Fig. 6. The microscopic fits yield r1 ~ 2.37 and r2 ~ 1.85, so the equality r1 = r2 is not supported by the measured variances. Since the theoretical binodals depend on r through Eqs. (16)-(19), please state the chosen r and provide a sensitivity analysis; otherwise the comparison in Fig. 6 is incompletely specified.
  3. [II.E] Section II.E, Eqs. (21) and (24): the hydrodynamic theory of the non-motile NRTSAIM predicts stable ordered homogeneous solutions for beta1 above beta_o, and this ordered region appears in the theoretical state diagrams Figs. 9(k-l). The microscopic simulations, however, show no ordered state for any nonzero JNR (Figs. 8(h)-(l) and the explicit statement in Section II.E). The authors acknowledge this discrepancy, but it contradicts the abstract's claim that the theory 'confirms the phase diagrams.' The abstract and phase diagrams should be modified to present the ordered region as a mean-field prediction not realized in the 2D simulations, and the validation claim should be restricted to the sectors where theory and simulations agree.
minor comments (5)
  1. [II.B] The same symbol gamma_i appears on both sides of the definition gamma_i = gamma_i exp(r_i/2rho); please use distinct notation, for example tilde-gamma_i, to avoid confusion.
  2. [II.C] The text refers to the fitted value 'r ~ 2.37' and later 'r2 ~ 1.85'; it would help to consistently use r1 and r2 throughout and to state explicitly which transition line each parameter controls.
  3. [Abstract and II.F] The stability of the HDPF state is restricted to sufficiently large diffusion (D slightly above 0.15 in the parameters shown in Fig. 10); the abstract should mention this qualification rather than presenting the HDPF state as unconditionally stable.
  4. [II.E and Fig. 7] The run-and-chase state is identified from snapshots and density profiles; a quantitative criterion, such as a band-propagation order parameter or a measured chase distance, would strengthen the phase classification.
  5. [IV.A] The statement that tmax/delta_t ~ 10^5-10^7 Monte Carlo steps is broad; please state the equilibration and steady-state criteria used for each phase diagram.

Circularity Check

2 steps flagged · score 6.0 of 10

The hydrodynamic 'confirmation' is partly circular: the parameters r1 and r2 are fitted to the simulated non-motile transition densities and then re-emerge as the theory's transition-line predictions, so those lines agree by construction.

  1. fitted input called prediction [Section II.C, 'Hydrodynamic theory of the TSAIM without species flip', Eq. (15), and the simulation fit in Fig. 2(c)]
    "The transition line can be fitted by the expression ρ∗ = r/(2β1 − 1), with r ≃ 2.37. ... The transition occurs when 2β1 − r1/ρ0 is larger than 1, i.e. for a density larger than ρ∗ = r1/(2β1 − 1) ... This expression is compatible with the transition line obtained from the simulations of the microscopic model and shown in Fig. 2(c)."

    The hydrodynamic parameter r1, defined in Eq. (11) as r1 = (2β1)^2 α_a, is used as the same value r that was obtained by fitting the microscopic simulations' non-motile transition density. Eq. (15) then yields ρ* = r1/(2β1−1), so the theory's non-motile transition line coincides with the simulation fit by construction. Calling this agreement a validation of the hydrodynamic theory is circular for this part of the phase diagram; the active-state binodals are additional output, but the non-motile 'confirmation' is the fitted input returned.

  2. fitted input called prediction [Section II.D, 'Hydrodynamic theory of the TSAIM with species flip', Eqs. (16) and (19), with the simulation fit in Fig. 5(c)]
    "The transition line can be fitted by the expression ρ∗2(T2) = r2/(2β2 − 1), with r2 ≃ 1.85. ... Here, we will restrain to the special case r1 = r2 = r ... From Eq. (16), we get a second-order phase transition ... for a density larger than ρ∗2 = r/(2β2 − 1) ... This expression is compatible with the transition line obtained from simulations of the microscopic model and shown in Fig. 5(c)."

    The species-flip transition density ρ*_2 is first fitted to the simulations with r2 ≈ 1.85, and then the hydrodynamic equations are restricted to r1 = r2 = r, so the predicted ρ*_2 = r/(2β2−1) reproduces the fitted line identically. The theory's 'compatibility' with Fig. 5(c) is therefore an identity rather than an independent check. This is a second instance of fitted input being presented as hydrodynamic confirmation.

full rationale

The paper's most valuable findings—the HDPF state in the reciprocal TSAIM and the run-and-chase state in the NRTSAIM—are established by direct particle simulations (Figs. 1–4 and 7–8) and would stand even if no hydrodynamic theory existed. Those results are not circular. The circularity is confined to the hydrodynamic-validation claim. In Sec. II.B the refined mean-field closure introduces r1 and r2 as 'two new parameters'; in Secs. II.C–D the non-motile transition densities are first fitted from the same simulations as ρ* = r/(2β1−1) with r ≈ 2.37 and ρ*_2 = r2/(2β2−1) with r2 ≈ 1.85, and the hydrodynamic equations are then solved with r1 = r2 = r. Because Eq. (15) gives ρ* = r1/(2β1−1) and Eq. (16) gives ρ*_2 = r/(2β2−1), the theory's agreement with the simulated transition lines is built in. The active-state binodals and HDPF profiles produced by the hydrodynamic equations are not forced by this fit and have some independent content, so the circularity is partial rather than total. In the non-reciprocal sector the hydrodynamic theory is not circular but is contradicted by the simulations: Eq. (21) and the stability analysis predict a stable ordered homogeneous state for β1 > βo, while the microscopic simulations show no ordered state for any nonzero JNR, as the authors explicitly concede in Sec. II.E. That independent failure, plus the fitted-input status of the reciprocal transition lines, means the abstract's blanket statement that the hydrodynamic theory 'validates our numerical simulations and confirms the phase diagrams' overstates what the theory establishes. No load-bearing self-citation was found: the one-species refined mean-field ansatz is attributed to Solon and Tailleur, and the self-cited TSVM and MIIP papers are used for context or mechanism, not to force the new phases.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The hydrodynamic theory rests on a refined mean-field ansatz that is an approximation, not a consequence of the microscopic rules. Its fluctuation variances are encoded in the free parameters r1 and r2, which are fit to the simulation transition lines. The species-flip rates assume an exponential Arrhenius form with local density normalization, standard in the field. No new physical entities are introduced; the claimed new phases are steady-state configurations of existing degrees of freedom.

free parameters (3)
  • r1 (fluctuation variance parameter for va) = ≈2.37 from ρ*=r/(2β1−1)
    Section II.C: the non-motile transition density is fit to r≈2.37; the same parameter enters the hydrodynamic equations (Eqs. 8-11) and is used to produce the active phase diagrams.
  • r2 (fluctuation variance parameter for m) = ≈1.85 from ρ*_2=r2/(2β2−1)
    Section II.D: fitted to the species order-disorder transition density in the species-flip case; used in the hydrodynamic theory with r1=r2=r.
  • r1 in non-reciprocal hydrodynamic theory = not stated
    Section II.E, Eq. (20): an r1 parameter appears in the non-reciprocal hydrodynamic equations. The text never gives its value, and the theory comparison uses a different density (ρ0=2.5) than the simulations (ρ0=4), suggesting the value is chosen or fitted rather than derived.
assumptions (3)
  • domain assumption Refined mean-field Gaussian ansatz for local order parameter fluctuations (σ²_m=α_m ρ, σ²_a=α_a ρ).
    Section II.B and Supplementary Notes 1-2. Needed to obtain phase-separated solutions in the hydrodynamic theory. Not derived from the microscopic rules; introduces free parameters r1, r2.
  • ad hoc to paper The fluctuation variances α_m and α_a (hence r1 and r2) are constants independent of density and temperature across the whole phase diagram.
    Section II.B defines r1 and r2 as single numbers; the hydrodynamic theory assumes one value applies for all densities and temperatures, an assumption not justified by the microscopic model.
  • domain assumption Spin-flip and species-flip rates have the local Arrhenius form with 1/ρ normalization (Eqs. 2 and 7).
    Standard in active Ising models (Solon-Tailleur); a modeling assumption for how alignment interactions translate into flipping probabilities.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Emergent complex phases in a discrete flocking model with reciprocal and non-reciprocal interactions." pith.science (2026). https://pith.science/paper/C2DUZX7J

@misc{pith2026241202501,
  author       = {Pith},
  title        = {Pith review of: Emergent complex phases in a discrete flocking model with reciprocal and non-reciprocal interactions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/C2DUZX7J}},
  note         = {Machine review of arXiv:2412.02501}
}
read the original abstract

There is growing interest in multi-species active matter systems with reciprocal and non-reciprocal interactions. While such interactions have been explored in continuous symmetry models, less is known about multi-species discrete-symmetry systems. To address this, we study the two-species active Ising model (TSAIM), a discrete counterpart of the two-species Vicsek model. Our investigation explores both inter-species reciprocal and non-reciprocal interactions, along with the possibility of species interconversion. In the reciprocal TSAIM, we observe the emergence of a high-density parallel flocking state, a feature not seen in previous flocking models. With species interconversion, the TSAIM corresponds to an active extension of the Ashkin-Teller model and exhibits rich state diagrams. In the non-reciprocal TSAIM, a run-and-chase dynamics emerge. We also find that the system is metastable due to droplet excitation and exhibits spontaneous motility-induced interface pinning. A hydrodynamic theory validates our numerical simulations and confirms the phase diagrams.

Figures

Figures reproduced from arXiv: 2412.02501 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: (f) exemplifies a highly efficient non-reciprocal con￾figuration, facilitating the B-particles to maintain the maximum distance from the chasing A-particles. The velocity c ∼ 1.96 of the flocking bands [Figs. 7(a–b)] is larger than the self-propulsion velocity of the p…
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p015_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p016_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p018_10.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Spatially patterned phases in a reaction-time-symmetry-broken model of flocking

    cond-mat.soft 2025-05 conditional novelty 7.0 of 10

    A Vicsek-like flocking model with index-ordered time delays exhibits a new phase deep in the ordered state: two counter-propagating density bands with opposite transverse velocities.

  2. Reentrant phase behavior in binary topological flocks with nonreciprocal alignment

    cond-mat.soft 2024-12 conditional novelty 7.0 of 10

    A small fraction of non-aligning dissenters in a Voronoi-neighbor flock produces reentrant traveling bands, including unusual bands that move through an ordered background.

Reference graph

Works this paper leans on

76 extracted references · 71 canonical work pages · cited by 2 Pith papers

  1. [1]

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liv- erpool, J.Prost, M.Rao, andR.A.Simha, Hydrodynamics of soft active matter, Rev. Mod. Phys.85, 1143 (2013)

  2. [2]

    Ramaswamy, The mechanics and statistics of active matter, Annu

    S. Ramaswamy, The mechanics and statistics of active matter, Annu. Rev. Condens. Matter Phys.1, 323 (2010); Active matter, J. Stat. Mech.: Theor. Exp. 054002 (2017)

  3. [3]

    de Magistris and D

    G. de Magistris and D. Marenduzzo,An introduction to the physics of active matter, Physica A418, 65 (2015)

  4. [4]

    Bechinger, R

    C. Bechinger, R. di Leonardo, H. Löwen, C. Reichhardt, G. Volpe, and G. Volpe, Active particles in complex and crowded environments, Rev. Mod. Phys.88, 045006 (2016)

  5. [5]

    M. R. Shaebani, A. Wysocki, R. G. Winkler, G. Gompper, and H. Rieger, Computational models for active matter, Nat. Rev. Phys.2, 181 (2020)

  6. [6]

    M. J. Bowick, N. Fakhri, M. C. Marchetti, and S. Ra- maswamy, Symmetry, thermodynamics, and topology in active matter, Phys. Rev. X12, 010501 (2022)

  7. [7]

    Strogatz,Sync: The Emerging Science of Spontaneous Order (Penguin, London, 2004)

    S. Strogatz,Sync: The Emerging Science of Spontaneous Order (Penguin, London, 2004)

  8. [8]

    Vicsek, A

    T. Vicsek, A. Czirók, E. Ben-Jacob, I. Cohen and O. Shochet, Novel Type of Phase Transition in a System of Self-Driven Particles, Phys. Rev. Lett.75, 1226 (1995)

Show all 76 references
  1. [9]

    Ginelli,The physics of the Vicsek model, Eur

    F. Ginelli,The physics of the Vicsek model, Eur. Phys. J. Special Topics225, 2099 (2016)

  2. [10]

    A. P. Solon, H. Chaté, and J. Tailleur,From phase to microphase separation in flocking models: The essential role of nonequilibrium fluctuations, Phys. Rev. Lett.114, 068101 (2015)

  3. [11]

    A. P. Solon, J. B. Caussin, D. Bartolo, H. Chaté, and J. Tailleur, Pattern formation in flocking models: A hydro- dynamic description, Phys. Rev. E92, 062111 (2015)

  4. [12]

    Toner and Y

    J. Toner and Y. Tu, Long-Range Order in a Two- Dimensional Dynamical XY Model: How Birds Fly To- gether, Phys. Rev. Lett.75, 4326 (1995)

  5. [13]

    Toner and Y

    J. Toner and Y. Tu,Flocks, herds, and schools: A quan- titative theory of flocking, Phys. Rev. E58, 4828 (1998)

  6. [14]

    Toner,Reanalysis of the hydrodynamic theory of fluid, 22 polar-ordered flocks, Phys

    J. Toner,Reanalysis of the hydrodynamic theory of fluid, 22 polar-ordered flocks, Phys. Rev. E86, 031918 (2012)

  7. [16]

    A.P.SolonandJ.Tailleur, Flocking with discrete symme- try: The two-dimensional active Ising model, Phys. Rev. E 92, 042119 (2015)

  8. [17]

    Mangeat, S

    M. Mangeat, S. Chatterjee, R. Paul, and H. Rieger, Flocking with a q-fold discrete symmetry: Band-to-lane transition in the active Potts model, Phys. Rev. E102, 042601 (2020)

  9. [18]

    Chatterjee, M

    S. Chatterjee, M. Mangeat, R. Paul, and H. Rieger, Flocking and reorientation transition in the 4-state active Potts model, EPL130, 66001 (2020)

  10. [19]

    Karmakar, S

    M. Karmakar, S. Chatterjee, M. Mangeat, H. Rieger, and R. Paul,Jamming and flocking in the restricted active Potts model, Phys. Rev. E108, 014604 (2023)

  11. [20]

    Chatterjee, M

    S. Chatterjee, M. Mangeat, and H. Rieger,Polar flocks with discretized directions: The active clock model ap- proaching the Vicsek model, EPL138, 41001 (2022)

  12. [21]

    Solon, H

    A. Solon, H. Chaté, J. Toner, and J. Tailleur,Suscepti- bility of polar flocks to spatial anisotropy, Phys. Rev. Lett. 128, 208004 (2022)

  13. [22]

    Codina, B

    J. Codina, B. Mahault, H. Chaté, J. Dobnikar, I. Pago- nabarraga, and X.-q. Shi,Small obstacle in a large polar flock, Phys. Rev. Lett.128, 218001 (2022)

  14. [23]

    Karmakar, S

    M. Karmakar, S. Chatterjee, R. Paul and H. Rieger,Con- sequence of anisotropy on flocking: the discretized Vicsek model, New J. Phys.26, 043023 (2024)

  15. [24]

    Besse, H

    M. Besse, H. Chaté, and A. Solon, Metastability of Constant-Density Flocks, Phys. Rev. Lett. 129, 268003 (2022)

  16. [25]

    Benvegnen, O

    B. Benvegnen, O. Granek, S. Ro, R. Yaacoby, H. Chaté, Y. Kafri, D. Mukamel, A. Solon, and J. Tailleur,Metasta- bility of discrete-symmetry flocks, Phys. Rev. Lett.131, 218301 (2023)

  17. [27]

    G. R. Graves and N.J. Gotelli,Assembly of avian mixed- species flocks in Amazonia, PNAS90, 1388, (1993)

  18. [28]

    A. J. W. Ward, T. M. Schaerf, A. L. J. Burns, J. T. Lizier, E. Crosato, M. Prokopenko, and M. M. Webster, Cohesion, order and information flow in the collective mo- tion of mixed-species shoals, R. Soc. open sci.5, 181132 (2018)

  19. [29]

    Ariel, O

    G. Ariel, O. Rimer, and E. Ben-Jacob,Order–disorder phase transition in heterogeneous populations of self- propelled particles, J. Stat. Phys.158, 579 (2015)

  20. [30]

    Peled, S

    S. Peled, S. D. Ryan, S. Heidenreich, M. Bär, G. Ariel, and A. Be’er,Heterogeneous bacterial swarms with mixed lengths, Phys. Rev. E103, 032413 (2021)

  21. [31]

    Zuo and Y

    W. Zuo and Y. Wu,Dynamic motility selection drives population segregation in a bacterial swarm, PNAS 117, 4693 (2020)

  22. [32]

    A. Jose, G. Ariel, and A. Be‘er, Physical characteris- tics of mixed-species swarming colonies, Phys. Rev. E105, 064404 (2022)

  23. [33]

    A. M. Menzel,Collective motion of binary self-propelled particle mixtures, Phys. Rev. E85, 021912 (2012)

  24. [34]

    S. R. McCandlish, A. Baskaran, and M. F. Hagan,Spon- taneous segregation of self-propelled particles with different motilities, Soft Matter8, 2527 (2012)

  25. [35]

    Mishra, K

    S. Mishra, K. Tunstrøm, I. D. Couzin, and C. Huepe, Collective dynamics of self-propelled particles with variable speed, Phys. Rev. E86, 011901 (2012)

  26. [36]

    Netzer, Y

    G. Netzer, Y. Yarom, and G. Ariel,Heterogeneous pop- ulations in a network model of collective motion, Physica A 530, 121550 (2019)

  27. [37]

    G. Book, C. Ingham, and G. Ariel,Modeling cooperating micro-organisms in antibiotic environment, Plos One12, e0190037 (2017)

  28. [38]

    Copenhagen, D

    K. Copenhagen, D. A. Quint, and A. Gopinathan,Self- organized sorting limits behavioral variability in swarms, Sci. Rep.6, 31808 (2016)

  29. [39]

    Khodygo, M

    V. Khodygo, M. T. Swain, and A. Mughal, Homoge- neous and heterogeneous populations of active rods in two- dimensional channels, Phys. Rev. E99, 022602 (2019)

  30. [40]

    P. K. Bera and A. K. Sood, Motile dissenters disrupt the flocking of active granular matter, Phys. Rev. E101, 052615 (2020)

  31. [41]

    N. A. Mecholsky, E. Ott, and T. M. Antonsen Jr.,Obsta- cle and predator avoidance in a model for flocking, Physica D 239, 988 (2010)

  32. [42]

    Sengupta, T

    A. Sengupta, T. Kruppa, and H. Löwen, Chemotactic predator-prey dynamics, Phys. Rev. E83, 031914 (2011)

  33. [44]

    Fruchart, R

    M. Fruchart, R. Hanai, P. B. Littlewood, and V. Vitelli, Non-reciprocal phase transitions, Nature592, 363 (2021)

  34. [45]

    D.Martin, D.Seara, Y.Avni, M.Fruchart, andV.Vitelli, The transition to collective motion in nonreciprocal active matter: coarse graining agent-based models into fluctuat- ing hydrodynamics, arXiv:2307.08251 (2024)

  35. [46]

    A. V. Ivlev, J. Bartnick, M. Heinen, C.-R. Du, V. Nosenko, and H. Löwen,Statistical mechanics where New- ton’s third law is broken, Phys. Rev. X5, 011035 (2015)

  36. [47]

    Z. You, A. Baskaran, and M. C. Marchetti,Nonreciproc- ity as a generic route to traveling states, PNAS117, 19767 (2020)

  37. [48]

    S. Saha, J. Agudo-Canalejo, and R. Golestanian,Scalar active mixtures: The nonreciprocal Cahn-Hilliard model, Phys. Rev. X10, 041009 (2020)

  38. [49]

    J. Chen, X. Lei, Y. Xiang, M. Duan, X. Peng, and H. P. Zhang,Emergent Chirality and Hyperuniformity in an Active Mixture with Nonreciprocal Interactions, Phys.Rev. Lett. 132, 118301 (2024)

  39. [50]

    K. L. Kreienkamp and S. H. L. Klapp,Dynamical struc- tures in phase-separating nonreciprocal polar active mix- tures, Phys. Rev. E110, 064135 (2024)

  40. [51]

    K. L. Kreienkamp and S. H. L. Klapp, Nonreciprocal Alignment Induces Asymmetric Clustering in Active Mix- tures, Phys. Rev. Lett.133, 258303 (2024)

  41. [52]

    Hong and S

    H. Hong and S. H. Strogatz,Kuramoto model of coupled oscillators with positive and negative coupling parameters: an example of conformist and contrarian oscillators, Phys. Rev. Lett.106, 054102 (2011)

  42. [53]

    Golomb and D

    D. Golomb and D. Hansel,The number of synaptic in- puts and the synchrony of large, sparse neuronal networks, Neural Comput.12, 1095 (2000)

  43. [54]

    Börgers and N

    C. Börgers and N. Kopell,Synchronization in networks of excitatory and inhibitory neurons with sparse, random connectivity, Neural Comput.15, 509 (2003)

  44. [55]

    Xiong, Y

    L. Xiong, Y. Cao, R. Cooper, W. J. Rappel, J. Hasty and L. Tsimring,Flower-like patterns in multi-species bacterial colonies, eLife9, e48885 (2020)

  45. [56]

    Theveneau, B

    E. Theveneau, B. Steventon, E. Scarpa, S. Garcia, X. 23 Trepat, A. Streit and R. Mayor,Chase-and-run between adjacent cell populations promotes directional collective migration, Nat. Cell Biol.15, 763 (2013)

  46. [57]

    C. H. Meredith, P. G. Moerman, J. Groenewold, Yu-Jen Chiu, W. K. Kegel, A. van Blaaderen and L. D. Zarzar, Predator–prey interactions between droplets driven by non- reciprocal oil exchange, Nat. Chem.12, 1136 (2020)

  47. [58]

    S. A.M. Loos andS. H.L. Klapp,Irreversibility, heat and information flows induced by non-reciprocal interactions, New J. Phys.22, 123051 (2020)

  48. [59]

    Durve and A

    M. Durve and A. Sayeed, First-order phase transition in a model of self-propelled particles with variable angular range of interaction, Phys. Rev. E93, 052115 (2016)

  49. [60]

    Durve, A

    M. Durve, A. Saha, and A. Sayeed,Active particle con- densation by non-reciprocal and time-delayed interactions, Eur. Phys. J. E41, 49 (2018)

  50. [61]

    Ashkin and E

    J. Ashkin and E. Teller,Statistics of Two-Dimensional Lattices with Four Components, Phys.Rev.64, 178(1943)

  51. [62]

    Fan,On critical properties of the Ashkin-Teller model, Physics Letters A39, 136 (1972)

    C. Fan,On critical properties of the Ashkin-Teller model, Physics Letters A39, 136 (1972)

  52. [63]

    Clifford and A

    P. Clifford and A. Sudbury,A model for spatial conflict, Biometrika60, 581 (1973)

  53. [64]

    R. A. Holley and T. M. Liggett, Ergodic theorems for weakly interacting infinite systems and the voter model, Ann. Probab.3, 643 (1975)

  54. [65]

    R. V. Ditzian, J. R. Banavar, G. S. Grest, and L. P. Kadanoff, Phase diagram for the Ashkin-Teller model in three dimensions, Phys. Rev. B22, 2542 (1980)

  55. [66]

    F. A. da Costa, M. J. de Oliveira, and S. R. Salinas, Asymmetric Baxter phase in the symmetric Ashkin-Teller model, Phys. Rev. B36, 7163 (1987)

  56. [68]

    Binder,Applications of Monte Carlo methods to sta- tistical physics, Rep

    K. Binder,Applications of Monte Carlo methods to sta- tistical physics, Rep. Prog. Phys.60, 487 (1997)

  57. [69]

    S. K. Baek, P. Minnhagen, and B. J. Kim, True and quasi-long-range order in the generalized q-state clock model, Phys. Rev. E80, 060101(R) (2009)

  58. [70]

    Y.Avni, M.Fruchart, D.Martin, D.Seara, andV.Vitelli, Nonreciprocal Ising model, Phys. Rev. Lett.134, 117103 (2025)

  59. [71]

    S. Basu, S. Majumder, S. Sutradhar, S. K. Das and R. Paul, Phase segregation in a binary fluid confined inside a nanopore, EPL116, 56003 (2016)

  60. [72]

    Emergent complex phases in a discrete flock- ing model with reciprocal and non-reciprocal interac- tions

    M. Mangeat, S. Chatterjee, J. D. Noh, and H. Rieger, Codes for “Emergent complex phases in a discrete flock- ing model with reciprocal and non-reciprocal interac- tions”, Zenodo, https://doi.org/10.5281/zenodo.14257074 (2024)

  61. [73]

    W. H. Press, Numerical recipes: The art of scientific computing (Cambridge University Press, Cambridge, UK, 2007)

  62. [74]

    Emergent complex phases in a discrete flocking model with reciprocal and non-reciprocal interactions

    R. Courant, K. Friedrichs, and H. Lewy, Über die partiellen Differenzengleichungen der mathematischen Physik, Math. Ann.100, 32–74 (1928). 24 Supplemental Material for “Emergent complex phases in a discrete flocking model with reciprocal and non-reciprocal interactions” Matthi...

  63. [75]

    A. P. Solon and J. Tailleur,Revisiting the Flocking Tran- sition Using Active Spins, Phys. Rev. Lett.111, 078101 (2013)

  64. [76]

    A. P. Solon and J. Tailleur,Flocking with discrete symme- try: The two-dimensional active Ising model, Phys. Rev. E 92, 042119 (2015)

  65. [77]

    Chatterjee, M

    S. Chatterjee, M. Mangeat, C. U. Woo, H. Rieger, and J. D. Noh, Flocking of two unfriendly species: The two- species Vicsek model, Phys. Rev. E107, 024607 (2023)

  66. [78]

    Emergent com- plex phases in a discrete flocking model with re- ciprocal and non-reciprocal interactions

    M. Mangeat, S. Chatterjee, J. D. Noh, and H. Rieger, Supplementary movies for “Emergent com- plex phases in a discrete flocking model with re- ciprocal and non-reciprocal interactions” , Zenodo, https://doi.org/10.5281/zenodo.14261008 (2024)

  67. [79]

    C. U. Woo and J. D. Noh,Motility-Induced Pinning in Flocking System with Discrete Symmetry, Phys. Rev. Lett. 133, 188301 (2024)

  68. [80]

    A. P. Solon, H. Chaté, and J. Tailleur, From phase to microphase separation in flocking models: The essential role of nonequilibrium fluctuations, Phys. Rev. Lett.114, 068101 (2015)

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.