{"id":"b0880cc9-67d0-4933-8760-b9291d394d5a","arxiv_id":"2607.07656","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":7,"one_line_summary":"A minimal single-field nematopolar model exhibits string-mediated defect interactions, (t/ln t)^{1/2} coarsening, and motility-induced charge symmetry breaking under self-advection.","lead":"This paper introduces a minimal single-field model for systems with both polar and nematic alignment, showing it produces depolarization strings, loop dynamics, and arrested coarsening under activity. A smart generalist might read it to understand how competing symmetries generate novel topological defect structures relevant to biological and synthetic active matter.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Arrested coarsening plateau lacks finite-size scaling verification; MICS stability arguments are purely qualitative without linear stability analysis.","rationale":"The reader correctly identifies that the analytical arguments rely on unspecified heuristic functions (g(k̄), h(k̄)) and that code is not shipped, both of which are real gaps. However, the reader's framing of the load-bearing concern as 'model minimality' (neglecting hydrodynamics, density, noise) is somewhat generic — this is a standard modeling choice in the dry active matter literature, and the paper is transparent about it. The more specific and consequential concern is the absence of finite-size scaling for the arrested coarsening plateau, which directly tests whether the paper's most novel claim (MICS + arrest) is genuine or a finite-size artifact. Without this check, the CONDITIONAL verdict is appropriate, but for a more precise reason than the reader states. The passive coarsening results are reasonably well-supported and do not require verdict adjustment. The paper's strengths include clear presentation, systematic parameter variation in the passive case, and internally consistent observations (defect counts matching correlation-length scaling). The weaknesses are the lack of analytical stability analysis for the active case, the single-k̄ exploration of MICS, and the missing finite-size control. None of these rise to the level of REJECT — the observations are real and interesting — but they prevent full ACCEPT.","tokens_in":20231,"tokens_out":3470,"duration_ms":332357,"concrete_test":"Run the active case (k̄=0.3, Λ=0.1) at three system sizes — N=256, N=512, N=1024 — and extract the plateau value of L_n(t). If the plateau height scales linearly with N, the arrest is a finite-size artifact. If it is size-independent (converges to the same value in units of ℓ_nat), the arrested coarsening claim is strengthened. Additionally, perform a linear stability analysis of the aster defect profile (Eq. B6 with q=+1) perturbed by the advection term Λ(p·∇)p to confirm that outward asters are linearly stable and inward asters unstable.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most novel claim — motility-induced charge symmetry breaking (MICS) and arrested coarsening — rests on two weak supports. First, the stability arguments in Section VI are entirely qualitative: the paper asserts that advection 'transports polarization from the core, promoting its depletion and thereby stabilizing' outward asters (for Λ>0), and that quatrefoil structures are 'incompatible with advection' because polarization is 'unevenly redistributed.' No linear stability analysis of defect configurations under self-advection is performed, so we cannot confirm that the observed defect selection is a genuine dynamical attractor rather than a transient preference at the chosen parameter values. Second, and more concretely, the evidence for arrested coarsening is the plateau in L_n(t) shown in Figure 6(d), but all active-case simulations use a single system size (N=512 or 1024 at fixed k̄=0.3). If the plateau value of L_n scales with system size L, then what appears as 'arrest' is simply finite-size saturation — the system has run out of room for further coarsening. The paper does not check whether the plateau height is independent of L. This is the standard control for distinguishing genuine dynamical arrest from finite-size effects, and its absence is the single most consequential gap. The passive coarsening claim ((t/ln t)^{1/2}) is comparatively well-supported by the scaling collapse and defect-count data in Figures 4–5, though the unspecified functions g(k̄) and h(k̄) in Equations 8–14 weaken the analytical rationalization of the controlled-configuration results without undermining the main coarsening claim.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This paper introduces a minimal single-field dry model for nematopolar systems, where a polar vector field evolves under relaxational dynamics with competing polar and nematic free energy contributions, plus a self-advection term for activity. The authors study phase ordering and defect dynamics in both passive and active regimes. In the passive case, they characterize depolarization strings and loops, finding non-monotonic string-mediated interactions and a coarsening law L(t)~(t/ln t)^{1/2}. In the active case, they report motility-induced charge symmetry breaking (MICS), where positive integer and negative half-integer defects coexist, and coarsening is arrested. The passive coarsening results are well-supported by scaling collapses and defect-count data. The active-case claims of MICS and arrested coarsening are interesting but rest on qualitative stability arguments and lack a finite-size scaling control.","tokens_in":20537,"tokens_out":1336,"duration_ms":219615,"significance":"The paper addresses a timely problem in active matter theory: the ordering dynamics of systems with competing polar and nematic symmetries. The single-field formulation is a genuinely minimal framework that cleanly separates the roles of competing symmetries and activity. The identification of distinct loop relaxation mechanisms (shrinking vs. rotational) with a phase diagram is a solid contribution. The passive coarsening law is supported by systematic scaling collapses and defect-count kinetics. The MICS phenomenon and arrested coarsening, if confirmed, represent a novel and falsifiable prediction for active nematopolar systems. The main weakness is that the most novel active-case claims lack the quantitative rigor of the passive-case results.","major_comments":[{"comment":"Section VI, Figure 6(d): The central claim of arrested coarsening rests on the plateau in L_n(t). However, all active-case simulations appear to use a single system size (N=512 or 1024 at fixed k̄=0.3). If the plateau value scales with system size L, the observed 'arrest' could simply be finite-size saturation. A finite-size scaling check — showing that the plateau height is independent of L — is the standard control for distinguishing genuine dynamical arrest from finite-size effects. Its absence is the most consequential gap in the manuscript.","section":null},{"comment":"Section VI: The stability arguments for MICS are entirely qualitative. The paper asserts that advection 'transports polarization from the core, promoting its depletion and thereby stabilizing' outward asters (for Λ>0), and that quatrefoil structures are 'incompatible with advection' because polarization is 'unevenly redistributed.' No linear stability analysis of defect configurations under self-advection is performed. Without such analysis, we cannot confirm that the observed defect selection is a genuine dynamical attractor rather than a transient preference at the chosen parameter values. At minimum, the authors should clarify the expected basin of attraction or provide evidence of robustness across a wider parameter range.","section":null},{"comment":"Section IV, Eqs. (7)-(14): The energetic arguments use unspecified functions g(k̄) and h(k̄), defined only by monotonicity properties, to fit qualitative trends. For instance, the 'prediction' of d_eq in Eq. (10) reduces to ℓ_nat / h(k̄), where h is an undefined function. While the scaling-level arguments are internally consistent, the claim that Eq. (10) 'predicts' the decreasing trend of d_eq with k̄ is circular: h(k̄) is defined to be increasing, so 1/h(k̄) decreases by construction. The authors should either derive h(k̄) explicitly or reframe these as consistency checks rather than predictions.","section":null}],"minor_comments":[{"comment":"Section II.A: The statement 'simulations are run for ~10^6 iterations' should specify the physical time range covered, as the iteration count alone is not directly interpretable without Δt.","section":null},{"comment":"Figure 2(b): The error bars (if any) on the annihilation velocity v are not discussed. Given that v is estimated from a ratio involving a single annihilation time, clarification of the statistical methodology and number of realizations would help.","section":null},{"comment":"Section V, Figure 5: The defect count curves are stated to be 'averaged over 5 independent runs,' while correlation functions in Figure 4 are 'averaged over 10 independent runs.' The reason for the different sample sizes is not given.","section":null},{"comment":"Section VI, Figure 6(d) inset: The statement that 'the dynamical scaling hypothesis remains valid when distances are rescaled using these plateauing L_n(t)' is unclear. If L_n(t) plateaus, rescaling by a constant does not test dynamic scaling in the usual sense. This point needs clarification.","section":null},{"comment":"Appendix B, Eq. (B8): The linearization leading to the exponential decay of δp neglects 1/r^n terms. It would help to state explicitly the range of r over which this approximation is expected to hold, and whether this range covers the defect separations studied in Section IV.","section":null},{"comment":"The term 'motility-induced charge symmetry breaking (MICS)' is introduced without discussion of whether analogous symmetry-breaking mechanisms have been reported in related active nematic or polar systems. A brief comparison would contextualize the novelty.","section":null}],"recommendation":"major_revision","confidential_remarks":"The reader's report and stress-test note correctly identify the finite-size scaling gap as the most consequential issue. This is a genuine concern: without demonstrating that the plateau in L_n(t) is system-size independent, the arrested coarsening claim cannot be distinguished from finite-size saturation. The qualitative nature of the MICS stability arguments is also a real weakness, though one that could potentially be addressed by a linear stability analysis or at least a more systematic parameter sweep. The circularity concern regarding g(k̄) and h(k̄) is valid but somewhat less severe — the scaling arguments are internally consistent and the authors do not overclaim quantitative predictive power. I would recommend major revision with the finite-size scaling check as a required fix. The passive-case results are largely sound and could form the basis of a solid paper even if the active-case claims need to be tempered."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The referee correctly identifies that the passive-case results are quantitatively well-supported while the active-case claims need additional controls. We address each major comment below and commit to revisions.","responses":[{"response":"The referee is correct that this is the most important control missing from the active-case analysis. We acknowledge that the current data cannot distinguish genuine dynamical arrest from finite-size saturation on the basis of a single system size alone. We will perform finite-size scaling: we will run active-case simulations at several system sizes (N = 512, 768, 1024, and 1536) at fixed k̄ = 0.3 and Λ = 0.1, and verify whether the plateau value of L_n(t) is independent of L. If the plateau height is system-size-independent, this confirms genuine arrest; if it scales with L, we will revise the claim accordingly. We will add a new figure showing L_n(t) for multiple system sizes and a collapse plot (L_n/L versus t/L^2 or similar). We agree this is essential and will include it in the revised manuscript.","revision_made":"yes","referee_comment":"Section VI, Figure 6(d): The central claim of arrested coarsening rests on the plateau in L_n(t). However, all active-case simulations appear to use a single system size (N=512 or 1024 at fixed k̄=0.3). If the plateau value scales with system size L, the observed 'arrest' could simply be finite-size saturation. A finite-size scaling check — showing that the plateau height is independent of L — is the standard control for distinguishing genuine dynamical arrest from finite-size effects. Its absence is the most consequential gap in the manuscript."},{"response":"We agree that the stability arguments as currently presented are qualitative and that a linear stability analysis would substantially strengthen the MICS claim. A full linear stability analysis of defect configurations under self-advection is technically involved — the defect core structure is not a simple perturbation around a uniform state — and we are not confident it can be completed within the revision timeframe. However, we can and will provide stronger numerical evidence of robustness: (1) we will show that MICS persists across a wider range of advection strengths Λ and balance parameters k̄, reporting the defect composition ratio n_def^{+1}/n_def^{-1/2} as a function of both parameters; (2) we will verify that the defect selection is not transient by showing that the plateau in defect counts is stable over times at least an order of magnitude longer than the transient; (3) we will test robustness to different initial conditions (including ordered initial states with imposed defect configurations). We will also soften the language from 'stabilizing' to 'dynamically selecting' where appropriate, and explicitly state that a linear stability analysis is deferred to future work. We acknowledge that without the analytical treatment, the MICS mechanism remains a numerically observed phenomenon rather than a rigorously established dynamical attractor.","revision_made":"partial","referee_comment":"Section VI: The stability arguments for MICS are entirely qualitative. The paper asserts that advection 'transports polarization from the core, promoting its depletion and thereby stabilizing' outward asters (for Λ>0), and that quatrefoil structures are 'incompatible with advection' because polarization is 'unevenly redistributed.' No linear stability analysis of defect configurations under self-advection is performed. Without such analysis, we cannot confirm that the observed defect selection is a genuine dynamical attractor rather than a transient preference at the chosen parameter values. At minimum, the authors should clarify the expected basin of attraction or provide evidence of robustness across a wider parameter range."},{"response":"The referee is correct. The functions g(k̄) and h(k̄) are not derived from first principles; they are defined by monotonicity properties chosen to be consistent with the observed physics. The claim that Eq. (10) 'predicts' the decreasing trend of d_eq with k̄ is indeed circular as stated, since h(k̄) is defined to be increasing. We will reframe the language throughout Section IV: these are scaling-level consistency checks that show the energetic framework is compatible with the observed trends, not first-principles predictions. We will replace 'predicts' with 'is consistent with' and add an explicit statement that the functional forms of g(k̄) and h(k̄) are not derived but are constrained by physical requirements (vanishing at k̄=0, monotonicity). We will also note that the non-trivial content of the argument lies not in the monotonicity of h itself but in the identification of the competing energy scales (interaction vs. string tension) whose balance determines d_eq, and in the fact that the natural-unit rescaling collapses the data. The referee's point is well-taken and we will adjust the manuscript accordingly.","revision_made":"yes","referee_comment":"Section IV, Eqs. (7)-(14): The energetic arguments use unspecified functions g(k̄) and h(k̄), defined only by monotonicity properties, to fit qualitative trends. For instance, the 'prediction' of d_eq in Eq. (10) reduces to ℓ_nat / h(k̄), where h is an undefined function. While the scaling-level arguments are internally consistent, the claim that Eq. (10) 'predicts' the decreasing trend of d_eq with k̄ is circular: h(k̄) is defined to be increasing, so 1/h(k̄) decreases by construction. The authors should either derive h(k̄) explicitly or reframe these as consistency checks rather than predictions."}],"tokens_in":20083,"tokens_out":1645,"duration_ms":75188,"standing_objections":["We cannot provide a full linear stability analysis of defect configurations under self-advection. The defect core structure involves non-perturbative spatial profiles of the polarization field, and a systematic linearization around such non-uniform backgrounds is a non-trivial calculation that is beyond the scope of this paper. We will provide the strongest numerical evidence we can (parameter scans, long-time stability checks, initial-condition robustness) but acknowledge that the analytical treatment remains open."]},"desk_editor":{"model":"glm-5.2","letter":"Here's my read on the Aprile et al. nematopolar paper. The bottom line: this is a solid numerical study with genuinely new phenomenology, but the most novel active-case claim needs a control that isn't there yet. It deserves a serious referee. What's new and good: the single-field dry nematopolar model is a clean construction — one vector field, competing polar and nematic free-energy terms, self-advection for activity. The passive results are the strongest part. The (t/ln t)^{1/2} coarsening law is well-supported by the scaling collapse in Figure 4 and the defect-count data in Figure 5, both consistent with point-defect-mediated coarsening in 2D. The controlled simulations of defect pairs and loops (Section IV) are nicely done — the non-monotonic annihilation velocity, finite equilibrium separation of like-charged defects, and the two loop-collapse mechanisms (shrinking vs. rotational) are all clearly presented with sensible energetic arguments. The loop phase diagram in Figure 3(g) is a useful summary. The depolarization strings and loops themselves are not new — Vafa & Doostmohammadi and others reported similar structures in two-field models — but the single-field context and the systematic characterization of their relaxation pathways is a genuine contribution. Now the soft spots. The stress-test concern about finite-size scaling for the arrested coarsening claim is the one that matters most. Figure 6(d) shows L_n(t) plateauing for several advection strengths, but all active-case simulations appear to use a single system size. If the plateau height scales with system size L, what looks like dynamical arrest is just the system running out of room. This is the standard control and its absence is conspicuous. The authors need to show the plateau is size-independent. Second, the MICS mechanism — preferential stabilization of +1 asters and -1/2 trefoils — is argued purely qualitatively. The flow-based stability arguments in Section VI are physically reasonable but there's no linear stability analysis of defect configurations under advection. I don't think this is fatal; the defect-count ratios in Figure 6(b) are suggestive. But calling it a dynamical attractor without checking stability more carefully is a stretch. On the unspecified functions g(k̄) and h(k̄): the reader flagged these as circular. I think that's slightly overstated. The functions are introduced as effective parametrizations of string morphology, defined by monotonicity properties, and used to rationalize qualitative trends. They're not pretending to be first-principles derivations. The arguments would be stronger with explicit forms, but the scaling relations for d_eq and the loop phase boundary are consistent with the data. This is a minor weakness, not a load-bearing one. Who benefits: researchers in active matter and soft-condensed matter who work on topological defects and phase ordering. The passive coarsening results are reliable; the active results are promising but need the finite-size check before they're conclusive. Recommendation: send to a serious referee. The paper is worth the effort — if the authors can show the arrest plateau is size-independent and tighten the MICS argument, this is a useful contribution.","headline":"Single-field nematopolar model reproduces depolarization strings and coarsening; active case shows MICS and arrested coarsening but lacks finite-size scaling check","tokens_in":21289,"tokens_out":742,"would_cite":true,"duration_ms":128609,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"One field, two symmetries: a minimal model for nematopolar ordering","keywords":[],"falsifier":"If adding density fluctuations, hydrodynamic coupling, or thermal noise to the model qualitatively changes the MICS mechanism or the arrested coarsening regime — for instance, by destabilizing the aster/trefoil defect pairs or restoring full coarsening — then the single-field dry description would be insufficient and the MICS phenomenon would be an artifact of the model's simplicity rather than a generic feature of active nematopolar systems.","tokens_in":20590,"feed_emoji":"🌀","tokens_out":883,"duration_ms":154965,"temperature":0.7,"pith_summary":"The paper proposes and numerically studies a minimal continuum model for nematopolar systems — materials where polar (head-tail asymmetric) and nematic (head-tail symmetric) alignment coexist and compete. The model uses a single vector field whose free energy contains both a polar elastic term (penalizing any deformation from parallel alignment) and a nematic elastic term (penalizing deviations from parallel-or-antiparallel alignment). Activity enters through a self-advection term in the spirit of dry active matter theory, where the polarization field simultaneously serves as the order parameter and the local velocity. The authors show that this single-field description reproduces the hallmark structures of nematopolar systems — depolarization strings connecting half-integer topological defects and closed depolarization loops separating domains of opposite polarization — and that the passive coarsening follows the (t/ln t)^{1/2} growth law known from two-dimensional systems with point-like defects. When self-advection is turned on, the system exhibits what the authors call motility-induced charge symmetry breaking (MICS): positive integer defects (outward asters) and negative half-integer defects (trefoils) are dynamically stabilized while their opposites are suppressed, leading to a persistent finite defect density and arrested coarsening.","feed_headline":"One field, two symmetries: a minimal model for nematopolar ordering","feed_subtitle":"Competing polar and nematic alignment in a single vector field reproduces string defects, loop dynamics, and activity-driven arrest — no","key_machinery":"The model's free energy (Eq. 1) combines a double-well potential fixing |p|=1, a polar elastic term proportional to k_p that penalizes all spatial gradients, and a nematic elastic term proportional to k_n built from the tensor P-hat = (pp^T - |p|^2 I / 2), which only penalizes deviations from parallel-or-antiparallel alignment. The ratio k_bar = k_p / k_n controls the balance between polar and nematic character. The dynamics (Eq. 3) couple relaxational dynamics (minimizing F) to a self-advection term Lambda (p · grad) p. The elementary objects are half-integer defects connected by depolarization strings (lines where |p| -> 0), closed loops of vanishing polarization, and — under advection — a","core_discovery":"The central discovery is that a single vector field with competing polar and nematic free-energy contributions plus self-advection is sufficient to reproduce the full nematopolar phenomenology: string-connected defect pairs, loop structures with two distinct relaxation pathways, (t/ln t)^{1/2} passive coarsening, and — crucially — a motility-induced charge symmetry breaking under activity that arrests coarsening by preferentially stabilizing positive integer and negative half-integer defects. The mechanism for MICS is dynamical rather than energetic: self-advection stabilizes defect morphologies whose polarization circulation is consistent with the advective transport direction, while defect","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Single-field model captures full nematopolar defect phenomenology","Self-advection breaks defect charge symmetry and arrests coarsening","Strings, loops, and arrested coarsening in a minimal nematopolar model","Activity stabilizes mixed integer and half-integer defects in nematopolars","Competing polar and nematic alignment yields non-monotonic defect interactions"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The claim that the model is 'minimal' and captures the essential physics rests on the assumption that a single vector field with purely relaxational dynamics plus self-advection is sufficient — deliberately neglecting hydrodynamic flows, density fluctuations, and thermal noise. Whether the observed motility-induced charge symmetry breaking and arrested coarsening are generic features of nematopolar systems or artifacts of this specific dry active matter limit depends on this","fun_headline_variants_meta":{"raw":{"variants":["Single-field model captures full nematopolar defect phenomenology","Self-advection breaks defect charge symmetry and arrests coarsening","Strings, loops, and arrested coarsening in a minimal nematopolar model","Activity stabilizes mixed integer and half-integer defects in nematopolars","Competing polar and nematic alignment yields non-monotonic defect interactions"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":715,"prompt_tokens":636,"completion_tokens":79,"prompt_tokens_details":null},"tokens_in":636,"tokens_out":79,"duration_ms":47931,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T03:32:06.662379+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If adding density fluctuations, hydrodynamic coupling, or thermal noise to the model qualitatively changes the MICS mechanism or the arrested coarsening regime — for instance, by destabilizing the aster/trefoil defect pairs or restoring full coarsening — then the single-field dry description would be insufficient and the MICS phenomenon would be an artifact of the model's simplicity rather than a generic feature of active nematopolar systems.","supporting_citations":[],"review_version":1}