{"id":"70354090-be14-4687-b027-9f23d019363d","arxiv_id":"2509.02386","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Self-propulsion in an active nematic model produces a non-monotonic ordering effect, with order, defect anti-hyperuniformity, and long-range vorticity correlations all peaking at an intermediate self-propulsion speed.","lead":"The paper adds a self-propulsion term to the standard equations for active nematic fluids and shows that an intermediate self-propulsion speed creates a special regime: stronger nematic order, defect positions that fluctuate wildly over long distances, and long-lived vorticity patterns. A generalist may care because nematically ordered, self-propelling systems such as cell monolayers and bacterial swarms are everywhere in biology, and this work identifies a single control kno","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"All V0-dependent results rest on the sharp, non-dynamical closure p = sign(n̂·v)n̂; the claimed optimum and scaling exponents are not shown to survive any regularization of this rule.","rationale":"The reader's weakest-assumption analysis identified the polarity assignment rule as the key vulnerability, largely on the grounds that real self-propelled particles carry an independent polar field. My concern sharpens this: the issue is not only external validity but also internal robustness. Because p is slaved to v through a discontinuous sign rule, the model is not a generic 'self-propulsive active nematic' but a specific, singular closure. The linear stability threshold (Eq. 9) is internally consistent—I re-derived the eigenvalue and threshold condition, and the algebra checks out—so the analytic core is not the weak point. The weak point is that the nonlinear results, which carry the paper's headline claims, are obtained with a single non-analytic polarity rule and no test of whether those results survive smoothing or a dynamical p. This is exactly the kind of assumption that can change exponents and phase boundaries in active matter models. I therefore do not propose rejection: the paper is a legitimate minimal model, and the authors state their assumption clearly. But the central claims should be presented as conditional on this closure until a robustness check is provided. Since the reader already returned CONDITIONAL, my analysis does not move the verdict; it reinforces the conditionality. The concrete test above is the minimal simulation sweep that would settle whether the concern lands.","tokens_in":11980,"tokens_out":18924,"duration_ms":226900,"concrete_test":"Keep Eqs. 1–6 and the numerical setup identical, but replace the hard assignment p = sign(n̂·v)n̂ with the smoothed rule p = n̂ tanh(β n̂·v) for β ∈ {1, 10, 100} (and compare with β → ∞). At ζ = 0.05, L = 1024, recompute the V0-dependence of the correlation length l_c (Fig. 3e), the structure-factor exponent α (Fig. 4b), and the vorticity exponent ν (Fig. 5b). If the peak remains at V0 ≈ 0.06 and α, ν converge to within the fitting uncertainty as β increases, the sharp flow-alignment closure is not the source of the claimed effects. If the peak shifts or α, ν vary systematically with β, the reported self-propulsion phenomenology is an artifact of the discontinuous closure, and the conclusions must be restricted to that specific rule.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The model's only polar ingredient is the advective term V0 p·∇Q, with p assigned instantaneously as the end of n̂ least deviating from the local flow (Methods/Fig. 1), i.e. p = sign(n̂·v)n̂. This is a discontinuous, non-dynamical closure: at n̂·v = 0 the assignment is undefined, and no equation of motion, persistence, or noise for p is introduced. All of the central nonlinear results—the order/correlation-length peak at V0 ≈ 0.06, defect anti-hyperuniformity with α ≈ −0.9, and vorticity decay with ν < 2—are therefore properties of this specific flow-alignment rule, not of an independent self-propulsion degree of freedom. If the rule is regularized (e.g. p = n̂ tanh(β n̂·v)) or if p is allowed to relax with its own dynamics, the optimum location and the extracted exponents can shift or disappear. The cited experiments [27–29] demonstrate flow alignment in particular biological settings, but they do not establish that intrinsic polarity is slaved to the local velocity on all scales, especially near topological defects where the sign of n̂·v changes. The paper reports no sensitivity analysis with respect to this closure, so the quantitative predictions are not yet anchored to the physical mechanism they claim to represent. This is a load-bearing modelling assumption, not merely an interpretation issue, because every headline claim in the abstract is conditioned on it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a minimal extension of the standard active nematic model by adding a self-advective term V0 p_k ∂_k Q_ij to the Beris–Edwards equation, with polarity p assigned instantaneously as the nematic axis closest to the local flow velocity. A linear stability analysis around a uniformly aligned state yields a closed expression for the critical self-propulsion speed V0* that suppresses the classical active nematic instability. Numerical simulations then report a non-monotonic dependence of nematic order on V0, with a peak near V0≈0.06, defect configurations becoming anti-hyperuniform (α≈−0.9, β≈0.7), long-range vorticity correlations (ν<2), and rotational symmetry breaking that is restored at larger V0. The paper interprets these results as evidence that self-propulsion can bridge active turbulence and ordered regimes.","tokens_in":12203,"tokens_out":3209,"duration_ms":40575,"significance":"If the results are robust, the paper offers a conceptually simple control parameter—self-propulsion speed—for tuning active nematics from turbulent to long-range-ordered states, and it connects to recent experiments on cell monolayers and bacterial colonies. The linear stability calculation is a genuine parameter-free prediction from the stated equations and is tested against simulations, which is a strength. The anti-hyperuniformity and vorticity scaling measurements are also internally cross-checked by independent estimators at V0=0.06. However, the central nonlinear results rest entirely on a sharp, non-dynamical polarity closure, and the reported scaling exponents lack error bars and fitting-range sensitivity analyses. These issues must be addressed before the quantitative claims can be considered established.","major_comments":[{"comment":"The model's only polar ingredient is the advective term V0 p·∇Q, with p assigned as the end of n̂ least deviating from the local flow, i.e. p = sign(n̂·v)n̂. This is a discontinuous, non-dynamical closure: it is undefined when n̂·v=0, and no equation of motion, persistence, or noise is introduced for p. Every V0-dependent headline result—the order/correlation-length peak at V0≈0.06, defect anti-hyperuniformity, and vorticity long-range order—is conditioned on this specific rule. The cited experiments [27–29] demonstrate flow alignment in particular settings, but they do not establish that intrinsic polarity is slaved to the local velocity on all scales, especially near topological defects. I request a sensitivity analysis with a regularized closure (e.g., p = n̂ tanh(β n̂·v)) or with a dynamical p (even a simple relaxational equation) to show that the predicted optimum and exponents do n","section":"§II, Eq. (1) and Fig. 1"},{"comment":"The exponents α in Fig. 4b are obtained by fitting S(q)∼q^α over |q|≤1/4, but no error bars, number of fitting points, or sensitivity to the fitting range are reported. Fig. 4a shows visible curvature and flattening for V0≥0.08, and the distinction between asymptotic anti-hyperuniformity and finite-range apparent scaling is exactly what such an analysis must demonstrate. The independent window-scaling estimate at V0=0.06 (β≈0.7, α≈−0.8) is welcome, but it is only given for two V0 values. I ask for confidence intervals on α and for explicit q_min sensitivity tests for at least V0=0.04, 0.06, 0.08, and 0.1.","section":"§III C, Fig. 4"},{"comment":"The energy-cascade exponent β and vorticity-decay exponent ν are extracted from insets without error bars and without stating the fitting intervals. The claim ν<2 (long-range vorticity order) depends on the small-q behaviour of the enstrophy spectrum / large-r decay of ⟨ω²(r)⟩; the insets appear to cover only a limited range and the reported values are plotted without uncertainties. Please provide the fit ranges, confidence intervals, and a demonstration that ν is stable under varying the fitting window. This is load-bearing for the 'anomalous long-range order' conclusion.","section":"§III D, Fig. 5"}],"minor_comments":[{"comment":"Typo: 'longitude' should be 'length' when defining L.","section":"§III A"},{"comment":"Reference [19] contains a typo: 'Physical Revew E' should be 'Physical Review E'.","section":"References"},{"comment":"The definition ρe = ⟨E/max(E)⟩ is ambiguous: max(E) over what set—time, space, or both? Please state the normalization explicitly.","section":"§III B, Fig. 3d"},{"comment":"The inset defines l∥ and l⊥ but the text does not state how these lengths are extracted. Please define them in the caption or text.","section":"§III D, Fig. 6"},{"comment":"The number of frames and the stationarity of the defect configuration used for the window counts are not stated. A brief note on temporal averaging and convergence would help.","section":"§III C, Fig. 4c"}],"recommendation":"major_revision","confidential_remarks":"The central modelling assumption is the weakest point. It is explicitly stated, which is good, but it is also the sole source of all V0-dependent nonlinear physics. The stress-test note is fair: without a regularized or dynamical polarity, the quantitative predictions are not yet anchored to the claimed physical mechanism. I do not see grounds for rejection, because the assumption can be tested within the manuscript's scope by adding a controlled regularization and repeating the key measurements."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: the linear stability analysis is the best part, and it checks out. I verified the Jacobian in Eq. A3, the eigenvalue expression in Eq. 8, and the threshold in Eq. 9; the algebra is correct. The paper then adds an advective self-propulsion term to the standard active nematic equations and shows in simulation a well-defined intermediate V0 window (peaking at V0≈0.06) where nematic order is maximal, defects become anti-hyperuniform with α≈−0.9, and vorticity correlations decay slower than 1/r^2. Those are new results and they give the paper a real experimental hook.\n\nThe soft spot is the polarity rule. p is not an independent field; it is slaved to the flow via p = sign(n·v)n at every time step. This rule is discontinuous at n·v=0, and there is no equation of motion, persistence, or noise for p. Consequently, every V0-dependent headline—the optimum, the exponents, the symmetry breaking—is a property of this particular closure, not of a genuine self-propulsion degree of freedom. The experiments cited for flow alignment are suggestive but they don't show that intrinsic polarity is locked to the local velocity everywhere, especially near topological defects where n·v changes sign. The paper does no sensitivity test against a regularized rule or a dynamical p, so the quantitative predictions are conditional. That is the main reason for a conditional verdict, not any fatal error.\n\nBeyond that, the exponent fits mostly have no error bars (only one α is quoted with uncertainty), and the text contains a 'V0∈[0.4,0.6]' interval that should clearly be '[0.04,0.06]' from context and from the later sentence. The defect detection protocol is not given, no code or data is shipped, and the 'flocking transition' that brackets the ordered window is never directly characterized. These are addressable issues.\n\nBottom line: this is a paper for active matter people interested in the polar–nematic crossover. It deserves a serious referee: the linear theory is solid, the simulation findings are suggestive, and the writing is clear. But the referee should insist on a robustness analysis of the polarity rule, error bars, and a clean-up of the typo. I'd take the results with a grain of salt until then.","headline":"The linear stability checks out, but the simulation results all hang on an unregularized flow-locked polarity rule; worth refereeing with major revision.","tokens_in":12846,"tokens_out":3301,"would_cite":false,"duration_ms":36888,"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":"Adding a self-advective term to the active nematic equations suppresses the nematic instability above a threshold speed, and at intermediate speeds it maximizes nematic order, turns defect arrangements anti-hyperuniform, and produces long-r","keywords":["active nematics","self-propulsion","active turbulence","nematic instability","topological defects","anti-hyperuniformity","long-range vorticity order","rotational symmetry breaking"],"falsifier":"A concrete test: run a simulation (or experiment) of self-propelled nematic rods in which the polarity is an independent field with its own relaxation, and measure the correlation length, defect structure factor, and vorticity exponent as a function of swimming speed. If the peak at V0≈0.06 disappears, the predicted optimum and the associated long-range order are artifacts of the flow-locking rule; if the peak persists, the slaving assumption is not the origin of the effect.","tokens_in":11713,"feed_emoji":"🌀","tokens_out":6443,"duration_ms":63896,"temperature":0.7,"pith_summary":"This paper argues that a minimal addition to the standard active nematic model—a term that lets each particle advect along a polarity direction—changes the physics of active nematics far beyond a perturbative shift. It shows, first, that self-propulsion above a critical speed V0* stabilizes the homogeneous nematic state, postponing the classical active instability. Second, at intermediate speeds (peaking near V0≈0.06 at activity ζ=0.05) the nematic order and correlation length reach a maximum, and the topological defects arrange into anti-hyperuniform patterns with giant number fluctuations, while vorticity correlations decay with an exponent ν<2, meaning true long-range order. Third, in this same window the system's rotational symmetry is spontaneously broken, and it is restored at larger speeds. If correct, self-propulsion is not just a perturbation but a quantitative control parameter that bridges active turbulence and flocking-like order.","feed_headline":"Self-propulsion gives active nematics long-range order","feed_subtitle":"At one intermediate speed, defects become anti-hyperuniform and vorticity order stretches across the system.","key_machinery":"The central object is the self-advective term V0 p_k ∂_k Q_ij in Equation (1), which breaks nematic symmetry by giving each nematic particle a polarity. The polarity p is not an independent field: it is slaved to the director n̂ and the flow velocity v, with p pointing along the end of the nematic axis that makes the smallest angle with the local flow. This term enters the linear stability analysis through the longitudinal mode, producing the closed-form threshold V0* (Equation 9) that separates the stable and unstable regions. In the nonlinear regime, the same term is responsible for the non-monotonic enhancement of order, the anti-hyperuniform defect configurations, and the broken rotation","core_discovery":"The central claim is that adding the self-advective term V0 p_k ∂_k Q_ij to the Beris–Edwards equation for the nematic tensor—with the polarity p chosen at each instant as the director end closest to the local flow—makes the homogeneous nematic state linearly stable when V0 exceeds the closed-form threshold V0* = (ΓK + η/ρ)√[(2+λ)/(2ΓKη) (ζ − q²K(2+λ)) − q²], verified by simulations. At lower but non-zero speeds, before the flocking transition, the same term enhances nematic order non-monotonically: at V0≈0.06 the elastic free energy density is minimal and the correlation length maximal. In this ordered window, topological defects exhibit anti-hyperuniform density fluctuations (structure fac","pith_inferences":["If the polarity-slaver assumption is relaxed to allow an independent polarity field with its own relaxation and contact alignment, the sharp optimum at V0≈0.06 might broaden or shift; a direct test would be to simulate the nematopolar models cited as [23, 24] with the same parameters.","The anti-hyperuniform defect packing at the order maximum suggests a connection to critical-like density fluctuations: measuring the compressibility or the structure factor of the defect gas at very small q could reveal a universality class shared with other active systems at a nonequilibrium critical point.","The resumption of rotational symmetry at higher V0 suggests that self-propulsion acts like a reentrant symmetry-breaking field; it may be worth testing whether this reentrance persists when activity ζ is varied along the stability boundary.","Because the model uses a single scalar V0, it yields a testable quantitative prediction for real systems: for a given suspension, the correlation length should first grow, peak, and then fall as swimming speed is increased, with the peak location set by the ratio ΓK/η times the V0* expression."],"forward_implications":["Above V0*, active nematic suspensions with self-propulsion can remain in a homogeneous, aligned state at activities that would otherwise drive spontaneous turbulence.","At intermediate speeds, the defect network becomes anti-hyperuniform, meaning density fluctuations grow faster than the Poisson law—observable as giant number fluctuations in experiment.","Vorticity correlations decay with an exponent smaller than the spatial dimension, implying long-range order in the flow and a symmetry-broken, anisotropic state.","The non-universal scaling exponents in the kinetic-energy spectrum provide a fingerprint that could be used to detect self-propulsion in experimental systems such as migrating cell monolayers or swarming bacteria.","Tuning self-propulsion can therefore be used as a design knob to control order, defect structure, and flow correlations in synthetic active-nematic materials."],"supporting_citations":[{"why":"Supplies the continuum active nematic framework (Beris–Edwards + Navier–Stokes) that the model extends.","marker":"[13]"},{"why":"Provides the linear stability analysis of the homogeneous nematic state that the paper generalizes.","marker":"[31]"},{"why":"Grounds the self-advective term V0 p·∂Q in the polar active fluid formalism.","marker":"[25]"},{"why":"Provides the hydrodynamic treatment and the giant number fluctuations context connecting to anti-hyperuniformity.","marker":"[26]"},{"why":"Defines hyperuniform/anti-hyperuniform scaling (S(q)~q^α, σ²_N~R^{d−α}) used to classify defect configurations.","marker":"[35]"},{"why":"Experimental evidence that Dictyostelium cells align polarity with flow, supporting the polarity assignment rule.","marker":"[27]"},{"why":"Experimental evidence of lamellipodia/intrinsic polarity alignment with force/flow in epithelial monolayers.","marker":"[28]"},{"why":"Shows Myxococcus xanthus aligns with flow direction, grounding the polarity rule for bacteria.","marker":"[29]"},{"why":"Baseline for the V0=0 kinetic energy spectrum that the paper compares against.","marker":"[38]"}],"fun_headline_variants":["Self-propulsion turns active nematic turbulence into order","Optimal self-propulsion maximizes order in active nematics","A little self-propulsion orders active nematic flow","Self-propulsion at moderate speed orders nematic vorticity","Self-propulsive nematic defects align at sweet spot speed"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is the polarity assignment rule: each particle is assumed to always polarize along the end of its nematic axis that makes the smallest angle with the local flow velocity, so the self-propulsion direction has no independent dynamics.","fun_headline_variants_meta":{"raw":{"variants":["Self-propulsion turns active nematic turbulence into order","Optimal self-propulsion maximizes order in active nematics","A little self-propulsion orders active nematic flow","Self-propulsion at moderate speed orders nematic vorticity","Self-propulsive nematic defects align at sweet spot speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000268,"raw_usage":{"total_tokens":1446,"prompt_tokens":729,"completion_tokens":717,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":473,"completion_tokens_details":{"reasoning_tokens":638}},"tokens_in":473,"tokens_out":717,"duration_ms":8309,"temperature":1.0,"reasoning_tokens":638,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:40:02.602982+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: run a simulation (or experiment) of self-propelled nematic rods in which the polarity is an independent field with its own relaxation, and measure the correlation length, defect structure factor, and vorticity exponent as a function of swimming speed. If the peak at V0≈0.06 disappears, the predicted optimum and the associated long-range order are artifacts of the flow-locking rule; if the peak persists, the slaving assumption is not the origin of the effect.","supporting_citations":[{"cited_title":"Local polar order controls mechanical stress and triggers layer formation in developing Myxococcus xanthus colonies","cited_arxiv_id":"2308.00368","evidence_quote":"Provides the linear stability analysis of the homogeneous nematic state that the paper generalizes."},{"cited_title":"2012 Polar patterns in active fluids.Soft Matter8, 129–139","cited_arxiv_id":null,"evidence_quote":"Grounds the self-advective term V0 p·∂Q in the polar active fluid formalism."},{"cited_title":"2021 Inertia Drives a Flocking Phase Transition in Viscous Active Fluids.Physical Review X 11, 031063","cited_arxiv_id":null,"evidence_quote":"Defines hyperuniform/anti-hyperuniform scaling (S(q)~q^α, σ²_N~R^{d−α}) used to classify defect configurations."},{"cited_title":"This behaviour is also shared for bacterium, such asMyxococcus xanthus, where each individual particle aligns with the flow direction [29]","cited_arxiv_id":null,"evidence_quote":"Experimental evidence that Dictyostelium cells align polarity with flow, supporting the polarity assignment rule."},{"cited_title":"2003 Shear flow-induced motility of Dictyostelium discoideum cells on solid substrate.Journal of Cell Science 116, 4331–4338","cited_arxiv_id":null,"evidence_quote":"Shows Myxococcus xanthus aligns with flow direction, grounding the polarity rule for bacteria."},{"cited_title":"2021 Local Number Fluctuations in Hyperuniform and Nonhyperuniform Systems: Higher-Order Moments and Distribution Functions","cited_arxiv_id":null,"evidence_quote":"Baseline for the V0=0 kinetic energy spectrum that the paper compares against."}],"review_version":1}