pith. sign in

arxiv: 2507.17108 · v2 · submitted 2025-07-23 · ❄️ cond-mat.soft · cond-mat.stat-mech

Defect-Mediated Aggregation and Motility-Induced Phase Separation in Self-Propelled Lattice-Gas Active XY Model

Shun Inoue , Satoshi Yukawa This is my paper

Pith reviewed 2026-05-19 03:35 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords active XY modelmotility-induced phase separationtopological defectsvortex chargeself-propelled particleslattice gasphase separation kineticsactive matter
0
0 comments X

The pith

Self-propulsion in the active XY model causes particles to aggregate around positive vortex defects, producing motility-induced phase separation whose relaxation time scales as the cube of system size.

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

The authors introduce a lattice-based active XY model that adds a self-propulsion bias to the usual spin interactions of the classical XY model. Numerical runs show that this bias drives particles to collect into dense clusters centered on topological defects that carry positive vortex charge, while defects with negative charge simply disappear. The clusters themselves enlarge through a two-stage exponential process whose characteristic time grows as the cube of the linear system size, echoing the scaling seen in equilibrium first-order phase transitions. The work therefore positions topological defects as the organizing agents that turn uniform active motion into macroscopic phase separation.

Core claim

In the proposed active XY model, self-propulsion induces motility-induced phase separation in which particles form clusters that nucleate and grow around topological defects possessing positive vortex charge. Defects carrying negative charge instead dissipate without forming stable aggregates. The temporal evolution of cluster size proceeds via two distinct exponential relaxation stages, and the associated characteristic time scales as τ ∼ L³ with system size L, a dependence that mirrors the kinetics of first-order phase separation in equilibrium systems.

What carries the argument

The self-propulsion parameter that imposes a directional bias on particle hops, acting together with the lattice-gas update rules to make positive-vortex-charge defects stable sites for particle accumulation.

If this is right

  • Positive vortex defects act as persistent nucleation centers for dense clusters while negative defects are transient and do not support aggregation.
  • Cluster growth occurs in two successive exponential stages whose duration is controlled by system size.
  • The L³ scaling of the relaxation time places the nonequilibrium clustering process in the same dynamical class as equilibrium first-order phase separation.
  • Topological defects therefore function as the microscopic agents that convert uniform self-propulsion into macroscopic phase separation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the same defect selectivity survives in off-lattice or hydrodynamic versions of the model, defect charge could serve as a general predictor of clustering locations across many active-matter systems.
  • The observed L³ scaling suggests that defect-mediated MIPS may exhibit finite-size effects that are experimentally measurable in microfluidic or colloidal active-matter setups of varying container size.
  • One could test whether adding explicit hydrodynamic interactions alters the two-stage relaxation or merely rescales the prefactor in the L³ law.

Load-bearing premise

The discrete lattice rules and the chosen form of self-propulsion bias are sufficient to capture the essential defect dynamics and clustering behavior of real active matter.

What would settle it

A set of larger-system simulations in which the measured relaxation time fails to grow proportionally to L³ or in which negative-charge defects remain stable and nucleate clusters would contradict the reported scaling and defect selectivity.

Figures

Figures reproduced from arXiv: 2507.17108 by Satoshi Yukawa, Shun Inoue.

