{"id":"4c033b63-934a-4cdf-8c0c-090dbf46cf5c","arxiv_id":"2608.05851","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A nonperturbative renormalization group calculation finds the ordered-state exponents of Malthusian flocks and a nonuniversal quasi-long-range-ordered phase whose endpoint differs from the BKT transition.","lead":"This paper derives the long-distance behavior of Malthusian flocks, collections of self-propelled particles that constantly die and reproduce. It finds a stable ordered state with specific scaling laws, plus a fragile intermediate phase with slow, power-law order that likely disappears when defects are included.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The QLRO claim rests on excluding defects via the spin-wave approximation; since the paper's own equilibrium overlay places its QLRO region past the BKT vortex-unbinding threshold, the physical phase remains unestablished until defects are included.","rationale":"I checked the fixed-point algebra in the smooth-regulator flow equations (8)-(11) with the anomalous dimensions (65)-(66). The Malthusian fixed point (bar_lambda^2 = 8/5, bar_mu_y = -bar_mu_3 = 1/5, bar_mu_x = 0) solves the equations with eta_x = 1/4, eta_y = 3/4, and the Gaussian line bar_lambda = bar_mu_3 = 0, bar_mu_x = bar_mu_y = c is indeed a line of fixed points whose lambda-eigenvalue is 1 - c. So the internal NPRG calculation is coherent. The weakest point is external: the QLRO line and endpoint are properties of the spin-wave (defect-free) truncation, and the authors' own comparison with equilibrium BKT indicates the predicted QLRO region is precisely where vortex unbinding would occur in the equilibrium XY model. Without a defect theory or a simulation including amplitude fluctuations, the central 'Quasi Long-Range Order in Malthusian Flocks' claim is not established for the physical system. Secondary issues, such as the apparent typo in sharp-regulator Eq. (79) (first term multiplies bar_mu_y rather than bar_mu_3) and the absence of the referenced Mathematica notebooks, support caution but are not the main load-bearing problem. The reader's weakest_assumption is the same one, so I agree; the CONDITIONAL verdict stands.","tokens_in":21775,"tokens_out":10597,"duration_ms":99589,"concrete_test":"Run a 2D lattice simulation of the full polar vector model with amplitude fluctuations and Malthusian turnover (e.g., discretized Supplemental Eq. (1), or an active XY model with birth/death and no fixed |p|) in the predicted QLRO basin, large D/mu and small lambda. Measure G(r) = <p(0) . p(r)>/p0^2 and vortex density rho_v(L). If G(r) decays exponentially or rho_v grows with L, defect unbinding destroys the claimed QLRO; if G(r) ~ r^{-2c} with c > 1 and rho_v -> 0, the spin-wave line survives. A cheaper analytical cross-check is to repeat the NPRG with a two-field ansatz including delta p and the potential U(|p|); if the attractive line and its endpoint disappear, they are artifacts of the spin-wave truncation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is the spin-wave approximation made in deriving the theta-only equation of motion, Eq. (1) and Supplemental Eq. (4): the polarity amplitude is fixed, delta p = 0, so theta is a smooth field and topological defects cannot occur. The paper explicitly states that the QLRO analysis 'crucially hinges' on this assumption and that nonperturbative amplitude fluctuations allow defect-pair creation. In the equilibrium sector the authors map their attractive line to the XY stiffness, bar_mu_x = 1/(4 pi K), overlay the BKT vortex-unbinding point at bar_mu_x = 1/8, and conclude that there is no overlap between the two lines, 'suggesting that QLRO is generically destabilized in Malthusian flocks, either by activity or defect fluctuations.' Their predicted QLRO regime is bar_mu_x > 1, which lies far beyond the equilibrium threshold; no mechanism for nonequilibrium suppression of vortex unbinding is calculated. Therefore the paper establishes QLRO only for a defect-free spin-wave model, not for the advertised Malthusian flock. The strong-coupling fixed point and exponents are internally consistent and match Chate-Solon, and I do not challenge them; but the title-level claim of quasi long-range order in Malthusian flocks, and the new critical point at bar_mu_x = 1, remain conditional on defect physics that is not treated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a nonperturbative renormalization group (NPRG) analysis of two-dimensional Malthusian flocks in the spin-wave approximation, in which the polarity amplitude is held fixed and topological defects are excluded. Using a second-order derivative expansion that respects rotational and chiral symmetry, the authors derive flow equations and identify a strong-coupling fixed point with exponents eta_x = 1/4, eta_y = 3/4, yielding chi = -1/4, zeta = 3/4, and z = 5/4, in agreement with the values predicted by Chate and Solon. They also report an attractive line of fixed points that would realize quasi long-range order for large noise or weak diffusion, ending at a critical point with bar_mu_x = 1 that they argue is distinct from the Berezinskii-Kosterlitz-Thouless universality class. The paper explicitly acknowledges that the quasi-long-range-order analysis relies crucially on the spin-wave approximation and that defect fluctuations are expected to destabilize this phase in the full Malthusian flock model.","tokens_in":22006,"tokens_out":18065,"duration_ms":153076,"significance":"If the strong-coupling fixed point is robust, the paper makes an important contribution by providing a nonperturbative derivation of the Malthusian flock exponents and by demonstrating that rotational invariance alone can reproduce the Chate-Solon scaling, contrary to earlier suggestions that additional symmetries are needed. The work is accompanied by detailed supplemental calculations and Mathematica notebooks, and the authors check two different regulators, which strengthens confidence in the internal consistency of the flow calculation. The quasi-long-range-order line is an interesting prediction for the defect-free spin-wave model and is independently reported by Sezik and Pruessner (Ref. [43]). However, the physical relevance of the quasi-long-range-order phase to actual Malthusian flocks remains conditional on defect physics that is not treated in the manuscript; the authors are candid about this limitation in the body of the paper.","major_comments":[{"comment":"Equation (14) is incorrect as written. For a Gaussian field theta, exp(-1/2 <[theta(0,r)-theta(0,0)]^2>) equals the order-parameter correlation <n(0,r) . n(0,0)>, not the mean-square difference <[n(0,r)-n(0,0)]^2>. The latter is 2 - 2<n(0,r) . n(0,0)> and tends to 2 at large r. Thus the power law (Lambda r)^(-2 bar_mu_x) should be assigned to the correlation function, not to the squared difference. The subsequent discussion of the critical end point at bar_mu_x = 1 (corresponding to an order-parameter correlation exponent of 2) should be stated consistently with this corrected identification.","section":"Results, Eq. (14)"},{"comment":"The title and abstract present quasi long-range order as a phase of Malthusian flocks, but the calculation is performed for the spin-wave model with delta p = 0 (Supplement Eq. (4)), and the text states that the QLRO analysis 'crucially hinges' on this approximation. Moreover, in the same section the authors use the equilibrium BKT mapping to conclude that 'QLRO is generically destabilized in Malthusian flocks, either by activity or defect fluctuations.' The manuscript should therefore explicitly frame the QLRO phase as a property of the defect-free spin-wave model, and adjust the title or abstract so that the conditional nature of this claim is not lost.","section":"Quasi long-range order and novel critical point / Discussion & Outlook"},{"comment":"The claim that the sharp cutoff leaves 'the quantitative values of the scaling exponents' unchanged is not supported by the numbers reported. The sharp-regulator fixed point is given as lambda* approx 0.95, mu_x* = 0, mu_y* approx 0.03, mu_3* approx 0.05. Substituting these values into Eqs. (80)-(81) gives eta_x approx 0.233 and eta_y approx 0.737, which differ by about 7% from the smooth-regulator values 1/4 and 3/4. The authors should report the sharp-regulator exponents explicitly and discuss the size of the regulator dependence, or soften the claim of quantitative agreement.","section":"Supplement: Flow equations with sharp regulator"}],"minor_comments":[{"comment":"The notation 'mu3x mu y' should be written as mu_x^3 mu_y to avoid ambiguity between the cube of mu_x and a coupling mu_3.","section":"Eq. (12)"},{"comment":"In the sharp-regulator flow equation for bar_mu_3, the first term on the right-hand side should read (1/2)(-3 eta_x - eta_y) bar_mu_3 rather than bar_mu_y.","section":"Supplement Eq. (79)"},{"comment":"There are several typographical errors: 'ressembling' should be 'resembling', 'destablized' should be 'destabilized', 'obtaines' should be 'obtains', and 'relavant' should be 'relevant'.","section":"Throughout"},{"comment":"The flow equations in the main text depend on eta_x and eta_y, but the explicit expressions for these anomalous dimensions appear only in the Supplement; including them in an appendix or in the main text would make the fixed-point solution more self-contained.","section":"Main text, Eqs. (8)-(11)"},{"comment":"The sentence 'the nonlinear response of such a state is then governed by the scaling exponents obtained in this work' is stronger than what has been computed; response functions have not been calculated here, so this statement should be softened.","section":"Discussion & Outlook"},{"comment":"Ref. [43] is listed as 'to be published'; if an arXiv number or a published reference is available, it should be provided.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about the spin-wave limitation, which is to its credit. The main issues are the incorrect identification of the correlation function in Eq. (14), the mismatch between the advertised QLRO phase and the conditional nature of the result, and an overstatement of regulator independence. These are fixable in revision and the central strong-coupling fixed point appears sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the ordered-phase calculation is real and worth engaging with; the QLRO claim is honestly labeled as conditional, and the paper's own equilibrium overlay undercuts it. If you read one thing, read the Discussion section.\n\nWhat's new: this is the first full second-order derivative expansion of the NPRG applied to the Malthusian flock equation, and it explicitly finds the strong-coupling fixed point. The exponents eta_x=1/4, eta_y=3/4 reproduce Chate-Solon without invoking extra symmetries, and the two-regulator comparison plus the DRG consistency check give me reasonable confidence the ordered-phase result is not an artifact of the truncation. I checked the fixed point algebra: bar_lambda*=sqrt(8/5), bar_mu_y*=-bar_mu_3*=1/5, bar_mu_x*=0 solves the flow equations with the stated anomalous dimensions. That part holds up.\n\nThe soft spots are in the QLRO claim. The derivation sets delta p=0, which excludes topological defects. The paper says the QLRO analysis 'crucially hinges' on this, and the equilibrium overlay places the attractive line at bar_mu_x>1, far past the BKT unbinding threshold at 1/8. The authors conclude QLRO is 'generically destabilized' in Malthusian flocks, which is honest but means the title claim is not established for the physical system. The new critical point at bar_mu_x=1 is real only in the spin-wave model until someone treats defects. I don't see a way around this within the paper; it's a genuine limitation, not a nitpick.\n\nMinor issues: no code or data is shipped despite references to Mathematica notebooks, which makes verification harder. There's a likely typo in Eq. (79) of the supplemental—the first term multiplies bar_mu_y instead of bar_mu_3, unlike the smooth-regulator Eq. (70). Also, the sharp-regulator fixed point has mu_3* with opposite sign to the smooth one; the authors note this but don't discuss what it means for regulator dependence.\n\nNovelty is moderate. The exponents were already in Ref. [28], and the QLRO line is independently in the unpublished Ref. [43]. The contribution here is the NPRG framework itself, which is genuinely new for this model and may be useful for future work.\n\nOverall: a careful calculation with an honest limitation. The ordered-phase result deserves referee time; the QLRO claim needs to be either defended against defects or reframed. I'd send it to review, with a request to address the defect issue and fix the supplemental typo.","headline":"Solid NPRG calculation that reproduces Chate-Solon exponents, but the QLRO phase is derived for a defect-free spin-wave model, not the physical Malthusian flock.","tokens_in":22613,"tokens_out":2643,"would_cite":true,"duration_ms":23757,"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 rotationally invariant, defect-free renormalization group calculation fixes the scaling exponents of two-dimensional Malthusian flocks at 1/4 and 3/4 and finds a quasi-long-range-ordered phase with a critical endpoint distinct from BKT.","keywords":["Malthusian flocks","active matter","nonperturbative renormalization group","quasi-long-range order","Berezinskii-Kosterlitz-Thouless transition","flocking","spin-wave approximation","universality class"],"falsifier":"A direct numerical test: simulate the spin-wave model (1) on a lattice and measure the equal-time orientation correlation function; the paper predicts algebraic decay with exponent $2\\bar\\mu_x^*$ that reaches 2 at the critical point, whereas exponential decay or a universal BKT exponent of 1/4 at unbinding would refute the claim.","tokens_in":21457,"feed_emoji":"🐦","tokens_out":9406,"duration_ms":86187,"temperature":0.7,"pith_summary":"Flocks whose constituents are continually created and destroyed—cells in a tissue, filaments in a cytoskeleton—are described by Malthusian flock hydrodynamics, but their two-dimensional phase diagram has resisted analytic treatment because nonlinearities cannot be handled perturbatively. This paper claims that a nonperturbative renormalization group calculation, rotationally invariant to second order in derivatives and excluding topological defects, resolves the problem. It explicitly finds the strong-coupling fixed point of the ordered phase with exponents $\\eta_x=1/4$ and $\\eta_y=3/4$ (hence $\\chi=-1/4$, $\\zeta=3/4$, $z=5/4$), recovering the prediction of Ref. [28] without invoking extra symmetries. It also uncovers an attractive line of fixed points that realizes quasi-long-range order in the spin-wave model, ending in a critical point similar to, but distinct from, the Berezinskii-Kosterlitz-Thouless transition. If correct, the universal scaling of ordered flocks with turnover is known, and the route into collective motion is controlled by self-advection alone.","feed_headline":"Birth-death flocks: scaling exponents fixed at 1/4 and 3/4","feed_subtitle":"A defect-free RG predicts universal scaling for 2D flocks with turnover and a quasi-ordered phase beyond BKT.","key_machinery":"The central object is the stochastic equation of motion for the orientation angle $\\theta$ obtained under the spin-wave approximation, Eq. (1): a single field theory with five parameters—self-advection $\\lambda$, isotropic and anisotropic diffusion $\\mu_1,\\mu_2,\\mu_3$, and noise amplitude $D$—whose form is completely fixed by rotational and chiral symmetry at second order in derivatives. The paper converts this Langevin equation into a response-field action and applies the exact nonperturbative renormalization group flow, truncating the effective average action to the same second-order derivative expansion. A mass-like regulator makes the loop integrals analytic, and a sharp cutoff confirms that the fixed-point structure and exponents do not depend on the regulator. The dimensionless couplings $\\bar\\lambda,\\bar\\mu_x,\\bar\\mu_y,\\bar\\mu_3$ then flow according to four coupled equations whose fixed points encode the universal physics.","core_discovery":"On its own terms, the paper establishes that the second-order, rotationally invariant derivative expansion of the effective action for the Goldstone mode of a two-dimensional Malthusian flock has three fixed-point structures. The stable strong-coupling Malthusian fixed point sits at $\\bar\\lambda^*=\\sqrt{8/5}$, $\\bar\\mu_y^*=-\\bar\\mu_3^*=1/5$, $\\bar\\mu_x^*=0$ and carries anomalous dimensions $\\eta_x=1/4$, $\\eta_y=3/4$, which translate into $\\chi=-1/4$, $\\zeta=3/4$, $z=5/4$; because the graphical correction to the self-advection $\\lambda$ vanishes at this fixed point, the exponents follow from rotational symmetry alone rather than an additional symmetry assumption. The paper also finds Toner's fixed point and, unexpectedly, an attractive line $\\bar\\lambda^*=0$, $\\bar\\mu_x^*=\\bar\\mu_y^*=1$, $\\bar\\mu_3^*=0$ that realizes quasi-long-range order, with a critical endpoint at $\\bar\\mu_x^*=1$ whose universal exponent differs from the BKT value $1/8$. The authors state plainly that the quasi-long-range-order part of the analysis depends on the spin-wave approximation and that defect fluctuations are not treated.","pith_inferences":["Editorial inference: the continuous line of fixed points implies a continuously varying quasi-long-range-order exponent; measuring the orientation correlation exponent as a function of noise or diffusion in a lattice simulation would provide a parameter-free check.","Editorial inference: because the critical point is controlled by $\\lambda$, biological collectives could switch from a weakly ordered responsive state to true collective migration by increasing cell motility rather than alignment strength, a functional consequence the authors mention but do not develop.","Editorial inference: if nonequilibrium effects modify defect binding, as reported for some active liquid crystals, the quasi-long-range-ordered phase could survive in real Malthusian flocks; the paper leaves this genuinely open, so a simulation with amplitude fluctuations would settle it.","Editorial inference: in the equilibrium limit the attractive line coincides with the XY spin-wave line and is unstable to defects; the paper's own comparison therefore suggests that purely equilibrium systems would melt the QLRO, while actively driven systems may stabilize it or replace it with a different non-BKT melting line."],"forward_implications":["If the central claim is correct, the long-distance scaling of ordered Malthusian flocks in two dimensions is universal, with exponents $\\eta_x=1/4$, $\\eta_y=3/4$ and derived exponents $\\zeta=3/4$, $z=5/4$.","In the spin-wave model without defects, quasi-long-range order exists and orientation correlations decay as a power law whose exponent is set by where the flow lands on the attractive line of fixed points.","The transition from quasi-long-range order to true long-range order is driven by the self-advection $\\lambda$ alone, so increasing motility can turn a globally disordered but responsive state into a collectively migrating one.","The critical endpoint is a new universality class: it resembles BKT in having a line of fixed points and a universal endpoint exponent, but the endpoint sits at $\\bar\\mu_x^*=1$, not the BKT value $1/8$.","Because the derivation uses only symmetries and conservation laws, the results should apply to any polar active system maintained at density homeostasis, including tissues and cytoskeletal networks."],"supporting_citations":[{"why":"Supplies the hydrodynamic Malthusian flock model that the paper analyzes.","marker":"[23]"},{"why":"Predicts the scaling exponents $\\eta_x=1/4$, $\\eta_y=3/4$ that the calculation reproduces without assuming extra symmetries.","marker":"[28]"},{"why":"Argues the nonlinearity structure is insufficiently constrained by symmetry; the paper shows rotational symmetry alone suffices at second order.","marker":"[29]"},{"why":"Independently obtains the quasi-long-range-ordered line and critical point with a field-theoretic RG, corroborating the QLRO finding.","marker":"[43]"},{"why":"Defines the Berezinskii-Kosterlitz-Thouless universality class used as the comparison for the new critical endpoint.","marker":"[44–46]"},{"why":"Raises the possibility that defect fluctuations destabilize even the true long-range-ordered phase, bounding the physical conclusions.","marker":"[26]"},{"why":"Provides the nonperturbative renormalization group formalism on which the calculation is based.","marker":"[34]"},{"why":"Gives the extension of the nonperturbative RG to nonequilibrium steady states, needed to treat the Langevin equation.","marker":"[48]"}],"fun_headline_variants":["Malthusian flocks: exponents pinned to 1/4 and 3/4","Quasi long-range order in birth-death flocks","Beyond BKT: new phase in Malthusian flocks","Defect-free RG: true and quasi order in flocks","Malthusian flock theory: new scaling universality"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the spin-wave approximation: the polarity amplitude is frozen, so the angle field is smooth and topological defects cannot form; if defect pairs unbind in the full model, the quasi-long-range-ordered phase may not exist and the ordered-phase exponents could also shift.","fun_headline_variants_meta":{"raw":{"variants":["Malthusian flocks: exponents pinned to 1/4 and 3/4","Quasi long-range order in birth-death flocks","Beyond BKT: new phase in Malthusian flocks","Defect-free RG: true and quasi order in flocks","Malthusian flock theory: new scaling universality"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000634,"raw_usage":{"total_tokens":2905,"prompt_tokens":908,"completion_tokens":1997,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":1909}},"tokens_in":524,"tokens_out":1997,"duration_ms":14253,"temperature":1.0,"reasoning_tokens":1909,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T22:22:38.997617+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical test: simulate the spin-wave model (1) on a lattice and measure the equal-time orientation correlation function; the paper predicts algebraic decay with exponent $2\\bar\\mu_x^*$ that reaches 2 at the critical point, whereas exponential decay or a universal BKT exponent of 1/4 at unbinding would refute the claim.","supporting_citations":[{"cited_title":"Chat´ e and A","cited_arxiv_id":null,"evidence_quote":"Predicts the scaling exponents $\\eta_x=1/4$, $\\eta_y=3/4$ that the calculation reproduces without assuming extra symmetries."},{"cited_title":"Sezik and G","cited_arxiv_id":null,"evidence_quote":"Independently obtains the quasi-long-range-ordered line and critical point with a field-theoretic RG, corroborating the QLRO finding."},{"cited_title":"Besse, H","cited_arxiv_id":null,"evidence_quote":"Raises the possibility that defect fluctuations destabilize even the true long-range-ordered phase, bounding the physical conclusions."},{"cited_title":"Canet, H","cited_arxiv_id":null,"evidence_quote":"Gives the extension of the nonperturbative RG to nonequilibrium steady states, needed to treat the Langevin equation."}],"review_version":1}