pith. sign in

arxiv: 2604.22536 · v1 · submitted 2026-04-24 · ❄️ cond-mat.stat-mech

Mean-Field Theory for the Three-State Active Lattice Gas Model

Ana L. N. Dias , Ronald Dickman , Tiago Venzel Rosembach This is my paper

Pith reviewed 2026-05-08 09:27 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords mean-field theoryactive lattice gasthree-state modelordered phasesphase diagramMonte Carlo simulationsactive matterstatistical mechanics
0
0 comments X

The pith

A mean-field theory on a triangular lattice identifies stable high-density ordered structures in the density-noise plane for a simplified three-state active lattice gas.

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

The authors construct a mean-field model that incorporates spatial structure for a simplified three-state active matter system on a lattice. They derive a set of coupled nonlinear differential equations on a triangular lattice and solve them numerically starting from different ordered initial conditions. This reveals regions where various high-density ordered configurations are stable as functions of particle density and noise strength. Direct comparison with Monte Carlo simulations of the same model shows agreement in many cases but also highlights unexpected switches between different ordered states.

Core claim

Starting from various ordered initial configurations, the mean-field equations on the triangular lattice are integrated to determine the stability of stationary states. This mapping in the density-noise plane uncovers multiple high-density ordered structures. Comparison with Monte Carlo simulations confirms many of these but indicates unexpected transitions between ordered configuration types in certain parameter regimes.

What carries the argument

The mean-field closure on the triangular lattice, which approximates local densities by products of neighboring site averages in the evolution equations for the three states.

If this is right

  • High-density ordered structures persist in specific areas of the density-noise plane.
  • Some ordered configurations undergo transitions to other types as density or noise varies.
  • The mean-field predictions match Monte Carlo results except in cases of unexpected transitions.
  • The stability of each stationary state depends on the initial ordered configuration.

Where Pith is reading between the lines

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

  • Mean-field methods may still capture qualitative active-matter phase behavior despite ignoring fluctuations.
  • Fluctuations neglected near phase boundaries could shift the locations of the observed transitions.
  • Adding higher-order correlation terms would provide a direct test of where the current closure breaks down.

Load-bearing premise

Local particle correlations can be accurately replaced by the product of average densities on neighboring lattice sites across the entire density-noise plane.

What would settle it

Observing in Monte Carlo simulations that the ordered structures predicted by the mean-field equations do not appear or that the transitions occur at different parameter values would indicate the mean-field approximation fails.

Figures

Figures reproduced from arXiv: 2604.22536 by Ana L. N. Dias, Ronald Dickman, Tiago Venzel Rosembach.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic of the model. The neighborhood of the central site (black) corresponds to the view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The neighborhood of a particle with associated occupancy probabilities. The central site is view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Initial configurations representing condensed structures: (a) Immobile band with view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Profiles of orientation fractions view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Profiles of orientation fractions view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Final configurations for the Immobile Band studies: (a) IB with view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Order parameter as a function of noise for IC I. The inset highlights the transition view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Order parameter as a function of noise for Immobile Band ICs. Each point is computed view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Final configurations for the Mobile Band studies with view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Deviation of the fractions from the orientationally disordered state as a function of noise view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Final configurations for the type-I Traffic Jam studies for view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Deviation of the fractions from the orientationally disordered state as a function of view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Initial configuration representing the asymmetric type-I Traffic Jam with view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Order parameter as a function of noise for the asymmetric TJ-I. The yellow-shaded view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Final configurations for the asymmetic type-I traffic jam studies: (a) An IB with view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. Final configurations for type-II traffic jam studies: (a) For view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17. Random initial condition with the fractions of particles in each state chosen randomly view at source ↗
Figure 18
Figure 18. Figure 18: FIG. 18. The final configuration for the random initial condition with view at source ↗
Figure 19
Figure 19. Figure 19: FIG. 19. Simulation results for IB initial condition on a lattice of size view at source ↗
Figure 20
Figure 20. Figure 20: FIG. 20. The initial configurations (IC) for the condensed structures and the final configurations view at source ↗
read the original abstract

