{"id":"4b1c95be-4e25-43ca-bc7f-0acd75942e26","arxiv_id":"2502.07054","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a nematic-aligning active-particle model, high individual persistence destroys global nematic order even when alignment interactions are strong, leaving only local order.","lead":"This paper uses a lattice model of self-propelled particles that align head-to-tail, and shows that the tendency of each particle to keep moving straight can prevent the whole population from lining up in one global direction. It matters because it identifies individual persistence, not just particle interactions, as a control knob for collective order in bacterial swarms and active materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Density-independence claim is unsupported: only ρ=0.2 is simulated, and the linear-stability argument for density independence does not cover the nonlinear ordered-phase stability that underlies the A-B transition.","rationale":"The reader's weakest_assumption is exactly the point that would sink the paper's central claim if wrong. The paper's own evidence is one density. The linear stability of the disordered state is not sufficient for density independence because the ordered-state stability, which controls the A-B boundary, depends on the self-consistent mean-field occupations f1 and f2; these are density-dependent, so the effective persistence threshold is not captured by the products αρ and βρ alone. The parameter mismatch in Fig. 7 strengthens the concern: the only ordered-state stability calculation is at β=0.8, not at the β=4 used for the simulated A-B transition, so the theoretical support for the persistence-induced frustration is weaker than claimed. Nevertheless, the raw simulation result at ρ=0.2 is plausible and the central qualitative behavior is likely correct; the issue is over-generalization. The paper should remain CONDITIONAL pending a density scan or a nonlinear density-collapse argument that directly addresses the ordered states, not just the homogeneous instability thresholds.","tokens_in":16867,"tokens_out":12844,"duration_ms":104752,"concrete_test":"Re-run the LGCA simulations at ρ=0.1 and ρ=0.4 (same L=120, periodic boundary conditions, same initial conditions and update rules) with β=4, scanning α from 0 to 10 in steps of 0.5. Measure SL after 2000 time steps, averaging over at least 100 independent runs. If the transition α* where SL drops from near 1 to ~0.5 shifts by more than ~20% relative to ρ=0.2, or if SL remains near 1 for α≥5 at ρ=0.4, then the density-independence claim and the generality of the central claim are refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central assertion—global nematic order is only possible for weakly persistent particles—is established by the phase diagram in Fig. 3, obtained solely at ρ=0.2 on a 120×120 lattice. Section III asserts that 'different densities lead to qualitatively the same stationary behavior,' citing a linear stability analysis in which the homogeneous-state thresholds depend on βρ and αρ. This argument covers only infinitesimal perturbations of the disordered state. The ordered nematic state Q+, whose persistence-induced destabilization is the stated mechanism for the A→B transition, is analyzed through the mean-field propagator in Eq. (14), which depends on the ordered-state occupations f1 and f2. These steady-state occupations are density-dependent solutions of Eq. (7), so the stability of Q+ against persistence is not a function of αρ and βρ alone. No simulations at any other density are reported, and no nonlinear stability argument is given. If at a still-dilute density (e.g., ρ=0.4) the A→B transition moves substantially in α, or if high-persistence populations retain global nematic order, the paper's headline claim as stated is false. Additionally, the Q+ stability computation in Fig. 7 uses β=0.8, whereas the simulated A→B transition in Fig. 4 uses β=4; thus the provided analysis does not even directly address the transition it claims to explain.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a lattice-gas cellular automaton (LGCA) model of self-propelled particles with nematic alignment, incorporating an independent parameter α for individual persistence of motion and β for the strength of nematic alignment interactions. Simulations at a single density ρ=0.2 on a 120×120 periodic lattice map a phase diagram in the (β, α) plane (Fig. 3), identifying a region of global nematic order (A), a region of partial or local nematic order (B) at high persistence and high alignment, a polar-ordered region (C) at high persistence and low alignment, and a disordered region. A mean-field linear stability analysis of the homogeneous state yields threshold conditions depending on the products βρ and αρ (Eq. (13)), and a stability analysis of a homogeneous nematic steady state Q+ (Eq. (14)) suggests that large persistence destabilizes global nematic order. The authors conclude that individual persistence can frustrate global nematic ordering while allowing local order, with qualitative comparisons to experiments on Myxococcus xanthus and other systems.","tokens_in":17149,"tokens_out":4869,"duration_ms":43751,"significance":"If the central claim holds, the paper offers a minimal, analytically tractable model showing that an individual-level motility trait—persistence—can control the global spatiotemporal organization of a collectively migrating population. The clean separation of persistence and alignment into two independent control parameters is a conceptual strength, and the combination of simulation with linear stability analysis makes the paper useful for further theoretical work. The mean-field derivations are mostly transparent and reproduce known instabilities as limiting cases. The claim is falsifiable and directly testable by additional simulations at other densities and by repeating the ordered-state stability analysis at the parameters of the simulated transition.","major_comments":[{"comment":"The statement in Section III that \"different densities lead to qualitatively the same stationary behavior\" is supported only by the linear stability of the homogeneous state, where thresholds depend on the products βρ and αρ. This argument does not cover the ordered-state stability that controls the A→B transition, which is the central evidence for the headline claim. The stability of Q+ in Eq. (14) depends on the steady-state occupations f1 and f2, which are density-dependent solutions of Eq. (7); therefore the persistence-induced destabilization of global nematic order is not a function of αρ and βρ alone. No simulations at other densities are reported, so the claim that high persistence prevents global nematic order has not been shown to hold beyond ρ=0.2. Additional simulations at, e.g., ρ=0.1 and ρ=0.4 would directly test the density independence of the A→B transition.","section":"Section III; Section IV, Eqs. (7) and (14)"},{"comment":"The stability analysis of the ordered state Q+ is performed at β=0.8 (Fig. 7), whereas the simulated A→B transition in Fig. 4 is obtained at β=4. Since the ordered-state occupations f1 and f2, and hence the stability of Q+ against persistence, depend on β, the provided analysis does not directly explain the transition that it claims to explain. The authors should repeat the eigenvalue calculation at β=4 or at least show how the instability threshold in α varies with β and compare it with the simulated phase boundary.","section":"Section IV, Fig. 7 vs. Fig. 4"},{"comment":"The mean-field derivation leading to Eq. (13) assumes a dilute system with negligible interference among reorientations at a node, but the simulated density ρ=0.2 corresponds to an average of 1.2 particles per node out of six channels. The validity of the dilute approximation at this density is not quantified. A direct comparison between the predicted instability thresholds (3αρ=1, 9βρ=1) and the first appearance of order in simulations would help establish that the linear analysis applies in the regime used for the phase diagram.","section":"Section IV, Eq. (13)"}],"minor_comments":[{"comment":"The caption contains a typo: \"incrased persistence\" should be \"increased persistence\".","section":"Section III, Fig. 4 caption"},{"comment":"The text states that patterns were characterized on \"a square lattice\", but the model is defined on a triangular (hexagonal) lattice; this is likely a typo and should be corrected.","section":"Section V"},{"comment":"The phrase \"the five different phase transition among the four observed steady states\" is unclear; the paper describes four distinct phases, not five transitions. Please revise the wording.","section":"Section V"},{"comment":"The phase diagram in Fig. 3 does not report error bars or the number of independent simulations used to compute the order parameters; adding this information would allow the reader to judge the statistical significance of the phase boundaries.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is in scope and the topic is suitable for the journal. The central idea is interesting and the model is well-defined. The main concern is the gap between the simulation evidence (single density) and the general claim, as well as the parameter mismatch in the theoretical explanation of the A→B transition. These are fixable with additional simulations and analysis, so I recommend major revision rather than rejection. No concerns about citation practice or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is a solid lattice-gas cellular automaton study that I think is worth a serious referee, with one clearly overreaching claim. The new thing here is a clean separation of single-particle persistence (alpha) from nematic alignment strength (beta) in a single model, and a phase diagram (Fig. 3) showing that high persistence destroys the globally ordered nematic phase (region A) and replaces it with a partially ordered or polar phase. That qualitative result—global nematic order only for weak persistence—is plausible and is backed by both simulation and a mean-field linear stability analysis. The mean-field treatment is competently done: the homogeneous-state instabilities are computed from the model's own transitions, and the ordered-state Q+ stability at low vs. high alpha (Fig. 7) is consistent with the A-to-B transition. The links to Myxococcus and other bacterial patterns are appropriately cautious, not overclaimed.\n\nThe soft spots are real but mostly addressable. The biggest one is the density-independence assertion in Sec. III (\"different densities lead to qualitatively the same stationary behavior\"). The support is the linear stability of the homogeneous state, where thresholds depend on alpha*rho and beta*rho. That covers infinitesimal perturbations of the disordered state only; it says nothing about the nonlinear stability of the ordered states, whose occupations f1 and f2 depend on rho in a nontrivial way. No simulations at other densities are shown. So the generalization is unsupported. If the A-B boundary shifts a lot with rho at dilute densities, the central claim would need qualification. I don't have a strong prior that it fails, but the paper should either simulate another density or soften the claim.\n\nSecond, the Q+ stability analysis in Fig. 7 uses beta=0.8, whereas the simulated A-B transition plotted in Fig. 4 uses beta=4. The analysis is qualitative, so this is not fatal, but it does mean the calculation is not directly checking the transition it's invoked to explain. A referee should ask whether the instability persists at beta=4. Third, Fig. 3 has no error bars and no code/data are shipped; that's a minor issue for a simulation paper, but it would help.\n\nThe derivation is not circular: the instabilities are computed from the model's own rules and compared with independent simulations. No experimental constants are fitted. The citation pattern is fine, including the authors' prior LGCA work, which is the correct formalism to build on.\n\nBottom line: the central physical claim is probably right, but as written the density-independence sentence is too strong and the stability analysis is one parameter set away from the claimed transition. With those fixed, this would be a solid contribution to active matter and biophysics. It deserves peer review, and I'd bring it to a reading group.","headline":"A clean LGCA model showing that high persistence frustrates global nematic order, but the density-independence claim is oversold and the ordered-state analysis is done off the simulated transition.","tokens_in":17662,"tokens_out":3456,"would_cite":true,"duration_ms":29684,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In a lattice model of nematically aligning self-propelled particles, this paper shows that high individual persistence prevents global nematic order even under strong alignment interactions, leaving only local order.","keywords":["active matter","nematic order","persistence of motion","lattice-gas cellular automaton","collective motion","mean-field stability analysis","orientational order","self-propelled particles"],"falsifier":"Run the same LGCA on a 120 by 120 lattice at densities $\\rho = 0.1$ and $\\rho = 0.4$ across the $(\\beta, \\alpha)$ parameter plane and measure the nematic order parameter $S_L$ in the stationary state; if at either density the region of high global nematic order extends to large $\\alpha$, or the transition from region A to region B disappears, then persistence-frustrated ordering is not a general dilute-regime property. A cheaper check would be to fix $\\beta = 4$, vary $\\alpha$ at both densities, and compare where $S_L$ drops.","tokens_in":16673,"feed_emoji":"🦠","tokens_out":7135,"duration_ms":60641,"temperature":0.7,"pith_summary":"This paper uses a lattice-gas cellular automaton (a discrete space-and-time model in which particles hop along six velocity channels on a triangular lattice) to ask how the persistence of an individual's own motion affects the collective order of a population that aligns nematically, i.e., head-to-tail rather than head-to-head. The central finding is that global nematic order only appears when particles are weakly persistent: raising the persistence parameter $\\alpha$ while keeping the alignment strength $\\beta$ high destroys the globally ordered nematic band state and leaves only partial, local nematic order. The paper establishes this by simulation of a 120 by 120 lattice at density $\\rho = 0.2$ and by a mean-field linear stability analysis of the homogeneous state. The practical point is that population-level order is not set by interactions alone; the intrinsic movement style of the individuals can act as an independent control parameter capable of frustrating or promoting a given collective phase.","feed_headline":"Stubborn motion blocks global order in head-to-tail crowds","feed_subtitle":"Simulations and mean-field theory show individual persistence can keep alignment local even when interactions are strong.","key_machinery":"The load-bearing object is the reorientation Hamiltonian of the LGCA, written as $H = H_{\\rm pers} + H_{\\rm align}$, where $H_{\\rm pers} = -\\alpha \\sum_{\\ell,j} (\\vec{c}_\\ell \\cdot \\vec{c}_j) s^I_\\ell s_j$ biases a particle to keep its current velocity channel and $H_{\\rm align} = -\\beta \\sum_{\\vec{r}' \\in N_{\\vec{r}}} \\sum_{\\ell,j} (\\vec{c}_\\ell \\cdot \\vec{c}_j)^2 s^I_\\ell s_j^{\\vec{r}'}$ biases reorientation toward the average nematic axis of neighboring nodes. The argument then proceeds through the mean-field Boltzmann propagator $\\Gamma_{\\ell,j}(\\vec{k})$, a matrix that maps Fourier modes of perturbations under the linearized mean-field dynamics; its eigenvalues decide whether small perturbations of a steady state grow. The key relations are the instability thresholds $3\\alpha\\rho = 1$ for flux (polar) modes and $9\\beta\\rho = 1$ for nematic modes, and the finding that large $\\alpha$ destabilizes the ordered nematic steady state $Q_+$. This is the mechanism by which individual persistence is shown to switch off global nematic order.","core_discovery":"On the paper's own terms, the discovery is a two-parameter phase structure for nematically aligning self-propelled particles. In the $(\\beta, \\alpha)$ plane, with $\\beta$ the sensitivity to nematic alignment and $\\alpha$ the persistence of individual motion, four regimes appear: a disordered random-walk state at low $\\beta$ and low $\\alpha$; a globally ordered nematic state ($S_L \\approx 1$, $S_F \\approx 0$) at high $\\beta$ and low $\\alpha$; a polar-ordered state at low $\\beta$ and high $\\alpha$; and a partially ordered, frustrated state at high $\\beta$ and high $\\alpha$ in which nematic bands form but no global nematic order is achieved. The linear stability analysis explains the frustration: the homogeneous disordered state is destabilized by nematic modes when $9\\beta\\rho = 1$ and by flux/polar modes when $3\\alpha\\rho = 1$, where $\\rho$ is density; at high $\\alpha$ the nematic steady state itself develops unstable polar and nematic modes, which is exactly the persistence-induced loss of global ordering seen in simulation.","pith_inferences":["Editorial inference: if the persistence-frustration mechanism is generic, then in engineered active materials one could toggle between global and local orientational order by changing the noise or persistence of individual agents, without changing the alignment interaction at all.","Editorial inference: the linear-stability density-independence argument does not by itself guarantee that the nonlinear phases (bands, clusters) are density-independent; simulating the phase diagram at two other dilute densities would be the direct test of whether the antagonism between $\\alpha$ and global nematic order survives outside the homogeneous-state analysis.","Editorial inference: because the persistence term is built from exponentially correlated Gaussian noise, the result may not carry over to power-law or L\\'evy persistent motion, which is observed in some bacteria; whether heavy-tailed persistence also frustrates nematic order is an open question the paper does not address.","Editorial inference: the model suggests a biological reading that the paper leaves implicit—populations that require global nematic sheets should keep their individual persistence low, whereas populations that tolerate streams and local order can be highly persistent."],"forward_implications":["At fixed dilute density, increasing individual persistence can push a population from a globally nematic-ordered band state into a partially ordered state, so persistence acts as a control parameter for the degree of collective order.","The linear stability thresholds $3\\alpha\\rho = 1$ and $9\\beta\\rho = 1$ imply that density and interaction strength enter only through their products, so the same phase behavior should reappear if both density and coupling are rescaled in the dilute regime.","When both persistence and nematic sensitivity are high, the ordered nematic steady state becomes unstable in the mean-field analysis, which explains the observed transition between the globally ordered region A and the partially ordered region B.","The model's patterns—a single nematic band at low persistence, a network of bands at intermediate values, and polar clusters at high persistence—mimic the range of structures seen in persistent and non-persistent bacterial populations, and this structural agreement motivates the biological reading of the result."],"supporting_citations":[{"why":"Provides the lattice-gas cellular automaton formalism and update rules on which the model is built.","marker":"[59]"},{"why":"Supplies the mean-field Boltzmann propagator and linear stability method used to derive the instability thresholds for flux and nematic modes.","marker":"[66]"},{"why":"Provides the continuum mean-field theory for velocity-alignment particles used for qualitative comparison of nematic and polar modes.","marker":"[22]"},{"why":"Derives cellular automaton models for time-correlated random walks, grounding the persistence Hamiltonian in exponentially correlated noise.","marker":"[61]"},{"why":"Shows how to extract cellular automaton transition rules from Langevin models, supplying the nematic interaction term.","marker":"[65]"},{"why":"Gives the finding, echoed here, that the nematic steady state is destabilized at large persistence, supporting the interpretation of region B.","marker":"[67]"},{"why":"Provides the experimental observations of persistent and non-persistent Myxococcus xanthus phenotypes whose sheet and stream patterns motivate the model.","marker":"[58]"}],"fun_headline_variants":["High persistence foils global nematic order even with strong alignment","Stubborn particles resist collective alignment when both traits are high","Persistence traps particles in local bands, blocking global order","Nematic alignment alone can't overcome individual persistence","Frustrated states emerge when persistence and alignment clash"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim rests on the assumption that the phase diagram measured at one density, $\\rho = 0.2$ on a 120 by 120 periodic lattice, is representative of the whole dilute regime; the density-independence argument given is a linear stability result that covers only infinitesimal perturbations of the homogeneous disordered state, not the nonlinear band and cluster phases at other densities.","fun_headline_variants_meta":{"raw":{"variants":["High persistence foils global nematic order even with strong alignment","Stubborn particles resist collective alignment when both traits are high","Persistence traps particles in local bands, blocking global order","Nematic alignment alone can't overcome individual persistence","Frustrated states emerge when persistence and alignment clash"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00063,"raw_usage":{"total_tokens":2926,"prompt_tokens":977,"completion_tokens":1949,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":1870}},"tokens_in":593,"tokens_out":1949,"duration_ms":13427,"temperature":1.0,"reasoning_tokens":1870,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T13:56:56.069732+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same LGCA on a 120 by 120 lattice at densities $\\rho = 0.1$ and $\\rho = 0.4$ across the $(\\beta, \\alpha)$ parameter plane and measure the nematic order parameter $S_L$ in the stationary state; if at either density the region of high global nematic order extends to large $\\alpha$, or the transition from region A to region B disappears, then persistence-frustrated ordering is not a general dilute-regime property. A cheaper check would be to fix $\\beta = 4$, vary $\\alpha$ at both densities, and compare where $S_L$ drops.","supporting_citations":[{"cited_title":"Deutsch and S","cited_arxiv_id":null,"evidence_quote":"Provides the lattice-gas cellular automaton formalism and update rules on which the model is built."},{"cited_title":"Mean- field analysis of a dynamical phase transition in a cellular automaton model for collective motion,","cited_arxiv_id":null,"evidence_quote":"Supplies the mean-field Boltzmann propagator and linear stability method used to derive the instability thresholds for flux and nematic modes."},{"cited_title":"A mean-field the- ory for self-propelled particles interacting by velocity alignment mechanisms,","cited_arxiv_id":null,"evidence_quote":"Provides the continuum mean-field theory for velocity-alignment particles used for qualitative comparison of nematic and polar modes."},{"cited_title":"Cellular automaton models for time- correlated random walks: derivation and analysis,","cited_arxiv_id":null,"evidence_quote":"Derives cellular automaton models for time-correlated random walks, grounding the persistence Hamiltonian in exponentially correlated noise."},{"cited_title":"Extracting cellular automaton rules from physical langevin equation models for single and col- lective cell migration,","cited_arxiv_id":null,"evidence_quote":"Shows how to extract cellular automaton transition rules from Langevin models, supplying the nematic interaction term."},{"cited_title":"Hydrodynamics of self-propelled hard rods,","cited_arxiv_id":null,"evidence_quote":"Gives the finding, echoed here, that the nematic steady state is destabilized at large persistence, supporting the interpretation of region B."},{"cited_title":"Directional reversals enable myxococ- cus xanthus cells to produce collective one-dimensional streams during fruiting-body formation,","cited_arxiv_id":null,"evidence_quote":"Provides the experimental observations of persistent and non-persistent Myxococcus xanthus phenotypes whose sheet and stream patterns motivate the model."}],"review_version":1}