{"id":"c9b669af-c225-47b4-9325-a854e0e2cd3c","arxiv_id":"2608.10844","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Triadic alignment interactions in a generalized Vicsek model cause a discontinuous order-disorder transition, unlike the continuous transition from pairwise alignment.","lead":"This paper adds three-particle alignment rules to the classic Vicsek model of flocking and finds that the transition to ordered motion becomes discontinuous instead of smooth. The result suggests that how individuals interact, not just noise or density, can control how abruptly collective order appears.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed discontinuous transition is not established by the reported finite-size data: the susceptibility peak scales as L^0.89, far below the L^2 volume scaling expected for a first-order transition, and Binder-cumulant deepening shows no saturation. Larger-system scaling tests are needed.","rationale":"The reader's weakest assumption was that the finite-size signatures at N up to 9000 indicate a true discontinuous transition in the thermodynamic limit. I agree that this is the load-bearing point, but I sharpen it: the reported χ_max ∼ L^0.89 scaling is not merely an absence of extrapolation; it is quantitatively inconsistent with the L^2 volume scaling required for a first-order transition with homogeneous coexisting phases. This makes the concern more specific than the reader's generic finite-size caveat. The Binder-cumulant minimum deepening linearly with N also lacks the saturation expected for a genuine first-order transition. The mean-field mechanism, in turn, contains a clear formal flaw: it replaces the angularly dependent triadic field magnitude by an averaged positive cubic term, discarding the sign flip that actually distinguishes the triadic rule from the pairwise rule. For these reasons, I do not find the central claim established as written, but the issue can in principle be settled by additional simulations at larger system sizes and by measuring the switching rate between putative coexisting phases. The reader's CONDITIONAL verdict remains appropriate, since the paper can be strengthened by such tests; I therefore do not change the verdict, though I emphasize that the current evidence is weaker than the abstract suggests.","tokens_in":9070,"tokens_out":13835,"duration_ms":146042,"concrete_test":"Run the α=0, ρ=10 HOVM for L=40, 60, and 90 (N=16000, 36000, 81000) with the same protocol, and fit χ_max versus L. For a discontinuous transition, the effective exponent must grow toward 2 (volume scaling); if it remains near 0.9, the first-order interpretation fails. Independently, run at least one very long trajectory (10^6–10^7 steps) at the apparent transition noise and count spontaneous switches between the ordered and disordered branches; if no switch occurs within ≳10^4 τ_max, the double-well free energy and hysteresis reflect metastability rather than equilibrium coexistence.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is that the finite-size signatures in Fig. 4 indicate a true discontinuous transition. This is undermined by the paper's own susceptibility scaling: Fig. 4(a) reports χ_max ∼ L^{0.89±0.03}, whereas a genuine first-order transition with two coexisting homogeneous phases must have χ_max ∼ L^d = L^2, because the variance of the order parameter between the two phases is O(1). The reported exponent is far smaller than volume scaling, suggesting that the simulated ensemble does not actually sample both phases on the measurement time scale. The Appendix autocorrelation analysis measures relaxation within a single phase (τ_max ∼ 10^3 steps), not the barrier-crossing rate between ordered and disordered states, so it cannot certify that the double-well free energy and hysteresis are equilibrium coexistence rather than metastability or initial-condition artifacts. The Binder-cumulant minimum deepening linearly with N without saturation (Fig. 4b) is also not the expected approach to a finite first-order value. Therefore, the central claim that triadic alignment changes the order of the transition is not established by the present data. Additionally, the mean-field derivation replaces ⟨cos(θ_i−Θ)⟩ = v_a inside a positive effective magnitude, erasing the sign-flipping character of the triadic field; the resulting cubic term is therefore an artifact of that replacement and does not correctly represent the mechanism.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a generalized Vicsek model in which a focal particle aligns with the mean velocity of its neighbors (pairwise term) and, additionally or exclusively, with a triadic field built from symmetrized contractions of neighbor velocities with the focal velocity (Eq. 1). A mixing parameter α interpolates between the standard Vicsek model (α=1) and the purely higher-order model (α=0). Agent-based simulations report that the α=0 limit shows an abrupt polarization drop, a sharp susceptibility peak, a negative Binder cumulant minimum, a double-well effective free energy, and hysteresis, which the authors interpret as a discontinuous order-disorder transition; the α=1 limit appears continuous at the same system sizes. A mean-field calculation with F(v_a)=αv_a+(1−α)v_a^3 is offered as the mechanism, and the paper concludes that the microscopic structure of the alignment interaction is an independent control parameter for the order of the transition.","tokens_in":9336,"tokens_out":10225,"duration_ms":95083,"significance":"The central question is interesting and timely: whether genuinely many-body alignment, rather than noise or density, can change the order of the collective-motion transition. If established, the claim would connect active-matter physics with the growing literature on higher-order interactions in complex systems. The paper is also commendable for depositing simulation data, source code, and plotting scripts on Zenodo, and for reporting several complementary observables rather than a single order-parameter curve. However, the quantitative finite-size scaling is not yet consistent with the standard first-order expectation, and the mean-field derivation contains a sign-averaging step that appears to be an artifact. Both issues are load-bearing for the stated conclusion, so the present version does not yet establish the discontinuous-transition claim.","major_comments":[{"comment":"The finite-size evidence is quantitatively inconsistent with a first-order interpretation. For a two-dimensional first-order transition between two coexisting phases, the susceptibility maximum should scale as the volume, χ_max ∼ L^2, because the order-parameter variance between the phases is O(1). The paper's own fit reports χ_max ∼ L^{0.89±0.03}, far below that benchmark. Likewise, the Binder cumulant minimum for a genuine first-order transition should approach a finite value as N→∞; the reported linear deepening of G_min with N (slope ≈ −3.26×10^{-5}) shows no saturation. Because no extrapolation to the thermodynamic limit is performed, the four signatures in Fig. 4 are also compatible with a continuous transition with strong finite-size effects or with metastability. Please provide larger-system data or a scaling collapse, state the expected L^2 benchmark explicitly, and test the N-dependence of the barrier between the two wells.","section":"Section III, Fig. 4(a) and Fig. 4(b)"},{"comment":"The mean-field derivation leading to F(v_a)=αv_a+(1−α)v_a^3 contains an invalid self-consistency step. From Eq. (1), the triadic contribution to the effective field for particle i is proportional to 2v_a^2 cos(θ_i−Θ) e_Θ. For particles with cos(θ_i−Θ)<0, this field points along −e_Θ, i.e., it reinforces the particle's current hemisphere rather than pulling it toward the global direction Θ. Replacing cos(θ_i−Θ) by its average v_a before taking the argument of u_i erases this sign-flipping character and produces a positive effective magnitude that does not represent the actual update rule. The angular distribution is therefore not the assumed von Mises distribution with concentration κ=F(v_a)/η, and the cubic term is an artifact of this substitution. The proposed mechanism for bistability is thus not established by the mean-field calculation; it needs to be derived from the full angular dynamics or presented explicitly as a separate phenomenological model with a justified form.","section":"Section IV, Eq. (1) and self-consistency condition"},{"comment":"The autocorrelation analysis in Appendix A measures relaxation of the global polarization within a single metastable branch, not the rate of barrier crossing between ordered and disordered states. At a genuine first-order transition, the inter-well switching time grows exponentially with N, so t_eq=10^4 with τ_max∼10^3 does not guarantee that the system samples both wells. Consequently, the double-well probability distribution in Fig. 4(c) and the hysteresis branches in Fig. 4(d) could reflect initial-condition dependence or metastability rather than equilibrium coexistence. Please test convergence of P(v_a) from ordered and disordered initial conditions at fixed η and report the inter-well switching rate as a function of N. In addition, the quasi-static sweep in Fig. 4(d) is not tested against sweep rate; hysteresis alone is a dynamical signature and should be shown to vanish in the appropriate limit.","section":"Appendix A, Fig. 4(c) and Fig. 4(d)"}],"minor_comments":[{"comment":"There is a typo: 'a discontinues phase transition' should read 'a discontinuous phase transition'.","section":"Conclusion"},{"comment":"The linear fit to G_min versus N is reported without error bars on the points or an uncertainty on the slope; please provide these to support the claimed linear deepening.","section":"Fig. 4(b) inset"},{"comment":"Reference [38] is a closely related higher-order Vicsek model, but the paper does not explain how the present triadic rule differs from it; please add a brief comparison.","section":"References, Ref. [38]"},{"comment":"The factor 2 arising from the two symmetrized contractions in Eq. (1) is not carried through to the mean-field expression F(v_a); this is immaterial for the bistability argument but should be stated explicitly to avoid confusion.","section":"Section II, Eq. (1)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The model is genuinely new: triadic alignment built from the symmetric tensor product of neighbor velocities contracted with the focal particle's velocity. That's a clean, non-pairwise rule, and the paper studies it with a reasonable protocol: fixed equilibration and sampling windows tied to measured autocorrelation times, 10^3 realizations, and several independent observables. Hysteresis, double-well free energy, Binder cumulant dip, and a growing susceptibility peak all point the same direction. I believe the authors when they say the triadic limit behaves differently from the pairwise one at these system sizes.\n\nThe soft spots are real but not all equally bad. The most obvious: reference [38] (León et al., 2025) already introduces a higher-order Vicsek model, and this paper never compares itself to it. That's a gap a referee would want closed before publication. The mean-field derivation is also sloppy—replacing cos(θ_i−Θ) by va inside the magnitude of the triadic field erases the sign dependence that is the actual mechanism, and the resulting cubic term is more an analogy than a derivation. The authors call it phenomenological, so I don't count it as a fatal flaw, but it shouldn't be presented as the mechanism without qualification.\n\nThe deeper worry is the thermodynamics. A true first-order transition in 2D with two coexisting homogeneous phases should show χ_max ~ L^2; the paper reports L^0.89. That's a large mismatch. The likely explanation is that the simulations don't sample both phases on the measurement time scale, and the hysteresis then reflects metastability, not equilibrium coexistence. The appendix's autocorrelation analysis measures relaxation within a phase, not barrier crossing between phases, so it doesn't resolve this. The Binder minimum deepening linearly with N without saturation is also not the standard approach to a finite first-order value. These issues mean the central claim—that triadic alignment changes the order of the transition—is plausible but not established. It needs larger system sizes and, ideally, a check for whether the coexisting phases are homogeneous or band-like, where the scaling expectations differ.\n\nThis paper deserves a serious referee. The model is interesting enough, and the simulations are clean enough, that a careful revision could make it solid. I'd send it to review, and I'd tell the authors to deal with [38], fix the MF, and either get larger-scale scaling data or soften their claim to 'first-order-like at accessible sizes.'\n\nFor a reading group: maybe, if your group likes active matter and higher-order interactions. I wouldn't cite it yet in my own work, because the core claim is still under-supported.","headline":"Plausible and worth a serious look, but the discontinuity claim rests on finite-size signatures that don't yet scale like a true first-order transition, and the paper misses a directly competing 2025 model.","tokens_in":9843,"tokens_out":3616,"would_cite":false,"duration_ms":34012,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that replacing pairwise alignment with a genuinely triadic alignment rule in the Vicsek model changes the order–disorder transition from continuous to discontinuous, making the microscopic interaction structure itself an…","keywords":["collective motion","Vicsek model","higher-order interactions","triadic alignment","discontinuous phase transition","active matter","Binder cumulant","mean-field theory"],"falsifier":"Extrapolate the deepest Binder-cumulant value $G_{\\min}$ to $N\\to\\infty$ for $\\alpha=0$ at $\\rho=10$: a genuine first-order transition requires $G_{\\min}$ to remain negative or keep decreasing, while a return toward zero would mark the dip as a finite-size artifact. A second check is to repeat the quasi-static noise sweep with an equilibration window much longer than the measured maximum autocorrelation time $\\tau_{\\max}$; if the forward and backward branches coincide, the apparent hysteresis is not true metastability.","tokens_in":8865,"feed_emoji":"🐦","tokens_out":13529,"duration_ms":112345,"temperature":0.7,"pith_summary":"This paper asks whether the order of the collective-motion phase transition can be changed by the structure of the alignment interaction alone, without tuning noise or density. It introduces a generalized Vicsek model with a single parameter $\\alpha$ that interpolates between standard pairwise alignment and a genuinely triadic rule coupling each particle to pairs of neighbors. Simulations show that the pure triadic limit produces a discontinuous order-disorder transition: hysteresis, a double-well free-energy landscape, and a Binder cumulant minimum that deepens with system size, while the standard pairwise model stays continuous at the same sizes. A mean-field analysis traces the discontinuity to a cubic term in the effective alignment strength. The paper concludes that microscopic interaction structure is an independent control parameter for the order of the phase transition in active matter.","feed_headline":"Triadic alignment makes Vicsek order-disorder transition discontinuous","feed_subtitle":"Triadic interactions bring hysteresis and phase coexistence; pairwise interaction remains continuous.","key_machinery":"The load-bearing object is the triadic alignment field of Eq. (1), $h_i(t)=(v_0^3|T_i(t)|)^{-1}\\sum_{(j,k)\\in T_i(t)}[(v_j\\cdot v_i)v_k+(v_k\\cdot v_i)v_j]$, which couples each focal particle $i$ to unordered pairs of neighbors inside its radius; each term involves three particles simultaneously and cannot be written as a sum of pairwise couplings. The hybrid model mixes this with the standard pairwise polarization through the effective field $u_i=\\alpha p_i+(1-\\alpha)h_i$. In the mean-field reduction, replacing neighbor velocities by the global polarization makes the pairwise field contribute linearly in the order parameter $v_a$ and the triadic field contribute a cubic term, giving effective alignment strength $F(v_a)=\\alpha v_a+(1-\\alpha)v_a^3$. With von Mises noise the self-consistency condition is $v_a=I_1(F(v_a)/\\eta)/I_0(F(v_a)/\\eta)$; the cubic term creates two stable branches separated by an unstable branch, and that bistability is the mechanism behind the discontinuous transition.","core_discovery":"The paper's central claim is that a genuinely triadic alignment interaction, in which each particle aligns with the joint effect of pairs of neighbors, changes the nature of the Vicsek order-disorder transition from continuous to discontinuous. In agent-based simulations at density $\\rho=10$ and system sizes up to $N=9000$, the pure higher-order model exhibits four independent signatures of a first-order transition: a susceptibility peak that grows and narrows with system size, a Binder cumulant minimum that becomes more negative with $N$, a double-well effective free energy near critical noise, and hysteresis under quasi-static noise sweeps. The standard pairwise model at the same sizes shows none of these. A mean-field calculation with von Mises angular noise yields a self-consistency condition whose effective alignment strength is $F(v_a)=\\alpha v_a+(1-\\alpha)v_a^3$; the cubic term produces two stable fixed-point branches separated by an unstable one in the $\\alpha=0$ limit, whereas $\\alpha=1$ has a single continuously vanishing branch. The authors conclude that the microscopic structure of the alignment interaction is an independent control parameter for the order of the phase transition in active matter.","pith_inferences":["Extension: the cubic self-consistency mechanism suggests that other triadic constructions will also produce bistability, so the discontinuity may be generic for many-body alignment rules rather than specific to Eq. (1).","Extension: in biological flocks where three-body correlations are measurable, polarization distributions should be bimodal near the transition if triadic interactions dominate; field data showing hysteresis would distinguish pairwise from higher-order control.","Extension: fixing $\\alpha=0$ and lowering density toward $\\rho=1$ should weaken or erase the discontinuity, separating the role of coordination number from the pure shape of the interaction.","Extension: repeating the finite-size analysis with larger $N$ and with Gaussian noise instead of uniform noise would show whether the discontinuity is tied to the specific noise distribution, since the mean-field uses von Mises statistics while the simulations use uniform noise."],"forward_implications":["Pure triadic alignment in the Vicsek model must show a first-order order-disorder transition, with phase coexistence, at densities where pairwise alignment is still continuous.","Collective order under triadic rules is harder to sustain because higher densities are needed, and once lost it breaks abruptly, so such systems would display sudden order collapse under small noise increases.","Mixing a small amount of pairwise alignment ($\\alpha>0$ but small) should broaden the transition and shift it to higher noise, so the transition order is tunable by the interaction structure itself.","The mean-field mechanism predicts that any local alignment rule whose effective strength is cubic in the order parameter will produce bistability, giving a general criterion for designing first-order active-matter transitions."],"supporting_citations":[{"why":"Defines the standard Vicsek model, the pairwise alignment rule, and the update protocol that the generalized model builds on and compares against.","marker":"[11]"},{"why":"Supplies the von Mises distribution whose moment ratio $I_1/I_0$ enters the mean-field self-consistency condition.","marker":"[40]"},{"why":"Provides the reference treatment of order-disorder transitions and autocorrelation decay in Vicsek-class models, used for the sampling protocol and transition characterization.","marker":"[18]"},{"why":"Establishes that higher-order interactions produce abrupt transitions and multistability in complex systems, the conceptual setting the paper extends to active matter.","marker":"[30]"},{"why":"Reviews irreducible higher-order couplings in networks, supporting the claim that triadic alignment cannot be decomposed into pairwise terms.","marker":"[31]"},{"why":"Presents a related higher-order Vicsek model with conformity interactions, giving the nearest prior construction that the present triadic rule generalizes or contrasts with.","marker":"[38]"},{"why":"Supplies the integrated-autocorrelation-time framework used to choose equilibration and sampling windows for the steady-state and hysteresis measurements.","marker":"[41]"}],"fun_headline_variants":["Triadic alignment flips Vicsek transition to first-order","Three-body alignment makes Vicsek transition discontinuous","Triadic interactions turn Vicsek order-disorder into a jump","Pairwise stays continuous, triadic jumps: interaction order matters","Triadic Vicsek model: first-order order-disorder transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on treating the finite-size signatures in Figure 4—deepening Binder minimum, double-well free energy, hysteresis at $N\\le 9000$—as evidence of a genuine first-order transition in the thermodynamic limit, and on trusting the explicitly phenomenological mean-field model that reproduces them only qualitatively.","fun_headline_variants_meta":{"raw":{"variants":["Triadic alignment flips Vicsek transition to first-order","Three-body alignment makes Vicsek transition discontinuous","Triadic interactions turn Vicsek order-disorder into a jump","Pairwise stays continuous, triadic jumps: interaction order matters","Triadic Vicsek model: first-order order-disorder transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001106,"raw_usage":{"total_tokens":4598,"prompt_tokens":921,"completion_tokens":3677,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":537,"completion_tokens_details":{"reasoning_tokens":3596}},"tokens_in":537,"tokens_out":3677,"duration_ms":27444,"temperature":1.0,"reasoning_tokens":3596,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:56:20.639991+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Extrapolate the deepest Binder-cumulant value $G_{\\min}$ to $N\\to\\infty$ for $\\alpha=0$ at $\\rho=10$: a genuine first-order transition requires $G_{\\min}$ to remain negative or keep decreasing, while a return toward zero would mark the dip as a finite-size artifact. A second check is to repeat the quasi-static noise sweep with an equilibration window much longer than the measured maximum autocorrelation time $\\tau_{\\max}$; if the forward and backward branches coincide, the apparent hysteresis is not true metastability.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the von Mises distribution whose moment ratio $I_1/I_0$ enters the mean-field self-consistency condition."},{"cited_title":"Chat´ e, F","cited_arxiv_id":null,"evidence_quote":"Provides the reference treatment of order-disorder transitions and autocorrelation decay in Vicsek-class models, used for the sampling protocol and transition characterization."},{"cited_title":"Boccaletti, P","cited_arxiv_id":null,"evidence_quote":"Reviews irreducible higher-order couplings in networks, supporting the claim that triadic alignment cannot be decomposed into pairwise terms."},{"cited_title":"Le´ on, R","cited_arxiv_id":"2512.19318","evidence_quote":"Presents a related higher-order Vicsek model with conformity interactions, giving the nearest prior construction that the present triadic rule generalizes or contrasts with."},{"cited_title":"Sokal, inFunctional Integration, NATO ASI Series, Vol","cited_arxiv_id":null,"evidence_quote":"Supplies the integrated-autocorrelation-time framework used to choose equilibration and sampling windows for the steady-state and hysteresis measurements."}],"review_version":1}