{"id":"70f3d5b4-c62a-472e-8b7d-98ad2d35bcc5","arxiv_id":"2506.12967","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A microscopic Landau-Ginzburg theory predicts the self-alignment strength at which an active crystal transitions from disordered motion to collective flocking, with a diverging correlation length.","lead":"This paper derives a simple formula that predicts when a dense crystal of self-propelled particles starts to move together in one direction, a phenomenon called flocking. It shows that this transition behaves like a classic second-order phase transition, the same family as magnetism.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Second-order claim rests on a rigid-lattice LG functional and finite-size simulations; no finite-size scaling or hysteresis check rules out density-coupling or fluctuation-induced first-order behavior.","rationale":"I read the paper's main claim as: for self-aligning active crystals, flocking is a second-order phase transition, and the microscopic theory predicts Bc=Pe^-1 with a diverging correlation length. The derivation of the transition point is clean, and the agreement of Bc=Pe^-1 with the phase diagram is a genuine success. However, the headline 'second-order' rests on two pillars: the rigid-lattice LG free energy and finite-size simulations of continuous order-parameter growth. The weakest link is not the algebra but the uncontrolled truncation of density and positional fluctuations, plus the absence of the standard checks (hysteresis, finite-size scaling) that distinguish a continuous transition from a weakly first-order one. The reader's weakest assumption identified the rigid-lattice approximation; I agree that this is central, but I would sharpen the test to the dynamical signatures of transition order. This supports the CONDITIONAL verdict without changing it.","tokens_in":9689,"tokens_out":9223,"duration_ms":120379,"concrete_test":"Run quasistatic B sweeps at Pe=10 for system sizes L=32, 64, and 128, ramping B upward and then downward through Bc=Pe^-1=0.1, and measure the polarization S(B) and susceptibility chi=L^2(<S^2>-<S>^2). If hysteresis is absent and Bc(L) extrapolates to Pe^-1 with a consistent finite-size scaling collapse, the second-order LG scenario survives; if hysteresis appears or Bc(L) drifts with L, the true transition is first-order or fluctuation-dominated. As a secondary check, repeat at a lower packing fraction such as Phi=0.9 to expose density-velocity coupling: a large shift in Bc or a discontinuous S would show that the F_i=0 rigid-lattice approximation is not inert.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—a second-order flocking transition with Bc=Pe^-1, |S| ~ sqrt(B-Bc), and correlation length lambda ~ (1-BPe)^-1/2—is obtained from the Landau-Ginzburg functional in Eq. (4), which in turn is derived under the rigid-lattice approximation F_i=0 in Eq. (2). This approximation removes positional and density degrees of freedom from the theory. The actual dynamics in Eq. (1) include density fluctuations, phonons, and possible lattice defects; if any of these couple to the local velocity field, the effective free energy is not simply Eq. (4), and the transition point or the order of the transition can change. The numerical evidence in Fig. 2(d) and Fig. 3(b) is consistent with a continuous transition, but there is no finite-size scaling analysis, no hysteresis sweep, and no error bars on the correlation length, so a weakly first-order transition or a system-size-dependent Bc is not excluded. The theory is mean-field-like in the sense that the sign change of the mass term guarantees a Mexican-hat potential in the homogeneous limit, but fluctuations in two dimensions can renormalize or even alter the predicted critical behavior. In addition, Eq. (7) contains an unstated lattice stiffness K, so the correlation-length divergence is not an independent quantitative prediction. If density fluctuations are irrelevant, the theory is a valuable exact mapping for the crystal phase; if they are relevant, the second-order claim is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a two-dimensional crystal of self-aligning active Brownian particles and proposes a microscopic mapping from the particle dynamics to a Landau-Ginzburg free energy for the velocity field. Under a rigid-lattice approximation (F_i=0), the authors derive an effective velocity evolution equation, Eq. (2), and from it the free-energy functional in Eq. (4), whose mass term changes sign at β_c/γ_r = 1/(v0 τ), i.e., B_c = Pe^{-1}. From this free energy they predict a mean-field polarization curve |S| ~ sqrt(B−B_c) and a correlation length λ ~ (1−BPe)^{-1/2}. Simulations of the full particle model show a smooth increase of the polarization and a growing correlation length, which the authors interpret as evidence that flocking in these crystals is a second-order phase transition. The paper claims to provide the first microscopic theory for self-alignment-induced flocking in dense active matter.","tokens_in":9951,"tokens_out":5799,"duration_ms":64737,"significance":"If the mapping is correct, the paper delivers a valuable and nontrivial result: a parameter-free prediction of the transition point, B_c = Pe^{-1}, that appears to match simulations, and a transparent connection between a microscopic active-particle model and a Landau-Ginzburg description. The derivation is elegant and the simulations are relevant to recent experimental and numerical observations of self-aligning granular particles and cell monolayers. The paper also makes an explicit and falsifiable prediction for the divergence of the correlation length. These strengths are substantial. However, the current manuscript does not yet provide all the evidence needed to support the strong claim that the transition is second order: the derivation is not self-contained, the correlation-length amplitude contains an unspecified stiffness K, and the numerical evidence lacks finite-size scaling and other checks that would rule out a weakly first-order transition.","major_comments":[{"comment":"The derivation of the effective velocity equation and of the Landau-Ginzburg free energy is deferred entirely to the Supplemental Material, which is not available with the manuscript. Since every subsequent prediction follows from Eq. (4), the main text should present the key steps of the derivation (or the SM should be provided for review), so that the mapping can be checked. This is a load-bearing point, not a presentation issue.","section":"Derivation of Eqs. (2) and (4)"},{"comment":"The correlation-length prediction in Eq. (7) contains the lattice stiffness K, whose numerical value is not given. The comparison in Fig. 3(b) therefore cannot be assessed as a parameter-free quantitative test. The authors should either compute K from the WCA potential and the lattice geometry, or state explicitly that K is fitted to the simulation data. This matters because the divergence of the correlation length is one of the two central pieces of evidence for the second-order claim.","section":"Eq. (7) and Fig. 3(b)"},{"comment":"The phase boundary in Fig. 2(a) is defined by the condition ⟨S⟩>0.5, which is an arbitrary threshold. For a continuous transition, the critical point is the onset of nonzero polarization in the thermodynamic limit, not the location of the 0.5 contour. Comparing the theoretical line B_c = Pe^{-1} to this threshold contour is not a rigorous test of Eq. (5); the authors should provide a finite-size extrapolation of the critical point, for example from Binder cumulants or from the crossing of the correlation length at different system sizes.","section":"Fig. 2(a) and Eq. (5)"},{"comment":"The numerical evidence for a second-order transition is incomplete. No finite-size scaling analysis, no hysteresis sweep, and no error bars on the correlation length are presented. A weakly first-order transition can produce a smooth-looking order-parameter curve and a growing correlation length in a finite box. To substantiate the classification as second order, the authors should examine the system-size dependence of ⟨S⟩ and ξ and check for discontinuities in the order parameter or in the Binder parameter across the transition.","section":"Sec. on spatial velocity correlations and Fig. 3"},{"comment":"The derivation of Eq. (2) assumes that particles sit at fixed lattice positions with balanced forces, F_i = 0. This rigid-lattice approximation removes density fluctuations, phonons, and positional disorder from the theory. Near the transition, if density or positional fluctuations couple to the velocity field, the effective free energy may differ from Eq. (4) and the transition point or its order could change. The manuscript should justify the irrelevance of these couplings in the crystal phase, or at least state this as an explicit limitation of the theory.","section":"Lattice approximation in Eq. (2)"}],"minor_comments":[{"comment":"There is a typo in the caption: 'self-alignment strenght' should be 'self-alignment strength'.","section":"Fig. 1 caption"},{"comment":"In the caption, 'P` eclet' should be 'Péclet' (the accent is misplaced).","section":"Fig. 2 caption"},{"comment":"The expression after the proportionality sign is written as sqrt(B−B_c/B), which is ambiguous; it should be sqrt((B−B_c)/B) to match the preceding formula.","section":"Eq. (6)"},{"comment":"In the caption for Fig. 3(b), the reference to Eq. (7) appears after the description of the critical points; it would be clearer to state explicitly whether the solid curves are the prediction of Eq. (7) with a specified K or a fit to the data.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper presents an appealing and potentially important microscopic mapping, and the transition-point prediction is attractive. However, the strength of the evidence does not yet justify the strong second-order claim. The derivation should be made available in the review process, the correlation-length amplitude should be either derived or clearly identified as a fit, and the numerical evidence should include finite-size scaling or at least a discussion of system-size dependence. I believe these issues are addressable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper's central result is a microscopic Landau-Ginzburg description for self-aligning active crystals, with a parameter-free prediction for the transition point, B_c = Pe^{-1}, and a continuous polarization curve |S| ~ sqrt(B-B_c). That part is genuinely new and likely correct. The second thing: the derivation leans on a rigid-lattice approximation (F_i=0) and the interaction curvature K is never specified in the main text, so the correlation-length prediction in Eq. (7) is not fully independent.\n\nWhat the paper does well. The self-alignment torque creates a sign-changing mass term in the effective velocity free energy; that's a real extension of the Marconi-Maggi/Fodor mappings. The transition point follows without fitting, and the simulations show the polarization rising continuously and the correlation length growing as B approaches Pe^{-1}. The comparison in Fig. 2(d) is honest about the additive constant. This is the best available theory for this system.\n\nWhere it's soft. The lattice assumption F_i=0 is load-bearing: it removes density and positional fluctuations from the theory. If those couple to the velocity field, the transition could be renormalized or even driven first-order. The simulations are consistent with a continuous transition, but there is no finite-size scaling analysis and no hysteresis sweep, so a weakly first-order transition is not excluded. Also, K is not given, so the amplitude of the divergence is not a prediction; the fit parameters a and xi in Fig. 3 effectively absorb that. The S>0.5 threshold for the transition line is arbitrary, though it does seem to land at B_c. All of these are fixable or at least quantifiable. The Supplemental Material likely contains the derivation and the K value; we only have the main text, so I can't verify that part.\n\nVerdict. This is a solid, useful paper for the active matter community, especially people working on granular polar disks and cell monolayers. The transition point prediction is concrete and testable. The second-order claim is plausible but not proven to the level the abstract suggests. I would send it to peer review: the referees should ask for the derivation details, the numerical value or fitting of K, and some finite-size analysis. The paper may need heavy revision, but it deserves referee time.","headline":"A plausible, parameter-free theory for the flocking transition in self-aligning active crystals; the central prediction B_c=Pe^{-1} is clean, but the rigid-lattice assumption and an unspecified stiffness K mean the second-order claim is not as settled as the abstract suggests.","tokens_in":10534,"tokens_out":2775,"would_cite":true,"duration_ms":30242,"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":"A self-aligning active crystal flocks through a continuous, second-order phase transition, and a Landau-Ginzburg free energy predicts where the order appears.","keywords":["active matter","flocking","self-alignment","Landau-Ginzburg theory","second-order phase transition","active crystals","velocity correlations","Péclet number"],"falsifier":"Measure the polarization and correlation length as functions of $B$ in a self-aligning crystal that is allowed to deform or host vacancies at a fixed Péclet number; a discontinuous jump in polarization, or a critical point displaced from $B_c = \\mathrm{Pe}^{-1}$ beyond simulation uncertainty, would show the fixed-lattice free energy does not capture the actual transition.","tokens_in":9451,"feed_emoji":"🐦","tokens_out":6150,"duration_ms":68224,"temperature":0.7,"pith_summary":"The paper claims that a dense two-dimensional crystal of particles whose orientation self-aligns to their velocity undergoes a genuine phase transition from a disordered to a flocking state, and that this transition is continuous, i.e. second order. A microscopic derivation maps the particle dynamics onto a Landau-Ginzburg free energy for the velocity field, in which the sign of the quadratic mass term flips at the critical self-alignment strength $\\beta_c/\\gamma_r = 1/(v_0 \\tau)$. This yields an explicit prediction for the transition point in terms of the Péclet number, $B_c = \\mathrm{Pe}^{-1}$, a mean-field polarization curve $|S| \\sim \\sqrt{B-B_c}$, and a correlation length that diverges as $(1-B/B_c)^{-1/2}$. If the mapping is correct, it provides a microscopic foundation for reading experiments on granular self-aligning particles and migrating cells as critical phenomena, and it distinguishes this mechanism from Vicsek-style velocity-alignment models.","feed_headline":"Flocking in active crystals is a second-order phase transition","feed_subtitle":"A Landau-Ginzburg free energy predicts the exact threshold and the diverging velocity correlations.","key_machinery":"The load-bearing object is a Landau-Ginzburg free-energy functional for the velocity field, Eq. (4): $F[v] = (3\\tau\\sigma^2 K/2\\gamma)(\\nabla v)^2 + (1 - v_0\\tau\\beta/\\gamma_r)|v|^2/2 + (\\tau\\beta/\\gamma_r)|v|^4/(4 v_0^2)$. The parameter $K$ is the curvature (second derivative) of the interaction potential on the lattice, and it turns the repulsive force between nearest neighbors into a discrete Laplacian coupling velocities. The sign of the mass term is controlled by the competition between the self-alignment time $\\gamma_r/\\beta$ and the persistence length $v_0 \\tau$; when the mass term goes negative, the potential becomes a Mexican hat, rotational symmetry breaks, and flocking sets in.","core_discovery":"The central discovery is that the velocity field of a self-aligning active crystal obeys a relaxational Landau-Ginzburg dynamics, Eq. (3), with the effective free energy of Eq. (4). The coefficient of $|v|^2$ changes sign exactly when $\\beta/\\gamma_r$ crosses $1/(v_0\\tau)$, so the free energy goes from a single well to a Mexican hat. The gradient term, produced by repulsive interactions acting as a discrete Laplacian on velocities, penalizes spatial variation and sets the correlation length. Simulations corroborate the predicted transition point, the square-root growth of polarization, and the divergence of the correlation length, which together establish that flocking in this system is a second-order phase transition.","pith_inferences":["If density fluctuations and lattice disorder are allowed to couple to the velocity field, the fixed-lattice prediction for the critical point might acquire corrections; a natural extension is to test whether $B_c = \\mathrm{Pe}^{-1}$ survives in deformable or defective crystals.","The mean-field square-root polarization and the $1/2$ correlation-length exponent suggest that flocking in this class may belong to a different universality class from the continuous but non-mean-field transitions of Vicsek-style models, a question a full renormalization-group analysis could settle.","The theory implies that velocity fluctuations should relax slowly near the transition; measuring the dynamic correlation time as a function of $B - B_c$ on granular particles would be a direct experimental test of the Landau-Ginzburg dynamics.","Because the crystal is used as a rigid background, applying the same mapping to active liquids or glasses, where positional order is weaker, may reveal whether the transition there remains second order or becomes first order."],"forward_implications":["The transition point is fixed by the dimensionless combination $B = \\beta\\sigma/\\gamma_r$ equaling $\\mathrm{Pe}^{-1}$, so increasing the persistence length lowers the self-alignment strength needed for flocking.","The polarization grows continuously as $|S| \\sim \\sqrt{B-B_c}$, a mean-field behavior reminiscent of the magnetization curve in the Ising model.","The correlation length of spatial velocity correlations diverges as $(1 - B/B_c)^{-1/2}$, making the system scale-free as the transition is approached.","No explicit velocity-alignment interaction between neighbors is required for flocking; repulsive interactions plus self-alignment suffice, which explains collective motion in polar granular particles and cell monolayers.","The free-energy structure provides a basis for building a hydrodynamic (Toner-Tu-type) description of self-aligning active matter."],"supporting_citations":[{"why":"Supplies the time-derivative mapping from position and orientation dynamics to effective velocity equations used to obtain Eq. (2).","marker":"[63, 64]"},{"why":"Provides the connected velocity-correlation analysis used to characterize the scale-free properties of the flocking state.","marker":"[69]"},{"why":"Gives experimental evidence that self-alignment induces a flocking transition in polar granular systems, the phenomenon the theory targets.","marker":"[52]"},{"why":"Provides the preceding numerical observation of a disordered-to-flocking transition in dense self-aligning systems that the theory explains.","marker":"[58]"},{"why":"Supplies the Ising and XY model analogy and the mean-field magnetization curve used to interpret the polarization behavior.","marker":"[29]"},{"why":"Provides the hydrodynamic description that the authors propose the free-energy theory can justify for self-aligning active matter.","marker":"[79]"}],"fun_headline_variants":["Active crystals flock via continuous transition","Self-aligning crystals show true phase transition","Flocking threshold predicted by Landau-Ginzburg","Diverging correlations signal flocking onset","Crystal flocking: a second-order phase transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation treats the crystal as a rigid, force-balanced lattice whose only slow variable is the velocity field, so density fluctuations, positional disorder, and the unspecified interaction curvature $K$ play no role.","fun_headline_variants_meta":{"raw":{"variants":["Active crystals flock via continuous transition","Self-aligning crystals show true phase transition","Flocking threshold predicted by Landau-Ginzburg","Diverging correlations signal flocking onset","Crystal flocking: a second-order phase transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000127,"raw_usage":{"total_tokens":1054,"prompt_tokens":826,"completion_tokens":228,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":442,"completion_tokens_details":{"reasoning_tokens":160}},"tokens_in":442,"tokens_out":228,"duration_ms":3235,"temperature":1.0,"reasoning_tokens":160,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:36:41.611163+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the polarization and correlation length as functions of $B$ in a self-aligning crystal that is allowed to deform or host vacancies at a fixed Péclet number; a discontinuous jump in polarization, or a critical point displaced from $B_c = \\mathrm{Pe}^{-1}$ beyond simulation uncertainty, would show the fixed-lattice free energy does not capture the actual transition.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives experimental evidence that self-alignment induces a flocking transition in polar granular systems, the phenomenon the theory targets."},{"cited_title":"Paoluzzi, D","cited_arxiv_id":null,"evidence_quote":"Provides the preceding numerical observation of a disordered-to-flocking transition in dense self-aligning systems that the theory explains."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Ising and XY model analogy and the mean-field magnetization curve used to interpret the polarization behavior."},{"cited_title":"Toner, Y","cited_arxiv_id":null,"evidence_quote":"Provides the hydrodynamic description that the authors propose the free-energy theory can justify for self-aligning active matter."}],"review_version":1}