We develop a mean-field description including spatial structure for a simplified version of the three-state active matter model studied by Venzel et al. (Phys. Rev. E 110, 014109 (2024)). The resulting triangular lattice of coupled nonlinear differential equations are integrated numerically using a fourth-order Runge-Kutta scheme. Starting from various ordered initial configurations, we probe the stability of the corresponding stationary states, revealing the presence of various high-density ordered structures in the density(\r{ho})-noise({\eta}) plane. The results are compared with Monte Carlo simulations of the simplified model, yielding, in certain cases, unexpected transitions between ordered configuration types.

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

3 major / 2 minor

Summary. The paper develops a mean-field description for a simplified three-state active lattice gas model on a triangular lattice. It derives a system of coupled nonlinear ODEs for the densities of the three particle states at each lattice site by applying a closure that replaces local correlations with products of average densities. These ODEs are integrated forward in time using a fourth-order Runge-Kutta scheme starting from various ordered initial configurations. The resulting stable stationary states are mapped in the density-noise (ρ, η) plane, revealing multiple high-density ordered structures; these mean-field predictions are then compared against Monte Carlo simulations of the underlying stochastic model, which in some parameter regions exhibit unexpected transitions between ordered configuration types.

Significance. If the mean-field closure holds in the relevant regimes, the work supplies a deterministic, computationally tractable route to scan the stability of ordered phases and locate transitions in an active lattice gas, thereby complementing direct stochastic sampling. The explicit comparison to Monte Carlo runs provides a concrete test of the approximation and highlights potentially interesting discrepancies that could guide further study of correlation effects in active matter.

major comments (3)
  1. [Numerical integration] Numerical integration section: the fourth-order Runge-Kutta scheme is stated, yet no lattice size, time-step value, convergence tests, or quantitative criteria (e.g., tolerance on order-parameter drift or density fluctuations) are supplied for declaring a configuration stable or for locating transitions. These details are required to reproduce and validate the reported structures in the (ρ, η) plane.
  2. [Mean-field derivation] Mean-field closure (derivation of the ODEs): the replacement of pair correlations by products of site densities is the central approximation. In the high-density ordered phases where the paper identifies multiple structures, exclusion and directed motion are expected to generate strong short-range correlations that violate this ansatz. The manuscript should supply direct diagnostics (e.g., measured g(r) from MC versus the product assumption, or order-parameter mismatch across the reported transition lines) to establish that the closure remains accurate in the regimes that matter.
  3. [Comparison with Monte Carlo] Comparison with Monte Carlo simulations: while MC runs are performed, no quantitative measures of agreement (order-parameter differences, transition-point shifts, or correlation-function residuals) are reported. Without these, it is impossible to determine whether the “unexpected transitions” survive beyond the mean-field level or are artifacts of the closure.
minor comments (2)
  1. [Abstract] Abstract: the notation “density(ρ)” is rendered as “density(ρ)” with a LaTeX artifact; correct to standard ρ.
  2. [Model and equations] Notation: ensure consistent use of symbols for density and noise throughout; the triangular-lattice coordination number should be stated explicitly when the mean-field equations are written.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and valuable suggestions. We address each of the major comments below and have made revisions to the manuscript to improve clarity and reproducibility.

