{"id":"943d48c4-2d0b-42c6-b351-93244b67e7fc","arxiv_id":"2506.03795","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Turnover generates outward flows from defect cores that overcome elastic attraction, stabilizing topological defects in compressible active polar fluids and producing lattices, foams, and vortex glasses.","lead":"This paper shows that turnover, the continuous assembly and disassembly of an active fluid, can stabilize topological defects in compressible polar active matter, producing states like foams, lattices, and vortex glasses. The result matters because stable defect patterns are thought to organize stress in biological materials such as the actin cortex and tissues.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The defect-stabilization mechanism vanishes when the density-polarity coupling chi=0: Eqs. (13)-(14) give zero turnover-driven radial flow, so the central claim is conditional on the positive chi assumed in Eq. (1).","rationale":"The reader's weakest-assumption analysis identified the positive density-polarity coupling as the key input, and I agree that this is the most load-bearing point. It is the only place in the argument where the central mechanism can fail within the model: if chi=0, the analytic core model gives zero radial flow, and Eq. (18) reduces to the usual defect-annihilation dynamics. The numerics support the conditional claim for chi=0.1, and I find no internal inconsistency in the derivation; the issue is external validity of a key modeling assumption. The other limitations noted by the reader (code not yet released, single-run phase diagram, absence of defect creation in the deterministic dynamics) are important for reproducibility and for the strength of the biological conclusions, but they do not by themselves falsify the stabilization mechanism. The proposed check is inexpensive and would settle whether the effect is generic turnover physics or a property of the chosen density-polarity coupling. For these reasons the reader's CONDITIONAL verdict remains appropriate; I would not move it to ACCEPT or REJECT without this check.","tokens_in":22246,"tokens_out":11321,"duration_ms":112442,"concrete_test":"Run the same ensemble used in Fig. 2(a) (rho0=0.6, tau=1, zeta_rho=4, L=10, 50 runs) for chi=0.01 and chi=0.001, keeping all other Table I parameters fixed, and plot the steady-state number of defects versus chi. If the defect number extrapolates to zero as chi -> 0, this confirms that epsilon proportional to chi (Eq. 14) is the sole source of the repulsive radial flow and that turnover only amplifies a density-polarity effect. A complementary analytical check is to solve the mass-balance equation (12) with chi=0; Eq. (14) already predicts epsilon=0, but the numerical test rules out any alternative turnover-induced repulsion mediated by zeta_rho. If the defect density remains finite for small chi, the mechanism is more robust than the chi-proportional formula suggests.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that turnover stabilizes topological defects. In the core model, the radial outflow that repels defects is v = tau^{-1} R epsilon / 2 (Eq. 13), with epsilon given by Eq. (14). Because epsilon is proportional to chi and vanishes at chi=0, without the density-polarity coupling chi > 0 in the free energy (1) there is no depletion of the defect core, no radial flow, and no repulsive contribution in the defect-pair balance (18) beyond the conventional elastic attraction. Thus turnover alone does not stabilize defects; the stabilization is produced by turnover acting on a density field depleted at defect cores because |p|^2 = rho/rho0 is imposed by Eq. (1). This coupling is an input assumption about actin-cortex physics, not derived or experimentally validated in the paper. The simulations use chi=0.1 (Table I), and the analytic treatment sets kappa=0 inside the core, leaving the core radius R as a free parameter, so the mechanism is not quantitatively pinned down independently of the numerics. The abstract statement that \"turnover readily leads to a stabilization of defects\" is therefore conditional: it holds for positive chi, and the sign and magnitude of this coupling are not established. If chi were absent or negative, the same construction would predict annihilation rather than stabilization.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a hydrodynamic theory of compressible active polar fluids that includes turnover (assembly/disassembly of the active component), density-dependent active stress, and a density-polarity coupling in the free energy. The central claim is that turnover stabilizes topological defects in the polar order field: density depletion at defect cores, combined with turnover-driven radial outflow, repels oppositely charged defects and prevents annihilation. The authors support this claim with long-time numerical solutions in two dimensions, identifying several asymptotic states (active foams, density waves, vortex glasses, and defect lattices) as functions of target density, turnover rate, and activity. They also present a simplified two-domain analytic model of an isolated defect and a defect-pair force balance, along with linear stability analyses of homogeneous polarized and isotropic states. The paper argues that turnover, rather than chaotic active turbulence or suppressed hydrodynamic interactions, is a generic mechanism for organizing defects in biological active matter.","tokens_in":22548,"tokens_out":4171,"duration_ms":45379,"significance":"If the central claim holds, the paper identifies a plausible and previously underappreciated mechanism by which topological defects can be stabilized in biological active matter, with potential relevance to actin-cortex organization and tissue morphogenesis. The strengths of the paper are its extensive numerical parameter sweeps, the explicit derivation of the hydrodynamic model, the transparent stability analysis of homogeneous states, and falsifiable predictions such as the critical activity threshold in Eq. (21) and the dependence of defect density on rho_0, zeta_rho, and tau. However, the significance is moderated by the fact that the stabilization mechanism is conditional on a positive density-polarity coupling chi and by the acknowledged free parameters (defect core radius R and interface thickness l) in the analytic defect-core model. The paper is a solid contribution to active-matter theory if these caveats are made explicit and the analytic model is presented as a scaling argument rather than a closed quantitative theory.","major_comments":[{"comment":"The radial outflow v(R-) = tau^{-1} R epsilon / 2 and the density contrast epsilon both vanish identically when the density-polarity coupling chi = 0, because epsilon in Eq. (14) is proportional to chi. Therefore the stabilization mechanism is not produced by turnover alone but by turnover acting on a density profile that is depleted at defect cores only because the equilibrium relation |p|^2 = rho/rho_0 is imposed by Eq. (1). The abstract's statement that 'turnover readily leads to a stabilization of defects' and the claim in Sec. III that turnover stabilizes defect pairs 'both in the presence or absence of active stress' should be qualified to state explicitly that positive chi is required. The paper would also benefit from a brief discussion of the physical evidence or modeling precedent for the sign and magnitude of chi, since the simulations use only chi = 0.1 (Table I) and the entire mechanism disappears for chi = 0.","section":"Abstract and Sec. III.A, Eqs. (13)-(14)"},{"comment":"The analytic defect-core calculation neglects the Frank free energy (kappa = 0), leaves the core radius R and the interface thickness l as free parameters, and the authors state that chemical and mechanical balance cannot be enforced simultaneously. As a consequence, Eq. (18) and the stabilization criterion tau^{-1} R^2 epsilon > kappa / bar{gamma} are scaling relations rather than closed quantitative predictions. This becomes load-bearing when Fig. 3 compares the numerical phase boundary to Eq. (21) using the relation a = 4 zeta_c_rho / 3 from Table I; this relation is an ad hoc modeling choice, not a derived or measured parameter. Please state clearly that R, l, and the a-zeta_c_rho relation are adjustable within the analytic model, and provide a sensitivity check showing how the predicted boundary in Fig. 3 changes under reasonable variations of R and l.","section":"Sec. III.A and Fig. 3"},{"comment":"The defect-pair balance uses the single-defect outflow v(r=d) as if it were the core-boundary value v(R-) = tau^{-1} R epsilon / 2. For the argument to work, the radial outflow must decay over a length comparable to R so that it is significant at d ~ R and negligible at d >> R, but the paper does not provide the radial decay profile or a derivation of v(r) outside the core. Without this profile, the statement that elastic attraction dominates at large d because radial flows are 'localized near the defect center' is only qualitative. A measurement of v(r) from the simulations, or a matched asymptotic solution, would substantially strengthen the analytic mechanism.","section":"Sec. III.B, Eq. (18)"}],"minor_comments":[{"comment":"The text says 'The code can be found at [51]', but Ref. [51] states 'Code will be made available upon publication.' Please provide a working repository link or change the wording to reflect the intended availability.","section":"Sec. IV.B and Ref. [51]"},{"comment":"There is a typo, 'perfromed', in the first paragraph of Sec. V; it should be 'performed'.","section":"Sec. V"},{"comment":"The discussion of the critical activity zeta_c_rho and the comparison in Fig. 3 would be easier to follow if the text explicitly stated that Eq. (21) is solved self-consistently when a = 4 zeta_c_rho / 3, and how this affects the shape of the predicted boundary.","section":"Sec. IV.A and Eq. (21)"},{"comment":"The caption states that each data point comes from a single numerical solution; given the known variability of defect counts, a brief note on the expected statistical error or a representative error bar would help the reader judge the phase boundaries.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed active-matter theory paper with a novel and plausible mechanism. The main risk is not the mathematics but the breadth of the central claim: the stabilization effect is conditional on a positive density-polarity coupling and on several free parameters in the analytic model. I recommend major revision rather than rejection because the numerical evidence is substantial and the required changes are clarifications and additional supporting analysis, not a fundamental reworking of the model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the useful news: this paper identifies a new mechanism for keeping topological defects alive in active polar fluids. Instead of relying on anisotropic active stress or friction, turnover (assembly/disassembly) combined with a density-polarity coupling depletes the active fluid at defect cores, generating a radial outflow that repels opposite-charged defects. The mechanism is supported by extensive numerical simulations across parameter space, a clean two-defect force balance (Eq. 18), and a linear stability analysis that rationalizes the observed phases. The survey of states—active foams, density waves, vortex glasses, defect lattices—is a solid piece of numerical phenomenology. I also give the authors credit for stating the main limitation themselves: without fluctuations the dynamics does not create defects, so the link to biological organization is indirect.\n\nWhere are the soft spots? The stabilizing effect is conditional on the sign of the density-polarity coupling chi in Eq. (1). The stress-test note is right that at chi=0 the radial flow vanishes, but the paper never claims otherwise; the mechanism is explicitly built on that coupling. The real question is whether chi>0 describes actin cortex physics, and that is not established experimentally. The analytic core model leaves the core radius R and interface thickness l as free parameters and drops Frank elasticity inside the core, so Eqs. (13)-(14) are a rationalization rather than a quantitative derivation. Fig. 3's phase diagram is based on single runs per point without error bars, and the code is promised but not yet available. These are addressable issues, not fatal ones.\n\nOverall, the central claim holds up as a statement about the model, and the paper is honest about what it does not do. It will be of most value to active matter theorists and people modeling the actin cortex or tissues. I would send it to peer review; with code release, ensemble statistics for the phase diagram, and a more careful phrasing of the abstract, it could be a solid publication.","headline":"A genuinely new turnover-based mechanism for stabilizing defects in active polar fluids, well-supported by numerics but conditional on the density-polarity coupling chi; deserves referee time.","tokens_in":23036,"tokens_out":2559,"would_cite":true,"duration_ms":24263,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Turnover stabilizes topological defects in compressible active polar fluids.","keywords":["active polar fluid","turnover","topological defects","defect stabilization","active foam","defect lattice","vortex glass","density-polarity coupling"],"falsifier":"A numerical experiment with the same equations but $\\chi=0$ should show no defect stabilization: opposite-charge pairs should annihilate for all turnover rates. Equivalently, in a simulation with $\\chi>0$, measure the density profile and velocity field around an isolated $+1$ defect: if the claim is correct, the density should dip at the core and the radial velocity should grow with turnover rate roughly as $\\tau^{-1}R\\varepsilon$, while a defect pair should settle at a finite separation that shrinks as $\\chi\\to0$.","tokens_in":22048,"feed_emoji":"🌀","tokens_out":5831,"duration_ms":56445,"temperature":0.7,"pith_summary":"The paper argues that continuous assembly and disassembly of a compressible active polar fluid—turnover—keeps topological defects from annihilating, so the fluid can sustain persistent populations of defects instead of relaxing to a uniform state. The mechanism is a density-polarity coupling: where polarity vanishes at a defect core, the equilibrium density is lower, and turnover then drives a radial outflow that pushes neighboring defects apart. The authors derive a velocity balance for a defect pair and confirm by simulation that the outflow defeats the elastic attraction at intermediate separations. If right, this gives a generic, anisotropic-stress-free route for biological matter such as the actin cortex to organize active stress around stable defect arrays, foams, and vortex glasses.","feed_headline":"Turnover stabilizes defects in compressible active polar fluids","feed_subtitle":"A density-polarity coupling makes defect cores push fluid outward, so opposite charges repel instead of annihilating.","key_machinery":"The load-bearing object is the density-polarity coupling $\\chi>0$ in the free energy $f = \\frac{a}{4}\\rho^4 + \\rho^2\\left[-\\frac{\\chi}{2}\\frac{\\rho}{\\rho_0}|p|^2 + \\frac{\\chi}{4}|p|^4 + \\frac{\\kappa}{2}(\\nabla p)^2\\right]$, which sets $|p|^2 = \\rho/\\rho_0$ in equilibrium so that defect cores ($p=0$) are density-depleted. Around such a core, mass conservation with turnover (source term $-\\tau^{-1}(\\rho-\\rho_0)$) produces a radial outflow; the two-domain core model gives the radial velocity $v(R^-)=\\tau^{-1}R\\varepsilon/2$ with $\\varepsilon = \\frac{3\\chi\\rho_0^2}{3\\chi\\rho_0^2+12\\rho_0^3(a\\rho_0-\\zeta_\\rho)+2\\xi\\tau^{-1}R\\ell}$. This outflow enters the defect-pair force balance $\\dot{d}=2v(r=d)-\\kappa/(d\\bar{\\gamma})$, the equation that determines when defects repel rather than annihilate.","core_discovery":"The central claim is that turnover stabilizes topological defect pairs in a compressible polar active fluid, with or without anisotropic active stress. In equilibrium the fluid's polarity magnitude obeys $|p|^2 = \\rho/\\rho_0$ for $\\chi>0$, so a defect core, where $p=0$, is depleted of active fluid. Turnover, modeled by the source term $-\\tau^{-1}(\\rho-\\rho_0)$, then sustains a radial Darcy-like outflow from the core with velocity $v(R^-)=\\tau^{-1}R\\varepsilon/2$ (Eq. 13), where $\\varepsilon$ is the density contrast set by stress balance and mass conservation (Eq. 14). For two opposite-charge defects the separation obeys $\\dot{d}=2v(r=d)-\\kappa/(d\\bar{\\gamma})$ (Eq. 18), so the outflow advects each defect away from the other while elastic interactions ($\\sim\\kappa/d$) pull them together; at intermediate distances the outflow wins, giving a stable nonzero separation. Long-time numerical solutions show that this stabilization organizes defects into active foams, density waves, vortex glasses, and spontaneously forming square or hexagonal defect lattices, with the phase selected by the turnover rate and target density, rationalized by a linear stability analysis of homogeneous states that turns the instability from type II to type I.","pith_inferences":["The short-range repulsion between defects could be viewed as an effective defect gas with a repulsive core set by the outflow zone; this suggests that collective defect statistics might follow from an effective interacting-particle model with a single length scale $R$.","Since defects are not spontaneously created in the noiseless theory, the model predicts that turnover regulates the fate of pre-existing defect populations; adding noise or defect-pair nucleation could turn it into a full theory of defect number selection.","The same density-polarity coupling could stabilize defects in active nematics with turnover if the nematic order parameter is density-dependent, so the mechanism may generalize beyond polar systems.","A concrete testable extension: in systems where turnover can be inhibited (e.g., by drug treatment), defect density should drop and pairs should annihilate, while increasing turnover should restore finite defect separations."],"forward_implications":["If turnover stabilizes defects without requiring anisotropic active stress, experimental systems with tunable turnover (reconstituted actomyosin, cell monolayers) should exhibit sustained defect populations controlled by assembly and disassembly rates.","Defect lattices in this model move at constant velocity without defect rearrangements, implying a genuinely self-propelled crystalline state of singularities.","The linear stability result predicts that turnover changes the onset of pattern formation from a type II to a type I instability, so the characteristic wavelength of the emerging pattern should be set by $\\tau$ and by the density diffusion coefficient.","Because the mechanism is independent of active stress anisotropy, it applies to both contractile and extensile situations as long as the density-polarity coupling is positive; anisotropic stress and flow alignment then only select defect subtypes or modify the pattern.","The phase sequence with increasing target density—foams, waves, vortex glass, lattices, uniform—can serve as a phase diagram for experiments mapping turnover rates."],"supporting_citations":[{"why":"Provides the two-component active-gel description of the actin cortex and the density-polarity coupling used as the starting hydrodynamic model.","marker":"[37]"},{"why":"Supplies the standard elastic interaction force $\\sim\\kappa/d$ between defects and the defect core energetics.","marker":"[16]"},{"why":"Establishes the defect-pair dynamics in active nematics whose anisotropic-stress-driven separation this paper contrasts with turnover-driven separation.","marker":"[47]"},{"why":"Gives the overdamped defect-position dynamics with effective friction coefficient that the pair velocity balance Eq. (18) uses.","marker":"[48]"},{"why":"Provides the generic constitutive equations of polar active gels underlying Eqs. (6)-(9).","marker":"[44]"},{"why":"Supplies the Cross-Hohenberg classification used to distinguish the type I instability found with turnover from the type II instability without it.","marker":"[49]"},{"why":"Prior work stabilizing defect lattices via substrate friction, the alternative mechanism this paper's turnover mechanism is compared against for lattice formation.","marker":"[30]"},{"why":"Shows friction-induced multidefect ordering, another competing route to defect lattices that turnover is shown not to require.","marker":"[31]"}],"fun_headline_variants":["Turnover stabilizes defect pairs in compressible active polar fluids","Defect lattices and foams arise from turnover in active fluids","Turnover repels defects, yielding active foams and vortex glasses","Density-polarity coupling stabilizes topological defects in active fluids","Turnover turns defect annihilation into stable defect phases"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole stabilization mechanism requires that polarity order grows with density (the coupling $\\chi>0$), so that a vanishing polarity at a defect core leaves the fluid locally depleted; without that density contrast, turnover produces no outward flow and nothing stops the defects from annihilating.","fun_headline_variants_meta":{"raw":{"variants":["Turnover stabilizes defect pairs in compressible active polar fluids","Defect lattices and foams arise from turnover in active fluids","Turnover repels defects, yielding active foams and vortex glasses","Density-polarity coupling stabilizes topological defects in active fluids","Turnover turns defect annihilation into stable defect phases"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000497,"raw_usage":{"total_tokens":2432,"prompt_tokens":940,"completion_tokens":1492,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":1417}},"tokens_in":556,"tokens_out":1492,"duration_ms":11956,"temperature":1.0,"reasoning_tokens":1417,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:55:01.739669+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A numerical experiment with the same equations but $\\chi=0$ should show no defect stabilization: opposite-charge pairs should annihilate for all turnover rates. Equivalently, in a simulation with $\\chi>0$, measure the density profile and velocity field around an isolated $+1$ defect: if the claim is correct, the density should dip at the core and the radial velocity should grow with turnover rate roughly as $\\tau^{-1}R\\varepsilon$, while a defect pair should settle at a finite separation that shrinks as $\\chi\\to0$.","supporting_citations":[{"cited_title":"Defect order in active nematics on a curved surface","cited_arxiv_id":"2002.06364","evidence_quote":"Provides the two-component active-gel description of the actin cortex and the density-polarity coupling used as the starting hydrodynamic model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the overdamped defect-position dynamics with effective friction coefficient that the pair velocity balance Eq. (18) uses."},{"cited_title":"Maroudas-Sacks and K","cited_arxiv_id":null,"evidence_quote":"Provides the generic constitutive equations of polar active gels underlying Eqs. (6)-(9)."},{"cited_title":"de Gennes and J","cited_arxiv_id":null,"evidence_quote":"Supplies the Cross-Hohenberg classification used to distinguish the type I instability found with turnover from the type II instability without it."},{"cited_title":"Blanch-Mercader, P","cited_arxiv_id":null,"evidence_quote":"Prior work stabilizing defect lattices via substrate friction, the alternative mechanism this paper's turnover mechanism is compared against for lattice formation."},{"cited_title":"Maroudas-Sacks, L","cited_arxiv_id":null,"evidence_quote":"Shows friction-induced multidefect ordering, another competing route to defect lattices that turnover is shown not to require."}],"review_version":1}