Figure 1
Figure 1. Figure 1: Schematic picture of particle movement rates and its orientation. The movement rate depends on the self-propulsion parameter ϵ and the spin orientation. Particles at each lattice site can move probabilistically in four directions: up, down, left, and right. The XY spin orientation is defined by an angle θ as in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Snapshots of the simulation for various values of ρ and ϵ. The density increases to the right, and the self-propulsion parameter increases downward in this grid layout. The hue color circle represents the particle orientation, corresponding to the color wheel and orientation shown in the upper left. The white color represents the empty sites. Increasing the self-propulsion parameter ϵ induces motility-indu… view at source ↗
Figure 3
Figure 3. Figure 3: Emerging topological defects. Due to the self-propulsion and the particle exclusion, the m = +1 topological defect shown on the leftmost figure tends to persist intuitively, whereas all instances of the m = −1 topological defects shown on the right two figures dissipate. Even among the m = +1 topological defects, the topological defect shown on the second left is less likely to serve as a nucleus for clust… view at source ↗
Figure 4
Figure 4. Figure 4: Cluster merging process. The system size is taken to be L = 240. ± signs represent topological defects with the vortex charge m = ±1. As time increases from the left figure to the right one, we observe that two clusters with m = +1 vortex charge merge in a rectangular region shown by a dashed line. In the center, the two clusters are just touching. Then, a topological defect with m = −1 charge emerges at t… view at source ↗
Figure 5
Figure 5. Figure 5: Relationship between the self-propulsion parameter ϵ and the total vorticity N at ρ = 0.5 and L = 60. The total vorticity is calculated every 2 500 MCSs after discarding the initial 500 MCSs for initial relaxation. The initial configuration is random, and the average vorticity is near zero. The horizontal red line indicates N = 0. The figure shows an average of over 100 samples. For ϵ = 0.0, corresponding … view at source ↗
Figure 6
Figure 6. Figure 6: Relationship between the density ρ and the total vorticity N with the self-propulsion ϵ = 1 and L = 60. The total vorticity is calculated every 2 500 MCSs after discarding the initial 500 MCSs for initial relaxation. The initial configuration is random, and the average vorticity is near zero. The horizontal red line indicates N = 0. The figure shows an average of over 100 samples. At ρ = 0 and ρ = 1, the t… view at source ↗
Figure 7
Figure 7. Figure 7: Heatmap of the total vorticity N on the ϵ − ρ plane. The heatmap is constructed from the data of the average of over 100 samples of the total vorticity calculated over MCSs from 250 000 to 500 000. The system size is taken to be L = 60. + + - [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The effect of the periodic boundary on the cluster merging process. The system size is taken to be L = 60. ± signs represent topological defects with the vortex charge m = ±1. As time increases from the left figure to the right one, we observe a single large cluster that connects itself across the boundary. Then, the m = −1 defect emerges at the connected point. After relaxation, the two topological defect… view at source ↗
Figure 9
Figure 9. Figure 9: Snapshots after the same elapsed time (100 000 MCSs) with ϵ = 1 for different system sizes L from the random initial configuration. In larger systems, the aggregation of clusters takes much longer MCSs. Comparing [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Time evolution of the cumulative distribution of cluster sizes for ρ = 0.5 for different system sizes L. The distributions are obtained by averaging over 100 samples. The horizontal axis represents the cluster size s, and the vertical axis represents the cumulative distribution P(s) defined in the text. The distributions at four different time steps, t = 1 000, 10 000, 100 000, and 1 000 000, are superimp… view at source ↗
Figure 10
Figure 10. Figure 10: shows the cumulative distribution of cluster sizes for the case of particle density ρ = 0.5, obtained by averaging over 100 samples for the several system sizes. In all system sizes, at early time steps of around t = 1 000, small clusters dominate and large clusters are scarcely observed. However, cluster merging and growth occur as time pro￾gresses, and the cumulative distribution shifts toward larger si… view at source ↗
Figure 11
Figure 11. Figure 11: Time evolution of the cumulative distribution of cluster sizes for ρ = 0.2 for different system sizes L. The distributions are obtained by averaging over 100 samples. The horizontal and the vertical axes, time steps, and system sizes are the same as shown in [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Plot of the area fraction of the maximum cluster size, ϕmax = ⟨S m⟩/L 2 , as a function of the particle density ρ. The cases for system sizes L = 60, 120, and 180 are compared. The result of the linear fitting, ϕmax ≃ 1.15ρ − 0.14, is also shown. ⟨S m⟩ is averaged over 100 samples after the relaxation of 2 000 000 MCSs. The error bar of each data point is smaller than the size of the markers. where ⟨S m⟩ … view at source ↗
Figure 13
Figure 13. Figure 13: Time evolution of the maximum cluster size. The horizontal axis represents the number of MCSs in units of 106 MCSs, and the vertical axis represents the area fraction of the maximum cluster, ϕmax = ⟨S m⟩/L 2 . The data for system sizes L = 60, 120, 180, and 240 are compared for several densities: (a) ρ = 0.3, (b) ρ = 0.4, (c) ρ = 0.5, and (d) ρ = 0.6. evolution, we employ ϕmax = ⟨S m⟩/L 2 , as defined in … view at source ↗
Figure 14
Figure 14. Figure 14: Scaling transformation of [PITH_FULL_IMAGE:figures/full_fig_p021_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Time evolution of the area fraction of the maximum cluster size ϕmax(t) for ρ = 0.5, with the time axis scaled by L 3 . (a) In the double-logarithmic plot, the data would form a straight line if a power-law relaxation is valid, but the results do not support a simple power-law. (b) In the semi-logarithmic plot, two distinct linear regions are observed, corresponding to a rapid initial relaxation and a slo… view at source ↗
Figure 16
Figure 16. Figure 16: Two exponential relaxation processes and snapshots of corresponding times. In rapid relaxation, the dynamics of small clusters are dominated. In contrast, the latter relaxation is caused by the growth of a single large cluster. Finally, we compare the two-stage relaxation for the case of other densities ρ = 0.3, ρ = 0.4, and ρ = 0.6 in [PITH_FULL_IMAGE:figures/full_fig_p023_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Time evolution of the area fraction of the maximum cluster size ϕmax(t), shown on a semi-logarithmic plot, for system sizes L = 60, 120, 180, and 240 under different density conditions: (a) ρ = 0.3, (b) ρ = 0.4, (c) ρ = 0.5, and (d) ρ = 0.6. In all cases, the data indicate that the growth is described as the sum of two exponential functions. Interestingly, when the time axis is normalized by L 3 , the tra… view at source ↗
read the original abstract

We propose an ``active XY model'' that incorporates key elements of both the classical XY model and the Vicsek model to study the role of topological defects in active matter systems. This model features self-propelled particles with XY spin degrees of freedom on a lattice and introduces a self-propulsion parameter controlling the directional bias of particle motion. Using numerical simulations, we demonstrate that self-propulsion induces motility-induced phase separation (MIPS), where particles aggregate into clusters around topological defects with positive vortex charge. In contrast, negative charge defects tend to dissipate. We analyze the evolution of these clusters and show that their growth follows a two-stage exponential relaxation process, with characteristic time scaling as $\tau \sim L^{3}$ with the system size $L$, reminiscent of first-order phase separation in equilibrium systems. Our results highlight the important role of topological defects in phase separations and clustering behavior in active systems, bridging nonequilibrium dynamics and equilibrium theory.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces a lattice-based active XY model combining elements of the classical XY model and the Vicsek model, controlled by a self-propulsion parameter. Numerical simulations show that self-propulsion induces motility-induced phase separation (MIPS), with particles aggregating into clusters around topological defects of positive vortex charge while negative-charge defects dissipate. Cluster growth is reported to follow a two-stage exponential relaxation whose characteristic time scales as τ ∼ L³ with system size L, analogous to first-order phase separation in equilibrium systems.

Significance. If the numerical observations prove robust, the work provides a useful lattice platform for examining how topological defects mediate clustering and phase separation in active matter. The reported L³ scaling for the relaxation time offers a concrete link to equilibrium first-order transition kinetics and could stimulate further studies of defect-driven nonequilibrium dynamics.

major comments (2)
  1. [Model Definition and Simulation Setup] Model and simulation sections: the assignment of local vorticity via plaquette angle sums and the subsequent classification into positive versus negative defects is central to the claim that only positive-charge defects nucleate stable MIPS clusters. No systematic variation of the vorticity cutoff, plaquette size, or coarse-graining length is described, nor is a comparison to an off-lattice or continuum implementation provided; without these controls the reported positive/negative distinction and the associated aggregation may contain lattice-scale artifacts.
  2. [Results on Cluster Growth and Scaling] Results on cluster evolution: the two-stage exponential fit and the τ ∼ L³ scaling are load-bearing for the analogy to equilibrium first-order phase separation. The manuscript does not report checks for robustness against finite-size effects, alternative initial conditions, or variations in the self-propulsion parameter range; these omissions leave open whether the scaling is an intrinsic feature of the active dynamics or sensitive to the chosen simulation protocol.
minor comments (2)
  1. [Abstract] Abstract: the phrase 'reminiscent of first-order phase separation in equilibrium systems' would benefit from a brief citation or explicit comparison to a specific equilibrium model (e.g., Ising or lattice gas).
  2. [Figures and Captions] Figure captions and text: ensure that error bars or standard deviations are shown for all reported scaling exponents and relaxation times.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and for recognizing the potential of our active XY model in studying defect-mediated dynamics in active matter. We address each major comment below and outline the revisions we will make to strengthen the manuscript.

read point-by-point responses

  1. Referee: Model and simulation sections: the assignment of local vorticity via plaquette angle sums and the subsequent classification into positive versus negative defects is central to the claim that only positive-charge defects nucleate stable MIPS clusters. No systematic variation of the vorticity cutoff, plaquette size, or coarse-graining length is described, nor is a comparison to an off-lattice or continuum implementation provided; without these controls the reported positive/negative distinction and the associated aggregation may contain lattice-scale artifacts.

    Authors: We agree that additional controls would enhance the robustness of our findings regarding the role of positive versus negative defects. In the revised version, we will include a systematic study varying the vorticity cutoff threshold (e.g., from 0.5 to 2.0 in units of 2π) and plaquette sizes. We will also perform simulations with different coarse-graining lengths to confirm that the aggregation around positive defects persists. While a full off-lattice comparison is beyond the scope of this lattice-focused study, we will add a brief discussion of how our results relate to continuum active matter models with topological defects. revision: yes

  2. Referee: Results on cluster evolution: the two-stage exponential fit and the τ ∼ L³ scaling are load-bearing for the analogy to equilibrium first-order phase separation. The manuscript does not report checks for robustness against finite-size effects, alternative initial conditions, or variations in the self-propulsion parameter range; these omissions leave open whether the scaling is an intrinsic feature of the active dynamics or sensitive to the chosen simulation protocol.

    Authors: We appreciate this point and will strengthen the results section accordingly. We will add data for multiple system sizes to better demonstrate finite-size scaling, include results from different initial conditions (e.g., random vs. clustered starts), and explore a broader range of self-propulsion parameters to show that the τ ∼ L³ scaling is robust. These additions will support the analogy to equilibrium first-order kinetics more convincingly. revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct simulation outputs on an explicitly defined lattice model

full rationale

The paper defines an active XY lattice-gas model with explicit self-propulsion bias and update rules, then reports observations from numerical simulations including defect charge assignment, cluster aggregation around positive vortices, negative defect dissipation, two-stage exponential growth, and τ ∼ L³ scaling. These quantities are computed directly from the simulated configurations and time series; no analytical derivation, parameter fitting to a target observable, or self-citation chain is invoked to obtain the central claims. The model and measurement protocols are self-contained and externally falsifiable by rerunning the simulations, satisfying the criteria for a non-circular finding.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The central claim rests on the new model definition and simulation outcomes; main additions are the self-propulsion bias and lattice interaction rules, with standard assumptions from XY and Vicsek models.

free parameters (1)
  • self-propulsion parameter
    Controls directional bias of particle motion; value chosen to induce the observed MIPS and defect effects.
axioms (2)
  • domain assumption Particles on lattice interact via nearest-neighbor XY spin alignment rules
    Standard background from classical XY model invoked to define the orientation dynamics.
  • domain assumption Self-propulsion adds a directional bias to particle hops without altering spin interactions
    Core modeling choice combining Vicsek motility with XY spins.
invented entities (1)
  • active XY model no independent evidence
    purpose: Hybrid lattice model to study topological defects in active matter
    New construct introduced to bridge equilibrium XY defects with nonequilibrium motility.

pith-pipeline@v0.9.0 · 5695 in / 1541 out tokens · 39176 ms · 2026-05-19T03:35:08.787723+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha: Reviews of Modern Physics 85 (2013) 1143

    work page 2013

  2. [2]

    Bechinger, R

    C. Bechinger, R. Di Leonardo, H. L ¨owen, C. Reichhardt, G. V olpe, and G. V olpe: Re- views of Modern Physics 88 (2016) 045006

    work page 2016

  3. [3]

    J. L. Toner: The Physics of Flocking: Birth, Death, and Flight in Active Matter (Cam- bridge University Press, Cambridge, 2024)

    work page 2024

  4. [4]

    C. W. Reynolds: ACM SIGGRAPH Computer Graphics 21 (1987) 25

    work page 1987

  5. [5]

    Shimoyama, K

    N. Shimoyama, K. Sugawara, T. Mizuguchi, Y . Hayakawa, and M. Sano: Physical Re- view Letters 76 (1996) 3870

    work page 1996

  6. [6]

    Niwa: Journal of Theoretical Biology 181 (1996) 47

    H.-S. Niwa: Journal of Theoretical Biology 181 (1996) 47

    work page 1996

  7. [7]

    Ito and N

    S. Ito and N. Uchida: Europhysics Letters 138 (2022) 17001

    work page 2022

  8. [8]

    Ito and N

    S. Ito and N. Uchida: Journal of the Physical Society of Japan 91 (2022) 064806

    work page 2022

  9. [9]

    M. S. Kumar and P. Philominathan: Biophysical Reviews 2 (2010) 13

    work page 2010

  10. [10]

    Nishiguchi, S

    D. Nishiguchi, S. Shiratani, K. A. Takeuchi, and I. S. Aranson: Proceedings of the Na- tional Academy of Sciences 122 (2025) e2414446122

    work page 2025

  11. [11]

    Kawaguchi, R

    K. Kawaguchi, R. Kageyama, and M. Sano: Nature 545 (2017) 327

    work page 2017

  12. [12]

    Sumino, K

    Y . Sumino, K. H. Nagai, Y . Shitaka, D. Tanaka, K. Yoshikawa, H. Chat´e, and K. Oiwa: Nature 483 (2012) 448

    work page 2012

  13. [13]

    Vicsek, A

    T. Vicsek, A. Czir ´ok, E. Ben-Jacob, I. Cohen, and O. Shochet: Physical Review Letters 75 (1995) 1226

    work page 1995

  14. [14]

    H. E. Stanley: Physical Review Letters 20 (1968) 589

    work page 1968

  15. [15]

    Tobochnik and G

    J. Tobochnik and G. Chester: Physical Review B 20 (1979) 3761

    work page 1979

  16. [16]

    J. M. Kosterlitz and D. J. Thouless: Journal of Physics C: Solid State Physics 6 (1973) 1181

    work page 1973

  17. [17]

    J. M. Kosterlitz: Journal of Physics C: Solid State Physics 7 (1974) 1046

    work page 1974

  18. [18]

    V . L. Berezinskii: Soviet Physics JETP32 (1971) 493

    work page 1971

  19. [19]

    V . L. Berezinskii: Soviet Physics JETP34 (1972) 610

    work page 1972

  20. [20]

    D. R. Nelson: Defects and Geometry in Condensed Matter Physics (Cambridge Univer- sity Press, Cambridge, 2002). 27/28 J. Phys. Soc. Jpn. FULL PAPERS

    work page 2002

  21. [21]

    Shankar, A

    S. Shankar, A. Souslov, M. J. Bowick, M. C. Marchetti, and V . Vitelli: Nature Reviews Physics 4 (2022) 380

    work page 2022

  22. [22]

    M. E. Cates and J. Tailleur: Annual Review of Condensed Matter Physics 6 (2015) 219

    work page 2015

  23. [23]

    Kourbane-Houssene, C

    M. Kourbane-Houssene, C. Erignoux, T. Bodineau, and J. Tailleur: Physical Review Letters 120 (2018) 268003

    work page 2018

  24. [24]

    Adachi and H

    K. Adachi and H. Nakano: Physical Review Research 6 (2024) 033234

    work page 2024

  25. [25]

    Nakano and K

    H. Nakano and K. Adachi: Physical Review Research 6 (2024) 013074

    work page 2024

  26. [26]

    Toner and Y

    J. Toner and Y . Tu: Physical Review Letters75 (1995) 4326

    work page 1995

  27. [27]

    Toner and Y

    J. Toner and Y . Tu: Physical Review E58 (1998) 4828

    work page 1998

  28. [28]

    Toner: Physical Review E 86 (2012) 031918

    J. Toner: Physical Review E 86 (2012) 031918

    work page 2012

  29. [29]

    Haldar, A

    A. Haldar, A. Sarkar, S. Chatterjee, and A. Basu: Physical Review E 108 (2023) L032101

    work page 2023

  30. [30]

    Haldar, A

    A. Haldar, A. Sarkar, S. Chatterjee, and A. Basu: Physical Review E108 (2023) 034114

    work page 2023

  31. [31]

    R. A. Fisher and F. Yates: Statistical tables for biological, agricultural and medical research (Edinburgh: Oliver and Boyd, 1963)

    work page 1963

  32. [32]

    D. E. Knuth: The Art of Computer Programming, volume 2 Seminumerical Algorithms (Addison-Wesley, Reading, Massachusetts, 1998) 3rd ed

    work page 1998

  33. [33]

    Hoshen and R

    J. Hoshen and R. Kopelman: Physical Review B 14 (1976) 3438

    work page 1976

  34. [34]

    Lifshitz and V

    I. Lifshitz and V . Slyozov: Journal of Physics and Chemistry of Solids19 (1961) 35

    work page 1961

  35. [35]

    P. M. Chaikin and T. C. Lubensky:Principles of Condensed Matter Physics (Cambridge University Press, 1995)

    work page 1995

  36. [36]

    Baldan: Journal of Materials Science 37 (2002) 2171

    A. Baldan: Journal of Materials Science 37 (2002) 2171

    work page 2002

  37. [37]

    A. G. Thompson, J. Tailleur, M. E. Cates, and R. A. Blythe: Journal of Statistical Me- chanics: Theory and Experiment 2011 (2011) P02029

    work page 2011

  38. [38]

    Adachi and K

    K. Adachi and K. Kawaguchi: arXiv preprint arXiv:2012.02517 (2020)

    work page arXiv 2012

  39. [39]

    Stenhammar, A

    J. Stenhammar, A. Tiribocchi, R. J. Allen, D. Marenduzzo, and M. E. Cates: Physical Review Letters 111 (2013) 145702. 28/28

    work page 2013