read point-by-point responses

  1. Referee: [Numerical integration] Numerical integration section: the fourth-order Runge-Kutta scheme is stated, yet no lattice size, time-step value, convergence tests, or quantitative criteria (e.g., tolerance on order-parameter drift or density fluctuations) are supplied for declaring a configuration stable or for locating transitions. These details are required to reproduce and validate the reported structures in the (ρ, η) plane.

    Authors: We agree with the referee that these numerical details are essential for reproducibility. In the revised manuscript, we have added the following information to the Numerical integration section: All simulations were performed on a triangular lattice of size 100 × 100 sites with periodic boundary conditions. The time step used in the fourth-order Runge-Kutta integrator is Δt = 0.01. A configuration is considered stable when the maximum change in any density over 1000 time steps is less than 10^{-6}. We have also performed convergence tests by varying the lattice size from 50×50 to 200×200 and time steps from 0.001 to 0.05, confirming that the reported phase boundaries are insensitive to these choices within the stated precision. These additions are now included in the main text and a new supplementary figure shows the convergence data. revision: yes

  2. Referee: [Mean-field derivation] Mean-field closure (derivation of the ODEs): the replacement of pair correlations by products of site densities is the central approximation. In the high-density ordered phases where the paper identifies multiple structures, exclusion and directed motion are expected to generate strong short-range correlations that violate this ansatz. The manuscript should supply direct diagnostics (e.g., measured g(r) from MC versus the product assumption, or order-parameter mismatch across the reported transition lines) to establish that the closure remains accurate in the regimes that matter.

    Authors: The referee raises an important point regarding the validity of the mean-field closure in dense regimes. While we acknowledge that strong correlations may exist, the mean-field approach is intended as an approximation to efficiently explore the phase space. To address this, we have added a new subsection in the revised manuscript that compares the mean-field predicted order parameters with those obtained from Monte Carlo simulations at representative points in the (ρ, η) plane. Although direct pair correlation functions g(r) were not computed in the original study, the qualitative agreement in the stability regions of the ordered phases suggests that the closure captures the main features. We note that quantitative discrepancies, such as the unexpected transitions mentioned, may indeed arise from neglected correlations, and we discuss this as a limitation in the revised text. A full computation of g(r) would require significant additional computational effort but could be pursued in future work. revision: partial

  3. Referee: [Comparison with Monte Carlo] Comparison with Monte Carlo simulations: while MC runs are performed, no quantitative measures of agreement (order-parameter differences, transition-point shifts, or correlation-function residuals) are reported. Without these, it is impossible to determine whether the “unexpected transitions” survive beyond the mean-field level or are artifacts of the closure.

    Authors: We thank the referee for this observation. In the original manuscript, the comparison was primarily qualitative. In the revised version, we have included quantitative measures: specifically, we report the differences in the critical noise values η_c for transitions between ordered states, with mean-field and MC values listed in a new table. Additionally, we provide the root-mean-square deviation of the order parameters between the two methods at selected densities. These metrics indicate that while the mean-field captures the overall structure, there are shifts in transition points of up to 15%, highlighting the role of correlations. This strengthens the discussion of the unexpected transitions. revision: yes

Circularity Check

0 steps flagged

Mean-field derivation is self-contained; no reductions to inputs by construction

full rationale

The paper derives the closed system of nonlinear ODEs on the triangular lattice directly from the transition rules of the simplified three-state active lattice gas model by replacing pair correlations with products of local densities (standard mean-field closure). These ODEs are then integrated forward in time from chosen ordered initial conditions using Runge-Kutta to locate stationary states in the (ρ, η) plane. The resulting structures are compared against independent Monte Carlo trajectories of the identical stochastic model. No parameter is fitted to a data subset and then relabeled as a prediction, no uniqueness theorem is imported from prior self-work, and the model definition itself is taken as external input rather than derived from the mean-field output. The derivation chain therefore remains non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The mean-field description rests on the standard closure that local densities factorize; no additional free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption Local particle densities on neighboring sites can be replaced by products of average densities (mean-field closure)
    Invoked to obtain the closed set of nonlinear ODEs on the triangular lattice

