{"id":"ac3d7063-fd7c-4b35-9a17-589f94830bd8","arxiv_id":"2604.00930","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper argues that the original arctan-based Vicsek update breaks O(2) symmetry and that an adaptive phase shift makes the flocking phase transition vanish.","lead":"The paper claims the original Vicsek model is not O(2)-symmetric because the arctangent angle update has a branch cut, and that an adaptively chosen global phase makes the famous flocking transition disappear. The claim rests on treating the unwrapped angle increment as physical, and the 'vanishing transition' is demonstrated only at one parameter set with a gauge choice that forces the branch cut onto the mean direction.","discovery_kind":"replication","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim collapses if the physical state is the velocity vector on S^1; the 2πn term in Eq. (3.7) is the identity on the circle, so the update is O(2)-equivariant and the adaptive-phase disorder is a branch-cut gauge artifact.","rationale":"The reader's weakest-assumption analysis and my own reading converge on the same issue: the paper's entire case rests on defining O(2) symmetry through the exact invariance of the unwrapped real angle increment Δθ_i. Physical observables in the Vicsek model are built from unit vectors e^{iθ_i}, for which the branch-cut wrapping 2πn(φ) in Eq. (3.7) is the group identity. Consequently, the original arctan-based update is O(2)-equivariant in the correct sense: rotating all input velocities rotates the output velocity distribution covariantly. The 'broken symmetry' identified in Section III is an artifact of lifting S^1 to R and then demanding more than the physical state warrants. The numerical claim that the phase transition vanishes is even weaker: the adaptive phase φ(t)=π−⟨θ_i(t)⟩ is chosen specifically to align the arctangent branch cut with the mean velocity direction, which forces the principal-value mean to jump by 2π and drives v_op to zero. That is a coordinate singularity, not a change in the underlying stochastic process. The paper also presents only a single trajectory at one parameter point, not a phase diagram, so the headline conclusion is unsupported even on its own terms. I find no other load-bearing concern that survives: the algebra in Section III is internally consistent but applied to an unphysical observable; the continuous-time discussion is tangential; and no independent verification (e.g., machine-checked proofs) is offered. The rejection stands, and no verdict adjustment is needed.","tokens_in":6650,"tokens_out":3683,"duration_ms":36969,"concrete_test":"Recompute Fig. 1 (upper) with φ(t)=π−⟨θ_i(t)⟩ but evaluate the order parameter after rotating the output velocities back by −φ(t) (i.e., use R(−φ(t))v_i(t) or equivalently θ_i(t)−φ(t) before wrapping). If the restored v_op(t) matches the φ(t)=0 curve at late times, the disappearance is purely a branch-cut gauge artifact. Additionally, scan η from 0.1 to 2.0 at fixed r_V=0.1 and report v_op(T) with and without the adaptive phase; if the adaptive-phase curve remains near zero for all η, the phase-transition claim would require a full phase diagram, not a single run.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central premise, introduced in Section III as Eq. (3.2), is that O(2) symmetry means exact invariance of the real-valued angle increment Δθ_i(t)=θ_i(t+Δt)-θ_i(t) under the global shift θ_i→θ_i+φ. This is the load-bearing assumption, and it is not the physically appropriate definition. The physical state is the velocity vector v_i(t)=v_abs(cos θ_i, sin θ_i) (Eq. 2.1b); angles are a lift to R, and θ and θ+2π describe the same velocity. In Eq. (3.7), the offending term is 2πn(φ), which is exactly the identity on S^1. The update is equivariant on the circle: arg(⟨e^{i(θ_j+φ)}⟩)=arg(⟨e^{iθ_j}⟩)+φ (mod 2π), and the noise term has the same distribution after rotation. Thus the 'absence of O(2) symmetry' is not a property of the model but of the chosen real-valued coordinate. The numerical demonstration in Fig. 1 (upper) is a direct artifact: φ(t)=π−⟨θ_i(t)⟩ places the principal-value branch cut on the instantaneous mean direction, so the wrapped mean flips to −π and the order parameter collapses regardless of the true dynamics. This is a gauge choice, not a physical disappearance of order. Moreover, the simulation uses a single parameter set (N=1600, L=8.0, η=0.75, r_V=0.1) and shows no phase diagram; it cannot support the statement that 'the phase transition vanishes'. The algebra in Section III is correct but applied to an unphysical observable, and the conclusion does not follow.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the original Vicsek model (Eq. (2.1)) lacks O(2) symmetry because the angular increment Δθ_i(t)=θ_i(t+Δt)-θ_i(t) is not invariant under global phase shifts θ_i→θ_i+φ when the mean direction is computed with the principal-value arctangent. The non-invariance arises through an additional 2πn(φ) term (Eq. (3.7)). The paper then presents a numerical simulation with an adaptive global phase φ(t)=π-⟨θ_i(t)⟩ (Eq. (4.1)) and claims that the order parameter collapses, so the phase transition reported in Vicsek et al. vanishes. It contrasts this with an arithmetic-mean variant (Eq. (2.6)) that is O(2)-symmetric and discusses a continuous-time limit (Section V).","tokens_in":7084,"tokens_out":3765,"duration_ms":35392,"significance":"If the central claim were correct, it would overturn the standard interpretation of the Vicsek transition as spontaneous O(2) symmetry breaking and would imply that the reported phase transition is a gauge artifact. The paper is clearly written, the algebra in Section III is internally consistent, and the distinction between the arctan and arithmetic-mean discretizations is a useful pedagogical point. The continuous-time derivation in Section V is also interesting. However, the load-bearing premise is a nonstandard definition of O(2) symmetry based on invariance of real-valued angular increments of an unwrapped coordinate, whereas the physical state is the velocity vector on S^1. Under the standard equivariant definition, the update is O(2)-equivariant and the adaptive-phase collapse is a branch-cut artifact of the principal-value arctangent. The numerical evidence is also limited to a single time series at one parameter point. The paper therefore does not establish the claimed absence of O(2) symmetry or the vanishing of the Vicsek phase transition.","major_comments":[{"comment":"The definition of O(2) symmetry as exact invariance of the real-valued increment Δθ_i(t) is not the physically appropriate one. The state of particle i is the velocity vector v_i=v_abs(cos θ_i, sin θ_i) (Eq. (2.1b)), so θ_i and θ_i+2π represent the same state. The offending term 2πn(φ) in Eq. (3.7) is the identity on S^1. On the circle the update is equivariant: arg(⟨e^{i(θ_j+φ)}⟩)=arg(⟨e^{iθ_j}⟩)+φ (mod 2π), and the noise distribution is rotationally invariant. Thus the 'absence of O(2) symmetry' is a property of the chosen unwrapped coordinate, not of the model. The central conclusion of the paper rests on this premise.","section":"Sec. III, Eqs. (3.2)-(3.7)"},{"comment":"The numerical demonstration is a direct branch-cut artifact. The adaptive phase φ(t)=π-⟨θ_i(t)⟩ (Eq. (4.1)) places the principal-value branch cut at the instantaneous mean direction, forcing the wrapped mean angle to jump and the order parameter to collapse regardless of the intrinsic dynamics. The paper itself acknowledges this: 'because a branch cut exists near θ=π.' Moreover, the simulation uses a single parameter set (N=1600, L=8.0, η=0.75, r_V=0.1) and shows no phase diagram; a single trajectory cannot support the statement that 'the phase transition reported in Ref. [1] vanishes completely.' At minimum one would need to vary η and r_V and demonstrate the absence of an ordered phase throughout the parameter plane.","section":"Sec. IV, Fig. 1"},{"comment":"The remark about periodic boundary conditions does not address the relevant symmetry. A global phase shift θ_i→θ_i+φ acts on the velocity directions and does not require rotating spatial positions, so the claimed incompatibility with periodic boundary conditions is irrelevant to the O(2)-equivariance of the update rule. If the authors intend to challenge the standard interpretation of the Vicsek transition, they must engage with the standard formulation in which the order parameter is a vector on S^1, rather than redefining symmetry through unwrapped angle increments.","section":"Sec. VI"}],"minor_comments":[{"comment":"The phrase 'the phase transition vanishes completely' is an overclaim relative to the evidence: Fig. 1 shows the time evolution of the order parameter at one parameter point (N=1600, L=8, η=0.75, r_V=0.1), not a phase diagram or a systematic parameter scan.","section":"Abstract and Sec. IV"},{"comment":"The caption uses φ(t)=0, −⟨θ_i(t)⟩, π−⟨θ_i(t)⟩, while the text refers to θ̄=π; please unify the notation. Also, the inset time range t∈[0,200] is not clearly labeled on the inset axes.","section":"Fig. 1"},{"comment":"Using the Frobenius norm ∥·∥_F for a vector is unconventional; the standard Euclidean norm is sufficient and clearer.","section":"Eq. (2.8)"},{"comment":"The continuous-time limit of the arithmetic-mean model is interesting but orthogonal to the main claim about O(2) symmetry. Consider moving it to an appendix or framing it explicitly as an independent contribution.","section":"Sec. V"}],"recommendation":"reject","confidential_remarks":"The paper is internally consistent but the central claim is a coordinate artifact: the physical state space is S^1, and the 2πn term in Eq. (3.7) is the identity on the circle. The numerical demonstration is a branch-cut gauge choice rather than a physical disappearance of order. This is a load-bearing error that cannot be fixed within the manuscript's current scope; the authors would need to reformulate the symmetry definition and provide a genuine parameter-space study to make the claim viable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, my take on arXiv:2604.00930 is that the central claim fails on a definitional point. The authors define O(2) symmetry as invariance of the real-valued angle increment Δθ_i(t) under a global shift θ_i→θ_i+φ, and then show that the principal-value arctangent introduces an extra 2πn(φ) term, which breaks that. The algebra is correct, but the premise is wrong. The physical state is the unit vector v_i=(cosθ_i, sinθ_i); angles are a lift. On the circle, θ and θ+2π are the same, so the 2πn term is the identity. The standard notion of O(2) equivariance for the update is v_i(t+Δt)=v_abs(cos(⟨θ_j⟩+Ξ_i), sin(...)), and rotating all θ_j by φ rotates the average direction by φ mod 2π; the noise distribution is unchanged. So the model is O(2)-equivariant on S^1. Their 'absence of symmetry' is a property of the coordinate, not the model.\n\nWhat is genuinely useful: the paper is a clear reminder that principal-value arctangent has branch cuts, and that wrapped angle updates can be implementation-sensitive. The arithmetic-mean variant, Eq. (2.6), is a sensible alternative, and the continuous-time discussion in Sec. V is fine. The algebra in Sec. III is correct as far as it goes. So the paper is not careless; it is conceptually misframed.\n\nThe soft spots are load-bearing. First, the definition of O(2) symmetry via real-valued increments is nonstandard and unphysical. Second, the adaptive phase φ(t)=π−⟨θ_i(t)⟩ deliberately places the branch cut on the instantaneous mean angle; any angle-based update will show a discontinuity there, and the order parameter collapse in Fig. 1(upper) is a gauge artifact, not a lost transition. Third, the numerical evidence is a single parameter set (N=1600, L=8.0, η=0.75, r_V=0.1) with no phase diagram, so the claim that 'the phase transition vanishes completely' is not supported. The remark about Ref. [15] and periodic boundaries is a side argument and does not rescue the main claim.\n\nWho this is for: someone benchmarking or implementing Vicsek-type models in angle coordinates might find the branch-cut warning useful. But the paper as written does not establish what it claims. I would not cite the central result. I would not send this to peer review as a substantive research paper; a serious editor would desk reject because the premise is unphysical. That said, the authors are earnest, the writing is clear, and there is a kernel of a useful numerical caution here. If they reframed it as a technical note on branch-cut handling and dropped the symmetry-breaking interpretation, it could become a worthwhile, much shorter contribution. My recommendation: decline in current form.","headline":"The paper's central claim is a coordinate artifact—the Vicsek update is O(2)-equivariant on the unit circle, and the adaptive-phase disorder is a gauge choice.","tokens_in":7546,"tokens_out":4255,"would_cite":false,"duration_ms":37331,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The original Vicsek model, with its principal-value arctangent update rule, is not invariant under global phase shifts, and its reported phase transition vanishes when the global phase is chosen adaptively.","keywords":["Vicsek model","O(2) symmetry","phase transition","active matter","circular mean","principal-value arctangent","branch cut","order parameter"],"falsifier":"Repeat the adaptive-phase simulation for the original arctan update with φ(t)=π−⟨θ_i(t)⟩, N=1600, L=8.0, v_abs=0.01, r_V=0.1, η=0.75, and Δt=1.0, and monitor the order parameter v_op(t). If v_op saturates toward 1 rather than staying near 0, the paper's central numerical claim is wrong; if it remains near 0, the claim is confirmed.","tokens_in":6480,"feed_emoji":"🐦","tokens_out":8220,"duration_ms":76617,"temperature":0.7,"pith_summary":"The paper revisits the 1995 Vicsek model and argues that its angle update rule—taking the mean direction as the principal-value arctangent of averaged sine and cosine, equivalently the argument of the averaged complex direction—is not invariant under a global rotation of all angles. At the level of the real-valued angular increment, a phase shift θ→θ+φ adds an extra 2πn(φ) term because the principal-value branch cut forces angles into (−π,π]. The authors conclude that the model lacks O(2) symmetry, counter to the standard story, and they show numerically that adaptively shifting the global phase—keeping the mean angle near π—destroys flocking even at low noise and large interaction radius. By contrast, the arithmetic-mean variant averages angles directly, is O(2)-symmetric, and its ordered phase is insensitive to the global phase. If correct, this would reframe a long-studied nonequilibrium phase transition as at least partly a coordinate artifact.","feed_headline":"Original Vicsek model breaks O(2) symmetry; its transition vanishes","feed_subtitle":"Adaptively choosing the global phase suppresses flocking at low noise, while the arithmetic-mean variant keeps the transition.","key_machinery":"The load-bearing object is the circular mean ⟨θ⟩ = arctan(⟨sin θ⟩/⟨cos θ⟩) = arg(⟨e^{iθ}⟩) used in the update rule. Its principal-value arctangent has a branch cut at ±π; under the global translation θ→θ+φ the increment Δθ_i acquires the extra term 2πn(φ). This branch-cut term is the entire mechanism: it breaks the invariance of the real-valued increment, and it is what the adaptive-phase simulation exploits by moving the branch cut through the flocking state. The arithmetic-mean update, Δθ_i = ⟨θ_j⟩−θ_i, has no branch cut and therefore is exactly invariant.","core_discovery":"The paper's central claim is that the original angle-based Vicsek update is not equivariant under the global transformation θ_i→θ_i+φ when the model is implemented as in the original paper. Because the mean direction is defined as arg(⟨e^{iθ_j}⟩) and computed with the principal value of arctan, a global shift changes the angular increment Δθ_i by 2πn(φ), a jump that is invisible when angles are treated modulo 2π but real when Δθ is read as a real number. The authors take this as a breaking of O(2) symmetry and demonstrate numerically that choosing the global phase adaptively—so that the branch cut sits near the mean angle—suppresses the order parameter entirely. In contrast, replacing the ci","pith_inferences":["A consequence the paper does not draw: because the branch-cut term is invisible to the unit-vector state e^{iθ}, the preference for the real-angle definition is what makes O(2) 'broken'; a reader who takes velocities as the physical state may regard the adaptive-phase simulation as a gauge choice rather than a change of physics.","Published numerical studies of the 'Vicsek model' may inherit the same coordinate dependence; rerunning their order-parameter and finite-size analyses under an adaptive global phase could identify which reported exponents are intrinsic and which are artifacts of the fixed branch cut.","The same diagnostic applies to other models using principal-value circular means, such as synchronization and active-rotation models: adaptively recentering the phase before the arctangent provides a simple test for hidden coordinate dependence."],"forward_implications":["The standard interpretation of the Vicsek transition as spontaneous O(2) symmetry breaking does not hold for the original arctan implementation; the symmetry is broken by the angle coordinate before any flocking order develops.","The order-disorder transition reported for the original model is phase-convention dependent: with an adaptive global phase, the system stays disordered even where low noise and large interaction radius should produce strong alignment.","The arithmetic-mean variant is the O(2)-symmetric version of the model, and its flocking transition is robust to global phase choices, making it a more reliable benchmark for the nonequilibrium transition.","The continuous-time limit, θ̇_i = κ⟨sin(θ_j−θ_i)⟩ + ξ, contains no branch cut, so the natural continuum description of alignment belongs to the arithmetic-mean family, not to the original update."],"fun_headline_variants":["Vicsek model's O(2) breaking destroys its transition","Original Vicsek update shatters O(2), erases phase transition","Global phase adaptation suppresses Vicsek order","O(2) asymmetry in Vicsek model: transition gone","Vicsek model breaks rotational symmetry, loses flocking"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The conclusion rests on treating the unwrapped real-valued angle difference Δθ_i(t) as the physical quantity that must be invariant under θ→θ+φ; if the physical state is instead the unit vector e^{iθ}, the extra 2πn term changes nothing and the original update is equivariant under global rotations.","fun_headline_variants_meta":{"raw":{"variants":["Vicsek model's O(2) breaking destroys its transition","Original Vicsek update shatters O(2), erases phase transition","Global phase adaptation suppresses Vicsek order","O(2) asymmetry in Vicsek model: transition gone","Vicsek model breaks rotational symmetry, loses flocking"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000417,"raw_usage":{"total_tokens":1950,"prompt_tokens":667,"completion_tokens":1283,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":411,"completion_tokens_details":{"reasoning_tokens":1200}},"tokens_in":411,"tokens_out":1283,"duration_ms":9637,"temperature":1.0,"reasoning_tokens":1200,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T16:56:52.170095+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the adaptive-phase simulation for the original arctan update with φ(t)=π−⟨θ_i(t)⟩, N=1600, L=8.0, v_abs=0.01, r_V=0.1, η=0.75, and Δt=1.0, and monitor the order parameter v_op(t). If v_op saturates toward 1 rather than staying near 0, the paper's central numerical claim is wrong; if it remains near 0, the claim is confirmed.","supporting_citations":[],"review_version":1}