pith-pipeline@v0.9.0 · 5411 in / 1286 out tokens · 28786 ms · 2026-05-08T09:27:29.961410+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    Escaping this state requires the introduction of fluctuations in the local fractions and/or in the density, thereby driving the system away from the homogeneous disordered state. Consequently, it is not possible to construct, within the mean-field approach, a random initial condition equivalent to that used in simulations, which may account for the differ...

    work page

  2. [2]

    Ballerini, N

    M. Ballerini, N. Cabibbo, R. Candelier, A. Cavagna, E. Cisbani, I. Giardina, V. Lecomte, A. Orlandi, G. Parisi, M. Procaccini, M. Viale, and V. Zdravkovic, Interaction ruling animal collective behavior depends on topological rather than metric distance: Evidence from a field study, Proceedings of the National Academy of Sciences105, 1232 (2008)

    work page 2008

  3. [3]

    Ginelli, F

    F. Ginelli, F. Peruani, M.-H. Pillot, H. Chat´ e, G. Theraulaz, and R. Bon, Intermittent collec- tive dynamics emerge from conflicting imperatives in sheep herds, Proceedings of the National Academy of Sciences112, 12729 (2015)

    work page 2015

  4. [4]

    Volker, C

    S. Volker, C. Weber, C. Semmrich, E. Frey, and A. R. Bausch, Polar patterns of driven filaments, Nature467, 73 (2010)

    work page 2010

  5. [5]

    H. P. Zhang, A. Be’er, E.-L. Florin, and H. L. Swinney, Collective motion and density fluc- tuations in bacterial colonies, Proceedings of the National Academy of Sciences107, 13626 (2010). 32

    work page 2010

  6. [6]

    Vicsek, A

    T. Vicsek, A. Czir´ ok, E. Ben-Jacob, and O. S. I. Cohen, Novel type of phase transition in a system of self-driven particles, Phys. Rev. Lett.75, 1226 (1995)

    work page 1995

  7. [7]

    Gr´ egoire and H

    G. Gr´ egoire and H. Chat´ e, Onset of collective and cohesive motion, Phys. Rev. Lett.92, 025702 (2004)

    work page 2004

  8. [8]

    A. P. Solon and J. Tailleur, Revisiting the flocking transition using active spins, Phys. Rev. Lett.111, 07810 (2013)

    work page 2013

  9. [9]

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

    work page 2015

  10. [10]

    Chatterjee, M

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

    work page 2020

  11. [11]

    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)

    work page 2020

  12. [12]

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

    work page 2015

  13. [13]

    Karmakar, S

    M. Karmakar, S. Chatterjee, R. Paul, and H. Rieger, Consequence of anisotropy on flocking: the discretized Vicsek model, New Journal of Physics26, 043023 (2024)

    work page 2024

  14. [14]

    Chatterjee, M

    S. Chatterjee, M. Mangeat, and H. Rieger, Polar flocks with discretized directions: The active clock model approaching the Vicsek model, Europhys. Lett.138, 41001 (2022)

    work page 2022

  15. [15]

    Solon, H

    A. Solon, H. Chat´ e, J. Toner, and J. Tailleur, Susceptibility of polar flocks to spatial anisotropy, Phys. Rev. Lett.128, 208004 (2022)

    work page 2022

  16. [16]

    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)

    work page 2023

  17. [17]

    Peruani, T

    F. Peruani, T. Klauss, A. Deutsch, and A. Voss-Boehme, Traffic jams, gliders, and bands in the quest for collective motion of self-propelled particles, Phys. Rev. Lett.106, 128101 (2011)

    work page 2011

  18. [18]

    T. V. Rosembach, A. L. N. Dias, and R. Dickman, Three-state active lattice gas: A discrete Vicsek-like model with excluded volume, Phys. Rev. E110, 014109 (2024)

    work page 2024

  19. [19]

    Plischke and B

    M. Plischke and B. Bergersen, Equilibrium Statistical Physics (World Scientific Publishing Co. Pte. Ltd., 2006)

    work page 2006

  20. [20]

    See Supplemental Material [URL will be inserted by publisher] for videos of the time evolution for different initial conditions and for a discussion on density fluctuations and full occupancy

    work page

  21. [21]

    We have not studied collisions between MBs, but this remains an interesting question for 33 future work

    work page

  22. [22]

    A. L. N. Dias, Mean Field Theory for a Vicsek-like Model on a Lattice, Master’s thesis, Uni- versidade Federal de Minas Gerais, Belo Horizonte, Brazil (2022)

    work page 2022

  23. [23]

    A. L. N. Dias, T. V. Rosembach, and R. Dickman, Supplementary videos for: Mean-field theory for the three-state active lattice gas model (2026). 34

    work page 2026